Ecological specialization in Trebouxia (Trebouxiophyceae) photobionts of Ramalina menziesii (Ramalinaceae) across six range-covering ecoregions of western North America

Metatranscriptomics reveals diversity of symbiotic interaction and mechanisms of carbon exchange in the marine cyanolichen Lichina pygmaea

Lichens are exemplar symbioses based upon carbon exchange between photobionts and their mycobiont hosts. Historically considered a two-way relationship, some lichen symbioses have been shown to contain multiple photobiont partners; however, the way in which these photobiont communities react to environmental change is poorly understood. Lichina pygmaea is a marine cyanolichen that inhabits rocky seashores where it is submerged in seawater during every tidal cycle. Recent work has indicated that L. pygmaea has a complex photobiont community including the cyanobionts Rivularia and Pleurocapsa. We performed rRNA-based metabarcoding and mRNA metatranscriptomics of the L. pygmaea holobiont at high and low tide to investigate community response to immersion in seawater. Carbon exchange in L. pygmaea is a dynamic process, influenced by both tidal cycle and the biology of the individual symbiotic components. The mycobiont and two cyanobiont partners exhibit distinct transcriptional responses to seawater hydration. Sugar-based compatible solutes produced by Rivularia and Pleurocapsa in response to seawater are a potential source of carbon to the mycobiont. We propose that extracellular processing of photobiont-derived polysaccharides is a fundamental step in carbon acquisition by L. pygmaea and is analogous to uptake of plant-derived carbon in ectomycorrhizal symbioses.

Circadian clock- and temperature-associated genes contribute to overall genomic differentiation along elevation in lichenized fungi

Circadian regulation is linked to local environmental adaptation, and many species with broad climatic niches display variation in circadian genes. Here, we hypothesize that lichenizing fungi occupying different climate zones tune their metabolism to local environmental conditions with the help of their circadian systems. We study two species of the genus Umbilicaria occupying similar climatic niches (Mediterranean and the cold temperate) in different continents. Using homology to Neurospora crassa genes, we identify gene sets associated with circadian rhythms (11 core, 39 peripheral genes) as well as temperature response (37 genes). Nucleotide diversity of these genes is significantly correlated with mean annual temperature, minimum temperature of the coldest month and mean temperature of the coldest quarter. Furthermore, we identify altitudinal clines in allele frequencies in several non-synonymous substitutions in core clock components, for example, white collar-like, frh-like and various ccg-like genes. A dN/dS approach revealed a few significant peripheral clock- and temperature-associated genes (e.g. ras-1-like, gna-1-like) that may play a role in fine-tuning the circadian clock and temperature-response machinery. An analysis of allele frequency changes demonstrated the strongest evidence for differentiation above the genomic background in the clock-associated genes in U. pustulata. These results highlight the likely relevance of the circadian clock in environmental adaptation, particularly frost tolerance, of lichens. Whether or not the fungal clock modulates the symbiotic interaction within the lichen consortium remains to be investigated. We corroborate the finding of genetic variation in clock components along altitude-not only latitude-as has been reported in other species.

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.

Regionally Varying Habitat Relationships in Lichens: The Concept and Evidence with an Emphasis on North-Temperate Ecosystems

Habitat ecology of lichens (lichen-forming fungi) involves diverse adaptations to stressful environments where lichens use specific habitat conditions. Field observations confirm that such habitat ‘preferences’ can vary significantly across species’ distribution ranges, sometimes revealing abrupt changes over short distances. We critically review and generalize such empirical evidence as broad ecological patterns, link these with the likely physiological mechanisms and evolutionary processes involved, and outline the implications for lichen conservation. Non-replicated correlative studies remain only suggestive because the data are frequently compromised by sampling bias and pervasive random errors; further noise is related to unrecognized cryptic species. Replicated evidence exists for three macroecological patterns: (a) regional limiting factors excluding a species from a part of its microhabitat range in suboptimal areas; (b) microhabitat shifts to buffer regionally adverse macroclimates; (c) substrate suitability changed by the chemical environment, notably air pollution. All these appear to be primarily buffering physiological challenges of the adverse conditions at the macrohabitat scale or, in favorable environments, coping with competition or predation. The roles of plasticity, adaptation, dispersal, and population-level stochasticity remain to be studied. Although lichens can inhabit various novel microhabitats, there is no evidence for a related adaptive change. A precautionary approach to lichen conservation is to maintain long-term structural heterogeneity in lichen habitats, and consider lichen ecotypes as potential evolutionarily significant units and a bet-hedging strategy for addressing the climate change-related challenges to biodiversity.

