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

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.

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.

Myrmecia, Not Asterochloris, Is the Main Photobiont of Cladonia subturgida (Cladoniaceae, Lecanoromycetes)

This study explores the diversity of photobionts associated with the Mediterranean lichen-forming fungus Cladonia subturgida. For this purpose, we sequenced the whole ITS rDNA region by Sanger using a metabarcoding method for ITS2. A total of 41 specimens from Greece, Italy, France, Portugal, and Spain were studied. Additionally, two specimens from Spain were used to generate four cultures. Our molecular studies showed that the genus Myrmecia is the main photobiont of C. subturgida throughout its geographic distribution. This result contrasts with previous studies, which indicated that the main photobiont for most Cladonia species is Asterochloris. The identity of Myrmecia was also confirmed by ultrastructural studies of photobionts within the lichen thalli and cultures. Photobiont cells showed a parietal chloroplast lacking a pyrenoid, which characterizes the species in this genus. Phylogenetic analyses indicate hidden diversity within this genus. The results of amplicon sequencing showed the presence of multiple ASVs in 58.3% of the specimens studied.

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.

Cryptic community structure and metabolic interactions among the heritable facultative symbionts of the pea aphid

Most insects harbour influential, yet non-essential heritable microbes in their hemocoel. Communities of these symbionts exhibit low diversity. But their frequent multi-species nature raises intriguing questions on roles for symbiont-symbiont synergies in host adaptation, and on the stability of the symbiont communities, themselves. In this study, we build on knowledge of species-defined symbiont community structure across US populations of the pea aphid, Acyrthosiphon pisum. Through extensive symbiont genotyping, we show that pea aphids' microbiomes can be more precisely defined at the symbiont strain level, with strain variability shaping five out of nine previously reported co-infection trends. Field data provide a mixture of evidence for synergistic fitness effects and symbiont hitchhiking, revealing causes and consequences of these co-infection trends. To test whether within-host metabolic interactions predict common versus rare strain-defined communities, we leveraged the high relatedness of our dominant, community-defined symbiont strains vs. 12 pea aphid-derived Gammaproteobacteria with sequenced genomes. Genomic inference, using metabolic complementarity indices, revealed high potential for cooperation among one pair of symbionts-Serratia symbiotica and Rickettsiella viridis. Applying the expansion network algorithm, through additional use of pea aphid and obligate Buchnera symbiont genomes, Serratia and Rickettsiella emerged as the only symbiont community requiring both parties to expand holobiont metabolism. Through their joint expansion of the biotin biosynthesis pathway, these symbionts may span missing gaps, creating a multi-party mutualism within their nutrient-limited, phloem-feeding hosts. Recent, complementary gene inactivation, within the biotin pathways of Serratia and Rickettsiella, raises further questions on the origins of mutualisms and host-symbiont interdependencies.

Fungal-Algal Association Drives Lichens' Mutualistic Symbiosis: A Case Study with Trebouxia-Related Lichens

Biotic and abiotic factors influence the formation of fungal-algal pairings in lichen symbiosis. However, the specific determinants of these associations, particularly when distantly related fungi are involved, remain poorly understood. In this study, we investigated the impact of different drivers on the association patterns between taxonomically diverse lichenized fungi and their trebouxioid symbiotic partners. We collected 200 samples from four biomes and identified 41 species of lichenized fungi, associating them with 16 species of trebouxioid green algae, of which 62% were previously unreported. The species identity of both the fungal and algal partners had the most significant effect on the outcome of the symbiosis, compared to abiotic factors like climatic variables and geographic distance. Some obviously specific associations were observed in the temperate zone; however, the nestedness value was lower in arid regions than in cold, polar, and temperate regions according to interaction network analysis. Cophylogenetic analyses revealed congruent phylogenies between trebouxioid algae and associated fungi, indicating a tendency to reject random associations. The main evolutionary mechanisms contributing to the observed phylogenetic patterns were "loss" and "failure to diverge" of the algal partners. This study broadens our knowledge of fungal-algal symbiotic patterns in view of Trebouxia-associated fungi.

Specialization patterns in symbiotic associations: A community perspective over spatial scales

Specialization, contextualized in a resource axis of an organism niche, is a core concept in ecology. In biotic interactions, specialization can be determined by the range of interacting partners. Evolutionary and ecological factors, in combination with the surveyed scale (spatial, temporal, biological, and/or taxonomic), influence the conception of specialization. This study aimed to assess the specialization patterns and drivers in the lichen symbiosis, considering the interaction between the principal fungus (mycobiont) and the associated Nostoc (cyanobiont), from a community perspective considering different spatial scales. Thus, we determined Nostoc phylogroup richness and composition of lichen communities in 11 Nothofagus pumilio forests across a wide latitudinal gradient in Chile. To measure specialization, cyanobiont richness, Simpson's and d' indices were estimated for 37 mycobiont species in these communities. Potential drivers that might shape Nostoc composition and specialization measures along the environmental gradient were analysed. Limitations in lichen distributional ranges due to the availability of their cyanobionts were studied. Turnover patterns of cyanobionts were identified at multiple spatial scales. The results showed that environmental factors shaped the Nostoc composition of these communities, thus limiting cyanobiont availability to establish the symbiotic association. Besides, specialization changed with the spatial scale and with the metric considered. Cyanolichens were more specialized than cephalolichens when considering partner richness and Simpson's index, whereas the d' index was mostly explained by mycobiont identity. Little evidence of lichen distributional ranges due to the distribution of their cyanobionts was found. Thus, lichens with broad distributional ranges either associated with several cyanobionts or with widely distributed cyanobionts. Comparisons between local and regional scales showed a decreasing degree of specialization at larger scales due to an increase in cyanobiont richness. The results support the context dependency of specialization and how its consideration changes with the metric and the spatial scale considered. Subsequently, we suggest considering the entire community and widening the spatial scale studied as it is crucial to understand factors determining specialization.

