Surveying lichen diversity in forests: A comparison of expert mapping and eDNA metabarcoding of bark surfaces.

Lichens are an important part of forest ecosystems, contributing to forest biodiversity, the formation of micro-niches and nutrient cycling. Assessing the diversity of lichenised fungi in complex ecosystems, such as forests, requires time and substantial skills in collecting and identifying lichens. The completeness of inventories thus largely depends on the expertise of the collector, time available for the survey and size of the studied area. Molecular methods of surveying biodiversity hold the promise to overcome these challenges. DNA barcoding of individual lichen specimens and bulk collections is already being applied; however, eDNA methods have not yet been evaluated as a tool for lichen surveys. Here, we assess which species of lichenised fungi can be detected in eDNA swabbed from bark surfaces of living trees in central European forests. We compare our findings to an expert floristic survey carried out in the same plots about a decade earlier. In total, we studied 150 plots located in three study regions across Germany. In each plot, we took one composite sample based on six trees, belonging to the species Fagussylvatica, Piceaabies and Pinussylvestris. The eDNA method yielded 123 species, the floristic survey 87. The total number of species found with both methods was 167, of which 48% were detected only in eDNA, 26% only in the floristic survey and 26% in both methods. The eDNA contained a higher diversity of inconspicuous species. Many prevalent taxa reported in the floristic survey could not be found in the eDNA due to gaps in molecular reference databases. We conclude that, currently, eDNA has merit as a complementary tool to monitor lichen biodiversity at large scales, but cannot be used on its own. We advocate for the further development of specialised and more complete databases.

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

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

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.

Biosynthetic gene cluster synteny: Orthologous polyketide synthases in Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata.

Lichens are symbiotic associations consisting of a photobiont (algae or cyanobacteria) and a mycobiont (fungus), which together generate a variety of unique secondary metabolites. To access this biosynthetic potential for biotechnological applications, deeper insights into the biosynthetic pathways and corresponding gene clusters are necessary. Here, we provide a comparative view of the biosynthetic gene clusters of three lichen mycobionts derived from Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata. In addition, we present a high-quality PacBio metagenome of Parmelia sulcata, from which we extracted the mycobiont bin containing 214 biosynthetic gene clusters. Most biosynthetic gene clusters in these genomes were associated with T1PKSs, followed by NRPSs and terpenes. This study focused on biosynthetic gene clusters related to polyketide synthesis. Based on ketosynthase homology, we identified nine highly syntenic clusters present in all three species. Among the four clusters belonging to nonreducing PKSs, two are putatively linked to lichen substances derived from orsellinic acid (orcinol depsides and depsidones, e.g., lecanoric acid, physodic acid, lobaric acid), one to compounds derived from methylated forms of orsellinic acid (beta orcinol depsides, e.g., atranorin), and one to melanins. Five clusters with orthologs in all three species are linked to reducing PKSs. Our study contributes to sorting and dereplicating the vast PKS diversity found in lichenized fungi. High-quality sequences of biosynthetic gene clusters of these three common species provide a foundation for further exploration into biotechnological applications and the molecular evolution of lichen substances.

Biosynthetic Potential of Hypogymnia Holobionts: Insights into Secondary Metabolite Pathways.

Lichens are symbiotic associations consisting of a photobiont (algae or cyanobacteria) and a mycobiont (fungus). They are known to produce a variety of unique secondary metabolites. To access this biosynthetic potential for biotechnological applications, deeper insights into the biosynthetic pathways and corresponding gene clusters are necessary. Here we provide a comprehensive view of the biosynthetic gene clusters of all organisms comprising a lichen thallus: fungi, green algae, and bacteria. We present two high-quality PacBio metagenomes, in which we identified a total of 460 biosynthetic gene clusters. Lichen mycobionts yielded 73-114 clusters, other lichen associated ascomycetes 8-40, green algae of the genus Trebouxia 14-19, and lichen-associated bacteria 101-105 clusters. The mycobionts contained mainly T1PKSs, followed by NRPSs, and terpenes; Trebouxia reads harbored mainly clusters linked to terpenes, followed by NRPSs and T3PKSs. Other lichen-associated ascomycetes and bacteria contained a mix of diverse biosynthetic gene clusters. In this study, we identified for the first time the biosynthetic gene clusters of entire lichen holobionts. The yet untapped biosynthetic potential of two species of the genus Hypogymnia is made accessible for further research.

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.

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

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

Lichen holobionts show compositional structure along elevation.

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

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

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

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.

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.

