Evolution of rDNA ITS1 and ITS2 sequences and RNA secondary structures within members of the fungal genera Grosmannia and Leptographium

Isolation, molecular characterization and first report of Dothiorella gregaria associated with fruit rot of walnuts of Jammu and Kashmir, India

Walnuts are known for their high levels of antioxidants, which are linked to various health benefits. However, challenges related to distribution and storage, as well as the risk of fungal infections, can affect the quality of walnut kernels. Fungal pathogens from the Botryosphaeriaceae family, including Dothiorella species and Diplodia species, can damage fruit and reduce its antioxidant content. To comprehend the cause of fruit rot in walnuts, Dothiorella gregaria isolates were studied using polyphasic methods, including multiple gene sequences and morphological identification, as well as analysis of polyphenol content and pathogenicity. The walnuts kernels purchased from market places of Jammu and Kashmir (J&K), India were observed to be affected by Dothiorella gregaria species causing the quality detoriation and decrease in polyphenol content thus undeniably with decreased antioxidant properties. D. gregaria Infected walnut kernels were having some brown and black spots and some were having white mycelial growth and however, most samples were asymptomatic. Pathogenicity testing revealed that the pathogen was able to develop all the symptoms under experimental conditions and the reisolated pathogen was morphologically similar to D. gregaria. The samples infected with this pathogen showed considerable decrease in polyphenol content, 10.9 ± 2.66 mgGAE/g (mean ± standard deviation) thus decreased antioxidant quality as compared to the samples which showed zero incidence of this pathogen, 52.50 ± 4.27 mgGAE/g (mean ± standard deviation). Furthermore, the pathogen was studied using polyphasic approach involving morphological, molecular and phylogenetic analysis. Combined nucleotide dataset of nuclear ribosomal ITS and tef1-α revealed that Dothiorella gregaria (NY6) formed a clade with Dothiorlla iberica (MAEC33), Dothiorella sarmentorium (MAEC28) and Dothiorella iberica (CAA905) strains with 83% bootstrap support. Besides, we observed six nucleotide changes, four were insertions or deletions and two were substitutions in the 502-bp region of the ITS rRNA gene when we compared our isolate to the most equivalent sequences submitted to NCBI GenBank. This is the first report of Dothiorella gregaria affecting walnuts purchased from various markets in J&K, India, causing fruit rot in walnuts after harvest. Given that local farmers store and export walnuts, it could pose an emerging threat to their livelihood. Thus, creating post-harvesting interventions for D. gregaria and knowing more about the fruit rot in walnuts can be benefited from morphological and molecular identification using several gene loci, genetic variability in the ITS rRNA gene, and total phenol analysis.

Polymorphisms of rDNA genes in Cyberlindnera yeast suggest birth-and-death evolution events

In eukaryotes, the ribosome machinery is encoded by repeats of the ribosomal RNA genes: 26/28S, 18S, 5.8S, and 5S, structured in tandem arrays and frequently homogenized within a genome. This homogenization is thought to be driven by concerted evolution, evolving as a unit, which contributes to its target as the species barcode in modern taxonomy. However, high heterogeneity of rDNA genes has been reported, including in Saccharomycotina yeasts. Here, we describe the polymorphisms and heterogeneity of D1/D2 domains (26S rRNA) and the intergenic transcribed spacer of a new yeast species with affinities to the genus Cyberlindnera and their evolution. Both regions are not homogenized, failing the prediction of concerted evolution. Phylogenetic network analysis of cloned sequences revealed that Cyberlindnera sp. rDNAs are diverse and evolved by reticulation rather than by bifurcating tree evolution model. Predicted rRNA secondary structures also confirmed structural differences, except for some conserved hairpin loops. We hypothesize that some rDNA is inactive within this species and evolves by birth-and-death rather than concerted evolution. Our findings propel further investigation into the evolution of rDNA genes in yeasts.

Taxonomy and phylogeny of the Leptographium olivaceum complex (Ophiostomatales, Ascomycota), including descriptions of six new species from China and Europe

