A diverse population of introns in the nuclear ribosomal genes of ericoid mycorrhizal fungi includes elements with sequence similarity to endonuclease-coding genes

Does One Size Fit All? Variations in the DNA Barcode Gaps of Macrofungal Genera

The nuclear ribosomal internal transcribed spacer (nrITS) region has been widely used in fungal diversity studies. Environmental metabarcoding has increased the importance of the fungal DNA barcode in documenting fungal diversity and distribution. The DNA barcode gap is seen as the difference between intra- and inter-specific pairwise distances in a DNA barcode. The current understanding of the barcode gap in macrofungi is limited, inhibiting the development of best practices in applying the nrITS region toward research on fungal diversity. This study examined the barcode gap using 5146 sequences representing 717 species of macrofungi from eleven genera, eight orders and two phyla in datasets assembled by taxonomic experts. Intra- and inter-specific pairwise distances were measured from sequence and phylogenetic data. The results demonstrate that barcode gaps are influenced by differences in intra- and inter-specific variance in pairwise distances. In terms of DNA barcode behavior, variance is greater in the ITS1 than ITS2, and variance is greater in both relative to the combined nrITS region. Due to the difference in variance, the barcode gaps in the ITS2 region are greater than in the ITS1. Additionally, the taxonomic approach of "splitting" taxa into numerous taxonomic units produces greater barcode gaps when compared to "lumping". The results show variability in the barcode gaps between fungal taxa, demonstrating a need to understand the accuracy of DNA barcoding in quantifying species richness. For taxonomic studies, variability in nrITS sequence data supports the application of multiple molecular markers to corroborate the taxonomic and systematic delineation of species.

Porcine fungal mock community analyses: Implications for mycobiome investigations

Introduction

The gut microbiome is an integral partner in host health and plays a role in immune development, altered nutrition, and pathogen prevention. The mycobiome (fungal microbiome) is considered part of the rare biosphere but is still a critical component in health. Next generation sequencing has improved our understanding of fungi in the gut, but methodological challenges remain. Biases are introduced during DNA isolation, primer design and choice, polymerase selection, sequencing platform selection, and data analyses, as fungal reference databases are often incomplete or contain erroneous sequences.

Methods

Here, we compared the accuracy of taxonomic identifications and abundances from mycobiome analyses which vary among three commonly selected target gene regions (18S, ITS1, or ITS2) and the reference database (UNITE - ITS1, ITS2 and SILVA - 18S). We analyze multiple communities including individual fungal isolates, a mixed mock community created from five common fungal isolates found in weanling piglet feces, a purchased commercial fungal mock community, and piglet fecal samples. In addition, we calculated gene copy numbers for the 18S, ITS1, and ITS2 regions of each of the five isolates from the piglet fecal mock community to determine whether copy number affects abundance estimates. Finally, we determined the abundance of taxa from several iterations of our in-house fecal community to assess the effects of community composition on taxon abundance.

Results

Overall, no marker-database combination consistently outperformed the others. Internal transcribed space markers were slightly superior to 18S in the identification of species in tested communities, but Lichtheimia corymbifera, a common member of piglet gut communities, was not amplified by ITS1 and ITS2 primers. Thus, ITS based abundance estimates of taxa in piglet mock communities were skewed while 18S marker profiles were more accurate. Kazachstania slooffiae displayed the most stable copy numbers (83-85) while L. corymbifera displayed significant variability (90-144) across gene regions.

Discussion

This study underscores the importance of preliminary studies to assess primer combinations and database choice for the mycobiome sample of interest and raises questions regarding the validity of fungal abundance estimates.

Analysis of small and large subunit rDNA introns from several ectomycorrhizal fungi species

The small (18S) and large (28S) nuclear ribosomal DNA (rDNA) introns have been researched and sequenced in a variety of ectomycorrhizal fungal taxa in this study, it is found that both 18S and 28S rDNA would contain introns and display some degree variation in size, nucleotide sequences and insertion positions within the same fungi species (Meliniomyces). Under investigations among the tested isolates, 18S rDNA has four sites for intron insertions, 28S rDNA has two sites for intron insertions. Both 18S and 28S rDNA introns among the tested isolates belong to group I introns with a set of secondary structure elements designated P1-P10 helics and loops. We found a 12 nt nucleotide sequences TACCACAGGGAT at site 2 in the 3'-end of 28S rDNA, site 2 introns just insert the upstream or the downstream of the12 nt nucleotide sequences. Afters sequence analysis of all 18S and 28S rDNA introns from tested isolates, three high conserved regions around 30 nt nucleotides (conserved 1, conserved 2, conserved 3) and identical nucleotides can be found. Conserved 1, conserved 2 and conserved 3 regions have high GC content, GC percentage is almost more than 60%. From our results, it seems that the more convenient host sites, intron sequences and secondary structures, or isolates for 18S and 28S rDNA intron insertion and deletion, the more popular they are. No matter 18S rDNA introns or 18S rDNA introns among tested isolates, complementary base pairing at the splicing sites in P1-IGS-P10 tertiary helix around 5'-end introns and exons were weak.