A Candidate Gene Cluster for the Bioactive Natural Product Gyrophoric Acid in Lichen-Forming Fungi.

Natural products of lichen-forming fungi are structurally diverse and have a variety of medicinal properties. Despite this, they have limited implementation in industry mostly because the corresponding genes are unknown for most of their natural products. Here, we implement a long-read sequencing and bioinformatic approach to identify the putative biosynthetic gene cluster of the bioactive natural product gyrophoric acid (GA). Using 15 high-quality genomes representing nine GA-producing species of the lichen-forming fungal genus Umbilicaria, we identify the most likely GA cluster and investigate the cluster gene organization and composition across the nine species. Our results show that GA clusters are promiscuous within Umbilicaria, and only three genes are conserved across species, including the polyketide synthase (PKS) gene. In addition, our results suggest that the same cluster codes for different, but structurally similar compounds, namely, GA, umbilicaric-, and hiascic acid, bringing new evidence that lichen metabolite diversity is also generated through regulatory mechanisms at the molecular level. Ours is the first study to identify the most likely GA cluster and, thus, provides essential information to open new avenues for biotechnological approaches to producing and modifying GA and similar lichen-derived compounds. GA PKS is the first tridepside PKS to be identified. IMPORTANCE The implementation of natural products in the pharmaceutical industry relies on the possibility of modifying the natural product (NP) pathway to optimize yields and pharmacological effects. Characterization of genes and pathways underlying natural product biosynthesis is a major bottleneck for exploiting the medicinal properties of the natural products. Genome mining is a promising and relatively cost- and time-effective approach to utilize unexplored NP resources for drug discovery. In this study, we identify the most likely gene cluster for the lichen-forming fungal depside gyrophoric acid in nine Umbilicaria species. This compound shows cytotoxic and antiproliferative properties against several cancer cell lines and is also a broad-spectrum antimicrobial agent. This information paves the way for generating GA analogs with modified properties by selective activation/deactivation of genes.

Lecanoric acid mediates anti-proliferative effects by an M phase arrest in colon cancer cells.

Lichen extracts containing, among other compounds, depsides such as evernic acid, atranorin, and lecanoric acid possess anti-proliferative effects. We aimed to identify lichen metabolites that are responsible for the observed anti-proliferative effects. We performed cytotoxicity, cell colony, cell cycle and apoptosis assays in various cell lines or primary immune cells. We analyzed several cell cycle proteins and apoptosis-related proteins to gain insights into the underlying mechanism. All depsides reduced the viability of the tested cell lines (HCT-116, HEK293T, HeLa, NIH3T3, RAW246.7) in a cell line-dependent manner with lecanoric acid being the most effective. Atranorin did not influence the cell cycle or colony formation in HCT-116 cells, but induced apoptosis in HCT-116 cells. Evernic acid showed no anti-proliferative effects. Lecanoric acid inhibited cell colony formation already at 0.03 µg/ml in HCT-116 cells and induced a G2 cell cycle block in several cell lines. Moreover, lecanoric acid arrested the cell cycle, presumably in the M phase, since expression of cyclin B1 and phosphorylated histone H3 was upregulated, whereas the inactive cyclin-dependent kinase 1 (CDK1) was reduced in HCT-116 cells. Most importantly, cell death induced by lecanoric acid was more prominent in cancer cells than in primary human immune and endothelial cells. In conclusion, lecanoric acid seems to mediate its anti-proliferative effects via arrest of cells in the M phase. Our data suggest lecanoric acid may be a potential new candidate for anti-cancer therapy, because it has anti-proliferative effects on cancer cell lines, and does not affect primary immune cells.

Genome mining as a biotechnological tool for the discovery of novel biosynthetic genes in lichens.