Symbiont-specific responses to environmental cues in a threesome lichen symbiosis

Photosymbiodemes are a special case of lichen symbiosis where one lichenized fungus engages in symbiosis with two different photosynthetic partners, a cyanobacterium and a green alga, to develop two distinctly looking photomorphs. We compared gene expression of thallus sectors of the photosymbiodeme-forming lichen Peltigera britannica containing cyanobacterial photobionts with thallus sectors with both green algal and cyanobacterial photobionts and investigated differential gene expression at different temperatures representing mild and putatively stressful conditions. First, we quantified photobiont-mediated differences in fungal gene expression. Second, because of known ecological differences between photomorphs, we investigated symbiont-specific responses in gene expression to temperature increases. Photobiont-mediated differences in fungal gene expression could be identified, with upregulation of distinct biological processes in the different morphs, showing that interaction with specific symbiosis partners profoundly impacts fungal gene expression. Furthermore, high temperatures expectedly led to an upregulation of genes involved in heat shock responses in all organisms in whole transcriptome data and to an increased expression of genes involved in photosynthesis in both photobiont types at 15 and 25°C. The fungus and the cyanobacteria exhibited thermal stress responses already at 15°C, the green algae mainly at 25°C, demonstrating symbiont-specific responses to environmental cues and symbiont-specific ecological optima.

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.

Photobiont Diversity in Lichen Symbioses From Extreme Environments

Fungal–algal relationships—both across evolutionary and ecological scales—are finely modulated by the presence of the symbionts in the environments and by the degree of selectivity and specificity that either symbiont develop reciprocally. In lichens, the green algal genus Trebouxia Puymaly is one of the most frequently recovered chlorobionts. Trebouxia species-level lineages have been recognized on the basis of their morphological and phylogenetic diversity, while their ecological preferences and distribution are still only partially unknown. We selected two cosmopolitan species complexes of lichen-forming fungi as reference models, i.e., Rhizoplaca melanophthalma and Tephromela atra, to investigate the diversity of their associated Trebouxia spp. in montane habitats across their distributional range worldwide. The greatest diversity of Trebouxia species-level lineages was recovered in the altitudinal range 1,000–2,500 m a.s.l. A total of 10 distinct Trebouxia species-level lineages were found to associate with either mycobiont, for which new photobionts are reported. One previously unrecognized Trebouxia species-level lineage was identified and is here provisionally named Trebouxia “A52.” Analyses of cell morphology and ultrastructure were performed on axenically isolated strains to fully characterize the new Trebouxia “A52” and three other previously recognized lineages, i.e., Trebouxia “A02,” T. vagua “A04,” and T. vagua “A10,” which were successfully isolated in culture during this study. The species-level diversity of Trebouxia associating with the two lichen-forming fungi in extreme habitats helps elucidate the evolutionary pathways that this lichen photobiont genus traversed to occupy varied climatic and vegetative regimes.

Phylogeny and Ecology of Trebouxia Photobionts From Bolivian Lichens

In the past few years, new phylogenetic lineages in Trebouxia were detected as a result of molecular approaches. These studies included symbiont selectivity in lichen communities, transects along altitudinal gradients at local and global scales and the photobiont diversity in local populations of lichen-forming fungal species. In most of these studies, phylogenetic and haplotype analyses based on the internal transcribed spacer (ITS) locus have continuously allowed the recognition of new monophyletic lineages, which suggests that still numerous undiscovered Trebouxia lineages can be hidden in lichens from unexplored areas, especially in the tropics. Here, we estimated the biodiversity of photobionts in Bolivian Andean vegetation and assessed their specificity. About 403 lichen samples representing 42 genera, e.g., Haematomma, Heterodermia, Hypotrachyna, Lecanora, Lepra, Leucodermia, Parmotrema, Pertusaria, Polyblastidium, and Usnea, containing Trebouxia photobionts, were analyzed. ITS ribosomal DNA (rDNA) and rbcL markers were used. We obtained Trebouxia sequences from Bolivian samples belonging to already described clades A, C, I, and S. Thirty-nine Trebouxia lineages were distinguished within these clades, while 16 were new. To reveal the structure of the community of Bolivian photobionts and their relationships with mycobionts, the comparative effects of climate, altitude, geographical distances, substrate, and habitat type, as well as functional traits of lichens such as growth forms, propagation mode and secondary metabolites, were analyzed. Furthermore, new Bolivian records were included in analysis on a global scale. In our study, the mycobiont genus or even species are the most important factors correlated with photobiont identity. Moreover, we revealed that the community of Bolivian photobionts is shaped by altitude.

