Super-resolution imaging with Pontamine Fast Scarlet 4BS enables direct visualization of cellulose orientation and cell connection architecture in onion epidermis cells
© Liesche et al.; licensee BioMed Central Ltd. 2013
Received: 23 October 2013
Accepted: 23 December 2013
Published: 28 December 2013
In plants, a complex cell wall protects cells and defines their shape. Cellulose fibrils form a multilayered network inside the cell-wall matrix that plays a direct role in controlling cell expansion. Resolving the structure of this network will allow us to comprehend the relationship of cellulose fibril orientation and growth.
The fluorescent dye Pontamine Fast Scarlet 4BS (PFS) was shown to stain cellulose with high specificity and could be used to visualize cellulose bundles in cell walls of Arabidopsis root epidermal cells with confocal microscopy. The resolution limit of confocal microscopy of some 200 nm in xy and 550 nm in z for green light, restricts the direct visualization of cellulose to relatively large bundles, whereas the structure of cellulose microfibrils with their diameter below 10 nm remains unresolved. Over the last decade, several so-called super-resolution microscopy approaches have been developed; in this paper we explore the potential of such approaches for the direct visualization of cellulose.
To ensure optimal imaging we determined the spectral properties of PFS-stained tissue. PFS was found not to affect cell viability in the onion bulb scale epidermis. We present the first super-resolution images of cellulose bundles in the plant cell wall produced by direct stochastic optical reconstruction microscopy (dSTORM) in combination with total internal reflection fluorescence (TIRF) microscopy. Since TIRF limits observation to the cell surface, we tested as alternatives 3D-structured illumination microscopy (3D-SIM) and confocal microscopy, combined with image deconvolution. Both methods offer lower resolution than STORM, but enable 3D imaging. While 3D-SIM produced strong artifacts, deconvolution gave good results. The resolution was improved over conventional confocal microscopy and the approach could be used to demonstrate differences in fibril orientation in different layers of the cell wall as well as particular cellulose fortifications around plasmodesmata.
Super-resolution light microscopy of PFS-stained cellulose fibrils is possible and the increased resolution over conventional approaches makes it a valuable tool for the investigation of the cell-wall structure. This is one step in method developments that will close the gap to more invasive techniques, such as atomic force and electron microscopy.
KeywordsFluorescent dye Cell wall Cellulose STORM Structured illumination Super-resolution microscopy TIRF Deconvolution
Cellulose microfibrils form the backbone of the complex cell wall of plant cells. The exact structure of the cellulose fibril network has important implications for our understanding of plant growth on a cellular level, as it has a direct role in controlling cell elongation [1–3]. In addition, precise knowledge of cell wall architecture facilitates the development of methods for efficient breakdown of plant cell walls in the production of biofuels [4, 5].
Pontamine Fast Scarlet 4BS (PFS), a fluorescent dye that binds cellulose with high specificity, has been used to visualize bundles of cellulose bundles in the cell wall of Arabidopsis root epidermis cells . With a confocal microscope, cellulose bundle orientation could be followed over time, providing the first direct proof of passive reorientation of cellulose bundles during cell elongation. Nevertheless, only the largest cellulose bundles could be visualized, as the resolution of conventional light microscopy is limited to around 250 nm in xy-direction. Cellulose microfibrils have a diameter below 10 nm [7, 8], but can be bundled into macrofibrils. Destructive methods that offer resolution on the nm scale such as scanning electron microscopy and atomic force microscopy indicate that macrofibrils do not necessarily show the same orientation as the majority of microfibrils [9, 10].
In order to achieve a more complete and dynamic understanding of the structure of the cellulose network in plant cell walls a resolution increase for light microscopy is required that allows showing the cell wall in its native state and its differentiation in living cells. The extensive preparation of invasive high resolution techniques, such as scanning electron microscopy (SEM) and atomic force microscopy (AFM) are not compatible with live microscopy and were discussed to have profound impact on the cellulose network structure [9, 11].
Recent developments in the field of light microscopy significantly increased the resolving capacity of fluorescence microscopes following different approaches (reviewed by Schermelleh et al. ). Basically, these approaches improve resolution by circumventing the diffraction-limits of microscopes as defined by Ernst Abbe . We used two of these so-called super-resolution microscopy techniques, 3D-structured illumination microscopy (3D-SIM) and direct stochastic optical reconstruction microscopy (dSTORM), to test visualization of PFS-stained cellulose fibrils. The third of the established super-resolution techniques, stimulated emission depletion microscopy (STED) is not compatible with green plant tissues, since the depletion laser wavelength is highly absorbed by chlorophyll in the sample. In comparison to wide field fluorescence microscopes, 3D-SIM and STORM improve the resolution of fluorescent signals by, respectively, 2 fold in xyz and by up to 10 fold in xy . In addition, we also tested the effects of deconvolution on confocal image stacks which leads to clearer micrographs by removing the distorting effects of the point-spread function. This function is an inherent factor of any optical system and can be compensated by appropriate algorithms.
Onion bulb scale epidermis, which is used in the present study, has been identified as an ideal model system for the investigation of how cells control growth . Knowing the exact structure of the cellulose network in the cell wall is therefore of immediate significance.
Results and discussion
Spectral properties and cell toxicity of Pontamine Fast Scarlet in situ
Direct stochastic optical reconstruction microscopy shows PFS-stained cellulose fibrils on the cell surface
Fundamentally, two factors are important for a good dSTORM dye: 1) The number of photons detected per dye molecule. This will determine the precision by which the dye molecules are localized. 2) The on-off duty cycle. Dempsey et al. defines the on-off duty cycle as: "the fraction of time a fluorophore spends in the on state". Even with perfectly localized dye molecules the duty cycle will influence the resolution of structures . In regard to these factors, PFS is not an ideal probe for dSTORM. The observed blinking events suggest that PFS, when bound to cellulose, generally has only a few switching cycles before irreversible photodestruction and has a long dark state, meaning that the duty cycle is very low. In addition, during the on-time of 5–10 ms only a relatively low amount of photons can be detected. In practice, this means that a large number of images have to be acquired in order to get a complete super-resolution image.
The method might be improved by employing an imaging buffer containing an oxygen-scavenging system and a primary thiol as reducing compound. Reducing conditions stabilize the off-state of fluorophores, thus preventing their irreversible photodestruction. At the same time, the buffer conditions induce spontaneous recovery of fluorophore molecules to their on-state. Optimization of thiol concentration and excitation laser intensity could lead to an increased number of localizations per frame and duty cycles. This should result in faster imaging and higher resolution .
The cellulose fibrils in Figure 3B have a diameter of around 100 nm with clear spacing between them and appear to branch at various positions. This structure on the cell surface can be significantly different from inner layers, which are generally expected to show narrower-spaced, more uniform fibrils . Nevertheless, the observed form and pattern reflect the structural changes achieved during cell expansion and can, as such, give indications of how this process happens. The dSTORM images obtained here could point to the possibility of increased aggregation of fibrils during restructuring, although data from other cell-wall layers obtained with a similar method would be necessary for comparison.
PFS is not suited for visualization of cellulose with 3D structured illumination
In 3D-SIM, patterned excitation light is used to take images at various phases, rotations and focus planes. From the resulting diffraction patterns, an image with higher resolution than wide-field images can be obtained via mathematical reconstruction in the Fourier space. Usually, this technique is well suited for imaging in plant tissue as it works with conventional probes and moderately thick specimens [21, 22]. It was, however, not possible to resolve cellulose fibrils in the outer cell wall of onion bulb scale epidermis cells clearly. The Additional file 1: Figure S1 shows a comparison of the PFS-stained outer wall epidermis in wide field and 3D-SIM mode with two different noise filter settings. While standard noise filtering results in loss of any structural information, minimal noise filtering results in typical rotation-symmetric 3D-SIM artefacts as seen in the center of Additional file 1: Figure S1. Only small parts of the image reveal a fibril-like structure. It cannot be ruled out that strongly ordered cell wall macrofibrils from other cell types, like sclerenchyma fibers, would be resolved by 3D-SIM. However, discrimination of artefacts from real structures appears difficult - in any case for the example used here.
The effects of different filter settings indicate that the background signal level in PSF-stained onion cells is too high. In this respect the 3D-SIM mode behaves like a wide-field microscope rather than a confocal microscope, where out-of focus fluorescence is strongly reduced already under recording.
Confocal imaging combined with image deconvolution can show the cellulose fibril network in 3D
While not usually considered as a super-resolution technique, image deconvolution, a mathematical operation to recover signal from blurring and noise, can also be used to enhance the resolution capacity of conventional microscopes beyond the diffraction limit. In comparison to STORM and 3D-SIM the gain in resolution is modest, but image deconvolution has the advantage of being applicable to confocal image data. This means that no special imaging hardware is required and that structures deep inside the specimen can be imaged due to the optical sectioning capacity of confocal microscopes.
Data can be acquired with point-scanning, as well as with spinning disc confocal systems. We used a spinning disc confocal, which offers a higher acquisition rate and lower photobleaching than a point confocal. After acquiring an image stack it is imported into adequate image processing software, of which the open source ImageJ and the commercial Huygens SVI are considered to give the best results . The critical factor for deconvolving images is the point spread function, which can be determined experimentally or theoretically. We chose theoretical calculation as the determination of experimental point-spread functions can be unreliable for thick specimen, like whole-tissue samples . Nevertheless, improvement in deconvolution efficiency might be realized if reference objects at the actual focal plane, for example fluorescent beads attached to the outside of a cell, are included during imaging and then used to calculate the sample-inherent point spread function.
Taken the optical resolution of some 200 by 550 nm (x by z) with the objective lens applied (N.A = 1.4) at the emission maximum of PFS (600 nm), the observed change of fibril orientation at a focal distance of 200 nm appears impossible. However, the strength of deconvolution lies exactly in its power to improve the 3-D resolution by removing the blurring effects of the zeppelin-shaped point-spread function. This technique would not allow to separate two points in the specimen lying on top of each other on the optical axis (= z-axis) within this distance, but is able to do so for the linear structures as the macrofibrils in the sample used here. This means that two microfibrils running parallel on top of each other with a distance of 200 nm in z cannot be separated, but those which are crisscrossing. Here the deconvolution algorithm is able to improve the contrast between the fiber in focus and the gap (no fluorescence) 200 nm below sufficiently for visualization.
In addition to cellulose orientation, deconvolution of images of PFS-stained onion epidermis cell wall also reveals features of the cell connection architecture (Figure 5). While raw data of co-staining with the callose-specific dye Anilin Blue suggests incongruent distribution domains of cellulose in callose in the cell wall (Figure 5A,B), this becomes much clearer on deconvolved images (Figure 5C, D), even showing the thin cell wall plate that remains in the pit fields. Interestingly, in 3D view or maximum intensity projections of deconvolved images, it becomes apparent that cellulose is deposited at much higher density in areas of 2 μm around callose compared to the rest of the wall (Figure 5F). This cellulose collar could potentially form a limit to the callose deposits that are generally enhanced by wounding in order to plug the plasmodesmata .
Wide-field fluorescence microscopy allows investigating cell walls in their native state, but further details of its substructure can be uncovered when improving the microscopic resolution power. Here we demonstrate that the cellulose-specific dye PFS is suited for use in TIRF microscopy and STORM, realizing the resolution of cellulose fibrils below 100 nm thickness. Nevertheless, since TIRF microscopy limits observation to the cell surface, no information can be gained from deeper cell wall layers. While 3D-SIM principally enables super-resolution imaging even more than 100 μm inside a biological specimen, the high signal density of PFS in the cell wall prevented a correct mathematical reconstruction of the fibrillar substructure in the onion epidermis. As an alternative to STORM and 3D-SIM we used conventional confocal microscopy with post-acquisition image deconvolution. This enables a more modest increase in resolution compared to the other approaches, but leverages the z-sectioning and deep imaging capacity of confocal microscopy. This approach yielded data that enabled us to assess the orientation of cellulose fibrils, which was different in two layers within the cell wall. It was also used in co-staining experiments with a dye for another cell wall component, which revealed new features of cell wall architecture around plasmodesmata between neighboring cells.
Future developments, like the adaptation of 3D-STORM  for plant tissues will, without doubt, help to further close the resolution gap between non-invasive light microscopy and invasive techniques like AFM and SEM. The results presented here show how a first step towards this goal can be realized with the super-resolution techniques that are currently available, demonstrating their potential to be a valuable tool for the investigation of cell wall structure of living plant cells at a resolution higher than what had been technically possible before.
Dyes and staining
Pontamine Fast Scarlet 4 BS was purchased from Aldrich Rare Chemicals Library (catalog no. S479896; Sigma-Aldrich, St. Louis, USA). It is now available under the synonym Direct Red. Stock solutions of 5 mg/ml in 0.1 M phosphate buffered saline (PBS) pH 7.2 can be stored at -4°C or -20°C for several weeks. The working solution was diluted to 0.1 – 0.01 mg/ml depending on the tissue and was always prepared fresh. Incubation time was usually 20 min, but longer in complex tissues. Staining was followed by brief washing in PBS.
Carboxyfluorescein-diacetate (Molecular Probes, Life Technologies Corporation, Carlsbad, USA) was used at 0.5 μg/ml in PBS. Tissue was stained in 10 to 15 minutes, followed by 5 minutes washing in phosphate buffer.
Aniline Blue (Sigma-Aldrich, St. Louis, USA) working solution was diluted to 50 μg/ml in PBS and staining time was 5 to 10 minutes.
Onion bulb scale epidermis was prepared from Allium cepa bulbs purchased in a supermarket. Using a razorblade 0.5 × 0.5 cm pieces were cut out from fresh scale and the epidermis peeled off from the inner side with tweezers.
A simultaneous excitation and emission scan (Lambda2-scan) was performed in situ on onion bulb scale epidermis cell walls stained with PFS using a Leica SP5-X confocal laser-scanning microscope (Leica Microsystems, Mannheim, Germany). A white light laser, calibrated to keep a constant output power over the available spectrum, was used as excitation source. Excitation scans were made every 10 nm between 470 nm and 620 nm. The AOBS-defined detection window of 10 nm width was shifted from excitation wavelength plus 15 nm (center position) to 745 nm for each excitation wavelength. Image analysis was performed in ImageJ.
Image acquisition and processing
For confocal imaging in the cell vitality test a Leica TCS SP2/MP (Leica Microsystems, Mannheim, Germany) was used. PFS and cFDA were excited by an Argon Laser at 488 nm. The following detection windows were defined: PFS: 560–605 nm, cFDA: 485–505 nm.
Super-resolution imaging was performed on a Zeiss Elyra PS.1 microscope (Zeiss, Jena, Germany) and an Andor iXion 860 EMCCD camera (Andor Technology, Belfast, Ireland) using a 561 nm solid state laser for excitation, a band pass filter with DM 561 nm, EM 570–650 nm.
For 3D-SIM 63 × 1.2 Zeiss Apochromat water-immersion and 63 × 1.4 Zeiss Apochromat oil-immersion objectives were used together with a 23 μm grating. Z-stacks were recorded with 5 phase-changes and 5 grating rotations for each section. The 3D-SIM wizard in the Zeiss Zen 2011 software was used for image processing with the following setting: SR frequency weighing 1, Baseline cut, theoretical PSF and noise filtering as specified in the image legend.
For STORM, the same microscope was used with a 100 × 1.46 Zeiss Apochromat oil-immersion objective. The TIRF angle was adjusted by focusing on the surface of the onion epidermis cells in epi-fluorescence mode, then switching to TIRF-mode and increasing the angle until cellulose fibers became visible (around 60°). For one image series, 20 000 images were acquired. These were processed inside the Zeiss Zen 2011 software using the PALM-wizard. The filter settings were chosen permissive compared to settings generally used when working with photoswitchable proteins with relatively low chi-square and background thresholds. False positives were filtered out by setting the minimum photon count to 100. From the vector map of localization precisions, the background-corrected point-spread functions were rendered with a localization precision of 75 nm. The on time and dark state interval of PFS fluorophores were determined in separate image series. Only a 128 × 128 pixel area of the camera sensor was read out in order to enable image acquisition times of about 3 ms per image.
Deconvolution of confocal images was done after recording optical sections with a spinning-disc confocal (Andor Revolution XD, 100 × 1.4 Olympus UPlanSApo objective, Andor iXion 897 EMCCD camera, 561 nm solid state laser, EX 564 nm, DM 561 nm, EM 570–650 nm). Deconvolution was performed in ImageJ using the plugins PSF Generator and DeconvolutionLab, both developed by Biomedical Imaging Group at the École Polytechnique de Lausanne . With PSF Generator, theoretical point spread functions were generated with the Born and Wolf algorithm , which takes account of numerical aperture, refractive indices, wavelength, pixel size and z-step spacing. The resulting point-spread function was used in the deconvolution algorithm developed by Richardson and Lucy  with 20 iterations. This Bayesian-based method is generally considered among the most appropriate for deconvolution of confocal data . The output data is compared to raw data without further image adjustments.
All microscopy was performed at the Center for Advanced Bioimaging (CAB) Denmark. We thank Michael Hansen, University of Copenhagen, for advice on super-resolution image processing.
- Baskin TI: Anisotropic expansion of the plant cell wall. Ann Rev Cell Dev Biol. 2005, 21: 203-222. 10.1146/annurev.cellbio.20.082503.103053.View ArticleGoogle Scholar
- Geitmann A: Mechanical modeling and structural analysis of the primary plant cell wall. Curr Opin Plant Biol. 2010, 13: 693-699. 10.1016/j.pbi.2010.09.017.PubMedView ArticleGoogle Scholar
- Suslov D, Verbelen JP: Cellulose orientation determines mechanical anisotropy in onion epidermis cell walls. J Exp Bot. 2006, 57: 2183-2192. 10.1093/jxb/erj177.PubMedView ArticleGoogle Scholar
- Himmel ME, Ding SY, Johnson DK, Adney WS, Nimlos MR, Brady JW, Foust TD: Biomass recalcitrance: engineering plants and enzymes for biofuels production. Science. 2007, 315: 804-807. 10.1126/science.1137016.PubMedView ArticleGoogle Scholar
- Keegstra K: Plant cell walls. Plant Phys. 2010, 154: 483-486. 10.1104/pp.110.161240.View ArticleGoogle Scholar
- Anderson CT, Carroll A, Akhmetova L, Somerville C: Real-time imaging of cellulose reorientation during cell wall expansion in Arabidopsis roots. Plant Phys. 2010, 152: 787-796. 10.1104/pp.109.150128.View ArticleGoogle Scholar
- Ha MA, Apperley DC, Evans BW, Huxham IM, Jardine WG, Vietor RJ, Reis D, Vian B, Jarvis MC: Fine structure in cellulose microfibrils: NMR evidence from onion and quince. Plant J. 1998, 16: 183-190. 10.1046/j.1365-313x.1998.00291.x.PubMedView ArticleGoogle Scholar
- Kennedy CJ, Cameron GJ, Sturcova A, Apperley DC, Altaner C, Wess TJ, Jarvis MC: Microfibril diameter in celery collenchyma cellulose: X-ray scattering and NMR evidence. Cellulose. 2007, 14: 235-246. 10.1007/s10570-007-9116-1.View ArticleGoogle Scholar
- Ding SY, Himmel ME: The maize primary cell wall microfibril: a new model derived from direct visualization. J Agr Food Chem. 2006, 54: 597-606. 10.1021/jf051851z.View ArticleGoogle Scholar
- McCann MC, Bush M, Milioni D, Sado P, Stacey NJ, Catchpole G, Defernez M, Carpita NC, Hofte H, Ulvskov P, et al: Approaches to understanding the functional architecture of the plant cell wall. Phytochemistry. 2001, 57: 811-821. 10.1016/S0031-9422(01)00144-3.PubMedView ArticleGoogle Scholar
- Sarkar P, Bosneaga E, Auer M: Plant cell walls throughout evolution: towards a molecular understanding of their design principles. J Exp Bot. 2009, 60: 3615-3635. 10.1093/jxb/erp245.PubMedView ArticleGoogle Scholar
- Schermelleh L, Heintzmann R, Leonhardt H: A guide to super-resolution fluorescence microscopy. J Cell Biol. 2010, 190: 165-175. 10.1083/jcb.201002018.PubMedPubMed CentralView ArticleGoogle Scholar
- Abbe E: A contribution to the theory of the microscope and the nature of microscopic vision. Proc Bristol Nat Soc. 1874, 1: 200-261.Google Scholar
- Suslov D, Verbelen JP, Vissenberg K: Onion epidermis as a new model to study the control of growth anisotropy in higher plants. J Exp Bot. 2009, 60: 4175-4187. 10.1093/jxb/erp251.PubMedView ArticleGoogle Scholar
- Hoch HC, Galvani CD, Szarowski DH, Turner JN: Two new fluorescent dyes applicable for visualization of fungal cell walls. Mycologia. 2005, 97: 580-588. 10.3852/mycologia.97.3.580.PubMedView ArticleGoogle Scholar
- Oparka KJ, Read N: The use of fluorescent probes for studies on living plant cells. Plant Cell Biology: A Practical Approach. Edited by: Harris N, Oparka KJ. Oxford: Oxford University Press; 1994, 27-50.Google Scholar
- Heilemann M, van de Linde S, Schuttpelz M, Kasper R, Seefeldt B, Mukherjee A, Tinnefeld P, Sauer M: Subdiffraction-resolution fluorescence imaging with conventional fluorescent probes. Angew Chem Int Edit. 2008, 47: 6172-6176. 10.1002/anie.200802376.View ArticleGoogle Scholar
- Shim SH, Xia CL, Zhong GS, Babcock HP, Vaughan JC, Huang B, Wang X, Xu C, Bi GQ, Zhuang XW: Super-resolution fluorescence imaging of organelles in live cells with photoswitchable membrane probes. Proc Natl Acad Sci U S A. 2012, 109: 13978-13983. 10.1073/pnas.1201882109.PubMedPubMed CentralView ArticleGoogle Scholar
- Dempsey GT, Vaughan JC, Chen KH, Bates M, Zhuang X: Evaluation of fluorophores for optimal performance in localization-based super-resolution imaging. Nat Methods. 2011, 8: 1027-1036. 10.1038/nmeth.1768.PubMedPubMed CentralView ArticleGoogle Scholar
- Vogelsang J, Cordes T, Forthmann C, Steinhauer C, Tinnefeld P: Controlling the fluorescence of ordinary oxazine dyes for single-molecule switching and superresolution microscopy. Proc Natl Acad Sci USA. 2009, 106: 8107-8112. 10.1073/pnas.0811875106.PubMedPubMed CentralView ArticleGoogle Scholar
- Fitzgibbon J, Bell K, King E, Oparka K: Super-resolution imaging of plasmodesmata using three-dimensional structured illumination microscopy. Plant Phys. 2010, 153: 1453-1463. 10.1104/pp.110.157941.View ArticleGoogle Scholar
- Linnik O, Liesche J, Tilsner J, Oparka KJ: Unraveling the structure of viral replication complexes at super-resolution. Front Plant Sci. 2013, 4: 6.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffa A, Garin N, Sage D: Comparison of deconvolution software in 3D microscopy. A user point of view. GIT Imaging Microsc. 2010, 1: 41-43.Google Scholar
- de Monvel JB, Le Calvez S, Ulfendahl M: Image restoration for confocal microscopy: improving the limits of deconvolution, with application to the visualization of the mammalian hearing organ. Biophys J. 2001, 80: 2455-2470. 10.1016/S0006-3495(01)76214-5.View ArticleGoogle Scholar
- Chan J: Microtubule and cellulose microfibril orientation during plant cell and organ growth. J Microsc. 2012, 247: 23-32. 10.1111/j.1365-2818.2011.03585.x.PubMedView ArticleGoogle Scholar
- Roberts A, Oparka K: Plasmodesmal structure and development. Plasmodesmata. Edited by: Oparka K. Oxford: Blackwell Publishing, 1-31. Annual Plant Reviews; 2005.View ArticleGoogle Scholar
- Huang B, Jones SA, Brandenburg B, Zhuang XW: Whole-cell 3D STORM reveals interactions between cellular structures with nanometer-scale resolution. Nat Methods. 2008, 5: 1047-1052. 10.1038/nmeth.1274.PubMedPubMed CentralView ArticleGoogle Scholar
- EPFL Biomedical imaging group algorithms. http://bigwww.epfl.ch/algorithms.html.
- Born M, Wolf E: Principles of Optics 7th edition. Cambridge: Cambridge University Press; 2002.Google Scholar
- Richardson W: Bayesian-based iterative method of image restoration. J Opt Soc Am. 1972, 62: 55. 10.1364/JOSA.62.000055.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.