Lichen holobionts show compositional structure along elevation.

Holobionts are dynamic ecosystems that may respond to abiotic drivers with compositional changes. Uncovering elevational diversity patterns within these microecosystems can further our understanding of community-environment interactions. Here, we assess how the major components of lichen holobionts-fungal hosts, green algal symbionts, and the bacterial community-collectively respond to an elevational gradient. We analyse populations of two lichen symbioses, Umbilicaria pustulata and U. hispanica, along an elevational gradient spanning 2100 altitudinal metres and covering three major biomes. Our study shows (i) discontinuous genomic variation in fungal hosts with one abrupt genomic differentiation within each of the two host species, (ii) altitudinally structured bacterial communities with pronounced turnover within and between hosts, and (iii) altitude-specific presence of algal symbionts. Alpha diversity of bacterial communities decreased with increasing elevation. A marked turnover in holobiont diversity occurred across two altitudinal belts: at 11°C-13°C average annual temperature (here: 800-1200 m a.s.l.), and at 7°C-9°C average annual temperature (here: 1500-1800 m a.s.l.). The two observed zones mark a clustering of distribution limits and community shifts. The three ensuing altitudinal classes, that is, the most frequent combinations of species in holobionts, approximately correspond to the Mediterranean, cool-temperate, and alpine climate zones. We conclude that multitrophic microecosystems, such as lichen holobionts, respond with concerted compositional changes to climatic factors that also structure communities of macroorganisms, for example, vascular plants.

Metatranscriptomics reveals contrasting effects of elevation on the activity of bacteria and bacterial viruses in soil.

Soil microbial diversity affects ecosystem functioning and global biogeochemical cycles. Soil bacterial communities catalyse a diversity of biogeochemical reactions and have thus sparked considerable scientific interest. One driver of bacterial community dynamics in natural ecosystems has so far been largely neglected: the predator-prey interactions between bacterial viruses (bacteriophages) and bacteria. To generate ground level knowledge on environmental drivers of these particular predator-prey dynamics, we propose an activity-based ecological framework to simultaneous capture community dynamics of bacteria and bacteriophages in soils. An ecological framework and specifically the analyses of community dynamics across latitudinal and elevational gradients have been widely used in ecology to understand community-wide responses of innumerable taxa to environmental change, in particular to climate. Here, we tested the hypothesis that the activity of bacteria and bacteriophages codeclines across an elevational gradient. We used metatranscriptomics to investigate bacterial and bacteriophage activity patterns at five sites across 400 elevational metres in the Swiss Alps in 2015 and 2017. We found that metabolic activity (transcription levels) of bacteria declined significantly with increasing elevation, but activity of bacteriophages did not. We showed that bacteriophages are consistently active in soil along the entire gradient, making bacteriophage activity patterns divergent from that of their putative bacterial prey. Future efforts will be necessary to link the environment-activity relationship to predator-prey dynamics, and to understand the magnitude of viral contributions to carbon, nitrogen and phosphorus cycling when infection causes bacterial cell death, a process that may represent an overlooked component of soil biogeochemical cycles.

Virus diversity in metagenomes of a lichen symbiosis (Umbilicaria phaea): complete viral genomes, putative hosts and elevational distributions.

Viruses can play critical roles in symbioses by initiating horizontal gene transfer, affecting host phenotypes, or expanding their host's ecological niche. However, knowledge of viral diversity and distribution in symbiotic organisms remains elusive. Here we use deep-sequenced metagenomic DNA (PacBio Sequel II; two individuals), paired with a population genomics approach (Pool-seq; 11 populations, 550 individuals) to understand viral distributions in the lichen Umbilicaria phaea. We assess (i) viral diversity in lichen thalli, (ii) putative viral hosts (fungi, algae, bacteria) and (iii) viral distributions along two replicated elevation gradients. We identified five novel viruses, showing 28%-40% amino acid identity to known viruses. They tentatively belong to the families Caulimoviridae, Myoviridae, Podoviridae and Siphoviridae. Our analysis suggests that the Caulimovirus is associated with green algal photobionts (Trebouxia) of the lichen, and the remaining viruses with bacterial hosts. We did not detect viral sequences in the mycobiont. Caulimovirus abundance decreased with increasing elevation, a pattern reflected by a specific algal lineage hosting this virus. Bacteriophages showed population-specific patterns. Our work provides the first comprehensive insights into viruses associated with a lichen holobiont and suggests an interplay of viral hosts and environment in structuring viral distributions.

