Unveiling microbial preservation under hyperacidic and oxidizing conditions in the Oligocene Rio Tinto deposit

1 Vol.:(0123456789) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports Unveiling microbial preservation under hyperacidic and oxidizing conditions in the Oligocene Rio Tinto deposit David C. Fernández‑Remolar1,9,10*, Daniel Carrizo2, Mourad Harir3, Ting Huang4, Ricardo Amils2,5, Philippe Schmitt‑Kopplin3,6, Laura Sánchez‑García2, David Gomez‑Ortiz7 & Per Malmberg8 The preservation of biosignatures on Mars is largely associated with extensive deposits of clays formed under mild early Noachian conditions (> 3.9 Ga). They were followed by widespread precipitation of acidic sulfates considered adverse for biomolecule preservation. In this paper, an exhaustive mass spectrometry investigation of ferric subsurface materials in the Rio Tinto gossan deposit (~ 25 Ma) provides evidence of well‑preserved molecular biosignatures under oxidative and acidic conditions. Time of flight secondary ion mass spectrometry (ToF–SIMS) analysis shows a direct association between physical‑templating biological structures and molecular biosignatures. This relation implies that the quality of molecular preservation is exceptional and provides information on microbial life formerly operating in the shallow regions of the Rio Tinto subsurface. Consequently, low‑pH oxidative environments on Mars could also record molecular information about ancient life in the same way as the Noachian clay‑rich deposits. The surface of Mars is widely acknowledged to be inhospitable for life because of the direct and indirect damaging effects of solar radiation1. Although there are a few bacterial forms that could potentially thrive under such extreme conditions2, their occurrence would be so low that the chances of detecting them would be negligible. Given these conditions, the subsurface areas of Mars have been claimed to host a potential biosphere as a consequence of its planetary evolution3,4. While the subsurface areas are considered a primary target for searching for extant or extinct life on Mars5, this is not the case for surface or subsurface materials affected or derived from the hydrochemical activity of the strongly oxidizing and acidic solutions, which mostly occurred during the Hesperian (< 3.7Ga)6. Although different authors report neutral to alkaline conditions for underground Mars7,8, the upper part of the planet interior has certainly been affected by the circulation and downwelling of low-pH and oxidizing solutions9–11. Under such conditions, interactions between the oxidizing and acidic aggressive solutions and the Martian crust should have enhanced mineralization from ion dissolution and the release of byproducts at the surface, such as different Al-bearing clays whose composition depends on the weathering degree12. Among other salts, different Fe-bearing compounds, such as sulfates and oxides, could subsequently precipitate underground13 by oversaturation, which would lead to mineralization by ferruginous materials infilling the substrate faults and cracks or cementing the porosity. During mineralization, different microbial forms thriving in the subsurface could eventually be entombed inside a saline mineral matrix by simple oversaturation and precipitation of the incoming solutions. The mineralization process could also be promoted by the microbial OPEN 1CEA, CNRS, IBS, Metalloproteins Unit, Université Grenoble Alpes, 38000 Grenoble, France. 2Centro de Astrobiología (INTA-CSIC), Madrid, Spain. 3Research Unit Analytical Biogeochemistry, Helmholtz Zentrum München, Neuherberg, Germany. 4State Key Laboratory of Lunar and Planetary Sciences, Macau University of Science and Technology, Macau, China. 5Centro de Biología Molecular Severo Ochoa, Universidad Autónoma de Madrid, Madrid, Spain. 6Chair of Analytical Food Chemistry, Technical University Munich, 85354 Freising-Weihenstephan, Germany. 7ESCET-Área de Geología, Universidad Rey Juan Carlos, 28933 Móstoles, Madrid, Spain. 8Chemistry and Chemical Engineering, Chalmers University of Technology, Kemivägen 10, 412 96 Gothenburg, Sweden. 9Present address: State Key Laboratory of Lunar and Planetary Sciences, Macau University of Science and Technology, Macau 999078, PR China. 10Present address: CNSA Macau Center for Space Exploration and Science, Macau 999078, PR China. *email: [email protected] 2 Vol:.(1234567890) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ activity itself through direct or indirect mechanisms that can trigger or even increase mineral precipitation on the cell surface14,15. Iron oxides have been claimed to preserve molecular biosignatures with implications for detecting traces of life on Mars16. However, acidic and oxidizing environments have currently been ruled out as ideal conditions for the preservation of molecular biosignatures. This is also the case for Martian materials that have resulted from alteration or precipitation due to the activity of such acidic and oxidizing fluids17,18. The reason lies in the fact that these solutions inherit a high oxidizing capacity that can be released through different biotic or abiotic pathways. If life arose on Mars, it could have been a very active agent for recycling inorganic and organic compounds by using inorganic electron acceptors (e.g., Fe3+ and SO42−) that are abundant in acidic fluids19. Although the solar X-ray and UV irradiation rates occurring in the Hesperian are unknown, different oxidants that would surely have been produced by incident radiation on the atmosphere and hydrosphere7,20 would have oxidized biomolecules and organic compounds without the participation of microbes. As the preservation potential of molecular biosignatures under acidic conditions is considered very low, the planetary community has ruled out the list of potential targets to search for molecular evidence of extinct life on Mars17. Therefore, deposits formed under mild conditions (neutral pH, low oxidizing to reducing conditions) have become first-class targets to find traces of living organisms on early Mars21. On Earth, records of ancient life activity can occur in the form of physical structures replicated by mineral precipitation22. Indeed, as the final diagenetic byproducts of acidic deposits are goethite (FeOOH) and hematite (Fe2O3), the preservation potential of mineralized structures is extremely high. The diagenetic processes affecting acidic materials consist of a combination of loss of sulfate ions and dehydration under oxic conditions that are controlled by the inflow of meteoric solutions after exhumation of materials from the initial acidic conditions23. Very recently, a diverse association of lipids has been detected in modern to recent materials formed under acidic and oxidizing conditions24, which suggests that organics can be preserved in the first stage of diagenesis. Indeed, the preservation potential of acidic deposits for organic compounds has been wrongly dismissed, as different organic compounds25, including evidence of peptide chains26, have been found in the 2-Ma upper terrace of Rio Tinto. In this paper, long-term preservation of complex molecular biosignatures in ancient Rio Tinto subsurface ferruginous deposits (~ 25Ma) is targeted. This task was accomplished by analyzing samples through state-ofthe-art time of flight secondary ion spectrometry (ToF–SIMS) with Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) and gas chromatograph mass spectrometry (G-MS). The identification of large molecular species suggesting that the acidic ferruginous materials have a high preservation potential agrees with the finding of peptide chains26 in the oldest Rio Tinto terrace, which also requires exceptional preservation conditions for its persistence. In this regard, preservation is favored in the Rio Tinto acidic deposits despite high water availability for mineral transformation. As this was not the case for sediments formed under low-pH conditions on Mars with limited water23, high preservation should be expected in those Hesperian materials. In this context, the extensive ferruginous deposits formed in acidic and oxidizing Hesperian environments should be considered first-class targets to look for ancient traces of life if any biosphere emerged at some time on the red planet. Geobiological settings. Gossan deposits occur in the Peña de Hierro area (Fig.1), which is considered one of the sources of the acidic waters of Rio Tinto12. Indirect rock dating using geomagnetic techniques27 suggests that the gossan materials in this location were formed through the acidic weathering of Carboniferous sulfide orebodies28 during the late Oligocene27. Such alteration coincides with a thermal climatic event29 that likely enhanced the acidic weathering of the orebodies. This was followed by a drop in the water table in response to subsequent climatic evolution during the Messinian saline crisis starting by 6 Ma27. The gossan originated from a complex direct and indirect chemolithotroph attack on pyrite (FeS2) orebodies30 through the following reactions: FeS2 (s) + H2O (l) + 3.5·O2 (aq) → Fe2+ (aq) + 2·SO42− (aq) + 2·H+ (aq) (direct biooxidation of sulfur). Fe2+ (aq) + 1/4·O2 (aq) + H+ (aq) → Fe3+ (aq) + 1/2·H2O (l) (biooxidation of released Fe2+). where O2 is used as an electron acceptor for the oxidation of FeS2 and free Fe2+. The ferric cation resulting from the microbial oxidation of Fe2+ becomes a secondary agent for pyrite oxidation through a subsequent abiotic process, where the sulfide anion is oxidized to sulfate30 to form Fe2+ again: FeS2 (s) + 14·Fe3+ (aq) + 8·H2O (l) → 15·Fe2+ (aq) + 2·SO42− (aq) + 16·H+ (aq). The solutions resulting from sulfide oxidation are eventually exposed to oversaturation, producing mainly acidic sulfates such as jarosite (e.g., H3O+Fe3(OH)6(SO4)2) and schwertmannite (Fe8O8(OH)6(SO4))23: 3·Fe3+ (aq) + M+ (aq) + 2·SO42+ (aq) + 6·H2O (l) ↔ MFe3(OH)6(SO4)2 (s) + 6·H+ (aq). where M = Na, K or H3O, and. 8·Fe3+ (aq) + SO42− (aq) + 14·H2O (l) ↔ Fe8O8(OH)6(SO4) (s) + 22·H+ (aq). The iron oxides that are the main mineral components of gossan and other underground ferric deposits in Peña del Hierro13 come from the diagenesis of primary acidic sulfates31 (Fig.2) when they are exposed to meteoric solutions with a circumneutral pH (~ 6–7). The process involves a combination of hydrolysis and dehydration that results in removing the sulfate groups to produce oxyhydroxides such as goethite (FeOOH) after dissolution-reprecipitation32. In addition, hematite (Fe2O3) can be directly formed from acidic sulfates by dehydration, which is accelerated by temperature31. Furthermore, hematite can also be formed by dehydration of nanocrystalline goethite32 or nanophase iron oxides23 whose composition is unclear. The diagenesis of materials formed under acidic conditions (Fig.2) is much more complex than a simple transformation of ferric sulfate to 3 Vol.:(0123456789) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ oxides. This diagenesis can involve different oxidation steps affecting primary ferrous sulfates (e.g., melanterite, FeSO4·7H2O) that precipitate in deeper areas under lower oxidizing conditions or through the recrystallization and hydrolysis of different sulfate phases by slight changes in the pH of the acidic solutions10,32. In the Peña del Hierro gossan, the abundance of hematite over other iron oxide species such as goethite13,23 suggests that the protolith was mainly composed of ferric sulfates that evolved to Fe2O3 during early diagenesis. Such a complex diagenesis that involves an interplay among dramatic pH and redox changes, hydrolysis and dehydration should have prevented the preservation of molecular compounds with a biological origin. However, the detection of abundant and diverse biomolecules in the Rio Tinto gossan shows that different preservation Figure1. Occurrence of the Upper Gossan deposits in the Peña de Hierro area. (a) Google Earth satellite image of the Upper Gossan unit topping the Peña de Hierro where the sample 010,109–1 was collected. (b) Outcrop photograph of the lower part of the Upper Gossan unit consisting of iron oxide deposits (FeOx) infilling cracks and voids left by the alteration of sulfide materials, which are limited by quartz-rich materials (Qz). (c) Sample 010,109was collected in the gossan material cropping out at a higher location in the Upper Gossan; it occurs as a lenticular body with internal massive fabric and topped with cm-sized popcorn and millimeter-thick filamentous structures. (d) Thin section of sample 010,109–1 around TA3 (see Fig.7) showing a laminated fabric not recognized through observation by the naked eye; this sample contains different 30-micron ovoid elements of likely fungal origin; because they correspond with the preservation of filamentous structures, the ovoid sections suggest that they have the same orientation. 4 Vol:.(1234567890) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ mechanisms have played essential roles in the fossilization of organic compounds that were produced by microbial activity in the subsurface of Peña de Hierro. Results To determine the preservation potential of ferruginous deposits that formed under acidic and oxidizing conditions, sample 010,109–1 (Fig.1b–d) collected from the top of the Upper Gossan unit in Peña de Hierro (~ 25Ma)13,22,33 was analyzed using different spectral techniques. Such analysis consisted of molecular mapping and characterization by means of three spectral techniques: FTICR-MS, GC–MS and ToF–SIMS. FTICR MS molecular fingerprinting. Ultrahigh-resolution Fourier transform ion cyclotron mass spectrometry was used to examine the extracted Upper Gossan fraction characteristic of the original preserved organic biogeochemical signature. Enabling highly accurate mass measurement and mass resolution, FTICR MS of terrestrial dissolved organic matter (TDOM) provides several thousand mass peaks for individual samples, of which several were assigned to C, H, N, O, and S compositional space34–36. The FTICR mass spectrum of the extracted Upper Gossan fraction is shown in Fig.3a. All mass signals detected correspond to singly charged ions and were assigned and classified as O, S, N-functionalized CHO, CHNO and CHOS molecular compositions (Fig.3b–f). In fact, these molecular series are visualized in a van Krevelen plot (Fig.3c), which is common to illustrate the molecular space of most complex natural samples34,35,37. The advantage of this plot is that it displays the assigned molecular species present in a sample rather than the computed average values. As shown in Fig.3c, the van Krevelen plot of the Upper Gossan extract exhibits regular smooth and extensive distributions over large ranges of unsaturation (high H/C ratio), down to high degrees of oxygenation (high O/C ratio). Condensed compounds were relevant (Fig.3c), had double bond equivalents (DBEs) higher than 3 with at least one or more phenyl rings, and were more abundant than saturated and unsaturated compounds with no phenyl ring (DBE ≤ 3) (Fig.S1). Moreover, the extracted gossan fraction revealed significant chemodiversity mainly attributed to CHO (39.0%), CHNO (51.0%), and CHOS (10.0%) molecular compositions (Fig.3d–f). While CHO-bearing molecules have shown the presence of saturated/unsaturated compounds (e.g., fatty acid-like, glycerins and oxygenated hydrocarbons CHOn>2), carbohydrate-like and carboxylic-rich alicyclic molecules (CRAM), including phytochemicals such as terpenes, flavonoids, lignins, polyphenol compounds and derivatives, as significant tracers of terrestrial organic matter (Fig.3d), CHNO-bearing molecules have shown enrichment in oxygenated compounds, peptide-like compounds and aliphatics (Fig.3e). Fewer CHNO-bearing molecules are localized Figure2. pH-Eh mineral stability diagram of different iron-bearing mineral phases and various clays. The blue arrows trace several diagenetic pathways starting in different Fe-bearing protoliths (e.g., sulfates and nanophase oxides), passing by an intermediate phase such as goethite (FeOOH), and ending in hematite (Fe2O3). The diagram shows that mineral maturation toward iron oxides requires hydrolysis and dehydration, which are controlled by changes in pH and Eh. The red arrows trace weathering pathways in the clays forming sedimentary deposits at the surface and subsurface of Rio Tinto. 5 Vol.:(0123456789) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ in the CRAM region (Fig.3e). Such enrichment in CHNO-bearing molecules suggests more preserved bioavailable organic matter and its derivatives, which are known as drivers of microbial community structure and fundamental metabolic strategies38–40. Conversely, CHOS-bearing molecules display a narrow range of hydrogen saturation H/C (0.55–2.2) evocative of reduced CHOS organic compounds and more oxygenated molecular species bounded by approximate O/C (0.25–0.75) and H/C (0.6–1.4) values characteristic of CRAM (Fig.3f). The appearance of CRAM and derivative compounds in the lignin-like and polyphenol-like regions on the van Krevelen plot could possibly arise from vegetation (i.e., vascular/higher plants and/or algal polyphenols). Indeed, the formation of organics by transitional organic matter involving monosaccharides, amino acids or fatty acids produces long-term humic organics via abiotic and/or humification reactions41. In parallel, the generation of insoluble, nonhydrolyzable and nonbiodegradable organics by various organisms may also be involved in the production and preservation processes of organic matter41,42. Overall, this large molecular diversity found in the Upper Gossan extract within and across the assigned molecular compositions is an example of the diverse biogeochemical availability of organic compound signatures, which is an acknowledged aspect of progressive organic matter preservation in the environment. GC–MS molecular screening. The analysis of the total lipid extract (i.e., including nonpolar, polar and acidic fractions) revealed the presence of diverse lipid families (i.e., normal and branched alkanes, including isoprenoids; alkenes; branched alkanols; steroids and hopanoids; and normal, branched, and monounsaturated fatty acids) (TableS1). The nonpolar fraction was mostly composed of straight chain or normal alkanes (n-alkanes) with 11 to 31 carbons that showed a bimodal distribution with maxima at C18-C21, C27 and C29. In conjunction with the n-alkanes, the isoprenoids pristane and phytane were also found (TableS1). The acidic fraction was composed of n-fatty acids with 11 to 30 carbons, monounsaturated fatty acids with 16, 18 and 20 carbons, branched fatty acids of iso/anteiso configurations (i.e., methyl groups at ultimate or penultimate positions) with 13 to 17 carbons (from i-/a-C13:0 to i/a-C17:0), and other minor compounds (TableS1). The n-fatty acid series showed a dominance of even-over-odd carbons, with maximum peaks at C24, C26 and C28. Meanwhile, the polar fraction was primarily composed of n-alkanols with 11 to 24 carbons and even-over-odd preference, glycerol compounds suggesting the occurrence of glycolipids, a lupane derivative (lupan-3-ol, acetate), and a series of sterols, including cholesterol β-sitosterol, stigmasterol, and an ergostanol derivate (TableS1). The observed lipid distribution patterns provide information about possible sources of organic matter in this sample, which are mainly related to bacteria, algae/phytoplankton and plants. The latter source is inferred from the detection in the sample of odd and high molecular weight (HMW) n-alkanes such as C25 for macrophytes Figure3. Main results obtained from the ESI(-)-FTICR-MS analysis. (a) Mass spectrum (150–950Da) of the Upper Gossan sample 010,109–1. Inset number in (a) represents the total number of assigned molecular members (i.e., CHO, CHNO and CHOS). (b) and (c) H/C ratio versus mass-to-charge (m/z) plot and van Krevelen diagram of all assigned molecular members, respectively (color-coded dots represent molecular members CHO, blue; CHNO, orange; and CHOS, green, while circular areas indicate relative mass peak intensity). Inset in (b) represents the frequency histogram of the exact mass differences corresponding to CH2, H2, On(1≤n≤5), NHn(1≤n≤3), and SOn(1≤n≤4) (error range 0.1ppm). (d), (e) and (f) represent compound classes included in CHO, CHNO and CHOS-molecular members, respectively. Highlighted areas in (d), (e) and (f) visualize the following compound classes: (1) Saturated/unsaturated compounds, (2) Amino acids-/peptide-like compounds, (3) Carbohydrate-like compounds, and (4) highly unsaturated and aromatic compounds. 6 Vol:.(1234567890) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ and C27, C29, and C31 for terrestrial plants43,44. In contrast, the detection of low molecular weight (LMW) n-alkanes (< C22) is associated with either bacterial (mostly even LMW n-alkanes) or phytoplanktonic (odd LMW n-alkanes) sources45. Pristane and phytane have as a primary source the degradation of chlorophyll pigments. Both isoprenoids mainly originate from phytol, the esterifying alcohol of phototrophic chlorophylls46, whose degradation gives rise to pristane or phytane in the presence or absence of oxygen, respectively47. Vegetal and microbial biomarkers are also found in the acidic and polar fractions, where LMW n-fatty acids and n-alkanols (i.e., < C20), as well as the hopanoid lupan-3-ol and acetate, are most likely associated with bacterial sources48, whereas HMW n-fatty acids and n-alkanols, as well as phytosterols (i.e., β-sitosterol or stigmasterol) can be related instead to macrophytes or land plants43,49. Sterols are ubiquitous constituents in all eukaryotic organisms, where stigmastanol or β-sitosterol is produced almost exclusively by higher plants50, whereas cholesterol is primarily produced by animals51. The presence of ergosterol can be related to fungal species52. Certain monounsaturated fatty acids detected in the sample have been associated with biomass from cyanobacteria or microalgae (i.e., 16:1 ω7, 18:1 ω9)53. Furthermore, iso-/anteisopairs of fatty acids (mostly 15:0 and 17:0) are typically associated with bacterial sources54 and are particularly abundant in sulfate-reducing bacteria55. ToF–SIMS molecular mapping. ToF–SIMS molecular imaging has allowed the characterization of different groups of molecules that are preferentially associated with distinctive morphologies, the sizes of which are generally smaller than 500μm (see Table1 where groups and morphologies are associated with compounds). Such morphologic groups show different fabrics and textural features when the target areas (TAs) 1 to 3 (TA1, TA2, and TA3) are analyzed (Figs.3, 4, and 5; Figs.S2–S5). The distribution of ions with biomolecular origin follows the fabric and textural pattern outlined by the main inorganic ions of acidic oxysalts and oxyhydroxides (Figs.S2 to S3), including Fe+ (m/z 55.93), Al+ (m/z 26.98), K+ (m/z 38.96), Na+ (m/z 22.99), SO2-/SO3- (m/z 63.96/79.96), and SiO2- (m/z 63.96), while PO2-/PO3- (m/z 62.96/78.96), and NO2-/NO3- (m/z 45.99/61.99) anions occur as microstructures smaller than 50 microns (Fig.S4). In sample 010,109–1, TA fabric and texture fit the distribution of inorganic and organic ions well (Figs.4, 5, 6, and 7; Figs.S2–S4), which mainly correspond to filamentous networks and isolated elements associated with organic compounds of biological origin (e.g., see Figs.4, 5, 6, and 7). In some cases, ToF–SIMS ion mapping is the only technique for recognizing the microstructures of samples when this is not possible using different microscopic techniques (e.g., light or scanning electron microscopes). Furthermore, the TOF–SIMS ioninduced secondary electron images also provided additional information on the microstructure of the sample Table 1. Identification of the morphological groups through the ToF–SIMS molecular distribution found in the Upper Gossan sample of Peña de Hierro. Group Microstructure Potential compounds Target area (TA)m/z positive ions m/z negative ions Group 1 Loosely to densely packed networks of < 50 micron thick filaments with strong mineralization shown by positive and negative sulfurbearing ions Lineal hydrocarbons fragmented and ionized from long lipids which form adducts with NH4+ likely sphingolipids Fatty acids coming from lipidic structures like ceramides TA1, TA2, TA3 Salt composition of mineralized structures SOx-bearing adducts of cholesterol-like sterols and lipids coming from the mineralized surface of fungal structures Group 2 Linear elements with an average thickness of > 80 microns disassociated from the other microstructures Likely terpenoid fragments sourced in plants TA1 Group 3 Ovoidal microstructures averaging 60 microns (TA1 and TA2); following larger microstructures and the mineral matrix in TA3 Fragments of amino acids and peptides. The occurrence of some ions in the TA3 matrix while others follow the biological microstructures suggest different a origin for them TA1, TA2, TA3 Group 4 Two coupled elongated nodular microstructures in which the elongation axes connect both nodule centers (TA1); occurring as an empty cast in TA3 Unknown origin n/a TA1, TA3 Group 5 600 micron-long elongated triangle-like morphology crosscutting the main structure fabric of Groups 1 and 2; occurring as an empty cast in TA3 Plant terpenoids like piperonal derivatives TA1, TA3 Group 6 Three 120-micron subcircular morphologies show an empty central void; occurring as an empty cast in TA3 Vanillic acid derivatives n/a TA1, TA3 Group 7 10to 50-micron micronodules following the main pattern appearing in Groups 3 and 4 Plant terpenoids TA1 Group 8 Structures with short branches with an average length of 200 microns that are transverse to the fabric of Group 1 n/a Tannins, or some polysaccharides TA1 Group 9 Densely packed network of 40-micron-long filament-like microstructures, which appear in TA2 along the sharp boundary that contacts Group 1; fits the primary microstructure found in TA3, while in the TA1 show a more homogeneous occurrence Glycerolipids and glycerophospholipids with possible occurrence of glycerophosphocholines. Fatty acids fragmented from the lipidic structure. Abundant in TA2 and TA3, and secondary in TA1 TA2, TA3 7 Vol.:(0123456789) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ Figure4. Occurrence of different biological microstructures in Target Area 1 (TA1) characterized by the imaging capabilities of ToF–SIMS as morphological Groups 1 to 8 through their molecular compositions. (a) Groups 1 to 8 outlined by the distribution of (a) positive and (b) negative fragments. The red rectangle traces a nodular microstructure in G3 occurring as an empty and complementary cast in G2, G5 and G6. The white rectangle encloses branched filamentous structures in G8 that complement part of the filamentous structures in G1 (see Fig.5) suggesting that the fungal filamentous structure show compositional variations. The white arrow points at a void (see Supplementary Fig.2) that is affecting some microstructures in G2 and G5. As discussed in the main text, the composition and structure of groups G1 to G8 suggest that the groups could come from the preservation of different biological structures such as filamentous fungi (G1 and G3), or plant fragments (G5 and G6). 8 Vol:.(1234567890) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ (Figs.S2–S4). By combining all the data sources retrieved via ToF–SIMS, the following morphological groups could be identified (Figs.4, 5, 6, and 7). Group 1 This group is characterized by networks of < 50-micron-thick filaments (Figs.4, 5, 6, and 7) showing strong mineralization by the enrichment in positive and negative sulfur-bearing ions such as Na3SO4+, K2NaSO3+, HS-, SO2-, and SO3at 164.92, 196.86, 32.98, 63.96, and 79.96 (Figs.S2–S5), respectively, as well as different cations such as Fe+ (55.94), K+ (38.97), Na+ (22.99) and Al+ (26.98) (Figs.S2 to S5). TA1 filamentous structures are thicker and looser (Fig.4) than filament networks found in TA2, which are thinner, overlap more tightly, and show 20-micron nodular ends (Fig.6). In TA3, the distribution of positive and negative ions outlines an elliptic microstructure that is perpendicular to straight laminar projections (Fig.7). Group 1 shows maximum intensities in different molecular groups, including a series of N-bearing compounds (e.g., m/z peaks at 268.30 and 282.31 corresponding to C18H38N+ and C19H40N+) and NH4+-bearing positive adducts found at 296.33, 494.56, 522.59 and 550.62 fitting C20H42N+, C34H72N+, C36H76N+, and C38H80N+, respectively. This group is also followed by a series of negative ions identified as fatty acids and/or ester carboxylates (e.g., m/z 197.15, 199.17, 227.20, or 255.23, corresponding to C12H21O2-, C12H23O2-, C14H27O2-, or C16H31O2-, respectively) and Figure5. Couple of ToF–SIMS red/green/blue (RGB) overlays for (a) G3, G6 and G5 and (b) G1, G3, and G8 resulting in two different microstructure associations. The first RGB association (a) could correspond to compositional variations in plants (G5, G6, and G3) or to plant fragments (G5, and G6) affected by microbial activity (G3). Interestingly, G4 appears as a compositional variation in the G3 structure. The second RGB association (b) could include three different groups, namely, G1, G3 and G8 that could correspond to a microbial structure such as a biofilm in which the biomass is dominated by microbial forms biomineralizing sulfate. G3 structures interlock empty complementary areas in G1, suggesting that these areas could be the remains of different biological entities, while G8 branched microstructures merge with some filamentous units in G1, suggesting compositional variations in some specific zones of the filamentous mass of G1. 9 Vol.:(0123456789) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ Figure6. Characterization of possible microbial structures preserved in Target Area 2 (TA2) by using the ToF–SIMS molecular mapping capabilities. The image analysis follows the same morphological groups G1 to G8 adding one new group G9. Two different main structures defined by G1 and G9 are limited by a sharp boundary between them. The G1 microstructure is characterized by a strong intensity of different Nand S-bearing compounds that outline a diffuse filament network, whose primary composition would be large lipids such as ceramides. However, G9 shows a more massive filamentous mass that is traced by a high intensity in acylglycerides and FA fragments. G3 occurs as an independent group of microstructures that interlock both G1 and G9 microstructures. The RGB overlay shows two main microbial structures (G1 and G9) in which empty cavities would be occupied by a third microbial unit defined by G3. Interestingly, the G1 microbial structure is strongly mineralized by sulfate while G3 and G9 microstructures are mineralized by ferruginous materials, suggesting that they were formed in two separate mineralization episodes. The microbial structure and composition of G1 and G9 suggest that these groups could correspond to mineralized biofilms. 16 Vol:.(1234567890) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ analysis of the gossan sample strongly supports organic production from multiple sources that include aquatic (i.e., macrophytes) and land plants, bacteria, algae/phytoplankton and fungi. Both FTICR-MS and GC–MS results are highlighted by the ToF–SIMS imaging capabilities, which have allowed the characterization of microbial communities through the distribution of distinctive biomolecular fragments. Such distributions distinguish different imaging groups (from 1 to 9), whose morphologies suggest that they correspond to parts of living organisms, microbial communities and/or biofilms92. These morphologies are observed in the molecular imaging of the different target areas that captures the occurrence of several biological structures. In fact, the G1 group in TA1 (Fig.4) shows thick filamentous structures with a high content of N-bearing positive fragments that can be attributed to sphingolipids. The ToF–SIMS spectral analysis of TA1 also shows high enrichments in different sets of positive ions that have been assigned to plain ceramides and glucosylceramides (e.g., cerebrosides and sphingomyelins). As such, an association of biomolecules is currently found in microbial fungi82,93 that occur as filament clusters (hyphae) displaying a magnitude comparable to that of the filamentous structures defined by N-bearing fragments as G1 (Fig.4). The filamentous structures are also well defined by high intensities of different inorganic SOn (2 ≤n ≤4) anions and organic sulfate adducts, which strongly suggests hyphae biomineralization, as commonly occurs in some modern Rio Tinto fungi14. The detection of high contents of sulfate adducts likely results from the binding between SO42 and HSO4anions with a positively charged fungal hyphae surface, which is partly formed by lipids. In fact, the occurrence of cholesteryl sulfate defined by a large peak (I > 400 to 1050 cps) at m/z 465.30 in the different target areas (Figs.4, 5, and 6) strongly suggests that the sulfate organic adducts result from the ionization of the biomineralized fungal surface by sulfate anions. Organic sulfate-bearing adducts are not found in biological structures other than the filamentous fungal structures, suggesting that the sulfate biomineralization is highly specific to some fungal microorganisms14. As discussed above, the different biological structures defined by morphological groups 1 to 9 are emphasized by the highly distinctive molecular composition distribution. The red/green/blue (RGB) overlay of the different morphological groups (Figs.5, 6, and 7) shows that most of the biological structures interlock with each other to form larger biological assemblages. Although the origin of some of the preserved structures is unknown, the morphology and composition of most of them strongly suggest that they correspond to an association of different fungal structures (e.g., G1, G3 and G7) that are associated with parts of other living forms, such as plants represented in both Groups 5 and 6 (Figs.4 and 5). As shown in the results, G1, G3 and G7 are formed by N-bearing fragments and sulfur-bearing adducts sourced in sphingolipids and biomineralized bodies by sulfate (Fig.2; TableS2), which corresponds to the preservation of fungal communities14,80,93, while G5 and G6 are defined by different ions fitting the fragmentation pattern of some phytochemicals, such as vanillic acid and piperonal56,57. In terms of biomass, the main fungal structure corresponds to G1, which is characterized by dozens of fragments (TableS2) that usually have a high intensity as well. This feature is very common in the different N-bearing m/z cations that are characterized as sphingolipid fragments and secondary anions (e.g., G8) that could have been produced by the fragmentation of large polysaccharides forming fungal biofilms94. The occurrence of exceptionally preserved biofilms with two different microbial forms is observed in TA2 (Figs.6 and 8; Figs.S10 and S11), where two different biological structures are found. They are distinguished by the molecular composition that defines G1, which is enriched in ceramides and sulfur adducts, versus G9, which is characterized by a high abundance of diacylglycerides, through a sharp boundary between the two groups (Fig.6). G3 interlocks G1 (Fig.6) in the same way as in TA1, which suggests that there are compositional changes in the fungal biofilm. Furthermore, in TA3, ToF–SIMS image analysis shows additional biological structures whose morphology and molecular composition suggest a fungal origin but also contain glycerolipids in addition to sphingolipids and SO4-mineralized compounds (Fig.7). As discussed above, the fungal structures defined in TA1 and TA2 show high concentrations of sphingolipids and SO4-mineralized compounds but depletion in Figure10. Fragmentation pattern of unknown compounds defining G8. The distribution of fragment masses suggests the presence of biological compounds with high contents of oxygen and/or nitrogen, which correspond with polysaccharides and polyketides. 17 Vol.:(0123456789) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ glycerolipids (Figs.4 and 6), which suggests that the biological structure found in TA3 corresponds to a different biological origin. The phytochemicals (e.g., vanillic acid and piperonal derivatives) that characterize some biological structures in TA1, such as G5 and G6 (Fig.4), occur in TA3 within the mineral matrix in which the biological structures are embedded (Figs.1 and 7). Such a distribution suggests that they were transported in the acidic solution after being released from plant fragments in more proximal and/or surficial areas where the plants were growing. They come together with a set of different PAHs (Fig.7; TableS3), which are recalcitrant compounds resulting from the photochemical process and are common in continental waters. This agrees with the parallel arrangement of the fungal structures (Fig.1), which suggests that they were oriented along an acidic flow with a high concentration in ferric oxysulfate polymers95. Indeed, the ToF–SIMS image analysis provides enough evidence in the form of chemical compositions to gain insight into the mineralization processes that led to the preservation of biological structures and species. This preservation was accomplished by means of different mineralizing episodes that started with sulfate biomineralization by the activity of some fungal and bacterial communities14. Sulfate biomineralization was followed by ferric-bearing mineralization, whose composition could have been mediated through different phases, including schwertmannite, hydronium jarosite, and nanophase goethite23. The mineralization in this case was a passive process in which biological structures were coated and infilled with ferric-bearing compounds whose composition was controlled by the solution hydrochemistry. The rock maturation ended in partial replacement of the primary phases by goethite and hematite (Fig.2). Such limited maturation is evidenced by high concentrations of sulfate and organic carbon, which were determined by ToF–SIMS analysis (Figs.4, 5, 6, and 7). In this context, the precipitation of the ferric and sulfate phases played an essential role in the preservation of the structure and composition of the biological elements in three ways: (1) building rigid mineralized structures resistant to sediment load or tectonic stress and chemical attack96,97, (2) limiting oxygen availability and thus inhibiting the degradation of biological compounds by a high concentration of cations98 that are very harmful for microbial life99, and (3) a fast decrease in water activity during mineral precipitation100,101. Conclusions The spectral and ToF–SIMS image analyses of the Upper Gossan sample show that hyperacidic and oxic conditions leading to ferric deposits have remarkable preservation potential in terms of structure and composition. Such a preservation potential is also supported by the occurrence of peptidic chains in 6to 2.2-Ma Rio Tinto terrace deposits26. The same low-pH and oxidizing conditions have been hypothesized for some regions on Mars during the Noachian-Hesperian transition6,7. In this period, different mineral phases, such as jarosite, schwertmannite, nanophase goethite, and other unknown Fe2+/Fe3+-bearing sulfate or oxide phases9,102,103, could have played essential roles in the mineralization of microbial remains if they were ever produced on Mars. The precipitation of ferric phases has been reported to follow fast saturation7, which is observed in the Rio Tinto Mars analog. The saturation of acidic solutions on Mars would have increased the concentrations of different cations, including toxic metals98,99, and decreased the water activity100, both of which greatly inhibit microbial activity. Indeed, the water activity decreases with temperature104, so the net preservation of biomolecular compounds would be greatly favored at the low temperatures on Mars. The deposit maturation to goethite and hematite would confer on them a higher resistance to deformation and weathering. High concentrations of ferric ions in rock and solutions105 would also have provided protection against UV radiation once the microbial remains were eventually produced on Mars. Although preservation in clay-rich deposits is currently associated with organic preservation21, the preservation potential can be greatly reduced when clay deposits are exposed to different maturation and alteration processes that can degrade organic compounds in later stages of rock diagenesis on Mars, such as those associated with enhanced weathering (Fig.S14), which likely occurred during the Noachian to Hesperian transition. The solutions involved in clay formation would have not had high concentrations of cations that could eventually mineralize the biological structures and protect them against UV radiation105 and lithostatic compression97. When mineralization does not occur, biomolecular compounds can be broken into smaller compounds by lithostatic compression. Indeed, such mild environments are associated with unfavorable mineralizing conditions, which implies low potential for preserving structures. Assuming that Mars was exposed to acidic episodes, there is a high possibility that if life emerged on the planet, its molecular preservation might have followed the same mechanisms to record traces of ancient life on the surface and in the interior of the planet. In fact, the search for life on Mars could be expanded to large and extensive acidic deposits formed billions of years ago. Material and methods Sample collection. A field campaign was conducted for surface sampling of the upper part of the Cenozoic gossan13 topping the Peña de Hierro Carboniferous basement (Fig.1a). This unit corresponds with the oldest alteration episode that formed the Upper Gossan materials through sulfide weathering by a water table that remained relatively stable during the late Oligocene-early Miocene27. Sample 010109–1. was collected in Peña de Hierro (GPS coordinates 37.724918°, 6.553717°) from the uppermost part of the iron oxide deposits corresponding to the Upper Gossan unit aged late Oligocene27. The collected materials consist of massive iron oxide strata-like deposits forming ~ 1-m-thick lenses that are associated with the central and highest part of the gossan (Fig.1a). Such materials currently vary downward or laterally to brecciated or tuffaceous oxide deposits (Fig.1b–c), showing that the gossan formation involves a complex interplay between surface and subsurface processes. In this regard, such massive deposits are currently delimited at the top by efflorescent pseudomorphs and coated filamentous structures (Fig.1c) that could result from a shallowing evolution from subsurface to surface conditions. The sample was collected directly from 18 Vol:.(1234567890) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ an outcrop showing few signs of alteration, and a geological hammer was used to prevent contamination that can be produced through drilling, which employs water as a coolant. The sample was carefully collected using nitrile gloves, covered by aluminum foil and introduced into a sterilized sampling bag to prevent contamination. Then, the sample was split into three parts; subsequently, each subsample surface was sterilized by igniting mass spectrometry-grade ethanol (sigma). Then, the three subsamples were delivered for spectrometry analysis. Observations under an optical microscope and a scanning electron microscope (SEM) of thin sections of the massive unit provided evidence for the occurrence of abundant biological fabrics that were templated by coating and mineralization (Fig.1c). Furthermore, the sample observations under the different microscope techniques also evaluated the origin of the biological traces that were found in the rock matrix (Figs.S14 and S15). Interestingly, scanning electron and energy dispersive spectroscopy (SEM–EDS) analysis suggested very slight, if any, alteration of the mineral matrix of the sample. Indeed, the lack of different cations, such as Ca, Mg and K, and the low content of Si (Figs.S15 and S16) indicate that no alteration occurred12. In this regard, Fernandez-Remolar etal.106 reported the formation of silica laminations. Such laminations might have formed if the gossan material had been exposed to subsequent alteration associated with the occurrence of secondary biomolecules that were younger than the rock. Indeed, the SEM–EDS images show the presence of colomorphic goethite (Figs. S15 and S16) covered by laminated hematite, which is the result of long sedimentary maturation of the primary iron oxysulfates and/or oxyhydroxides13,23. Furthermore, the EDS analysis implies that goethite is associated with traces of S (see spectrum 3 in both Figs.S15 and S16) that likely came from a primary mineral composition of sulfates. Interestingly, SEM–EDS also targeted goethite with microstructures (Fig.S16) that were very similar to the filamentous microstructures shown in Figs.4 and 5. The SEM–EDS image shows that the filamentous goethite is enclosed inside a hematite-rich lamina (Fig.S16), which indicates an old age and supports a primary origin for the preserved biological structures. Sample preparation. Sample analysis was performed using different state-ofthe-art analytical techniques, including time of flight secondary ion mass spectrometry (ToF–SIMS) and Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS)26,36. These methods were supported by the more conventional technique of gas chromatography coupled with mass spectrometry (GC–MS), all of which have access to molecular databases with many compounds. In this regard, GC–MS is a technique that provides essential information to confirm compound identification. The preparation for the sample analyses using these three distinctive spectrometry techniques followed two different processes: simple thin sectioning and organic extraction from powdered samples. Thin sections of sample 010109–1 were prepared for mineral and textural analysis under an optical microscope and a Scanning Electron Microscope (SEM) coupled with an Electron Dispersive Spectroscopy (EDS) probe, and molecular surface analysis by ToF-SIMS107. They were polished an extra time by using a 0.3μm alumina paste to reduce surface imperfections that produce frequent interference during surface analysis. On the other hand, as FTICR-MS and GC–MS analyses of compounds require organic extraction, sample 010,109–1 was initially crushed under sterile conditions to prevent contamination. Microscope analysis. The inner structure of sample thin sections was observed under an upright trinocular microscope Motic BA400 coupled to a MoticamPro 282a camera. Optical observations were complemented with backscattered electron imagery using a Jeol JSM-5600 LV Scanning Electron Microscope (SEM). In addition, semiquantitative chemical microanalysis was performed with a Cambridge INCAx SIGHT EDS (Electron Dispersive Spectroscopy) to associate microstructure and mineral substrate with elemental composition. GC–MS and FTICR‑MS analyses of organic extracts. Before lipid extraction, the sample was externally cleaned, and the outer layer was removed to avoid cross contamination with allochthonous material from the surface. In the first step, an external layer approximately 0.5cm thick was removed with a solvent-cleaned (methanol and dichloromethane) saw and discarded. Subsequently, the remaining inner sample was externally cleaned with a solution of approximately 100ml of hydroxide potassium (6% in Milli-Q water) while being gently shaken for 15min in a glass beaker. After rinsing with Milli-Q water, the sample was then cleaned with nitric acid (6%) under the same conditions (shaking for 15min). Finally, the sample was rinsed with Milli-Q water and ethanol (> 99.5% purity) and then dried in an oven (45°C). Once dried, the sample was ground with a solvent-cleaned mortar and pestle for lipid analysis. Approximately 100g of sample 010,109–1 was prepared for organic-solvent extraction for GC–MS analysis (Fig.S13). The extraction was performed with a mixture of dichloromethane/methanol (DCM/MeOH, 3:1, v/v) for 24h with a Soxhlet apparatus. Internal standards (tetracosane-D50, myristic acid-D27, and 2-hexadecanol) were added prior to the extraction. The total lipid extract was concentrated by rotary evaporation to approximately 2ml. After this process, activated copper was added, and the sample was allowed to stand overnight for elemental sulfur removal. The extract was separated into three fractions of different polarities (nonpolar, polar, and acidic). In the first step, a Bond-elute column (bond phase NH2, 500mg, 40µm particle size) was used to separate the neutral and acidic fractions by eluting first with 15ml DCM/2-propanol (2:1, v/v) and subsequently with 15ml diethyl ether with 2% acetic acid. In the second step, the neutral fraction was further separated into nonpolar and polar fractions by using an Al2O3 column (activated, neutral, 0.05–0.15mm particle size) and eluting with 5ml of hexane/DCM (9:1, v/v) and with 4ml of DCM/methanol (1:1, v/v), respectively. The acidic fraction was derivatized with BF3 in methanol, and the polar fraction was derivatized with N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA). Together with the sample, a blank was extracted (i.e., an empty cellulose thimble), concentrated, fractionated and analyzed in parallel as a negative control. All glassware and sorbents were heated in a furnace at 550°C for 24h. All solvents were provided by Sigma Aldrich (Spain) and were of analytical grade (> 99.5% purity). GC–MS analysis was performed to characterize the different nonpolar, acidic, and polar fractions. This analysis used a 6850 GC system and a 5975 VL MSD with a triple axis detector (Agilent Technologies). The recovery 19 Vol.:(0123456789) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ of spectral data was obtained using full scan mode (m/z 50 to 650) that operated with electron ionization at 70eV. The analytes were injected (1µl) into an HP-5MS column (30m × 0.25mm i.d. × 0.25µm film thickness) with He as a carrier gas (1.1ml min-1). For the nonpolar fraction, the oven temperature was programmed from 50 to 130°C at 20°C min−1 and then to 300°C at 6°C min-1 (held for 20min). For the acidic fraction, the oven temperature was programmed from 70 to 130°C at 20°C min−1 and then to 300°C at 10°C min−1 (held for 10min). For the polar fraction, the oven temperature program was the same as that for the acidic fraction, except that the oven was held for 15min at 300°C. The injector temperature was 290°C, the transfer line was 300°C, and the MS source was 240°C. The identification of compounds was based on the comparison of retention time and mass spectra with the National Institute of Standards and Technology (NIST) mass spectral database and reference materials: n-alkanes (C10 to C40), fatty acid methyl esters (FAMEs; C8 to C24), n-alkanols (C10, C14, C18, and C20), and branched isoprenoids (2,6,10-trimethyl-docosane, crocetane, pristane, phytane, squalane, and squalene). All chemicals and standards were supplied by Sigma Aldrich and Chiron. The GC–MS data were acquired and processed by Agilent MSD ChemStation software, which performed molecular identification though NIST library matches of peak spectra. To obtain the molecular information for sample 010,109–1, a solid-phase extraction was performed on the supernatant as described elsewhere36,37, and the extract was used for FTICR-MS analysis. This analysis was conducted with a SolariX FTICR mass spectrometer (Bruker Daltonik GmbH, Germany) equipped with a 12T superconducting magnet (Magnex Scientific Inc., GB) and an APOLLO II electrospray ionization source (Bruker Daltonik GmbH) in negative ionization mode. The injection was performed using a microliter pump at a liquid flow rate of 120ml h-1. Both sheath gas and curtain gas consisted of nitrogen. A source heater temperature of 200°C was maintained to ensure rapid solvent evaporation in the ionized droplets. The FTICR mass spectrometer was precalibrated using known arginine clusters. The mass spectrum was zero filled to a processing size of four megawords using 750 scans within a mass range of m/z 147.4 to 1000 and internal calibration with an appropriate reference mass list, reaching accuracy values lower than 100ppb. A blank sample was run before sample injection and did not show any clusters of contamination. Data treatment was performed using Compass Data Analysis 4.0 (Bruker Daltonics, Bremen, Germany). The exact masses of the detected m/z ions were defined, and their molecular formulas were calculated using a software tool written in-house109. The final formulas were branched into groups containing CHO, CHNO and CHOS molecular compositions, which were used to reconstruct the group-selective mass spectrum. ToF–SIMS analysis. ToF–SIMS is an analytic surface technique that characterizes the distribution of positive and negative ions in a sample. For this reason, it can provide essential information about the preservation of organic compounds, establishing a direct association between the mineral matrix and the occurrence of molecular biosignatures in the sample108. ToF–SIMS has a high detection level for different compounds, including contaminants that can potentially be attached to the thin section surface. To prevent this problem, the thin sections of sample 010109–1 were cleaned by surface sputtering with a 100 nA 500 eV oxygen ion beam for 3 s. The molecular surface analyses were performed by using a ToF–SIMS IV (ION-TOF, Münster, Germany) system under a pressure of 5·10–9mbar. The thin section was sputtered with a pulsed bismuth liquid metal ion source (Bi32+) at 50keV. The ion beam was operated with a 20ns pulse width and 0.3 pA pulsed ion current for a dosage lower than 5·109 ions·cm-2, below the threshold level of the static SIMS conditions (1·1011 ions·cm−2). The released secondary ions were detected with a reflectron time-of-flight analyzer, a multichannel plate (MCP), and a time-to-digital converter (TDC). Charge neutralization was achieved with a low energy (20eV) electron flood gun. Secondary ion spectra were obtained from a randomly rasterized square surface area with a 500-μm side within the thin section surface. Secondary ions were extracted with a 2kV accelerating voltage and were postaccelerated to 10keV kinetic energy just before hitting the detector. Furthermore, the observations of microfabric and microtextural features of the inspected areas for each sample were accomplished by the ToF–SIMS ability to collect ion-induced secondary electron images. This was the first step that was taken during the ToF–SIMS analytic process that searched for differences in fabric and textural properties along the polished surface of the sample. As a result, different target areas showing characteristic fabrics and textures were selected in sample 010,109–1 (Figs.2, 3, 4, 5, 6, and 7; Figs.S2 to S3). The characterization of the fabric and texture through the imaging capabilities of ToF–SIMS was essential to recognize the association among fabric/texture, mineralized microstructures, and distribution of preserved biological compounds. Mass spectral acquisition and image analysis, including the overload of red, green and blue (RGB) images for molecular mapping, were performed with ION-TOF Surface Lab software (version 6.8). The analysis of the molecular fragmentation pattern was completed through the Chemical online tool109 and the open source mass spectrometry tool Mmass110. Chemspider, METLIN and LIPID MAPS Structure databases111–113 were used as sources of information to identify molecular fragments and compounds. Received: 14 April 2021; Accepted: 11 October 2021 References 1. Hassler, D. M. et al. Mars’ surface radiation environment measured with the mars science laboratory’s curiosity rover. Science 343, 1244797. https:// doi. org/ 10. 1126/ scien ce. 12447 97 (2014). 2. Leuko, S. et al. On the stability of deinoxanthin exposed to mars conditions during a long-term space mission and implications for biomarker detection on other planets. Front. Microbiol. 8, 1680–1680. https:// doi. org/ 10. 3389/ fmicb. 2017. 01680 (2017). 3. Michalski, J. R. et al. The Martian subsurface as a potential window into the origin of life. Nat. Geosci. 11, 21–26. https:// doi. org/ 10. 1038/ s415610170015-2 (2018). 20 Vol:.(1234567890) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ 4. Orosei, R. et al. Radar evidence of subglacial liquid water on Mars. Science (2018). 5. Des Marais, D. J. et al. The NASA astrobiology roadmap. Astrobiology 8, 715–730. https:// do i. org/ 10. 1089/ ast. 2008. 0819 (2008). 6. Ehlmann, B. L. & Edwards, C. S. Mineralogy of the martian surface. Annu. Rev. Earth Planet. Sci. 42, 291–315. https:// doi. org/ 10. 1146/ annur evearth060313055024 (2014). 7. Hurowitz, J. A., Fischer, W. W., Tosca, N. J. & Milliken, R. E. Origin of acidic surface waters and the evolutionof atmospheric chemistry on early Mars. Nat. Geosci. 3, 323–326. https:// doi. org/ 10. 1038/ ngeo8 31 (2010). 8. Michalski, J. R. & Niles, P. B. Deep crustal carbonate rocks exposed by meteor impact on Mars. Nat. Geosci. 3, 751. https:// doi. org/ 10. 1038/ ngeo9 71 (2010). 9. Lichtenberg, K. A. et al. Stratigraphy of hydrated sulfates in the sedimentary deposits of Aram Chaos. Mars. J. Geophys. Res. Planets https:// doi. org/ 10. 1029/ 2009J E0033 53 (2010). 10. Fernández-Remolar, D. C. et al. The environment of early Mars and the missing carbonates. Meteorit. Planet. Sci. 46, 1447–1469. https:// doi. org/ 10. 1111/j. 19455100. 2011. 01238.x (2011). 11. Wray, J. J. et al. Orbital evidence for more widespread carbonate-bearing rocks on Mars. J. Geophys. Res. Planets ht tps:// doi. org/ 10. 1002/ 2015J E0049 72 (2016). 12. Fernández-Remolar, D. C. et al. Río Tinto sedimentary mineral assemblages: A terrestrial perspective that suggests some formation pathways of phyllosilicates on Mars. Icarus 211, 114–138. https:// doi. org/ 10. 1016/j. icarus. 2010. 09. 008 (2011). 13. Fernández-Remolar, D. C. et al. Underground habitats in the Río Tinto basin: A model for subsurface life habitats on mars. Astrobiology 8, 1023–1047. https:// doi. org/ 10. 1089/ ast. 2006. 0104 (2008). 14. Oggerin, M., Tornos, F., Rodríguez, N., De Moral, C. & Sánchez-Román, M. Specific jarosite biomineralization by Purpureocillium lilacinum, an acidophilic fungi isolated from Río Tinto. Environ. Microbiol. 15, 2228–2237 (2013). 15. Knoll, A. H., Canfield, D. E. & Konhauser, K. O. Fundamentals of geobiology. 480 (Blackwell Publishing Ltd, 2012). 16. Floyd, M. A. M., Williams, A. J., Grubisic, A. & Emerson, D. Metabolic processes preserved as biosignatures in iron-oxidizing microorganisms: Implications for biosignature detection on mars. Astrobiology 19, 40–52. https:// doi. org/ 10. 1089/ as t. 2017. 1745 (2018). 17. Sumner, D. Y. Poor preservation potential of organics in Meridiani Planum hematite-bearing sedimentary rocks. J. Geophys. Res. Planets https:// doi. org/ 10. 1029/ 2004J E0023 21 (2004). 18. Summons, R. E. et al. Preservation of Martian organic and environmental records: Final report of the mars biosignature working group. Astrobiology 11, 157–181. https:// doi. org/ 10. 1089/ ast. 2010. 0506 (2011). 19. Amils, R. & Fernández-Remolar, D. C. Río Tinto: A geochemical and mineralogical terrestrial analogue of mars. Life 4, 511–534 (2014). 20. Lasne, J. et al. Oxidants at the surface of mars: A review in light of recent exploration results. Astrobiology 16, 977–996. https:// doi. org/ 10. 1089/ ast. 2016. 1502 (2016). 21. Ehlmann, B. L. et al. Clay minerals in delta deposits and organic preservation potential on Mars. Nat. Geosci. 1, 355. ht tps:// doi. org/ 10. 1038/ ngeo2 07 (2008). 22. Fernández-Remolar, D. C. & Knoll, A. H. Fossilization potential of iron-bearing minerals in acidic environments of Rio Tinto, Spain: Implications for Mars exploration. Icarus 194, 72–85 (2008). 23. Fernández-Remolar, D. C., Morris, R. V., Gruener, J. E., Amils, R. & Knoll, A. H. The Rio Tinto Basin, Spain: Mineralogy, sedimentary geobiology, and implications for interpretation of outcrop rocks at Meridiani Planum Mars. Earth Planet. Sci. Lett. 240, 149–167 (2005). 24. Tan, J., Lewis, J. M. T. & Sephton, M. A. The fate of lipid biosignatures in a mars-analogue sulfur stream. Sci. Rep. 8, 7586. https:// doi. org/ 10. 1038/ s4159801825752-7 (2018). 25. Preston, L. J. et al. The preservation and degradation of filamentous bacteria and biomolecules within iron oxide deposits at Rio Tinto Spain. Null 9, 233–249. https:// doi. org/ 10. 1111/j. 14724669. 2011. 00275.x (2011). 26. Colín-García, M. et al. Detection of peptidic sequences in the ancient acidic sediments of Río Tinto, Spain. Origins Life Evolut. Biospheres 41, 523–527. https:// doi. org/ 10. 1007/ s110840119258-x (2011). 27. Essalhi, M. et al. A case study of the internal structures of gossans and weathering processes in the Iberian Pyrite Belt using magnetic fabrics and paleomagnetic dating. Miner. Deposita 46, 981–999. https:// doi. org/ 10. 1007/ s001260110361-8 (2011). 28. Leistel, J. M. et al. The volcanic-hosted massive sulphide deposits of the Iberian Pyrite Belt. Miner. Deposita 33, 2–30 (1998). 29. Zachos, J., Pagani, M., Sloan, L., Thomas, E. & Billups, K. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science 292, 686. https:// doi. org/ 10. 1126/ scien ce. 10594 12 (2001). 30. Sand, W., Gehrke, T., Jozsa, P.-G. & Schippers, A. (Bio)chemistry of bacterial leaching–direct vs indirect bioleaching. Hydrometallurgy 59, 159–175 (2001). 31. Barrón, V., Torrent, J. & Greenwood, J. P. Transformation of jarosite to hematite in simulated Martian brines. Earth Planet. Sci. Lett. 251, 380–385. https:// doi. org/ 10. 1016/j. epsl. 2006. 09. 022 (2006). 32. Tosca, N. J., McLennan, S. M., Dyar, M. D., Sklute, E. C. & Michel, F. M. Fe oxidation processes at Meridiani Planum and implications for secondary Fe mineralogy on Mars. J. Geophys. Res. Planets https:// doi. org/ 10. 1029/ 2007J E0030 19 (2008). 33. Fernández-Remolar, D. C., Banerjee, N., Gómez-Ortiz, D., Izawa, M. & Amils, R. A mineralogical archive of the biogeochemical sulfur cycle preserved in the subsurface of the Río Tinto system. Am. Miner. 103, 394–411 (2018). 34. Hertkorn, N., Harir, M. & Schmitt-Kopplin, P. Nontarget analysis of Murchison soluble organic matter by high-field NMR spectroscopy and FTICR mass spectrometry. Magn. Reson. Chem. 53, 754–768. https:// doi. org/ 10. 1002/ mrc. 4249 (2015). 35. Schmitt-Kopplin, P. et al. High molecular diversity of extraterrestrial organic matter in Murchison meteorite revealed 40 years after its fall. Proc. Nat. Acad. Sci. 107, 2763. https:// doi. org/ 10. 1073/ pnas. 09121 57107 (2010). 36. Fernández-Remolar, D. C. et al. Molecular preservation in haliteand perchlorate-rich hypersaline subsurface deposits in the Salar Grande basin (Atacama Desert, Chile): Implications for the search for molecular biomarkers on Mars. J. Geophys. Res. Biogeosci. 118, 922–939. https:// doi. org/ 10. 1002/ jgrg. 20059 (2013). 37. Handle, F. et al. Tracking aging of bitumen and its saturate, aromatic, resin, and asphaltene fractions using high-field fourier transform ion cyclotron resonance mass spectrometry. Energy Fuels 31, 4771–4779. https:// doi. org/ 10. 1021/ acs. energ yfuels. 6b033 96 (2017). 38. McCarren, J. et al. Microbial community transcriptomes reveal microbes and metabolic pathways associated with dissolved organic matter turnover in the sea. Proc. Nat. Acad. Sci. 107, 16420. https:// doi. org/ 10. 1073/ pnas. 10107 32107 (2010). 39. Landa, M., Blain, S., Christaki, U., Monchy, S. & Obernosterer, I. Shifts in bacterial community composition associated with increased carbon cycling in a mosaic of phytoplankton blooms. ISME J. 10, 39–50. https:// doi. org/ 10. 1038/ ism ej . 2015. 105 (2016). 40. D’Andrilli, J., Smith, H. J., Dieser, M. & Foreman, C. M. Climate driven carbon and microbial signatures through the last ice age. Geochem. Perspect. Lett. 4, 29–34. https:// doi. org/ 10. 7185/ geoch emlet. 1732 (2017). 41. Burdige, D. J. Preservation of organic matter in marine sediments: Controls, mechanisms, and an imbalance in sediment organic carbon budgets?. Chem. Rev. 107, 467–485. https:// doi. org/ 10. 1021/ cr050 347q (2007). 42. Gelin, F. et al. Distribution of aliphatic, nonhydrolyzable biopolymers in marine microalgae. Org. Geochem. 30, 147–159. https:// doi. org/ 10. 1016/ S01466380(98) 00206-X (1999). 43. Eglinton, G. & Hamilton, R. J. Leaf epicuticular waxes. Science 156, 1322. https:// doi. org/ 10. 1126/ scien ce. 156. 3780. 1322 (1967). 44. Ficken, K. J., Li, B., Swain, D. L. & Eglinton, G. An n-alkane proxy for the sedimentary input of submerged/floating freshwater aquatic macrophytes. Org. Geochem. 31, 745–749. https:// doi. org/ 10. 1016/ S01466380(00) 00081-4 (2000). 21 Vol.:(0123456789) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ 45. Grimalt, J. & Albaigés, J. Sources and occurrence of C12-C22n-alkane distributions with even carbon-number preference in sedimentary environments. Geochim. Cosmochim. Acta 51, 1379–1384. https:// doi. org/ 10. 1016/ 00167037(87) 90322-X (1987). 46. Didyk, B. M., Simoneit, B. R. T., Brassell, S. C. & Eglinton, G. Organic geochemical indicators of Paleo-environmental conditions of sedimentation. Nature 272, 216–222 (1978). 47. Peters, K. E., Walters, C. C. & Moldowan, J. M. The Biomarker Guide—Part II—Biomarkers and Isotopes in Petroleum Exploration and Earth History. (Cambridge University Press, 2005). 48. O’Leary, W. M. in Microbial Composition: Carbohydrates, Lipids, and Minerals Vol. IV CRC Handbook of Microbiology (eds A I Laskin & H A Lechevalier) 391–434 (CRC Press, 1982). 49. Feng, X. & Simpson, M. J. The distribution and degradation of biomarkers in Alberta grassland soil profiles. Org. Geochem. 38, 1558–1570. https:// doi. org/ 10. 1016/j. orgge ochem. 2007. 05. 001 (2007). 50. Goad, J. & Akishisa, T. Analysis of Sterols. (Blackie Academic & Professional, 1997). 51. Volkman, J. K. A review of sterol markers for marine and terrigenous organic matter. Org. Geochem. 9, 83–99. https:// doi. org/ 10. 1016/ 01466380(86) 90089-6 (1986). 52. Weete, J. D., Abril, M. & Blackwell, M. Phylogenetic distribution of fungal sterols. PLoS ONE 5, e10899. https:// do i. org/ 10. 1371/ journ al. pone. 00108 99 (2010). 53. Ahlgren, G., Gustafsson, I.-B. & Boberg, M. Fatty acid content and chemical composition of freshwater microalgae1. J. Phycol. 28, 37–50. https:// doi. org/ 10. 1111/j. 00223646. 1992. 00037.x (1992). 54. Kaneda, T. Isoand anteiso-fatty acids in bacteria: Biosynthesis, function, and taxonomic significance. Microbiol. Rev. 55, 288–302 (1991). 55. Langworthy, T. A., Holzer, G., Zeikus, J. G. & Tornabene, T. G. Isoand anteiso-branched glycerol diethers of the thermophilic anaerobe thermodesulfotobacterium commune. Syst. Appl. Microbiol. 4, 1–17. https:// doi. org/ 10. 1016/ S07232020(83) 80029-0 (1983). 56. Doi, S., Hashimoto, Y., Tomita, C., Kumano, T. & Kobayashi, M. Discovery of piperonal-converting oxidase involved in the metabolism of a botanical aromatic aldehyde. Sci. Rep. 6, 38021. https:// doi. org/ 10. 1038/ srep3 8021 (2016). 57. Zhou, X. & Wu, F. Vanillic acid changed cucumber (Cucumis sativus L.) seedling rhizosphere total bacterial, Pseudomonas and Bacillus spp. communities. Sci. Rep. 8, 4929. https:// doi. org/ 10. 1038/ s4159801823406-2 (2018). 58. Hwang, I. H., Swenson, D. C., Gloer, J. B. & Wicklow, D. T. Pestaloporonins: Caryophyllene-derived sesquiterpenoids from a fungicolous isolate of pestalotiopsis sp. Org. Lett. 17, 4284–4287. https:// doi. org/ 10. 1021/ acs. orgle tt. 5b020 80 (2015). 59. Gonçalves Tavares, D., Viana Lessa Barbosa, B., Lopes Ferreira, B., Ferreira Duarte, R. & Gomes Cardoso, W. Antioxidant activity and phenolic compounds of the extract from pigment-producing fungi isolated from Brazilian caves. Biocatalysis Agric. Biotechnol. 16, 148–154. https:// doi. org/ 10. 1016/j. bcab. 2018. 07. 031 (2018). 60. Leefmann, T. et al. Biomarker imaging of single diatom cells in a microbial mat using time-of-flight secondary ion mass spectrometry (ToF-SIMS). Org. Geochem. 57, 23–33. https:// doi. org/ 10. 1016/j. orgge ochem. 2013. 01. 005 (2013). 61. Perkins, M. C. et al. Surface morphology and chemistry of Prunus laurocerasus L. leaves: a study using X-ray photoelectron spectroscopy, time-of-flight secondary-ion mass spectrometry, atomic-force microscopy and scanning-electron microscopy. Planta 221, 123–134. https:// doi. org/ 10. 1007/ s004250041417-0 (2005). 62. Vichi, S., Cortés-Francisco, N., Romero, A. & Caixach, J. Direct chemical profiling of olive (Olea europaea) fruit epicuticular waxes by direct electrospray-ultrahigh resolution mass spectrometry. J. Mass Spectrom. 50, 558–566. https:// doi. org/ 10. 1002/ jms. 3562 (2015). 63. Weete, J. D. in Chemistry and biochemistry of natural waxes Vibrational spectra and structure (ed P E Kolattukudy) 349–418 (Elsevier Scientific, 1976). 64. Huang, X., Xue, J. & Guo, S. Long chain n-alkanes and their carbon isotopes in lichen species from western Hubei Province: Implication for geological records. Front. Earth Sci. 6, 95–100. https:// doi. org/ 10. 1007/ s117070120300-8 (2012). 65. Stephan, T. TOF-SIMS in cosmochemistry. Planet. Space Sci. 49, 859–906. https:// doi. org/ 10. 1016/ S00320633(01) 00037-X (2001). 66. Stephan, T., Jessberger, E. K., Heiss, C. H. & Rost, D. TOF-SIMS analysis of polycyclic aromatic hydrocarbons in Allan Hills 84001. Meteorit. Planet. Sci. 38, 109–116. https:// doi. org/ 10. 1111/j. 19455100. 2003. tb010 49.x (2003). 67. Gadd, G. M. in Encyclopedia of Geobiology (eds Joachim Reitner & Volker Thiel) 416–432 (Springer Netherlands, 2011). 68. López-Mondéjar, R., Zühlke, D., Becher, D., Riedel, K. & Baldrian, P. Cellulose and hemicellulose decomposition by forest soil bacteria proceeds by the action of structurally variable enzymatic systems. Sci. Rep. 6, 25279. https:// doi. org/ 10. 1038/ srep2 5279 (2016). 69. Killops, S. D. & Massoud, M. S. Polycyclic aromatic hydrocarbons of pyrolytic origin in ancient sediments: evidence for Jurassic vegetation fires. Org. Geochem. 18, 1–7. https:// doi. org/ 10. 1016/ 01466380(92) 90137-M (1992). 70. Heim, C., Sjövall, P., Lausmaa, J., Leefmann, T. & Thiel, V. Spectral characterisation of eight glycerolipids and their detection in natural samples using time-of-flight secondary ion mass spectrometry. Rapid Commun. Mass Spectrom. 23, 2741–2753. https:// doi. org/ 10. 1002/ rcm. 4183 (2009). 71. Fernández-Remolar, D. C. et al. Formation of iron-rich shelled structures by microbial communities. J. Geophys. Res. Biogeosci. 120, 147–168. https:// doi. org/ 10. 1002/ 2014J G0027 45 (2015). 72. Harwood, J. L. & Russell, N. J. Lipids in plants and microbes. 162 (George Allen & Unwin, 1984). 73. Sánchez-García, L. et al. Molecular biomarkers in the subsurface of the Salar Grande (Atacama, Chile) evaporitic deposits. Biogeochemistry 140, 31–52. https:// doi. org/ 10. 1007/ s105330180477-3 (2018). 74. Deelai, S., Suetrong, S., Damrianant, S., Unagul, P. & Sakkayawong, N. Isolation and identification of native lower fungi for polyunsaturated fatty acid (PUFA) production in Thailand, and the effect of carbon and nitrogen sources on growth and production. Afr. J. Biotech. 14, 1449–1460 (2015). 75. Imai, H., Ohnishi, M., Kinoshita, M., Kojima, M. & Ito, S. Structure and distribution of cerebroside containing unsaturated hydroxy fatty acids in plant leaves. Biosci. Biotechnol. Biochem. 59, 1309–1313. https:// doi. org/ 10. 1271/ bbb. 59. 1309 (1995). 76. Killops, S. D. & Killops, V. Introduction to organic chemistry (Blackwell Science, 2005). 77. Van Dyk, M. S., Kock, J. L. F. & Botha, A. Hydroxy long-chain fatty acids in fungi. World J. Microbiol. Biotechnol. 10, 495–504. https:// doi. org/ 10. 1007/ BF003 67653 (1994). 78. Kinderlerer, J. L. Biotransformations of lauric acid oils by fungi. J. Chem. Technol. Biotechnol. 55, 400–402. https:// doi. org/ 10. 1002/ jctb. 28055 0419 (1992). 79. Vaezian, B., Anderton, C. R. & Kraft, M. L. Discriminating and imaging different phosphatidylcholine species within phaseseparated model membranes by principal component analysis of TOF-secondary ion mass spectrometry images. Anal. Chem. 11, 10006–10014 (2010). 80. Singh, A. & Del Poeta, M. Sphingolipidomics: An important mechanistic tool for studying fungal pathogens. Front. Microbiol. 7, 501 (2016). 81. Qu, F., Zhang, H., Zhang, M. & Hu, P. Sphingolipidomic profiling of rat serum by UPLC-Q-TOF-MS: Application to rheumatoid arthritis study. Molecules https:// doi. org/ 10. 3390/ molec ules2 30613 24 (2018). 82. Barreto-Bergter, E., Sassaki, G. L. & de Souza, L. M. Structural analysis of fungal cerebrosides. Front. Microbiol. 2, 239–239. https:// doi. org/ 10. 3389/ fmicb. 2011. 00239 (2011). 22 Vol:.(1234567890) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ 83. Siljeström, S., Parenteau, M. N., Jahnke, L. L. & Cady, S. L. A comparative ToF-SIMS and GC–MS analysis of phototrophic communities collected from an alkaline silica-depositing hot spring. Org. Geochem. 109, 14–30. https:// doi. org/ 10. 1016/j. orgge ochem. 2017. 03. 009 (2017). 84. Kawecki, M. & Bernard, L. Database of proteinogenic amino acid reference spectra for Bismuth-cluster ToF-SIMS I. Negative polarity. Surface Sci. Spectra 25, 015001. https:// doi. org/ 10. 1116/1. 50249 87 (2018). 85. Kawecki, M. & Bernard, L. Database of proteinogenic amino acid reference spectra for Bismuth-cluster ToF-SIMS II. Positive polarity. Surface Sci. Spectra 25, 015002. https:// doi. org/ 10. 1116/1. 50249 88 (2018). 86. Leefmann, T. et al. Spectral characterization of ten cyclic lipids using time-of-flight secondary ion mass spectrometry. Rapid Commun. Mass Spectrom. 27, 565–581. https:// doi. org/ 10. 1002/ rcm. 6483 (2013). 87. Banoub, J. et al. A critique on the structural analysis of lignins and application of novel tandem mass spectrometric strategies to determine lignin sequencing. J. Mass Spectrom. 50, 5–48. https:// doi. org/ 10. 1002/ jms. 3541 (2015). 88. Janusz, G. et al. Lignin degradation: microorganisms, enzymes involved, genomes analysis and evolution. FEMS Microbiol. Rev. 41, 941–962. https:// doi. org/ 10. 1093/ femsre/ fux049 (2017). 89. Godden, B., Ball, A. S., Helvenstein, P., Mccarthy, A. J. & Penninckx, M. J. Towards elucidation of the lignin degradation pathway in actinomycetes. Microbiology 138, 2441–2448. https:// doi. org/ 10. 1099/ 00221 287138112441 (1992). 90. Eglinton, G., Logan, G. A., Ambler, R. P., Boon, J. J. & Perizonius, W. R. K. Molecular preservation [and discussion]. Philosophical Trans. Biol. Sci. 333, 315–328 (1991). 91. Volkman, J. K. Sterols in microorganisms. Appl. Microbiol. Biotechnol. 60, 495–506 (2003). 92. Escudero, C., Vera, M., Oggerin, M. & Amils, R. Active microbial biofilms in deep poor porous continental subsurface rocks. Sci. Rep. 8, 1538. https:// doi. org/ 10. 1038/ s4159801819903-z (2018). 93. Del Poeta, M., Nimrichter, L., Rodrigues, M. L. & Luberto, C. Synthesis and biological properties of fungal glucosylceramide. PLoS Pathog 10, e1003832. https:// doi. org/ 10. 1371/ journ al. ppat. 10038 32 (2014). 94. Mitchell, K. F., Zarnowski, R. & Andes, D. R. Fungal super glue: The biofilm matrix and its composition, assembly, and functions. PLoS Pathog 12, e1005828. https:// doi. org/ 10. 1371/ journ al. ppat. 10058 28 (2016). 95. Fernández-Remolar, D., Rodríguez, N., Gómez, F. & Amils, R. The geological record of an acidic environment driven by iron hydrochemistry: The Tinto River system. J. Geophys. Res. 10, 8. https:// doi. org/ 10. 1029/ 2002J E0019 18 (2003). 96. Briggs, D. E. G. The role of decay and mineralization in the preservation of soft-bodied fossils. Annu. Rev. Earth Planet. Sci. 31, 275–301. https:// doi. org/ 10. 1146/ annur ev. earth. 31. 100901. 144746 (2003). 97. Einsele, G. Sedimentary Basins. Evolution, facies and sedimentary budget. (Springer-Verlag, 1992). 98. Bird, L. J., Coleman, M. L. & Newman, D. K. Iron and Copper Act Synergistically To Delay Anaerobic Growth of Bacteria. Appl. Environ. Microbiol. 79, 3619. https:// doi. org/ 10. 1128/ AEM. 0394412 (2013). 99. Igiri, B. E. et al. Toxicity and bioremediation of heavy metals contaminated ecosystem from tannery wastewater: A review. J. Toxicol. 2018, 2568038. https:// doi. org/ 10. 1155/ 2018/ 25680 38 (2018). 100. Tosca, N. J., Knoll, A. H. & McLennan, S. M. Water activity and the challenge for life on early mars. Science 320, 1204. https:// doi. org/ 10. 1126/ scien ce. 11554 32 (2008). 101. Stevenson, A. et al. Is there a common water-activity limit for the three domains of life?. ISME J. 9, 1333–1351. https:// doi. org/ 10. 1038/ ismej. 2014. 219 (2015). 102. Klingelhöfer, G. et al. Jarosite and hematite at meridiani planum from opportunity’s mossbauer spectrometer. Science 306, 1740–1745 (2004). 103. Fernández-Remolar, D. et al. The Tinto River, an extreme acidic environment under control of iron, as an analog of the Terra Meridiani hematite site of Mars. Planet. Space Sci. 52, 239–248 (2004). 104. Syamaladevi, R. M. et al. Influence of water activity on thermal resistance of microorganisms in low-moisture foods: A review. Comprehens. Rev. Food Sci. Food Safety 15, 353–370. https:// doi. org/ 10. 1111/ 15414337. 12190 (2016). 105. Gómez, F., Aguilera, A. & Amils, R. Soluble ferric iron as an effective protective agent against UV radiation: Implications for early life. Icarus 191, 352–359. https:// doi. org/ 10. 1016/j. icarus. 2007. 04. 008 (2007). 106. Fernández-Remolar, D. C. et al. Preservation of underground microbial diversity in ancient subsurface deposits (> 6 Ma) of the Rio Tinto basement. Microorganisms https:// doi. org/ 10. 3390/ micro organ isms9 081592 (2021). 107. Fernández-Remolar, D., Harir, M., Carrizo, D., Schmitt-Kopplin, P. & Amils, R. Formation of shale-hosted massive sulfide deposits in the iberian pyrite belt (Tharsis, Spain). J. Geophys. Res. 123, 1017–1040 (2018). 108. Sodhi, R. N. S. Time-of-flight secondary ion mass spectrometry (TOF-SIMS)—versatility in chemical and imaging surface analysis. Analyst 129, 483–487. https:// doi. org/ 10. 1039/ B4026 07C (2004). 109. Patiny, L. & Borel, A. ChemCalc: A building block for tomorrow’s chemical infrastructure. J. Chem. Inf. Model. 53, 1223–1228. https:// doi. org/ 10. 1021/ ci300 563h (2013). 110. Strohalm, M., Kavan, D., Novák, P., Volný, M. & Havlíček, V. mMass 3: A cross-platform software environment for precise analysis of mass spectrometric data. Anal. Chem. 82, 4648–4651. https:// doi. org/ 10. 1021/ ac100 818g (2010). 111. Pence, H. E. & Williams, A. ChemSpider: An online chemical information resource. J. Chem. Educ. 87, 1123–1124. https:// doi. org/ 10. 1021/ ed100 697w (2010). 112. Smith, C. A. et al. METLIN: A metabolite mass spectral database. Ther Drug Monit 27, 747–751 (2005). 113. Fahy, E., Sud, M., Cotter, D. & Subramaniam, S. LIPID MAPS online tools for lipid research. Nucleic Acids Res. 35, W606–W612. https:// doi. org/ 10. 1093/ nar/ gkm324 (2007). Author contributions D.C.F.R., D.G.O., T.H. and R.A. performed fieldwork and sample collection. D.C.F.R., D.C., M.H., and L.S.G. conceived and performed analyses. D.C.F.R., D.C., M.H., P.S.K., and P.M. analyzed the data. The paper was written by D.C.F.R., D.C., and M.H. with contribution with the other Co-As. Funding Funding was provided by Investissements d’Avenir program (Grant No. ANR-15-IDEX-02), European Regional Development Fund (Grant No. RYC2018-023943-I), Horizon 2020 Framework Programme (Grant No. RYC-2014-19446). Competing interests The authors declare no competing interests. Additional information Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s4159802100730-8. 23 Vol.:(0123456789) Scientific Reports | (2021) 11:21543 | https://doi.org/10.1038/s41598-021-00730-8 www.nature.com/scientificreports/ Correspondence and requests for materials should be addressed to D.C.F.-R. 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:// creat iveco mmons. org/ licen ses/ by/4. 0/. © The Author(s) 2021, corrected publication 2021