rRNA Pseudogenes in Filamentous Ascomycetes as Revealed by Genome Data

Elucidation of intragenomic variation of ribosomal DNA sequences in the enigmatic fungal genus Ceraceosorus, including a newly described species Ceraceosorus americanus

Multicopy nuclear ribosomal DNA (rDNA) genes have been used as markers for fungal identification for three decades. The rDNA sequences in a genome are thought to be homogeneous due to concerted evolution. However, intragenomic variation of rDNA sequences has recently been observed in many fungi, which may make fungal identification and species abundance estimation based on these loci problematic. Ceraceosorus is an enigmatic genus in the smut lineage Ustilaginomycotina for which very limited distribution data exist. Our previous research demonstrated intragenomic variation in the internal transcribed spacer (ITS1-5.8S-ITS2) region of two Ceraceosorus species. In this study, we described the fourth known species of Ceraceosorus, C. americanus, isolated from an asymptomatic rosemary leaf collected in Louisiana, USA. This is the first report of this genus in the Americas. We then selected all four known Ceraceosorus species, plus exemplar smut fungi representing all major lineages of subphylum Ustilaginomycotina, to examine sequence heterogeneity in three regions of the rDNA repeat (partial 18S, ITS, and partial 28S regions). Three methods were used: PCR-cloning-Sanger sequencing, targeted amplicon high-throughput sequencing, and whole-genome shotgun high-throughput sequencing. Our results show that Ceraceosorus is the only sampled fungal genus in Ustilaginomycotina with significant intragenomic variation at the ITS, with up to 25 nucleotide variant sites in the ITS1-5.8S-ITS2 region and 2.6% divergence among analyzed ITS haplotypes. We found many conflicting patterns across the three detection methods, with up to 27 conflicting variant sites recorded from a single individual. At least 40% of the conflicting patterns are possibly due to PCR-cloning-sequencing errors, as the corresponding variant sites were not observed in the other detection methods. Based on our data and the literature, we evaluated the characteristics and advantages/disadvantages of each detection method. Finally, a model for how intragenomic variation in the rDNA copies within a genome may arise is presented.

Whole genome sequencing of Wilsonomyces carpophilus, an incitant of shot hole disease in stone fruits: insights into secreted proteins of a necrotrophic fungal repository

BACKGROUND

Shot hole is one of the important fungal diseases in stone fruits viz., peach, plum, apricot and cherry caused by Wilsonomyces carpophilus and almond among nut crops. Fungicides significantly decrease the disease. Pathogenicity studies proved a wide host range of the pathogen infecting all stone fruits and almond among the nut crops, however, the mechanism underlying host-pathogen interaction is still unknown. Molecular detection of the pathogen using polymerase chain reaction (PCR) based simple sequence repeat (SSR) markers is also unknown due to the unavailability of the pathogen genome.

METHODS AND RESULTS

We examined the morphology, pathology and genomics of the Wilsonomyces carpophilus. Whole genome sequencing of the W. carpophilus was carried out by Illumina HiSeq and PacBio high throughput sequencing plate-forms through hybrid assembly. Constant selection pressure alters the molecular mechanism of the pathogen causing disease. The studies revealed that the necrotrophs are more lethal with a complex pathogenicity mechanism and little-understood effector repositories. The different isolates of necrotrophic fungus W. carpophilus causing shot hole in stone fruits namely peach, plum, apricot and cherry, and almonds among the nut crops showed a significant variation in their morphology, however, the probability value (p = 0.29) suggests in-significant difference in the pathogenicity. Here, we reported draft genome of W. carpophilus of size 29.9 Mb (Accession number: PRJNA791904). A total of 10,901 protein-coding genes were predicted, including heterokaryon incompatibility genes, cytochrome-p450 genes, kinases, sugar transporters among others. We found 2851 simple sequence repeats (SSRs), tRNAs, rRNAs and pseudogenes in the genome. The most prominent proteins showing necrotrophic lifestyle of the pathogen were hydrolases, polysaccharide-degrading enzymes, esterolytic, lipolytic, and proteolytic enzymes accounted for 225 released proteins. Among the 223 fungal species, top-hit species distribution revealed the majority of hits against the Pyrenochaeta species followed by Ascochyta rabiei and Alternaria alternata.

CONCLUSION

Draft genome of W. carpophilus is 29.9 Mb based on Illumina HiSeq and PacBio hybrid assembly. The necrotrophs are more lethal with a complex pathogenicity mechanism. A significant variation in morphology was observed in different pathogen isolates. A total of 10,901 protein-coding genes were predicted in the pathogen genome including heterokaryon incompatibility, cytochrome-p450 genes, kinases and sugar transporters. We found 2851 SSRs, tRNAs, rRNAs and pseudogenes, and prominent proteins showing necrotrophic lifestyle such as hydrolases, polysaccharide-degrading enzymes, esterolytic, lipolytic and proteolytic enzymes. The top-hit species distribution were against the Pyrenochaeta spp. followed by Ascochyta rabiei.

Improving Environmental DNA Sensitivity for Dreissenid Mussels by Targeting Tandem Repeat Regions of the Mitochondrial Genome