The Leptographiumolivacea complex encompasses species in the broadly defined genus Leptographium (Ophiostomatales, Ascomycota) that are generally characterized by synnematous conidiophores. Most species of the complex are associates of conifer-infesting bark beetles in Europe and North America. The aims of this study were to reconsider the delineation of known species, and to confirm the identity of several additional isolates resembling L.olivacea that have emerged from recent surveys in China, Finland, Poland, Russia, and Spain. Phylogenetic analyses of sequence data for five loci (ACT, TUB, CAL, ITS2-LSU, and TEF-1α) distinguished 14 species within the complex. These included eight known species (L.cucullatum, L.davidsonii, L.erubescens, L.francke-grosmanniae, L.olivaceum, L.olivaceapini, L.sagmatosporum, and L.vescum) and six new species (herein described as L.breviuscapum, L.conplurium, L.pseudoalbum, L.rhizoidum, L.sylvestris, and L.xiningense). New combinations are provided for L.cucullatum, L.davidsonii, L.erubescens, L.olivaceum, L.olivaceapini, L.sagmatosporum and L.vescum. New Typifications: Lectotypes are designated for L.olivaceum, L.erubescens and L.sagmatosporum. Epitypes were designated for L.olivaceapini and L.sagmatosporum. In addition to phylogenetic separation, the synnematous asexual states and ascomata with almost cylindrical necks and prominent ostiolar hyphae, distinguish the L.olivaceum complex from others in Leptographium.

Phylogenetic and morphological analyses of the mycoparasitic genus Piptocephalis

The Piptocephalidaceae (Zoopagales, Zoopagomycota) contains three genera of mycoparasitic, haustoria-forming fungi: Kuzuhaea, Piptocephalis, and Syncephalis. Although the species in this family are diverse and ubiquitous in soil and dung, they are among the least studied fungi. Co-cultures of Piptocephalis and their hosts are relatively easy to isolate from soil and dung samples across the globe, making them a good model taxon for the order Zoopagales. This study focuses on the systematics of the genus Piptocephalis. Despite the fact that there are approximately 40 described Piptocephalis species, there are no modern taxonomic or molecular phylogenetic treatments of this group. Minimal sequence data are available, and relatively little is known about the true diversity or biogeography of the genus. Our study addresses two aspects: Piptocephalis systematics and analyses of the length and inter- and infraspecific variation of the nuc rDNA internal transcribed spacer (ITS1-5.8S-ITS2 = ITS) region. First, we generated a large subunit (28S) nuc rDNA phylogeny and evaluated several morphological characters by testing their correlation with the phylogeny using Bayesian Tip-association Significance testing (BaTS). We found monophyly of Piptocephalis species identified based on morphological traits, but morphological character states were not conserved across clades, suggesting that there have been multiple gains and losses of morphological characters. We also found that Kuzhuaea is nested within Piptocephalis. Second, we amplified the ITS from many Piptocephalis isolates, created a sequence alignment, and measured the lengths using the software ITSx. Piptocephalis species had ITS regions that were longer than the average for most Dikarya but were similar in length to those of the related genus Syncephalis.

Evaluation of the ribosomal DNA internal transcribed spacer (ITS), specifically ITS1 and ITS2, for the analysis of fungal diversity by deep sequencing

The nuclear ribosomal DNA internal transcribed spacer (ITS) has been widely used to assess the fungal composition in different environments by deep sequencing. To evaluate the ITS in the analysis of fungal diversity, comparisons of the clustering and taxonomy generated by sequencing with different portions of the whole fragment were conducted in this study. For a total of 83,120 full-length ITS sequences obtained from the UNITE database, it was found that, on average, ITS1 varied more than ITS2 within the kingdom Fungi; this variation included length and GC content variations and polymorphisms, with some polymorphisms specific to particular fungal groups. The taxonomic accuracy for ITS was higher than that for ITS1 or ITS2. The commonly used operational taxonomic unit (OTU) for evaluating fungal diversity and richness assigned several species to a single OTU even with clustering at 99.00% sequence similarity. The clustering and taxonomic capacities did not differ between ITS1 and ITS2. However, the OTU commonality between ITS1 and ITS2 was very low. To test this observation further, 219,741 pyrosequencing reads, including 39,840 full-length ITS sequences, were obtained from 10 soil samples and were clustered into OTUs. The pyrosequencing results agreed with the results of the in silico analysis. ITS1 might overestimate the fungal diversity and richness. Analyses using ITS, ITS1 and ITS2 yielded several different taxa, and the taxonomic preferences for ITS and ITS2 were similar. The results demonstrated that ITS2 alone might be a more suitable marker for revealing the operational taxonomic richness and taxonomy specifics of fungal communities when the full-length ITS is not available.

Ophiostomatoid fungi associated with pines infected by Bursaphelenchus xylophilus and Monochamus alternatus in China, including three new species

