Phylogeography of Ramalina menziesii, a widely distributed lichen-forming fungus in western North America

Phylogeographic analyses of an epiphytic foliose lichen show multiple dispersal events westward from the Hengduan Mountains of Yunnan into the Himalayas

Lobaria pindarensis is an endemic species of the Himalayas and the Hengduan Mountains. Little information is available on the phylogeography genetics and colonization history of this species or how its distribution patterns changed in response to the orographic history of the Himalayas and Hengduan Mountains. Based on samples covering a major part of the species' distribution range, we used 443 newly generated sequences of nine loci for molecular coalescent analyses in order to reconstruct the evolutionary history of L. pindarensis, and to reconstruct the species' ancestral phylogeographic distributions using Bayesian binary MCMC analyses. The results suggest that current populations originated from the Yunnan region of the Hengduan Mountains in the middle Pliocene, and that the Himalayas of Bhutan were colonized by a lineage that diverged from Yunnan ca. 2.72 Ma. The analysis additionally indicates that the Nepal and Xizang areas of the Himalayas were colonized from Yunnan as well, and that there was later a second dispersal event from Yunnan to Bhutan. We conclude that the change in climate and habitat related to the continuous uplift of the Himalayas and the Hengduan Mountains in the late Pliocene and middle Pleistocene influenced the geographic distribution pattern of L. pindarensis.

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.

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.

How do genes flow? Identifying potential dispersal mode for the semi-aquatic lichen Dermatocarpon luridum using spatial modelling and photobiont markers

Background

Landscape genetics is an interdisciplinary field that combines tools and techniques from population genetics with the spatially explicit principles from landscape ecology. Spatial variation in genotypes is used to test hypotheses about how landscape pattern affects dispersal in a wide range of taxa. Lichens, symbiotic associations between mycobionts and photobionts, are an entity for which little is known about their dispersal mechanism. Our objective was to infer the dispersal mechanism in the semi-aquatic lichen Dermatocarpon luridum using spatial models and the spatial variation of the photobiont, Diplosphaera chodatii. We sequenced the ITS rDNA and the β-actin gene regions of the photobiont and mapped the haplotype spatial distribution in Payuk Lake. We subdivided Payuk Lake into subpopulations and applied four spatial models based on the topography and hydrology to infer the dispersal mechanism.

Results

Genetic variation corresponded with the topography of the lake and the net flow of water through the waterbody. A lack of isolation-by-distance suggests high gene flow or dispersal within the lake. We infer the dispersal mechanism in D. luridum could either be by wind and/or water based on the haplotype spatial distribution of its photobiont using the ITS rDNA and β-actin markers.

Conclusions

We inferred that the dispersal mechanism could be either wind and/or water dispersed due to the conflicting interpretations of our landscape hypotheses. This is the first study to use spatial modelling to infer dispersal in semi-aquatic lichens. The results of this study may help to understand lichen dispersal within aquatic landscapes, which can have implications in the conservation of rare or threatened lichens.

The fruticose genera in the Ramalinaceae (Ascomycota, Lecanoromycetes): their diversity and evolutionary history

We present phylogenetic analyses of the fruticose Ramalinaceae based on extensive collections from many parts of the world, with a special focus on the Vizcaíno deserts in north-western Mexico and the coastal desert in Namibia. We generate a four-locus DNA sequence dataset for accessions of Ramalina and two additional loci for Niebla and Vermilacinia. Four genera are strongly supported: the subcosmopolitan Ramalina, the new genus Namibialina endemic to SW Africa, and a duo formed by Niebla and Vermilacinia, endemic to the New World except the sorediate V. zebrina that disjunctly occurs in Namibia. The latter three genera are restricted to coastal desert and chaparral where vegetation depends on moisture from ocean fog. Ramalina is subcosmopolitan and much more diverse in its ecology. We show that Ramalina and its sister genus Namibialina diverged from each other at c. 48 Myrs, whereas Vermilacinia and Niebla split at c. 30 Myrs. The phylogeny of the fruticose genera remains unresolved to their ancestral crustose genera. Species delimitation within Namibialina and Ramalina is rather straightforward. The phylogeny and taxonomy of Vermilacinia are fully resolved, except for the two youngest clades of corticolous taxa, and support current taxonomy, including four new taxa described here. Secondary metabolite variation in Niebla generally coincides with major clades which are comprised of species complexes with still unresolved phylogenetic relationships. A micro-endemism pattern of allopatric species is strongly suspected for both genera, except for the corticolous taxa within Vermilacinia. Both Niebla and saxicolous Vermilacinia have chemotypes unique to species clades that are largely endemic to the Vizcaíno deserts. The following new taxa are described: Namibialina gen. nov. with N. melanothrix (comb. nov.) as type species, a single new species of Ramalina (R. krogiae) and four new species of Vermilacinia (V. breviloba, V. lacunosa, V. pustulata and V. reticulata). The new combination V. granulans is introduced. Two epithets are re-introduced for European Ramalina species: R. crispans (= R. peruviana auct. eur.) and R. rosacea (= R. bourgeana auct. p.p). A lectotype is designated for Vermilacinia procera. A key to saxicolous species of Vermilacinia is presented.