Climate-specific biosynthetic gene clusters in populations of a lichen-forming fungus.

Natural products can contribute to abiotic stress tolerance in plants and fungi. We hypothesize that biosynthetic gene clusters (BGCs), the genomic elements that underlie natural product biosynthesis, display structured differences along elevation gradients. We analysed biosynthetic gene variation in natural populations of the lichen-forming fungus Umbilicaria pustulata. We collected a total of 600 individuals from the Mediterranean and cold-temperate climates. Population genomic analyses indicate that U. pustulata contains three clusters that are highly differentiated between the Mediterranean and cold-temperate populations. One entire cluster is exclusively present in cold-temperate populations, and a second cluster is putatively dysfunctional in all cold-temperate populations. In the third cluster variation is fixed in all cold-temperate populations due to hitchhiking. In these two clusters the presence of consistent allele frequency differences among replicate populations/gradients suggests that selection rather than drift is driving the pattern. We advocate that the landscape of fungal biosynthetic genes is shaped by both positive and hitchhiking selection. We demonstrate, for the first time, the presence of climate-associated BGCs and BGC variations in lichen-forming fungi. While the associated secondary metabolites of the candidate clusters are presently unknown, our study paves the way for targeted discovery of natural products with ecological significance.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata.

BACKGROUND: Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution. RESULTS: A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting. CONCLUSIONS: The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.

Biosynthetic Gene Content of the 'Perfume Lichens' Evernia prunastri and Pseudevernia furfuracea.

Lichen-forming fungi produce a vast number of unique natural products with a wide variety of biological activities and human uses. Although lichens have remarkable potential in natural product research and industry, the molecular mechanisms underlying the biosynthesis of lichen metabolites are poorly understood. Here we use genome mining and comparative genomics to assess biosynthetic gene clusters and their putative regulators in the genomes of two lichen-forming fungi, which have substantial commercial value in the perfume industry, Evernia prunastri and Pseudevernia furfuracea. We report a total of 80 biosynthetic gene clusters (polyketide synthases (PKS), non-ribosomal peptide synthetases and terpene synthases) in E. prunastri and 51 in P. furfuracea. We present an in-depth comparison of 11 clusters, which show high homology between the two species. A ketosynthase (KS) phylogeny shows that biosynthetic gene clusters from E. prunastri and P. furfuracea are widespread across the Fungi. The phylogeny includes 15 genomes of lichenized fungi and all fungal PKSs with known functions from the MIBiG database. Phylogenetically closely related KS domains predict not only similar PKS architecture but also similar cluster architecture. Our study highlights the untapped biosynthetic richness of lichen-forming fungi, provides new insights into lichen biosynthetic pathways and facilitates heterologous expression of lichen biosynthetic gene clusters.

Environment and host identity structure communities of green algal symbionts in lichens.

An understanding of how biotic interactions shape species' distributions is central to predicting host-symbiont responses under climate change. Switches to locally adapted algae have been proposed to be an adaptive strategy of lichen-forming fungi to cope with environmental change. However, it is unclear how lichen photobionts respond to environmental gradients, and whether they play a role in determining the fungal host's upper and lower elevational limits. Deep-coverage Illumina DNA metabarcoding was used to track changes in the community composition of Trebouxia algae associated with two phylogenetically closely related, but ecologically divergent fungal hosts along a steep altitudinal gradient in the Mediterranean region. We detected the presence of multiple Trebouxia species in the majority of thalli. Both altitude and host genetic identity were strong predictors of photobiont community assembly in these two species. The predominantly clonally dispersing fungus showed stronger altitudinal structuring of photobiont communities than the sexually reproducing host. Elevation ranges of the host were not limited by the lack of compatible photobionts. Our study sheds light on the processes guiding the formation and distribution of specific fungal-algal combinations in the lichen symbiosis. The effect of environmental filtering acting on both symbiotic partners appears to shape the distribution of lichens.