The activity of the pine wood nematodeBursaphelenchusxylophilusleads to extremely serious economic, ecological and social losses in East Asia. The nematode causes pine wilt disease, which is currently regarded as the most important forest disease in China. The pathogenic nematode feeds on dendrocola fungi to complete its cycle of infection. As the vector of the nematode, the Japanese pine sawyer (Monochamusalternatus) also carries dendrocola fungi. Pine woods, infected byB.xylophilusand tunnelled byM.alternatus, are also inhabited by ophiostomatoid fungi. These fungi are well known for their association with many bark and ambrosia beetles. They can cause sapstain and other serious tree diseases. The aims of our study were to investigate and identify the ophiostomatoid communities associated with the epidemic pine wood nematode and the pine sawyer inPinusmassonianaandP.thunbergiiforests, which are the main hosts of the pine wood nematode in China. Two hundred and forty strains of ophiostomatoid fungi were isolated from nematode and sawyer–infected trees in the coastal Shandong and Zhejiang Provinces, representing newly and historically infected areas, respectively. Six ophiostomatoid species were identified on the basis of morphological, physiological and molecular data. For the latter, DNA sequences of the internal transcribed spacer (ITS1–5.8S–ITS2) region and partial b-tubulin gene were examined. The ophiostomatoid species included one known species,Ophiostomaips, three novel species, viz.Ophiostomaalbumsp. nov.,Ophiostomamassonianasp. nov.andSporothrixzhejiangensissp. nov.and two species whose identities are still uncertain, Ophiostomacf.deltoideosporum and Graphilbumcf.rectangulosporium, due to the paucity of the materials obtained. The ophiostomatoid community was dominated byO.ips. This study revealed that a relatively high species diversity of ophiostomatoid fungi are associated with pine infected byB.xylophilusandM.alternatusin China.

Ophiostomatoid fungi associated with pines infected by Bursaphelenchusxylophilus and Monochamusalternatus in China, including three new species

The activity of the pine wood nematode Bursaphelenchusxylophilus leads to extremely serious economic, ecological and social losses in East Asia. The nematode causes pine wilt disease, which is currently regarded as the most important forest disease in China. The pathogenic nematode feeds on dendrocola fungi to complete its cycle of infection. As the vector of the nematode, the Japanese pine sawyer (Monochamusalternatus) also carries dendrocola fungi. Pine woods, infected by B.xylophilus and tunnelled by M.alternatus, are also inhabited by ophiostomatoid fungi. These fungi are well known for their association with many bark and ambrosia beetles. They can cause sapstain and other serious tree diseases. The aims of our study were to investigate and identify the ophiostomatoid communities associated with the epidemic pine wood nematode and the pine sawyer in Pinusmassoniana and P.thunbergii forests, which are the main hosts of the pine wood nematode in China. Two hundred and forty strains of ophiostomatoid fungi were isolated from nematode and sawyer-infected trees in the coastal Shandong and Zhejiang Provinces, representing newly and historically infected areas, respectively. Six ophiostomatoid species were identified on the basis of morphological, physiological and molecular data. For the latter, DNA sequences of the internal transcribed spacer (ITS1-5.8S-ITS2) region and partial b-tubulin gene were examined. The ophiostomatoid species included one known species, Ophiostomaips, three novel species, viz. Ophiostomaalbum sp. nov., Ophiostomamassoniana sp. nov. and Sporothrixzhejiangensis sp. nov. and two species whose identities are still uncertain, Ophiostomacf.deltoideosporum and Graphilbumcf.rectangulosporium, due to the paucity of the materials obtained. The ophiostomatoid community was dominated by O.ips. This study revealed that a relatively high species diversity of ophiostomatoid fungi are associated with pine infected by B.xylophilus and M.alternatus in China.

Birth-and-Death Evolution and Reticulation of ITS Segments of Metschnikowia andauensis and Metschnikowia fructicola rDNA Repeats

