Early Ongoing Speciation of Ogataea uvarum Sp. Nov. Within the Grape Ecosystem Revealed by the Internal Variability Among the rDNA Operon Repeats

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.

Yeast Starter Culture Identification to Produce of Red Wines with Enhanced Antioxidant Content

Grape variety, quality, geographic origins and phytopathology can influence the amount of polyphenols that accumulate in grape tissues. Polyphenols in wine not only shape their organoleptic characteristics but also significantly contribute to the positive impact that this beverage has on human health. However, during the winemaking process, the total polyphenol content is substantially reduced due to the adsorption onto yeast wall polymers and subsequent lees separation. Despite this, limited information is available regarding the influence of the yeast starter strain on the polyphenolic profile of wine. To address this issue, a population consisting of 136 Saccharomyces cerevisiae strains was analyzed to identify those with a diminished ability to adsorb polyphenols. Firstly, the reduction in concentration of polyphenolic compounds associated to each strain was studied by assaying Total Phenolic Content (TPC) and Trolox Equivalent Antioxidant Capacity (TEAC) in the wines produced by micro-scale must fermentation. A total of 29 strains exhibiting a TPC and TEAC reduction ≤ 50%, when compared to that detected in the utilized grape must were identified and the nine most-promising strains were further validated by larger-scale vinification. Physico-chemical analyses of the resulting wines led to the identification of four strains, namely ITEM6920, ITEM9500, ITEM9507 and ITEM9508 which showed, compared to the control wine, a TPC and TEAC reduction ≤ 20 in the produced wines. They were denoted by a significant (p < 0.05) increased amount of anthocyanin, quercetin and trans-coutaric acid, minimal volatile acidity (<0.2 g/L), absence of undesirable metabolites and a well-balanced volatile profile. As far as we know, this investigation represents the first clonal selection of yeast strains aimed at the identifying "functional" fermentation starters, thereby enabling the production of regional wines with enriched polyphenolic content.

MinION Sequencing of Yeast Mock Communities To Assess the Effect of Databases and ITS-LSU Markers on the Reliability of Metabarcoding Analysis

Microbial communities play key roles both for humans and the environment. They are involved in ecosystem functions, maintaining their stability, and provide important services, such as carbon cycle and nitrogen cycle. Acting both as symbionts and as pathogens, description of the structure and composition of these communities is important. Metabarcoding uses ribosomal DNA (rDNA) (eukaryotic) or rRNA gene (prokaryotic) sequences for identification of species present in a site and measuring their abundance. This procedure requires several technical steps that could be source of bias producing a distorted view of the real community composition. In this work, we took advantage of an innovative "long-read" next-generation sequencing (NGS) technology (MinION) amplifying the DNA spanning from the internal transcribed spacer (ITS) to large subunit (LSU) that can be read simultaneously in this platform, providing more information than "short-read" systems. The experimental system consisted of six fungal mock communities composed of species present at various relative amounts to mimic natural situations characterized by predominant and low-frequency species. The influence of the sequencing platform (MinION and Illumina MiSeq) and the effect of different reference databases and marker sequences on metagenomic identification of species were evaluated. The results showed that the ITS-based database provided more accurate species identification than LSU. Furthermore, a procedure based on a preliminary identification with standard reference databases followed by the production of custom databases, including only the best outputs of the first step, is proposed. This additional step improved the estimate of species proportion of the mock communities and reduced the number of ghost species not really present in the simulated communities. IMPORTANCE Metagenomic analyses are fundamental in many research areas; therefore, improvement of methods and protocols for the description of microbial communities becomes more and more necessary. Long-read sequencing could be used for reducing biases due to the multicopy nature of rDNA sequences and short-read limitations. However, these novel technologies need to be assessed and standardized with controlled experiments, such as mock communities. The interest behind this work was to evaluate how long reads performed identification and quantification of species mixed in precise proportions and how the choice of database affects such analyses. Development of a pipeline that mitigates the effect of the barcoding sequences and the impact of the reference database on metagenomic analyses can help microbiome studies go one step further.

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.

When barcoding fails: genome chimerisation (admixing) and reticulation obscure phylogenetic and taxonomic relationships

DNA barcoding is based on the premise that the barcode sequences can distinguish individuals (strains) of different species because their sequence variation between species exceeds that within species. The primary barcodes used in fungal and yeast taxonomy are the ITS segments and the LSU (large subunit) D1/D2 domain of the homogenized multicopy rDNA repeats. The secondary barcodes are conserved segments of protein‐encoding genes, which usually have single copies in haploid genomes. This study shows that the analysis of barcode sequences fails to reconstruct accurate species trees and differentiate species when the organisms have chimeric genomes composed of admixed mosaics of different origins. It is shown that the type strains of 10 species of the pulcherrima clade of the ascomycetous yeast genus Metschnikowia cannot be differentiated with standard barcodes because their intragenomic diversity is comparable to or even higher than the interstrain diversity. The analysis of a large group of genes of the sequenced genomes of the clade and the viability and segregation of the hybrids of ex‐type strains indicate that the high intragenomic barcode differences can be attributed to admixed genome structures. Because of the mosaic structures of the genomes, the rDNA repeats do not form continuous arrays and thus cannot be homogenized. Since the highly diverse ITS and D1/D2 sequences of the type strains form a continuous pool including pseudogenes, the evolution of their rDNA appears to involve reticulation. The secondary barcode sequences and the nonbarcode genes included in the analysis show incongruent phylogenetic relationships among the type strains, which can also be attributed to differences in the phylogenetic histories of the genes.