Antimicrobial-producing Pseudoalteromonas from the marine environment of Panama shows a high phylogenetic diversity and clonal structure.

Pseudoalteromonas is a genus of marine bacteria often found in association with other organisms. Although several studies have examined Pseudoalteromonas diversity and their antimicrobial activity, its diversity in tropical environments is largely unexplored. We investigated the diversity of Pseudoalteromonas in marine environments of Panama using a multilocus phylogenetic approach. Furthermore we tested their antimicrobial capacity and evaluated the effect of recombination and mutation in shaping their phylogenetic relationships. The reconstruction of clonal relationships among 78 strains including 15 reference Pseudoalteromonas species revealed 43 clonal lineages, divided in pigmented and non-pigmented strains. In total, 39 strains displayed moderate to high activity against Gram-positive and Gram-negative bacteria and fungi. Linkage disequilibrium analyses showed that the Pseudoalteromonas strains of Panama have a highly clonal structure and that, although present, recombination is not frequent enough to break the association among alleles. This clonal structure is in contrast to the high rates of recombination generally reported for aquatic and marine bacteria. We propose that this structure is likely due to the symbiotic association with marine invertebrates of most strains analyzed. Our results also show that there are several putative new species of Pseudoalteromonas in Panama to be described.

Adaptive differentiation coincides with local bioclimatic conditions along an elevational cline in populations of a lichen-forming fungus.

BACKGROUND: Many fungal species occur across a variety of habitats. Particularly lichens, fungi forming symbioses with photosynthetic partners, have evolved remarkable tolerances for environmental extremes. Despite their ecological importance and ubiquity, little is known about the genetic basis of adaption in lichen populations. Here we studied patterns of genome-wide differentiation in the lichen-forming fungus Lasallia pustulata along an altitudinal gradient in the Mediterranean region. We resequenced six populations as pools and identified highly differentiated genomic regions. We then detected gene-environment correlations while controlling for shared population history and pooled sequencing bias, and performed ecophysiological experiments to assess fitness differences of individuals from different environments. RESULTS: We detected two strongly differentiated genetic clusters linked to Mediterranean and temperate-oceanic climate, and an admixture zone, which coincided with the transition between the two bioclimates. High altitude individuals showed ecophysiological adaptations to wetter and more shaded conditions. Highly differentiated genome regions contained a number of genes associated with stress response, local environmental adaptation, and sexual reproduction. CONCLUSIONS: Taken together our results provide evidence for a complex interplay between demographic history and spatially varying selection acting on a number of key biological processes, suggesting a scenario of ecological speciation.

Fungal-algal association patterns in lichen symbiosis linked to macroclimate.

Both macroclimate and evolutionary events may influence symbiont association and diversity patterns. Here we assess how climatic factors and evolutionary events shape fungal-algal association patterns in the widely distributed lichen-forming fungal genus Protoparmelia. Multilocus phylogenies of fungal and algal partners were generated using 174 specimens. Coalescent-based species delimitation analysis suggested that 23 fungal hosts are associating with 20 algal species. Principal component analysis (PCA) was performed to infer how fungal-algal association patterns varied with climate. Fungi associated with one to three algal partners whereas algae accepted one to five fungal partners. Both fungi and algae were more specific, associating with fewer partners, in the warmer climates. Interaction with more than one partner was more frequent in cooler climates for both the partners. Cophylogenetic analyses suggest congruent fungal-algal phylogenies. Host switch was a more common event in warm climates, whereas failure of the photobiont to diverge with its fungal host was more frequent in cooler climates. We conclude that both environmental factors and evolutionary events drive fungal and algal evolution in Protoparmelia. The processes leading to phylogenetic congruence of fungi and algae are different in different macrohabitats in our study system. Hence, closely related species inhabiting diverse habitats may follow different evolutionary pathways.

Taxonomic revision and species delimitation of coccoid green algae currently assigned to the genus Dictyochloropsis (Trebouxiophyceae, Chlorophyta).

