Heterogeneity of intron presence or absence in rDNA genes of the lichen species Physcia aipolia and P. stellaris

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.

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.

IMA Genome - F13: Draft genome sequences of Ambrosiella cleistominuta, Cercospora brassicicola, C. citrullina, Physcia stellaris, and Teratosphaeria pseudoeucalypti

Draft genomes of the fungal species Ambrosiella cleistominuta, Cercospora brassicicola, C. citrullina, Physcia stellaris, and Teratosphaeria pseudoeucalypti are presented. Physcia stellaris is an important lichen forming fungus and Ambrosiella cleistominuta is an ambrosia beetle symbiont. Cercospora brassicicola and C. citrullina are agriculturally relevant plant pathogens that cause leaf-spots in brassicaceous vegetables and cucurbits respectively. Teratosphaeria pseudoeucalypti causes severe leaf blight and defoliation of Eucalyptus trees. These genomes provide a valuable resource for understanding the molecular processes in these economically important fungi.

Characterizing the ribosomal tandem repeat and its utility as a DNA barcode in lichen-forming fungi

Background

Regions within the nuclear ribosomal operon are a major tool for inferring evolutionary relationships and investigating diversity in fungi. In spite of the prevalent use of ribosomal markers in fungal research, central features of nuclear ribosomal DNA (nrDNA) evolution are poorly characterized for fungi in general, including lichenized fungi. The internal transcribed spacer (ITS) region of the nrDNA has been adopted as the primary DNA barcode identification marker for fungi. However, little is known about intragenomic variation in the nrDNA in symbiotic fungi. In order to better understand evolution of nrDNA and the utility of the ITS region for barcode identification of lichen-forming fungal species, we generated nearly complete nuclear ribosomal operon sequences from nine species in the Rhizoplaca melanophthalma species complex using short reads from high-throughput sequencing.

Results

We estimated copy numbers for the nrDNA operon, ranging from nine to 48 copies for members of this complex, and found low levels of intragenomic variation in the standard barcode region (ITS). Monophyly of currently described species in this complex was supported in phylogenetic inferences based on the ITS, 28S, intergenic spacer region, and some intronic regions, independently; however, a phylogenetic inference based on the 18S provided much lower resolution. Phylogenetic analysis of concatenated ITS and intergenic spacer sequence data generated from 496 specimens collected worldwide revealed previously unrecognized lineages in the nrDNA phylogeny.

Conclusions

The results from our study support the general assumption that the ITS region of the nrDNA is an effective barcoding marker for fungi. For the R. melanophthalma group, the limited amount of potential intragenomic variability in the ITS region did not correspond to fixed diagnostic nucleotide position characters separating taxa within this species complex. Previously unrecognized lineages inferred from ITS sequence data may represent undescribed species-level lineages or reflect uncharacterized aspects of nrDNA evolution in the R. melanophthalma species complex.

Heterogeneity of intron presence/absence in Olifantiella sp. (Bacillariophyta) contributes to the understanding of intron loss

Although hypotheses have been proposed and developed to interpret the origins and functions of introns, substantial controversies remain about the mechanism of intron evolution. The availability of introns in the intermediate state is quite helpful for resolving this debate. In this study, a new strain of diatom (denominated as DB21-1) was isolated and identified as Olifantiella sp., which possesses multiple types of 18S rDNAs (obtained from genomic DNA; lengths ranged from 2,056 bp to 2,988 bp). Based on alignments between 18S rDNAs and 18S rRNA (obtained from cDNA; 1,783 bp), seven intron insertion sites (IISs) located in the 18S rDNA were identified, each of which displayed the polymorphism of intron presence/absence. Specific primers around each IIS were designed to amplify the introns and the results indicated that introns in the same IIS varied in lengths, while terminal sequences were conserved. Our study showed that the process of intron loss happens via a series of successive steps, and each step could derive corresponding introns under intermediate states. Moreover, these results indicate that the mechanism of genomic deletion that occurs at DNA level can also lead to exact intron loss.

"Cryptic" group-I introns in the nuclear SSU-rRNA gene of Verticillium dahliae

