www.nature.com/scientificreports open A correlation analysis of Light Microscopy and X‑ray Microct imaging methods applied to archaeological plant remains’ morphological attributes visualization cristina Marilin calo1*, Marcia A. Rizzutto1, Sandra M. carmello‑Guerreiro2, carlos S. B. Dias3, Jennifer Watling4, Myrtle p. Shock5, carlos A. Zimpel6, Laura p. furquim4, francisco pugliese4,7 & eduardo G. neves4 in this work, several attributes of the internal morphology of drupaceous fruits found in the archaeological site Monte Castelo (Rondonia, Brazil) are analyzed by means of two different imaging methods. The aim is to explore similarities and differences in the visualization and analytical properties of the images obtained via High Resolution Light Microscopy and X‑ray micro‑computed tomography (X‑ray Microct) methods. Both provide data about the three‑layered pericarp (exo‑, meso- and endocarp) of the studied exemplars, defined by cell differentiation, vascularisation, cellular contents, presence of sclerenchyma cells and secretory cavities. However, it is possible to identify a series of differences between the information that can be obtained through each of the methods. These variations are related to the definition of contours and fine details of some characteristics, their spatial distribution, size attributes, optical properties and material preservation. the results obtained from both imaging methods are complementary, contributing to a more exhaustive morphological study of the plant remains. X‑ray Microct in phase‑contrast mode represents a suitable non‑ destructive analytic technique when sample preservation is required. The identification of the taxonomic provenance and preservation conditions of plant remains recovered from archaeological deposits is at the basis of archaeobotanical research. The whole procedure involves a detailed description of observable characteristics of the morphological phenotype of the sample, together with a comparative analysis of modern and ancient exemplars of related taxa. On the other hand, plant remains from archaeological sites are not always recognizable through their external morphological characteristics, for reasons that include the presence of similar traits between different species, morphological changes induced by taphonomic processes, and other issues which hinder the identification of specimens. In these cases, the analysis of the internal characteristics of the specimen can increase the amount of morphological data used to determine the anatomic and taxonomic provenance of ancient plant remains. In general, Archaeobotany has a wide palette of microscopy techniques that allow for various descriptions of both external and internal morphological features of archaeological plant samples1,2. In fact, it is possible to obtain 1 Laboratory of Archaeometry and Applied Sciences To Cultural Heritage (LACAPC), Institute of Physics, University of Sao Paulo, Sao Paulo, Brazil. 2Department of Plant Biology, Institute of Biology, University of Campinas (UNICAMP), Campinas, Brazil. 3Brazilian Synchrotron Light National Laboratory (LNLS), Brazilian Centre for Research in Energy and Materials (CNPEM), Campinas, Brazil. 4Museum of Archaeology and Ethnology (MAE), University of Sao Paulo, Sao Paulo, Brazil. 5Federal University of Western Para (UFOPA), Santarem, Brazil. 6Department of Archaeology, Federal University of Rondonia Foundation (UNIR), Porto Velho, Brazil. 7Institute of Geosciences, University of Brasilia, Brasília, Brazil. *email: marilincalo@usp.br Scientific RepoRtS | (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 1 Vol.:(0123456789) www.nature.com/scientificreports/ very detailed high resolution images of studied objects using some optical or electronic microscopy methods. Nevertheless, these images are limited to regions near the surface of the sample and/or to thin transparent or semitransparent samples. Internal structure and organization of thick and opaque objects are accessible only using a combination of microscopy, staining and microtome techniques, all of these having a necessary impact on sample integrity3,4. By contrast, no destructive sample preparation procedure is required for high resolution X-ray MicroCT. Volume reconstruction and virtual slicing prevent manipulation damage and the loss of material and information. As with staining, digital segmentation allows distinguishing and analyzing any feature separately, based on the image contrast. As a whole, the use of highly detailed 3D X-ray MicroCT images to study the internal structure of materials and objects is a growing and promising area in scientific research. Non-destructiveness is one of the most remarkable advantages of this method and its relevance increases when samples are rare, non-reproducible or very fragile. For this reason, X-ray MicroCT has been suggested as a suitable and useful procedure for applications on historical, artistic and culturally valuable objects made from diverse materials5. Several ethnographic and historical items manufactured on raw plant materials have been studied by X-Ray MicroCT i.e.6–8. Archaeological application of this technique on plant remains includes the analysis of contaminants in charred remains dated by radiocarbon9, dendrochronology10, wood characterization i.e.11,12 and identification of carpological remains13,14. However, so far, no studies have been conducted to assess the extent to which the X-ray MicroCT method can provide, by itself, the appropriate data for addressing issues of plant remains identification. In other words, if it is possible to obtain morphological information that is quantitatively and qualitatively comparable to that accessible by more regularly used light microscopy methods This study seeks to shed light on how three-dimensional data on the internal structure of plant remains contribute to the morphological examination of plant remains and what challenges arise from their application in Archaeobotany. This article is based on a previous X-ray MicroCT analysis of non-charred plant remains from the archaeological site Monte Castelo (Rondonia, Brazil). The results of that study were especially relevant to question the hypothesis of the caryopsis (maize-type) morphology, externally observed in several exemplars. In contrast, the technique allowed the identification of their drupaceous anatomy, and the suggestion of some morphological relationships with the Anacardiaceae15. Those same specimens were imaged using X-ray MicroCT and High Resolution Light Microscopy (HRLM) to specify the presence of a three-layered pericarp, vascular bundles and differentiated cellular contents. This essentially methodological contribution is aimed to explore variations in the visualization and analytical scope of both techniques, focusing on these morphological attributes. Methodology and materials the samples. This work compares the results obtained from the application of two alternative imaging techniques, HRLM and X-ray MicroCT, to the analysis of internal morphological characteristics of three plant remains samples. The specimens SMC_01, SMC_10 and SMC_11 were recovered at the archaeological site Monte Castelo, located in the floodplain of the Middle Guapore River, in the southwestern Brazilian Amazon. The site is an artificial mound with dimensions of ca. 120 m by 100 m and a height of 6.3 m, constructed in an area of ecological transition between periodically flooded savannahs and high-ground, evergreen tropical forest. Archaeological excavations reveal an uninterrupted sequence of occupation events showing evidence of several domestic activities dating from ca. 6,000 years BP16–19 (Fig. 1). The samples were selected from a large set of dried, non-charred, small drupe fruits (Fig. 2a). This botanical material comes from soil collections made after the definition of the cultural layers and, therefore, have secure provenience. These fruits were found concentrated in the most recent stratigraphic package of the shellmound that is related to the Bacabal Archaeological Phase, although they are also present in other archeological layers. This Phase consists predominately of layers of dark soil and shells, where richly decorated pottery, animal bones, plant materials, and human burials have been found15,18–21 (Fig. 2b). Examined with the naked eye and low magnification optics, the samples have globular to sub-globular shapes with pointy ends and rather irregular surfaces. The average maximum length and diameter of the studied exemplars are 4.9 mm and 4.6 mm respectively, while volume ranges between 33.3 to 40.4 mm3. Longitudinal lines or furrows are present, extending from the pointy end to approximately the medial region. Opposite to the pointy end, the rounded contour is only interrupted by a marked protuberance. imaging methods. The methodology is focused on the application of two different imaging methods X-ray Computed Microtomography (X-ray MicroCT) and High Resolution Light Microscopy (HRLM), whose performance for visualization and analysis are compared on the basis of a specific set of drupe morphological characteristics (see “Referenced internal morphological attributes of drupes” section). For the first imaging method, tomographic 3D images at a micrometer scale of samples SMC_01 and SMC_10 were reconstructed based on the X-ray MicroCT technique in propagated-based phase-contrast mode at the IMX Beamline of the Brazilian Synchrotron Light Laboratory (LNLS-CNPEM). The technical features of the IMX Beamline include the detector PCO pco.2000 area camera, 0.82 μm effective pixel size, 2048 × 2048 pixel, 14-bit CCD cooled camera. For both experiments an X-ray energy range (pink-beam) 4 to 14 keV was selected. The exposure time/ projection varied from 0.4 s in SMC_01 to 0.5 s in SMC_10. The sample-to-detector distance and the angular range were established at 260 mm and 360°, respectively. The effective pixel size was 0.82 μm for a total of 2048 projections on each microtomographied sample. A 350 μm Si filter was used in both cases. These parameters were set based on a previous test imaging experiment, performed on the same exemplar SMC_01 at the IMX Beamline, and another 10 specimens imaged using industrial microtomographers (Calo et al. 2019). They were Scientific RepoRtS | Vol:.(1234567890) (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 2 www.nature.com/scientificreports/ Figure 1. The location of the Monte Castelo archaeological site in the Middle Guapore River basin area (adapted)18 (QGis version 3.8, URL: https://www.qgis.org). selected in order to optimize resolution and contrast for the minimum necessary scanned area of the sample to reveal internal specific characteristics and the global time requirements to complete the experiment (delimited by the synchrotron experiments dynamic). X-ray MicroCT images result from the application of mathematical algorithms22 to X-ray experimental data, obtaining a sequence of 2D images (slices) and stacked together to reconstruct a 3D map of the object’s X-ray attenuation coefficient (transmission or absorption tomography). If certain experimental conditions are fulfilled, then phase shifts of the transmitted X-rays can also be detected. Absorption and phase modulation are effects that occur on X-rays crossing any kind of materials. The use of phase sensitive imaging techniques has the advantage of enhancing the visibility of low absorption materials as well as accentuating the edges between materials with small differences in refraction index. This requires a highly spatially-coherent X-ray beam provided mostly by synchrotron light sources5,23,24, but also by microfocus X-ray sources25. This study used the X-ray propagation based—phase contrast MicroCT via synchrotron radiation26,27. The in-house software Raft, developed at LNLS-CNPEM28–30, was applied to reconstruct data acquired from the IMX Beamline. The Paganin filter (sometimes considered a PCI projection restoration)31, regularly implemented in the IMX Beamline reconstruction software for edge enhancement of low contrast samples, were applied to projections acquired in this study. Visualization, segmentation and 3D morphological analysis of microtomographic images were achieved using CT Analyzer 1.15.4.0, developed by Bruker32 and the FIJI distribution version 1.52p of the image processing software ImageJ33–35 with the plug-ins MorphoLibJ36, 3DViewer37 and Volume Viewer38. All figures were edited using GNU Image Manipulation Program (GIMP) in the version 2.10.2 (https://www.gimp.org/). For the HRLM imaging method, 2D photomicrography of the internal morphological characteristics of the sample SMC_11 were examined and captured using two different microscopes. The first of them was an Olympus System Microscope model BX51 coupled to an Olympus Microscope Digital Camera Model DP71 in the Plant Anatomy Laboratory of the Biology Institute, University of Campinas. The second was an Erns Leitz GmbH microscope model Ortoplan-Pol coupled to a Kodak Digital Camera model DC4800 in the Laboratory of Complex Fluids of the Institute of Physics, University of Sao Paulo. As non-destructiveness (including damage or modification) is considered a relevant issue of the imaging method for plant remains, only one specimen was subjected to the HRLM method of sample preparation in the Plant Anatomy Laboratory of the Biology Institute, University of Campinas. The specimen SMC_11 was previously fixed in FAA (formaldehyde, acetic acid, 50% ethanol) for 24 h39. Then, it was dehydrated in an ethanol Scientific RepoRtS | (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 3 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 2. (a) Photomicrography of the sample SMC_10 (scale in mm) (Photo: Pedro Campos—Laboratory of Archaeometry and Applied Sciences to Cultural Heritage, Institute of Physics, University of Sao Paulo); (b) scheme of the north profile of the Monte Castelo shellmound and radiocarbon dates (cal. 2-sigma) (adapted)18. series and embedded in hydroxyethyl methacrylate resin (Historesin Leica)40. A total of 96 longitudinal slices of 8.0 µm thick were obtained from one of the halves of the lengthwise divided sample, using a Microm HM340E rotary microtome. The slices were stained with 0.05% Toluidine Blue in sodium acetate buffer pH 4.741 and fixed to microscope slides. All slices were mounted with water for examination. Referenced internal morphological attributes of drupes The internal characteristics considered here for the analysis of the fruit remains comprises the three-layered pericarp (endocarp, mesocarp and exocarp) defined by cell differentiation, vascularisation, presence of cellular contents, sclerenchyma cells, cavities and ducts. This set of descriptive variables is based on specific studies on fruit anatomy and classification42–44. In general, drupes can be distinguished from other fruits by their fleshy mesocarp and hard endocarp tissues. The peach could be considered a typical example of a drupe42,43. The exocarp (strictu sensu) is usually uniseriate and, in many fruits, represents the main protective layer of the pericarp. Notwithstanding, multiseriate exocarps are very seldom found in several fruits. This layer of tissue may be composed of regular epidermal cells of isodiametric or elongated shapes in surface view, but palisade cells could appear in some drupes. The outer tangential walls of the epidermal cells are generally thick and covered by a well-developed cuticle. Stomata, lenticels and trichomes are frequently present. Anthocyanins are regularly reported as contents in the cell sap of epidermal cells, producing the red, purple, blue and black colors of the fruit. Calcium oxalate crystals, raphides, siliceous impregnations, cutines, suberines and waxes also occur in specialized epidermis cells42,45. The endocarp (strictu sensu) is often composed of smaller cells than the exocarp, especially in fleshy fruits, and they are frequently elongated. Drupes develop a distinctive hard multilayered endocarp formed by sclerenchymatic tissue. Differently from parenchymatous endocarps, its main function is to serve as supporting tissue and to protect the seed from injuries. The sclerenchyma is formed by thick secondary walls and strongly lignified cells denominated sclereids (short cells), fibers (long cells) and fiber-sclereids. Several types of sclereids could be present in the endocarp structure. The most common are brachysclereids, or stone cells, roughly isodiametric or somewhat elongated cells; macrosclereids, represented by elongated columnar cells; and osteosclereids or bone-cells, also columnar but with enlarged ends42,45. As crystals are often found in the vicinity of sclerenchymatic cells, a crystal layer frequently accompanies the endocarp in drupes42. Inorganic deposits on plants consist mostly of calcium salts and anhydrides of silica which occurs in many plant families. Although calcium carbonate crystals are rarely present, the crystalline form of calcium oxalate is very common. They appear as prismatic, raphide, druse and styloid crystal morphologies. In some tissues, calcium oxalate crystals arise in cells that resemble adjacent, crystal-free cells. In others, the crystals Scientific RepoRtS | Vol:.(1234567890) (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 4 www.nature.com/scientificreports/ Figure 3. General views of the three-layered structure of the fruit as seen by HRLM image of a longitudinal microtome section from SMC_11 (a); virtual X-ray MicroCT transversal section from SMC_01 (b) and SMC_10 (c) idem for longitudinal section (d, e) (ImageJ-FIJI). are formed in idioblasts, specialized secretory cells with a distinctive form (vesicles, sacs, or tubes). On the other hand, silica compounds are more present in cell walls or form silica bodies or phytoliths in the lumen of cells45. Besides crystals, idioblasts can also produce oils, mucilage and tannins42 and can be distributed all over the pericarp. Other secretory structures in drupes are the cavities and ducts, formed by several specialized cells which release those secondary metabolites to the intercellular space. This characteristic is very conspicuous in the mesocarp of drupes of the Anacardiceae family i.e.46–52, also referred to as lacunar mesocarp44. The parenchyma plays an essential part in the formation of the fleshy mesocarp of drupes. It generally consists of large rounded or elongated (in a radial direction, generally the inner mesocarp) cells with thin walls and large vacuole. In ripe drupes, these are usually distinguished by their juiciness. The cell sap of parenchyma may be enriched with anthocyanin and sugar, together with tannins, acids and other contents. Crystals of different shapes may also occur and appear in specially preformed layers. Small intercellular spaces may eventually develop between the cells. As the parenchyma represents the edible part of drupes, it functions as a storage tissue in maturing fruits; it may contain starch, sugar, oil, fats, etheric oil (deposited in idioblasts, special glands or oil ducts) and latex in laticifers. Collenchymatic tissue is also found frequently in the fleshy mesocarp of drupes. Scientific RepoRtS | (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 5 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 4. Detail of cell morphologies of the exocarp and endocarp layers: (a, b) HRLM from sample SMC_11; X-ray MicroCT transversal section; (c, d) from sample SMC_01; and (e, f) from sample SMC_10 (ImageJ-FIJI). This small-celled and thick walled cells tissue occupies the periphery of the pericarp or the vicinity of vascular bundles. Its primary function is mechanical support42,43. Vascular bundles are the main transport channels of water and nutrients in the fruit. They form an anastomosed network in the drupe’s pericarp running through the mesocarp. The vascular supply of fruit varies according to the type of ovary from which it originates. The fleshy fruits develop an additional bundle system and the whole finer bundles are more numerous in fleshy tissues than in dry and/or stony parts45. Scientific RepoRtS | Vol:.(1234567890) (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 6 www.nature.com/scientificreports/ Figure 5. Volume reconstruction of (a) seed-embryo; (b) endocarp; (c) mesocarp; (d) exocarp; (e) complete composed volume; (f) longitudinally sectioned, from sample SMC_10 (CT Analyzer; ImageJ-FIJI, plug-in: 3D Viewer). Results pericarp. Both HRLM and X-ray MicroCT longitudinal sections show a three-layered internal structure marked by straight differences between cellular components (Fig. 3). In HLRM images, the exocarp corresponds to a multiseriate layer stained in brown-green. Cells are small and slightly regular with prevailing isodiametric shapes and some contents. No palisade cells were observed. A fine cuticle layer is present on the outer cells in some regions of the exocarp. The mesocarp appears as the most voluminous tissue, formed by large, irregular parenchymatic cells. Primary cell walls are well preserved (Fig. 4a). The endocarp is a multiseriate layer formed by mostly regular, small and highly lignified sclerenchymatic cells (brachysclereids) stained in green–blue (Fig. 4b). No structure related to seed-embryo tissues was observed in HRLM images from the studied exemplar. In X-ray MicroCT images, virtual longitudinal and transversal sections allow one to distinguish the threelayered pericarp in terms of the differences between the morphologies of their composing cells (Fig. 3b,c), very similar to those observed in HRLM images. In sample SMC_01, the limits between the layers are barely defined, resulting in a gradate transition from one cell morphology to another. Any relevant contrast level differences based on variations of the density values of layers are evident (Fig. 4c–f). This character makes the segmentation of images less evident, even when non-algorithmic procedures are used. A membranous tissue remnant of the seed-embryo structure is clearly observable inside the central cavity. It appears weakly attached to the inner side of the endocarp (Fig. 5) (Supplementary Video S1 online). presence of cellular contents. Most of the exocarp cells show the presence of brown-reddish stained contents in the HRLM images. They could be associated with phenolic contents and/or pigments, commonly present in the skin of fruits (Fig. 4a). However, this characteristic is not observed in the X-ray MicroCT images, at least for the experimental parameters set for this study. HRLM also shows crystal contents in the cells of the endocarp. Some of these cells contain calcium oxalate crystals which exhibit the typical extinction cross under polarized light (Fig. 6a,b). These particles were also evidenced in X-ray MicroCT images, which exhibit several white and compacted particles in the inner layer of some cells. Their high X-ray attenuation coefficient makes them easily distinguishable from tissues and other particles in the samples because they are comparatively denser and their composition is very different from surrounding material (Fig. 6c–e). A total of 19 particles were segmented from a random section of the endocarp of sample SMC_01, most of them measuring between 1.50E+07 μm3 and 6.40 E+07 μm3 (Fig. 7a,b). The Scientific RepoRtS | (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 7 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 6. (a) Longitudinal section HRLM image from the inner layer region showing the thick walled cells, some of which contains crystals; (b) the same image evidencing the crystalline character of the particles under polarized light; (c) X-ray MicroCT reconstructed volume from a section of the endocarp of the sample SMC_01; the lines in blue, green and red indicate the position of the slices at the XY, XZ and YZ sections, respectively and their images appear on the left side; (d, e) the slices appear repositioned in order to reveal the localization of the particles inside the thick walled cells; the yellow lines in the volume indicate the clipping plane entering/cutting the volume; (f) reconstructed volumes from particles 1 and 2 (scale in μm) (CT Analyzer; ImageJ-FIJI, plug-ins: 3D Viewer and Volume Viewer). reconstructed volumes of some of them suggest crystal-like multi-faceted structures defined by edges and angles (Fig. 6f). Notwithstanding, the resolution of the X-ray microtomographies in this study is not high enough to confirm this morphology. X-ray MicroCT images also show a second type of dense particle present in the mesocarp and inside the cavity, next to the open pointy region of the samples. They do not seem to be directly associated with intracellular contents and their 3D morphological analysis shows a total of 111 particles which are bigger than the crystal ones. The volume of most of these particles is between 6.6 E+05 μm3 and 5.4 E+09 μm3, while others can reach 1.2 E+11 μm3 (Fig. 7c,d). Their Sauter diameter (the diameter of the sphere with the same volume/area ratio of a given particle) ranges between 5.6 E+05 μm3and 2.9 E+06 μm3. They moreover lack the multi-faceted morphology and might be better described as amorphous concretion-like structures (Fig. 8). The complete list of volume Scientific RepoRtS | Vol:.(1234567890) (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 8 www.nature.com/scientificreports/ Figure 7. Volume and Sauter diameter distribution of particles in the endocarp (a, b) and mesocarp (c, d) of the sample SMC_01 (CT Analyzer). and Sauter diameter values, as well as other measured variables (‘Surface’, ‘Volume-equivalent sphere diameter’ and ‘Surface-equivalent sphere diameter’) can be accessed in the Supplementary Table S2 online. Vascularisation. Several segments of the vascular system of the fruit can be observed in HRLM slices as vascular elements or fragments of bundles showing spiral thickened cell walls (Fig. 9a–c). They are present abundantly in the fleshy mesocarp layer, and are also morphologically distinguishable in X-ray MicroCT longitudinal virtual slices (Fig. 9d,e). The volumetric reconstruction of the more conspicuous bundles through the mesocarp (at least 12) reveals their distribution in the region near to the endocarp (Fig. 9f–i). cavities and ducts. Some microtomed longitudinal sections show several cavities with sub-circular contours which is compatible to those previously described for several species of the Anacardiaceae. The diameter of these structures ranges from 100 µm to 250 µm and some of them still obtrude with contents (Fig. 10a–d). They were observed less frequently in the distal half of the sample. Nonetheless, the microtomographied exemplars showed some of them in this position in both longitudinal and transversal virtual sections of the mesocarp (Fig. 10e–j). Contents are not evident in 2D virtual slices. Their size ranges from 80 µm to 100 µm in transversal and longitudinal diameter, according to their 3D sub-spherical/drop morphology (Figs. 10k and 11a). Another type of structure with contents was also observed in the Micro-CT images, appearing as open spaces along the mesocarp, where some substances persist, probably due to their secretory function (Fig. 11b). conclusions The high resolution images obtained by both techniques were able to evidence each of the studied morphological traits in the archaeobotanical samples from Monte Castelo. HRLM and X-ray MicroCT are able to recognize the three differentiated layers of tissues composing the pericarp of the fruit and some of its most relevant histological characteristics: the presence of vascular tissues in the mesocarp; the existence of mineral particles within and between the tissues; and the occurrence of secretory structures. However, it is possible to point out some differences between the information obtained through each method. Scientific RepoRtS | (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 9 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 8. 2D Visualization of dense amorphous particles in a virtual transversal section from sample SMC_01 (a) Volume reconstruction of amorphous particles in SMC_01 placed in the orthogonal section of the whole sample volume; (b) Volume reconstruction of the same sample showing the localization of three amorphous particles and their respective 3D image; (c–e) (CT Analyzer; ImageJ-FIJI, plug-ins: 3D Viewer and Volume Viewer). The clear contours and detailed definition of the objects in the HRLM images are among their most notable characteristics. The contrast between the different colors and color shades induced by staining reinforces these effects. HRLM allowed the visualization of very fine structures in vascular tissues, cell walls and contents, which were either barely or not visible in the X-ray MicroCT images. Scientific RepoRtS | Vol:.(1234567890) (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 10 www.nature.com/scientificreports/ Figure 9. Vascular elements in microtomed longitudinal sections of SMC_11 (a–c); Vascular bundles as seen in X-ray MicroCT virtual longitudinal sections from samples SMC_01 (d) and SMC_10 (e); Volume reconstruction of a part of the vascularisation system in the region of the mesocarp-endocarp edge and their emplacement in exemplars SMC_01 (f, g) and SMC_10 (h, i) (CT Analyzer; ImageJ-FIJI, plug-ins: MorphoLibJ and 3D Viewer). Similarly, the use of polarized light coupled to HRLM enhances the analysis of mineral particles, giving a clear idea about their crystalline nature. This was not possible using microtomography. Despite this, the reconstructed volumes of the particles show differences between faceted and non-faceted morphologies, suggesting that a morphological analysis of crystalline minerals could be achieved focusing on appropriate X-ray MicroCT experimental settings. Indeed, the possibility of performing 3D reconstructions of samples and visualizing their internal structure in a three-dimensional space allows the user to reveal the position and distribution of vascular bundles, mineral particles and secretory structures. Moreover, dynamic and non-destructive virtual slicing addresses the full range of slicing planes describing characteristics in a way that is not available to HRLM static slicing. The non-destructive character of X-ray MicroCT helped to describe the presence of remnant tissues inside the central cavity of the samples which are absent in the microscopy images and probably represent the seed-embryo structures. As shown by microtomography, these structures are weakly tied to the rest of the fruit and could Scientific RepoRtS | (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 11 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 10. Sub-circular cavities observed in the microtomed longitudinal section from SMC_11 (a–d); Virtual X-ray MicroCT sections of cavities from sample SMC_01, suggesting their sub-circular/drop morphology (e–j); Volume reconstruction of one of these X-ray MicroCT imaged cavities (k) (CT Analyzer; ImageJ-FIJI, plug-in: 3D Viewer). be detached and lost from the microtomed-for-microscopy sample. Likewise, some differences in the average size of the secretory structures were observed when measured by X-ray MicroCT or HRLM images. This effect could be linked to distortions originating in sample preparation procedures, including sample dehydration and resin embedding, meaning all quantitative comparisons between specimens should be made under the same imaging method53. Results obtained from both imaging methods are complementary in terms of a more exhaustive analysis of morphological characteristics of plant remains. However, X-ray MicroCT represents a suitable alternative analytic technique when sample preservation is required, for example when it is a rare, fragile specimen, or it will later be intended for further analysis. Based on the results of this study, the application of this analytical technique presents two types of challenge. The first is of a technical nature, and has to do with the limitations of the technique in defining contours and subtle definition of details when applied to biological materials. Beyond the loss of information that this may represent for the study of plant remains, its resolution depends fundamentally on the technological development Scientific RepoRtS | Vol:.(1234567890) (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 12 www.nature.com/scientificreports/ Figure 11. Volume reconstruction placing the cavities (in red) in the pericarp of sample SMC_01 (grey) with the endocarp highlighted in blue (a); Section on the volume reconstruction of the sample SMC_10 showing contents in a mesocarp secretory structure (b) (CTAnalyzer; ImageJ-FIJI, plug-ins: 3D Viewer and Volume Viewer). of X-ray MicroCT techniques, new instruments and even more advanced X-ray microscopy modalities such as ptychography. The second challenge is a methodological one, linked to some basic procedures of archaeobotanical studies. A good part of Archaeobotany research is focused on taxonomic identification of archaeological plant remains. This task is based primarily on morphological comparative analysis between old specimens and contemporary plant reference material. These references may be previous bibliographical studies, modern described plant exemplars itself or, as is generally the case, a combination of both. However, the analytical parameters used by morphology and plant anatomy to describe different plant taxa have been constructed based on two-dimensional, light or SEM microscopy images. This means that a large number of the variables used to describe and compare plant material, relies on images of a different nature than three-dimensional X-ray MicroCT. It is therefore necessary to consider providing three-dimensional descriptions of reference plant material to archaeobotanical studies using this technique. Nonetheless, it is worth mentioning that virtual slicing would offer two-dimensional images equivalent to Light Microscopy, but properly related to volume, which could pave the way for new taxonomical standards. The present study has shed light on some differences observed between analytic images based on the variations of the light reflected by the object (HRLM) from those based on the variations of the X-ray attenuation coefficient (and phase shifts) of the irradiated object (X-ray phase-contrast MicroCT). A more detailed analysis of these differences, their impact on the recognition and comparability of specific morphological characteristics, and the production of three-dimensional reference material, will be valuable in assisting the implementation of X-ray MicroCT to the study of archaeological plant remains. This research benefitted from resources at the Brazilian Synchrotron Light Laboratory (LNLS), an open national facility operated by the Brazilian Centre for Research in Energy and Materials (CNPEM) for the Brazilian Ministry for Science, Technology, Innovations and Communications (MCTIC). The IMX Beamline staff is acknowledged for the assistance during the experiments. We are also thankful to the Laboratory of Complex Fluids of the Institute of Physics (IF), University of Sao Paulo (USP), for microscopy sessions, especially to Dennys Reis for technical support with microscopy images. Received: 9 April 2020; Accepted: 30 June 2020 References 1. Jones, M. K. & Colledge, S. M. Archaeobotany and the transition to agriculture. In Handbook of Archaeological Sciences (eds Brothwell, D. R. & Pollard, A. M.) 393–401 (Wiley, New York, 2001). 2. Weiss, E. & Kislev, M. E. Plant remains as a tool for reconstruction of the past environment, economy, and society: Archaeobotany in Israel. Isr. J. Earth Sci. 56, 163–173 (2007). 3. Landis, E. N. & Keane, D. T. X-ray microtomography. Mater. Charact. 61, 1305–1316 (2010). 4. Stuppy, W. H., Maisano, J. A., Colbert, M. W., Rudall, P. J. & Rowe, T. B. Three-dimensional analysis of plant structure using highresolution X-ray computed tomography. Trends Plant Sci. 8, 2–6 (2003). Scientific RepoRtS | (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 13 Vol.:(0123456789) www.nature.com/scientificreports/ 5. Mancini, L., Tromba, G. & Zanini, F. Structural microanalysis with synchrotron radiation: archaeometric applications at Elettra. J. Neutron Res. 14, 75–79 (2006). 6. Fioravanti, M. et al. Non-invasive wood identification of historical musical bows. IAWA J. 38, 285–296 (2017). 7. Haneca, K. et al. X-ray sub-micron tomography as a tool for the study of archaeological wood preserved through the corrosion of metal objects. Archaeometry 54, 893–905 (2012). 8. McNoughton, A., Smith, C. A. & Lowe, B. J. Micro-computed tomography for plant identification in artefacts. in Proceedings of the ICOM-CC 17th Triennial Conference, 1606. (2014). 9. Bird, M. I., Ascough, P. L., Young, I. M., Wood, C. V. & Scott, A. C. X-ray microtomographic imaging of charcoal. J. Archaeol. Sci. 35, 2698–2706 (2008). 10. Okochi, T., Hoshino, Y., Fujii, H. & Mitsutani, T. Nondestructive tree-ring measurements for Japanese oak and Japanese beech using micro-focus X-ray computed tomography. Dendrochronologia 24, 155–164 (2007). 11. Dreossi, D. et al. Synchrotron radiation Micro-tomography: a non-invasive tool for th characterization of archaeological wood. In Wood Science for Conservation of Cultural Heritage (ed. Uzielli, L.) (Firenze University Press, Firenze, 2010). 12. Mizuno, S., Torizu, R. & Sugiyama, J. Wood identification of a wooden mask using synchrotron X-ray microtomography. J. Archaeol. Sci. 37, 2842–2845 (2010). 13. Coubray, S., Zech-Matterne, V. & Mazurier, A. The earliest remains of a Citrus fruit from a western Mediterranean archaeological context? A microtomographic-based re-assessment. C. R. Palevol 9, 277–282 (2010). 14. Murphy, C. & Fuller, D. Q. Seed coat thinning during horsegram (Macrotyloma uniflorum) domestication documented through synchrotron tomography of archaeological seeds. Sci. Rep. 7, 5369 (2017). 15. Calo, C. M. et al. Study of plant remains from a fluvial shellmound (Monte Castelo, RO, Brazil) using the X-ray MicroCT imaging technique. J. Archaeol. Sci. Rep. 26, 101902 (2019). 16. Miller, E. T. Algumas Culturas Ceramistas, do Noroeste do Pantanal do Guaporé à Encosta e Altiplano Sudoeste do Chapadão dos Parecis. Origem, Difusão/Migração e Adaptação – do Noroeste da América do Sul ao Brasil. Rev. Bras. Linguística Antropológica 5, 335–383 (2013). 17. Miller, E. Th. Pesquisas Arqueológicas no Pantanal do Guaporé-RO, Brasil: A Seqüência Seriada da Cerâmica da Fase Bacabal. in Arqueologia Interpretativa O Método Quantitativo para Estabelecimento de Sequências Cerâmicas: Estudos de Caso (ed. Meggers, B.) (UNITINS, 2009). 18. Pugliese, F. A. A história indígena profunda do sambaqui Monte Castelo: um ensaio sobre a longa duração da cerâmica e das paisagens no sudoeste amazônico. Doctoral Thesis, Museu de Arqueologia e Etnologia, Universidade de São Paulo (2018). https ://doi.org/10.11606/T.71.2019.tde-02102018-162628. 19. Pugliese, F. A., Zimpel, C. A. & Neves, E. G. What do Amazonian Shellmounds Tell Us About the Long-Term Indigenous History of South America? in Encyclopedia of Global Archaeology 1–25 (Springer International Publishing, 2019). https://doi.org/10.1007/9783-319-51726-1_3030-1. 20. Furquim, L. P. Arqueobotânica e Mudanças Socioeconômicas durante o Holoceno Médio no Sudoeste da Amazônia. Master Dissertation, Museu de Arquelogia e Etnologia, Universidade de São Paulo (2018). https://doi.org/10.11606/D.71.2019.tde-30112 018-102517. 21. Zimpel, C. A. A fase Bacabal e seus correlatos arqueológicos no sudoeste da Amazônia. Doctoral Thesis, Museu de Arqueologia e Etnologia, Universidade de São Paulo (2018). https://doi.org/10.11606/T.71.2019.tde-29112018-144329. 22. Stock, S. R. Microcomputed Tomography: Methodology and Applications (CRC Press, Boca Raton, 2008). 23. Baruchel, J. et al. Phase imaging using highly coherent X-rays: radiography, tomography, diffraction topography. J. Synchrotron Radiat. 7, 196–201 (2000). 24. Snigirev, A., Snigireva, I., Kohn, V., Kuznetsov, S. & Schelokov, I. On the possibilities of x-ray phase contrast microimaging by coherent high-energy synchrotron radiation. Rev. Sci. Instrum. 66, 5486–5492 (1995). 25. Wilkins, S. W., Gureyev, T. E., Gao, D., Pogany, A. & Stevenson, A. W. Phase-contrast imaging using polychromatic hard X-rays. Nature 384, 335 (1996). 26. Cloetens, P. et al. Observation of microstructure and damage in materials by phase sensitive radiography and tomography. J. Appl. Phys. 81, 5878–5886 (1997). 27. Cloetens, P., Barrett, R., Baruchel, J., Guigay, J.-P. & Schlenker, M. Phase objects in synchrotron radiation hard x-ray imaging. J. Phys. Appl. Phys. 29, 133–146 (1996). 28. Koshev, N., Helou, E. S. & Miqueles, E. X. Fast Backprojection Techniques for High Resolution Tomography (2016). 29. Miqueles, E. X., Archilha, N. L., Anjos, M. R. D., Westfahl Jr., H. & Helou, E. S. Phase-Retrieval as a Regularization Problem. ArXiv170205092 Math (2017). 30. Miqueles, E. X., Rinkel, J., O’Dowd, F. & Bermúdez, J. S. V. Generalized Titarenko’s algorithm for ring artefacts reduction. J. Synchrotron. Radiat. 21, 1333–1346 (2014). 31. Paganin, D., Mayo, S. C., Gureyev, T. E., Miller, P. R. & Wilkins, S. W. Simultaneous phase and amplitude extraction from a single defocused image of a homogeneous object. J. Microsc. 206, 33–40 (2002). 32. Bruker. Skyscan softwares updates. Bruker microCT | Software Downloads https://bruker-microct.com/products/downloads.htm (2018). 33. Rueden, C. T. et al. Image J2: ImageJ for the next generation of scientific image data. BMC Bioinformatics 18, 529 (2017). 34. Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012). 35. Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012). 36. Legland, D., Arganda-Carreras, I. & Andrey, P. MorphoLibJ: integrated library and plugins for mathematical morphology with ImageJ. Bioinformatics 32, 3532–3534 (2016). 37. Schmid, B., Schindelin, J., Cardona, A., Longair, M. & Heisenberg, M. A high-level 3D visualization API for Java and ImageJ. BMC Bioinformatics 11, 274 (2010). 38. Barthel, K. U. Volume Viewer. ImageJ - Image Processing and Analysis in Java https://imagej.nih.gov/ij/plugins/volume-viewer.html (2005). 39. Johansen, D. A. Plant Microtechnique (McGraw-Hill, New York, 1940). 40. Gerrits, P. O. & Smid, L. A new, less toxic polymerization system for the embedding of soft tissues in glycol methacrylate and subsequent preparing of serial sections. J. Microsc. 132, 81–85 (1983). 41. O’Brien, T. P., Feder, N. & McCully, M. E. Polychromatic staining of plant cell walls by toluidine blue O. Protoplasma 59, 368–373 (1964). 42. Roth, I. Fruits of Angiosperms (Gebr, Borntraeger, 1977). 43. Spjut, R. W. A systematic treatment of fruit types. Mem NY Bot Gard USA 70, 1–182 (1994). 44. Barroso, G. M., Morim, M. P., Peixoto, A. L. & Ichaso, C. L. F. Frutos e Sementes. Morfología aplicada à sistemática das Dicotiledôneas. (UFV - Universidade Federal de Viçosa, 1999). 45. Evert, R. F. Esau’s Plant Anatomy: Meristems, Cells, and Tissues of the Plant Body: Their Structure, Function, and Development (Wiley, New York, 2006). 46. Carmello-Guerreiro, S. M. & Paoli, A. A. S. Anatomy of the pericarp and seed-coat of Lithraea molleoides (Vell.) Engl. (Anacardiaceae) with taxonomic notes. Braz. Arch. Biol. Technol. 48, 599–610 (2005). 47. Herrera, F. et al. Fruit morphology and anatomy of the spondioid anacardiaceae. Bot. Rev. 84, 315–393 (2018). Scientific RepoRtS | Vol:.(1234567890) (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 14 www.nature.com/scientificreports/ 48. Lacchia, A. P. S. & Guerreiro, S. M. C. Ultra-structural aspects of secretory canals in vegetative and reproductive organs of Anacardiaceae. Acta Bot. Bras. 23, 376–388 (2009). 49. Li, X., Baskin, J. M. & Baskin, C. C. Pericarp ontogeny and anatomy in Rhus aromatica Ait. and R. glabra L. (Anacardiaceae). J. Torrey Bot. Soc. 126, 279 (1999). 50. Machado, S. R. & Carmello-Guerreiro, S. M. Estrutura e desenvolvimento de canais secretores em frutos de Schinus terebinthifolius Raddi (Anacardiaceae). Acta Bot Bras 7 (2001). 51. Tölke, E. E. A. D. et al. Pericarp ontogeny of Tapirira guianensis Aubl. (Anacardiaceae) reveals a secretory endocarp in young stage. Acta Bot. Bras. 31, 319–329 (2017). 52. Wannan, B. S. & Quinn, C. J. Pericarp structure and generic affinities in the Anacardiaceae. Bot. J. Linn. Soc. 102, 225–252 (1990). 53. Dias, C. S. B. et al. Comparative analysis of sample preparation protocols of soft biological tissues for morphometric studies using synchrotron-based X-ray microtomography. J. Synchrotron Radiat. 26, 1–11 (2019). Acknowledgements This research benefitted from resources at the Brazilian Synchrotron Light Laboratory (LNLS), an open national facility operated by the Brazilian Centre for Research in Energy and Materials (CNPEM) for the Brazilian Ministry for Science, Technology, Innovations and Communications (MCTIC). The IMX Beamline staff is acknowledged for the assistance during the experiments. We are also thankful to the Laboratory of Complex Fluids of the Institute of Physics (IF), University of Sao Paulo (USP), for microscopy sessions, especially to Dennys Reis for technical support with microscopy images. Author contributions C.M.C. performed the X-Ray MicroCT experiments and the H.R.L.M. slices examination, as well as the processing and analysis of the X-ray MicroCT data. She wrote the main manuscript text and prepared the figures. M.A.R. and C.S.B.D. contributed to the design and conduct of both x-ray MicroCT experiments and helped with image processing. S.M.C.G. conducted the sample preparation procedures for H.R.L.M. and intervened in their morphological analysis and microscope images production. L.P.F., J.W., M.P.S. and C.A.Z. were involved in the excavation of Monte Castelo, the recovery of the samples and the discussion of their archaeological significance. L.P.F. also performed the preliminary analysis of their external macroscopic morphology. J.W. revised the final english version of the main text. F.A.P. and E.G.N. coordinate the archaeological site excavation and the study of the samples’ provenance contexts and significance. All authors have contributed to the discussion and revision of the manuscript. funding Resources for this study was provided by the São Paulo Research Foundation (FAPESP) Fellowships 201612867-7 (C. M. Calo) and 2017-25157-0 (J. Watling), FAPESP Research Grant 2017-11817-9 (E. G. Neves) and the National Council for Scientific and Technological Development (CNPq) Research Project 307179/2013-3 (E. G. Neves). competing interests The authors declare no competing interests. Additional information Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-71726-z. Correspondence and requests for materials should be addressed to C.M.C. Reprints and permissions information is available at www.nature.com/reprints. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2020 Scientific RepoRtS | (2020) 10:15105 | https://doi.org/10.1038/s41598-020-71726-z 15 Vol.:(0123456789)