Coccoid green algae traditionally classified in Dictyochloropsis have a complex, reticulate chloroplast, when mature, without a pyrenoid. They occupy remarkably diverse ecological niches as free-living organisms or in association with lichen-forming fungi and were recently shown to form two distinct lineages within Trebouxiophyceae. We used a polyphasic approach to revise the taxonomy of the genus. Based on phylogenetic analysis of the 18S rRNA gene, and detailed morphological investigation using comparative conventional light and confocal microscopy, we have assigned these lineages to two genera, Dictyochloropsis and Symbiochloris gen. nov. We have reconsidered the diagnostic generic features as follows: Dictyochloropsis comprises only free-living algae with a reticulate chloroplast, forming lobes in a parallel arrangement at some ontogenetic stages, and which reproduce only by means of autospores. This agrees with Geitler's original diagnosis of Dictyochloropsis, but not with the later emendation by Tschermak-Woess. Consequently, the species of Dictyochloropsis sensu Tschermak-Woess are assigned to Symbiochloris, with new combinations proposed. Symbiochloris encompasses free-living and/or lichenized algae with lobed chloroplasts and that reproduce by forming zoospores characterized by two subapical isokont flagella that emerge symmetrically near the flattened apex. In addition, using coalescent-based approaches, morphological characters and secondary structure of ITS transcripts, we inferred species boundaries and taxonomic relationships within the newly proposed genera. Two species of Dictyochloropsis and nine species of Symbiochloris are delimited, including the newly described species D. asterochloroides, S. handae, S. tropica, and S. tschermakiae. Our results further support the non-monophyly of autosporine taxa within Trebouxiophyceae.

Evolution of complex symbiotic relationships in a morphologically derived family of lichen-forming fungi.

We studied the evolutionary history of the Parmeliaceae (Lecanoromycetes, Ascomycota), one of the largest families of lichen-forming fungi with complex and variable morphologies, also including several lichenicolous fungi. We assembled a six-locus data set including nuclear, mitochondrial and low-copy protein-coding genes from 293 operational taxonomic units (OTUs). The lichenicolous lifestyle originated independently three times in lichenized ancestors within Parmeliaceae, and a new generic name is introduced for one of these fungi. In all cases, the independent origins occurred c. 24 million yr ago. Further, we show that the Paleocene, Eocene and Oligocene were key periods when diversification of major lineages within Parmeliaceae occurred, with subsequent radiations occurring primarily during the Oligocene and Miocene. Our phylogenetic hypothesis supports the independent origin of lichenicolous fungi associated with climatic shifts at the Oligocene-Miocene boundary. Moreover, diversification bursts at different times may be crucial factors driving the diversification of Parmeliaceae. Additionally, our study provides novel insight into evolutionary relationships in this large and diverse family of lichen-forming ascomycetes.

Molecular phylogeny and symbiotic selectivity of the green algal genus Dictyochloropsis s.l. (Trebouxiophyceae): a polyphyletic and widespread group forming photobiont-mediated guilds in the lichen family Lobariaceae.

Dictyochloropsis s.l. is an ecologically important, common but little-studied genus of green algae. Here, we examined the diversity and host selectivity of algae attributed to this genus at both species-to-species and species-to-community levels. We conducted a molecular investigation of 15 cultured strains and several lichen photobionts, using 18S rRNA, rbcL and ITS sequence data. We further used seven alga-specific microsatellite markers to study algal sharing among fungi of the family Lobariaceae in two populations in Madeira and Taiwan (454 lichens). We found that the genus Dictyochloropsis s.l. is polyphyletic. Dictyochloropsis clade 1 comprises only free-living algae whereas Dictyochloropsis clade 2 includes lichenized algae as well as free-living algae. Fungal selectivity towards algae belonging to Dictyochloropsis clade 2 is high. Selectivity varies geographically, with photobionts being restricted to a single region. Finally, we showed that Dictyochloropsis clade 2 individuals are shared among different fungal hosts in communities of lichens of the Lobariaceae. As for other green algal lineages, there is a high amount of cryptic diversity in Dictyochloropsis. Furthermore, co-evolution between Dictyochloropsis clade 2 algae and representatives of the Lobariaceae is manifested at the community level, with several unrelated fungal species being horizontally connected by shared photobiont clones.