Disentangling arbuscular mycorrhizal fungi and bacteria at the soil-root interface

Arbuscular mycorrhizal fungi (AMF) are essential components of the plant root mycobiome and are found in approximately 80% of land plants. As obligate plant symbionts, AMF harbor their own microbiota, both inside and outside the plant root system. AMF-associated bacteria (AAB) possess various functional traits, including nitrogen fixation, organic and inorganic phosphate mobilization, growth hormone production, biofilm production, enzymatic capabilities, and biocontrol against pathogen attacks, which not only contribute to the health of the arbuscular mycorrhizal symbiosis but also promote plant growth. Because of this, there is increasing interest in the diversity, functioning, and mechanisms that underlie the complex interactions between AMF, AAB, and plant hosts. This review critically examines AMF-associated bacteria, focusing on AAB diversity, the factors driving richness and community composition of these bacteria across various ecosystems, along with the physical, chemical, and biological connections that enable AMF to select and recruit beneficial bacterial symbionts on and within their structures and hyphospheres. Additionally, potential applications of these bacteria in agriculture are discussed, emphasizing the potential importance of AMF fungal highways in engineering plant rhizosphere and endophyte bacteria communities, and the importance of a functional core of AAB taxa as a promising tool to improve plant and soil productivity. Thus, AMF and their highly diverse bacterial taxa represent important tools that could be efficiently explored in sustainable agriculture, carbon sequestration, and reduction of greenhouse gas emissions related to nitrogen fertilizer applications. Nevertheless, future studies adopting integrated multidisciplinary approaches are crucial to better understand AAB functional diversity and the mechanisms that govern these tripartite relationships.

Metabarcoding of Antarctic Lichens from Areas with Different Deglaciation Times Reveals a High Diversity of Lichen-Associated Communities

Lichens have developed numerous adaptations to optimise their survival under harsh abiotic stress, colonise different substrates, and reach substantial population sizes and high coverage in ice-free Antarctic areas, benefiting from a symbiotic lifestyle. As lichen thalli represent consortia with an unknown number of participants, it is important to know about the accessory organisms and their relationships with various environmental conditions. To this end, we analysed lichen-associated communities from Himantormia lugubris, Placopsis antarctica, P. contortuplicata, and Ramalina terebrata, collected from soils with differing deglaciation times, using a metabarcoding approach. In general, many more Ascomycete taxa are associated with the investigated lichens compared to Basidiomycota. Given our sampling, a consistently higher number of lichen-associated eukaryotes are estimated to be present in areas with deglaciation times of longer than 5000 years compared to more recently deglaciated areas. Thus far, members of Dothideomycetes, Leotiomycetes, and Arthoniomycetes have been restricted to the Placopsis specimens from areas with deglaciation times longer than 5000 years. Striking differences between the associated organisms of R. terebrata and H. lugubris have also been discovered. Thus, a species-specific basidiomycete, Tremella, was revealed for R. terebrata, as was a member of Capnodiales for H. lugubris. Our study provides further understanding of the complex terricolous lichen-associated mycobiome using the metabarcoding approach. It also illustrates the necessity to extend our knowledge of complex lichen symbiosis and further improve the coverage of microbial eukaryotes in DNA barcode libraries, including more extended sampling.

Lichen ecophysiology in a changing climate

Lichens are one of the most iconic and ubiquitous symbioses known, widely valued as indicators of environmental quality and, more recently, climate change. Our understanding of lichen responses to climate has greatly expanded in recent decades, but some biases and constraints have shaped our present knowledge. In this review we focus on lichen ecophysiology as a key to predicting responses to present and future climates, highlighting recent advances and remaining challenges. Lichen ecophysiology is best understood through complementary whole-thallus and within-thallus scales. Water content and form (vapor or liquid) are central to whole-thallus perspectives, making vapor pressure differential (VPD) a particularly informative environmental driver. Responses to water content are further modulated by photobiont physiology and whole-thallus phenotype, providing clear links to a functional trait framework. However, this thallus-level perspective is incomplete without also considering within-thallus dynamics, such as changing proportions or even identities of symbionts in response to climate, nutrients, and other stressors. These changes provide pathways for acclimation, but their understanding is currently limited by large gaps in our understanding of carbon allocation and symbiont turnover in lichens. Lastly, the study of lichen physiology has mainly prioritized larger lichens at high latitudes, producing valuable insights but underrepresenting the range of lichenized lineages and ecologies. Key areas for future work include improving geographic and phylogenetic coverage, greater emphasis on VPD as a climatic factor, advances in the study of carbon allocation and symbiont turnover, and the incorporation of physiological theory and functional traits in our predictive models.

Fungal Host Affects Photosynthesis in a Lichen Holobiont