Group-I introns are widespread--though irregularly distributed--in eukaryotic organisms, and they have been extensively used for discrimination and phylogenetic analyses. Within the Verticillium genus, which comprises important phytopathogenic fungi, a group-I intron was previously identified in the SSU-rRNA (18S) gene of only V. longisporum. In this work, we aimed at elucidating the SSU-located intron distribution in V. dahliae and other Verticillium species, and the assessment of heterogeneity regarding intron content among rDNA repeats of fungal strains. Using conserved PCR primers for the amplification of the SSU gene, a structurally similar novel intron (sub-group IC1) was detected in only a few V. dahliae isolates. However, when intron-specific primers were used for the screening of a diverse collection of Verticillium isolates that originally failed to produce intron-containing SSU amplicons, most were found to contain one or both intron types, at variable rDNA repeat numbers. This marked heterogeneity was confirmed with qRT-PCR by testing rDNA copy numbers (varying from 39 to 70 copies per haploid genome) and intron copy ratios in selected isolates. Our results demonstrate that (a) IC1 group-I introns are not specific to V. longisporum within the Verticillium genus, (b) V. dahliae isolates of vegetative compatibility groups (VCGs) 4A and 6, which bear the novel intron at most of their rDNA repeats, are closely related, and (c) there is considerable intra-genomic heterogeneity for the presence or absence of introns among the ribosomal repeats. These findings underline that distributions of introns in the highly heterogeneous repetitive rDNA complex should always be verified with sensitive methods to avoid misleading conclusions for the phylogeny of fungi and other organisms.

Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota) and Pleistocene population expansions

Background

Factors promoting diversification in lichen symbioses remain largely unexplored. While Pleistocene events have been important for driving diversification and affecting distributions in many groups, recent estimates suggest that major radiations within some genera in the largest clade of macrolichens (Parmeliaceae, Ascomycota) vastly predate the Pleistocene. To better understand the temporal placement and sequence of diversification events in lichens, we estimated divergence times in a common lichen-forming fungal genus, Melanohalea, in the Northern Hemisphere. Divergence times were estimated using both concatenated gene tree and coalescent-based multilocus species tree approaches to assess the temporal context of major radiation events within Melanohalea. In order to complement our understanding of processes impacting genetic differentiation, we also evaluated the effects of Pleistocene glacial cycles on population demographics of distinct Melanohalea lineages, differing in reproductive strategies.

Results

We found that divergence estimates, from both concatenated gene tree and coalescent-based multilocus species tree approaches, suggest that diversification within Melanohalea occurred predominantly during the Miocene and Pliocene, although estimated of divergence times differed by up to 8.3 million years between the two methods. These results indicate that, in some cases, taxonomically diagnostic characters may be maintained among divergent lineages for millions of years. In other cases, similar phenotypic characters among non-sister taxa, including reproductive strategies, suggest the potential for convergent evolution due to similar selective pressures among distinct lineages. Our analyses provide evidence of population expansions predating the last glacial maximum in the sampled lineages. These results suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species, albeit with apparently different demographic histories between sexually and vegetatively reproducing lineages.

Conclusions

Our results contribute to the understanding of how major changes during the Miocene and Pliocene have been important in promoting diversification within common lichen-forming fungi in the northern Hemisphere. Additionally, we provide evidence that glacial oscillations have influenced current population structure of broadly distributed lichenized fungal species throughout the Holarctic.

Factors that affect large subunit ribosomal DNA amplicon sequencing studies of fungal communities: classification method, primer choice, and error

Nuclear large subunit ribosomal DNA is widely used in fungal phylogenetics and to an increasing extent also amplicon-based environmental sequencing. The relatively short reads produced by next-generation sequencing, however, makes primer choice and sequence error important variables for obtaining accurate taxonomic classifications. In this simulation study we tested the performance of three classification methods: 1) a similarity-based method (BLAST + Metagenomic Analyzer, MEGAN); 2) a composition-based method (Ribosomal Database Project naïve bayesian classifier, NBC); and, 3) a phylogeny-based method (Statistical Assignment Package, SAP). We also tested the effects of sequence length, primer choice, and sequence error on classification accuracy and perceived community composition. Using a leave-one-out cross validation approach, results for classifications to the genus rank were as follows: BLAST + MEGAN had the lowest error rate and was particularly robust to sequence error; SAP accuracy was highest when long LSU query sequences were classified; and, NBC runs significantly faster than the other tested methods. All methods performed poorly with the shortest 50-100 bp sequences. Increasing simulated sequence error reduced classification accuracy. Community shifts were detected due to sequence error and primer selection even though there was no change in the underlying community composition. Short read datasets from individual primers, as well as pooled datasets, appear to only approximate the true community composition. We hope this work informs investigators of some of the factors that affect the quality and interpretation of their environmental gene surveys.

DNA barcoding of lichenized fungi demonstrates high identification success in a floristic context

• Efforts are currently underway to establish a standard DNA barcode region for fungi; we tested the utility of the internal transcribed spacer (ITS) of nuclear ribosomal DNA for DNA barcoding in lichen-forming fungi by sampling diverse species across eight orders. • Amplification of the ITS region (ITS1-5.8S-ITS2) was conducted for 351 samples, encompassing 107, 55 and 28 species, genera and families, respectively, of lichenized fungi. We assessed the ability of the entire ITS vs the ITS2 alone to discriminate between species in a taxonomic dataset (members of the genus Usnea) and a floristic dataset. • In the floristic dataset, 96.3% of sequenced samples could be assigned to the correct species using ITS or ITS2; a barcode gap for ITS is present in 92.1% of species. Although fewer species have a barcode gap in the taxonomic dataset (73.3% with ITS and 68.8% with ITS2), up to 94.1% of samples were assigned to the correct species using BLAST. • While discrimination between the most closely related species will remain challenging, our results demonstrate the potential to identify a high percentage of specimens to the correct species, and the remainder to the correct genus, when using DNA barcoding in a floristic context.