Insights into intrathalline genetic diversity of the cosmopolitan lichen symbiotic green alga Trebouxia decolorans Ahmadjian using microsatellite markers.

Trebouxia decolorans is a widespread and common symbiotic green alga that is found in association with different species of lichen-forming fungi. By applying T. decolorans-specific microsatellite markers, we investigated the within-thallus diversity of T. decolorans in thalli of Xanthoria parietina and Anaptychia ciliaris. We found several algal strains in most of the thalli of both hosts. High genetic differentiation among thalli suggests that algal diversity is generated de novo via mutation in both fungal hosts. Rarefied allelic richness of the algae was higher in thalli of X. parietina. Our results indicate that in X. parietina intrathalline algal diversity is additionally created by environmental uptake of algae either at the start of the symbiotic association or during the lifetime of the thallus. This study indicates that promiscuous host-symbiont associations in lichen symbioses with Trebouxia spp. may be more common than currently recognized.

Integrating coalescent and phylogenetic approaches to delimit species in the lichen photobiont Trebouxia.

The accurate assessment of species boundaries in symbiotic systems is a prerequisite for the study of speciation, co-evolution and selectivity. Many studies have shown the high genetic diversity of green algae from the genus Trebouxia, the most common photobiont of lichen-forming fungi. However, the phylogenetic relationships, and the amount of cryptic diversity of these algae are still poorly understood, and an adequate species concept for trebouxiophycean algae is still missing. In this study we used a multifaceted approach based on coalescence (GMYC, STEM) and phylogenetic relationships to assess species boundaries in the trebouxioid photobionts of the lichen-forming fungus Lasallia pustulata. We further investigated whether putative species of Trebouxia found in L. pustulata are shared with other lichen-forming fungi. We found that L. pustulata is associated with at least five species of Trebouxia and most of them are shared with other lichen-forming fungi, showing different patterns of species-to-species and species-to-community interactions. We also show that one of the putative Trebouxia species is found exclusively in association with L. pustulata and is restricted to thalli from localities with Mediterranean microclimate. We suggest that the species delimitation method presented in this study is a promising tool to address species boundaries within the heterogeneous genus Trebouxia.

Biotic interactions outweigh abiotic factors as drivers of bark microbial communities in Central European forests.

Bark surfaces are extensive areas within forest ecosystems, which provide an ideal habitat for microbial communities, through their longevity and seasonal stability. Here we provide a comprehensive account of the bark surface microbiome of living trees in Central European forests, and identify drivers of diversity and community composition. We examine algal, fungal, and bacterial communities and their interactions using metabarcoding on samples from over 750 trees collected in the Biodiversity Exploratories in northern, central, and southern Germany. We show that mutual biotic influence is more important than the abiotic environment with regard to community composition, whereas abiotic conditions and geography are more important for alpha diversity. Important abiotic factors are the relative humidity and light availability, which decrease the algal and bacterial alpha diversity but strongly increase fungal alpha diversity. In addition, temperature is important in shaping the microbial community, with higher temperature leading to homogeneous communities of dominant fungi, but high turnover in bacterial communities. Changes in the community dissimilarity of one organismal group occur in close relation to changes in the other two, suggesting that there are close interactions between the three major groups of the bark surface microbial communities, which may be linked to beneficial exchange. To understand the functioning of the forest microbiome as a whole, we need to further investigate the functionality of interactions within the bark surface microbiome and combine these results with findings from other forest habitats such as soil or canopy.

Phylogenomics reveals the evolutionary origins of lichenization in chlorophyte algae.

Mutualistic symbioses have contributed to major transitions in the evolution of life. Here, we investigate the evolutionary history and the molecular innovations at the origin of lichens, which are a symbiosis established between fungi and green algae or cyanobacteria. We de novo sequence the genomes or transcriptomes of 12 lichen algal symbiont (LAS) and closely related non-symbiotic algae (NSA) to improve the genomic coverage of Chlorophyte algae. We then perform ancestral state reconstruction and comparative phylogenomics. We identify at least three independent gains of the ability to engage in the lichen symbiosis, one in Trebouxiophyceae and two in Ulvophyceae, confirming the convergent evolution of the lichen symbioses. A carbohydrate-active enzyme from the glycoside hydrolase 8 (GH8) family was identified as a top candidate for the molecular-mechanism underlying lichen symbiosis in Trebouxiophyceae. This GH8 was acquired in lichenizing Trebouxiophyceae by horizontal gene transfer, concomitantly with the ability to associate with lichens fungal symbionts (LFS) and is able to degrade polysaccharides found in the cell wall of LFS. These findings indicate that a combination of gene family expansion and horizontal gene transfer provided the basis for lichenization to evolve in chlorophyte algae.