Transposable Elements in the Genome of the Lichen-Forming Fungus Umbilicaria pustulata and Their Distribution in Different Climate Zones along Elevation

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host's ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.

Turnover of Lecanoroid Mycobionts and Their Trebouxia Photobionts Along an Elevation Gradient in Bolivia Highlights the Role of Environment in Structuring the Lichen Symbiosis

Shifts in climate along elevation gradients structure mycobiont-photobiont associations in lichens. We obtained mycobiont (lecanoroid Lecanoraceae) and photobiont (Trebouxia alga) DNA sequences from 89 lichen thalli collected in Bolivia from a ca. 4,700 m elevation gradient encompassing diverse natural communities and environmental conditions. The molecular dataset included six mycobiont loci (ITS, nrLSU, mtSSU, RPB1, RPB2, and MCM7) and two photobiont loci (ITS, rbcL); we designed new primers to amplify Lecanoraceae RPB1 and RPB2 with a nested PCR approach. Mycobionts belonged to Lecanora s.lat., Bryonora, Myriolecis, Protoparmeliopsis, the "Lecanora" polytropa group, and the "L." saligna group. All of these clades except for Lecanora s.lat. occurred only at high elevation. No single species of Lecanoraceae was present along the entire elevation gradient, and individual clades were restricted to a subset of the gradient. Most Lecanoraceae samples represent species which have not previously been sequenced. Trebouxia clade C, which has not previously been recorded in association with species of Lecanoraceae, predominates at low- to mid-elevation sites. Photobionts from Trebouxia clade I occur at the upper extent of mid-elevation forest and at some open, high-elevation sites, while Trebouxia clades A and S dominate open habitats at high elevation. We did not find Trebouxia clade D. Several putative new species were found in Trebouxia clades A, C, and I. These included one putative species in clade A associated with Myriolecis species growing on limestone at high elevation and a novel lineage sister to the rest of clade C associated with Lecanora on bark in low-elevation grassland. Three different kinds of photobiont switching were observed, with certain mycobiont species associating with Trebouxia from different major clades, species within a major clade, or haplotypes within a species. Lecanoraceae mycobionts and Trebouxia photobionts exhibit species turnover along the elevation gradient, but with each partner having a different elevation threshold at which the community shifts completely. A phylogenetically defined sampling of a single diverse family of lichen-forming fungi may be sufficient to document regional patterns of Trebouxia diversity and distribution.

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.

Macroecological diversification and convergence in a clade of keystone symbionts

Lichens are classic models of symbiosis, and one of the most frequent nutritional modes among fungi. The ecologically and geographically widespread lichen-forming algal (LFA) genus Trebouxia is one of the best-studied groups of LFA and associates with over 7000 fungal species. Despite its importance, little is known about its diversification. We synthesized twenty years of publicly available data by characterizing the ecological preferences of this group and testing for time-variant shifts in climatic regimes over a distribution of trees. We found evidence for limited shifts among regimes, but that disparate lineages convergently evolved similar ecological tolerances. Early Trebouxia lineages were largely forest specialists or habitat generalists that occupied a regime whose extant members occur in moderate climates. Trebouxia then convergently diversified in non-forested habitats and expanded into regimes whose modern representatives occupy wet-warm and cool-dry climates. We rejected models in which climatic diversification slowed through time, suggesting climatic diversification is inconsistent with that expected under an adaptive radiation. In addition, we found that climatic and vegetative regime shifts broadly coincided with the evolution of biomes and associated or similar taxa. Together, our work illustrates how this keystone symbiont from an iconic symbiosis evolved to occupy diverse habitats across the globe.

New lineages of photobionts in Bolivian lichens expand our knowledge on habitat preferences and distribution of Asterochloris algae

