Skip to main content
Download PDF

The symbiotic alga Trebouxia fuels a coherent soil ecosystem on the landscape scale in the Atacama Desert

Environmental Microbiome volume 19, Article number: 59 (2024) Cite this article

Abstract

Biocrusts represent associations of lichens, green algae, cyanobacteria, fungi and other microorganisms, colonizing soils in varying proportions of principally arid biomes. The so-called grit crust represents a recently discovered type of biocrust situated in the Coastal Range of the Atacama Desert (Chile) made of microorganisms growing on and in granitoid pebbles, resulting in a checkerboard pattern visible to the naked eye on the landscape scale. This specific microbiome fulfills a broad range of ecosystem services, all probably driven by fog and dew-induced photosynthetic activity of mainly micro-lichens. To understand its biodiversity and impact, we applied a polyphasic approach on the phototrophic microbiome of this biocrust, combining isolation and characterization of the lichen photobionts, multi-gene phylogeny of the photobionts and mycobionts based on a direct sequencing and microphotography approach, metabarcoding and determination of chlorophylla+b contents. Metabarcoding showed that yet undescribed lichens within the Caliciaceae dominated the biocrust together with Trebouxia as the most abundant eukaryote in all plots. Together with high mean chlorophylla+b contents exceeding 410 mg m−2, this distinguished the symbiotic algae Trebouxia as the main driver of the grit crust ecosystem. The trebouxioid photobionts could be assigned to the I (T. impressa/gelatinosa) and A (T. arboricola) clades and represented several lineages containing five potential species candidates, which were identified based on the unique phylogenetic position, morphological features, and developmental cycles of the corresponding isolates. These results designate the grit crust as the only known coherent soil layer with significant landscape covering impact of at least 440 km2, predominantly ruled by a single symbiotic algal genus.

Introduction

One of the most abundant complex and natural microbial systems are biological soil crusts (biocrusts), which occur globally on 12% of the terrestrial surface, especially in water-limited environments where vascular plants are restricted [1]. They are considered pioneer communities in harsh environments such as polar- and non-polar deserts because they are composed of a predominantly photoautotrophic, poikilohydric group of organisms called cryptogams (meaning ‘hidden reproduction’): green algae, cyanobacteria, bryophytes as well as lichens together with heterotrophic bacteria and fungi [2]. All together, these organisms can colonize the first millimeters of soil and concatenate particles with their extracellular polymeric substances, filaments, and root-like structures, resulting in a physical, biotic crust [3]. Depending on soil characteristics, biome and also microbial community composition biocrusts can be divided into several types according to functional groups and/or the dominant organism [4,5,6]. Biocrusts from polar regions, for example, are often dominated by bryophytes [7,8,9], and those of the Negev and Sahara Desert by cyanobacteria [10, 11] while lichens occurring in small patches are the predominant biocrust organisms of the Colorado Plateau [12], or in the Namib Desert [13]. All of these biocrust types have significant and multifactorial influences on their environment by fulfilling several ecosystem services such as erosion prevention [14], but also modulating nitrogen-, phosphorus- and carbon cycles [15]. In addition, biocrusts exert strong influences on hydrological processes in drylands by modifying numerous soil properties that affect water retention and movement in soils such as infiltration, runoff, soil moisture and evaporation [16]. As ecosystem engineers, they can be pioneers in the colonization of inert substrates, where they play a crucial role in ecological succession by promoting initial soil formation or priming [17, 18], often resulting in the subsequent establishment of vascular plants [19, 20]. All of these important ecosystem services are the result of the extremotolerant character of cryptogams which allows them to withstand extreme temperature amplitudes and low precipitation through desiccation tolerance [21,22,23], but also their complex organismic [24, 25] and symbiotic nature spanning domains, kingdoms, and phyla [26, 27].

Among such symbiotic associations of biocrusts are the lichens. In general, they represent an iconic example of symbiosis where green algae and/or cyanobacteria (photobionts) team up with a fungus (mycobiont), and include other symbiotic microorganisms, which significantly contribute to the symbiosis (e.g., [28,29,30]). The distribution and habitat preferences of lichens are shaped by the requirements of all symbionts forming the holobiont [31], and according to this the symbionts allow their hosts to live in habitats from which they would be otherwise excluded. Their complexity has recently been acknowledged by redefining lichens as self-sustaining complex micro-ecosystems [32], albeit the role of other microorganisms—apart from the mycobiont and photobiont—are not demonstrably known [33, 34]. The debate about this terminology exemplifies how crucial the understanding of such mycobiont-photobiont relationships are, adding another diversity level to biocrusts. However, lichens and especially their photobionts associated with biocrusts have rarely been investigated [8, 35], which might be because such lichens often occur in isolated patches due to their mainly squamulose growth forms. This isolated appearance is the case for Psora decipiens from the Colorado Plateau [12] or Teloschistes capensis from the Namib Desert [13], while, for example, cyanobacteria fill the lichen-free spaces between them resulting in a much higher cohesive coverage than the single lichens. Consequently, a patchy presence of lichens within biocrusts also limits the appearance as well as the functional role of the lichen photobionts. Such lichen photobionts are for example members of the green algal genus Trebouxia, one of the most abundant lichen photobionts [36,37,38], which has rarely been found free living [39,40,41,42]. The taxonomic system for the diverse photobiont genus Trebouxia was recently updated, dividing the genus into four main lineages: A- (T. arboricola), C- (T. corticola), I- (T. impressa / T. gelatinosa) and S-clade (T. simplex / T. jamesii). This new framework is based on historical taxonomic work on isolated Trebouxia-species [43,44,45,46,47] in combination with modern multilocus phylogenies and morpho-anatomical features such as chloroplast architecture [48,49,50]. Following this approach, species-level lineages within each of the four major clades can also be designated, including potential species candidates, based on the phylogenetic position and morphology of representative isolates [51, 52].

So far, such a comprehensive approach has not been carried out for any biocrust environment and provides an interesting opportunity to evaluate the role of photobionts in a novel and unique lichen-dominated grit crust environment from the Atacama Desert [53]. This microbiome has been described as a transitional stage between a saxicolous community (growing on rocks) and a biocrust because the microorganisms are colonizing granitoid pebbles called grits of about 6 mm in diameter and not a fine substrate equivalent to soil [3, 53]. The microorganisms creating the grit crust reach comparably high densities, resulting in a blackish coloration of the grits, so that the biocrust’s extension appears as a checkerboard pattern on the landscape scale covering about 80% of the 440 km2 encompassing National Park Pan de Azucar [53]. There, abiotic conditions are not as harsh as in hyper-arid parts of the Atacama Desert, but radiation is still extreme, and water availability is restricted to fog and dew [54]. The biodiversity of the grit crust has not yet been fully understood, but important roles in CO2 fixation [53], bioweathering, initial pedogenesis [55], implications for astrobiology [56], and nature conservation [57], have been attributed mainly to lichens. Those lichens can best be described as crustose- to squamulose micro-lichens forming thalli not bigger than a few millimeters where often one specimen colonizes only a single grit, and they have a reduced set of morphological features (e.g., rare formation of apothecia), all of which exacerbates detailed studies. This is in contrast to any other known biocrust type or environment where lichens mostly form macroscopic, fruticose and/or large squamulose thalli in patches [58,59,60], while the micro-lichens of the grit crust form a coherent crust on the landscape scale [53]. In addition, insights into the ecology and phylogeny of lichen mycobionts and their photobionts from South America are rare [61, 62], but this holds particularly true for the Atacama Desert, which is considered a lichen hotspot, not only due to its fog oases [63,64,65] but also to its striking aridity gradient that correlates with latitude [66] and elevation [67]. Up to now, detailed insights into the microbial diversity of the grit crust are missing including information on symbiotic interactions such as between lichen mycobionts and their algal photobionts. It also remains an open question whether a single photobiont population or a diverse set of different photobionts or -lineages have been successful under the unique climate regime of the location leading to the broad, cohesive extent of the grit crust.

For this reason, we provide a polyphasic approach on grit crust lichens in order to answer questions about the organismic composition of the biocenosis, the lichen’s host-photobiont specificity, and the ecological role of their photobionts in the environment. We combine a Sanger-sequencing approach for taxonomic gene regions of myco- and photobionts based on low biomass picked directly from single lichens growing on single grits, which we micro-photographed, along with isolation of the photobionts. Those small-scale analyses were underpinned by metabarcoding and chlorophylla+b content determination alongside eleven locations of the grit crust community distributed over the 440 km2 spanning National Park Pan de Azucar. Due to the often described outstanding character of the grit crust and the great dominance of green algal lichens we combined classic biocrust research methods (metabarcoding, chlorophyll analysis) with analyses designed specifically for lichens and their trebouxioid photobionts (isolation, morphological studies of life cycles and chloroplasts, multigene-phylogenies). This unusual methodological combination helped us frame the grit crust's microbiome, to uncover the identity of the major lichen groups and their photobiont lineages, including the designation of unique potential photobiont species candidates. On the other hand, we intend to uncover the role of the grit crust’s lichen photobionts as the prime organism on the landscape scale amongst other biocrust-associated primary producers, which has so far never been exposed.

Material and methods

Sampling location

The National Park Pan de Azúcar, situated within the arid expanse of the Atacama Desert in northern Chile, spans approximately 440 km2 and encompasses a diverse range of ecological landscapes. The National Park Pan de Azúcar is a crucial site for scientific inquiry and preservation efforts in the context of arid and coastal ecosystems (Fig. 1A), which was highlighted during the research project EarthShape. As a result, geological and climatological information on the National Park can now be accessed [15, 25, 68].

Fig. 1
figure 1

Grit crust landscape, sSampling locations and Chlorophylla+b m−2 content. A close-up photography of the grit crust in the foreground and the patchy distribution in the landscape in the background. B map of the National Park Pan de Azucar including the 11 sampling locations where plots were established for metabarcoding and chlorophyll determination. C Chlorophylla+b m−2 content of each plot

Full size image

Sampling

In total, 11 homogeneously colonized plots with 1 m2 area each were established in the National Park Pan de Azucar (Fig. 1A, B) within comparable flat areas covered with densely colonized grit stones which appeared black [53, 55]. The corners of the plots were marked with long metal screws and the sides of the plots were arranged according to the cardinal points.

Within each of these plots the top centimeter of a 10 cm2 area was removed with a sterile brush and collected in plastic tubes for chlorophyll determination.

A second representative 10 cm2 area within the same plot was sampled as described above for a metabarcoding approach. Those tubes were filled with LifeGuard Soil Preservation Solution (Qiagen, Hilden, Germany) directly after sampling in order to preserve them.

Additionally, three samples of 50 g top-soil layer, each, were collected with a spatula directly beneath the 1 m2 plots and stored in plastic containers for the subsequent direct analyses of the lichens and their photobionts and the isolation of the latter.

The samples for chlorophyll determination and the taxonomic investigations of the lichens were stored air-dried at room temperature in the dark, while the preserved samples were stored at − 20 °C in the freezer until the beginning of the study.

Chlorophyll extraction and quantification