Spatio-temporal formation of the genetic diversity in the Mediterranean dwelling lichen during the Neogene and Quaternary epochs

Genetic patterns of lichenized fungi often display a mosaic-like and difficult to interpret structure blurring their evolutionary history. The genetic diversity and phylogeographic pattern of a mycobiont of the predominantly Mediterranean dwelling lichen Solenopsora candicans were investigated on the base of extensive sampling (361 individuals, 77 populations) across its entire distribution range. We tested whether the genetic pattern of S. candicans mirrors paleoclimatic and paleogeological events in the Mediterranean and adjacent regions. The divergence time estimates indicated a Tertiary origin for S. candicans, with formation of intraspecific diversity initiated in the Late Miocene. The distribution of the most divergent haplotypes, mostly of a pre-Pleistocene origin, was restricted to the eastern or western extremities of the Mediterranean exhibiting Kiermack disjunction. The population genetic diversity analyses indicated multiple diversity centres and refugia for S. candicans across the entire Mediterranean Basin. While the south Mediterranean regions harboured both the Tertiary and Quaternary born diversity, conforming to the 'cumulative refugia' paradigm, the Apennine and Balkan Peninsulas in the north hosted mostly younger Pleistocene haplotypes and lineages. The recent population expansion of S. candicans might have occurred in the middle Pleistocene with a population burst in the Apennine and Balkan peninsulas. The presence of unique haplotypes in Central Europe indicates the existence of extra-Mediterranean microrefugia. This study presents the first comprehensive lichen phylogeography from the Mediterranean region and simultaneously reports for the first time the glacial survival of a warm-adapted lichen in the temperate zone.

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.

Evolutionary lessons from California plant phylogeography

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

A worldwide phylogeography of the whiteworm lichens Thamnolia reveals three lineages with distinct habitats and evolutionary histories

Thamnolia is a lichenized fungus with an extremely wide distribution, being encountered in arctic and alpine environments in most continents. In this study, we used molecular markers to investigate the population structure of the fungal symbiont and the associated photosynthetic partner of Thamnolia. By analyzing molecular, morphological, and chemical variation among 253 specimens covering the species distribution range, we revealed the existence of three mycobiont lineages. One lineage (Lineage A) is confined to the tundra region of Siberia and the Aleutian Islands, a second (Lineage B) is found in the high alpine region of the Alps and the Carpathians Mountains, and a third (Lineage C) has a worldwide distribution and covers both the aforementioned ecosystems. Molecular dating analysis indicated that the split of the three lineages is older than the last glacial maximum, but the distribution ranges and the population genetic analyses suggest an influence of last glacial period on the present‐day population structure of each lineage. We found a very low diversity of Lineage B, but a higher and similar one in Lineages A and C. Demographic analyses suggested that Lineage C has its origin in the Northern Hemisphere, possibly Scandinavia, and that it has passed through a bottleneck followed by a recent population expansion. While all three lineages reproduce clonally, recombination tests suggest rare or past recombination in both Lineages A and C. Moreover, our data showed that Lineage C has a comparatively low photobiont specificity, being found associated with four widespread Trebouxia lineages (three of them also shared with other lichens), while Lineages A and B exclusively harbor T. simplex s. lat. Finally, we did not find support for the recognition of taxa in Thamnolia based on either morphological or chemical characters.

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

Background

Many fungal species occur across a variety of habitats. Particularly lichens, fungi forming symbioses with photosynthetic partners, have evolved remarkable tolerances for environmental extremes. Despite their ecological importance and ubiquity, little is known about the genetic basis of adaption in lichen populations. Here we studied patterns of genome-wide differentiation in the lichen-forming fungus Lasallia pustulata along an altitudinal gradient in the Mediterranean region. We resequenced six populations as pools and identified highly differentiated genomic regions. We then detected gene-environment correlations while controlling for shared population history and pooled sequencing bias, and performed ecophysiological experiments to assess fitness differences of individuals from different environments.

