Sequencing genomes from mixed DNA samples - evaluating the metagenome skimming approach in lichenized fungi.

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.

Circumscription of the genus Lepra, a recently resurrected genus to accommodate the "Variolaria"-group of Pertusaria sensu lato (Pertusariales, Ascomycota).

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.

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.

Potential and pitfalls of eukaryotic metagenome skimming: a test case for lichens.

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.

Meta-analysis of deep-sequenced fungal communities indicates limited taxon sharing between studies and the presence of biogeographic patterns.

High-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.