The internal transcribed spacer (ITS) region (ITS1, 5.8S rDNA, and ITS2) separates the genes coding for the SSU 18S and the LSU 26S genes in the rDNA units which are organized into long tandem arrays in the overwhelming majority of fungi. As members of a multigenic family, these units are subject of concerted evolution, which homogenizes their sequences. Exceptions have been observed in certain groups of plants and in a few fungal species. In our previous study we described exceptionally high degree of sequence diversity in the D1/D2 domains of two pulcherrimin-producing Metschnikowia (Saccharomycotina) species which appeared to evolve by reticulation. The major goals of this study were the examination of the diversity of the ITS segments and their evolution. We show that the ITS sequences of these species are not homogenized either, differ from each other by up to 38 substitutions and indels which have dramatic effects on the predicted secondary structures of the transcripts. The high intragenomic diversity makes the D1/D2 domains and the ITS spacers unsuitable for barcoding of these species and therefore the taxonomic position of strains previously assigned to them needs revision. By analyzing the genome sequence of the M. fructicola type strain, we also show that the rDNA of this species is fragmented, contains pseudogenes and thus evolves by the birth-and-death mechanism rather than by homogenisation, which is unusual in yeasts. The results of the network analysis of the sequences further indicate that the ITS regions are also involved in reticulation. M. andauensis and M. fructicola can form interspecies hybrids and their hybrids segregate, providing thus possibilities for reticulation of the rDNA repeats.

CONSTAX: a tool for improved taxonomic resolution of environmental fungal ITS sequences

BACKGROUND

One of the most crucial steps in high-throughput sequence-based microbiome studies is the taxonomic assignment of sequences belonging to operational taxonomic units (OTUs). Without taxonomic classification, functional and biological information of microbial communities cannot be inferred or interpreted. The internal transcribed spacer (ITS) region of the ribosomal DNA is the conventional marker region for fungal community studies. While bioinformatics pipelines that cluster reads into OTUs have received much attention in the literature, less attention has been given to the taxonomic classification of these sequences, upon which biological inference is dependent.

RESULTS

Here we compare how three common fungal OTU taxonomic assignment tools (RDP Classifier, UTAX, and SINTAX) handle ITS fungal sequence data. The classification power, defined as the proportion of assigned OTUs at a given taxonomic rank, varied among the classifiers. Classifiers were generally consistent (assignment of the same taxonomy to a given OTU) across datasets and ranks; a small number of OTUs were assigned unique classifications across programs. We developed CONSTAX (CONSensus TAXonomy), a Python tool that compares taxonomic classifications of the three programs and merges them into an improved consensus taxonomy. This tool also produces summary classification outputs that are useful for downstream analyses.

CONCLUSIONS

Our results demonstrate that independent taxonomy assignment tools classify unique members of the fungal community, and greater classification power is realized by generating consensus taxonomy of available classifiers with CONSTAX.

Effectiveness of ITS and sub-regions as DNA barcode markers for the identification of Basidiomycota (Fungi)

BACKGROUND

Fungi are among the most abundant and diverse organisms on Earth. However, a substantial amount of the species diversity, relationships, habitats, and life strategies of these microorganisms remain to be discovered and characterized. One important factor hindering progress is the difficulty in correctly identifying fungi. Morphological and molecular characteristics have been applied in such tasks. Later, DNA barcoding has emerged as a new method for the rapid and reliable identification of species. The nrITS region is considered the universal barcode of Fungi, and the ITS1 and ITS2 sub-regions have been applied as metabarcoding markers. In this study, we performed a large-scale analysis of all the available Basidiomycota sequences from GenBank. We carried out a rigorous trimming of the initial dataset based in methodological principals of DNA Barcoding. Two different approaches (PCI and barcode gap) were used to determine the performance of the complete ITS region and sub-regions.

RESULTS

For most of the Basidiomycota genera, the three genomic markers performed similarly, i.e., when one was considered a good marker for the identification of a genus, the others were also; the same results were observed when the performance was insufficient. However, based on barcode gap analyses, we identified genomic markers that had a superior identification performance than the others and genomic markers that were not indicated for the identification of some genera. Notably, neither the complete ITS nor the sub-regions were useful in identifying 11 of the 113 Basidiomycota genera. The complex phylogenetic relationships and the presence of cryptic species in some genera are possible explanations of this limitation and are discussed.

CONCLUSIONS

Knowledge regarding the efficiency and limitations of the barcode markers that are currently used for the identification of organisms is crucial because it benefits research in many areas. Our study provides information that may guide researchers in choosing the most suitable genomic markers for identifying Basidiomycota species.

Determining threshold values for barcoding fungi: lessons from Cortinarius (Basidiomycota), a highly diverse and widespread ectomycorrhizal genus

