NGS barcode sequencing in taxonomy and diagnostics, an application in "Candida" pathogenic yeasts with a metagenomic perspective

Targeted sequencing analysis pipeline for species identification of human pathogenic fungi using long-read nanopore sequencing

Among molecular-based techniques for fungal identification, Sanger sequencing of the primary universal fungal DNA barcode, the internal transcribed spacer (ITS) region (ITS1, 5.8S, ITS2), is commonly used in clinical routine laboratories due to its simplicity, universality, efficacy, and affordability for fungal species identification. However, Sanger sequencing fails to identify mixed ITS sequences in the case of mixed infections. To overcome this limitation, different high-throughput sequencing technologies have been explored. The nanopore-based technology is now one of the most promising long-read sequencing technologies on the market as it has the potential to sequence the full-length ITS region in a single read. In this study, we established a workflow for species identification using the sequences of the entire ITS region generated by nanopore sequencing of both pure yeast isolates and mocked mixed species reads generated with different scenarios. The species used in this study included Candida albicans (n = 2), Candida tropicalis (n = 1), Nakaseomyces glabratus (formerly Candida glabrata) (n = 1), Trichosporon asahii (n = 2), Pichia kudriavzevii (formerly Candida krusei) (n = 1), and Cryptococcus neoformans (n = 1). Comparing various methods to generate the consensus sequence for fungal species identification, the results from this study indicate that read clustering using a modified version of the NanoCLUST pipeline is more sensitive than Canu or VSEARCH, as it classified species accurately with a lower abundance cluster of reads (3% abundance compared to 10% with VSEARCH). The modified NanoCLUST also reduced the number of classified clusters compared to VSEARCH, making the subsequent BLAST+ analysis faster. Subsampling of the datasets, which reduces the size of the datasets by approximately tenfold, did not significantly affect the identification results in terms of the identified species name, percent identity, query coverage, percentage of reads in the classified cluster, and the number of clusters. The ability of the method to distinguish mixed species within sub-populations of large datasets has the potential to aid computer analysis by reducing the required processing power. The herein presented new sequence analysis pipeline will facilitate better interpretation of fungal sequence data for species identification.

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.

Next-Generation Sequencing Applications for the Study of Fungal Pathogens

Next-generation sequencing (NGS) has become a widely used technology in biological research. NGS applications for clinical pathogen detection have become vital technologies. It is increasingly common to perform fast, accurate, and specific detection of clinical specimens using NGS. Pathogenic fungi with high virulence and drug resistance cause life-threatening clinical infections. NGS has had a significant biotechnological impact on detecting bacteria and viruses but is not equally applicable to fungi. There is a particularly urgent clinical need to use NGS to help identify fungi causing infections and prevent negative impacts. This review summarizes current research on NGS applications for fungi and offers a visual method of fungal detection. With the development of NGS and solutions for overcoming sequencing limitations, we suggest clinicians test specimens as soon as possible when encountering infections of unknown cause, suspected infections in vital organs, or rapidly progressive disease.

Plague reservoir species throughout the world

Plague has been known since ancient times as a re-emerging infectious disease, causing considerable socioeconomic burden in regional hotspots. To better understand the epidemiological cycle of the causative agent of the plague, its potential occurrence, and possible future dispersion, one must carefully consider the taxonomy, distribution, and ecological requirements of reservoir-species in relation either to natural or human-driven changes (e.g. climate change or urbanization). In recent years, the depth of knowledge on species taxonomy and species composition in different landscapes has undergone a dramatic expansion, driven by modern taxonomic methods such as synthetic surveys that take into consideration morphology, genetics, and the ecological setting of captured animals to establish their species identities. Here, we consider the recent taxonomic changes of the rodent species in known plague reservoirs and detail their distribution across the world, with a particular focus on those rodents considered to be keystone host species. A complete checklist of all known plague-infectable vertebrates living in plague foci is provided as a Supporting Information table.

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.

Population Dynamics and Yeast Diversity in Early Winemaking Stages without Sulfites Revealed by Three Complementary Approaches

