Microclimatic differentiation of gene pools in theLobaria pulmonariasymbiosis in a primeval forest landscape

The microbiome of the lichen Lobaria pulmonaria varies according to climate on a subcontinental scale

The Lobaria pulmonaria holobiont comprises algal, fungal, cyanobacterial and bacterial components. We investigated L. pulmonaria's bacterial microbiome in the adaptation of this ecologically sensitive lichen species to diverse climatic conditions. Our central hypothesis posited that microbiome composition and functionality aligns with subcontinental-scale (a stretch of ~1100 km) climatic parameters related to temperature and precipitation. We also tested the impact of short-term weather dynamics, sampling season and algal/fungal genotypes on microbiome variation. Metaproteomics provided insights into compositional and functional changes within the microbiome. Climatic variables explained 41.64% of microbiome variation, surpassing the combined influence of local weather and sampling season at 31.63%. Notably, annual mean temperature and temperature seasonality emerged as significant climatic drivers. Microbiome composition correlated with algal, not fungal genotype, suggesting similar environmental recruitment for the algal partner and microbiome. Differential abundance analyses revealed distinct protein compositions in Sub-Atlantic Lowland and Alpine regions, indicating differential microbiome responses to contrasting environmental/climatic conditions. Proteins involved in oxidative and cellular stress were notably different. Our findings highlight microbiome plasticity in adapting to stable climates, with limited responsiveness to short-term fluctuations, offering new insights into climate adaptation in lichen symbiosis.

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.

Population genetics and biogeography of the lungwort lichen in North America support distinct Eastern and Western gene pools

PREMISE

Populations of species with large spatial distributions are shaped by complex forces that differ throughout their ranges. To maintain the genetic diversity of species, genepool-based subsets of widespread species must be considered in conservation assessments.

METHODS

The population genetics of the lichenized fungus Lobaria pulmonaria and its algal partner, Symbiochloris reticulata, were investigated using microsatellite markers to determine population structure, genetic diversity, and degree of congruency in eastern and western North America. Data loggers measuring temperature and humidity were deployed at selected populations in eastern North America to test for climatic adaptation. To better understand the role Pleistocene glaciations played in shaping population patterns, a North American, range-wide species distribution model was constructed and hindcast to 22,000 years before present and at 500-year time slices from then to the present.

RESULTS

The presence of two gene pools with minimal admixture was supported, one in the U.S. Pacific Northwest and one in eastern North America. Western populations were significantly more genetically diverse than eastern populations. There was no evidence for climatic adaptation among eastern populations, though there was evidence for range-wide adaptation to evapotranspiration rates. Hindcast distribution models suggest that observed genetic diversity may be due to a drastic Pleistocene range restriction in eastern North America, whereas a substantial coastal refugial area is inferred in the west.

CONCLUSIONS

Taken together the results show different, complex population histories of L. pulmonaria in eastern and western North America, and suggest that conservation planning for each gene pool should be considered separately.

Deep divergence between island populations in lichenized fungi

Macaronesia is characterized by a high degree of endemism and represents a noteworthy system to study the evolutionary history of populations and species. Here, we compare the population-genetic structure in three lichen-forming fungi, the widespread Lobaria pulmonaria and two Macaronesian endemics, L. immixta and L. macaronesica, based on microsatellites. We utilize population genetic approaches to explore population subdivision and evolutionary history of these taxa on the Canary Islands, Madeira, Azores, and the western Iberian Peninsula. A common feature in all species was the deep divergence between populations on the Azores, a pattern expected by the large geographic distance among islands. For both endemic species, there was a major split between archipelagos. In contrast, in the widespread L. pulmonaria, divergent individuals were distributed across multiple archipelagos, suggesting a complex evolutionary history involving repeated migration between islands and mainland.

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.

Low genetic differentiation between apotheciate Usnea florida and sorediate Usnea subfloridana (Parmeliaceae, Ascomycota) based on microsatellite data

Accurate species delimitation has a pivotal role in conservation biology, and it is especially important for threatened species where decisions have political and economic consequences. Finding and applying appropriate character sets and analytical tools to resolve interspecific relationships remains challenging in lichenized fungi. The main aim of our study was to re-assess the species boundaries between Usnea subfloridana and Usnea florida, which have been phylogenetically indistinguishable until now, but are different in reproductive mode and ecological preferences, using fungal-specific simple sequence repeats (SSR), i.e. microsatellite markers. Bayesian clustering analysis, discriminant analysis of principal components (DAPC), minimal spanning network (MSN), and principal component analysis (PCA) failed to separate U. florida and U. subfloridana populations. However, a low significant differentiation between the two taxa was observed across all populations according to AMOVA results. Also, analysis of shared haplotypes and statistical difference in clonal diversity (M) supported the present-day isolation between the apotheciate U. florida and predominantly sorediate U. subfloridana. Our results do not provide a clear support either for the separation of species in this pair or the synonymization of U. florida and U. subfloridana. We suggest that genome-wide data could help resolve the taxonomic question in this species pair.

Distribution of Foliicolous Lichen Strigula and Genetic Structure of S. multiformis on Jeju Island, South Korea