We studied the biodiversity of Asterochloris photobionts found in Bolivian lichens to better understand their global spatial distribution and adaptation strategies in the context of a worldwide phylogeny of the genus. Based on nuclear ITS rDNA, the chloroplast rbcL gene and the actin type I gene we reconstructed a phylogenetic tree that recovered nine new Asterochloris lineages, while 32 Bolivian photobiont samples were assigned to 12 previously recognized Asterochloris lineages. We also show that some previously discovered Asterochloris photobiont species and lineages may occur in a broader spectrum of climatic conditions, and mycobiont species and photobionts may show different preferences along an altitude gradient. To reveal general patterns of of mycobiont specificity towards the photobiont in Asterochloris, we tested the influence of climate, altitude, geographical distance and effects of symbiotic partner (mycobiont) at the species level of three genera of lichen forming fungi: Stereocaulon, Cladonia and Lepraria. Further, we compared the specificity of mycobionts towards Asterochloris photobionts in cosmopolitan, Neotropical, and Pantropical lichen forming fungi. Interestingly, cosmopolitan species showed the lowest specificity to their photobionts, but also the lowest haplotype diversity. Neotropical and Paleotropical mycobionts, however, were more specific.

Contrasting co-occurrence patterns of photobiont and cystobasidiomycete yeast associated with common epiphytic lichen species

The popular dual definition of lichen symbiosis is under question with recent findings of additional microbial partners living within the lichen body. Here we compare the distribution and co-occurrence patterns of lichen photobiont and recently described secondary fungus (Cyphobasidiales yeast) to evaluate their dependency on lichen host fungus (mycobiont). We sequenced the nuclear internal transcribed spacer (ITS) strands for mycobiont, photobiont, and yeast from six widespread northern hemisphere epiphytic lichen species collected from 25 sites in Switzerland and Estonia. Interaction network analyses and multivariate analyses were conducted on operational taxonomic units based on ITS sequence data. Our study demonstrates the frequent presence of cystobasidiomycete yeasts in studied lichens and shows that they are much less mycobiont-specific than the photobionts. Individuals of different lichen species growing on the same tree trunk consistently hosted the same or closely related mycobiont-specific Trebouxia lineage over geographic distances while the cystobasidiomycete yeasts were unevenly distributed over the study area - contrasting communities were found between Estonia and Switzerland. These results contradict previous findings of high mycobiont species specificity of Cyphobasidiales yeast at large geographic scales. Our results suggest that the yeast might not be as intimately associated with the symbiosis as is the photobiont.

Cophylogenetic patterns in algal symbionts correlate with repeated symbiont switches during diversification and geographic expansion of lichen-forming fungi in the genus Sticta (Ascomycota, Peltigeraceae)

Species in the fungal genus Sticta form symbiotic associations primarily with either green algae or cyanobacteria, but tripartite associations or photosymbiodemes involving both types of photobionts occur in some species. Sticta is known to associate with green algae in the genus Symbiochloris. However, previous studies have shown that algae from other genera, such as Heveochlorella, may also be suitable partners for Sticta. We examined the diversity of green algal partners in the genus Sticta and assessed the patterns of association between the host fungus and its algal symbiont. We used multi-locus sequence data from multiple individuals collected in Australia, Cuba, Madagascar, Mauritius, New Zealand, Reunion and South America to infer phylogenies for fungal and algal partners and performed tests of congruence to assess coevolution between the partners. In addition, event-based methods were implemented to examine which cophylogenetic processes have led to the observed association patterns in Sticta and its green algal symbionts. Our results show that in addition to Symbiochloris, Sticta associates with green algae from the genera Chloroidium, Coccomyxa, Elliptochloris and Heveochlorella, the latter being the most common algal symbiont associated with Sticta in this study. Geography plays a strong role in shaping fungal-algal association patterns in Sticta as mycobionts associate with different algal lineages in different geographic locations. While fungal and algal phylogenies were mostly congruent, event-based methods did not find any evidence for cospeciation between the partners. Instead, the association patterns observed in Sticta and associated algae, were largely explained by other cophylogenetic events such as host-switches, losses of symbiont and failure of the symbiont to diverge with its host. Our results also show that tripartite associations with green algae evolved multiple times in Sticta.

Formally described species woefully underrepresent phylogenetic diversity in the common lichen photobiont genus Trebouxia (Trebouxiophyceae, Chlorophyta): An impetus for developing an integrated taxonomy