Corals and lichens are iconic examples of photosynthetic holobionts, i.e., ecological and evolutionary units resulting from the tightly integrated association of algae and prokaryotic microbiota with animal or fungal hosts, respectively. While the role of the coral host in modulating photosynthesis has been clarified to a large extent in coral holobionts, the role of the fungal host in this regard is far less understood. Here, we address this question by taking advantage of the recent discovery of highly specific fungal-algal pairings corresponding to climatically adapted ecotypes of the lichen-forming genus Umbilicaria. Specifically, we compared chlorophyll a fluorescence kinetics among lichen thalli consisting of different fungal-algal combinations. We show that photosynthetic performance in these lichens is not only driven by algal genotype, but also by fungal host species identity and intra-host genotype. These findings shed new light on the closely intertwined physiological processes of fungal and algal partners in the lichen symbiosis. Indeed, the specific combinations of fungal and algal genotypes within a lichen individual-and the resulting combined functional phenotype-can be regarded as a response to the environment. Our findings suggest that characterizing the genetic composition of both eukaryotic partners is an important complimentary step to understand and predict the lichen holobiont's responses to environmental change.

Local adaptation in fungi

In this review, I explore the pervasive but underappreciated role of local adaptation in fungi. It has been difficult historically to study local adaptation in fungi because of the limited understanding of fungal species and their traits, but new hope has been offered with technological advances in sequencing. The filamentous nature of fungi invalidates some assumptions made in evolution because of their ability to exist as multinucleate entities with genetically different nuclei sharing the same cytoplasm. Many insights on local adaptation have come from studying fungi, and much of the empirical evidence gathered about local adaptation in the context of host-pathogen interactions comes from studying fungal virulence genes, drug resistance, and environmental adaptation. Together, these insights paint a picture of the variety of processes involved in fungal local adaptation and their connections to the unusual cell biology of Fungi (multinucleate, filamentous habit), but there is much that remains unknown, with major gaps in our knowledge of fungal species, their phenotypes, and the ways by which they adapt to local conditions.

Identification and expression of functionally conserved circadian clock genes in lichen-forming fungi

Lichen-forming fungi establish stable symbioses with green algae or cyanobacteria. Many species have broad distributions, both in geographic and ecological space, making them ideal subjects to study organism-environment interactions. However, little is known about the specific mechanisms that contribute to environmental adaptation in lichen-forming fungi. The circadian clock provides a well-described mechanism that contributes to regional adaptation across a variety of species, including fungi. Here, we identify the putative circadian clock components in phylogenetically divergent lichen-forming fungi. The core circadian genes (frq, wc-1, wc-2, frh) are present across the Fungi, including 31 lichen-forming species, and their evolutionary trajectories mirror overall fungal evolution. Comparative analyses of the clock genes indicate conserved domain architecture among lichen- and non-lichen-forming taxa. We used RT-qPCR to examine the core circadian loop of two unrelated lichen-forming fungi, Umbilicaria pustulata (Lecanoromycetes) and Dermatocarpon miniatum (Eurotiomycetes), to determine that the putative frq gene is activated in a light-dependent manner similar to the model fungus Neurospora crassa. Together, these results demonstrate that lichen-forming fungi retain functional light-responsive mechanisms, including a functioning circadian clock. Our findings provide a stepping stone into investigating the circadian clock in the lichen symbiosis, e.g. its role in adaptation, and in synchronizing the symbiotic interaction.

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.

The Role of Photobionts as Drivers of Diversification in an Island Radiation of Lichen-Forming Fungi

Speciation in oceanic islands has attracted the interest of scientists since the 19th century. One of the most striking evolutionary phenomena that can be studied in islands is adaptive radiation, that is, when a lineage gives rise to different species by means of ecological speciation. Some of the best-known examples of adaptive radiation are charismatic organisms like the Darwin finches of the Galapagos and the cichlid fishes of the great African lakes. In these and many other examples, a segregation of the trophic niche has been shown to be an important diversification driver. Radiations are known in other groups of organisms, such as lichen-forming fungi. However, very few studies have investigated their adaptive nature, and none have focused on the trophic niche. In this study, we explore the role of the trophic niche in a putative radiation of endemic species from the Macaronesian Region, the Ramalina decipiens group. The photobiont diversity was studied by Illumina MiSeq sequencing of the ITS2 region of 197 specimens spanning the phylogenetic breadth and geographic range of the group. A total of 66 amplicon sequence variants belonging to the four main clades of the algal genus Trebouxia were found. Approximately half of the examined thalli showed algal coexistence, but in most of them, a single main photobiont amounted to more than 90% of the reads. However, there were no significant differences in photobiont identity and in the abundance of ITS2 reads across the species of the group. We conclude that a segregation of the trophic niche has not occurred in the R. decipiens radiation.

Choosing the Right Life Partner: Ecological Drivers of Lichen Symbiosis

Lichens are an iconic example of symbiotic systems whose ecology is shaped by the requirements of the symbionts. Previous studies suggest that fungal (mycobionts) as well as photosynthesizing (phycobionts or cyanobionts) partners have a specific range of acceptable symbionts that can be chosen according to specific environmental conditions. This study aimed to investigate the effects of climatic conditions and mycobiont identity on phycobiont distribution within the lichen genera Stereocaulon, Cladonia, and Lepraria. The study area comprised the Canary Islands, Madeira, Sicily, and the Aeolian Islands, spanning a wide range of climatic conditions. These islands are known for their unique and diverse fauna and flora; however, lichen phycobionts have remained unstudied in most of these areas. In total, we genetically analyzed 339 lichen samples. The phycobiont pool differed significantly from that outside the studied area. Asterochloris mediterranea was identified as the most abundant phycobiont. However, its distribution was limited by climatic constraints. Other species of Asterochloris and representatives of the genera Chloroidium, Vulcanochloris, and Myrmecia were also recovered as phycobionts. The selection of symbiotic partners from the local phycobiont pool was driven by mycobiont specificity (i.e., the taxonomic range of acceptable partners) and the environmental conditions, mainly temperature. Interestingly, the dominant fungal species responded differently in their selection of algal symbionts along the environmental gradients. Cladonia rangiformis associated with its phycobiont A. mediterranea in a broader range of temperatures than Stereocaulon azoreum, which favors other Asterochloris species along most of the temperature gradient. Stereocaulon vesuvianum associated with Chloroidium spp., which also differed in their temperature optima. Finally, we described Stereocaulon canariense as a new endemic species ecologically distinct from the other Stereocaulon species on the Canary Islands.