Habitat and tree species identity shape aboveground and belowground fungal communities in central European forests.

Introduction: Trees interact with fungi in mutualistic, saprotrophic, and pathogenic relationships. With their extensive aboveground and belowground structures, trees provide diverse habitats for fungi. Thus, tree species identity is an important driver of fungal community composition in forests. Methods: Here we investigate how forest habitat (bark surface vs. soil) and tree species identity (deciduous vs. coniferous) affect fungal communities in two Central European forests. We assess differences and interactions between fungal communities associated with bark surfaces and soil, in forest plots dominated either by Fagus sylvatica, Picea abies, or Pinus sylvestris in two study regions in southwestern and northeastern Germany. Results: ITS metabarcoding yielded 3,357 fungal amplicon sequence variants (ASVs) in the northern and 6,088 in the southern region. Overall, soil communities were 4.7 times more diverse than bark communities. Habitat type explained 48-69% of the variation in alpha diversity, while tree species identity explained >1-3%. NMDS ordinations showed that habitat type and host tree species structured the fungal communities. Overall, few fungal taxa were shared between habitats, or between tree species, but the shared taxa were highly abundant. Network analyses, based on co-occurrence patterns, indicate that aboveground and belowground communities form distinct subnetworks. Discussion: Our study suggests that habitat (bark versus soil) and tree species identity are important factors structuring fungal communities in temperate European forests. The aboveground (bark-associated) fungal community is currently poorly known, including a high proportion of reads assigned to "unknown Ascomycota" or "unknown Dothideomycetes." The role of bark as a habitat and reservoir of unique fungal diversity in forests has been underestimated.

Gene abundance linked to climate zone: Parallel evolution of gene content along elevation gradients in lichenized fungi.

Introduction: Intraspecific genomic variability affects a species' adaptive potential toward climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many of these species have remarkable environmental tolerances, and often form populations across different climate zones. Here, we combine comparative and population genomics to assess the presence and absence of genes in high and low elevation genomes of two lichenized fungi of the genus Umbilicaria. Methods: The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked, which genes increase or decrease in abundance within populations along elevation. Results: In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinase superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents and found that the number of this protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations. Discussion: We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.

Physiological Plasticity as a Strategy to Cope with Harsh Climatic Conditions: Ecophysiological Meta-Analysis of the Cosmopolitan Moss Ceratodon purpureus in the Southern Hemisphere.

Determining the physiological tolerance ranges of species is necessary to comprehend the limits of their responsiveness under strong abiotic pressures. For this purpose, the cosmopolitan moss Ceratodon purpureus (Hedw.) Brid. is a good model due to its wide geographical distribution throughout different biomes and habitats. In order to disentangle how this species copes with stresses such as extreme temperatures and high radiation, we designed a meta-analysis by including the main photosynthetic traits obtained by gas exchange measurements in three contrasting habitats from the Southern Hemisphere. Our findings highlight that traits such as respiration homeostasis, modulation of the photosynthetic efficiency, adjustment of the optimal temperature, and switching between shade and sun-adapted forms, which are crucial in determining the responsiveness of this species. In fact, these ecophysiological traits are in concordance with the climatic particularities of each habitat. Furthermore, the photosynthetic trends found in our study point out how different Livingston Island (Maritime Antarctica) and Granite Harbour (Continental Antarctica) are for plant life, while the population from the Succulent Karoo Desert (South Africa) shares traits with both Antarctic regions. Altogether, the study highlights the high resilience of C. purpureus under abrupt climate changes and opens new perspectives about the wide spectrum of physiological responses of cryptogams to cope with climate change scenarios.