Lichens provide valuable systems for studying symbiotic interactions. In lichens, these interactions are frequently described in terms of availability, selectivity and specificity of the mycobionts and photobionts towards one another. The lichen-forming, green algal genus Trebouxia Puymaly is among the most widespread photobiont, associating with a broad range of lichen-forming fungi. To date, 29 species have been described, but studies consistently indicate that the vast majority of species-level lineages still lack formal description, and new, previously unrecognized lineages are frequently reported. To reappraise the diversity and the evolutionary relationships of species-level lineages in Trebouxia, we assembled DNA sequence data from over 1600 specimens, compiled from a range of sequences from previously published studies, axenic algal cultures, and lichens collected from poorly sampled regions. From these samples, we selected representatives of the currently known genetic diversity in the lichenized Trebouxia and inferred a phylogeny from multi-locus sequence data (ITS, rbcL, cox2). We demonstrate that the current formally described species woefully underrepresent overall species-level diversity in this important lichen-forming algal genus. We anticipate that an integrative taxonomic approach, incorporating morphological and physiological data from axenic cultures with genetic data, will be required to establish a robust, comprehensive taxonomy for Trebouxia. The data presented here provide an important impetus and reference dataset for more reliably characterizing diversity in lichenized algae and in using lichens to investigate the evolution of symbioses and holobionts.

Expanding the mutualistic niche: parallel symbiont turnover along climatic gradients

Keystone mutualisms, such as corals, lichens or mycorrhizae, sustain fundamental ecosystem functions. Range dynamics of these symbioses are, however, inherently difficult to predict because host species may switch between different symbiont partners in different environments, thereby altering the range of the mutualism as a functional unit. Biogeographic models of mutualisms thus have to consider both the ecological amplitudes of various symbiont partners and the abiotic conditions that trigger symbiont replacement. To address this challenge, we here investigate 'symbiont turnover zones'--defined as demarcated regions where symbiont replacement is most likely to occur, as indicated by overlapping abundances of symbiont ecotypes. Mapping the distribution of algal symbionts from two species of lichen-forming fungi along four independent altitudinal gradients, we detected an abrupt and consistent β-diversity turnover suggesting parallel niche partitioning. Modelling contrasting environmental response functions obtained from latitudinal distributions of algal ecotypes consistently predicted a confined altitudinal turnover zone. In all gradients this symbiont turnover zone is characterized by approximately 12°C average annual temperature and approximately 5°C mean temperature of the coldest quarter, marking the transition from Mediterranean to cool temperate bioregions. Integrating the conditions of symbiont turnover into biogeographic models of mutualisms is an important step towards a comprehensive understanding of biodiversity dynamics under ongoing environmental change.

Evolutionary conservation of within-family biodiversity patterns

The tendency for species to retain their ancestral biological properties has been widely demonstrated, but the effect of phylogenetic constraints when progressing from species to ensemble-level properties requires further assessment. Here we test whether community-level patterns (environmental shifts in local species richness and turnover) are phylogenetically conserved, assessing whether their similarity across different families of lichens, insects, and birds is dictated by the relatedness of these families. We show a significant phylogenetic signal in the shape of the species richness-elevation curve and the decay of community similarity with elevation: closely related families share community patterns within the three major taxa. Phylogenetic influences are partly explained by similarities among families in conserved traits defining body plan and interactions, implying a scaling of phylogenetic effects from the organismal to the community level. Consequently, the phylogenetic signal in community-level patterns informs about how the historical legacy of a taxon and shared responses among related taxa to similar environments contribute to community assembly and diversity patterns.

A case study on the re-establishment of the cyanolichen symbiosis: where do the compatible photobionts come from?

BACKGROUND AND AIMS

In order to re-establish lichen symbiosis, fungal spores must first germinate and then associate with a compatible photobiont. To detect possible establishment limitations in a sexually reproducing cyanolichen species, we studied ascospore germination, photobiont growth and photobiont association patterns in Pectenia plumbea.

METHODS

Germination tests were made with ascospores from 500 apothecia under different treatments, and photobiont growth was analysed in 192 isolates obtained from 24 thalli. We determined the genotype identity [tRNALeu (UAA) intron] of the Nostoc cyanobionts from 30 P. plumbea thalli from one population. We also sequenced cyanobionts of 41 specimens of other cyanolichen species and 58 Nostoc free-living colonies cultured from the bark substrate.

KEY RESULTS