Biodiversity Patterns and Ecological Preferences of the Photobionts Associated With the Lichen-Forming Genus Parmelia

The worldwide, ecologically relevant lichen-forming genus Parmelia currently includes 41 accepted species, of which the Parmelia sulcata group (PSULgp) and the Parmelia saxatilis group (PSAXgp) have received considerable attention over recent decades; however, phycobiont diversity is poorly known in Parmelia s. lat. Here, we studied the diversity of Trebouxia microalgae associated with 159 thalli collected from 30 locations, including nine Parmelia spp.: P. barrenoae, P. encryptata, P. ernstiae, P. mayi, P. omphalodes, P. saxatilis, P. serrana, P. submontana, and P. sulcata. The mycobionts were studied by carrying out phylogenetic analyses of the nrITS. Microalgae genetic diversity was examined by using both nrITS and LSU rDNA markers. To evaluate putative species boundaries, three DNA species delimitation analyses were performed on Trebouxia and Parmelia. All analyses clustered the mycobionts into two main groups: PSULgp and PSAXgp. Species delimitation identified 13 fungal and 15 algal species-level lineages. To identify patterns in specificity and selectivity, the diversity and abundance of the phycobionts were identified for each Parmelia species. High specificity of each Parmelia group for a given Trebouxia clade was observed; PSULgp associated only with clade I and PSAXgp with clade S. However, the degree of specificity is different within each group, since the PSAXgp mycobionts were less specific and associated with 12 Trebouxia spp., meanwhile those of PSULgp interacted only with three Trebouxia spp. Variation-partitioning analyses were conducted to detect the relative contributions of climate, geography, and symbiotic partner to phycobiont and mycobiont distribution patterns. Both analyses explained unexpectedly high portions of variability (99 and 98%) and revealed strong correlations between the fungal and algal diversity. Network analysis discriminated seven ecological clusters. Even though climatic conditions explained the largest proportion of the variation among these clusters, they seemed to show indifference relative to climatic parameters. However, the cluster formed by P. saxatilis A/P. saxatilis B/Trebouxia sp. 2/Trebouxia sp. S02/Trebouxia sp. 3A was identified to prefer cold-temperate as well as humid summer environments.

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.

Macroclimatic conditions as main drivers for symbiotic association patterns in lecideoid lichens along the Transantarctic Mountains, Ross Sea region, Antarctica

Lecideoid lichens as dominant vegetation-forming organisms in the climatically harsh areas of the southern part of continental Antarctica show clear preferences in relation to environmental conditions (i.e. macroclimate). 306 lichen samples were included in the study, collected along the Ross Sea coast (78°S–85.5°S) at six climatically different sites. The species compositions as well as the associations of their two dominant symbiotic partners (myco- and photobiont) were set in context with environmental conditions along the latitudinal gradient. Diversity values were nonlinear with respect to latitude, with the highest alpha diversity in the milder areas of the McMurdo Dry Valleys (78°S) and the most southern areas (Durham Point, 85.5°S; Garden Spur, 84.5°S), and lowest in the especially arid and cold Darwin Area (~ 79.8°S). Furthermore, the specificity of mycobiont species towards their photobionts decreased under more severe climate conditions. The generalist lichen species Lecanora fuscobrunnea and Lecidea cancriformis were present in almost all habitats, but were dominant in climatically extreme areas. Carbonea vorticosa, Lecidella greenii and Rhizoplaca macleanii were confined to milder areas. In summary, the macroclimate is considered to be the main driver of species distribution, making certain species useful as bioindicators of climate conditions and, consequently, for assessing the consequences of climate change.

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.

Opening the Gap: Rare Lichens With Rare Cyanobionts - Unexpected Cyanobiont Diversity in Cyanobacterial Lichens of the Order Lichinales

