Ecological successions throughout the desiccation of Tirez lagoon (Spain) as an astrobiological time-analog for wet-to-dry transitions on Mars
The ecological baseline in TirezThe geology and the climate of the Tirez region favored the generation and maintenance of a type of hypersaline habitat characterized by extreme seasonality: the sulfate-chlorine waters, with sodium and magnesium cations, showed significant seasonal variations15. The alkaline pH, the low oxidant value for the redox potential of the water column and the highly reduced sediments imposed extreme conditions (see Table 1 and Supplementary Information for details). This extreme seasonality requires to define a valid representative ecological baseline to compare the ecology of the lagoon between 2002 and 2021 and, in this way, set the basis to proposing our model of ecological succession with increasing dryness as a “time-analog” for early Mars. Taxonomic data from 2002 is a snapshot of the community during one season, so we include in our discussion the results presented by Montoya et al. (2013) from a sample campaign carried out in 2005, because they16 analyzed both water and sediment and during both the wet and dry seasons.We consider here only the results obtained by Montoya et al.16 by gene cloning, since those obtained by isolation and sequencing are not comparable. At the level of large groups, no major seasonal differences were observed: Pseudomonadota, followed by Bacteroidetes, were the dominant phyla, in both water and sediments, and both in the dry and the rainy seasons; although Alphaproteobacteria was the dominant class in water, while Gammaproteobacteria was dominant in sediments (in both dry and rainy seasons). With respect to the archaeal domain, all the identified sequences were affiliated to Halobacteriales order, mainly Halorubrum (water) and Halobacterium (sediment), both within Halobacteriaceae family. We can consider these results presented in Montoya et al.16 as the “ecological baseline” for Tirez, however taken with a grain of salt, because only 43 bacterial and 35 archaeal sequences, including rainy and dry seasons and water and sediments, were considered for analysis.Prokaryotic diversity in 2002As can be expected for an extreme environment, the bacterial diversity detected in 2002 was low, although we cannot exclude the possibility that this may reflect the limitation of DNA sequencing techniques at the time. 59% of the obtained clones in the then-wet sediments corresponded to the Malaciobacter genus. Malaciobacter (previously named17 Arcobacter) is an aerotolerant Epsilonproteobacteria. Species within this genus are moderately halophilic, e.g., M. halophilus, capable to grow in up to 4% NaCl. Even though the role that Malaciobacter can play in the environment is not known, it seems to thrive in aquatic systems, like sewage, with a high organic matter content17: e.g., M. canalis, M. cloacae, or M. defluvii.After Malaciobacter-like clones, the next most numerous group belongs to the phylum Bacillota (27% of the sequenced clones; Table 2). Under stressful environmental conditions, members of the genus Virgibacillus produce endospores, a very useful property in an extreme and variable environment (ionic strength, temperature, light intensity), easy to compare with early Mars. Endospores facilitate species survival, allowing them to overcome drastic negative changes, like dry periods, and to germinate when the conditions are favorable again. The closest identified species was the halotolerant V. halodenitrificans, but with low homology, not far from other halotolerant (e.g., V. dokdonensis) or halophilic (e.g., V. marismortui) species within the same genus. The other Gram-positive clones belong to the order Clostridiales. These clones cluster in two taxonomic units related with the strictly anaerobic genus Tissierella.Despite the abundance of Pseudomonadota, their biodiversity was very low, reduced to only two genera within the Epsilon- and Delta-proteobacteria. Six sequences affiliated to Deltaproteobacteria, and clustered in one OTU (salB38, similarity 96.6% with Desulfotignum), were retrieved. Its presence in anaerobic media rich in sulfates (Table 1) seems reasonable. In fact, sulfate-reducing activity was detected using a specific enrichment assay.Finally, one taxon belonging to the phylum Spirochaetota (previously named Spirochaetes) was identified. The presence of Spirochaetota in this system is not strange because members of the genus Spirochaeta are very often found in mud and anaerobic marine environments rich in sulfates18. Moreover, the closest species to SalB63, although with a low similarity of 87%, was Spirochaeta bajacaliforniensis, a spirochete isolated19 from a microbial mat in Laguna Figueroa (Baja California), an extensive hypersaline lagoon with high gypsum content, very similar, although much bigger, than Tirez lagoon.The diversity within the domain Archaea was very low in 2002. The phylogenetic analysis of 96 clones indicate that they correspond to one specie belonging to the obligate halophile genus Methanohalophilus. Their high similarity (99.3%) with several species of Methanohalophilus, such as M. portucalensis (isolated from sediments of a solar saltern in Portugal), M. mahii (isolated from sediments of the Great Salt Lake), or M. halophilus (isolated20 from a cyanobacterial mat at Hamelin Pool, Australia), makes impossible its adscription to any particular species level. Methanohalophilus is strictly methylotrophic, which is consistent with this environment, given that the methylotrophic methanogenesis pathway, non-competitive at low-salt conditions, is predominant at high saline concentrations21. We further confirmed methanogenic activity in Tirez by the measurement of methane by gas chromatography in enrichment cultures.Prokaryotic diversity in context of other studies between 2002 and 2021It was challenging to establish a timeline for the succession of the populations involved, because the scarcity of data harvested and published so far from Tirez. However, combining our results with the few data available in Montoya et al.16 and Preston et al.22 on samplings carried out on 2005 and 2017, respectively, we can see a clear predominance of the phylum Pseudomonadota: Epsilonproteobacteria, i.e. Arcobacter-like, and Deltaproteobacteria, mainly sulfate-reducing bacteria (this work, sampling 2002), and Gammaproteobacteria16 when Tirez maintained a water film, to eventually a final predominance of Gammaproteobacteria, e.g. Chromatiales and Pseudomonadales, in the dry Tirez (this work, 2020 sampling). The Bacillales order has remained widely represented both in the wet and dry Tirez.Regarding the archaeal domain, the few references available (Refs.16,22; this work) confirm that the members of the Halobacteriaceae family are well adapted to both the humid and dry ecosystems of Tirez, being predominant in both conditions. Preston et al.22 found that the second most abundant group of archaea in the dry sediments of Tirez was the Methermicoccaceae family, within the Methanosarcinales order, Methanomicrobia class. Taking into account the results obtained in the dry Tirez (Preston et al.22; and this work, sample 2020), the methanogenic archaea have decreased drastically through time, probably due to salt stress and the competition with sulfate-reducing bacteria.Prokaryotic diversity in 2021From a metabolic point of view, most of the bacteria present today in the sediment are chemoorganotrophs, anaerobes, and halophilic or halotolerant. Scarce information is available about the predominant OTU, Candidate Division OP1. The OP1 division was one of the main bacterial phyla in a sulfur-rich sample in the deepest analyzed samples from the Red Sea sediments under brine pools23. In addition, the phylogenetically related Candidate division KB1 has been observed in deep-sea hypersaline anoxic basins at Orca Basin (Gulf of Mexico), and other hypersaline environments24. Eight of the nine genera identified show coverage greater than 1% of the sequences: i.e., Rubinisphaera, Halothiobacillus, Thiohalophilus, Anaerobacillus/Halolactibacillus, Halomonas, Halothermothrix, and Aliifodinibius are halophilic or halotolerant genera13,25.Regarding archaea, our analyses reveal archaeal groups that seem to thrive in sediments from extreme environments, e.g., marine brine pools/deep water anoxic basins or hypersaline lakes. The most abundant OTU, Thermoplasmata KTK4A, was found prominent and active in the sediment of Lake Strawbridge, a hypersaline lake in Western Australia26, and in soda-saline lakes in China27. The creation of a Candidatus Haloplasmatales, a novel order to include KTK4A-related Thermoplasmata, has been proposed27. On the other hand, both in the aforementioned soda-saline lakes in China27 and in a sulfur-rich section of the sediments from below the Red Sea brine pools23, retrieved sequences were assigned to Marine Benthic Groups B, D, and E. Finally, in the section of nitrogen-rich sediments from the aforementioned Red Sea brine pools, the unclassified lineage ST-12K10A represented the most abundant archaeal group. In the Tirez Lagoon sediment after desiccation, all Methanomicrobia readings belonged to this group.The significance and implications of an ecosystem characterized in 2021 by high diversity, high inequality, and lack of isolated representatives, resides in that Tirez is today an ecosystem in which many (most) of the species/OTUs present are dormant, and they do not play any metabolic role. Hence the high percentage of raretons, greater than 80% for both bacteria and archaea, which are actually present in the lagoon but with only one or two copies each. Only those species adapted to the conditions imposed by the extreme environment are able to actually thrive, and consequently only a few species carry out all the metabolic activity. We conclude that the microbiota in Tirez today represents an ecosystem with a high resilience capacity in the face of environmental changes that may occur.We want to clearly highlight that the technique available in 2002 to study the microbiota of the Tirez lagoon only allowed to obtain a low-resolution image, but that was the state-of-the-art procedure at the time, and the Tirez lagoon cannot be sampled again with the conditions back in 2002, which no longer exist and are not expected to return. Although we have kept in storage several samples of water and sediment from the 2002 Tirez lagoon, it is reasonable to assume that those laboratory microcosms would have chemically and microbiologically changed during the last 20 years, and as such no longer represent reliable replicas of the original lagoon, so we cannot use them for the purposes of this work. Therefore, we are aware that any comparisons of the 2002 laboratory results with the much more robust results obtained by Illumina in 2021 need to be taken with a grain of salt. With all the precautions required, in a high-level, first-order comparison, the most noticeable difference between 2002 and 2021 is a drastic change in the microbial Tirez population. Only some OTUs within Bacillales (Virgibacillus/Anaerobacillus), sulfate-reducing Deltaproteobacteria (Desulfotignum/Desulfobacteraceae-Desulfovibrio), and Spirochaetes are shared among the 2002 and 2021 samples. This comparison is enough for the purposes of this work, as we are interested in the evolution of the lagoon system as a whole to establish a “time-analog” with the wet-to-dry transition on early Mars, and not in the particular outcome of each and every OTU in Tirez. With the results at hand, we conclude that, since 2002, the lacustrine microbiota has shifted to one more adapted to the extreme conditions in the dry sediments, derived from the gradual and persistent desiccation concluding ca. 7 years ago (i.e., completely desiccated in 2015), such as lack of light, absence of oxygen, and lack of water availability. This shift has likely been triggered because organisms that were originally in the lagoon but at low abundance in 2002 became dominant as they were better adapted to desiccation, and because the incoming of new microorganisms transported by birds or wind28.Lipid biomarkers analysis of the desiccated lake sedimentsThe analysis of cell membrane-derived lipid compounds on the dry lake sediments at present allow to provide another perspective of the microbial communities inhabiting the Tirez lagoon, by contributing additional information about the ecosystem and depositional environment. It is important to note that, analyzing only the 2021 lake sediments, we cannot differentiate between lipidic biomarkers of the microorganisms inhabiting Tirez in 2002 and before from those left behind by the microorganisms living in the dried sediments today. Instead, the analyses of lipid biomarkers provide clues about the different microorganisms that have populated Tirez through time, including both older communities inhabiting the former aqueous system and also younger communities better adapted to the present dry conditions. Thus, the lipid biomarkers analysis can be considered as a time-integrative record of the microbial community inhabiting Tirez during the last decades.Based on the molecular distribution of lipid biomarkers, the presence of gram-positive bacteria was inferred from the relative abundance of the monounsaturated alkanoic acid C18:1[ω9], or iso/anteiso pairs of alkanoic acids from 12 to 17 carbons29 with dominance of i/a-C15:0 and i/a-C17:0 (Fig. 3B). In contrast, generally ubiquitous alkanoic acids such as C16:1[ω7], C18:1[ω7], or C18:2[ω6] suggested a provenance rather related to gram-negative bacteria30. The combined detection of the i/a-C15:0 and i/a-C17:0 acids, with dominance of the iso over the anteiso congeners, together with other biomarkers such as the mid-chain branched 10Me16:0, the monounsaturated C17:1, or the cyclopropyl Cy17:0 and Cy19:0 acids, may be associated with a community of SRB31 in today´s dry sediments of Tirez. Specifically, most of those alkanoic acids have been found in a variety of Deltaproteobacteria and/or Bacteroidota (previously named Bacteroidetes). The presence of archaea was deduced from the detection of prominent peaks of archaeol in the polar fraction32 (Fig. 3C), as well as squalene and relatives (dihydrosqualene and tetrahydrosqualene) in the apolar fraction33 (Fig. 3A). Squalene and a variety of unsaturated derivatives are present in the neutral lipid fractions of many archaea with high abundances in saline lakes34. The relative abundance of autotrophs over heterotrophs35 can be estimated by the ratio of the autotrophically-related pristane and phytane over the both autotrophically- and heterotrophically-produced n-C17 and n-C18 alkanes ([Pr + Ph]/[n-C17 + n-C18]). A ratio of 0.56 in the Tirez sediments suggest the presence of a relevant proportion of heterotrophs in the ancient lacustrine system.Furthermore, the lipid biomarkers analysis was able to detect compounds specific of additional microbial sources, such as cyanobacteria36 (n-C17, C17:1, or 7Me-C15 and 7Me-C17), microalgae and/or diatoms (phytosterols37; or C20:5, and C22:6 alkanoic acids30), and other photoautotrophs (phytol and potentially degradative compounds such as pristane and phytane31). A relatively higher preservation of the cell-membrane remnants (i.e., lipids) compared to the DNA-composing nucleic acids may contribute to explain the lack of detection of cyanobacteria, diatoms and microalgae, and other phototrophs by DNA analysis (a deficit in our results shared with Montoya et al.16, and Preston et al.22). Although abundant in higher plants38, sterols such as those detected here (i.e., the sterols campesterol, stigmasterol, and β-sitosterol, as well as ergosterol) are also major sterols in some microalgal classes37 (such as Bacillariophyceae, Chrysophyceae, Euglenophyceae, Eustigmatophyceae, Raphidophyceae, Xanthophyceae, and Chlorophyceae), cyanobacteria (β-sitosterol), and fungi (ergosterol39).The carbon isotopic composition of lipid biomarkers provides a rapid screening of the carbon metabolism in a system, by recognizing the principal carbon fixation pathways used by autotrophs. The range of δ13C values measured in the Tirez sediments (from − 33.9 to − 16.1‰) denotes a mixed use of different carbon assimilation pathways, involving mostly the reductive pentose phosphate (a.k.a. Calvin–Benson–Bassham or just Calvin) cycle (from − 19 to − 30‰), and in lesser extent the reductive acetyl-CoA (a.k.a. Wood–Ljungdahl) pathway (from − 28 to − 44‰), and/or the reverse tricarboxylic acid (rTCA) cycle (from − 12 to − 21‰).The lipids synthesized by microorganisms using the Calvin or reductive acetyl-CoA pathway are typically depleted relative to the bulk biomass, particularly those produced via de latter pathway. In the dry Tirez sediments, the majority of the lipid compounds are more depleted in 13C than the bulk biomass (Fig. 4). In particular, the branched alkane DiMeC18 (Fig. 4A) and the SRB-indicative 10Me16:0 acid (Fig. 4B) showed the most depleted δ13C values and suggested the use of the reductive acetyl-CoA pathway. The rest of lipid compounds showed isotopic signatures (from − 16.1 to − 31.4‰) compatible with the prevalence of the Calvin pathway. These values may directly reflect the autotrophic activity of microorganisms fixing carbon via the Calvin cycle or heterotrophic activity of microorganisms growing on their remnants. Thus, the saturated and linear alkyl chains of lipids (i.e., n-alkanes, n-alkanoic acids, and n-alkanols) showing the most negative δ13C values (e.g., alkanes n-C17 and C17:1; or acid C18:1[ω7]) reflect prokaryotic sources of Calvin-users autotrophs (e.g., cyanobacteria or purple sulfur bacteria), while the rest of compounds with slightly less negative δ13C values instead stem from the autotrophic activity of eukaryotes also users of the Calvin cycle (unsaturated fatty acids and sterols) or from the metabolism of heterotrophs such as SRB (iso/anteiso-, other branched, and cyclopropyl fatty acids) and haloarchaea (isoprenoids, phytanol, and archaeol). All in all, the compound-specific isotope composition of the dry sediments in the today´s Tirez lagoon may indirectly reflect the prevailing autotrophic mechanisms in the present lacustrine system of Tirez, by showing isotopic signatures of secondary lipids similar to their carbon source40.In addition, the use of a number of lipid molecular ratios or proxies allow further characterization of the lacustrine ecosystem and depositional environment. For example, the average chain length of the n-alkanes (24.1) suggests a relevant presence of eukaryotic biomass in the lacustrine sediments, since long-chained alkanes ( > C20) are known to originate from epicuticular leaf waxes in higher plants41. Highlighting the relevance of eukaryotes and their ecological roles is one of the major contributions of this work, because previous studies on the microbial ecology of hypersaline environments have been focused primarily on prokaryotes42.The proportion of odd n-alkanes of high molecular chain (i.e., n-C27, n-C29, and n-C31) over even n-alkanes of low molecular chain (i.e., n-C15, n-C17, and n-C19) provides an estimate of the relative abundance of terrigenous over aqueous biomass43, which in Tirez is TAR = 1.8. The Paq index may also be used to differentiate the proportion of terrigenous versus aquatic (emergent and submerged) plant biomass44. A Paq of 0.3 in the Tirez sediments from 2021 supported the relative abundance of land plants. Finally, the depositional environment in the lacustrine system of Tirez may be also characterized analyzing the ratio of pristane over phytane (Pr/Ph), which is higher than 1 when phytol degrades to pristane under oxic conditions45. Assuming that both isoprenoids in the Tirez sediments derived from phytol31, according to their similarly depleted δ13C (Fig. 4A), we can conclude that the sediments in the Tirez lagoon were deposited under predominantly oxic conditions (i.e., Pr/Ph ratio of 1.1).In summary, the lipid biomarkers study revealed useful information about the depositional environment and lacustrine ecosystem, including the presence of active or past autotrophic metabolisms involving prokaryotes (e.g., cyanobacteria and purple sulfur bacteria) and eukaryotes (plants, diatoms and other microalgae), as well as heterotrophic metabolisms of likely SRB and haloarchaea growing on Calvin-users exudates. These results are quite in agreement with the microbial community previously reported16 in sediments from the wet and dry seasons: abundant Gammaproteobacteria and Alphaproteobacteria, together with Algae and Cyanobacteria, dinoflagellates and filamentous fungi, Bacillota, Actinomycetota (previously named Actinomycetes), and a halophilic sulfate-reducing Deltaproteobacteria.Tirez as the first astrobiological “time-analog” for early MarsEarly Mars most likely had a diversity of environments in terms of pH, redox conditions, geochemistry, temperature, and so on. Field research in terrestrial analog environments contribute to understand the habitability of this diversity of environments on Mars in the past, because terrestrial analogues are places on Earth characterized by environmental, mineralogical, geomorphological, or geochemical conditions similar to those observed on present or past Mars9. Therefore, so far analogs have been referred to terrestrial locations closely similar to any of the geochemical environments that have been inferred on Mars, i.e., they are “site-analogs” that represent snapshots in time: one specific environmental condition at a very specific place and a very specific time. Because of this, each individual field analog site cannot be considered an adequate representation of the changing martian environmental conditions through time. Here we introduce the concept of astrobiological “time-analog”, referred to terrestrial analogs that may help understand environmental transitions and the related possible ecological successions on early Mars. In this sense, they should be “time-resolved analogs”: dynamic analog environments where we can analyze changes over time. To the best of our knowledge, this is the first study that looks at the environmental microbiology of a Mars astrobiological analog site over a significant and long period of change, and try to understand the ecological successions to put them in the context of martian environmental evolution.As Mars lost most of its surface water at the end of the Hesperian5,9,12, this wet-to-dry global transition can be considered the major environmental perturbation in the geological history of Mars, and therefore merits to be the first one to be assigned a “time-analog” for its better understanding and characterization. The drying of Mars was probably a stepwise process, characterized by multiple transitions between drier and wetter environments12,47, and therefore the seasonal fluctuations and eventual full desiccation of Tirez represent a suitable analog to better understand possible ecological transitions during the global desiccation of most of the Mars’s surface before the Amazonian (beginning 3.2 Ga).To introduce Tirez as the first Mars astrobiological “time-analog” of the wet-to-dry transition on early Mars, the objective of this study was threefold: first, we wanted to identify the dominant prokaryotic microorganisms in the active Tirez lagoon 20 years ago, a unique hypersaline ecosystem with an ionic composition different from that of marine environments, and therefore potentially analogous to ancient saline lacustrine environments on Mars during the Noachian and into the Hesperian46,47. Our results provide a preliminary basis to hypothesize how the microbial communities on the Noachian Mars could have developed in salty environments with dramatically fluctuating water availability. The requirement to deal with important variations in ionic strength and water availability, involving at times the complete evaporation of the water, could have represented additional constraints48 for microorganisms on early Mars.The second objective of this investigation was the identification of the microbial community inhabiting the desiccated Tirez sediments today, after all the water was lost, as a potential analog to desiccated basins on Mars at the end of the Hesperian1,3,4,47. Our results suggest that hypothetical early microbial communities on early Mars, living with relative abundance of liquid water during the Noachian, would have been forced to adapt to increasingly desiccating surface environments, characterized by extreme conditions derived from the persistent dryness and lack of water availability. Our investigation in Tirez suggest that hypothetical microorganisms at the end of the Hesperian would have needed to evolve strategies similar to those of microorganisms on Earth adapted to living at very low water activity49, to thrive in the progressively desiccating sediments.And the third objective of this investigation was the identification of the lipidic biomarkers left behind by the microbial communities in Tirez, as a guide to searching and identifying the potential leftovers of a hypothetical ancient biosphere on Mars. Lipids (i.e., fatty acids and other biosynthesized hydrocarbons) are structural components of cell membranes bearing recognized higher resistance to degradation relative to other biomolecules, thus with potential to reconstruct paleobiology in a broader temporal scale than more labile molecules50. Our results reinforce the notion that lipidic biomarkers should be preferred targets in the search for extinct and/or extant life on Mars precisely because they are so recalcitrant. More