Ericoid mycorrhizal symbiosis: theoretical background and methods for its comprehensive investigation

Despite decades of intensive research (especially from 1970s to 1990s), the ericoid mycorrhizal (ErM) hair root is still largely terra incognita and this simplified guide is intended to revive and promote the study of its mycobiota. Basic theoretical knowledge on the ErM symbiosis is summarized, followed by practical advices on Ericaceae root sample collection and handling, microscopic observations and photo-documentation of root fungal colonization, mycobiont isolation, maintenance and identification and resynthesis experiments with ericoid plants. The necessity of a proper selection of the root material and its surface sterilization prior to mycobiont isolation is stressed, together with the need of including suitable control treatments in inoculation experiments. The culture-dependent approach employing plating of single short (~ 2 mm) hair root segments on nutrient media is substantiated as a useful tool for characterization of Ericaceae root-associated fungal communities; it targets living mycelium and provides metabolically active cultures that can be used in physiological experiments and taxonomic studies, thus providing essential reference material for culture-independent approaches. On the other hand, it is stressed that not every mycobiont isolated from an ericoid hair root necessarily represent an ErM fungus. Likewise, not every intracellular hyphal coil formed in the Ericaceae rhizodermis necessarily represents the ErM symbiosis. Taxonomy of the most important ericoid mycobionts is updated, mutualism in the ErM symbiosis is briefly discussed from the mycobiont perspective, and some interesting lines of possible future research are highlighted.

Comparison of the performance of ITS1 and ITS2 as barcodes in amplicon-based sequencing of bioaerosols

This paper presents the performance of two eukaryotic genomic ribosomal regions, ITS1 and ITS2, in describing fungal diversity in aerosol samples using amplicon-based High-Throughput Sequencing (HTS). Composting sites, biomethanization facilities, and dairy farms, all affected by the presence of fungi, were visited to collect air samples. The amplicon-based HTS approach is a target enrichment method that relies on the amplification of a specific target using particular primers before sequencing. Thus, the results are highly dependent on the quality of amplification. For this reason, the authors of this paper used a shotgun metagenomic approach to compare its outcome with the amplicon-based method. Indeed, shotgun metagenomic does not rely on any amplification prior to sequencing, because all genes are sequenced without a specific target. In addition, culture methods were added to the analyses in biomethanization and dairy farms samples to validate their contribution to fungal diversity of aerosols. The results obtained are unequivocal towards ITS1 outperformance to ITS2 in terms of richness, and taxonomic coverage. The differential abundance analysis did demonstrate that some taxa were exclusively detected only by ITS2, and vice-versa for ITS1. However, the shotgun metagenomic approach showed a taxonomic profile more resembling to ITS1 than ITS2. Based on these results, neither of the barcodes evaluated is perfect in terms of distinguishing all species. Using both barcodes offers a broader view of the fungal aerosol population. However, with the actual knowledge, the authors strongly recommend using ITS1 as a universal fungal barcode for quick general analyses of diversity and when limited financial resources are available, primarily due its ability to capture taxonomic profiles similar to those obtained using the shotgun metagenomic. The culture comparison with amplicon-based sequencing showed the complementarity of both approaches in describing the most abundant taxa.

Integrated evolution of ribosomal RNAs, introns, and intron nurseries

The initial components of ribosomes first appeared more than 3.8 billion years ago during a time when many types of RNAs were evolving. While modern ribosomes are complex molecular machines consisting of rRNAs and proteins, they were assembled during early evolution by the association and joining of small functional RNA units. Introns may have provided the means to ligate many of these pieces together. All four classes of introns (group I, group II, spliceosomal, and archaeal) are present in many rRNA gene loci over a broad phylogenetic range. A survey of rRNA intron sequences across the three major life domains suggests that some of the classes of introns may have diverged from one another within rRNA gene loci. Analyses of rRNA sequences revealed self-splicing group I and group II introns are present in ancestral regions of the SSU (small subunit) and LSU (large subunit), whereas spliceosomal and archaeal introns appeared in sections of the rRNA that evolved later. Most classes of introns increased in number for approximately 1 billion years. However, their frequencies are low in the most recently evolved regions added to the SSU and LSU rRNAs. Furthermore, many of the introns appear to have been in the same locations for billions of years, suggesting an ancient origin for these sequences. In this Perspectives paper, I reviewed and analyzed rRNA intron sequences, locations, structural characteristics, and splicing mechanisms; and suggest that rRNA gene loci may have served as evolutionary nurseries for intron formation and diversification.

An amplicon-based sequencing approach for the study of aeromycology

Not available.

Accurate Estimation of Fungal Diversity and Abundance through Improved Lineage-Specific Primers Optimized for Illumina Amplicon Sequencing