The last decades of research led to a change in understanding of lichens that are now seen as self-sustaining micro-ecosystems, harboring diverse microbial organisms in tight but yet not fully understood relationships. Among the diverse interdependencies, the relationship between the myco- and photobiont is the most crucial, determining the shape, and ecophysiological properties of the symbiotic consortium. Roughly 10% of lichens associate with cyanobacteria as their primary photobiont, termed cyanolichens. Up to now, the diversity of cyanobionts of bipartite lichens resolved by modern phylogenetic approaches is restricted to the filamentous and heterocytous genera of the order Nostocales. Unicellular photobionts were placed in the orders Chroococcales, Pleurocapsales, and Chroococcidiopsidales. However, especially the phylogeny and taxonomy of the Chroococcidiopsidales genera remained rather unclear. Here we present new data on the identity and phylogeny of photobionts from cyanolichens of the genera Gonohymenia, Lichinella, Peccania, and Peltula from a broad geographical range. A polyphasic approach was used, combining morphological and cultivation-depending characteristics (microscopy, staining techniques, life cycle observation, baeocyte motility, and nitrogen fixation test) with phylogenetic analyses of the 16S rRNA and 16S–23S ITS gene region. We found an unexpectedly high cyanobiont diversity in the cyanobacterial lichens of the order Lichinales, including two new genera and seven new species, all of which were not previously perceived as lichen symbionts. As a result, we describe the novel unicellular Chroococcidiopsidales genera Pseudocyanosarcina gen. nov. with the species Pseudocyanosarcina phycocyania sp. nov. (from Peltula clavata, Australia) and Compactococcus gen. nov. with the species Compactococcus sarcinoides sp. nov. (from Gonohymenia sp., Australia) and the new Chroococcidiopsidales species Aliterella compacta sp. nov. (from Peltula clavata, Australia), Aliterella gigantea sp. nov. (from Peltula capensis; South Africa), Sinocapsa ellipsoidea sp. nov. (from Peccania cerebriformis, Austria), as well as the two new Nostocales species Komarekiella gloeocapsoidea sp. nov. (from Gonohymenia sp., Czechia) and Komarekiella globosa sp. nov. (from Lichinella cribellifera, Canary Islands, Spain). Our study highlights the role of cyanolichens acting as a key in untangling cyanobacterial taxonomy and diversity. With this study, we hope to stimulate further research on photobionts, especially of rare cyanolichens.

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.

Microbial Communities of Cladonia Lichens and Their Biosynthetic Gene Clusters Potentially Encoding Natural Products

Lichens have been widely used in traditional medicine, especially by indigenous communities worldwide. However, their slow growth and difficulties in the isolation of lichen symbionts and associated microbes have hindered the pharmaceutical utilisation of lichen-produced compounds. Advances in high-throughput sequencing techniques now permit detailed investigations of the complex microbial communities formed by fungi, green algae, cyanobacteria, and other bacteria within the lichen thalli. Here, we used amplicon sequencing, shotgun metagenomics, and in silico metabolomics together with compound extractions to study reindeer lichens collected from Southern Finland. Our aim was to evaluate the potential of Cladonia species as sources of novel natural products. We compared the predicted biosynthetic pathways of lichen compounds from isolated genome-sequenced lichen fungi and our environmental samples. Potential biosynthetic genes could then be further used to produce secondary metabolites in more tractable hosts. Furthermore, we detected multiple compounds by metabolite analyses, which revealed connections between the identified biosynthetic gene clusters and their products. Taken together, our results contribute to metagenomic data studies from complex lichen-symbiotic communities and provide valuable new information for use in further biochemical and pharmacological studies.

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.

Thallus Growth Stage and Geographic Origin Shape Microalgal Diversity in Ramalina farinacea Lichen Holobionts

Lichen symbioses are microecosystems hosting many other living organisms besides the two major lichen symbionts (i.e., lichenized fungi [the mycobiont] and green microalgae or cyanobacteria [the photobiont]). Recent investigations evidenced that other fungi, non-photosynthetic bacteria, and microalgae co-inhabit within the lichen thalli, but their diversity and their roles are still underinvestigated. Here we present an ad hoc stratified sampling design and in-depth Illumina paired-end metabarcoding approach to explore microalgal diversity in lichen thalli of the model species Ramalina farinacea from different ecologies. Lichen thalli were surveyed according to three different sizes, and different thallus parts were considered for molecular, bioinformatics, and community diversity analyses. The results revealed that microalgal diversity strongly depends on the growth stage of the thalli, the geographic area, and the habitat type. The results also show that microalgal diversity does not vary along the thallus branches (lacinias)-that is, it does not correlate with the apical growth and founder effects-and that there is no balanced co-presence of two main photobionts as previously established in R. farinacea. The sampling design performed here minimizes bias in the assessment of photobiont diversity in lichens and is proposed to be reliable and applicable to further study microalgal diversity in lichen symbioses.

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.

Exploring syntenic conservation across genomes for phylogenetic studies of organisms subjected to horizontal gene transfers: A case study with Cyanobacteria and cyanolichens

Understanding the evolutionary history of symbiotic Cyanobacteria at a fine scale is essential to unveil patterns of associations with their hosts and factors driving their spatiotemporal interactions. As for bacteria in general, Horizontal Gene Transfers (HGT) are expected to be rampant throughout their evolution, which justified the use of single-locus phylogenies in macroevolutionary studies of these photoautotrophic bacteria. Genomic approaches have greatly increased the amount of molecular data available, but the selection of orthologous, congruent genes that are more likely to reflect bacterial macroevolutionary histories remains problematic. In this study, we developed a synteny-based approach and searched for Collinear Orthologous Regions (COR), under the assumption that genes that are present in the same order and orientation across a wide monophyletic clade are less likely to have undergone HGT. We searched sixteen reference Nostocales genomes and identified 99 genes, part of 28 COR comprising three to eight genes each. We then developed a bioinformatic pipeline, designed to minimize inter-genome contamination and processed twelve Nostoc-associated lichen metagenomes. This reduced our original dataset to 90 genes representing 25 COR, which were used to infer phylogenetic relationships within Nostocales and among lichenized Cyanobacteria. This dataset was narrowed down further to 71 genes representing 22 COR by selecting only genes part of one (largest) operon per COR. We found a relatively high level of congruence among trees derived from the 90-gene dataset, but congruence was only slightly higher among genes within a COR compared to genes across COR. However, topological congruence was significantly higher among the 71 genes part of one operon per COR. Nostocales phylogenies resulting from concatenation and species tree approaches based on the 90- and 71-gene datasets were highly congruent, but the most highly supported result was obtained when using synteny, collinearity, and operon information (i.e., 71-gene dataset) as gene selection criteria, which outperformed larger datasets with more genes.