Not a single fungal ascospore germinated and none of the photobiont isolates produced motile hormogonia. Genetic analyses revealed that P. plumbea shares Nostoc genotypes with two other cyanolichen species of the same habitat, but these photobionts were hardly present in the bark substrate.

CONCLUSIONS

Due to the inability of both symbionts to thrive independently, the establishment of P. plumbea seems to depend on Dendriscocaulon umhausense, the only cyanolichen species in the same habitat that reproduces asexually and acts as a source of appropriate cyanobionts. This provides support to the hypothesis about facilitation among lichens.

Comparing Sanger sequencing and high-throughput metabarcoding for inferring photobiont diversity in lichens

The implementation of HTS (high-throughput sequencing) approaches is rapidly changing our understanding of the lichen symbiosis, by uncovering high bacterial and fungal diversity, which is often host-specific. Recently, HTS methods revealed the presence of multiple photobionts inside a single thallus in several lichen species. This differs from Sanger technology, which typically yields a single, unambiguous algal sequence per individual. Here we compared HTS and Sanger methods for estimating the diversity of green algal symbionts within lichen thalli using 240 lichen individuals belonging to two species of lichen-forming fungi. According to HTS data, Sanger technology consistently yielded the most abundant photobiont sequence in the sample. However, if the second most abundant photobiont exceeded 30% of the total HTS reads in a sample, Sanger sequencing generally failed. Our results suggest that most lichen individuals in the two analyzed species, Lasallia hispanica and L. pustulata, indeed contain a single, predominant green algal photobiont. We conclude that Sanger sequencing is a valid approach to detect the dominant photobionts in lichen individuals and populations. We discuss which research areas in lichen ecology and evolution will continue to benefit from Sanger sequencing, and which areas will profit from HTS approaches to assessing symbiont diversity.

Oligocene origin and drivers of diversification in the genus Sticta (Lobariaceae, Ascomycota)

A major challenge to evolutionary biologists is to understand how biodiversity is distributed through space and time and across the tree of life. Diversification of organisms is influenced by many factors that act at different times and geographic locations but it is still not clear which have a significant impact and how drivers interact. To study diversification, we chose the lichen genus Sticta, by sampling through most of the global range and producing a time tree. We estimate that Sticta originated about 30 million years ago, but biogoegraphic analysis was unclear in estimating the origin of the genus. Furthermore, we investigated the effect of dispersal ability finding that Sticta has a high dispersal rate, as collections from Hawaii showed that divergent lineages colonized the islands at least four times. Symbiont interactions were investigated using BiSSE to understand if green-algal or cyanobacterial symbiont interactions influenced diversification, only to find that the positive results were driven almost completely by Type I error. On the other hand, another BiSSE analysis found that an association with Andean tectonic activity increases the speciation rate of species.

Sharing of photobionts in sympatric populations of Thamnolia and Cetraria lichens: evidence from high-throughput sequencing

In this study, we explored the diversity of green algal symbionts (photobionts) in sympatric populations of the cosmopolitan lichen-forming fungi Thamnolia and Cetraria. We sequenced with both Sanger and Ion Torrent High-Throughput Sequencing technologies the photobiont ITS-region of 30 lichen thalli from two islands: Iceland and Öland. While Sanger recovered just one photobiont genotype from each thallus, the Ion Torrent data recovered 10-18 OTUs for each pool of 5 lichen thalli, suggesting that individual lichens can contain heterogeneous photobiont populations. Both methods showed evidence for photobiont sharing between Thamnolia and Cetraria on Iceland. In contrast, our data suggest that on Öland the two mycobionts associate with distinct photobiont communities, with few shared OTUs revealed by Ion Torrent sequencing. Furthermore, by comparing our sequences with public data, we identified closely related photobionts from geographically distant localities. Taken together, we suggest that the photobiont composition in Thamnolia and Cetraria results from both photobiont-mycobiont codispersal and local acquisition during mycobiont establishment and/or lichen growth. We hypothesize that this is a successful strategy for lichens to be flexible in the use of the most adapted photobiont for the environment.

Clair LL, Williams L, Wedin M, Lumbsch HT. Multiple, Distinct Intercontinental Lineages but Isolation of Australian Populations in a Cosmopolitan Lichen-Forming Fungal Taxon, Psora decipiens (Psoraceae, Ascomycota)