While high-throughput sequencing methods are revolutionizing fungal ecology, recovering accurate estimates of species richness and abundance has proven elusive. We sought to design internal transcribed spacer (ITS) primers and an Illumina protocol that would maximize coverage of the kingdom Fungi while minimizing nontarget eukaryotes. We inspected alignments of the 5.8S and large subunit (LSU) ribosomal genes and evaluated potential primers using PrimerProspector. We tested the resulting primers using tiered-abundance mock communities and five previously characterized soil samples. We recovered operational taxonomic units (OTUs) belonging to all 8 members in both mock communities, despite DNA abundances spanning 3 orders of magnitude. The expected and observed read counts were strongly correlated (r = 0.94 to 0.97). However, several taxa were consistently over- or underrepresented, likely due to variation in rRNA gene copy numbers. The Illumina data resulted in clustering of soil samples identical to that obtained with Sanger sequence clone library data using different primers. Furthermore, the two methods produced distance matrices with a Mantel correlation of 0.92. Nonfungal sequences comprised less than 0.5% of the soil data set, with most attributable to vascular plants. Our results suggest that high-throughput methods can produce fairly accurate estimates of fungal abundances in complex communities. Further improvements might be achieved through corrections for rRNA copy number and utilization of standardized mock communities.

IMPORTANCE

Fungi play numerous important roles in the environment. Improvements in sequencing methods are providing revolutionary insights into fungal biodiversity, yet accurate estimates of the number of fungal species (i.e., richness) and their relative abundances in an environmental sample (e.g., soil, roots, water, etc.) remain difficult to obtain. We present improved methods for high-throughput Illumina sequencing of the species-diagnostic fungal ribosomal marker gene that improve the accuracy of richness and abundance estimates. The improvements include new PCR primers and library preparation, validation using a known mock community, and bioinformatic parameter tuning.

Evaluation of the authenticity of a highly novel environmental sequence from boreal forest soil using ribosomal RNA secondary structure modeling

