- Open Access
Nano-mechanical characterization of the wood cell wall by AFM studies: comparison between AC- and QI™ mode
Plant Methods volume 13, Article number: 60 (2017)
Understanding the arrangement and mechanical properties of wood polymers within the plant cell wall is the basis for unravelling its underlying structure–property relationships. As state of the art Atomic Force Microscopy (AFM) has been used to visualize cell wall layers in contact resonance- and amplitude controlled mode (AC) on embedded samples. Most of the studies have focused on the structural arrangement of the S2 layer and its lamellar structure.
In this work, a protocol for AFM is proposed to characterize the entire cell wall mechanically by quantitative imaging (QI™) at the nanometer level, without embedding the samples. It is shown that the applied protocol allows for distinguishing between the cell wall layers of the compound middle lamella, S1, and S2 of spruce wood based on their Young’s Moduli. In the transition zone, S12, a stiffness gradient is measured.
The QI™ mode pushes the limit of resolution for mechanical characterization of the plant cell wall to the nanometer range. Comparing QI™- against AC images reveals that the mode of operation strongly influences the visualization of the cell wall.
Wood is a unique biological material with several levels of structural hierarchy from the nanometer level to the macroscale . In particular the organization on the nanoscale, within the cell wall, is still under debate [2, 3]. The cell wall has been studied with numerous high resolution techniques like Transmission Electron Microscopy, TEM [4, 5], Scanning Near-field Optical Microscopy (SNOM) , or Atomic Force Microscopy (AFM) [3, 7,8,9,10]. A comparison of these studies shows that the observed nanostructural organization of the cell wall is strongly influenced by the used technique, and the sample preparation. Most of the previously cited studies have focused on the structural arrangement of the S2 layer and its lamellar organization, but also the other layers like the S1 and S3 are of importance for the mechanical behavior of the entire cell wall composite. TEM measurements on stained samples elucidated a transition zone between the S1 and S2 layers, called S12, where the lignin concentration decreases [11,12,13,14]. It is well known that in case of intra-wall failure, cell walls rupture predominately between the S1 and the S2-layer , however mechanical characterization of this mechanically highly relevant zone is missing.
State of the art imaging of secondary cell walls by AFM is in Resonant Contact mode [15, 16] and more common Amplitude Controlled mode (AC mode, an intermittent contact mode) of embedded cells [7, 17, 18]. Recent technical developments in the field of AFM allow to conduct a mechanical characterization in addition to topography studies, without sample embedding. Such multichannel AFM studies have first been applied for the characterization of primary walls [19, 20], and little has been done on secondary cell walls yet. Peakforce Quantitative Nanomechanics (PFQNM) has been used to image bamboo fiber cell walls  and to characterize the shape of cellulose microfibrils in primary plant cell walls in water using a very sharp probe with a curvature below 5 nm . Adhesion force mapping was conducted to measure the inactivation of a freshly cut cell wall surface, and to study the influence of surface roughness and tip geometry [22, 23]. Recently, Arnould et al.  gave a proof of concept of high resolution force mapping on embedded flax fibers by studying the mechanical gradients using nanoindentation as a reference and additionally Muraille et al.  applied the protocol on poplar fiber cell walls, however without providing specific structural details.
The purpose of the present paper is to systematically explore the feasibility of using the Quantitative Imaging mode (QI™ mode), a force spectroscopy mode, which records a force–distance (FD) curve in every pixel, as a tool to characterize the wood cell wall with a specific focus on the transition from the middle lamella to the S2 layer on a cell wall cross-section of spruce. To be able to compare the structures imaged in QI™ mode to literature data, the cell wall was also scanned in AC mode. Thereby, structural information was revealed with high resolution, helping to better understand the underlying structure–property relationships of wood cell walls.
Spruce cube preparation
The cross-section of an air dried spruce cube (5 × 5 × 5 mm3) was prepared by means of a two-step polishing process. The plane of the sectioning was oriented perpendicular to the longitudinal fiber axis of the wood. A microtome with a steel knife (RM2255, Leica) was used to smoothen the surface under wet conditions. Following the protocol of Keplinger et al. , an ultramicrotome (Ultracut, Reichert-Jung) equipped with a Diatome Histo diamond knife was used to polish. The typical root-mean-square roughness of a 1 × 1 µm2 cross-section was less than 1.4 nm. The microfibril angle (MFA) was measured on a spruce sample, cut out of the same strip of wood to gain samples from longitudinally matched positions, by wide angle X-ray scattering on the longitudinal-transverse plane. A MFA of 6° was recorded, which points to a mature wood probe.
Quantitative imaging AFM
AFM imaging was performed using a NanoWizard 4 (JPK Instruments AG) in QI™ mode under controlled climatic conditions (temperature 20 °C, humidity 65%). As a cantilever, a non-contact cantilever (NCHR, Nano World, resonance frequency 320 kHz) with a silicon probe was used. The cantilever was calibrated with the contact-free method for a beam shaped cantilever by giving the environmental conditions and cantilever dimensions . The force constant was calibrated to be 30 N m−1 (±1 N m−1, n = 5) and the deflection sensitivity was determined to be 24 nm V−1 (±1 nm V−1, n = 5). The setpoint of the measurement was defined according to the cantilever stiffness (60 nN) and the z-length (50 nm) and the pixel time (12 ms) were set fixed, to ensure a similar velocity for each measurement. This resulted in an extend rate of 62.5 kHz (the extend rate controls the speed of the movement on the extend part of the FD curve), and an extend speed of 10.42 µm s−1.
The software extension Advanced QI™ mode was used to have full access to all FD curves. The mapping resolution was chosen to be 256 × 256 pixels and the scan size was set between 10 × 10 µm2 and 1 × 1 µm2 (theoretical resolution limit 3.9 nm). The AFM was operated in z-closed loop, therefore the nonlinearity and hysteresis of the piezo were corrected during the movement (channel: height measured) and FD curves had a constant speed. The baseline was adjusted to correct for any changes of the setpoint.
Prior and after a cell scan, the tip resolution at 0° and 90° scan direction was checked using a test specimen (Product No. 628-AFM tip and resolution test specimen, Pelcro, TED PELLA INC.). The tip radius was fitted according to the stiffness of polystyrene (2 GPa) after scanning a polymer test specimen (PS-LDPE-GS, Veeco Metrology Group) (Fig. 1).
The data was analyzed in the JPK image processing software (JPK Instruments AG). A line fit was applied to correct the height measure image. To avoid characteristic shadows around high objects, the areas with elevated features were excluded from the line fit by a region of interest. As a second step possible default lines were replaced with the average between the adjacent lines. Force curve batch processing was performed in the following order on the extend curve: (1) Calibration of V-deflection (sensitivity and spring constant). (2) Smoothing of force data. (3) Baseline subtraction for offset and tilt by defining the fit range from 100% to the snap in, normally around 10%. (4) Contact point determination: Calculates vertical tip position, corrects the height signal for the cantilever deflection. (5) The Young’s Modulus can be calculated from the slope of the force curve by applying a contact mechanics model. As suggested by the suppliers, in PFQNM (Bruker) the Young’s Modulus is determined with the Derjaguin–Müller–Toporov (DMT) model of elastic contact on the retract curve , whereas in QI™ the extend curve is fitted to avoid any influence from plastic deformation. The Young’s Modulus was calculated with the DMT model, assuming a Poisson’s ratio of 0.4 (adapted from Gibson and Ashby ), and the respective radius of the cantilever. The quality of the fit was inspected visually for a typical FD curve. Due to the high amount of curves generated within one image, the fit cannot be optimized for each curve, but the weight of a false processed curve is therefore limited. The results were plotted in a histogram revealing the distribution of the data.
To remove cutting artefacts, the image data can be corrected with a two-dimensional fast Fourier transform (FFT) using Gwyddion 2.47.
AC mode AFM
In AC mode the mapping resolution was increased to 512 × 512 pixels, taking the equivalent time as the QI™ image. The full z-range of 15 µm was used. For the calibration of the cantilever the setpoint amplitude was selected to be around 70%. The approach was performed with constant velocity and baseline update at the starting point. The gain parameters were optimized in respect to the offset of the trace and retrace line in the oscilloscope. The channels were post-processed according to the height measured image described for the QI™ mode.
Results and discussion
Imaging cell walls in QI™ mode
By using the QI™ mode it is possible to scan over a whole unembedded cell through the lumen, because for FD curves the cantilever only moves the z-length that was set, in this case 200 nm (Fig. 2). Although the resolution of such a large scan (18 µm × 31 µm) is reduced in x-direction to 35 nm and in y-direction to 60 nm, stiffness changes between the stiff S2 and the softer middle lamella and S1 layer regions can clearly be seen in the Young’s Modulus image (Fig. 2b). Furthermore, the quality of the cut could be judged based on the visible inclination at the lumen/cell wall interface in the height image (Fig. 2a, black arrow).
The middle lamella reveals no clear transition to the adjacent primary wall. Therefore, the middle lamella and both adjacent primary walls are termed compound middle lamella, CML. The height images of the CML appear different depending on the scanning direction, which can be perpendicular, or in line to the cell wall layers (Fig. 3). This is a typical artefact in AFM imaging when scanning oriented structures , and needs to be considered when comparing different positions in the cell wall. Interestingly, the corresponding Young’s Modulus images are not affected by the scanning direction.
Comparison of QI™- and AC mode
In QI™ mode the algorithm of the tip motion measures a FD curve in every pixel, with a defined setpoint . Thereby, besides topological information the sample can be characterized mechanically with a high spatial resolution at high speed. There is no xy movement during the FD curve recording, which ensures a measurement under constant velocity. In AC mode the scanning cantilever is oscillating at a kilohertz range frequency, and as an additional channel the phase image can be displayed. The lock-in-amplifier measures a phase shift between the drive signal and the cantilever movement in dependency of the tip-sample interaction, including mechanical information, adhesion, and dissipation of cantilever energy . In order to compare the imaging mode QI™ to the state of the art applied AC mode, we choose to scan with both modes the area of transition from the compound middle lamella to the S2.
Figure 4a shows an overview image of a cell corner and Fig. 4b–d display a zoom into the transition zone from the middle lamella to the S2-layer. The vertical straight lines in the images of the S2 are cutting artefacts that result from imperfections of the diamond knife. The CML is buildup of an isotropic structure consisting mainly of lignin, that is assumed to organize in a self-assembly process [31, 32]. After the CML follows a narrow zone of approximately 100 nm, most presumably the S1. While the S1 has a comparable denser structure, the S2 appears as a woven network. The typical lamella structure of the S2 cannot be clearly visualized because, (1) cutting artefacts may overly the lamella, (2) the scanned area might be too close to the CML, and/or (3) microfibrils are strictly parallel aligned .
The Young’s Modulus of the S1 was determined by selecting a rectangle of maximum size inside the cell wall layer and calculating the root mean square value (Fig. 4g). Due to its large MFA, typically 70–90°, the S1 possesses a lower stiffness than the S2 (Fig. 4e) . Surprisingly, the S1 stiffness is with 1.3 GPa lower than the CML (data obtained from a region of interest comprising approximately 5625 FD curves). This might be explained by the different textures of the loose CML and the denser S1 layer, leading to artificially higher peak values in the CML. This indicates that due to the specific surface- cantilever tip interactions the structural patterns of the cell wall layers affect the stiffness values and that the stiffness ratios between the cell wall layers can be altered.
In the corresponding stiffness images and histograms it can be seen, how the Young’s Modulus changes from a narrow Gaussian distribution mean value of around 2.5 ± 0.6 GPa in the CML (Fig. 4j) to a broader distribution around 3.2 ± 0.8 GPa in the S2 layer (Fig. 4h, for a comparison the values are summarized in Fig. 4k, data obtained from one image comprising approximately 65536 FD curves). The histogram distribution of the S2 might be larger, because of being more heterogeneous in terms of biopolymer composition compared to the middle lamella, which is mainly composed of lignin . The above-mentioned transition zone between the S1 and S2, called S12, where the lignin concentration decreases [11, 12], can be visualized in the QI™ mode over a length of 2 µm by an increase in Young’s Modulus from the S1 to the S2 (Fig. 4f,i). Depending on the depth information measured by AFM, the increase in Young’s Modulus can arise from (1) the stiff cellulose microfibrils oriented with a small MFA, (2) the change in the surface structure due to a change in MFA, and/or (3) the change in composition of the different cell wall layers. Further experiments on composite model systems, with defined MFA on a similar length scale than wood, are required to make a certain statement about the parameters influencing the measured Young’s Modulus. Certainly, very low stiffness values for the wood cell wall, respectively its individual layers were measured in comparison to mechanical data obtained by tensile tests [34,35,36] and nanoindentation [37, 38]. The entirely different loading conditions and test geometries in micro- and macroscopic tests do not allow for a direct comparison with the AFM data. Since the cell wall stiffness values obtained by nanoindentation and those derived from micro- and macroscopic tests are in the same range, a focus is laid on a comparison with values measured by nanoindentation, which was recently discussed by Arnould and Arinero .
The contact mechanism of the probe of a nanoindenter considerably differs from an AFM tip in terms of geometry and indentation depth, which results in different interactions with the wood surface, because of the anisotropy of wood [39, 40]. Comparing stiffness values obtained by different AFM modes, Arnould et al.  applied PFQNM on a flax fiber with a setpoint of 200 nN on the retrace curve. For the S1 and CML they calculated a Young’s Modulus of 7 GPa. For the S2 they obtained values from 13 to 18 GPa, which is comparable to the values they measured with a nanoindenter. Muraille et al.  applied PFQNM on a transverse section of poplar plant fiber cell wall with a setpoint of 600 nN. The indentation moduli of the CML, S1 and S2 were averaged 17, 21, and 26 GPa. In our measurement the penetration depth is very small as at a setpoint of 60 nN, the penetration is around 5 nm. Thereby, only surface properties are determined and no plastic deformation takes place. The setpoints mentioned in the two studies before are three to ten times higher than in our measurement, therefore it is assumed that also the penetration depth was larger. As wood is a viscoelastic material, the scan speed of the measurements also needs to be taken into account, when comparing the obtained values. In QI™ mode the supplier suggests to use the trace curve, whereas in PFQNM it is suggested to use the retrace curve, for the fitting routine of the contact model. In our measurements, the slope of the retrace curve was too steep to be fitted. The trace curves could be fitted with the Hertz-  and the DMT model , with the DMT model giving higher stiffness values. The roughness of the microtome polished surface was very small. For the CR-AFM and the PFQNM measurement of Arnould et al.  the Young’s Modulus increases from the CML to the S2 by a factor of two. Although the values in our measurements were smaller, the same tendency was found.
To compare the QI™ mode with the AC mode, Fig. 5 displays the transition from the CML to the S2 imaged in AC mode. The height images in the two modes are comparable (Fig. 5b–d). The phase image shows lignin in the CML as ellipsoid structures and the adjacent S1 layer possess a low phase signal (Fig. 5e–g). In phase contrast images, typically darker areas can be correlated with regions of lower stiffness . By applying a FFT the directed cutting artefacts can be selectively removed from the phase image, thereby the structural change from the S1 layer to the S2 layer becomes more obvious (Fig. 5h–j). The transition zone S12 shows a parallel orientation to the CML (white arrow), whereas the S2 has a granular structure. The transition zone visualized in the phase image is smaller than the one from the Young’s Modulus image, as the Modulus image is more sensitive for detecting changes in mechanical properties.
Fahlén and Salmén  visualized individual cellulose aggregates ranging from 10 to 30 nm in the S2 layer by phase imaging oriented with a regularity that was interpreted as an additional lignin pattern. Here, we fitted the structures in the FFT of the S2 and they lay in the same size range compared to the one observed by Fahlén and Salmén . The lignin structures of the CML had an ellipsoid shape ranging from 30 to 60 nm diameter. Due to the influence of the cutting artefacts a clear lamellar structure in the S2 could not be detected.
Comparison fast Fourier transform
In Fig. 6 two enlarged FFT corrected images of the transition between the CML and S2 scanned in QI™- (Fig. 6a) and AC mode (Fig. 6b) can be seen. After FFT processing the variability of the Young’s Modulus values could be reduced to two colors (the lighter, the stiffer), thereby the gradient information is lost. One can only distinguish the CML from the S2. The CML appears as a porous network, whereas the S2 seems to be dense. The phase image has a twice as high resolution, therefore it appears sharper. The CML gives the impression of a negative of the Young’s Modulus image and the S2 has a granular appearance. The phase shift is proportional to the stiffness, but phase contrast interpretation is not as straight forward like the analysis of FD curves, due to the contributions from contact area, viscoelastic properties, and capillary forces .
In conclusion, QI™ mode gives the opportunity to mechanically characterize the different secondary cell wall layers on the nanometer level by obtaining FD curves that can be analyzed with a mechanical model. Although too low stiffness values are measured, it outperforms mechanical characterization by nanoindentation in terms of resolution, and provides a more distinct image of stiffness distribution than what can be obtained by phase contrast imaging in AC mode.
The imaging mode QI™ is on the one hand very robust, as it can scan over whole unembedded cells, and on the other hand it is sensitive enough to detect small changes, as can be seen from the visualization of the transition zone S12. The nanostructure of wooden cell walls is just at the beginning to be characterized mechanically. Comparative studies are needed to unravel the influence of never dried wood, different cell types, or different species on the mechanics. Further insights into the cell wall assembly can be used to understand the localization of chemical modifications within wood beyond the resolution limit of commonly used techniques , or to set up reliable computational models on the microscale [43,44,45].
Atomic Force Microscopy
contact resonance AFM
- FD curve:
fast Fourier transform
Peakforce Quantitative Nanomechanics
Fratzl P, Weinkamer R. Nature’s hierarchical materials. Prog Mater Sci. 2007;52:1263–334.
Donaldson LA. A three-dimensional computer model of the tracheid cell wall as a tool for interpretation of wood cell wall ultrastructure. IAWA J. 2001;22:213–33.
Zimmermann T, Thommen V, Reimann P, Hug HJ. Ultrastructural appearance of embedded and polished wood cell walls as revealed by atomic force microscopy. J Struct Biol. 2006;156:363–9.
Ruel K, Barnould F, Goring DAI. Lamellation in the S2 layer of softwood tracheids as demonstrated by scanning transmission electron microscopy. Wood Sci Technol. 1978;12:287–91.
Kerr AJ, Goring DAI. The ultrastructural arrangement of the wood cell wall. Cellul Chem Technol. 1975;9:563–73.
Keplinger T, Konnerth J, Aguié-Béghin V, Rüggeberg M, Gierlinger N, Burgert I. A zoom into the nanoscale texture of secondary cell walls. Plant Methods. 2014;10:1.
Fahlén J, Salmén L. On the lamellar structure of the tracheid cell wall. Plant Biol. 2002;4:339–45.
Schwarze FWMR, Engels J. Cavity formation and the exposure of peculiar structures in the secondary wall (S2) of tracheids and fibres by wood degrading Basidiomycetes. Holzforschung. 1998;52:117–23.
Zimmermann T, Eckstein JSD. Rasterelektronenmikroskopische Untersuchung an Zugbruchflächen von Fichtenholz. Holz Roh Werkst. 1994;52:223–9.
Muraille L, Aguié-Béghin V, Chabbert B, Molinari M. Bioinspired lignocellulosic films to understand the mechanical properties of lignified plant cell walls at nanoscale. Sci Rep. 2017;7:1–11.
Fromm J, Rockel B, Lautner S, Windeisen E, Wanner G. Lignin distribution in wood cell walls determined by TEM and backscattered SEM techniques. J Struct Biol. 2003;143:77–84.
Brändström J, Bardage SL, Daniel G, Nilsson T. The structural organisation of the S1 cell wall layer of norway spruce tracheids. IAWA J. 2003;24:27–40.
Reza M, Ruokolainen J, Vuorinen T. Out-of-plane orientation of cellulose elementary fibrils on spruce tracheid wall based on imaging with high-resolution transmission electron microscopy. Planta. 2014;240:565–73.
Donaldson LA. Cell wall fracture properties in relation to lignin distribution and cell dimensions among three genetic groups of radiate pine. Wood Sci Technol. 1995;29:51–63.
Clair B, Arinero R, Lévèque G, Ramonda M, Thibaut B. Imaging the mechanical properties of wood cell wall layers by atomic force modulation microscopy. IAWA J. 2003;24:223–30.
Nair SS, Wang S, Hurley DC. Nanoscale characterization of natural fibers and their composites using contact-resonance force microscopy. Compos A Appl Sci Manuf. 2010;41:624–31.
Fahlén J, Salmén L. Cross-sectional structure of the secondary wall of wood fibers as affected by processing. J Mater Sci. 2003;38:119–26.
Fahlén J, Salmén L. Pore and matrix distribution in the fiber wall revealed by atomic force microscopy and image analysis. Biomacromol. 2005;6:433–8.
Zhang T, Zheng Y, Cosgrove DJ. Spatial organization of cellulose microfibrils and matrix polysaccharides in primary plant cell walls as imaged by multichannel atomic force microscopy. Plant J. 2016;85:179–92.
Peaucelle A, Braybrook SA, Le Guillou L, Bron E, Kuhlemeier C, Hofte H. Pectin-induced changes in cell wall mechanics underlie organ initiation in Arabidopsis. Curr Biol. 2011;21:1720–6.
Ren D, Wang H, Yu Z, Wang H, Yu Y. Mechanical imaging of bamboo fiber cell walls and their composites by means of peakforce quantitative nanomechanics (PQNM) technique. Holzforschung. 2015;69:975–84.
Frybort S, Obersriebnig M, Muller U, Gindl-Altmutter W, Konnerth J. Variability in surface polarity of wood by means of AFM adhesion force mapping. Colloid Surf A. 2014;457:82–7.
Jin X, Kasal B. Adhesion force mapping on wood by atomic force microscopy: influence of surface roughness and tip geometry. R Soc Open Sci. 2016;3:160248.
Arnould O, Siniscalco D, Bourmaud A, Le Duigou A, Baley C. Better insight into the nano-mechanical properties of flax fibre cell walls. Ind Crop Prod. 2017;97:224–8.
Sader JE, Chon JWM, Mulvaney P. Calibration of rectangular atomic force microscope cantilevers. Rev Sci Instrum. 1999;70:3967–9.
Derjaguin BV, Muller VM, Toporov YP. Effect of contact deformations on adhesion of particles. J Colloid Interface Sci. 1975;53:314–26.
Gibson LJ, Ashby MF. Cellular solids. Structure and properties. Cambridge: Cambridge University Press; 2001.
Tsukruk VV, Bliznyuk VN, Visser D, Campbell AL, Bunning TJ, Adams WW. Electrostatic deposition of polyionic monolayers on charged surfaces. Macromolecules. 1997;30:6615–25.
Chopinet L, Formosa C, Rols MP, Duval RE, Dague E. Imaging living cells surface and quantifying its properties at high resolution using AFM in QITM mode. Micron. 2013;48:26–33.
Tsukruk VV, Singamani S. Scanning probe microscopy of soft matter: fundamentals and practice. 1st ed. London: Wiley; 2012.
Salmén L. Wood morphology and properties from molecular perspectives. Ann For Sci. 2015;72:679–84.
Salmén L, Olsson A-M, Stevanic JS, Simonovic J, Radotic K. Structural organisation of the wood polymers in the wood fibre structure. Bioresources. 2012;7:512–32.
Fengel D, Wegener G. Wood chemistry, ultrastructure, reactions. Berlin: Walter de Gruyter; 1984.
Gindl W, Schöberl T. The significance of the elastic modulus of wood cell walls obtained from nanoindentation measurements. Compos A Appl Sci Manuf. 2004;35:1345–9.
Page DH, El-Hosseiny F, Winkler K. Behaviour of single wood fibres under axial tensile strain. Nature. 1971;3:229–52.
Burgert I, Keckes J, Frühmann K, Fratzl P, Tschegg SE. A comparison of two techniques for wood fibre isolation- evaluation by tensile tests on single fibres with different microfibril angle. Plant Biol. 2002;4:9–12.
Jäger A, Hofstetter K, Buksnowitz C, Gindl-Altmutter W, Konnerth J. Identification of stiffness tensor components of wood cell walls by means of nanoindentation. Compos A Appl Sci Manuf. 2011;42:2101–9.
Wimmer R, Lucas BN, Tsui TY, Oliver WC. Longitudinal hardness and Young’s modulus of spruce tracheid secondary walls using nanoindentation technique. Wood Sci Technol. 1997;31:131–41.
Arnould O, Arinero R. Towards a better understanding of wood cell wall characterisation with contact resonance atomic force microscopy. Compos A Appl Sci Manuf. 2015;74:69–76.
Gindl W, Gupta HS, Schöberl T, Lichtenegger HC, Fratzl P. Mechanical properties of spruce wood cell walls by nanoindentation. Appl Phys A. 2004;79:2069–73.
Sneddon N. The relation between load and penetration in the axisymmetric boussinesq problem for a punch of arbitrary profile. Int J Eng Sci. 1985;3:47–57.
Gierlinger N, Keplinger T, Harrington M. Imaging of plant cell walls by confocal Raman microscopy. Nat Protoc. 2012;7:1694–708.
Barber NF, Meylan BA. The anisotropic shrinkage of wood: A theoretical model. Holzforschung. 1964;18:146–56.
Mishnaevsky L, Qing H. Micromechanical modelling of mechanical behaviour and strength of wood: state-of-the-art review. Comp Mater Sci. 2008;44:363–70.
Yamamoto H, Kojima Y. Properties of cell wall constituents in relation to longitudinal elasticity of wood. Wood Sci Technol. 2002;36:55–74.
KC designed the experiment and analyzed data. KC, TK and IB co-wrote the paper. All authors discussed results and commented on the manuscript. All authors read and approved the final manuscript.
The authors thank Thomas Schnider for cutting the wood samples. Valuable discussion with John Berg is acknowledged.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.