Multiple drivers shape the spatial distribution of species, including dispersal capacity, niche incumbency, climate variability, orographic barriers, and plate tectonics. However, biogeographic patterns of fungi commonly do not fit conventional expectations based on studies of animals and plants. Fungi, in general, are known to occur across exceedingly broad, intercontinental distributions, including some important components of biological soil crust communities (BSCs). However, molecular data often reveal unexpected biogeographic patterns in lichenized fungal species that are assumed to have cosmopolitan distributions. The lichen-forming fungal species Psora decipiens is found on all continents, except Antarctica and occurs in BSCs across diverse habitats, ranging from hot, arid deserts to alpine habitats. In order to better understand factors that shape population structure in cosmopolitan lichen-forming fungal species, we investigated biogeographic patterns in the cosmopolitan taxon P. decipiens, along with the closely related taxa P. crenata and P. saviczii. We generated a multi-locus sequence dataset based on a worldwide sampling of these taxa in order to reconstruct evolutionary relationships and explore phylogeographic patterns. Both P. crenata and P. decipiens were not recovered as monophyletic; and P. saviczii specimens were recovered as a monophyletic clade closely related to a number of lineages comprised of specimens representing P. decipiens. Striking phylogeographic patterns were observed for P. crenata, with populations from distinct geographic regions belonging to well-separated, monophyletic lineages. South African populations of P. crenata were further divided into well-supported sub-clades. While well-supported phylogenetic substructure was also observed for the nominal taxon P. decipiens, nearly all lineages were comprised of specimens collected from intercontinental populations. However, all Australian specimens representing P. decipiens were recovered within a single well-supported monophyletic clade consisting solely of Australian samples. Our study supports up to 10 candidate species-level lineages in P. decipiens, based on genealogical concordance and coalescent-based species delimitation analyses. Our results support the general pattern of the biogeographic isolation of lichen-forming fungal populations in Australia, even in cases where closely related congeners have documented intercontinental distributions. Our study has important implications for understanding factors influencing diversification and distributions of lichens associated with BSC.

Novel climates reverse carbon uptake of atmospherically dependent epiphytes: Climatic constraints on the iconic boreal forest lichen Evernia mesomorpha

PREMISE OF THE STUDY

Changing climates are expected to affect the abundance and distribution of global vegetation, especially plants and lichens with an epiphytic lifestyle and direct exposure to atmospheric variation. The study of epiphytes could improve understanding of biological responses to climatic changes, but only if the conditions that elicit physiological performance changes are clearly defined.

METHODS

We evaluated individual growth performance of the epiphytic lichen Evernia mesomorpha, an iconic boreal forest indicator species, in the first year of a decade-long experiment featuring whole-ecosystem warming and drying. Field experimental enclosures were located near the southern edge of the species' range.

KEY RESULTS

Mean annual biomass growth of Evernia significantly declined 6 percentage points for every +1°C of experimental warming after accounting for interactions with atmospheric drying. Mean annual biomass growth was 14% in ambient treatments, 2% in unheated control treatments, and -9% to -19% (decreases) in energy-added treatments ranging from +2.25 to +9.00°C above ambient temperatures. Warming-induced biomass losses among persistent individuals were suggestive evidence of an extinction debt that could precede further local mortality events.

CONCLUSIONS

Changing patterns of warming and drying would decrease or reverse Evernia growth at its southern range margins, with potential consequences for the maintenance of local and regional populations. Negative carbon balances among persisting individuals could physiologically commit these epiphytes to local extinction. Our findings illuminate the processes underlying local extinctions of epiphytes and suggest broader consequences for range shrinkage if dispersal and recruitment rates cannot keep pace.

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.

Evolutionary lessons from California plant phylogeography

Phylogeography documents the spatial distribution of genetic lineages that result from demographic processes, such as population expansion, population contraction, and gene movement, shaped by climate fluctuations and the physical landscape. Because most phylogeographic studies have used neutral markers, the role of selection may have been undervalued. In this paper, we contend that plants provide a useful evolutionary lesson about the impact of selection on spatial patterns of neutral genetic variation, when the environment affects which individuals can colonize new sites, and on adaptive genetic variation, when environmental heterogeneity creates divergence at specific loci underlying local adaptation. Specifically, we discuss five characteristics found in plants that intensify the impact of selection: sessile growth form, high reproductive output, leptokurtic dispersal, isolation by environment, and the potential to evolve longevity. Collectively, these traits exacerbate the impact of environment on movement between populations and local selection pressures-both of which influence phylogeographic structure. We illustrate how these unique traits shape these processes with case studies of the California endemic oak, Quercus lobata, and the western North American lichen, Ramalina menziesii Obviously, the lessons we learn from plant traits are not unique to plants, but they highlight the need for future animal, plant, and microbe studies to incorporate its impact. Modern tools that generate genome-wide sequence data are now allowing us to decipher how evolutionary processes affect the spatial distribution of different kinds of genes and also to better model future spatial distribution of species in response to climate change.