Natural products (NPs) and their derivatives are a major contributor to modern medicine. Historically, microorganisms such as bacteria and fungi have been instrumental in generating drugs and lead compounds because of the ease of culturing and genetically manipulating them. However, the ever-increasing demand for novel drugs highlights the need to bioprospect previously unexplored taxa for their biosynthetic potential. Next-generation sequencing technologies have expanded the range of organisms that can be explored for their biosynthetic content, as these technologies can provide a glimpse of an organism's entire biosynthetic landscape, without the need for cultivation. The entirety of biosynthetic genes can be compared to the genes of known function to identify the gene clusters potentially coding for novel products. In this study, we mine the genomes of nine lichen-forming fungal species of the genus Umbilicaria for biosynthetic genes, and categorize the biosynthetic gene clusters (BGCs) as "associated product structurally known" or "associated product putatively novel". Although lichen-forming fungi have been suggested to be a rich source of NPs, it is not known how their biosynthetic diversity compares to that of bacteria and non-lichenized fungi. We found that 25%-30% of biosynthetic genes are divergent as compared to the global database of BGCs, which comprises 1,200,000 characterized biosynthetic genes from plants, bacteria, and fungi. Out of 217 BGCs, 43 were highly divergant suggesting that they potentially encode structurally and functionally novel NPs. Clusters encoding the putatively novel metabolic diversity comprise polyketide synthases (30), non-ribosomal peptide synthetases (12), and terpenes (1). Our study emphasizes the utility of genomic data in bioprospecting microorganisms for their biosynthetic potential and in advancing the industrial application of unexplored taxa. We highlight the untapped structural metabolic diversity encoded in the lichenized fungal genomes. To the best of our knowledge, this is the first investigation identifying genes coding for NPs with potentially novel properties in lichenized fungi.

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.

Depside and Depsidone Synthesis in Lichenized Fungi Comes into Focus through a Genome-Wide Comparison of the Olivetoric Acid and Physodic Acid Chemotypes of Pseudevernia furfuracea.

Primary biosynthetic enzymes involved in the synthesis of lichen polyphenolic compounds depsides and depsidones are non-reducing polyketide synthases (NR-PKSs), and cytochrome P450s. However, for most depsides and depsidones the corresponding PKSs are unknown. Additionally, in non-lichenized fungi specific fatty acid synthases (FASs) provide starters to the PKSs. Yet, the presence of such FASs in lichenized fungi remains to be investigated. Here we implement comparative genomics and metatranscriptomics to identify the most likely PKS and FASs for olivetoric acid and physodic acid biosynthesis, the primary depside and depsidone defining the two chemotypes of the lichen Pseudevernia furfuracea. We propose that the gene cluster PF33-1_006185, found in both chemotypes, is the most likely candidate for the olivetoric acid and physodic acid biosynthesis. This is the first study to identify the gene cluster and the FAS likely responsible for olivetoric acid and physodic acid biosynthesis in a lichenized fungus. Our findings suggest that gene regulation and other epigenetic factors determine whether the mycobiont produces the depside or the depsidone, providing the first direct indication that chemotype diversity in lichens can arise through regulatory and not only through genetic diversity. Combining these results and existing literature, we propose a detailed scheme for depside/depsidone synthesis.

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.

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.

Unraveling the Pharmacological Potential of Lichen Extracts in the Context of Cancer and Inflammation With a Broad Screening Approach.

Lichen-forming fungi are symbiotic organisms that synthesize unique natural products with potential for new drug leads. Here, we explored the pharmacological activity of six lichen extracts (Evernia prunastri, Pseudevernia furfuracea, Umbilicaria pustulata, Umbilicaria crustulosa, Flavoparmelia caperata, Platismatia glauca) in the context of cancer and inflammation using a comprehensive set of 11 functional and biochemical in vitro screening assays. We assayed intracellular Ca2+ levels and cell migration. For cancer, we measured tumor cell proliferation, cell cycle distribution and apoptosis, as well as the angiogenesis-associated proliferation of endothelial cells (ECs). Targeting inflammation, we assayed leukocyte adhesion onto ECs, EC adhesion molecule expression, as well as nitric oxide production and prostaglandin (PG)E2 synthesis in leukocytes. Remarkably, none of the lichen extracts showed any detrimental influence on the viability of ECs. We showed for the first time that extracts of F. caperata induce Ca2+ signaling. Furthermore, extracts from E. prunastri, P. furfuracea, F. caperata, and P. glauca reduced cell migration. Interestingly, F. caperata extracts strongly decreased tumor cell survival. The proliferation of ECs was significantly reduced by E. prunastri, P. furfuracea, and F. caperata extracts. The extracts did not inhibit the activity of inflammatory processes in ECs. However, the pro-inflammatory activation of leukocytes was inhibited by extracts from E. prunastri, P. furfuracea, F. caperata, and P. glauca. After revealing the potential biological activities of lichen extracts by an array of screening tests, a correlation analysis was performed to evaluate particular roles of abundant lichen secondary metabolites, such as atranorin, physodic acid, and protocetraric acid as well as usnic acid in various combinations. Overall, some of the lichen extracts tested in this study exhibit significant pharmacological activity in the context of inflammation and/or cancer, indicating that the group lichen-forming fungi includes promising members for further testing.

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.