Extracting the GEMs: Genotype, Environment, and Microbiome Interactions Shaping Host Phenotypes

One of the fundamental tenets of biology is that the phenotype of an organism (Y) is determined by its genotype (G), the environment (E), and their interaction (GE). Quantitative phenotypes can then be modeled as Y = G + E + GE + e, where e is the biological variance. This simple and tractable model has long served as the basis for studies investigating the heritability of traits and decomposing the variability in fitness. The importance and contribution of microbe interactions to a given host phenotype is largely unclear, nor how this relates to the traditional GE model. Here we address this fundamental question and propose an expansion of the original model, referred to as GEM, which explicitly incorporates the contribution of the microbiome (M) to the host phenotype, while maintaining the simplicity and tractability of the original GE model. We show that by keeping host, environment, and microbiome as separate but interacting variables, the GEM model can capture the nuanced ecological interactions between these variables. Finally, we demonstrate with an in vitro experiment how the GEM model can be used to statistically disentangle the relative contributions of each component on specific host phenotypes.

Heavy-metal pollution induces changes in the genetic composition and anatomical properties of photobionts in pioneer lichens colonising post-industrial habitats

Certain lichens are effective colonisers of polluted sites. However, little is known about the tolerance of photobionts and the degree of mycobiont selectivity to photobionts relative to metal pollution. The present study recognises the genetic and anatomical diversity of Asterochloris photobionts in epigeic lichens, i.e. Cladonia cariosa, C. rei, and Diploschistes muscorum, in relation to a wide spectrum of soil pollution. In accordance with phylogenetic analysis, photobionts were clustered in 7 moderately- to well-supported clades, including 19 haplotypes. The mycobionts of all studied lichens demonstrated a low level of selectivity and were capable of associating with various Asterochloris lineages. This tendency was also expressed by the frequent (~25%) occurrence of multiple algal genotypes in a single thallus. This indicates that identified Asterochloris lineages are generally tolerant to heavy-metal pollution, and the low level of selectivity of mycobionts enables them to select the most suitable and/or available partner. The trend of increasing incidence of certain Asterochloris lineages and decreasing frequency of others along with increasing soil pollution was observed. This proves the superior adaptation of some photobionts to polluted sites. Such symbiotic plasticity constitute an adaptive feature necessary for the successful colonisation. High number of haplotypes at polluted sites could be the result of multiple introduction events from different areas during the initial stages of spontaneous succession. Regardless of the genetic pattern, Asterochloris cells were considerably smaller, and the density and compaction of cells in the algal layer were higher, in lichen specimens from polluted sites, indicating that photobiont characteristics may be closely dependent on heavy-metal pollution.

Extensive photobiont sharing in a rapidly radiating cyanolichen clade

Recent studies have uncovered remarkable diversity in Dictyonema s.lat. basidiolichens, here recognized as subtribe Dictyonemateae. This group includes five genera and 148 species, but hundreds more await description. The photobionts of these lichens belong to Rhizonema, a recently resurrected cyanobacterial genus known by a single species. To further investigate photobiont diversity within Dictyonemateae, we generated 765 new cyanobacterial sequences from 635 specimens collected from 18 countries. The ITS barcoding locus supported the recognition of 200 mycobiont (fungal) species among these samples, but the photobiont diversity was comparatively low. Our analyses revealed three main divisions of Rhizonema, with two repeatedly recovered as monophyletic (proposed as new species), and the third mostly paraphyletic. The paraphyletic lineage corresponds to R. interruptum and partnered with mycobionts from all five genera in Dictyonemateae. There was no evidence of photobiont-mycobiont co-speciation, but one of the monophyletic lineages of Rhizonema appears to partner predominantly with one of the two major clades of Cora (mycobiont) with samples collected largely from the northern Andes. Molecular clock estimations indicate the Rhizonema species are much older than the fungal species in the Dictyonemateae, suggesting that these basidiolichens obtained their photobionts from older ascolichen lineages and the photobiont variation in extant lineages of Dictyonemateae is the result of multiple photobiont switches. These results support the hypothesis of lichens representing "fungal farmers," in which diverse mycobiont lineages associate with a substantially lower diversity of photobionts by sharing those photobionts best suited for the lichen symbiosis among multiple and often unrelated mycobiont lineages.

Lichen-forming fungi in postindustrial habitats involve alternative photobionts

Mycobionts of many lichen genera appear to demonstrate strong selectivity in the choice of algal partner. The biological properties of a photobiont and its availability in an environment significantly determine the habitat requirements of lichens. Flexibility in photobiont choice extends the ecological amplitude of lichens; therefore, it may constitute an important adaptive strategy for colonization of extreme habitats. The photobiont inventory of the three epigeic lichens most resistant to soil pollution, i.e., Cladonia cariosa, C. rei, and the hyperaccumulator Diploschistes muscorum, was examined to verify whether and to what extent algal composition depends on the type of habitat and substrate enrichment with heavy metals. Photobionts Asterochloris and Trebouxia were identified in the studied lichen species; however, the presence of Trebouxia was directly related to anthropogenic sites with technogenic substrates, and the proportion of lichen specimens with these algae clearly depended on the level of heavy-metal soil pollution and the habitat type. The total number of algal haplotypes increased with increasing soil pollution, and the richness was associated more with soil pollution than with a given lichen species. Additionally, a large number of lichen individuals bearing multiple algal genotypes at polluted sites were recorded. Although Cladonia lichens were previously thought to be restricted to Asterochloris, they are able to start the relichenization process with Trebouxia under specific habitat conditions and to establish a stable association with these algae when colonization of disturbed sites takes place. Comparative analysis of the internal transcribed spacer (ITS) rDNA sequences revealed as many as 13 haplotypes of Trebouxia, and phylogenetic analysis grouped them into two different clades. Such a high level of genetic diversity indicates that Trebouxia is well adapted to metal pollution and could be an alternative photosynthetic partner for certain lichens, especially in polluted sites.