In order to estimate the chlorophylla+b content of the grit crust, the top centimeter of a 10 cm2 area within the plots was removed with a brush and collected. For the chlorophyll extraction the dimethylsulfoxide (DMSO) method after Ronen and Galun [69] was applied which has been shown to be most suitable for biocrust samples [70].

DNA extraction, amplification and sequencing of metabarcoding samples

Prior to DNA extraction, the LifeGuard solution was removed by centrifuging the samples at 2500 g for 5 min and discarding the supernatant. Due to colonization also of inner structures of the grit stones the samples were homogenized by beat-beating with steel beads (1 mm and 2.3 mm diameter) and a powder mill (Retsch, Haan, Germany) according to the DNeasy PowerSoil kit (Qiagen, Hilden, Germany), which was used for DNA extraction. The V4 region of the 18S rRNA gene (~ 390–430 bp) was amplified for eukaryotes using the universal eukaryotic V4 primers TAReuk454FWD1 and TAReukRev3 [71] and the ITS2 (~ 300–400 bp) region for fungi using the primer pair ITS3_KYO2/ITS4 [72]. The PCRs according to the conditions for each gene region, library preparation following the Tagsteady protocol [73] and sequencing was carried out by Fasteris (Plan-les-Ouates, Switzerland) using an Illumina MiSeq platform (2 × 300 bp). The full bioinformatic analyses were automatically processed using DeltaMP v0.6 [74] by high-performance computing (“Centre de Calcul de la Faculté des Sciences”, University of Neuchâtel, Switzerland). Raw Illumina reads were demultiplexed using cutadapt v4.1 [75] allowing up to 2 mismatches on the barcode sequences (8-nt long plus 2 to 4-nt long heterogeneity spacers at 5’-end). After stripping the primer sequences at the 5′-end with cutadapt (up to 3 mismatches allowed for 18S rRNA datasets and 5 for ITS dataset), only the amplified biological sequences were further analyzed. These first step produced four sets of reads, that is paired libraries (R1 and R2) for each orientation (i.e. R1 library starting either with the forward or the reverse primer). End of reads were stripped due to low quality. An automatic algorithm optimized this process in order to keep 75% of raw demultiplexed reads while maintaining an expected error on the remaining stripped reads below two and enough length to allow for pair end assembly. In the 18S rRNA dataset, forward orientation libraries were trimmed at 265 (R1) and 190 nt (R2) and reverse orientation libraries were trimmed at 270 (R1) and 225 nt (R2). In the ITS dataset, forward orientation libraries were trimmed at 260 (R1) and 170 nt (R2) and reverse orientation libraries were trimmed at 235 (R1) and 180 nt (R2). Amplicon sequence variants (ASV) were called using the DADA2 R package v1.26 [76]. Error models were built separately for each library and orientation, with turning the pool option on. Paired reads in R1 and R2 libraries were assembled with the dada2 mergePairs function, with a minimum overlap length of 10 nucleotides and a maximum of 1 mismatch in the overlap region. ASV originally with the reverse primer in the R1 library were reverse complemented using the seqinr package v4.2-16 [77]. The presence of chimera was double-checked first with the dada2 function removeBimeraDenovo and with the UCHIME de-novo algorithm as implemented in VSEARCH v2.22.1 [78, 79]. The databases used to identify ASVs were PR2 v5.0.0 [80] for the 18S dataset and UNITE v9.0 (UNITE [81, 82]) for the ITS dataset. The UNITE database was augmented with sequences of isolates produced in this study (53 Lecanoromycetes and 49 Chlorophyta, see below). The taxonomy was assigned to each ASV using the VSEARCH option–usearch_global. When multiple best matches were found in the database, a consensus taxonomy was resolved at a threshold of 60%. The final ASV matrices contain 2.308.569 18S reads and 2.493.453 ITS reads distributed over 4.401 18S ASVs and 18.421 ITS ASVs. Statistical analyses were performed using R (v4.2.3) [83] and the package vegan (v2.5-7) [84], if not specified otherwise. For the 18S dataset ASVs not belonging to protist (Metazoa, Fungi and Embryophyta) were removed, as well as the ones having a percentage of identity below 70% with the database reference sequences and those that have less than 3 reads and that are present in only in one sample. For the ITS dataset, ASVs with a percentage of identity below 60% and that have less than 3 reads and that are present in only one sample were removed in addition to ITS sequences assigned to Chlorophyta in order to focus on fungi. General ecological functions of fungi were assigned based on the most updated literature for the corresponding ITS-derived ASVs consolidating the FungalTraits database [85]. Metabarcoding data were submitted to the European Nucleotide Archive (ENA) under the project code PRJEB72845 https://www.ebi.ac.uk/ena/browser/view/PRJEB72845.

Isolation and cultivation of Lichen photobionts

Single grit stones with lichens growing attached to them were inspected under a binocular stereoscope (Stemi 508, Carl Zeiss Microscopy GmbH, Jena, Germany) under sterile conditions while small parts of the lichens were taken with a sterile metal needle and transferred into sterile PCR tubes filled with 200 µl Bold’s Basal Medium (BBM, [86]). The inoculated PCR tubes were subsequently streaked out on 12-well plates filled with solidified BBM (0.9%), sealed with parafilm and incubated in a culture room for four weeks at 17 °C, 30 µmol photons m−2 s−1 and a light:dark cycle of 18:6 h. Afterwards the growth of the colonies was inspected under the binocular stereoscope and transferred to fresh BBM plates until unialgal cultures were obtained. As a final step the 23 unialgal cultures were transferred to culture flasks filled with liquid BBM.

Photo documentation and sequencing from samples and isolates

