BackgroundMany 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.ResultsWe 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.ConclusionsTaken 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 materialThe online version of this article (doi:10.1186/s12862-017-0929-8) contains supplementary material, which is available to authorized users.

SummaryHigh-throughput amplicon sequencing gives new insights into fungal community ecology. Massively generated molecular data lead to the discovery of vast fungal diversity. However, it is unclear to what extent operational taxonomic units (OTUs) overlap among independent studies, because no comparative studies exist.We compared fungal diversity based on the internal transcribed spacer (ITS1) region among 10 published studies. Starting from the raw 454 pyrosequencing data, we used a uniform pipeline to prune the reads. We investigated fungal richness and taxonomic composition among phyllosphere and soil fungal communities, as well as biogeographic signals in the data.We did not find globally distributed OTUs, even when comparing fungal communities from similar habitats (phyllosphere or soil). This suggests that high local fungal diversity scales up to high global diversity. The most OTU-rich classes in the phyllosphere were Dothideomycetes (21%) and Sordariomycetes (14%), and in the soil were Sordariomycetes (13%) and Agaricomycetes (12%). The richness estimates suggest the presence of undiscovered fungal diversity even in deeply sequenced study systems.The small number of OTUs shared among studies indicates that globally distributed taxa and habitat generalists may be rare. Latitudinal diversity decline and distance decay relationships suggest the presence of biogeographic patterns similar to those in plants and animals.

The metagenome skimming approach, i.e. low coverage shotgun sequencing of multi-species assemblages and subsequent reconstruction of individual genomes, is increasingly used for in-depth genomic characterization of ecological communities. This approach is a promising tool for reconstructing genomes of facultative symbionts, such as lichen-forming fungi, from metagenomic reads. However, no study has so far tested accuracy and completeness of assemblies based on metagenomic sequences compared to assemblies based on pure culture strains of lichenized fungi. Here we assembled the genomes of Evernia prunastri and Pseudevernia furfuracea based on metagenomic sequences derived from whole lichen thalli. We extracted fungal contigs using two different taxonomic binning methods, and performed gene prediction on the fungal contig subsets. We then assessed quality and completeness of the metagenome-based assemblies using genome assemblies as reference which are based on pure culture strains of the two fungal species. Our comparison showed that we were able to reconstruct fungal genomes from uncultured lichen thalli, and also cover most of the gene space (86–90%). Metagenome skimming will facilitate genome mining, comparative (phylo)genomics, and population genetics of lichen-forming fungi by circumventing the time-consuming, sometimes unfeasible, step of aposymbiotic cultivation.

Whole-genome shotgun sequencing of multispecies communities using only a single library layout is commonly used to assess taxonomic and functional diversity of microbial assemblages. Here, we investigate to what extent such metagenome skimming approaches are applicable for in-depth genomic characterizations of eukaryotic communities, for example lichens. We address how to best assemble a particular eukaryotic metagenome skimming data, what pitfalls can occur, and what genome quality can be expected from these data. To facilitate a project-specific benchmarking, we introduce the concept of twin sets, simulated data resembling the outcome of a particular metagenome sequencing study. We show that the quality of genome reconstructions depends essentially on assembler choice. Individual tools, including the metagenome assemblers Omega and MetaVelvet, are surprisingly sensitive to low and uneven coverages. In combination with the routine of assembly parameter choice to optimize the assembly N50 size, these tools can preclude an entire genome from the assembly. In contrast, MIRA, an all-purpose overlap assembler, and SPAdes, a multisized de Bruijn graph assembler, facilitate a comprehensive view on the individual genomes across a wide range of coverage ratios. Testing assemblers on a real-world metagenome skimming data from the lichen Lasallia pustulata demonstrates the applicability of twin sets for guiding method selection. Furthermore, it reveals that the assembly outcome for the photobiont Trebouxia sp. falls behind the a priori expectation given the simulations. Although the underlying reasons remain still unclear, this highlights that further studies on this organism require special attention during sequence data generation and downstream analysis.

Lasallia hispanica (Frey) Sancho & A. Crespo is one of three Lasallia species occurring in central-western Europe. It is an orophytic, photophilous Mediterranean endemic which is sympatric with the closely related, widely distributed, highly clonal sister taxon L. pustulata in the supra- and oro-Mediterranean belts. We sequenced the genome of L. hispanica from a multispore isolate. The total genome length is 41·2 Mb, including 8488 gene models. We present the annotation of a variety of genes that are involved in protein secretion, mating processes and secondary metabolism, and we report transposable elements. Additionally, we compared the genome of L. hispanica to the closely related, yet ecologically distant, L. pustulata and found high synteny in gene content and order. The newly assembled and annotated L. hispanica genome represents a useful resource for future investigations into niche differentiation, speciation and microevolution in L. hispanica and other members of the genus.

Pertusarialean lichens include more than 300 species belonging to several independent phylogenetic lineages. Only some of these phylogenetic clades have been comprehensively sampled for molecular data, and formally described as genera. Here we present a taxonomic treatment of a group of pertusarialean lichens formerly known as “Pertusaria amara-group”, “Monomurata-group”, or “Variolaria-group”, which includes widespread and well-known taxa such as P. amara, P. albescens, or P. ophthalmiza. We generated a 6-locus data set with 79 OTUs representing 75 species. The distinction of the Variolaria clade is supported and consequently, the resurrection of the genus Lepra is followed. Thirty-five new combinations into Lepra are proposed and the new species Lepra austropacifica is described from mangroves in the South Pacific. Lepra is circumscribed to include species with disciform ascomata, a weakly to non-amyloid hymenial gel, strongly amyloid asci without clear apical amyloid structures, containing 1 or 2, single-layered, thin-walled ascospores. Chlorinated xanthones are not present, but thamnolic and picrolichenic acids occur frequently, as well as orcinol depsides. Seventy-one species are accepted in the genus. Although the distinction of the genus from Pertusaria is strongly supported, the relationships of Lepra remain unresolved and the genus is tentatively placed in Pertusariales incertae sedis.

Abstract:We cultured Parmelina carporrhizans and P. quercina in Corn Meal Agar and 0.2% glucose Malt Yeast Agar for 160 days. Chemosyndrome of natural thalli and mycobiont cultures were analyzed by HPLC. Lecanoric acid, atranorin, chloratranorin and ergosterol were detected in P. carporrhizans thalli, while lecanoric acid, chloratranorin and aliphates were found in P. quercina thalli. The secondary methabolites pattern between thalli and mycobiont culture was completely different in both species. Both species secreted the phenalenone myeloconone C in culture media and was also detected in P. quercina mycobiont aggregates. Interestingly, the phenolic compounds produced by the mycobiont culture of P. carporrhizans are related to those produced by natural thallus by the same biosynthetic pathway, while the chemosyndrome of P. quercina mycobiont implies switch of biosynthetic pathway from acetate-polymalonate pathway to shikimic acid pathway, with pulvinic acid as major compound of mycobiont culture. The role of Myelochonone C, confluentic acid and pulvinic acid produced by mycobiont culture is discussed as possible adaptive vantage in field as photoprotective agent or as byproduct result of stressing artificial culture conditions.

No abstract