The number of sequences from both formally described taxa and uncultured environmental DNA deposited in the International Nucleotide Sequence Databases has increased substantially over the last two decades. Although the majority of these sequences represent authentic gene copies, there is evidence of DNA artifacts in these databases as well. These include lab artifacts, such as PCR chimeras, and biological artifacts such as pseudogenes or other paralogous sequences. Sequences that fall in basal positions in phylogenetic trees and appear distant from known sequences are particularly suspect. Phylogenetic analyses suggest that a novel sequence type (NS1) found in two boreal forest soil clone libraries belongs to the fungal kingdom but does not fall unambiguously within any known phylum. We have evaluated this sequence type using an array of secondary-structure analyses. To our knowledge, such analyses have never been used on environmental ribosomal sequences. Ribosomal secondary structure was modeled for four rRNA loci (ITS1, 5.8S, ITS2, 5' LSU). These models were analyzed for the presence of conserved domains, conserved nucleotide motifs, and compensatory base changes. Minimal free energy (MFE) foldings and GC contents of sequences representing the major fungal clades, as well as NS1, were also compared. NS1 displays secondary rRNA structures consistent with other fungi and many, but not all, conserved nucleotide motifs found across eukaryotes. However, our analyses show that many other authentic sequences from basal fungi lack more of these conserved motifs than does NS1. Together our findings suggest that NS1 represents an authentic gene copy. The methods described here can be used on any rRNA-coding sequence, not just environmental fungal sequences. As new-generation sequencing methods that yield shorter sequences become more widely implemented, methods that evaluate sequence authenticity should also be more widely implemented. For fungi, the adjacent 5.8S and ITS2 loci should be prioritized. This region is not only suited to distinguishing between closely related species, but it is also more informative in terms of expected secondary structure.

Ericoid mycorrhizal root fungi and their multicopper oxidases from a temperate forest shrub

Ericoid mycorrhizal fungi (ERM) may specialize in capturing nutrients from their host's litter as a strategy for regulating nutrient cycles in terrestrial ecosystems. In spite of their potential significance, we know little about the structure of ERM fungal communities and the genetic basis of their saprotrophic traits (e.g., genes encoding extracellular enzymes). Rhododendron maximum is a model ERM understory shrub that influences the nutrient cycles of montane hardwood forests in the southern Appalachians (North Carolina, USA). We sampled ERM roots of R. maximum from organic and mineral soil horizons and identified root fungi by amplifying and sequencing internal transcribed spacer (ITS) ribosomal DNA (rDNA) collected from cultures and clones. We observed 71 fungal taxa on ERM roots, including known symbionts Rhizoscyphus ericae and Oidiodendron maius, putative symbionts from the Helotiales, Chaetothyriales, and Sebacinales, ectomycorrhizal symbionts, and saprotrophs. Supporting the idea that ERM fungi are adept saprotrophs, richness of root-fungi was greater in organic than in mineral soil horizons. To study the genetic diversity of oxidative enzymes that contribute to decomposition, we amplified and sequenced a portion of genes encoding multicopper oxidases (MCOs) from ERM ascomycetes. Most fungi possessed multiple copies of MCO sequences with strong similarities to known ferroxidases and laccases. Our findings indicate that R. maximum associates with a taxonomically and ecologically diverse fungal community. The study of MCO gene diversity and expression may be useful for understanding how ERM root fungi regulate the cycling of nutrients between the host plant and the soil environment.

Evolution of small putative group I introns in the SSU rRNA gene locus of Phialophora species

BACKGROUND

Group I introns (specifically subgroup IC1) are common in the nuclear ribosomal RNA genes of fungi. While most range in length from more than 200 to nearly 1800 nucleotides (nt) in length, several small putative (or degenerate) group I introns have been described that are between 56 and 81 nt. Although small, previously we demonstrated that the PaSSU intron in the rRNA small subunit gene of Phialophora americana isolate Wang 1046 is capable of in vitro splicing using a standard group I intron pathway, thus qualifying it as a functional ribozyme.

FINDINGS

Here, we describe eight short putative group I introns, ranging in length from 63 to 75 nt, in the rRNA small subunit genes of Phialophora isolates, a fungal genus that ranges from saprobic to pathogenic on plants and animals. All contain putative pairing regions P1, P7, and P10, as well as a pairing region formed between the middle of the intron and part of the 3' exon. The other pairing regions common in the core of standard group I introns are absent. However, parts of the 3' exon may aid in the stabilization of these small introns. Although the eight putative group I introns were from at least three species of Phialophora, phylogenetic analysis indicated that the eight are monophyletic. They are also monophyletic with the small introns of two lichen-forming fungi, Porpidia crustulata and Arthonia lapidicola.

CONCLUSIONS

The small putative group I introns in Phialophora have common features that may represent group I introns at their minima. They appear to have a single origin as indicated by their monophyly in phylogenetic analyses.

Occurrence and characteristics of group 1 introns found at three different positions within the 28S ribosomal RNA gene of the dematiaceous Phialophora verrucosa: phylogenetic and secondary structural implications

BACKGROUND

Group 1 introns (ribozymes) are among the most ancient and have the broadest phylogenetic distribution among the known self-splicing ribozymes. Fungi are known to be rich in rDNA group 1 introns. In the present study, five sequences of the 28S ribosomal RNA gene (rDNA) regions of pathogenic dematiaceous Phialophora verrucosa were analyzed using PCR by site-specific primers and were found to have three insertions, termed intron-F, G and H, at three positions of the gene. We investigated the distribution of group 1 introns in this fungus by surveying 34 strains of P. verrucosa and seven strains of Phialophora americana as the allied species.

RESULTS

Intron-F's (inserted at L798 position) were found in 88% of P. verrucosa strains, while intron-G's (inserted at L1921) at 12% and intron-H's (inserted at L2563) at 18%. There was some correlation between intron distribution and geographic location. In addition, we confirmed that the three kinds of introns are group 1 introns from results of BLAST search, alignment analysis and Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR). Prediction of secondary structures and phylogenetic analysis of intron sequences identified introns-F and G as belonging to subgroup IC1. In addition, intron-H was identified as IE.

CONCLUSION

The three intron insertions and their insertion position in the 28S rDNA allowed the characterization of the clinical and environmental isolates of P. verrucosa and P. americana into five genotypes. All subgroups of introns-F and G and intron-H were characterized and observed for the first time in both species.

Taxonomy of Nephroselmis viridis sp. nov. (Nephroselmidophyceae, Chlorophyta), a sister marine species to freshwater N. olivacea

The genus Nephroselmis (Nephroselmidophyceae), which had been placed in the Prasinophyceae, is one of the primitive green flagellates that are important to our understanding of the early evolution of green plants. We studied a new species of Nephroselmis isolated from Japan, Fiji and South Africa. This species has been known for a long time as undescribed species 'N. viridis.' N. viridis possesses some ultrastructural characters shared with only the freshwater type species N. olivacea, including a disc-like structure beneath the pyrenoid and bipolar spiny body scales with 1-5-8-5-1 spines. Molecular phylogenetic analysis based on 18S rDNA also supports a sister relationship between N. viridis and N. olivacea. However, N. viridis is distinguishable from N. olivacea by the shape of its starch sheath, its scales, its pigment composition and its habitat. In this paper, we designate the formal description of N. viridis sp. nov. We also describe variability in the 18S rDNA introns of various N. viridis strains. This detailed study of N. viridis provides some insights into the evolution of Nephroselmis.

A bioinformatics pipeline for sequence-based analyses of fungal biodiversity

The internal transcribed spacer (ITS) is the locus of choice with which to characterize fungal diversity in environmental samples. However, methods to analyze ITS datasets have lagged behind the capacity to generate large amounts of sequence information. Here, we describe our bioinformatics pipeline to process large fungal ITS sequence datasets, from raw chromatograms to a spreadsheet of operational taxonomic unit (OTU) abundances across samples. Steps include assembling of reads originating from one clone, identifying primer "barcodes" or "tags," trimming vectors and primers, marking low-quality base calls and removing low-quality sequences, orienting sequences, extracting the ITS region from longer amplicons, and grouping sequences into OTUs. We expect that the principles and tools presented here are relevant to datasets arising from ever-evolving new technologies.

Molecular characterization of the symbionts associated with marine nematodes of the genus Robbea

Marine nematodes that carry sulfur-oxidizing bacteria on their cuticle (Stilbonematinae, Desmodoridae) migrate between oxidized and reduced sand layers thereby supplying their symbionts with oxygen and sulfide. These symbionts, in turn, constitute the worms' major food source. Due to the accessibility, abundance and relative simplicity of this association, stilbonematids may be useful to understand symbiosis establishment. Nevertheless, only the symbiont of Laxus oneistus has been found to constitute one single phylotype within the Gammaproteobacteria. Here, we characterized the symbionts of three yet undescribed nematodes that were morphologically identified as members of the genus Robbea. They were collected at the island of Corsica, the Cayman Islands and the Belize Barrier Reef. The surface of these worms is covered by a single layer of morphologically undistinguishable bacteria. 18S rDNA-based phylogenetic analysis showed that all three species belong to the Stilbonematinae, although they do not form a distinct cluster within that subfamily. 16S rDNA-based analysis of the symbionts placed them interspersed in the cluster comprising the sulfur-oxidizing symbionts of L. oneistus and of marine gutless oligochaetes. Finally, the presence and phylogeny of the aprA gene indicated that the symbionts of all three nematodes can use reduced sulfur compounds as an energy source.

Molecular phylogeny suggests a single origin of insect symbiosis in the Pucciniomycetes with support for some relationships within the genus Septobasidium

In the Pucciniomycetes, a class of fungi that includes the plant pathogenic rust fungi, insect parasitism is restricted to a single family, the Septobasidiaceae. The Septobasidiaceae form a variety of symbioses with scale insects and have remained largely unstudied since the 1930s. Transitions between plant and animal parasitism and between mutualism and parasitism cannot be fully addressed in the Basidiomycota without a clear phylogenetic hypothesis for the Septobasidiales. Here, molecular phylogenetic methods were applied to understand the origin of scale insect parasitism, test the monophyly of the order Septobasidiales, and evaluate the infrageneric concepts in the largest genus of scale insect parasites, Septobasidium. DNA sequence data from rRNA genes were used to infer higher-level relationships within the Pucciniomycetes, and data from translation elongation factor 1-alpha (tef1) were added for phylogenetic inference within the Septobasidiaceae. Data from tef1 revealed different intron arrangements within Septobasidium, but the molecule did not provide much additional phylogenetically informative data. Likelihood-model-based phylogenetic analyses of 44 Pucciniomycotina taxa provided moderate support for a single origin of insect parasitism. Within the Septobasidiaceae, there was little or no support for a monophyletic Septobasidium, and well-resolved subclades of Septobasidium species contradict previous morphological delimitations of groups within the genus.

Variability of nuclear SSU-rDNA group introns within Septoria species: incongruence with host sequence phylogenies

We report structural features and distribution patterns of 26 different group I introns located at three distinct nucleotide positions in nuclear small subunit ribosomal DNA (SSU-rDNA) of 10 Septoria and 4 other anamorphic species related to the teleomorphic genus Mycosphaerella. Secondary structure and sequence characteristics assigned the introns to the common IC1 and IE groups. Intron distribution patterns and phylogenetic relationships strongly suggested that some horizontal transfer events have occurred among the closely related fungal species sampled. To test this hypothesis, we used a comparative approach of intron- and rDNA-based phylogenies through MP- and ML-based topology tests. Our results showed two statistically well-supported major incongruences between the intron and the equivalent internal transcribed spacer (ITS) tree comparisons made. Such absence of a co-evolutive history between group I introns and host sequences is discussed relatively to the intron structures, the mechanisms of intron movement, and the biology of the Mycosphaerella pathogenic fungi.

Soil fungal communities underneath willow canopies on a primary successional glacier forefront: rDNA sequence results can be affected by primer selection and chimeric data

Soil fungal communities underneath willow canopies that had established on the forefront of a receding glacier were analyzed by cloning the polymerase chain reaction (PCR)-amplified partial small subunit (18S) of the ribosomal (rRNA) genes. Congruence between two sets of fungus-specific primers targeting the same gene region was analyzed by comparisons of inferred neighbor-joining topologies. The importance of chimeric sequences was evaluated by Chimera Check (Ribosomal Database Project) and by data reanalyses after omission of potentially chimeric regions at the 5'- and 3'-ends of the cloned amplicons. Diverse communities of fungi representing Ascomycota, Basidiomycota, Chytridiomycota, and Zygomycota were detected. Ectomycorrhizal fungi comprised a major component in the early plant communities in primary successional ecosystems, as both primer sets frequently detected basidiomycetes (Russulaceae and Thelephoraceae) forming mycorrhizal symbioses. Various ascomycetes (Ophiostomatales, Pezizales, and Sordariales) of uncertain function dominated the clone libraries amplified from the willow canopy soil with one set of primers, whereas the clone libraries of the amplicons generated with the second primer set were dominated by basidiomycetes. Accordingly, primer bias is an important factor in fungal community analyses using DNA extracted from environmental samples. A large proportion (>30%) of the cloned sequences were concluded to be chimeric based on their changing positions in inferred phylogenies after omission of possibly chimeric data. Many chimeric sequences were positioned basal to existing classes of fungi, suggesting that PCR artifacts may cause frequent discovery of new, higher level taxa (order, class) in direct PCR analyses. Longer extension times during the PCR amplification and a smaller number of PCR cycles are necessary precautions to allow collection of reliable environmental sequence data.

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.

Widespread association between the ericoid mycorrhizal fungus Rhizoscyphus ericae and a leafy liverwort in the maritime and sub-Antarctic

A recent study identified a fungal isolate from the Antarctic leafy liverwort Cephaloziella varians as the ericoid mycorrhizal associate Rhizoscyphus ericae. However, nothing is known about the wider Antarctic distribution of R. ericae in C. varians, and inoculation experiments confirming the ability of the fungus to form coils in the liverwort are lacking. Using direct isolation and baiting with Vaccinium macrocarpon seedlings, fungi were isolated from C. varians sampled from eight sites across a 1875-km transect through sub- and maritime Antarctica. The ability of an isolate to form coils in aseptically grown C. varians was also tested. Fungi with 98-99% sequence identity to R. ericae internal transcribed spacer (ITS) region and partial large subunit ribosomal (r)DNA sequences were frequently isolated from C. varians at all sites sampled. The EF4/Fung5 primer set did not amplify small subunit rDNA from three of five R. ericae isolates, probably accounting for the reported absence of the fungus from C. varians in a previous study. Rhizoscyphus ericae was found to colonize aseptically-grown C. varians intracellularly, forming hyphal coils. This study shows that the association between R. ericae and C. varians is apparently widespread in Antarctica, and confirms that R. ericae is at least in part responsible for the formation of the coils observed in rhizoids of field-collected C. varians.

High diversity of endophytic fungi from the pharmaceutical plant, Heterosmilax japonica Kunth revealed by cultivation-independent approach

Heterosmilax japonica Kunth is well recognized for its diuretic effects in China. However, little is known about its endophytic fungi. In this study, microbial communities inhabiting the stems of H. japonica in spring and summer were investigated by light microscopy and cultivation-independent approaches, such as RFLP analysis and sequencing of rDNA ITS library. Molecular phylogenetic analysis showed that a broad spectrum of fungi, including Mycosphaerella, Phomopsis, Aureobasidium, Cladosporium, Glomerella, Botryosphaeria, Guignardia, is able to colonize the plants internally. Particularly, several rDNA sequences determined in this study like YJ4-61 are not specifically affiliated with any currently documented fungal sequences in the public database. Several sequence types, such as YJ4-9 and YJ4-70, are significantly similar to some uncultured environmental samples. Furthermore, our result also showed that the samples collected in spring harbored more abundant endophytic populations than that in summer, implying a seasonal fluctuation for the endophytes in H. japonica.

Phylogenetic analyses suggest reverse splicing spread of group I introns in fungal ribosomal DNA

BACKGROUND

Group I introns have spread into over 90 different sites in nuclear ribosomal DNA (rDNA) with greater than 1700 introns reported in these genes. These ribozymes generally spread through endonuclease-mediated intron homing. Another putative pathway is reverse splicing whereby a free group I intron inserts into a homologous or heterologous RNA through complementary base-pairing between the intron and exon RNA. Reverse-transcription of the RNA followed by general recombination results in intron spread. Here we used phylogenetics to test for reverse splicing spread in a taxonomically broadly sampled data set of fungal group I introns including 9 putatively ancient group I introns in the rDNA of the yeast-like symbiont Symbiotaphrina buchneri.

RESULTS

Our analyses reveal a complex evolutionary history of the fungal introns with many cases of vertical inheritance (putatively for the 9 introns in S. buchneri) and intron lateral transfer. There are several examples in which introns, many of which are still present in S. buchneri, may have spread through reverse splicing into heterologous rDNA sites. If the S. buchneri introns are ancient as we postulate, then group I intron loss was widespread in fungal rDNA evolution.

CONCLUSION

On the basis of these results, we suggest that the extensive distribution of fungal group I introns is at least partially explained by the reverse splicing movement of existing introns into ectopic rDNA sites.

Divergent histories of rDNA group I introns in the lichen family Physciaceae

The wide but sporadic distribution of group I introns in protists, plants, and fungi, as well as in eubacteria, likely resulted from extensive lateral transfer followed by differential loss. The extent of horizontal transfer of group I introns can potentially be determined by examining closely related species or genera. We used a phylogenetic approach with a large data set (including 62 novel large subunit [LSU] rRNA group I introns) to study intron movement within the monophyletic lichen family Physciaceae. Our results show five cases of horizontal transfer into homologous sites between species but do not support transposition into ectopic sites. This is in contrast to previous work with Physciaceae small subunit (SSU) rDNA group I introns where strong support was found for multiple ectopic transpositions. This difference in the apparent number of ectopic intron movements between SSU and LSU rDNA genes may in part be explained by a larger number of positions in the SSU rRNA, which can support the insertion and/or retention of group I introns. In contrast, we suggest that the LSU rRNA may have fewer acceptable positions and therefore intron spread is limited in this gene.

The spread of LAGLIDADG homing endonuclease genes in rDNA

Group I introns that encode homing endonuclease genes (HEGs) are highly invasive genetic elements. Their movement into a homologous position in an intron-less allele is termed homing. Although the mechanism of homing is well understood, the evolutionary relationship between HEGs and their intron partners remains unclear. Here we have focused on the largest family of HEGs (encoding the protein motif, LAGLIDADG) to understand how HEGs and introns move in rDNA. Our analysis shows the phylogenetic clustering of HEGs that encode a single copy of the LAGLIDADG motif in neighboring, but often evolutionarily distantly related, group I introns. These endonucleases appear to have inserted into existing introns independent of ribozymes. In contrast, our data support a common evolutionary history for a large family of heterologous introns that encode HEGs with a duplicated LAGLIDADG motif. This finding suggests that intron/double-motif HEG elements can move into heterologous sites as a unit. Our data also suggest that a subset of the double-motif HEGs in rDNA originated from the duplication and fusion of a single-motif HEG encoded by present-day ribozymes in LSU rDNA.

Nuclear large subunit rDNA group I intron distribution in a population of Beauveria bassiana strains: phylogenetic implications

Four group I introns, designated Bb1, Bb2, Bb3 and Bb4, were identified in the entomopathogenic hyphomycete Beauveria bassiana. Sequence analyses of these introns verified that they were invariably inserted at specific target sequences after conserved positions Ec2563, Ec2449, Ec2066 and Ec1921 of the large nuclear subunit (LSU) rDNA 3'-end. Secondary structure modelling confirmed that Bb1 and Bb3 belonged to subgroup IE while Bb2 and Bb4 belonged to subgroup IC1. Intron presence, distribution and size-variation were studied in a population of 125 B. bassiana strains using site-specific primers. Nucleotide sequences and secondary structures were compared and showed considerable variations usually at P1, P6 and P9 helices, but concomitantly, high homology between members of the same site-specific group. Intron distribution studies revealed that few (7.2%) strains were intron-less, most contained one (28%), two (48%) or three (16%) introns, while only one strain contained all four introns. Bb4-like introns (Ec1921) were the most abundant (86.4%), whereas the other three introns were evenly represented (ca 30%) in the B. bassiana population. Analysis of intron genotype distribution indicated a tenuous dependence upon geographic origin or insect host species. Phylogenetic analysis of all B. bassiana LSU introns and their close relatives from other entomopathogenic fungi showed a strong correlation between specific insertion sites and intron subgroups, fully supported by corresponding clades, suggesting common ancestry of the site specific LSU introns.

Structural features and evolutionary considerations of group IB introns in SSU rDNA of the lichen fungus Teloschistes

Different species of the lichen-forming ascomycete fungus Teloschistes were found to contain group IB introns at position S1506 in the small subunit ribosomal RNA gene. We have characterized the structural organization and phylogeny of the Teloschistes introns Tco.S1506, Tla.S1506, and Tvi.S1506. Common features to all the introns are a small size, a compact RNA structure, and an atypical catalytic ribozyme core sequence motif. Variations in intron sizes, due to sequence extensions in the P1 and P8 loop segments, were observed in different species and isolates. Phylogenetic analyses based on the ITS1-5.8S-ITS2 region as well as the introns show that the Teloschistes S1506 introns represent a distinct evolutionary isolated cluster among the nuclear group I introns. Furthermore, introns from different lineages of Teloschistes villosus appear not strictly vertically inherited probably due to horizontal transfer in one of the lineages.

The evolution of homing endonuclease genes and group I introns in nuclear rDNA

Group I introns are autonomous genetic elements that can catalyze their own excision from pre-RNA. Understanding how group I introns move in nuclear ribosomal (r)DNA remains an important question in evolutionary biology. Two models are invoked to explain group I intron movement. The first is termed homing and results from the action of an intron-encoded homing endonuclease that recognizes and cleaves an intronless allele at or near the intron insertion site. Alternatively, introns can be inserted into RNA through reverse splicing. Here, we present the sequences of two large group I introns from fungal nuclear rDNA, which both encode putative full-length homing endonuclease genes (HEGs). Five remnant HEGs in different fungal species are also reported. This brings the total number of known nuclear HEGs from 15 to 22. We determined the phylogeny of all known nuclear HEGs and their associated introns. We found evidence for intron-independent HEG invasion into both homologous and heterologous introns in often distantly related lineages, as well as the "switching" of HEGs between different intron peripheral loops and between sense and antisense strands of intron DNA. These results suggest that nuclear HEGs are frequently mobilized. HEG invasion appears, however, to be limited to existing introns in the same or neighboring sites. To study the intron-HEG relationship in more detail, the S943 group I intron in fungal small-subunit rDNA was used as a model system. The S943 HEG is shown to be widely distributed as functional, inactivated, or remnant ORFs in S943 introns.

Soil fungal community assembly in a primary successional glacier forefront ecosystem as inferred from rDNA sequence analyses

•  Soil fungal community assembly in nonvegetated areas on the forefront of a receding glacier was analyzed by cloning of the PCR-amplified partial small subunit (18S) of the ribosomal RNA genes (rDNA) from soil DNA samples. •  Fungal sequences obtained from areas adjacent to the present glacier terminus (young substrate) represented three fungal phyla, whereas those obtained adjacent to the terminal moraine (old substrate) were distributed among Ascomycetes and Hymenomycetes. The cloned sequences from both substrates represented mainly filamentous ascomycetes or basidiomycetes with a likely affinity to Agaricales. Unexpected biotrophic fungi with affinities to Taphrinomycetes, Urediniomycetes (the rust fungi) and Ustilaginomycetes (the smut fungi) plus an unknown, likely chytridiomycetous group were detected exclusively in the young substrates. •  These observations of biotrophic fungi are attributed to an aerially deposited, dormant spore bank, which may also be present in the older substrate but is masked by larger active mycelial biomass. •  This study underlines the importance of stochastic events and airborne spore deposition in the assembly of early fungal communities.

Molecular diversity and phylogenetic affinities of symbiotic root-associated ascomycetes of the Helotiales in burnt and metal polluted habitats

•  The diversity and phylogenetic affinities of symbiotic root-associated ascomycetes of the Helotiales are reported here based on ITS1-5.8S-ITS2 (internal transcribed spacer, ITS) nrDNA sequences. •  Mycobionts were obtained from roots of ericoid plants and grasses and from Piceirhiza bicolorata ectomycorrhizas (pbECM) on conifers and hardwoods, predominantly in burnt and metal-polluted habitats. The mycobionts were sequenced through the ITS and compared with sequences of known helotialean taxa. •  We recognized 132 fungal ITS-sequences with affinity to the Helotiales, of which 75% (54 different ITS-genotypes) grouped within the Hymenoscyphus ericae aggregate including Phialophora finlandia. This aggregate showed stronger affinity to members of the Hyaloscyphaceae and Dermateaceae than to Hymenoscyphus fructigenus (genus-type species; Helotiaceae). Most of the pbECM mycobionts grouped with P. finlandia, although some grouped with H. ericae. Two genotypes co-occurred in ericoid and ectomycorrhizal roots. •  The H. ericae aggregate may be referable to a generic unit, and includes a diverse group of closely related, more or less darkly pigmented, root-associated ascomycetes where the borders between intra- and interspecific ITS-sequence variation, as well as different mycorrhizal and nonmycorrhizal root-symbioses, remains unclear.

Evolutionary dynamics of multiple group I introns in nuclear ribosomal RNA genes of endoparasitic fungi of the genus Cordyceps

A large number of group I introns were discovered in coding regions of small and large subunits of nuclear ribosomal RNA genes (SSU rDNA and LSU rDNA) in ascomycetous fungi of the genus CORDYCEPS: From 28 representatives of the genus, we identified in total 69 group I introns which were inserted at any of four specific sites in SSU rDNA and four specific sites in LSU rDNA. These group I introns reached sizes of up to 510 bp, occurred in up to eight sites in the same organism, and belonged to either subgroup IB3 or subgroup IC1 based on their sequence and structure. Introns inserted at the same site were closely related to each other among Cordyceps fungi, whereas introns inserted at different sites were phylogenetically distinct even in the same species. Mapped on the host phylogeny, the group I introns were generally not restricted to a particular lineage, but, rather, widely and sporadically distributed among distinct lineages. When the phylogenetic relationships of introns inserted at the same site were compared with the phylogeny of their hosts, the topologies were generally significantly congruent to each other. From these results, the evolutionary dynamics of multiple group I introns in Cordyceps fungi was inferred as follows: (1) most of the group I introns were already present at the eight sites in SSU and LSU rDNAs of the ancestor of the genus Cordyceps; (2) the introns have principally been immobile and vertically transmitted throughout speciation and diversification of Cordyceps fungi, which resulted in the phylogenetic congruence between the introns at the same site and their hosts; (3) in the course of vertical transmission, the introns have repeatedly been lost in a number of lineages independently, which has led to the present sporadic phylogenetic distribution of the introns; and (4) a few acquisitions of new introns, presumably through horizontal transmission, were identified in the evolutionary history of the genus Cordyceps, while no transpositions were detected. Losses of group I introns in SSU rDNA have occurred at least 27 times in the evolutionary course of the 28 Cordyceps members.

Role of Horizontal Gene Transfer in the Evolution of Fungi

Although evidence for horizontal gene transfer (HGT) in eukaryotes remains largely anecdotal, literature on HGT in fungi suggests that it may have been more important in the evolution of fungi than in other eukaryotes. Still, HGT in fungi has not been widely accepted because the mechanisms by which it may occur are unknown, because it is usually not directly observed but rather implied as an outcome, and because there are often equally plausible alternative explanations. Despite these reservations, HGT has been justifiably invoked for a variety of sequences including plasmids, introns, transposons, genes, gene clusters, and even whole chromosomes. In some instances HGT has also been confirmed under experimental conditions. It is this ability to address the phenomenon in an experimental setting that makes fungi well suited as model systems in which to study the mechanisms and consequences of HGT in eukaryotic organisms.