A direct PCR approach was applied to the micro-lichens growing on the grit stones forming the grit crust of the Atacama Desert. Single stones with lichens were selected under a binocular stereoscope, hydrated with sterile ddH2O and photographed using a 3D 4 K stereo microscope (VHX-7000, Keyence Deutschland GmbH, Neu-Isenburg, Germany). Subsequently a sterile needle was taken to pick small proportions of the lichens which were transferred to PCR tubes filled with 20 µl lysis buffer of the Platinum Direct PCR Universal Master Mix-Kit (Thermo Fisher Scientific Inc). The combination of photographing and extraction of biomass for molecular work allowed us to correctly assign DNA sequences of mycobionts and their photobionts to the morphology of the lichens. The DNA extraction of the biomass was carried out as described by the manufacturer including the setup of the mastermix and slightly modified [87]. For amplification of the ITS1 gene of the mycobiont the primers ITS1f [88] and LR3 [89] were used, for the green algal photobionts the nuclear internal transcribed spacer of the ribosomal DNA (ITS, ITS1, 5.8S, ITS2 was amplified using the primers ITS1T and ITS4T [90] with the identical lysate as DNA template.

In order to increase phylogenetic resolution for the isolated Trebouxia strains, rbcL was amplified using the primers rbcL1T and rbcL2T [91].

Success of amplification was checked by gel electrophoresis. PCR products that contained sufficient DNA concentration were cleaned using the NucleoSpin® Gel and PCR Clean-up Kit (Marchery Nagel, New England, Canada) and sent to Genewiz (Göttingen, Germany) for Sanger sequencing.

Phylogeny of mycobionts and photobionts

The generated sequences were merged and the primers were trimmed using the software Geneious Prime (2022.01.1). Sequences for individual loci were assigned based on the NCBI GenBank dataset using the Basic Local Alignment Search Tool (BLAST) of Mega 11 [92]. Then, multiple sequence alignments were performed using the multiple sequence alignment method (Muscle) algorithm (Edgar 2004) implemented in Mega 11. The final datasets of the alignments consisted of 143 sequences of 983 bp comprising phylogeny of the ITS1 gene region of the mycobionts and a 2241 bp sequences comprising the concatenated ITS1, 5.8S, ITS2 and rbcL phylogeny of the photobionts derived from direct PCR on lichen material and the 23 isolates. For this dataset sequences from the study of Kosecka et al. [61] were also incorporated (Suppl. Tab1), since here sequences covering the same gene regions of photobionts of lichens from South America (Bolivia) were generated. We also incorporated sequences generated by Muggia et al. [50]. For each dataset the evolutionary model that was best suited to the database used was selected on the basis of the lowest Akaike Information Criterion (AIC) value and calculated in Mega 11 which was GTR + G + I for all. The phylogenetic reconstructions were calculated using NGPhylogeny.fr [93] with 500 bootstrap replications.

Besides Maximum Likelihood analyses (ML) also Bayesian Inference (BI) was calculated for each of these trees using Mr. Bayes 3.2.1 for the later analysis [94] as described in Kosecka et al. [61]. The tree topologies obtained by the BI method did not contradict the ML trees, thus only the ML trees are shown. Branches with bootstrap support and posterior probabilities > 80% were considered to be strongly supported and are marked by an asterisk.

Whenever possible, we used the numbering system for the designated photobiont clades according to Muggia et al. [50] and Kosecka et al. [61], but for clades which were newly designated during this study we used our own system. Phylogenetic trees were finally depicted in iTOL [95] and the resulting figures were edited in BioRender.

We also chose to prepare a phylogenetic reconstruction of the metabarcoding-generated ITS sequences on the ASV level which were assigned to Trebouxia instead of using the 18S ASV sequences because microalgae can show ecophysiological differences even when they are phylogenetically close, which is best displayed in ITS sequences for Trebouxia [50, 96]. Those sequences were combined with the full ITS-covering dataset from Muggia et al. [50] and the sequences generated during this study. The alignments were calculated using Mafft [93], manually trimmed in Mega 11 leading to a final dataset of 7.470 sequences of 561 bp length covering the ITS2 gene region. The phylogenetic tree was constructed using the GTR + G substitution model set to 500 replications at NGPhylogeny.fr [93].

Microscopy

The 3D 4 K stereo microscope (VHX-7000, Keyencnce, Neu-Isenburg, Germany) was used as described above to capture in situ images during the growth of the isolated Trebouxia photobionts after 6 weeks on BBM agar plates.

In addition, light- and DIC microscopy were applied to take detailed images of the isolated Trebouxia photobionts using a BX51 (Olympus, Tokyo) coupled with a camera (MicroLive, Bremen, Germany) and MicroLive 5 software (MicroLive, Bremen, Germany).

Documentation of the developmental stages of the photobionts within the lichens were also made. This was carried out by removing a hydrated lichen thallus from the grit stones with tweezers under a binocular stereoscope and transferring it to a drop of water on a glass object slide. The thallus was afterwards slightly squeezed by pressing a cover slide on the object which was then investigated under the microscope.

To depict the chloroplast structure of the isolated Trebouxia confocal laser scanning microscopy (CLSM) was carried out using a Leica SP8 with an HC PL APO × 63/1.40 oil CS2 objective. Cells were immobilized on 1% (w/v) agar pads and excited at 458 nm. Fluorescence emission was recorded from 650 to 750 nm. Images were analyzed with Leica Application Suite X (version 3.5.2) and FIJI [97].

Results

Chlorophylla+b content as a proxy for photosynthetic biomass

On average 420 mg Chla+b per m2 were detected with a variability of 20%. In general, plots closer to the coast had higher Chla+b contents (e.g.: 1 = 672 mg m−2, 4 = 879 mg m−2) than others in distance to the coast (e.g.: 2 = 276 mg m−2, 3 = 261 mg m−2, 7 = 21 mg m−2) (Fig. 1C).

Metabarcoding of eukaryotes and fungi associated with the grit crust

The eukaryotic community of the grit crust microbiome was dominated by Chlorophyta (18S, 50%) and Fungi (ITS, 20–60%) (Fig. 2), whereby the latter contained mainly lichenized fungi referring to lichens as the dominant organismic group. Among Chlorophyta, more than 80% of all 18S-derived ASVs could be assigned solely to the genus Trebouxia but other genera such as Stichococcus, Myrmecia, Elliptochloris and Diplosphaera were detected (Fig. 2).

Fig. 2
figure 2

Metabarcoding. A Relative abundance of fungi-assigned ITS covering ASVs. B Relative abundance of Lecanoromycetes-assigned ITS covering ASVs. C Relative abundance of eukaryote-assigned 18S covering ASVs. D Relative abundance of Chlorophyta-assigned 18S covering ASVs

Full size image

More than 50% of all ITS-derived ASVs of the fungal microbiome were assigned to Lecanoromycetes (lichenized fungi), followed by less than 5% of each, Tremellomycetes (lichenicolous, parasitic fungi), Dothideomycetes (lichenized endophytes) and Eurotiomycetes (lichenized saprobiontic) (Fig. 2). Among the Lecanoromycetes the predominant order was Caliciales (60%) but also members of other lichens were detected such as Acarosporales, Pertusariales, Lecanorales and Teloschistales (Fig. 2).

Both targeted gene regions, ITS and 18S rRNA, covered algae and fungi but it turned out that, for example, Lecanoromycetes were underrepresented in the 18S rRNA dataset (61 ASVs) while the diversity of eukaryotic algae in general was underrepresented in the ITS dataset. For these reasons characterization of eukaryotic algae of the grit crust was best represented by the 18S dataset while fungal communities show a more realistic representation based on the ITS dataset, which has been commonly accepted [98].

Phylogeny of ASVs assigned to Trebouxia

The metabarcoding of the ITS2 gene region resulted in 6.406 ASVs assigned to Trebouxia. 96.08% of all grit crust Trebouxia ASVs could be assigned to the A clade, while 3.92% fell into the I clade. No ASVs of the C and S clade were detected. Within the I clade, the majority of ASV derived from the grit crust formed a separate cluster (Fig. 3, box and arrow) in distance to other sequences from epiphytic lichens (Fig. 3, arrow) sharing the same ecosystem. The grit crust derived ASVs of the A clade were distributed across the whole clade.

Fig. 3
figure 3

Phylogeny of the 6.406 Trebouxia-assigned ITS metabarcoding ASVs. The phylogenetic reconstruction represents the ITS2 covering sequences derived from metabarcoding of the grit crust and also includes the dataset of Muggia et al. [50], as well as the sequences created during this study based on the direct approach

Full size image

Phylogenies of mycobionts, their uncultured photobionts and photobiont isolates

The combination of microphotography of lichens on the grit stones and direct sequencing of the corresponding myco- and photobionts allowed for the determination of the phylogenetic position of 41 uncultured Trebouxia-photobionts and additional 23 isolates (Fig. 4) as well as 69 mycobiont sequences (Figs. 5, 6).

Fig. 4
figure 4

Concatenated Trebouxia phylogeny covering the ITS1, 5.8S, ITS2 and rbcL gene regions of sequences derived from direct sequencing of lichen material and isolates. A four major Trebouxia clades A (green), I (blue), C (orange) and S (red) after [50] and their assigned predominant pyrenoid and chloroplast structure reproduced after Bordenave et al. [48]. B, C phylogenetic reconstruction of the I clade (blue) and the A clade (green), respectively. Novel lineages are highlighted in color while already detected lineages are marked in gray. Stars mark potential species candidates for each novel lineage based on its phylogenetic positioning and morphological traits. Bold sequences without color boxes are derived from Trebouxia photobionts of epiphytic lichens sharing the same habitat but are not part of the grit crust. Orange circles indicate the number of the genus the corresponding mycobiont can be assigned to referring to Figs. 5, 6. Branches supported > 80% by Bayesian Inference (BI) and Maximum Likelihood bootstrap (ML) are labeled by an asterisk

Full size image
Fig. 5
figure 5

Grit Crust mycobiont phylogeny covering the ITS1 gene region derived from direct sequencing of lichen material. Orange circles number the mycobiont genera while green and blue circles indicate the lineage of the corresponding photobiont. Microphotographs show the exact origin of the biomass from which the mycobiont sequences and the corresponding photobiont sequences were generated. Scale bars are 500 µm. Branches supported > 80% by Bayesian Inference (BI) and Maximum Likelihood bootstrap (ML) are labeled by an asterisk

Full size image
Fig. 6
figure 6

Grit Crust mycobiont phylogeny of Caliciaceae covering the ITS1 gene region derived from direct sequencing of lichen material. Orange circles number the mycobiont genera while green and blue circles indicate the lineage of the corresponding photobiont. Microphotographs show the exact origin of the biomass from which the mycobiont sequences and the corresponding photobiont sequences were generated. Scale bars are 1000 µm. Branches supported > 80% by Bayesian Inference (BI) and Maximum Likelihood bootstrap (ML) are labeled by an asterisk

Full size image

In detail, five uncultured Trebouxia photobionts and one isolate were assigned to the I clade (T. impressa/gelatinosa) based on the latest phylogenetic reconstruction of this genus [50] including the gene regions ITS1, 5.8S, ITS2 and rbcL (Fig. 4B). Together, those six sequences formed their own clade (I42) and other well-defined clades of photobiont sequences from the epiphytic lichens Everniopsis trulla (I40; [64]), Cenocosia cava, C. excorticola, Heterodermia adunca and H. follmannii (I02; [57, 65, 99]) sharing the same habitat were also detected within the I clade (Fig. 4B, bold).

Additionally, 36 uncultured Trebouxia photobionts and 22 isolates were assigned to the A clade (T. arboricola) (Fig. 4C). Within this A clade ten new clades were identified and a few sequences fell within three clades previously detected based on Trebouxia-sequences from Bolivia by Kosecka et al. [61] (A07, A53, A18). The ten new clades also comprised the photobionts isolated from Acarospora conafii (A22; [64]) but were not related to photobionts isolated from the terricolous lichen Placidium sp. (A04), the epiphytic Ramalina thrausta (A11) which both share the same habitat [64] or sequences generated from Caloplaca lichens from other areas of the Atacama Desert (A19; [63]).

Based on the corresponding lichen material used to generate the photobiont phylogenies 69 ITS1 sequences were generated (Figs. 5, 6), which were assigned to eight different clusters within the Caliciaceae representing potential new genera (Fig. 5) and to the genera Caloplaca, Wetmoreana, Candelariella, Acarospora conafii, Aspicilia, and one undetermined, potential new genus within the Pertusariales (Fig. 5). In some cases biomass and therewith DNA was extracted from black lichen prothallus, apothecia and lichen thalli of the same lichen in order to test if these structures belonged to the same species, which was always the case (e.g.: Fig. 5, FPBH2, FPB1a, FPB1b).

Whenever the sequencing of the myco- and the photobiont of the same lichen was successful, their affiliation was indicated in the corresponding phylogenetic tree in order to get preliminary insights into the photobiont-host specificity (Fig. 4B, C orange circles; Figs. 5, 6 green/blue circles). Based on these results a rather low photobiont-host specificity was detected, where lichens of several genera (Acarospora conafii, Aspicilia, Caliciaceae cluster 7 and 8) shared Trebouxia photobionts of a single clade (A22), but also that Trebouxia-photobionts of different clades (I32, A23) teamed up with lichen mycobionts of a single clade (Caliciaceae cluster 6).

Trebouxia species candidates and developmental stages

Potential species candidates within the designated lineages were evaluated based on morphological criteria of the isolates such as cell morphology, developmental stages in the lichen and in culture, chloroplast structure and their phylogenetic position deduced from the multigene phylogeny. Following this, five isolated species candidates could be identified: 19C (I clade), 1F, 1K, 17G and 17D2 (A clade) (Fig. 7).

Fig. 7
figure 7

Potential Trebouxia-species candidates and their developmental cycles according to own observations after Friedl [44]. Top left shows the grit crust stones with the exact origin (encircled) of each isolate with the number of the isolate in the right corner as well as the Trebouxia cluster in the colored circle. Mid row shows microscopic images of each isolate. Right row shows 4 k microscopic images of the cultures grown on agar. Bottom: Broken lines: reproduction of Trebouxia within lichen thalli as interpreted from developmental stages seen in squashed preparations of algal layers of lichens. Cell cycle A: (1) zoospore; (2) juvenile vegetative cell (autospore); (3) adult vegetative cell; (4) protoplast division, chloroplast already divided during autosporangium (type 4) formation; (5) autosporangium (type 4) containing few autospores compressed within the sporangial wall; (6) autosporangium with autospores that begin to dissociate and where new protoplast divisions occur while the autospores that begin to dissociate and where new protoplast divisions occur while the autospores still adhere within the sporangial wall; (7) some autospores of stage 6 are developed into sporangia with numerous daughter cells which either escape as nonmotile autospores (a, autosporangium) or as free-swimming zoospores (z, zoosporangium) while other autospores are developed into an autosporangia (a) with few cells where new protoplast divisions occur; (8) autosporangium of 6 dissociates into single cells. Cell cycle B: (1) zoospore; (2) juvenile vegetative cell; (3) adult vegetative cell; (4) protoplast division, chloroplast already divided into four during autosporangium (type 4.1) formation; (5) cell containing numerous protoplasts (type 5.1 autosporangium); (6) development from stage 5.1 either into a sporangium with numerous nonmotile autospores (a) or into a zoosporangium (z). Microscopic images show all cell types / developmental stages

Full size image

Isolate 19C shared a unique cluster (here named I32) with other sequences of uncultured Trebouxia-photobionts from the grit crust (Fig. 4B, C). The developmental stages of 19C correspond to cycle B according to developmental life cycles after Friedl [44] (Fig. 7), because irregular, spherical autosporangia only of the type 5.1 were formed in culture and in the lichen thallus. The cells were equally rounded, up to 23 µm in diameter, and showed a structured cell surface and a weakly lobed chloroplast filling up 2/3 of the cell, with distinct nucleus and often several, up to 10 µm large vacuoles (Fig. 7, top row). Zoosporangia reached up to 30 µm and were frequently detected. On agar plates a rather flat biofilm was produced with dispersed cells. The chloroplast was divided into almost regularly formed, rounded lobes (Fig. 8).

Fig. 8
figure 8

Confocal laser scanning microscopy (CLSM) showing the chloroplast structure of the Trebouxia sp. isolates. Gray images show the Trebouxia isolate and the assigned Trebouxia cluster in the colored circle. Arrows mark the cell which was visualized in detail

Full size image

Isolate 1F had an isolated position (here named A15) within the multigene phylogeny (Fig. 4) and was further distinguished by highly aggregating cells of irregular shape, comparable to short pseudo-filaments resulting from incomplete deliberation of autosporangia of type 5 (Fig. 7), which was also observed in the lichen thallus. As a result 1F was assigned to the A life cycle (Fig. 7), but did not show zoosporangia and consequently also no zoospores in culture. The irregular, rounded cells had a diameter of 12 to 20 µm with a chloroplast which was fringed at the surface and filled the full cell (Fig. 7, second row). The isolate formed irregular colonies made of aggregated cells on agar which resulted in several layers and a conical structure. Confocal laser scanning microscopy revealed a lobed chloroplast with an irregularly formed fringed partitioning, sometimes forming small conical protrusions (Fig. 8).

Together with an uncultured Trebouxia from the grit crust the sequence of isolate 1K joined one other Trebouxia sequence (AV077; cluster A07) defined by Voytsekhovich and Beck [47] from Lecanora albescens from Ukraine. Cells of the isolate were round, up to 32 µm in diameter, with a broadly lobed, centered chloroplast filling up to 2/3 of the cell. An irregular structure on the cell surface was always visible. Developmental stages could be assigned to the B cycle with round autosporangia exclusively of the 5.1 type (Fig. 7) and up to 40 µm in diameter, which contained hundreds of round autospores (Fig. 7, third row). Zoosporangia and zoospores were not detected. Regular, conical to rounded colonies were formed on agar with loosely attached cells. The chloroplast was divided into irregular, thin lobes which extended into fringes near the cell periphery (Fig. 8).

The isolate 17D2 joined the cluster A24 which was only formed by Trebouxia-sequences generated from lichens of the grit crust (Fig. 4). The morphology of 17D2 was identical to all remaining isolates of uniquely defined clusters within the A lineage (A14, A1, A22, A23, A25, A28), thus this specific isolate was randomly chosen as a potential species candidate reflecting the morphology but not the genetic position of the named clusters/isolates. Cells of these isolates were rounded, only slightly aggregated, up to 22 µm in diameter with a centered chloroplast which filled most of the cell (Fig. 7, fourth row). Developmental stages followed cycle A (Fig. 7) with the formation of round, regular autosporangia of 25 µm diameter containing densely pressed autospores of angular shape. Even in young autosporangia the autospores mimicked the morphology of adult cells, including the shape of the chloroplast. Zoosporangia and biflagellate zoospores without stigma were rarely detected in culture. During growth on agar plates those isolates developed irregular colonies where the different developmental stages were separated, reflecting the age of the cells (Fig. 7, fourth row). The chloroplast was sometimes shifted and shallowly lobed. The chloroplast’s outer surface was made of small, regular protrusions resembling the ‘typical’ T. arboricola chloroplast morphology (Fig. 8).

Isolate 17G clustered together with a single sequence (A53) from a trebouxioid photobiont of Myrolecis spp. from Bolivia (Fig. 4) generated by Medeiros et al. [62]. The cells of the isolate grew densely aggregated made of cell packages. The cells themselves were irregularly ovoid to spherical with a diameter of 12 µm (Fig. 7, bottom row). Chloroplasts fully filled the cells and had a coarsely structured surface. The isolate followed life cycle A in cultivation and in the lichen thallus with autosporangia of type 5 (Fig. 7). The autosporangia were irregularly formed and comprised only a few cells with autospores that mostly stayed attached to the other autospores forming small cell packages. This was also reflected in the growth on agar plates where single colonies were formed which accumulated, resulting in round, dense structures reaching several mm. Zoosporangia and zoospores were not observed. Confocal laser scanning microscopy showed deep incisions of the chloroplast leading to irregular, thick lobes (Fig. 8).

Discussion

Our holistic approach, applied to an unexplored area, enabled the characterization of biodiversity through different morphological and sequencing methods which, in combination, unveiled the hidden and impressive biodiversity of the photobionts as the integral part of the grit crust. This is in contrast to the majority of related studies on biocrusts, as they focus most often on the isolation or analysis of single phyla, thereby overlooking such vivid biodiversity. We demonstrated that an approach which spanned from microscopic observations on single micro-lichens to data generated from samples distributed over the 440 km2 wide National Park Pan de Azucar designated the symbiotic algal genus Trebouxia as the main primary producer of the grit crust, thriving under the unique conditions of the coastal Atacama Desert.

The grit crust lichen community: saxicolous or biocrust associated?

Based on previous and ongoing work, it has been suspected that free living cyanobacteria and green algae only minorly contribute to the grit crust, while lichens dominate both in terms of biomass and function [53, 55]. It is not unusual for lichens to be a dominant component of biocrusts, some authors for example expect between 500 to 1000 species to be involved [100], albeit their biodiversity is usually low and restricted to a few species [101, 102]. Other studies using amplicon sequencing techniques on biocrusts only focused on certain phyla [7, 11, 103], did not detect interactions between Trebouxiophyceae and fungi [104], or only defined single lichen species as contributing locally [24]. For the grit crust this is a different situation compared to all other biocrust types known worldwide. This is due to the diverse and dominant group of lichens assigned to the family Caliciaceae with at least 8 genetically different clades as part of the community, complemented by genera such as Candelariella, Acarospora, Wetmoreana, Caloplaca, Aspicilia and Pertusaria (Fig. 5). This can be explained by the fact that the grits, being small rocks, favor a putative saxicolous lichen community over a soil crust-associated community. A detailed observation of these lichens shows that this is an oversimplification as all Acarosporales and a few Teloschistales ASVs can be assigned to Acarospora conafii and an undescribed species related to Wetmoreana genus. Both are strictly terricolous lichens from the same location growing only between the grits in the fine substrate and never on them [64, 99]. In addition, the core group of lichens contributing to the grit crust are lichens of the family Caliciaceae, which usually grow on diverse substrates such as bark, dead wood or rocks [105]. The rather graduated ecology of the non-epiphytic lichens from this habitat thus strengthens the definition of the grit crust as a transitional stage between a saxicolous community and a biological soil crust as previously outlined [3, 53]. Some members of the Caliciaceae, particularly those of the type genus, Calicium, are characterized by the presence of a mazaedium. This is an accumulation of loose, maturing spores covering the surface of the fruiting body, which is usually black, and either sits atop a long thin stalk, or rests on (or is sometimes immersed within) the surface of the lichen substrate. Consequently this results in a passive spore dispersal, a rare feature amongst Ascomycota [106]. However, many other Caliciaceae species generate spores in an apothecium, which typically resembles a flattened black disc (discoid, sessile and lecidiene ascomata) which is the case for all Caliciaceae-lichens of the grit crust (Fig. 4). In addition, most Caliciaceae-lichens of the grit crust have a well-developed prothallus made of a venous network of melanized, black hyphae stretching out several millimeters across the grit stones (Fig. 4; see also [53, 56]). In combination, the unique phylogenetic position of the Caliciaceae grit crust lichens based on the ITS gene region, the distinction of single, well supported clusters, the mutual absence of a mazaedium, their prothallus network and their extremotolerant ecology, will certainly justify the description of new genera and species. Most Caliciaceae species have a cosmopolitan distribution, although they are especially predominant in temperate and tropical montane areas, where they tend to avoid competition with other lichens and often grow in microhabitats uninhabited by other lichens [107]. This could explain the dominance of a unique Caliciaceae lichen subset thriving in such a habitat, which has been stable with an on-going aridity at least since the mid Miocene, 16–12 million years ago [108, 109], or even longer [110], leading to diversification over time.

Alongside the development of such an isolated lichen community are probably adaptations to the abiotic conditions such as the extreme radiation and the restricted but cyclic water input via probably fog and dew, which can be disadvantageous to other biocrusts [111, 112]. Also within the National Park Pan de Azucar differences in abiotic conditions affecting the grit crust microbiome exist, albeit mainly concerning water availability such as increased amount of fog water deposition close to the coast areas (climate data given for location 4 in [64] compared to topologically more enclosed inland areas (climate data given for location 7 in [68]). Since especially water availability is well known to shape lichen (-photobiont) communities this factor might explain diversity patterns but this is part of ongoing research in the area addressing for example fog water collection [113, 114].

However, adaptations of the organisms can on the one hand be supported, for example, by the melanized prothallus network of the grit crust lichens, which has been reported as a protection mechanism against UV radiation in other fungi and lichens from comparable ecosystems [115, 116]. On the other hand, ecophysiological adaptations of the grit crust community mirroring the abiotic conditions of the ecosystem have already been shown, identifying a comparably narrow ecological niche for optimum net primary productivity above 10 °C during the day, with no decrease in net photosynthesis up to at least 20 °C, 0.25 mm water content and a comparably high light compensation point (LCP) of 460 µmol m−2 s−1 [53]. Especially the exceptionally low water content required for optimum photosynthesis compared to other biocrust lichens is remarkable, since, for example, 6 mm were reported for Acarospora schleicheri from the fog zones of the Namib Desert [117] as well as 0.5–1 mm for lichen-dominated biocrusts from the Sonoran Desert [118]. However, the ecophysiological performance of lichens are most likely controlled, at least to a certain degree, by the lichen photobiont providing the carbon source for the whole lichen holobiont [119], supporting the theory that isolated lichen communities do not exist without an equally unique photobiont population [52, 61]. The low amount of water required for optimum photosynthesis of the grit crust might also be the key factor determining the dominance of green algae. In general, green algae and chlorolichens require minimal amounts of water for the activation of their photosynthesis and can do so even based on high relative air humidity alone, while cyanobacteria and cyanolichens depend on liquid water and higher amounts of water compared to green algae [120].

Trebouxia: the dominant organism fueling the ecosystem of the grit crust

The lichen-forming algal genus Trebouxia is one of the best-studied symbiotic organisms associating with over 7.000 fungal species found in virtually all biomes linked to a diversification which took place during mid‐Paleozoic to mid‐Mesozoic [121, 122]. Getting closer to understanding this diversification process is a major challenge that was recently approached by defining four major Trebouxia-lineages A, C, I and S [48, 50]. According to this, certain lineages could be detected in a certain ecosystem, such as Icelandic cetrarioid lichens sharing photobionts of the S-clade [91], but a rather mixed subset was detected in most cases, especially when extreme ecosystems were considered [52, 61]. However, it has been proposed that the crown I-clade occupies warmer, as well as cooler and drier, habitats, i.e., Eastern Europe, while clade C occupies hot and humid climates in partially or exclusively forested habitats [123]. Moreover, the stem A-clade is dominant in only or partially forested habitats, but subsequently expanded in the Cenozoic Era to occupy regimes characterized by cooler and drier habitats—loosely coinciding with or after clade S invasion to similar climatic regimes [123]. Our study now sheds a new light on this concept since the Atacama Desert represents a climatically long-term stable, extreme, deforested, cool and dry environment from which hardly any lichen photobiont data exist. According to this, we found that the majority of grit crust associated lichens shared photobionts of the A-clade (96%, Fig. 3), while only some hosted an isolated subset of Trebouxia photobionts of the I-clade (4%; Fig. 3). This supports the findings of Nelson et al. [27], since on the one hand the cooler and drier habitat preference assigned to the I-clade is comparable to those found in the coastal zone of the Atacama Desert. On the other hand, the A-clade assigned Trebouxia members found in the grit crust might reflect the proposed photobiont subset after their expansion in the Cenozoic Era. However, Nelsen et al. [123] also argue that their concept does not give the full picture since some discrepancies occur (e.g. [124, 125]), which is probably due to an uneven dataset in public databases consolidated for the study as most recently outlined based on Bolivian lichens [61].

There, trebouxioid photobionts of all four clades including 16 new lineages among them were identified from 403 lichens of the Bolivian Andean vegetation providing the first study covering South American lichens and their photobionts according to the updated concepts. This data set also helped us to untangle the photobiont lineages of the National Park Pan de Azucar as great overlap could be found such as, for example, the assignment of photobionts of Heterodermia / Polyplastidium / Leucodermia and Everniopsis trulla from Bolivia and epiphytic relatives from the Atacama Desert sharing photobionts of the I-clade (Fig. 4; [61, 65]). However, those I-clade members from epiphytic lichens are only distantly related to those 4% of photobionts from the grit crust lichens and also formed a separated, unique lineages within all other I-clade photobionts (Fig. 3). The remaining 96% of the photobiont community of the grit crust could be assigned to the A-clade including photobionts from Acarospora conafii, Wetmoreana sp. and Buellia-related Caliciaceae lichens, all of which have formerly been assigned to this lineage form the same habitat [64, 99], but also from other habitats [126, 127].

Within each of those major A and I clades we could detect several new clades based on a multi-gene phylogeny such as I32 or A24 (Fig. 4), which are both the two most abundant representatives for the grit crust including the potential species candidates based on the isolates 19C (I32) and 17D2 (A24). The combination of the traditional, but still reliable, morphological comparison between isolated Trebouxia strains (Fig. 4) with the Trebouxia sequences obtained from direct amplification from thallus and Sanger sequencing (Figs. 7, 8), allowed us to demonstrate the genetic- but also morphological plasticity of the heterogeneous genus Trebouxia, which is reflected for example, in those two isolates. All potential species candidates showed significant differences in their cell sizes, chloroplast architecture, the formation of autosporangia, growth on agar and developmental stages in culture and in the intact lichen thallus (Figs. 7, 8). All of these characteristics were also identical for additional isolates sharing the same phylogeny-based clade. We also found identical phenotypic characteristics of isolates but with well-supported diverging genotypes (1B (A21); 17D2 (A24); 17C1 (A28); Fig. 4), showcasing that relying on either molecular analyses or morphological observations alone can lead to misinterpretations, as already proposed [50]. However, the unique set of isolated Trebouxia strains generated during this study allows a future formal description of the potential species candidates based on larger multi-gene phylogenies, which include for example, the cox2 locus [50].

Implications for the ecosystem

Our metabarcoding analyses highlight the presence of other green algae such as Myrmecia, or Diplosphaera (Fig. 2), both of which can be photobionts of lichens and/or free-living [128, 129]. We can exclude them from being lichenized since they were never found during the direct sequencing of grit crust lichen material. In addition, the metabarcoding derived 18S / ITS data did not indicate the presence of any lichen-forming fungi associated with e.g. Diplosphaera, such as Verrucaria or Dermatocarpon (Eurotiomycetes). However, it can not be ruled out that they might be the photobiont partner of epiphytic lichens of the environment or of grit crust lichen taxa which occur only very rarely. In addition, the metabarcoding approach we have chosen here provides a compromise between a comparable environmental study using standard primers and a lichen myco-/photobiont specific analyses due to biases introduced by primer selectivity [130]. The dominance of the green algal photobionts amongst the grit crust microbiome in general was also reflected in the high chla+b contents exceeding 800 mg m−2 at some plots with frequent fog influence close to the coast (Fig. 1B) and an average of 420 mg chla+b m−2 +- 20%. These values are comparable to well developed mesic biocrusts, e.g. from Germany (600–800 mg chla+b m−2, [70]), but exceed those found in comparable ecosystems such as the Negev or Namib Desert (< 150 mg chla+b m−2, [5, 131]). In general, the chlorophyll content of biocrusts undergoes seasonal fluctuations [5], and even higher contents can be expected for the grit crust during the winter season, where the fog water input and condensation strongly increase compared to the summer season [54, 64], during which the sampling for this study took place.

Those high chlorophyll contents can also reflect the productivity of a biocrust-ecosystem e.g. in terms of metabolic activity, net primary production and other ecosystem services which designates the grit crust as an active and principal driver of the habitat, which is otherwise mostly devoid of vegetation. More specifically, we showed that considerably the lichenized Trebouxia spp. substantially account for this, because this organism is responsible for the presence of its mycobiont partner, and therewith various other microorganisms all of which are sustained by the carbon provided by photosynthesis of the symbiotic alga.

In the context of biocrust research it has often been shown that water availability, water source and temperature as well as the resulting overall biome determine the most dominant organism(s) of the underlying biocrust [3], Fig. 9). Wherever water availability is not limiting, chlorolichens and mosses dominate, such as in temperate regions [132], but when cold temperatures dominate and desiccation occurs via frost, also cyanobacteria co-occur [9, 133]. However, wherever water availability is so limiting, that fog and dew become the major water sources, cyanobacteria disappear [134]. The grit crust as such extends the scale by defining cold deserts with low water availability and predominantly fog and/or dew as the main water source dominated by green algal lichens and green algae (Fig. 9). This theory is supported by ‘the enigmatic absence of cyanobacterial biocrusts’, which has recently been reported from the fog belt of the Namib Desert [134].

Fig. 9
figure 9

Graphical scheme of dominant biocrust-assigned cryptogams according to biomes. This scheme roughly illustrates which cryptogams are dominant in biocrusts according to temperature, water availability and -source. The grit crust of the Atacama Desert represents a biocrust type dominated by chlorolichens due to low water availability restricted to fog and / or dew as well as comparably low temperatures

Full size image

Our findings corroborate the globally unique character of the microbiome by designating the grit crust as the only known coherent soil layer with significant landscape covering impact, substantially ruled by diverse lineages of clade A of the symbiotic algal genus Trebouxia only.

Outlook

For the first time, the symbiotic green alga Trebouxia has been shown to be the dominant primary producer of a vast area represented by the grit crust ecosystem in the coastal range of the Atacama Desert. This gives a new outlining to our ecological understanding of biocrusts which are dominated by lichens. These insights will allow future studies to target the resulting consequences for associated microorganisms such as a more detailed analyses of the lichen microbiome itself (e.g. amplicon sequencing of the lichen microbiome), metabolic exchanges from the trebouxioid photobiont to other microorganisms (e.g. metabolomics of the lichen microbiome) or direct interactions between the photobionts and the lithomatrix of the stones (e.g. bioweathering). Ecophysiological experiments with the isolated photobionts will now be possible in order to determine to which extent a certain Trebouxia-lineage or -species influences the overall carbon budget or the ecological niche of the grit crust lichens. In addition, the generated data can also help to detect and compare grit crust microbiomes outside of the National Park Pan de Azucar supported by approaches which aim to identify relations between certain mycobiont-/photobiont populations and spatio-temporal distribution patterns of those.

Availability of data materials

All generated sequences were submitted to NCBI GenBank (https://www.ncbi.nlm.nih.gov/) and can be found under the accession numbers shown in Supl. Table 1. Metabarcoding data were submitted to the European Nucleotide Archive (ENA) under the project code PRJEB72845 https://www.ebi.ac.uk/ena/browser/view/PRJEB72845.

References

  1. Rodriguez-Caballero E, Belnap J, Büdel B, Crutzen PJ, Andreae MO, Pröschl U, Weber B. Dryland photoautotrophic soil surface communities endangered by global change. Nat Geosci. 2018;11:185–9.

    Article  CAS  Google Scholar 

  2. Belnap J, Büdel B, Lange OL. Biological soil crusts: characteristics and distribution. Biolo Soil Crusts Struct Funct Manag. 2003;8:3–30.

    Google Scholar 

  3. Weber B, Belnap J, Büdel B, Antoninka AJ, Barger NN, Chaudhary VB, Bowker MA. What is a biocrust? A refined, contemporary definition for a broadening research community. Biol Rev. 2022;97(5):1768–85.

    Article  PubMed  Google Scholar 

  4. Belnap J. The world at your feet: desert biological soil crusts. Front Ecol Environ. 2003;1(4):181–9.

    Article  Google Scholar 

  5. Büdel B, Darienko T, Deutschewitz K, Dojani S, Friedl T, Mohr KI, Weber B. Southern African biological soil crusts are ubiquitous and highly diverse in drylands, being restricted by rainfall frequency. Microb Ecol. 2009;57:229–47.

    Article  PubMed  Google Scholar 

  6. Lázaro R, Cantón Y, Solé-Benet A, Bevan J, Alexander R, Sancho LG, Puigdefábregas J. The influence of competition between lichen colonization and erosion on the evolution of soil surfaces in the Tabernas badlands (SE Spain) and its landscape effects. Geomorphology. 2008;102(2):252–66.

    Article  Google Scholar 

  7. Rippin M, Lange S, Sausen N, Becker B. Biodiversity of biological soil crusts from the Polar Regions revealed by metabarcoding. FEMS Microbiol Ecol. 2018;94(4):fiy06.

    Article  Google Scholar 

  8. Williams L, Colesie C, Ullmann A, Westberg M, Wedin M, Büdel B. Lichen acclimation to changing environments: photobiont switching versus climate-specific uniqueness in Psora decipiens. Ecol Evol. 2017;7(8):2560–74.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Williams L, Borchhardt N, Colesie C, Baum C, Komsic-Buchmann K, Rippin M, Büdel B. Biological soil crusts of arctic Svalbard and of Livingston Island, Antarctica. Polar Biol. 2017;40:399–411.

    Article  Google Scholar 

  10. Hagemann M, Henneberg M, Felde VJ, Drahorad SL, Berkowicz SM, Felix-Henningsen P, Kaplan A. Cyanobacterial diversity in biological soil crusts along a precipitation gradient, Northwest Negev Desert, Israel. Microb Ecol. 2015;70:219–30.

    Article  PubMed  Google Scholar 

  11. Mehda S, Muñoz-Martín MÁ, Oustani M, Hamdi-Aïssa B, Perona E, Mateo P. Microenvironmental conditions drive the differential cyanobacterial community composition of biocrusts from the Sahara desert. Microorganisms. 2021;9(3):487.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Torres-Cruz TJ, Howell AJ, Reibold RH, McHugh TA, Eickhoff MA, Reed SC. Species-specific nitrogenase activity in lichen-dominated biological soil crusts from the Colorado Plateau, USA. Plant Soil. 2018;429:113–25.

    Article  CAS  Google Scholar 

  13. Pfiz M, Loris K, Erb E, Wirth V, Küppers M. Changing patterns of lichen growth form distributions within the lichen fields of the Central Namib. Biodivers S Afr. 2010;2:33–7.

    Google Scholar 

  14. Garcia-Pichel F, Wojciechowski MF. The evolution of a capacity to build supra-cellular ropes enabled filamentous cyanobacteria to colonize highly erodible substrates. PLoS ONE. 2009;4(11):e7801.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Baumann K, Jung P, Samolov E, Lehnert LW, Büdel B, Karsten U, Leinweber P. Biological soil crusts along a climatic gradient in Chile: richness and imprints of phototrophic microorganisms in phosphorus biogeochemical cycling. Soil Biol Biochem. 2018;127:286–300.

    Article  CAS  Google Scholar 

  16. Eldridge DJ, Reed S, Travers SK, Bowker MA, Maestre FT, Ding J, Zhao Y. The pervasive and multifaceted influence of biocrusts on water in the world’s drylands. Glob Change Biol. 2020;26(10):6003–14.

    Article  Google Scholar 

  17. Wang B, Huang Y, Li N, Yao H, Yang E, Soromotin AV, An S. Initial soil formation by biocrusts: nitrogen demand and clay protection control microbial necromass accrual and recycling. Soil Biol Biochem. 2022;167:108607.

    Article  CAS  Google Scholar 

  18. Xu H, Zhang Y, Shao X, Liu N. Soil nitrogen and climate drive the positive effect of biological soil crusts on soil organic carbon sequestration in drylands: a meta-analysis. Sci Total Environ. 2022;803:150030.

    Article  PubMed  CAS  Google Scholar 

  19. DeFalco LA, Detling JK, Tracy CR, Warren SD. Physiological variation among native and exotic winter annual plants associated with microbiotic crusts in the Mojave Desert. Plant Soil. 2001;234:1–14.

    Article  CAS  Google Scholar 

  20. Miralles I, Lázaro R, Sánchez-Marañón M, Soriano M, Ortega R. Biocrust cover and successional stages influence soil bacterial composition and diversity in semiarid ecosystems. Sci Total Environ. 2020;709:134654.

    Article  PubMed  CAS  Google Scholar 

  21. Green TA, Proctor MC. Physiology of photosynthetic organisms within biological soil crusts: their adaptation, flexibility, and plasticity. Biol Soil Crusts Organ Princip Drylands. 2016;6:347–81.

    Article  Google Scholar 

  22. Lange OL, Belnap J, Reichenberger H. Photosynthesis of the cyanobacterial soil-crust lichen Collema tenax from arid lands in southern Utah, USA: role of water content on light and temperature responses of CO2 exchange. Funct Ecol. 1998;12(2):195–202.

    Article  Google Scholar 

  23. Raggio J, Pintado A, Vivas M, Sancho LG, Büdel B, Colesie C, Green TGA. Continuous chlorophyll fluorescence, gas exchange and microclimate monitoring in a natural soil crust habitat in Tabernas badlands, Almería, Spain: progressing towards a model to understand productivity. Biodivers Conserv. 2014;23:1809–26.

    Article  Google Scholar 

  24. Pombubpa N, Pietrasiak N, De Ley P, Stajich JE. Insights into dryland biocrust microbiome: geography, soil depth and crust type affect biocrust microbial communities and networks in Mojave Desert, USA. FEMS Microbiol Ecol. 2020;96(9):fiaa125.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Samolov E, Baumann K, Büdel B, Jung P, Leinweber P, Mikhailyuk T, Glaser K. Biodiversity of algae and cyanobacteria in biological soil crusts collected along a climatic gradient in Chile using an integrative approach. Microorganisms. 2020;8(7):1047.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Nelson C, Garcia-Pichel F. Beneficial cyanosphere heterotrophs accelerate establishment of cyanobacterial biocrust. Appl Environ Microbiol. 2021;87(20):e01236-e1321.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Nelson C, Giraldo-Silva A, Garcia-Pichel F. A symbiotic nutrient exchange within the cyanosphere microbiome of the biocrust cyanobacterium, Microcoleus vaginatus. ISME J. 2021;15(1):282–92.

    Article  PubMed  CAS  Google Scholar 

  28. Cardinale M, Puglia AM, Grube M. Molecular analysis of lichen-associated bacterial communities. FEMS Microbiol Ecol. 2006;57(3):484–95.

    Article  PubMed  CAS  Google Scholar 

  29. Grimm M, Grube M, Schiefelbein U, Zühlke D, Bernhardt J, Riedel K. The lichens’ microbiota, still a mystery? Front Microbiol. 2021;12:623839.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Spribille T, Tuovinen V, Resl P, Vanderpool D, Wolinski H, Aime MC, McCutcheon JP. Basidiomycete yeasts in the cortex of ascomycete macrolichens. Science. 2016;353(6298):488–92.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Vančurová L, Muggia L, Peksa O, Řídká T, Škaloud P. The complexity of symbiotic interactions influences the ecological amplitude of the host: a case study in Stereocaulon (lichenized Ascomycota). Mol Ecol. 2018;27(14):3016–33.

    Article  PubMed  Google Scholar 

  32. Hawksworth DL, Grube M. Lichens redefined as complex ecosystems. New Phytol. 2020;227:1281.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Hawksworth DL, Grube M. Reflections on lichens as ecosystems. New Phytol. 2024;241(3):972–3.

    Article  PubMed  Google Scholar 

  34. Sanders WB. The disadvantages of current proposals to redefine lichens. New Phytol. 2024;241(3):969–71.

    Article  PubMed  Google Scholar 

  35. Moya P, Molins A, Chiva S, Bastida J, Barreno E. Symbiotic microalgal diversity within lichenicolous lichens and crustose hosts on Iberian Peninsula gypsum biocrusts. Sci Rep. 2020;10(1):14060.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Ahmadjian V. Trebouxia: reflections on a perplexing and controversial lichen photobiont. In: Symbiosis: mechanisms and model systems. Dordrecht: Springer; 2001. p. 373–83.

    Chapter  Google Scholar 

  37. Rambold G, Friedl T, Beck A. Photobionts in lichens: possible indicators of phylogenetic relationships? Bryologist. 1998;8:392–7.

    Article  Google Scholar 

  38. Schweiger AH, Ullmann GM, Nürk NM, Triebel D, Schobert R, Rambold G. Chemical properties of key metabolites determine the global distribution of lichens. Ecol Lett. 2022;25(2):416–26.

    Article  PubMed  Google Scholar 

  39. Bubrick P, Galun M, Frensdorff A. Observations on free-living Trebouxia De Puymalyand Pseudotrebouxia Archibald, and evidence that both symbionts from Xanthoria parietina (L.) Th. Fr. can be found free-living in nature. New Phytol. 1984;97(3):455–62.

    Article  Google Scholar 

  40. Sanders WB, Masumoto H. Lichen algae: the photosynthetic partners in lichen symbioses. Lichenologist. 2021;53(5):347–93.

    Article  Google Scholar 

  41. Veselá V, Malavasi V, Škaloud P. A synopsis of green-algal lichen symbionts with an emphasis on their free-living lifestyle. Phycologia. 2024;6:1–22.

    Google Scholar 

  42. Yung CC, Chan Y, Lacap DC, Pérez-Ortega S, de Los Rios-Murillo A, Lee CK, Pointing SB. Characterization of chasmoendolithic community in Miers valley, McMurdo dry valleys, Antarctica. Microb Ecol. 2014;68:351–9.

    Article  PubMed  Google Scholar 

  43. Archibald PA. Trebouxia de Pulmaly (Chlorophyceae, Chlorococcales) and Pseudotrebouxia gen. nov. (Chlorophyceae, Chlorosarcinales). Phycologia. 1975;14(3):125–37.

    Article  Google Scholar 

  44. Friedl T. New aspects of the reproduction by autospores in the lichen alga Trebouxia (Microthamniales, Chlorophyta). Arch Protistenk. 1993;143(1–3):153–61.

    Article  Google Scholar 

  45. Friedl T, Besendahl A, Pfeiffer P, Bhattacharya D. The distribution of group I introns in lichen algae suggests that lichenization facilitates intron lateral transfer. Mol Phylogenet Evol. 2000;14(3):342–52.

    Article  PubMed  CAS  Google Scholar 

  46. Helms G. Taxonomy and symbiosis in associations of Physciaceae and Trebouxia; 2003.

  47. Voytsekhovich A, Beck A. Lichen photobionts of the rocky outcrops of Karadag massif (Crimean Peninsula). Symbiosis. 2016;68(1–3):9–24.

    Article  Google Scholar 

  48. Bordenave CD, Muggia L, Chiva S, Leavitt SD, Carrasco P, Barreno E. Chloroplast morphology and pyrenoid ultrastructural analyses reappraise the diversity of the lichen phycobiont genus Trebouxia (Chlorophyta). Algal Res. 2022;61:102561. https://doi.org/10.1016/j.algal.2021.102561.

    Article  Google Scholar 

  49. Leavitt SD, Kraichak E, Nelsen MP, Altermann S, Divakar PK, Alors D, Lumbsch T. Fungal specificity and selectivity for algae play a major role in determining lichen partnerships across diverse ecogeographic regions in the lichen-forming family Parmeliaceae (Ascomycota). Mol Ecol. 2015;24(14):3779–97.

    Article  PubMed  Google Scholar 

  50. Muggia L, Nelsen MP, Kirika PM, Barreno E, Beck A, Lindgren H, et al. Formally described species woefully underrepresent phylogenetic diversity in the common lichen photobiont genus Trebouxia (Trebouxiophyceae, Chlorophyta): an impetus for developing an integrated taxonomy. Mol Phylog Evol. 2020;149:106821. https://doi.org/10.1016/j.ympev.2020.106821.

    Article  Google Scholar 

  51. Barreno E, Muggia L, Chiva S, Molins A, Bordenave C, García-Breijo F, Moya P. Trebouxia lynnae sp. nov (former Trebouxia sp. TR9): biology and biogeography of an epitome lichen symbiotic microalga. Biology. 2022;11(8):1196.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. De Carolis R, Cometto A, Moya P, Barreno E, Grube M, Tretiach M, Muggia L. Photobiont diversity in lichen symbioses from extreme environments. Front Microbiol. 2022;13:809804.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Jung P, Baumann K, Emrich D, Springer A, Felde VJ, Dultz S, Leinweber P. Lichens bite the dust—a bioweathering scenario in the Atacama Desert. Iscience. 2020;23:11.

    Article  Google Scholar 

  54. Lehnert LW, Thies B, Trachte K, Achilles S, Osses P, Baumann K, Bendix J. A case study on fog/low stratus occurrence at Las Lomitas, Atacama Desert (Chile) as a water source for biological soil crusts. Aerosol Air Qual Res. 2018;18(1):254–326.

    Article  Google Scholar 

  55. Jung P, Baumann K, Lehnert LW, Samolov E, Achilles S, Schermer M, Büdel B. Desert breath—How fog promotes a novel type of soil biocenosis, forming the coastal Atacama Desert’s living skin. Geobiology. 2020;18(1):113–24.

    Article  PubMed  CAS  Google Scholar 

  56. Jung P, Lehnert LW, Bendix J, Lentendu G, Grube M, Alfaro FD, Lakatos M. The grit crust: a poly-extremotolerant microbial community from the Atacama Desert as a model for astrobiology. Front Astron Space Sci. 2022;9:1052278.

    Article  Google Scholar 

  57. Jung P, van den Brink L, Alvarado JLG. Chilean biocrusts: an impetus for national park management. Trends Microbiol. 2023a;31:544–7.

    Article  PubMed  CAS  Google Scholar 

  58. Bowker MA, Maestre FT, Mau RL. Diversity and patch-size distributions of biological soil crusts regulate dryland ecosystem multifunctionality. Ecosystems. 2013;16:923–33.

    Article  CAS  Google Scholar 

  59. Eldridge DJ, Leys JF. Wind dispersal of the vagant lichen Chondropsis semiviridis in semi-arid Eastern Australia. Aust J Bot. 1999;47(2):157–64.

    Article  Google Scholar 

  60. Rogers RW. Soil surface lichens in arid and subarid south-eastern Australia. III. The relationship between distribution and environment. Aust J Bot. 1972;20(3):301–16.

    Article  Google Scholar 

  61. Kosecka M, Kukwa M, Jabłońska A, Flakus A, Rodriguez-Flakus P, Ptach Ł, Guzow-Krzemińska B. Phylogeny and ecology of Trebouxia photobionts from Bolivian lichens. Front Microbiol. 2022;13:779784.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Medeiros ID, Mazur E, Miadlikowska J, Flakus A, Rodriguez-Flakus P, Pardo-De la Hoz CJ, Lutzoni F. Turnover of lecanoroid mycobionts and their Trebouxia photobionts along an elevation gradient in Bolivia highlights the role of environment in structuring the lichen symbiosis. Front Microbiol. 2021;12:774839.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Castillo RV, Beck A. Photobiont selectivity and specificity in Caloplaca species in a fog-induced community in the Atacama Desert, northern Chile. Fungal Biol. 2012;116(6):665–76.

    Article  Google Scholar 

  64. Jung P, Emrich D, Briegel-Williams L, Schermer M, Weber L, Baumann K, Büdel B. Ecophysiology and phylogeny of new terricolous and epiphytic chlorolichens in a fog oasis of the Atacama Desert. Microbiol Open. 2019;8(10):e894.

    Article  CAS  Google Scholar 

  65. Jung P, Baumann K, Emrich D, Schermer M, Eckhardt KU, Jandl G, Lakatos M. The dark side of orange: multiorganismic continuum dynamics within a lichen of the Atacama Desert. Mycologia. 2023;116:1–15.

    Google Scholar 

  66. Carvajal DE, Loayza AP, Squeo FA. Functional diversity and spatial association analyses at different spatial scales reveal no changes in community assembly processes along an aridity gradient in the Atacama Desert. Sci Rep. 2023;13(1):19905.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  67. Díaz FP, Latorre C, Carrasco-Puga G, Wood JR, Wilmshurst JM, Soto DC, Gutiérrez RA. Multiscale climate change impacts on plant diversity in the Atacama Desert. Glob Change Biol. 2019;25(5):1733–45.

    Article  Google Scholar 

  68. Bernhard N, Moskwa LM, Schmidt K, Oeser RA, Aburto F, Bader MY, Kühn P. Pedogenic and microbial interrelations to regional climate and local topography: new insights from a climate gradient (arid to humid) along the Coastal Cordillera of Chile. CATENA. 2018;170:335–55.

    Article  CAS  Google Scholar 

  69. Ronen R, Galun M. Pigment extraction from lichens with dimethylsulfoxide (DMSO) and estimation of chlorophyll degradation. Environ Exp Bot. 1984;24(3):239–45. https://doi.org/10.1016/0098-8472(84)90004-2.

    Article  CAS  Google Scholar 

  70. Caesar J, Tamm A, Ruckteschler N, Leifke AL, Weber B. Revisiting chlorophyll extraction methods in biological soil crusts–methodology for determination of chlorophyll a and chlorophyll a+ b as compared to previous methods. Biogeosciences. 2018;15(5):1415–24.

    Article  CAS  Google Scholar 

  71. Stoeck T, Bass D, Nebel M, Christen R, Jones MDM, Breiner H-W, Richards TA. Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water. Mol Ecol. 2010;19:326.

    Article  Google Scholar 

  72. Toju H, Tanabe AS, Yamamoto S, Sato H. High-coverage ITS primers for the DNA-based identification of ascomycetes and basidiomycetes in environmental samples. PloS One. 2012;7(7):e40863

  73. Carøe C, Bohmann K. Tagsteady: a metabarcoding library preparation protocol to avoid false assignment of sequences to samples. Mol Ecol Resour. 2020;20(6):1620–31. https://doi.org/10.1111/1755-0998.13227.

    Article  PubMed  CAS  Google Scholar 

  74. Lentendu G. DeltaMP, a flexible, reproducible and resource efficient metabarcoding amplicon pipeline for HPC (0.6) [Bash]. 2023. https://github.com/lentendu/DeltaMP.

  75. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2. https://doi.org/10.14806/ej.17.1.200.

    Article  Google Scholar 

  76. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3. https://doi.org/10.1038/nmeth.3869.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  77. Charif D, Lobry JR. SeqinR 1.0-2: a contributed package to the R project for statistical computing devoted to biological sequences retrieval and analysis. In: Bastolla U, Porto M, Roman HE, Vendruscolo M, editors. Structural approaches to sequence evolution: molecules, networks, populations. New York: Springer; 2007. p. 207–32. https://doi.org/10.1007/978-3-540-35306-5_10.

    Chapter  Google Scholar 

  78. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200. https://doi.org/10.1093/bioinformatics/btr381.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  79. Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584. https://doi.org/10.7717/peerj.2584.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Guillou L, Bachar D, Audic S, Bass D, Berney C, Bittner L, Christen R. The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy. Nucleic Acids Res. 2012;41(D1):D597–604.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Abarenkov K, Nilsson RH, Larsson KH, Taylor AF, May TW, Frøslev TG, Kõljalg U. The UNITE database for molecular identification and taxonomic communication of fungi and other eukaryotes: sequences, taxa and classifications reconsidered. Nucleic Acids Res. 2024;52(D1):D791–7.

    Article  PubMed  Google Scholar 

  82. UNITE Community. UNITE USEARCH/UTAX release for eukaryotes. Version 25.07.2023. UNITE Community. 2022. https://doi.org/10.15156/BIO/2483924.

  83. R. Core Team. R: A language and environment for statistical computing, 2013. R Foundation for Statistical Computing, Vienna, Austria. 2014.

  84. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin P, O’Hara RB, Simpson GL, Solymos P. Vegan: community ecology package. R Package Version 2.5-7; 2021. p. 2020.

  85. Pölme S, Abarenkov K, Nilsson RH, Lindahl BD, Clemmensen KE, Kauserud H, Tedersoo L. FungalTraits: a user friendly traits database of fungi and fungus-like stramenopiles, vol 105, 116; 2021.

  86. Bischoff HW, Bold HC, Physiological studies: IV. In Some soil algae from enchanted rock and related algal species. Austin: University of Texas: 1963; p. 6318.

  87. Jung P, Briegel-Williams L, Werner L, Jost E, Schultz M, Nürnberg DJ, Lakatos M. A direct PCR approach with low-biomass insert opens new horizons for molecular sciences on cryptogam communities. Appl Environ Microbiol 2024;90(3):e00024-24

  88. Gardes M, Bruns TD. ITS primers with enhanced specificity for basidiomycetes-application to the identification of mycorrhizae and rusts. Mol Ecol. 1993;2(2):113–8. https://doi.org/10.1111/j.1365-294X.1993.tb00005.x.

    Article  PubMed  CAS  Google Scholar 

  89. Friedl T, Rokitta C. Species relationships in the Lichen Alga Trebouxia (Chlorophyta, Trebouxiophyceae): molecular phylogenetic analyses of nuclear-encoded large subunit rRNA gene sequences. Symbiosis. 1997;3:125–48.

    Google Scholar 

  90. Kroken S, Taylor JW. Phylogenetic species, reproductive mode, and specificity of the green alga Trebouxia forming lichens with the fungal genus Letharia. Bryol. 2000;103:645–60.

    Article  CAS  Google Scholar 

  91. Xu M, De Boer H, Olafsdottir ES, Omarsdottir S, Heidmarsson S. Phylogenetic diversity of the lichenized algal genus Trebouxia (Trebouxiophyceae, Chlorophyta): a new lineage and novel insights from fungal-algal association patterns of Icelandic cetrarioid lichens (Parmeliaceae, Ascomycota). Bot J Linn Soc. 2020;194(4):460–8.

    Article  Google Scholar 

  92. Tamura K, Stecher G, Kumar S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol Biol Evol. 2021;38(7):3022-3027. https://doi.org/10.1093/molbev/msab120

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  93. Lemoine F, Correia D, Lefort V, Doppelt-Azeroual O, Mareuil F, Cohen-Boulakia S, Gascuel O. NGPhylogeny. fr: new generation phylogenetic services for non-specialists. Nucleic Acids Res. 2019;47(W1):W260–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  94. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003;19(12):1572-1574

    Article  PubMed  CAS  Google Scholar 

  95. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucl Acids Res 2021:49(W1):W293-W296

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  96. Sadowsky A, Mettler-Altmann T, Ott S. Metabolic response to desiccation stress in strains of green algal photobionts (Trebouxia) from two Antarctic lichens of southern habitats. Phycologia. 2016;55:703–14. https://doi.org/10.2216/15-127.1.

    Article  CAS  Google Scholar 

  97. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9:676–82.

    Article  PubMed  CAS  Google Scholar 

  98. Tedersoo L, Bahram M, Zinger L, Nilsson RH, Kennedy PG, Yang T, Mikryukov V. Best practices in metabarcoding of fungi: from experimental design to results. Mol Ecol. 2022;31(10):2769–95.

    Article  PubMed  Google Scholar 

  99. Jung P, Werner L, Briegel-Williams L, Emrich D, Lakatos M. Roccellinastrum, Cenozosia and Heterodermia: ecology and phylogeny of fog lichens and their photobionts from the coastal Atacama Desert. MycoKeys. 2023;98:317.

    Article  PubMed  PubMed Central  Google Scholar 

  100. Nimis P, Martellos S. Keys to the lichens of Italy. I. Terricolous species. In: Goliardiche, editor; 2004. pp. 1–341.

  101. Green TA, Pintado A, Raggio J, Sancho LG. The lifestyle of lichens in soil crusts. Lichenologist. 2018;50(3):397–410.

    Article  Google Scholar 

  102. Rosentreter R, Eldridge DJ, Westberg M, Williams L, Grube M. Structure, composition, and function of biocrust lichen communities. Biol Soil Crusts Organ Principle Drylands. 2016;26:121–38.

    Article  Google Scholar 

  103. Zhang XC, Li JY, Liu JL, Yuan CX, Li YN, Liu BR, Yan XF. Temporal shifts in cyanobacterial diversity and their relationships to different types of biological soil crust in the southeastern Tengger Desert. Rhizosphere. 2021;17:100322.

    Article  Google Scholar 

  104. Pushkareva E, Elster J, Holzinger A, Niedzwiedz S, Becker B. Biocrusts from Iceland and Svalbard: Does microbial community composition differ substantially? Front Microbiol. 2022;13:1048522.

    Article  PubMed  PubMed Central  Google Scholar 

  105. Prieto M, Wedin M. Phylogeny, taxonomy and diversification events in the Caliciaceae. Fungal Divers. 2017;82(1):221–38.

    Article  Google Scholar 

  106. Prieto M, Baloch E, Tehler A, Wedin M. Mazaedium evolution in the Ascomycota (Fungi) and the classification of mazaediate groups of formerly unclear relationship. Cladistics. 2013;29(3):296–308.

    Article  PubMed  Google Scholar 

  107. Selva SB. The calicioid lichens and fungi of the Acadian Forest ecoregion of northeastern North America, I. New species and range extensions. Bryologist. 2013;116(3):248–56.

    Article  Google Scholar 

  108. Dunai TJ, López GAG, Juez-Larré J. Oligocene-Miocene age of aridity in the Atacama Desert revealed by exposure dating of erosion-sensitive landforms. Geology. 2005;33(4):321–4.

    Article  Google Scholar 

  109. Sun T, Bao H, Reich M, Hemming SR. More than ten million years of hyper-aridity recorded in the Atacama Gravels. Geochim Cosmochim Acta. 2018;227:123–32.

    Article  CAS  Google Scholar 

  110. Hartley AJ, Chong G, Houston J, Mather AE. 150 million years of climatic stability: evidence from the Atacama Desert, northern Chile. J Geol Soc. 2005;162(3):421–4.

    Article  Google Scholar 

  111. Belnap J, Phillips SL, Miller ME. Response of desert biological soil crusts to alterations in precipitation frequency. Oecologia. 2004;141:306–16.

    Article  PubMed  Google Scholar 

  112. Reed SC, Coe KK, Sparks JP, Housman DC, Zelikova TJ, Belnap J. Changes to dryland rainfall result in rapid moss mortality and altered soil fertility. Nat Clim Change. 2012;2(10):752. https://doi.org/10.1038/nclimate1596.

    Article  CAS  Google Scholar 

  113. del Río C, Lobos F, Siegmund A, Tejos C, Osses P, Huaman Z, García JL. GOFOS, ground optical fog observation system for monitoring the vertical stratocumulus-fog cloud distribution in the coast of the Atacama Desert, Chile. J Hydrol. 2021;597:126190.

    Article  Google Scholar 

  114. Lobos-Roco F, Suárez F, Aguirre-Correa F, Keim K, Aguirre I, Vargas C, del Río C. Understanding inland fog and dew dynamics for assessing potential non-rainfall water use in the Atacama. J Arid Environ. 2024;221:105125.

    Article  Google Scholar 

  115. Gessler NN, Egorova AS, Belozerskaya TA. Melanin pigments of fungi under extreme environmental conditions. Appl Biochem Microbiol. 2014;50:105–13.

    Article  CAS  Google Scholar 

  116. Mafole TC, Solhaug KA, Minibayeva FV, Beckett RP. Occurrence and possible roles of melanic pigments in lichenized ascomycetes. Fungal Biol Rev. 2019;33(3–4):159–65.

    Article  Google Scholar 

  117. Lange OL, Meyer A, Zellner H, Heber U. Photosynthesis and water relations of lichen soil crusts: field measurements in the coastal fog zone of the Namib Desert. Funct Ecol. 1994;8(2):253–64.

    Article  Google Scholar 

  118. Büdel B, Vivas M, Lange OL. Lichen species dominance and the resulting photosynthetic behavior of Sonoran Desert soil crust types (Baja California, Mexico). Ecol Process. 2013;2:1–9.

    Article  Google Scholar 

  119. Stanton DE, Ormond A, Koch NM, Colesie C. Lichen ecophysiology in a changing climate. Am J Bot. 2023;110(2):e16131.

    Article  PubMed  Google Scholar 

  120. Lange OL, Kilian E, Ziegler H. Water vapor uptake and photosynthesis of lichens: performance differences in species with green and blue-green algae as phycobionts. Oecologia. 1986;71(1):104–10.

    Article  PubMed  CAS  Google Scholar 

  121. Nelsen MP, Lücking R, Boyce CK, Lumbsch HT, Ree RH. No support for the emergence of lichens prior to the evolution of vascular plants. Geobiology. 2020;18(1):3–13.

    Article  PubMed  Google Scholar 

  122. Nelsen MP, Leavitt SD, Heller K, Muggia L, Lumbsch HT. Contrasting patterns of climatic niche divergence in Trebouxia—a Clade of lichen-forming algae. Front Microbiol. 2022;13:64.

    Article  Google Scholar 

  123. Nelsen MP, Leavitt SD, Heller K, Muggia L, Lumbsch HT. Macroecological diversification and convergence in a clade of keystone symbionts. FEMS Microbiol Ecol. 2021;97(6):fiab072.

    Article  PubMed  CAS  Google Scholar 

  124. Molins A, Chiva S, Calatayud Á, Marco F, García-Breijo F, Armiñana JR, et al. Multidisciplinary approach to describe Trebouxia diversity within lichenized fungi Buellia zoharyi from the Canary Islands. Symbiosis. 2020;82:19–34. https://doi.org/10.1007/s13199-020-00722-8.

    Article  CAS  Google Scholar 

  125. Singh G, Kukwa M, Dal Grande F, Łubek A, Otte J, Schmitt I. A glimpse into genetic diversity and symbiont interaction patterns in lichen communities from areas with different disturbance histories in Białowieża forest. Poland Microorgan. 2019;7:335. https://doi.org/10.3390/microorganisms7090335.

    Article  CAS  Google Scholar 

  126. Pérez-Ortega S, Ortiz-Álvarez R, Green TA, de los Ríos A. Lichen myco- and photobiont diversity and their relationships at the edge of life (McMurdo Dry Valleys, Antarctica). FEMS Microbiol Ecol. 2012;82:429–48. https://doi.org/10.1007/s00300-020-02754-8.

    Article  PubMed  Google Scholar 

  127. Wagner M, Bathke AC, Cary SC, Green TA, Junker RR, Trutschnig W, et al. Myco-and photobiont associations in crustose lichens in the McMurdo Dry Valleys (Antarctica) reveal high differentiation along an elevational gradient. Polar Biol. 2020;43:1967–83. https://doi.org/10.1007/s00300-020-02754-8.

    Article  Google Scholar 

  128. Pino-Bodas R, Blázquez M, de los Ríos A, Pérez-Ortega S. Myrmecia, not asterochloris, is the main photobiont of Cladonia Subturgida (Cladoniaceae, Lecanoromycetes). J Fungi. 2023;9(12):1160.

    Article  CAS  Google Scholar 

  129. Thüs H, Muggia L, Pérez-Ortega S, Favero-Longo SE, Joneson S, O’Brien H, Gueidan C. Revisiting photobiont diversity in the lichen family Verrucariaceae (Ascomycota). Eur J Phycol. 2011;46(4):399–415.

    Article  Google Scholar 

  130. Wright B, Clair LLS, Leavitt SD. Is targeted community DNA metabarcoding suitable for biodiversity inventories of lichen-forming fungi? Ecol Ind. 2019;98:812–20.

    Article  CAS  Google Scholar 

  131. Kidron GJ, Li XR, Jia RL, Gao YH, Zhang P. Assessment of carbon gains from biocrusts inhabiting a dunefield in the Negev Desert. Geoderma. 2015;253:102–10.

    Article  Google Scholar 

  132. Corbin JD, Thiet RK. Temperate biocrusts: mesic counterparts to their better-known dryland cousins. Front Ecol Environ. 2020;18(8):456–64.

    Article  Google Scholar 

  133. Jung P, Briegel-Williams L, Simon A, Thyssen A, Büdel B. Uncovering biological soil crusts: carbon content and structure of intact Arctic, Antarctic and alpine biological soil crusts. Biogeosciences. 2018;15(4):1149–60.

    Article  CAS  Google Scholar 

  134. Kidron GJ. The enigmatic absence of cyanobacterial biocrusts from the Namib fog belt: Do dew and fog hold the key? Flora. 2019;257:151416.

    Article  Google Scholar 

Download references

Acknowledgements

The authors want to thank the CONAF (Corporación Nacional Forestal) rangers of the National Park Pan de Azúcar for their ongoing support during field work and interpretation of the ecosystem. We also want to thank BioRender and the assistance of the Core Facility BioSupraMol supported by the German Research Council for confocal microscopy.

Funding

Open Access funding enabled and organized by Projekt DEAL. PJ, LW and EJ were funded by the German Research Council (Grit Life JU 3228/1-1), LBW by internal funding of the University of Applied Sciences Kaiserslautern and by the Federal Ministry of Education and Research (W2V-Strategy2Value, 03WIR4502A). ML was funded by Ministry of Science and Health Rhineland-Palatinate (724-0116#2021/004-1501 15405) and by the Federal Ministry of Education and Research (W2V-Strategy2Value, 03WIR4502A) as well as RB (W2V-NovelBioChem, 03WIR4505B). GL was funded by the Swiss National Science Foundation (grant n° 182531).

Author information

Authors and Affiliations

  1. Department of Integrative Biotechnology, University of Applied Sciences Kaiserslautern, Pirmasens, Germany

    Patrick Jung, Rebekah Brand, Laura Briegel-Williams, Lina Werner, Emily Jost & Michael Lakatos

  2. Laboratory of Soil Biodiversity, Université de Neuchâtel, Neuchâtel, Switzerland

    Guillaume Lentendu

  3. Soil Science and Environment Group, Changins, HES-SO University of Applied Sciences and Arts Western Switzerland, Nyon, Switzerland

    David Singer

  4. Institute for Experimental Physics, Freie Universität Berlin, Berlin, Germany

    Rujuta Athavale & Dennis J. Nürnberg

  5. Dahlem Centre of Plant Sciences, Freie Universität Berlin, Berlin, Germany

    Dennis J. Nürnberg

  6. GEMA Center for Genomics, Ecology and Environment, Universidad Mayor, Santiago, Chile

    Fernando D. Alfaro

  7. Biology, Rhineland-Palatinate Technical University Kaiserslautern Landau, Kaiserslautern, Germany

    Burkhard Büdel

Authors
  1. Patrick Jung
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Rebekah Brand
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Laura Briegel-Williams
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Lina Werner
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Emily Jost
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Guillaume Lentendu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. David Singer
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Rujuta Athavale
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Dennis J. Nürnberg
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Fernando D. Alfaro
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Burkhard Büdel
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Michael Lakatos
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

PJ, LBW and ML took samples, generated funding and supervised the work. PJ wrote the manuscript and prepared the figures while EJ and LW conducted lab work. GL performed bioinformatic analyses as well as DS and RB performed the statistical analyses of the metabarcoding dataset. RA and DJN performed confocal microscopy and analyzed images. All authors contributed and edited the manuscript.

Corresponding author

Correspondence to Patrick Jung.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jung, P., Brand, R., Briegel-Williams, L. et al. The symbiotic alga Trebouxia fuels a coherent soil ecosystem on the landscape scale in the Atacama Desert. Environmental Microbiome 19, 59 (2024). https://doi.org/10.1186/s40793-024-00601-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40793-024-00601-5

Keywords

Environmental Microbiome

ISSN: 2524-6372