Strigula is a pantropic foliicolous lichen living on the leaf surfaces of evergreen broadleaf plants. In South Korea, Strigula is the only genus of foliicolous lichen recorded from Jeju Island. Several Strigula species have been recorded, but the ecology of Strigula in South Korea has been largely unexplored. This study examined the distribution and genetic structure of Strigula on Jeju Island. The distribution was surveyed and the influence of environmental factors (e.g., elevation, forest availability, and bioclimate) on the distribution was analyzed using a species distribution modeling analysis. In addition, the genetic variations and differentiation of Strigula multiformis populations were analyzed using two nuclear ribosomal regions. The distribution of Strigula was largely restricted to a small portion of forest on Jeju Island, and the forest availability was the most important factor in the prediction of potential habitats. The genetic diversity and differentiation of the S. multiformis population were found to be high and were divided according to geography. On the other hand, geographic and environmental distance did not explain the population differentiation. Distribution and population genetic analysis suggested that the available habitat and genetic exchange of Strigula on Jeju Island are limited by the lack of available forest in the lowlands.

Microsatellite based genetic diversity of the widespread epiphytic lichen Usnea subfloridana (Parmeliaceae, Ascomycota) in Estonia: comparison of populations from the mainland and an island

Understanding the distribution of genetic patterns and structure is an essential target in population genetics and, thereby, important for conservation genetics. The main aim of our study was to investigate the population genetics of Usnea subfloridana, a widespread lichenised fungus, focusing on a comparison of genetic variation of its populations amongst three geographically remote and disconnected regions, in order to determine relationships amongst environmental data, variation in lichen secondary chemistry and microsatellite data in genotyped populations. In all, 928 Usnea thalli from 17 populations were genotyped using seven specific fungal microsatellite markers. Different measures of genetic diversity (allelic richness, private allelic richness, Nei's unbiased genetic diversity and clonal diversity) were calculated and compared between lichen populations. Our results revealed a low genetic differentiation of U. subfloridana populations amongst three distant areas in Estonia and also a high level of gene flow. The results support suggestion of the long-range vegetative dispersal of subpendulous U. subfloridana via symbiotic propagules (soralia, isidia or fragments of thalli). Our study has also provided evidence that environmental variables, including mean annual temperature and geographical longitude, shape the genetic structure of U. subfloridana populations in Estonia. Additionally, a weak but statistically significant correlation between lichen chemotypes and microsatellite allele distribution was found in genotyped specimens.

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.

Climate-woodland effects on population genetics for two congeneric lichens with contrasting reproductive strategies

Genetic variation is expected to be influenced by the interaction between reproductive mode and dispersal traits on the one hand and environmental and habitat setting affecting establishment success on the other. We evaluated how environmental/habitat setting affects population genetic variation (i.e. variation in genetic diversity and structure) when regulated by contrasting dispersal traits. We used fungus-specific microsatellite markers to examine genetic diversity and structure of two closely related epiphytic lichen fungi that differ in their primary reproductive mode: Nephroma laevigatum (sexually reproducing, n = 191, 10 microsatellites) and N. parile (asexually reproducing, n = 182, 12 microsatellites), along a steep climatic gradient in Scotland. Despite their reproductive differences, we found a high proportion of clones in both species and a background pattern of genetic structure related to climatic gradients. We also demonstrated that woodland connectivity, rather than geographic distance, explained genetic diversity in both species. Environmental/habitat setting, modulated by the reproductive mode of the species, affects genetic diversity and structure, but the putative dissimilarity in their reproductive mode is less important than has been previously assumed. We reinforce the importance of protecting highly connected populations, positioned along a gradient capturing the segregation of gene pool differences in response to climatic variation.

First genome-wide analysis of the endangered, endemic lichen Cetradonia linearis reveals isolation by distance and strong population structure

PREMISE OF THE STUDY

Lichenized fungi are evolutionarily diverse and ecologically important, but little is known about the processes that drive their diversification and genetic differentiation. Distributions are often assumed to be wholly shaped by ecological requirements rather than dispersal limitations. Furthermore, although asexual and sexual reproductive structures are observable, the lack of information about recombination rates makes inferences about reproductive strategies difficult. We investigated the population genomics of Cetradonia linearis, a federally endangered lichen in the southern Appalachians of eastern North America, to test the relative contributions of environmental and geographic distance in shaping genetic structure, and to characterize the mating system and genome-wide recombination.

METHODS

Whole-genome shotgun sequencing was conducted to generate data for 32 individuals of C. linearis. A reference genome was assembled, and reads from all samples were aligned to generate a set of single-nucleotide polymorphisms for further analyses.

KEY RESULTS

We found evidence for low rates of recombination and for isolation by distance, but not for isolation by environment. The species is putatively unisexual, given that only one mating-type locus was found. Hindcast species distribution models and the distribution of genetic diversity support C. linearis having a larger range during the Last Glacial Maximum in the southern portion of its current extent.

CONCLUSIONS

Our findings contribute to the understanding of factors that shape genetic diversity in C. linearis and in fungi more broadly. Because all populations are highly genetically differentiated, the extirpation of any population would mean the loss of unique genetic diversity; therefore, our results support the continued conservation of this species.

Strong specificity and network modularity at a very fine phylogenetic scale in the lichen genus Peltigera

Identifying the drivers and evolutionary consequences of species interactions is a major goal of community ecology. Network-based analyses can provide mathematical tools to detect non-random patterns of interactions, and potentially help predicting the consequences of such patterns on evolutionary dynamics of symbiotic systems. Here, we characterize the structure of a lichen network at a very fine phylogenetic scale, by identifying the photosynthetic partners (i.e., cyanobacteria of the genus Nostoc) of lichenized fungi belonging to a monophyletic section of a single genus (i.e., section Polydactylon of the genus Peltigera), worldwide. Even at such a fine phylogenetic scale, we found that interactions were highly modular and anti-nested, indicating strong preferences in interactions. When considering local Peltigera communities, i.e., datasets at small spatial scales with only a slightly broader phylogenetic range, interactions remained modular but were asymmetric, with generalist Nostoc partners interacting with specialized Peltigera species. This asymmetry was not detected with our global spatial scale dataset. We discuss these results in the light of lichen community assembly, and explore how such interaction patterns may influence coevolution in lichens and the evolutionary stability of the mutualism in general.

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

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

Discovery of long-distance gamete dispersal in a lichen-forming ascomycete

Accurate estimates of gamete and offspring dispersal range are required for the understanding and prediction of spatial population dynamics and species persistence. Little is known about gamete dispersal in fungi, especially in lichen-forming ascomycetes. Here, we estimate the dispersal functions of clonal propagules, gametes and ascospores of the epiphytic lichen Lobaria pulmonaria. We use hierarchical Bayesian parentage analysis, which integrates genetic and ecological information from multiannual colonization and dispersal source data collected in a large, old-growth forest landscape. The effective dispersal range of gametes is several hundred metres to kilometres from potential paternal individuals. By contrast, clonal propagules disperse only tens of metres, and ascospores disperse over several thousand metres. Our study reveals the dispersal distances of individual reproductive units; clonal propagules, gametes and ascospores, which is of great importance for a thorough understanding of the spatial dynamics of ascomycetes. Sexual reproduction occurs between distant individuals. However, whereas gametes and ascospores disperse over long distances, the overall rate of colonization of trees is low. Hence, establishment is the limiting factor for the colonization of new host trees by the lichen in old-growth landscapes.

Multi-gene phylogeny of the genus Lobaria: Evidence of species-pair and allopatric cryptic speciation in East Asia

PREMISE OF THE STUDY

Accurate species delimitation has critical implications for ecological and conservation studies. The lichen genus Lobaria is widely distributed in old-growth forests. Particularly in East Asia, this genus includes many rare and poorly known taxa that are circumscribed as morpho- or chemospecies, as well as species-pairs.

METHODS

To critically examine the relationships between species identified via morphological and chemical criteria, phylogenetic species recognition (PSR) was applied to the genus Lobaria. Morphological and chemical patterns of 87 individuals were examined and three independent nuclear loci were sequenced. The East Asian L. meridionalis-group was additionally studied using split decomposition and haplotype network analysis.

KEY RESULTS

The genus Lobaria and most of its species were strongly supported statistically. Split decomposition and haplotype networks suggest complex evolutionary histories of species within the East Asian L. meridionalis-group.

CONCLUSIONS

Phylogenetic analyses confirmed the monophyly of the genus Lobaria, including L. anomala. Within Lobaria, three major clades were found. These clades associate with different photobionts and comprise 18 known species and 5 undescribed species. Several chemical compounds were found to be neither stable nor invariant characters. Some taxa of the L. meridionalis-group appear to be monophyletic but remain as allopatric cryptic species. In three clades, this study found evidence for diversification processes between isidiate and nonisidiate specimens (species-pair). These findings are discussed in the context of evolutionary hypotheses for speciation processes.

The threatened epiphytic lichen Lobaria pulmonaria in the Iberian Peninsula: genetic diversity and structure across a latitudinal gradient

The current genetic diversity and structure of a species plays a marked role in the species' future response to environmental changes. Identification of the factors that might ensure the long-term viability of populations along its distribution area is therefore important for conserving biodiversity. In this work, infraspecific genetic diversity and structure of the threatened lichen Lobaria pulmonaria was investigated along a latitudinal gradient, spanning the Spanish latitudinal range of L. pulmonaria. Eighteen populations in Northern, Central, and Southern Spain were analysed using six specific fungal microsatellites of L. pulmonaria. Genetic diversity indices were calculated and compared among populations. Genetic differentiation was assessed using AMOVA and Bayesian methods. Additionally, a redundancy analysis was used to estimate the relative importance of environmental factors on the genetic variation among populations. Annual precipitation was the only factor affecting the genetic diversity probably through its influence on population and thallus size of L. pulmonaria, and significantly higher levels of genetic diversity were detected in southern populations. Isolation by distance was not significant, being environmental variables most important factors controlling genetic variation in L. pulmonaria populations.