The recent genetic revolution through the analysis of aquatic environmental DNA (eDNA) has become a powerful tool for improving the detection of rare and/or invasive species. For the majority of eDNA studies, genetic assays are designed to target mitochondrial genes commonly referred to as “barcode” regions. However, unlike the typical structure of an animal mitochondrial genome, those for the invasive zebra and quagga mussels are greatly expanded with large extended tandem repeat regions. These sections of repeated DNA can appear hundreds of times within the genome compared to a single copy for the mitochondrial barcode genes. This higher number of target copies per mitochondrial genome presents an opportunity to increase eDNA assay sensitivity for these species. Therefore, we designed and evaluated new eDNA assays to target the extended repeat sections for both zebra and quagga mussels. These assays lower the limit of detection of genomic DNA by 100-fold for zebra mussels and 10-fold for quagga mussels. Additionally, these newly developed assays provided longer durations of detection during degradation mesocosm experiments and greater sensitivity for eDNA detection from water samples collected across western Lake Erie compared to standard assays targeting mitochondrial genes. This work illustrates how understanding the complete genomic structure of an organism can improve eDNA analysis.

Caveats of the internal transcribed spacer region as a barcode to resolve species boundaries in Diaporthe

Species in Diaporthe are largely reported as important plant pathogens. Identification of species in this genus has been complemented by morphological and molecular features. However, one important factor delaying this process is the struggle to formulate robust species concepts to create adequate international phytosanitary measures. Regardless of the wide use of the internal transcribed spacer (ITS) rDNA region, established as the primary DNA barcode for fungi, the tendency for intraspecific variation has been reported, misleading interpretation of phylogenetic analyses. Therefore, the present study aimed to illustrate, using specific examples, how the ITS region may be problematic for species delimitation. We showed that the ITS region is highly variable, with strains of Diaporthe malorum and Diaporthe novem falling into more than one clade, which if analyzed on their own, would be likely recognized as distinct taxa. Divergent ITS paralogs were also proven to coexist within the genome of D. novem. We also suggest that ITS may have escaped from concerted evolution or has undergone a duplication event. Furthermore, this study reports for the first time the existence of a putative hybrid in the genus Diaporthe. Our findings offer new clues towards the intraspecific and intragenomic variation in the ITS region, raising questions about its value for barcoding, i.e., identifying species in the genus Diaporthe. Therefore, we recommend that the ITS region be analyzed cautiously and always compared for congruence prior to description of novel taxa.

RIP mutated ITS genes in populations of Ophiocordyceps sinensis and their implications for molecular systematics

Different hypotheses have been proposed to interpret the observed unusual ITS (internal transcribed spacer) sequences in Ophiocordyceps sinensis. The coexistence of diverged ITS paralogs in a single genome was previously shown by amplifying the ITS region from mono-ascospore isolates using specific primers designed for different ITS paralog groups. Among those paralogs, are AT-biased ITS sequences which were hypothesized to result from repeat-induced point mutation (RIP). This is a process that detects and mutates repetitive DNA and frequently leads to epigenetic silencing, and these mutations have been interpreted as pseudogenes. Here we investigate the occurrence and frequency of ITS pseudogenes in populations of O. sinensis using large-scale sampling, and discusses the implications of ITS pseudogenes for fungal phylogenetic and evolutionary studies. Our results demonstrate a wide distribution of ITS pseudogenes amongst different geographic populations, and indicate how ITS pseudogenes can contribute to the reconstruction of the evolutionary history of the species.

Unambiguous identification of fungi: where do we stand and how accurate and precise is fungal DNA barcoding?

True fungi (Fungi) and fungus-like organisms (e.g. Mycetozoa, Oomycota) constitute the second largest group of organisms based on global richness estimates, with around 3 million predicted species. Compared to plants and animals, fungi have simple body plans with often morphologically and ecologically obscure structures. This poses challenges for accurate and precise identifications. Here we provide a conceptual framework for the identification of fungi, encouraging the approach of integrative (polyphasic) taxonomy for species delimitation, i.e. the combination of genealogy (phylogeny), phenotype (including autecology), and reproductive biology (when feasible). This allows objective evaluation of diagnostic characters, either phenotypic or molecular or both. Verification of identifications is crucial but often neglected. Because of clade-specific evolutionary histories, there is currently no single tool for the identification of fungi, although DNA barcoding using the internal transcribed spacer (ITS) remains a first diagnosis, particularly in metabarcoding studies. Secondary DNA barcodes are increasingly implemented for groups where ITS does not provide sufficient precision. Issues of pairwise sequence similarity-based identifications and OTU clustering are discussed, and multiple sequence alignment-based phylogenetic approaches with subsequent verification are recommended as more accurate alternatives. In metabarcoding approaches, the trade-off between speed and accuracy and precision of molecular identifications must be carefully considered. Intragenomic variation of the ITS and other barcoding markers should be properly documented, as phylotype diversity is not necessarily a proxy of species richness. Important strategies to improve molecular identification of fungi are: (1) broadly document intraspecific and intragenomic variation of barcoding markers; (2) substantially expand sequence repositories, focusing on undersampled clades and missing taxa; (3) improve curation of sequence labels in primary repositories and substantially increase the number of sequences based on verified material; (4) link sequence data to digital information of voucher specimens including imagery. In parallel, technological improvements to genome sequencing offer promising alternatives to DNA barcoding in the future. Despite the prevalence of DNA-based fungal taxonomy, phenotype-based approaches remain an important strategy to catalog the global diversity of fungi and establish initial species hypotheses.