Using High-Throughput Amplicon Sequencing to Evaluate Intragenomic Variation and Accuracy in Species Identification of Cordyceps Species

While recent sequencing technologies (third generation sequencing) can successfully sequence all copies of nuclear ribosomal DNA (rDNA) markers present within a genome and offer insights into the intragenomic variation of these markers, high intragenomic variation can be a source of confusion for high-throughput species identification using such technologies. High-throughput (HT) amplicon sequencing via PacBio SEQUEL I was used to evaluate the intragenomic variation of the ITS region and D1-D2 LSU domains in nine Cordyceps species, and the accuracy of such technology to identify these species based on molecular phylogenies was also assessed. PacBio sequences within strains showed variable level of intragenomic variation among the studied Cordyceps species with C. blackwelliae showing greater variation than the others. Some variants from a mix of species clustered together outside their respective species of origin, indicative of intragenomic variation that escaped concerted evolution shared between species. Proper selection of consensus sequences from HT amplicon sequencing is a challenge for interpretation of correct species identification. PacBio consensus sequences with the highest number of reads represent the major variants within a genome and gave the best results in terms of species identification.

What Is the Best Lens? Comparing the Resolution Power of Genome-Derived Markers and Standard Barcodes

Fungal species delimitation was traditionally carried out with multicopy ribosomal RNA (rRNA) genes, principally for their ease of amplification. Since the efficacy of these markers has been questioned, single-copy protein-encoding genes have been proposed alone or in combination for Multi-Locus Sequence Typing (MLST). In this context, the role of the many sequences obtained with Next-Generation Sequencing (NGS) techniques, in both genomics and metagenomics, further pushes toward an analysis of the efficacy of NGS-derived markers and of the metrics to evaluate the marker efficacy in discriminating fungal species. This paper aims at proposing MeTRe (Mean Taxonomic Resolution), a novel index that could be used both for measuring marker efficacy and for assessing the actual resolution (i.e., the level of separation) between species obtained with different markers or their combinations. In this paper, we described and then employed this index to compare the efficacy of two rRNAs and four single-copy markers obtained from public databases as both an amplicon-based approach and genome-derived sequences. Two different groups of species were used, one with a pathogenic species of Candida that was characterized by relatively well-separated taxa, whereas the other, comprising some relevant species of the sensu stricto group of the genus Saccharomyces, included close species and interspecific hybrids. The results showed the ability of MeTRe to evaluate marker efficacy in general and genome-derived markers specifically.

Single Strain High-Depth NGS Reveals High rDNA (ITS-LSU) Variability in the Four Prevalent Pathogenic Species of the Genus Candida

Ribosomal RNA in fungi is encoded by a series of genes and spacers included in a large operon present in 100 tandem repeats, normally in a single locus. The multigene nature of this locus was somehow masked by Sanger sequencing, which produces a single sequence reporting the prevalent nucleotide of each site. The introduction of next generation sequencing led to deeper knowledge of the individual sequences (reads) and therefore of the variants between the same DNA sequences located in different tandem repeats. In this framework, NGS sequencing of the rDNA region was used to elucidate the extent of intra- and inter-genomic variation at both the strain and species level. Specifically, the use of an innovative NGS technique allowed the high-throughput high-depth sequencing of the ITS1-LSU D1/D2 amplicons of 252 strains belonging to four opportunistic yeast species of the genus Candida. Results showed the presence of a large extent of variability among strains and species. These variants were differently distributed throughout the analyzed regions with a higher concentration within the Internally Transcribed Spacer (ITS) region, suggesting that concerted evolution was not able to totally homogenize these sequences. Both the internal variability and the SNPs between strain can be used for a deep typing of the strains and to study their ecology.

Homoplasy as an Auxiliary Criterion for Species Delimitation

Homoplasy is a sort of noise in phylogenetic reconstructions, due to the accumulation of backmutations, convergent evolution and horizontal gene transfer (HGT), which is considered the major trigger of homoplasy in microorganism for its massive presence. It is also known that homoplasy increases with the complexity of the tree with both real and simulated data. In this paper, we analyzed the variation of homoplasy with the two widely used taxonomic markers ITS and LSU in four taxonomic models characterized by differences in the intra-specific distances. An algorithm (HomoDist) was developed to analyze the homoplasy index (HI) variation upon addition of a single element (strain or species) in increasing distance from a starting element. This algorithm allows to follow changes of the consistency index (CI), complementary to the HI, with the increase of the number of taxa and with the increase of the distance among elements. Results show that homoplasy increases-as expected-with the number of taxa, but also as a function of the overall distance among species, often with an almost linear relationship between distance and HI. No HI change was observed in trees with few taxa spanning through short distances, indicating that this noise is not prohibitive in this context, although the analysis of the ratio between HI and distance can be recommended as a criterion for tree acceptance. The absence of large changes of the HI within the species, and its increase when new species are added by HomoDist, suggest that homoplasy variation can be used as an auxiliary test in distance-based species delimitation with any type of marker.