Nowadays, the use of sulfur dioxide (SO2) during the winemaking process is a controversial societal issue. In order to reduce its use, various alternatives are emerging, in particular bioprotection by adding yeasts, with different impacts on yeast microbiota in early winemaking stages. In this study, quantitative-PCR and metabarcoding high-throughput sequencing (HTS) were combined with MALDI-TOF-MS to monitor yeast population dynamic and diversity in the early stages of red winemaking process without sulfites and with bioprotection by Torulaspora delbrueckii and Metschnikowia pulcherrima addition. By using standard procedures for yeast protein extraction and a laboratory-specific database of wine yeasts, identification at species level of 95% of the isolates was successfully achieved by MALDI-TOF-MS, thus confirming that it is a promising method for wine yeast identification. The different approaches confirmed the implantation and the niche occupation of bioprotection leading to the decrease of fungal communities (HTS) and Hanseniaspora uvarum cultivable population (MALDI-TOF MS). Yeast and fungi diversity was impacted by stage of maceration and, to a lesser extent, by bioprotection and SO2, resulting in a modification of the nature and abundance of the operational taxonomic units (OTUs) diversity.

DNA barcoding of phytopathogens for disease diagnostics and bio-surveillance

DNA barcoding has proven to be a versatile tool for plant disease diagnostics in the genomics era. As the mass parallel and next generation sequencing techniques gained importance, the role of specific barcodes came under immense scrutiny. Identification and accurate classification of phytopathogens need a universal approach which has been the main application area of the concept of barcode. The present review entails a detailed description of the present status of barcode application in plant disease diagnostics. A case study on the application of Internal Transcribed Spacer (ITS) as barcode for Aspergillus and Fusarium spp. sheds light on the requirement of other potential candidates as barcodes for accurate identification. The challenges faced while barcoding novel pathogens have also been discussed with a comprehensive outline of integrating more recent technologies like meta-barcoding and genome skimming for detecting plant pathogens.

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.

FoodOmics as a new frontier to reveal microbial community and metabolic processes occurring on table olives fermentation

Table olives are considered the most widespread fermented food in the Mediterranean area and their consumption is expanding all over the world. This fermented vegetable can be considered as a natural functional food thanks to their high nutritional value and high content of bioactive compounds that contribute to the health and well-being of consumers. The presence of bioactive compounds is strongly influenced by a complex microbial consortium, traditionally exploited through culture-dependent approaches. Recently, the rapid spread of omics technologies has represented an important challenge to better understand the function, the adaptation and the exploitation of microbial diversity in different complex ecosystems, such as table olives. This review provides an overview of the potentiality of omics technologies to in depth investigate the microbial composition and the metabolic processes that drive the table olives fermentation, affecting both sensorial profile and safety properties of the final product. Finally, the review points out the role of omics approaches to raise at higher sophisticated level the investigations on microbial, gene, protein, and metabolite, with huge potential for the integration of table olives composition with functional assessments.

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.

Lost in translation: Population genomics and long-read sequencing reveals relaxation of concerted evolution of the ribosomal DNA cistron

Concerted evolution of the ribosomal DNA array has been studied in numerous eukaryotic taxa, yet is still poorly understood. rDNA genes are repeated dozens to hundreds of times in the eukaryotic genome (Eickbush and Eickbush, 2007) and it is believed that these arrays are homogenized through concerted evolution (Zimmer et al., 1980; Dover, 1993) preventing the accumulation of intragenomic, and intraspecific, variation. However, numerous studies have reported rampant intragenomic and intraspecific variation in the rDNA array (Ganley and Kobayashi, 2011; Naidoo et al., 2013; Hughes and Petersen, 2001; Lindner and Banik, 2011; Li et al., 2013; Lindner et al., 2013; Hughes et al., 2018), contradicting our current understanding of concerted evolution. The internal transcribed spacers (ITS) of the rDNA cistron are the most commonly used DNA barcoding region in Fungi (Schoch et al., 2012), and rely on concerted evolution to homogenize the rDNA array leading to a "barcode gap" (Puillandre et al., 2012). Here we show that in Boletus edulis Bull., ITS intragenomic variation persists at low allele frequencies throughout the rDNA array, this variation does not correlate with genomic relatedness between populations, and rDNA genes may not evolve in a strictly concerted fashion despite the presence of unequal recombination and gene conversion. Under normal assumptions, heterozygous positions found in ITS sequences represent hybridization between populations, yet through allelic mapping of the rDNA array we found numerous heterozygous alleles to be stochastically introgressed throughout, presenting a dishonest signal of gene flow. Moreover, despite the signal of gene flow in ITS, our organisms were highly inbred, indicating a disconnect between true gene flow and barcoding signals. In addition, we show that while the mechanisms of concerted evolution are ongoing in pseudo-heterozygous individuals, they are not fully homogenizing the ITS array. Concerted evolution of the rDNA array may insufficiently homogenize the ITS gene, allowing for misleading signals of gene flow to persist, vastly complicating the use of the ITS locus for DNA barcoding in Fungi.