Gene flow and natural selection shape spatial patterns of genes in tree populations: implications for evolutionary processes and applications

A central question in evolutionary biology is how gene flow and natural selection shape geographic patterns of genotypic and phenotypic variation. My overall research program has pursued this question in tree populations through complementary lines of inquiry. First, through studies of contemporary pollen and seed movement, I have studied how limited gene movement creates fine-scale genetic structure, while long-distance gene flow promotes connectivity. My collaborators and I have provided new tools to study these processes at a landscape scale as well as statistical tests to determine whether changes in landscape conditions or dispersal vectors affect gene movement. Second, my research on spatial patterns of genetic variation has investigated the interacting impacts of geography and climate on gene flow and selection. Third, using next-generation genomic tools, I am now studying genetic variation on the landscape to find initial evidence of climate-associated local adaptation and epigenetic variation to explore its role in plant response to the climate. By integrating these separate lines of inquiry, this research provides specific insight into real-world mechanisms shaping evolution in tree populations and potential impacts of landscape transformation and climate change on these populations, with the prospective goal of contributing to their management and conservation.

Distribution patterns of haplotypes for symbionts from Umbilicaria esculenta and U. muehlenbergii reflect the importance of reproductive strategy in shaping population genetic structure

BACKGROUND

The diversity of lichen fungal components and their photosynthetic partners reflects both ecological and evolutionary factors. In present study, molecular investigations of the internal transcribed spacer of the nuclear ribosomal DNA (ITS nrDNA) region were conducted to analyze the genetic diversity of Umbilicaria esculenta and U. muehlenbergii together with their associated green algae.

RESULT

It was here demonstrated that the reproductive strategy is a principal reason for fungal selectivity to algae. U. muehlenbergii, which disperses via sexual spores, exhibits lower selectivity to its photosynthetic partners than U. esculenta, which has a vegetative reproductive strategy. The difference of genotypic diversity (both fungal and algal) between these two Umbilicaria species is low, although their nucleotide diversity can vary greatly.

CONCLUSIONS

The present study illustrates that lichen-forming fungi with sexual reproductive strategies are less selective with respect to their photobionts; and reveals that both sexual and vegetative reproduction allow lichens to generate similar amounts of diversity to adapt to the environments. The current study will be helpful for elucidating how lichens with different reproductive strategies adapt to changing environments.

Biogeography and Genetic Structure in Populations of a Widespread Lichen (Parmelina tiliacea, Parmeliaceae, Ascomycota)

The genetic diversity and population structure of the foliose lichenized fungus Parmelina tiliacea has been analyzed through its geographical range, including samples from Macaronesia (Canary Islands), the Mediterranean, and Eurosiberia. DNA sequences from the nuclear ribosomal internal transcribed spacer, the mitochondrial large subunit ribosomal RNA gene, and the translation elongation factor 1-α were used as molecular markers. The haplotypes of the three markers and the molecular variance analyses of multilocus haplotypes showed the highest diversity in the Canary Islands, while restricted haplotypes occurred at high frequencies in Mediterranean coastal samples. The multilocus haplotypes formed three unevenly distributed clusters (clusters 1-3). In the Canary Islands all the haplotypes were present in a similar proportion, while the coastal Mediterranean sites had almost exclusively haplotypes of cluster 3; cluster 2 predominated in inland Mediterranean sites; and cluster 1 was more abundant in central and northern Europe (Eurosiberian area). The distribution of clusters is partially explained by climatic factors, and its interaction with local spatial structure, but much of the variation remains unexplained. The high frequency of individuals in the Canary Islands with haplotypes shared with other areas suggests that could be a refugium of genetic diversity, and the high frequency of individuals of the Mediterranean coastal sites with restricted haplotypes indicates that gene flow to contiguous areas may be restricted. This is significant for the selection of areas for conservation purposes, as those with most genetic variation may reflect historical factors and biological properties of the species.