Dissociation and horizontal transmission of codispersing lichen symbionts in the genus Lepraria (Lecanorales: Stereocaulaceae)

Lichenized fungi of the genus Lepraria lack ascomata and conidiomata, and symbionts codisperse by soredia. Here, it is determined whether algal symbionts associated with Lepraria are monophyletic, and whether fungal and algal phylogenies are congruent, both of which are indicative of a long-term, continuous association between symbionts. The internal transcribed spacer (ITS) and part of the actin type I locus were sequenced from algae associated with Lepraria, and the fungal ITS and mitochondrial small subunit (mtSSU) were sequenced from fungal symbionts. Phylogenetic analyses tested for monophyly of algal symbionts and congruence between algal and fungal phylogenies. Algae associated with Lepraria were not monophyletic, and identical algae associated with different Lepraria individuals and species. Algal and fungal phylogenies were not congruent, suggesting a lack of strict codiversification. This study suggests that associations between symbionts are not strictly maintained over evolutionary time. The ability to switch partners may provide benefits similar to genetic recombination, which may have helped this lineage persist.

Phylogenetic comparison of protein-coding versus ribosomal RNA-coding sequence data: A case study of the Lecanoromycetes (Ascomycota)

The resolving power and statistical support provided by two protein-coding (RPB1 and RPB2) and three ribosomal RNA-coding (nucSSU, nucLSU, and mitSSU) genes individually and in various combinations were investigated based on maximum likelihood bootstrap analyses on lichen-forming fungi from the class Lecanoromycetes (Ascomycota). Our results indicate that the optimal loci (single and combined) to use for molecular systematics of lichen-forming Ascomycota are protein-coding genes (RPB1 and RPB2). RPB1 and RPB2 genes individually were phylogenetically more efficient than all two- and three-locus combinations of ribosomal loci. The 3rd codon position of each of these two loci provided the most characters in support of phylogenetic relationships within the Lecanoromycetes. Of the three ribosomal loci we used in this study, mitSSU contributed the most to phylogenetic analyses when combined with RPB1 and RPB2. Except for the mitSSU, ribosomal genes were the most difficult to recover because they often contain many introns, resulting in PCR bias toward numerous and intronless co-extracted contaminant fungi (mainly Dothideomycetes, Chaetothyriomycetes, and Sordariomycetes in the Ascomycota, and members of the Basidiomycota), which inhabit lichen thalli. Maximum likelihood analysis on the combined five-locus data set for 82 members of the Lecanoromycetes provided a well resolved and well supported tree compared to existing phylogenies. We confirmed the monophyly of three recognized subclasses in the Lecanoromycetes, the Acarosporomycetidae, Ostropomycetidae, and Lecanoromycetideae; the latter delimited as monophyletic for the first time, with the exclusion of the family Umbilicariaceae and Hypocenomyce scalaris. The genus Candelariella (formerly in the Candelariaceae, currently a member of the Lecanoraceae) represents the first evolutionary split within the Lecanoromycetes, before the divergence of the Acarosporomycetidae. This study provides a foundation necessary to guide the selection of loci for future multilocus phylogenetic studies on lichen-forming and allied ascomycetes.

Patterns of group I intron presence in nuclear SSU rDNA of the Lichen family Parmeliaceae

Group I introns are commonly reported within nuclear SSU ribosomal DNA of eukaryotic micro-organisms, especially in lichen-forming fungi. We have studied the primary and secondary structure of 70 new nuclear SSU rDNA group I introns of Parmeliaceae (Ascomycota: Lecanorales) and compared them with those available in databases, covering more than 60 species. The analyzed samples of Parmeliaceae fell into two groups, one having an intron at the 1506 site and another lacking this one but having another at the 1516 or 1521 position. Introns at the 1521 position seem to be transposed from 1516 sites. Introns at the 1516 position were similar in structure to ones previously reported at this site and known from other lecanoralean fungi, while those at the 1506 position showed structural differences and no similar introns are known from related fungi. The study of the distribution of group I introns within a large monophyletic ensemble of fungi has revealed an unexpected correlation between intron types and ecological and geographical parameters. The introns at the 1516 position occurred in mainly arctic, boreal, and temperate lichens, while those at position 1506 were present in mainly tropical and subtropical to oceanic mild-temperate taxa. Further, the 1516 introns occurred in genera with few distributed species that could represent older taxa, while the 1506 ones were mainly in species-rich genera that could be of recent speciation, as many species have wide distribution areas. The transition between two different environments has been accompanied by a change in introns gained and lost.