Myco- and photobiont associations in crustose lichens in the McMurdo Dry Valleys (Antarctica) reveal high differentiation along an elevational gradient

Climatically extreme regions such as the polar deserts of the McMurdo Dry Valleys (78° S) in Continental Antarctica are key areas for a better understanding of changes in ecosystems. Therefore, it is particularly important to analyze and communicate current patterns of biodiversity in these sensitive areas, where precipitation mostly occurs in form of snow and liquid water is rare. Humidity provided by dew, clouds, and fog are the main water sources, especially for rock-dwelling crustose lichens as one of the most common vegetation-forming organisms. We investigated the diversity and interaction specificity of myco-/photobiont associations of 232 crustose lichen specimens, collected along an elevational gradient (171–959 m a.s.l.) within the McMurdo Dry Valleys. The mycobiont species and photobiont OTUs were identified by using three markers each (nrITS, mtSSU, RPB1, and nrITS, psbJ-L, COX2). Elevation, positively associated with water availability, turned out to be the key factor explaining most of the distribution patterns of the mycobionts. Pairwise comparisons showed Lecidea cancriformis and Rhizoplaca macleanii to be significantly more common at higher elevations and Carbonea vorticosa and Lecidea polypycnidophora at lower elevations. Lichen photobionts were dominated by the globally distributed Trebouxia OTU, Tr_A02 which occurred at all habitats. Network specialization resulting from myco-/photobiont bipartite network structure varied with elevation and associated abiotic factors. Along an elevational gradient, the spatial distribution, diversity, and genetic variability of the lichen symbionts appear to be mainly influenced by improved water relations at higher altitudes.

Symbiotic microalgal diversity within lichenicolous lichens and crustose hosts on Iberian Peninsula gypsum biocrusts

This study analyses the interactions among crustose and lichenicolous lichens growing on gypsum biocrusts. The selected community was composed of Acarospora nodulosa, Acarospora placodiiformis, Diploschistes diacapsis, Rhizocarpon malenconianum and Diplotomma rivas-martinezii. These species represent an optimal system for investigating the strategies used to share phycobionts because Acarospora spp. are parasites of D. diacapsis during their first growth stages, while in mature stages, they can develop independently. R. malenconianum is an obligate lichenicolous lichen on D. diacapsis, and D. rivas-martinezii occurs physically close to D. diacapsis. Microalgal diversity was studied by Sanger sequencing and 454-pyrosequencing of the nrITS region, and the microalgae were characterized ultrastructurally. Mycobionts were studied by performing phylogenetic analyses. Mineralogical and macro- and micro-element patterns were analysed to evaluate their influence on the microalgal pool available in the substrate. The intrathalline coexistence of various microalgal lineages was confirmed in all mycobionts. D. diacapsis was confirmed as an algal donor, and the associated lichenicolous lichens acquired their phycobionts in two ways: maintenance of the hosts' microalgae and algal switching. Fe and Sr were the most abundant microelements in the substrates but no significant relationship was found with the microalgal diversity. The range of associated phycobionts are influenced by thallus morphology.

Disease swamps molecular signatures of genetic-environmental associations to abiotic factors in Tasmanian devil (Sarcophilus harrisii) populations

Landscape genomics studies focus on identifying candidate genes under selection via spatial variation in abiotic environmental variables, but rarely by biotic factors (i.e., disease). The Tasmanian devil (Sarcophilus harrisii) is found only on the environmentally heterogeneous island of Tasmania and is threatened with extinction by a transmissible cancer, devil facial tumor disease (DFTD). Devils persist in regions of long-term infection despite epidemiological model predictions of species' extinction, suggesting possible adaptation to DFTD. Here, we test the extent to which spatial variation and genetic diversity are associated with the abiotic environment (i.e., climatic variables, elevation, vegetation cover) and/or DFTD. We employ genetic-environment association analyses using 6886 SNPs from 3287 individuals sampled pre- and post-disease arrival across the devil's geographic range. Pre-disease, we find significant correlations of allele frequencies with environmental variables, including 365 unique loci linked to 71 genes, suggesting local adaptation to abiotic environment. The majority of candidate loci detected pre-DFTD are not detected post-DFTD arrival. Several post-DFTD candidate loci are associated with disease prevalence and were in linkage disequilibrium with genes involved in tumor suppression and immune response. Loss of apparent signal of abiotic local adaptation post-disease suggests swamping by strong selection resulting from the rapid onset of DFTD.

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.

Differential Colonization and Succession of Microbial Communities in Rock and Soil Substrates on a Maritime Antarctic Glacier Forefield