Results

We detected two strongly differentiated genetic clusters linked to Mediterranean and temperate-oceanic climate, and an admixture zone, which coincided with the transition between the two bioclimates. High altitude individuals showed ecophysiological adaptations to wetter and more shaded conditions. Highly differentiated genome regions contained a number of genes associated with stress response, local environmental adaptation, and sexual reproduction.

Conclusions

Taken together our results provide evidence for a complex interplay between demographic history and spatially varying selection acting on a number of key biological processes, suggesting a scenario of ecological speciation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0929-8) contains supplementary material, which is available to authorized users.

Fungal and cyanobacterial gene expression in a lichen symbiosis: Effect of temperature and location

Organisms have evolved different cellular mechanisms to deal with environmental stress, primarily through complex molecular mechanisms including protein refolding and DNA repair. As mutualistic symbioses, lichens offer the possibility of analyzing molecular stress responses in a particularly tight interspecific relationship. We study the widespread cyanolichen Peltigera membranacea, a key player in carbon and nitrogen cycling in terrestrial ecosystems at northern latitudes. We ask whether increased temperature is reflected in mRNA levels of selected damage control genes, and do the response patterns show geographical associations? Using real-time PCR quantification of 38 transcripts, differential expression was demonstrated for nine cyanobacterial and nine fungal stress response genes (plus the fungal symbiosis-related lec2 gene) when the temperature was increased from 5 °C to 15 °C and 25 °C. Principle component analysis (PCA) revealed two gene groups with different response patterns. Whereas a set of cyanobacterial DNA repair genes and the fungal lec2 (PC1 group) showed an expression drop at 15 °C vs. 5 °C, most fungal candidates (PC2 group) showed increased expression at 25 °C vs. 5 °C. PC1 responses also correlated with elevation. The correlated downregulation of lec2 and cyanobacterial DNA repair genes suggests a possible interplay between the symbionts warranting further studies.

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

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

Evolutionary and demographic history of the Californian scrub white oak species complex: an integrative approach

Understanding the factors promoting species formation is a major task in evolutionary research. Here, we employ an integrative approach to study the evolutionary history of the Californian scrub white oak species complex (genus Quercus). To infer the relative importance of geographical isolation and ecological divergence in driving the speciation process, we (i) analysed inter- and intraspecific patterns of genetic differentiation and employed an approximate Bayesian computation (ABC) framework to evaluate different plausible scenarios of species divergence. In a second step, we (ii) linked the inferred divergence pathways with current and past species distribution models (SDMs) and (iii) tested for niche differentiation and phylogenetic niche conservatism across taxa. ABC analyses showed that the most plausible scenario is the one considering the divergence of two main lineages followed by a more recent pulse of speciation. Genotypic data in conjunction with SDMs and niche differentiation analyses support that different factors (geography vs. environment) and modes of speciation (parapatry, allopatry and maybe sympatry) have played a role in the divergence process within this complex. We found no significant relationship between genetic differentiation and niche overlap, which probably reflects niche lability and/or that multiple factors, have contributed to speciation. Our study shows that different mechanisms can drive divergence even among closely related taxa representing early stages of species formation and exemplifies the importance of adopting integrative approaches to get a better understanding of the speciation process.

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

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

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

• Premise of the study: Many lichens exhibit extensive ranges spanning several ecoregions. It has been hypothesized that this wide ecological amplitude is facilitated by fungal association with locally adapted photobiont strains.• Methods: We studied the identity and geographic distribution of photobionts of the widely distributed North American lichen Ramalina menziesii based on rbcL (chloroplast DNA) and nuclear ribosomal ITS DNA sequences. To test for ecological specialization, we associate photobiont genotypes with local climate and phorophyte.• Key results: Of the photobiont lineages of R. menziesii, 94% belong to a clade including Trebouxia decolorans. The remaining are related to T. jamesii. The photobionts showed (1) significant structure according to ecoregion and phorophyte species and (2) genetic associations with phorophyte species and climate.• Conclusions: Geography, climate, and ecological specialization shape genetic differentiation of lichen photobionts. One great advantage of independent dispersal of the fungus is symbiotic association with locally adapted photobiont strains.