Intragenomic polymorphisms in the ITS region of high-quality genomes of the Hypoxylaceae (Xylariales, Ascomycota)

The internal transcribed spacer (ITS) region of the ribosomal DNA (rDNA) has been established (and is generally accepted) as a primary “universal” genetic barcode for fungi for many years, but the actual value for taxonomy has been heavily disputed among mycologists. Recently, twelve draft genome sequences, mainly derived from type species of the family Hypoxylaceae (Xylariales, Ascomycota) and the ex-epitype strain of Xylaria hypoxylon have become available during the course of a large phylogenomic study that was primarily aimed at establishing a correlation between the existing multi-gene-based genealogy with a genome-based phylogeny and the discovery of novel biosynthetic gene clusters encoding for secondary metabolites. The genome sequences were obtained using combinations of Illumina and Oxford nanopore technologies or PacBio sequencing, respectively, and resulted in high-quality sequences with an average N50 of 3.2 Mbp. While the main results will be published concurrently in a separate paper, the current case study was dedicated to the detection of ITS nrDNA copies in the genomes, in an attempt to explain certain incongruities and apparent mismatches between phenotypes and genotypes that had been observed during previous polyphasic studies. The results revealed that all of the studied strains had at least three copies of rDNA in their genomes, with Hypoxylon fragiforme having at least 19 copies of the ITS region, followed by Xylaria hypoxylon with at least 13 copies. Several of the genomes contained 2–3 copies that were nearly identical, but in some cases drastic differences, below 97% identity were observed. In one case, ascribable to the presence of a pseudogene, the deviations of the ITS sequences from the same genome resulted in only ca. 90% of overall homology. These results are discussed in the scope of the current trends to use ITS data for species recognition and segregation of fungi. We propose that additional genomes should be checked for such ITS polymorphisms to reassess the validity of this non-coding part of the fungal DNA for molecular identification.

Intragenomic internal transcribed spacer 2 variation in a genus of parasitoid wasps (Hymenoptera: Braconidae): implications for accurate species delimitation and phylogenetic analysis

A recent DNA barcoding study of Australian microgastrines (Hymenoptera: Braconidae) sought to use next-generation sequencing of the cytochrome c oxidase subunit 1 (COI) barcoding gene region, the wingless (WG) gene and the internal transcribed spacer 2 (ITS2) to delimit molecular species in a highly diverse group of parasitic wasps. Large intragenomic distances between ITS2 variants, often larger than the average interspecific variation, caused difficulties in using ITS2 for species delimitation in both threshold and tree-based approaches, and the gene was not included in the reported results of the previous DNA barcoding study. We here report on the intragenomic, and the intra- and interspecies, variation in ITS2in the microgastrine genus Diolcogasterto further investigate the value of ITS2as a marker for species delimitation and phylogenetics of the Microgastrinae. Distinctive intragenomic variant patterns were found in different species of Diolcogaster, with some species possessing a single major variant, and others possessing many divergent variants. Characterizing intragenomic variation of ITS2is critical as it is a widely used marker in hymenopteran phylogenetics and species delimitation, and large intragenomic distances such as those found in this study may obscure phylogenetic signal.

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.

Considerations and consequences of allowing DNA sequence data as types of fungal taxa

Nomenclatural type definitions are one of the most important concepts in biological nomenclature. Being physical objects that can be re-studied by other researchers, types permanently link taxonomy (an artificial agreement to classify biological diversity) with nomenclature (an artificial agreement to name biological diversity). Two proposals to amend the International Code of Nomenclature for algae, fungi, and plants (ICN), allowing DNA sequences alone (of any region and extent) to serve as types of taxon names for voucherless fungi (mainly putative taxa from environmental DNA sequences), have been submitted to be voted on at the 11th International Mycological Congress (Puerto Rico, July 2018). We consider various genetic processes affecting the distribution of alleles among taxa and find that alleles may not consistently and uniquely represent the species within which they are contained. Should the proposals be accepted, the meaning of nomenclatural types would change in a fundamental way from physical objects as sources of data to the data themselves. Such changes are conducive to irreproducible science, the potential typification on artefactual data, and massive creation of names with low information content, ultimately causing nomenclatural instability and unnecessary work for future researchers that would stall future explorations of fungal diversity. We conclude that the acceptance of DNA sequences alone as types of names of taxa, under the terms used in the current proposals, is unnecessary and would not solve the problem of naming putative taxa known only from DNA sequences in a scientifically defensible way. As an alternative, we highlight the use of formulas for naming putative taxa (candidate taxa) that do not require any modification of the ICN.