Different distance-based threshold selection approaches were used to assess and compare use of the internal transcribed spacer (ITS) region to distinguish among 901 Cortinarius species represented by >3000 collections. Sources of error associated with genetic markers and selection approaches were explored and evaluated using MOTUs from genus and lineage based-alignments. Our study indicates that 1%-2% more species can be distinguished by using the full-length ITS barcode as compared to either the ITS1 or ITS2 regions alone. Optimal threshold values for different picking approaches and genetic marker lengths inferred from a subset of species containing major lineages ranged from 97.0% to 99.5% sequence similarity using clustering optimization and UNITE SH, and from 1% to 2% sequence dissimilarity with CROP. Errors for the optimal cutoff ranged from 0% to 70%, and these can be reduced to a maximum of 22% when excluding species lacking a barcode gap. A threshold value of 99% is suitable for distinguishing species in the majority of lineages in the genus using the entire ITS region but only 90% of the species could be identified using just the ITS1 or ITS2 region. Prior identification of species, lacking barcode gaps and their subsequent separate analyses, maximized the accuracy of threshold approaches.

ITS1, 5.8S and ITS2 secondary structure modelling for intra-specific differentiation among species of the Colletotrichum gloeosporioides sensu lato species complex

The Colletotrichum gloeosporioides species complex is among the most destructive fungal plant pathogens in the world, however, identification of member species which are of quarantine importance is impacted by a number of factors that negatively affect species identification. Structural information of the rRNA marker may be considered to be a conserved marker which can be used as supplementary information for possible species identification. The difficulty in using ITS rDNA sequences for identification lies in the low level of sequence variation at the intra-specific level and the generation of artificially-induced sequence variation due to errors in polymerization of the ITS array during DNA replication. Type and query ITS sequences were subjected to sequence analyses prior to generation of predicted consensus secondary structures, including the pattern of nucleotide polymorphisms and number of indel haplotypes, GC content, and detection of artificially-induced sequence variation. Data pertaining to structure stability, the presence of conserved motifs in secondary structures and mapping of all sequences onto the consensus C. gloeosporioides sensu stricto secondary structure for ITS1, 5.8S and ITS2 markers was then carried out. Motifs that are evolutionarily conserved among eukaryotes were found for all ITS1, 5.8S and ITS2 sequences. The sequences exhibited conserved features typical of functional rRNAs. Generally, polymorphisms occurred within less conserved regions and were seen as bulges, internal and terminal loops or non-canonical G–U base-pairs within regions of the double stranded helices. Importantly, there were also taxonomic motifs and base changes that were unique to specific taxa and which may be used to support intra-specific identification of members of the C. gloeosporioides sensu lato species complex.

An open source chimera checker for the fungal ITS region