High-Throughput Rapid and Inexpensive Assay for Quantitative Determination of Low Cell-Density Yeast Cultures

A procedure for microbial cell density determination with a high-throughput densitometric assay was developed to allow a precise quantification of both free and sessile cells, such as those of a biofilm, with a large range from low to high cell densities. Densitometry was chosen because it allows fast, rapid and cost-effective measures; it is non-disruptive; and has an easy learning curve. The method setup, and the further validation, was carried out with strains of Candida albicans, C. tropicalis and C. parapsilosis. Equations were developed at the level of the single strains, of the three species and finally a general one applicable to all three species. In the cross validation, with strains absent from the training set, the method was shown to be robust and flexible. The best results were obtained with species specific equations, although the global equation performed almost as well in terms of correlation between real and estimated density values. In all cases, a correlation around 0.98 between effective and predicted density was obtained with figures ranging from 10² to 10⁸ cells mL⁻¹. The entire analytical part of the procedure can be accomplished with a MS Excel macro provided free of charge.

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

A yeast strain was isolated during a study on vineyard-associated yeast strains from Apulia in Southern Italy. ITS and LSU D1/D2 rDNA sequences showed this strain not to belong to any known species and was described as the type strain of Ogataea uvarum sp. nov., a close relative of O. philodendri. Several secondary peaks appeared in the sequences, suggesting internal heterogeneity among the copies of the rDNA. This hypothesis was tested by sequencing single clones of the marker region. The analyses showed different levels of variability throughout the operon with differences between the rRNA encoding genes and the internally transcribed regions. O. uvarum and O. philodendri share high frequency variants, i.e., variants frequently found in many clones, whereas there is a large variability of the low frequency polymorphisms, suggesting that the mechanism of homogenization is more active with the former than with the latter type of variation. These findings indicate that low frequency variants are detected in Sanger sequencing as secondary peaks whereas in Next Generation Sequencing (NGS) of metagenomics DNA would lead to an overestimate of the alpha diversity. For the first time in our knowledge, this investigation shed light on the variation of the copy number of the rDNA cistron during the yeast speciation process. These polymorphisms can be used to investigate on the processes occurring in these taxonomic markers during the separation of fungal species, it being a genetic process highly frequent in the complex microbial ecosystem existing in grape, must and wine.

Sequence-based nomenclature: a reply to Thines et al. and Zamora et al. and provisions for an amended proposal "from the floor" to allow DNA sequences as types of names

We reply to two recently published, multi-authored opinion papers by opponents of sequence-based nomenclature, namely Zamora et al. (IMA Fungus9: 167-175,2018) and Thines et al. (IMA Fungus9: 177-183, 2018). While we agree with some of the principal arguments brought forward by these authors, we address misconceptions and demonstrate that some of the presumed evidence presented in these papers has been wrongly interpreted. We disagree that allowing sequences as types would fundamentally alter the nature of types, since a similar nature of abstracted features as type is already allowed in the Code (Art. 40.5), namely an illustration. We also disagree that there is a high risk of introducing artifactual taxa, as this risk can be quantified at well below 5 %, considering the various types of high-throughput sequencing errors. Contrary to apparently widespread misconceptions, sequence-based nomenclature cannot be based on similarity-derived OTUs and their consensus sequences, but must be derived from rigorous, multiple alignment-based phylogenetic methods and quantitative, single-marker species recognition algorithms, using original sequence reads; it is therefore identical in its approach to single-marker studies based on physical types, an approach allowed by the Code. We recognize the limitations of the ITS as a single fungal barcoding marker, but point out that these result in a conservative approach, with "false negatives" surpassing "false positives"; a desirable feature of sequence-based nomenclature. Sequence-based nomenclature does not aim at accurately resolving species, but at naming sequences that represent unknown fungal lineages so that these can serve as a means of communication, so ending the untenable situation of an exponentially growing number of unlabeled fungal sequences that fill online repositories. The risks are outweighed by the gains obtained by a reference library of named sequences spanning the full array of fungal diversity. Finally, we elaborate provisions in addition to our original proposal to amend the Code that would take care of the issues brought forward by opponents to this approach. In particular, taking up the idea of the Candidatus status of invalid, provisional names in prokaryote nomenclature, we propose a compromise that would allow valid publication of voucherless, sequence-based names in a consistent manner, but with the obligate designation as "nom. seq." (nomen sequentiae). Such names would not have priority over specimen- or culture-based names unless either epitypified with a physical type or adopted for protection on the recommendation of a committee of the International Commission on the Taxonomy of Fungi following evaluation based on strict quality control of the underlying studies based on established rules or recommendations.