Glacier forefields provide a unique chronosequence to assess microbial or plant colonization and ecological succession on previously uncolonized substrates. Patterns of microbial succession in soils of alpine and subpolar glacier forefields are well documented but those affecting high polar systems, including moraine rocks, remain largely unexplored. In this study, we examine succession patterns in pioneering bacterial, fungal and algal communities developing on moraine rocks and soil at the Hurd Glacier forefield (Livingston Island, Antarctica). Over time, changes were produced in the microbial community structure of rocks and soils (ice-free for different lengths of time), which differed between both substrates across the entire chronosequence, especially for bacteria and fungi. In addition, fungal and bacterial communities showed more compositional consistency in soils than rocks, suggesting community assembly in each niche could be controlled by processes operating at different temporal and spatial scales. Microscopy revealed a patchy distribution of epilithic and endolithic lithobionts, and increasing endolithic colonization and microbial community complexity along the chronosequence. We conclude that, within relatively short time intervals, primary succession processes at polar latitudes involve significant and distinct changes in edaphic and lithic microbial communities associated with soil development and cryptogamic colonization.

Morphology and secondary chemistry in species recognition of Parmelia omphalodes group - evidence from molecular data with notes on the ecological niche modelling and genetic variability of photobionts

To evaluate the importance of morphological and chemical characters used in the recognition of species within the Parmelia omphalodes group, we performed phylogenetic, morphological and chemical analyses of 335 specimens, of which 34 were used for molecular analyses. Phylogenetic analyses, based on ITS rDNA sequences, show that P. pinnatifida is distinct from P. omphalodes and the most important difference between those species is the development of pseudocyphellae. In P. pinnatifida, they are mostly marginal and form white rims along lobes margins, but laminal pseudocyphellae can develop in older parts of thalli and are predominantly connected with marginal pseudocyphellae. In contrast, in P. omphalodes laminal pseudocyphellae are common and are predominantly not connected to marginal pseudocyphellae. Chemical composition of secondary lichen metabolites in both analysed species is identical and therefore this feature is not diagnostic in species recognition. Few samples of P. discordans, species morphologically similar to P. omphalodes and P. pinnatifida, were also included in the analyses and they are nested within the clade of P. omphalodes, despite the different chemistry (protocetraric acid present versus salazinic acid in P. omphalodes). All taxa of the P. omphalodes group occupy similar niches, but their potential distributions are wider than those currently known. The absence of specimens in some localities may be limited by the photobiont availability. Parmelia omphalodes and P. pinnatifida are moderately selective in photobiont choice as they form associations with at least two or three lineages of Trebouxia clade S. Parmelia pinnatifida, as well as P. discordans are associated with Trebouxia OTU S02 which seems to have a broad ecological amplitude. Other lineages of Trebouxia seem to be rarer, especially Trebouxia sp. OTU S04, which is sometimes present in P. pinnatifida. This study indicates the importance of extensive research including morphology, chemistry and analysis of molecular markers of both bionts in taxonomical studies of lichens.

A Glimpse into Genetic Diversity and Symbiont Interaction Patterns in Lichen Communities from Areas with Different Disturbance Histories in Białowieża Forest, Poland

Anthropogenic disturbances can have strong impacts on lichen communities, as well as on individual species of lichenized fungi. Traditionally, lichen monitoring studies are based on the presence and abundance of fungal morphospecies. However, the photobionts, as well photobiont mycobiont interactions also contribute to the structure, composition, and resilience of lichen communities. Here we assess the genetic diversity and interaction patterns of algal and fungal partners in lichen communities along an anthropogenic disturbance gradient in Białowieża Forest (Poland). We sampled a total of 224 lichen thalli in a protected, a managed, and a disturbed area of the forest, and sequenced internal transcribed spacer (ITS) ribosomal DNA (rDNA) of both, fungal and algal partners. Sequence clustering using a 97% similarity threshold resulted in 46 fungal and 23 green algal operational taxonomic units (OTUs). Most of the recovered photobiont OTUs (14 out of 23) had no similar hit in the NCBI-BLAST search, suggesting that even in well studied regions, such as central Europe, a lot of photobiont diversity is yet undiscovered. If a mycobiont was present at more than one site, it was typically associated with the same photobiont OTU(s). Generalist species, i.e., taxa that associate with multiple symbiont partners, occurred in all three disturbance regimes, suggesting that such taxa have few limitations in colonizing or persisting in disturbed areas. Trebouxia jamesii associated with 53% of the fungal OTUs, and was generally the most common photobiont OTU in all areas, implying that lichens that associate with this symbiont are not limited by the availability of compatible photobionts in Central European forests, regardless of land use intensity.

Climate Change, Bioclimatic Models and the Risk to Lichen Diversity

This paper provides an overview of bioclimatic models applied to lichen species, supporting their potential use in this context as indicators of climate change risk. First, it provides a brief summary of climate change risk, pointing to the relevance of lichens as a topic area. Second, it reviews the past use of lichen bioclimatic models, applied for a range of purposes with respect to baseline climate, and the application of data sources, statistical methods, model extents and resolution and choice of predictor variables. Third, it explores additional challenges to the use of lichen bioclimatic models, including: 1. The assumption of climatically controlled lichen distributions, 2. The projection to climate change scenarios, and 3. The issue of nonanalogue climates and model transferability. Fourth, the paper provides a reminder that bioclimatic models estimate change in the extent or range of a species suitable climate space, and that an outcome will be determined by vulnerability responses, including potential for migration, adaptation, and acclimation, within the context of landscape habitat quality. The degree of exposure to climate change, estimated using bioclimatic models, can help to inform an understanding of whether vulnerability responses are sufficient for species resilience. Fifth, the paper draws conclusions based on its overview, highlighting the relevance of bioclimatic models to conservation, support received from observational data, and pointing the way towards mechanistic approaches that align with field-scale climate change experiments.

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.