The internal transcribed spacer (ITS) region of the nuclear ribosomal repeat unit holds a central position in the pursuit of the taxonomic affiliation of fungi recovered through environmental sampling. Newly generated fungal ITS sequences are typically compared against the International Nucleotide Sequence Databases for a species or genus name using the sequence similarity software suite blast. Such searches are not without complications however, and one of them is the presence of chimeric entries among the query or reference sequences. Chimeras are artificial sequences, generated unintentionally during the polymerase chain reaction step, that feature sequence data from two (or possibly more) distinct species. Available software solutions for chimera control do not readily target the fungal ITS region, but the present study introduces a blast-based open source software package (available at http://www.emerencia.org/chimerachecker.html) to examine newly generated fungal ITS sequences for the presence of potentially chimeric elements in batch mode. We used the software package on a random set of 12 300 environmental fungal ITS sequences in the public sequence databases and found 1.5% of the entries to be chimeric at the ordinal level after manual verification of the results. The proportion of chimeras in the sequence databases can be hypothesized to increase as emerging sequencing technologies drawing from pooled DNA samples are becoming important tools in molecular ecology research.

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.

CloVR-ITS: Automated internal transcribed spacer amplicon sequence analysis pipeline for the characterization of fungal microbiota

BACKGROUND

Besides the development of comprehensive tools for high-throughput 16S ribosomal RNA amplicon sequence analysis, there exists a growing need for protocols emphasizing alternative phylogenetic markers such as those representing eukaryotic organisms.

RESULTS

Here we introduce CloVR-ITS, an automated pipeline for comparative analysis of internal transcribed spacer (ITS) pyrosequences amplified from metagenomic DNA isolates and representing fungal species. This pipeline performs a variety of steps similar to those commonly used for 16S rRNA amplicon sequence analysis, including preprocessing for quality, chimera detection, clustering of sequences into operational taxonomic units (OTUs), taxonomic assignment (at class, order, family, genus, and species levels) and statistical analysis of sample groups of interest based on user-provided information. Using ITS amplicon pyrosequencing data from a previous human gastric fluid study, we demonstrate the utility of CloVR-ITS for fungal microbiota analysis and provide runtime and cost examples, including analysis of extremely large datasets on the cloud. We show that the largest fractions of reads from the stomach fluid samples were assigned to Dothideomycetes, Saccharomycetes, Agaricomycetes and Sordariomycetes but that all samples were dominated by sequences that could not be taxonomically classified. Representatives of the Candida genus were identified in all samples, most notably C. quercitrusa, while sequence reads assigned to the Aspergillus genus were only identified in a subset of samples. CloVR-ITS is made available as a pre-installed, automated, and portable software pipeline for cloud-friendly execution as part of the CloVR virtual machine package (http://clovr.org).

CONCLUSION

The CloVR-ITS pipeline provides fungal microbiota analysis that can be complementary to bacterial 16S rRNA and total metagenome sequence analysis allowing for more comprehensive studies of environmental and host-associated microbial communities.

ITS1 versus ITS2 as DNA metabarcodes for fungi

The nuclear ribosomal Internal Transcribed Spacer ITS region is widely used as a DNA metabarcoding marker to characterize the diversity and composition of fungal communities. In amplicon pyrosequencing studies of fungal diversity, one of the spacers ITS1 or ITS2 of the ITS region is normally used. In this methodological study we evaluate the usability of ITS1 vs. ITS2 as a DNA metabarcoding marker for fungi. We analyse three data sets: two comprising ITS1 and ITS2 sequences of known taxonomic affiliations and a third comprising ITS1 and ITS2 environmental amplicon pyrosequencing data. Clustering analyses of sequences with known taxonomy using the bioinformatics pipeline ClustEx revealed that a 97% similarity cut-off represent a reasonable threshold for estimating the number of known species in the data sets for both ITS1 and ITS2. However, no single threshold value worked well for all fungi at the same time within the curated UNITE database, and we found that the Operational Taxonomic Unit (OTU) concept is not easily translated into the level of species because many species are distributed over several clusters. Clustering analyses of the 134 692 ITS1 and ITS2 pyrosequences using a 97% similarity cut-off revealed a high similarity between the two data sets when it comes to taxonomic coverage. Although some groups are under- or unrepresented in the two data sets due to, e.g. primer mismatches, our results indicate that ITS1 and ITS2 to a large extent yield similar results when used as DNA metabarcodes for fungi.

Nucleotide polymorphism in the 5.8S nrDNA gene and internal transcribed spacers in Phakopsora pachyrhizi viewed from structural models

The assessment of nucleotide polymorphisms in environmental samples of obligate pathogens requires DNA amplification through the polymerase chain reaction (PCR) and bacterial cloning of PCR products prior to sequencing. The drawback of this strategy is that it can give rise to false polymorphisms owing to DNA polymerase misincorporation during PCR or bacterial cloning. We investigated patterns of nucleotide polymorphism in the internal transcribed spacer (ITS) region for Phakopsora pachyrhizi, an obligate biotrophic fungus that causes the Asian soybean rust. Field-collected samples of P. pachyrhizi were obtained from all major soybean production areas worldwide, including Brazil and the United States. Bacterially-cloned, PCR products were obtained using a high fidelity DNA polymerase. A total of 370 ITS sequences that were subjected to an array of complementary sequence analyses, which included analyses of secondary structure stability, the pattern of nucleotide polymorphisms, GC content, and the presence of conserved motifs. The sequences exhibited features of functional rRNAs. Overall, polymorphisms took place within less conserved motives, such as loops and bulges; alternatively, they gave rise to non-canonical G-U pairs within conserved regions of double stranded helices. We discuss the usefulness of structural analyses to filter out putative 'suspicious' bacterially cloned ITS sequences, thus keeping artificially-induced sequence variation to a minimum.

Low evolutionary diversification in a widespread and abundant uncultured protist (MAST-4)

Recent culture-independent studies of marine planktonic protists have unveiled a large diversity at all phylogenetic scales and the existence of novel groups. MAST-4 represents one of these novel uncultured lineages, and it is composed of small (~2 μm) bacterivorous eukaryotes that are widely distributed in marine systems. MAST-4 accounts for a significant fraction of the marine heterotrophic flagellates at the global level, playing key roles in the marine ecological network. In this study, we investigated the diversity of MAST-4, aiming to assess its limits and structure. Using ribosomal DNA (rDNA) sequences obtained in this study (both pyrosequencing reads and clones with large rDNA operon coverage), complemented with GenBank sequences, we show that MAST-4 is composed of only five main clades, which are well supported by small subunit and large subunit phylogenies. The differences in the conserved regions of the internal transcribed spacers 1 and 2 (ITS1 and ITS2) secondary structures strongly suggest that these five clades are different biological species. Based on intraclade divergence, ITS secondary structures and comparisons of ITS1 and ITS2 trees, we did not find evidence of more than one species within clade A, whereas as many as three species might be present within other clades. Overall, the genetic divergence of MAST-4 was surprisingly low for an organism with a global population size estimated to be around 10(24), indicating a very low evolutionary diversification within the group.

Ribosomal DNA polymorphisms in the yeast Geotrichum candidum

The dimorphic yeast Geotrichum candidum (teleomorph: Galactomyces candidus) is commonly used to inoculate washed-rind and bloomy-rind cheeses. However, little is known about the phylogenetic lineage of this microorganism. We have sequenced the complete 18S, 5.8S, 26S ribosomal RNA genes and their internal transcribed spacers (ITS1) and ITS2 regions (5126 nucleotides) from 18 G. candidum strains from various environmental niches, with a focus on dairy strains. Multiple sequence alignments revealed the presence of 60 polymorphic sites, which is generally unusual for ribosomal DNA (rDNA) within a given species because of the concerted evolution mechanism. This mechanism drives genetic homogenization to prevent the divergent evolution of rDNA copies within individuals. While the polymorphisms observed were mainly substitutions, one insertion/deletion (indel) polymorphism was detected in ITS1. No polymorphic sites were detected downstream from this indel site, that is, in 5.8S and ITS2. More surprisingly, many sequence electrophoregrams generated during the sequencing of the rDNA had dual peaks, suggesting that many individuals exhibited intragenomic rDNA variability. The ITS1-5.8S-ITS2 regions of four strains were cloned. The sequence analysis of 68 clones revealed 32 different ITS1-5.8S-ITS2 variants within these four strains. Depending on the strain, from four to twelve variants were detected, indicating that multiple rDNA copies were present in the genomes of these G. candidum strains. These results contribute to the debate concerning the use of the ITS region for barcoding fungi and suggest that community profiling techniques based on rDNA should be used with caution.

A close-up view on ITS2 evolution and speciation - a case study in the Ulvophyceae (Chlorophyta, Viridiplantae)

BACKGROUND

The second Internal Transcriber Spacer (ITS2) is a fast evolving part of the nuclear-encoded rRNA operon located between the 5.8S and 28S rRNA genes. Based on crossing experiments it has been proposed that even a single Compensatory Base Change (CBC) in helices 2 and 3 of the ITS2 indicates sexual incompatibility and thus separates biological species. Taxa without any CBC in these ITS2 regions were designated as a 'CBC clade'. However, in depth comparative analyses of ITS2 secondary structures, ITS2 phylogeny, the origin of CBCs, and their relationship to biological species have rarely been performed. To gain 'close-up' insights into ITS2 evolution, (1) 86 sequences of ITS2 including secondary structures have been investigated in the green algal order Ulvales (Chlorophyta, Viridiplantae), (2) after recording all existing substitutions, CBCs and hemi-CBCs (hCBCs) were mapped upon the ITS2 phylogeny, rather than merely comparing ITS2 characters among pairs of taxa, and (3) the relation between CBCs, hCBCs, CBC clades, and the taxonomic level of organisms was investigated in detail.

RESULTS

High sequence and length conservation allowed the generation of an ITS2 consensus secondary structure, and introduction of a novel numbering system of ITS2 nucleotides and base pairs. Alignments and analyses were based on this structural information, leading to the following results: (1) in the Ulvales, the presence of a CBC is not linked to any particular taxonomic level, (2) most CBC 'clades' sensu Coleman are paraphyletic, and should rather be termed CBC grades. (3) the phenetic approach of pairwise comparison of sequences can be misleading, and thus, CBCs/hCBCs must be investigated in their evolutionary context, including homoplasy events (4) CBCs and hCBCs in ITS2 helices evolved independently, and we found no evidence for a CBC that originated via a two-fold hCBC substitution.

CONCLUSIONS

Our case study revealed several discrepancies between ITS2 evolution in the Ulvales and generally accepted assumptions underlying ITS2 evolution as e.g. the CBC clade concept. Therefore, we developed a suite of methods providing a critical 'close-up' view into ITS2 evolution by directly tracing the evolutionary history of individual positions, and we caution against a non-critical use of the ITS2 CBC clade concept for species delimitation.

Evolutionary dynamics of the mS952 intron: a novel mitochondrial group II intron encoding a LAGLIDADG homing endonuclease gene

Examination of the mitochondrial small subunit ribosomal RNA (rns) gene of five species of the fungal genus Leptographium revealed that the gene has been invaded at least once at position 952 by a group II intron encoding a LAGLIDADG homing endonuclease gene. Phylogenetic analyses of the intron and homing endonuclease sequences indicated that each element in Leptographium species forms a single clade and is closely related to the group II intron/homing endonuclease gene composite element previously reported at position 952 of the mitochondrial rns gene of Cordyceps species and of Cryphonectria parasitica. The results of an intron survey of the mt rns gene of Leptographium species superimposed onto the phylogenetic analysis of the host organisms suggest that the composite element was transmitted vertically in Leptographium lundbergii. However, its stochastic distribution among strains of L. wingfieldii, L. terebrantis, and L. truncatum suggests that it has been horizontally transmitted by lateral gene transfer among these species, although the random presence of the intron may reflect multiple random loss events. A model is proposed describing the initial invasion of the group II intron in the rns gene of L. lundbergii by a LAGLIDADG homing endonuclease gene and subsequent evolution of this gene to recognize a novel DNA target site, which may now promote the mobility of the intron and homing endonuclease gene as a composite element.

Diversity and decay ability of basidiomycetes isolated from lodgepole pines killed by the mountain pine beetle

When lodgepole pines ( Pinus contorta Douglas ex Louden var. latifolia Engelm. ex S. Watson) that are killed by the mountain pine beetle ( Dendroctonus ponderosae ) and its fungal associates are not harvested, fungal decay can affect wood and fibre properties. Ophiostomatoids stain sapwood but do not affect the structural properties of wood. In contrast, white or brown decay basidiomycetes degrade wood. We isolated both staining and decay fungi from 300 lodgepole pine trees killed by mountain pine beetle at green, red, and grey stages at 10 sites across British Columbia. We retained 224 basidiomycete isolates that we classified into 34 species using morphological and physiological characteristics and rDNA large subunit sequences. The number of basidiomycete species varied from 4 to 14 species per site. We assessed the ability of these fungi to degrade both pine sapwood and heartwood using the soil jar decay test. The highest wood mass losses for both sapwood and heartwood were measured for the brown rot species Fomitopsis pinicola and the white rot Metulodontia and Ganoderma species. The sap rot species Trichaptum abietinum was more damaging for sapwood than for heartwood. A number of species caused more than 50% wood mass losses after 12 weeks at room temperature, suggesting that beetle-killed trees can rapidly lose market value due to degradation of wood structural components.

Rapid identification and detection of pine pathogenic fungi associated with mountain pine beetles by padlock probes

Fifteen million hectares of pine forests in western Canada have been attacked by the mountain pine beetle (Dendroctonus ponderosae; MPB), leading to devastating economic losses. Grosmannia clavigera and Leptographium longiclavatum, are two fungi intimately associated with the beetles, and are crucial components of the epidemic. To detect and discriminate these two closely related pathogens, we utilized a method based on ligase-mediated nucleotide discrimination with padlock probe technology, and signal amplification by hyperbranched rolling circle amplification (HRCA). Two padlock probes were designed to target species-specific single nucleotide polymorphisms (SNPs) located at the inter-generic spacer 2 region and large subunit of the rRNA respectively, which allows discrimination between the two species. Thirty-four strains of G. clavigera and twenty-five strains of L. longiclavatum representing a broad geographic origin were tested with this assay. The HRCA results were largely in agreement with the conventional identification based on morphology or DNA-based methods. Both probes can also efficiently distinguish the two MPB-associated fungi from other fungi in the MPB, as well as other related fungi in the order Ophiostomatales. We also tested this diagnostic method for the direct detection of these fungi from the DNA of MPB. A nested PCR approach was used to enrich amplicons for signal detection. The results confirmed the presence of these two fungi in MPB. Thus, the padlock probe assay coupled with HRCA is a rapid, sensitive and reproducible method for the identification and detection of these ophiostomatoid fungi.