Dadaist2: A Toolkit to Automate and Simplify Statistical Analysis and Plotting of Metabarcoding Experiments.

The taxonomic composition of microbial communities can be assessed using universal marker amplicon sequencing. The most common taxonomic markers are the 16S rDNA for bacterial communities and the internal transcribed spacer (ITS) region for fungal communities, but various other markers are used for barcoding eukaryotes. A crucial step in the bioinformatic analysis of amplicon sequences is the identification of representative sequences. This can be achieved using a clustering approach or by denoising raw sequencing reads. DADA2 is a widely adopted algorithm, released as an R library, that denoises marker-specific amplicons from next-generation sequencing and produces a set of representative sequences referred to as 'Amplicon Sequence Variants' (ASV). Here, we present Dadaist2, a modular pipeline, providing a complete suite for the analysis that ranges from raw sequencing reads to the statistics of numerical ecology. Dadaist2 implements a new approach that is specifically optimised for amplicons with variable lengths, such as the fungal ITS. The pipeline focuses on streamlining the data flow from the command line to R, with multiple options for statistical analysis and plotting, both interactive and automatic.

Diverse Lineages of Candida albicans Live on Old Oaks.

The human pathogen Candida albicans is considered an obligate commensal of animals, yet it is occasionally isolated from trees, shrubs, and grass. We generated genome sequence data for three strains of C. albicans that we isolated from oak trees in an ancient wood pasture, and compared these to the genomes of over 200 clinical strains. C. albicans strains from oak are similar to clinical C. albicans in that they are predominantly diploid and can become homozygous at the mating locus through whole-chromosome loss of heterozygosity. Oak strains differed from clinical strains in showing slightly higher levels of heterozygosity genome-wide. Using phylogenomic analyses and in silico chromosome painting, we show that each oak strain is more closely related to strains from humans and other animals than to strains from other oaks. The high genetic diversity of C. albicans from old oaks shows that they can live in this environment for extended periods of time.

Saccharomyces jurei sp. nov., isolation and genetic identification of a novel yeast species from Quercus robur.

Two strains, D5088T and D5095, representing a novel yeast species belonging to the genus Saccharomyces were isolated from oak tree bark and surrounding soil located at an altitude of 1000 m above sea level in Saint Auban, France. Sequence analyses of the internal transcribed spacer (ITS) region and 26S rRNA D1/D2 domains indicated that the two strains were most closely related to Saccharomyces mikatae and Saccharomyces paradoxus. Genetic hybridization analyses showed that both strains are reproductively isolated from all other Saccharomyces species and, therefore, represent a distinct biological species. The species name Saccharomyces jurei sp. nov. is proposed to accommodate these two strains, with D5088T (=CBS 14759T=NCYC 3947T) designated as the type strain.

Apiotrichum terrigenum sp. nov., a soil-associated yeast found in both the UK and mainland Europe.

Five arthroconidium-producing yeast strains representing a novel Trichosporon-like species were independently isolated from the UK, Hungary and Norway. Two strains (Bio4T and Bio21) were isolated from biogas reactors used for processing grass silage, with a third strain (S8) was isolated from soil collected at the same UK site. Two additional strains were isolated in mainland Europe, one from soil in Norway (NCAIM Y.02175) and the other from sewage in Hungary (NCAIM Y.02176). Sequence analyses of the D1/D2 domains of the LSU rRNA gene and internal transcribed spacer (ITS) region indicated that the novel species belongs to the recently reinstated genus Apiotrichum and is most closely related to Apiotrichum scarabaeorum, a beetle-associated species first found in South Africa. Despite having similar physiological characteristics, the two species can be readily distinguished from one another by ITS sequencing. The species name Apiotrichum terrigenum sp. nov. is proposed to accommodate these strains, with Bio4T (=CBS 11373T=NCYC 3540T) designated as the type strain. The Mycobank deposit number is MB817431.

Prevalence and Dynamics of Ribosomal DNA Micro-heterogeneity Are Linked to Population History in Two Contrasting Yeast Species.

Despite the considerable number and taxonomic breadth of past and current genome sequencing projects, many of which necessarily encompass the ribosomal DNA, detailed information on the prevalence and evolutionary significance of sequence variation in this ubiquitous genomic region are severely lacking. Here, we attempt to address this issue in two closely related yet contrasting yeast species, the baker's yeast Saccharomyces cerevisiae and the wild yeast Saccharomyces paradoxus. By drawing on existing datasets from the Saccharomyces Genome Resequencing Project, we identify a rich seam of ribosomal DNA sequence variation, characterising 1,068 and 970 polymorphisms in 34 S. cerevisiae and 26 S. paradoxus strains respectively. We discover the two species sets exhibit distinct mutational profiles. Furthermore, we show for the first time that unresolved rDNA sequence variation resulting from imperfect concerted evolution of the ribosomal DNA region follows a U-shaped allele frequency distribution in each species, similar to loci that evolve under non-concerted mechanisms but arising through rather different evolutionary processes. Finally, we link differences between the shapes of these allele frequency distributions to the two species' contrasting population histories.

Identification of furfural resistant strains of Saccharomyces cerevisiae and Saccharomyces paradoxus from a collection of environmental and industrial isolates.

BACKGROUND: Fermentation of bioethanol using lignocellulosic biomass as a raw material provides a sustainable alternative to current biofuel production methods by utilising waste food streams as raw material. Before lignocellulose can be fermented, it requires physical, chemical and enzymatic treatment in order to release monosaccharides, a process that causes the chemical transformation of glucose and xylose into the cyclic aldehydes furfural and hydroxyfurfural. These furan compounds are potent inhibitors of Saccharomyces fermentation, and consequently furfural tolerant strains of Saccharomyces are required for lignocellulosic fermentation. RESULTS: This study investigated yeast tolerance to furfural and hydroxyfurfural using a collection of 71 environmental and industrial isolates of the baker's yeast Saccharomyces cerevisiae and its closest relative Saccharomyces paradoxus. The Saccharomyces strains were initially screened for growth on media containing 100 mM glucose and 1.5 mg ml(-1) furfural. Five strains were identified that showed a significant tolerance to growth in the presence of furfural, and these were then screened for growth and ethanol production in the presence of increasing amounts (0.1 to 4 mg ml(-1)) of furfural. CONCLUSIONS: Of the five furfural tolerant strains, S. cerevisiae National Collection of Yeast Cultures (NCYC) 3451 displayed the greatest furfural resistance and was able to grow in the presence of up to 3.0 mg ml(-1) furfural. Furthermore, ethanol production in this strain did not appear to be inhibited by furfural, with the highest ethanol yield observed at 3.0 mg ml(-1) furfural. Although furfural resistance was not found to be a trait specific to any one particular lineage or population, three of the strains were isolated from environments where they might be continually exposed to low levels of furfural through the ongoing natural degradation of lignocelluloses, and would therefore develop elevated levels of resistance to these furan compounds. Thus, these strains represent good candidates for future studies of genetic variation relevant to understanding and manipulating furfural resistance and in the development of tolerant ethanologenic yeast strains for use in bioethanol production from lignocellulose processing.

Kazachstania yasuniensis sp. nov., an ascomycetous yeast species found in mainland Ecuador and on the Galápagos.

Seven strains representing a novel yeast species belonging to the genus Kazachstania were found at several collection sites on both mainland Ecuador (Yasuní National Park) and the Galápagos (Santa Cruz Island). Two strains (CLQCA 20-132(T) and CLQCA 24SC-045) were isolated from rotten wood samples, two further strains (CLQCA 20-280 and CLQCA 20-348) were isolated from soil samples, and three strains (CLQCA 20-198, CLQCA 20-374 and CLQCA 20-431) were isolated from decaying fruits. Sequence analyses of the D1/D2 domains of the LSU rRNA gene and ribosomal internal transcribed spacer (ITS) region indicated that the novel species is most closely related to Kazachstania servazzii and Kazachstania unispora. Although the strains could not be distinguished from one another based upon their differing geographical origins, they could be differentiated according to their isolation source (fruit, soil or wood) by ITS sequencing. The species name Kazachstania yasuniensis sp. nov. is proposed to accommodate these strains, with CLQCA 20-132(T) ( = CBS 13946(T) = NCYC 4008(T)) designated the type strain.

Ribosomal DNA sequence heterogeneity reflects intraspecies phylogenies and predicts genome structure in two contrasting yeast species.

The ribosomal RNA encapsulates a wealth of evolutionary information, including genetic variation that can be used to discriminate between organisms at a wide range of taxonomic levels. For example, the prokaryotic 16S rDNA sequence is very widely used both in phylogenetic studies and as a marker in metagenomic surveys and the internal transcribed spacer region, frequently used in plant phylogenetics, is now recognized as a fungal DNA barcode. However, this widespread use does not escape criticism, principally due to issues such as difficulties in classification of paralogous versus orthologous rDNA units and intragenomic variation, both of which may be significant barriers to accurate phylogenetic inference. We recently analyzed data sets from the Saccharomyces Genome Resequencing Project, characterizing rDNA sequence variation within multiple strains of the baker's yeast Saccharomyces cerevisiae and its nearest wild relative Saccharomyces paradoxus in unprecedented detail. Notably, both species possess single locus rDNA systems. Here, we use these new variation datasets to assess whether a more detailed characterization of the rDNA locus can alleviate the second of these phylogenetic issues, sequence heterogeneity, while controlling for the first. We demonstrate that a strong phylogenetic signal exists within both datasets and illustrate how they can be used, with existing methodology, to estimate intraspecies phylogenies of yeast strains consistent with those derived from whole-genome approaches. We also describe the use of partial Single Nucleotide Polymorphisms, a type of sequence variation found only in repetitive genomic regions, in identifying key evolutionary features such as genome hybridization events and show their consistency with whole-genome Structure analyses. We conclude that our approach can transform rDNA sequence heterogeneity from a problem to a useful source of evolutionary information, enabling the estimation of highly accurate phylogenies of closely related organisms, and discuss how it could be extended to future studies of multilocus rDNA systems. [concerted evolution; genome hydridisation; phylogenetic analysis; ribosomal DNA; whole genome sequencing; yeast].

Wickerhamomyces arborarius f.a., sp. nov., an ascomycetous yeast species found in arboreal habitats on three different continents.

Five strains representing a novel yeast species belonging to the genus Wickerhamomyces were independently isolated from Ecuador, Taiwan and the USA. One strain (CLQCA 10-161(T)) was isolated from the white flower of an unidentified plant species collected in the Maquipucuna cloud forest reserve, near Quito, in Ecuador. A second strain (GY7L12) was isolated from the leaf of a Chinese sumac or nutgall tree (Rhus chinensis 'roxburghiana') collected in the Taoyuan mountain area, Kachsiung, in Taiwan. Three additional strains (A543, A546 and A563) were isolated from two species of wood-boring beetle (Xyleborus glabratus and Xyleborinus saxeseni) collected near Clyo, Georgia, USA. Analysis of the D1/D2 domains of the LSU rRNA gene indicated that the novel species belongs to the genus Wickerhamomyces, and is most closely related to Wickerhamomyces sydowiorum, an insect-associated species predominantly found in South Africa. The North American and Taiwanese strains have identical internal transcribed spacer (ITS) sequences and can be distinguished from the Ecuadorian strain based on a single nucleotide substitution in the ITS1 region. The species name of Wickerhamomyces arborarius f.a., sp. nov. is proposed to accommodate these strains, with CLQCA 10-161(T) ( = CBS 12941(T) = NCYC 3743(T)) designated the type strain.

Candida ecuadorensis sp. nov., an ascomycetous yeast species found in two separate regions of Ecuador.

In the course of an on-going study aimed at cataloguing the natural yeast biodiversity found in Ecuador, two strains (CLQCA 13-025 and CLQCA 20-004(T)) were isolated from samples of cow manure and rotten wood collected in two separate provinces of the country (Orellana and Bolívar). These strains were found to represent a novel yeast species based on the sequences of their D1/D2 domain of the large-subunit (LSU) rRNA gene and their physiological characteristics. Phylogenetic analysis based on LSU D1/D2 sequences revealed this novel species to belong to the Metschnikowia clade and to be most closely related to Candida suratensis, a species recently discovered in a mangrove forest in Thailand. The species name of Candida ecuadorensis sp. nov. is proposed to accommodate these strains, with strain CLQCA 20-004(T) (=CBS 12653(T) = NCYC 3782(T)) designated as the type strain.

Candida theae sp. nov., a new anamorphic beverage-associated member of the Lodderomyces clade.

Four strains representing a novel yeast species belonging to the genus Candida were independently isolated in Taiwan and Ecuador. Two strains (G17(T) and G31) were isolated in Taiwan, by pellet precipitation from plastic-bottled tea drinks produced in Indonesia, while two additional strains (CLQCA 10-049 and CLQCA 10-062) were recovered from ancient chicha fermentation vessels found in tombs in Quito, Ecuador. These four strains were morphologically, and phylogenetically identical to each other. No sexual reproduction was observed on common sporulation media. Large-subunit (LSU) rRNA gene sequence analysis revealed the four strains to belong to the Lodderomyces clade, closely related to members of the Candida parapsilosis species complex. The four strains, which have identical LSU D1/D2 sequences, differ from their closest phylogenetic neighbors, Candida orthopsilosis and Candida parapsilosis, by 6-9 nt substitutions, respectively. Physiologically, the four strains are similar to Candida parapsilosis, although they can be distinguished from their closest relative by the assimilation of arbutin, nitrite, and creatine. The Indonesian and Ecuadorian strain sets can also be distinguished from one another based on ITS sequencing, differing by 4 substitutions in ITS1 and 1 single nucleotide indel in ITS2. Collectively, the results indicate that the four strains represent a previously unrecognized species of Candida. The name Candida theae sp. nov. is proposed to accommodate these strains, with G-17(T) (BCRC 23242(T)=CBS 12239(T)=ATCC MYA-4746(T)) designated as the type strain.

Candida northwykensis sp. nov., a novel yeast isolated from the gut of the click beetle Melanotus villosus.

Two yeast morphotypes, BET 4(T) and BET 7, were isolated from the gut of click beetle Melanotus villosus. Click beetles were collected from the decaying timber within the woodlands of North Wyke Research, South West England, UK (latitude, 50°46'29″N; longitude, 3°55'23″W). Morphotype BET 7 was identified as Debaryomyces hansenii, and the other morphotype, BET 4(T), was found to differ from Priceomyces castillae and Priceomyces haplophilus, its closest phylogenetic neighbours, by 5.0% with respect to the nucleotide sequence of the D1/D2 domain of the large-subunit (LSU) rRNA gene, and by 8.0% with respect to the ribosomal internal-transcribed spacer (ITS) region. BET 4(T) also differ from P. castillae and P. haplophilus in a number of different phenotypic characteristics. Thus, based on the unique nucleotide sequences of its D1/D2 domain and ITS region, its physiological characteristics and an inability to sporulate, strain BET 4(T) is assigned the status of a new species of Candida, for which the name Candida northwykensis sp. nov., is proposed. The type strain is BET 4(T) (NCYC 3525(T) = CBS 11370(T)).

Saturnispora quitensis sp. nov., a yeast species isolated from the Maquipucuna cloud forest reserve in Ecuador.

A single strain, CLQCA-10-114(T), representing a novel yeast species belonging to the genus Saturnispora was isolated from the fruit of an unidentified species of bramble (Rubus sp.), collected from the Maquipucuna cloud forest reserve, near Quito, in Ecuador. Sequence analyses of the D1/D2 domains of the large-subunit rRNA gene and ribosomal internal transcribed spacer region indicated that the novel species is most closely related to the recently described species Saturnispora gosingensis, isolated from the fruiting body of a mushroom collected in Taiwan, and Saturnispora hagleri, a Drosophila-associated yeast found in Brazil. The name Saturnispora quitensis sp. nov. is proposed to accommodate this strain; the type strain is CLQCA-10-114(T) (=CBS 12184(T)=NCYC 3744(T)).

TURNIP: tracking unresolved nucleotide polymorphisms in large hard-to-assemble regions of repetitive DNA sequence.

UNLABELLED: TURNIP comprises a suite of Perl scripts and modules that facilitates the resolution of microheterogeneity within hard-to-assemble repetitive DNA sequences. TURNIP was originally developed for the Saccharomyces Genome Resequencing Project (SGRP) within which the ribosomal DNA (rDNA) of 36 strains of S.cerevisiae were analysed to investigate the occurrence of potential polymorphisms. Here, 'partially resolved SNPs', or pSNPs, as well as indels, were found to be far more prevalent than previously suspected. More generally, the TURNIP software ascertains degrees of variation between large tandem repeats within a single locus, offering insights into mechanisms of genome stability and gene conversion in any organism for which genome sequence data are available. AVAILABILITY: The TURNIP source code, results files and online help are available at http://www.ncyc.co.uk/software/turnip.html.

Cryptococcus shivajii sp. nov.: a novel basidiomycetous yeast isolated from biogas reactor.

Five yeast morphotypes were isolated from biogas reactors at North Wyke Research, Okehampton, UK. Out of the five morphotypes, four were identified as known species. In contrast, the fifth morphotype strain, Bio10(T), was found to differ from Bullera dendrophila and Kwoniella mangroviensis, its closest phylogenetic neighbours, by 2.6-2.9% with respect to the nucleotide sequence of the D1/D2 domain of the 26S rRNA gene and by 5.6-6.2% with respect to the internal transcribed spacer 1 (ITS1)-5.8S rRNA gene-ITS2 region. Bio10(T) also differs from these two species by a number of phenotypic characteristics. Thus, based on the phenotypic differences and phylogenetic analysis, strain Bio10(T) is assigned the status of a new species of Cryptococcus, for which the name Cryptococcus shivajii sp. nov. is proposed. The type strain is Bio10(T) (NCYC 3541(T) = CBS 11374(T)).

Candida carvajalis sp. nov., an ascomycetous yeast species from the Ecuadorian Amazon jungle.

In the course of a yeast biodiversity survey of different ecological habitats found in Ecuador, two yeast strains (CLQCA 20-011(T) and CLQCA20-014) were isolated from samples of rotten wood and fallen leaf debris collected at separate sites in the central region of the Ecuadorian Amazonia. These strains were found to represent a novel yeast species based on the sequences of their D1/D2 domain of the large-subunit (LSU) rRNA gene and their physiological characteristics. Phylogenetic analysis based on LSU D1/D2 sequences revealed this novel species to be most closely related to Candida asparagi, Candida fructus, Candida musae and two as yet undescribed Candida species, with the six yeast taxa collectively forming a distinct species group within the Clavispora clade. The species name of Candida carvajalis sp. nov. is proposed to accommodate these strains, with CLQCA 20-011(T) (NCYC 3509(T), CBS 11361(T)) designated as the type strain.

Population genomics of domestic and wild yeasts.

Since the completion of the genome sequence of Saccharomyces cerevisiae in 1996 (refs 1, 2), there has been a large increase in complete genome sequences, accompanied by great advances in our understanding of genome evolution. Although little is known about the natural and life histories of yeasts in the wild, there are an increasing number of studies looking at ecological and geographic distributions, population structure and sexual versus asexual reproduction. Less well understood at the whole genome level are the evolutionary processes acting within populations and species that lead to adaptation to different environments, phenotypic differences and reproductive isolation. Here we present one- to fourfold or more coverage of the genome sequences of over seventy isolates of the baker's yeast S. cerevisiae and its closest relative, Saccharomyces paradoxus. We examine variation in gene content, single nucleotide polymorphisms, nucleotide insertions and deletions, copy numbers and transposable elements. We find that phenotypic variation broadly correlates with global genome-wide phylogenetic relationships. S. paradoxus populations are well delineated along geographic boundaries, whereas the variation among worldwide S. cerevisiae isolates shows less differentiation and is comparable to a single S. paradoxus population. Rather than one or two domestication events leading to the extant baker's yeasts, the population structure of S. cerevisiae consists of a few well-defined, geographically isolated lineages and many different mosaics of these lineages, supporting the idea that human influence provided the opportunity for cross-breeding and production of new combinations of pre-existing variations.

Repetitive sequence variation and dynamics in the ribosomal DNA array of Saccharomyces cerevisiae as revealed by whole-genome resequencing.

Ribosomal DNA (rDNA) plays a key role in ribosome biogenesis, encoding genes for the structural RNA components of this important cellular organelle. These genes are vital for efficient functioning of the cellular protein synthesis machinery and as such are highly conserved and normally present in high copy numbers. In the baker's yeast Saccharomyces cerevisiae, there are more than 100 rDNA repeats located at a single locus on chromosome XII. Stability and sequence homogeneity of the rDNA array is essential for function, and this is achieved primarily by the mechanism of gene conversion. Detecting variation within these arrays is extremely problematic due to their large size and repetitive structure. In an attempt to address this, we have analyzed over 35 Mbp of rDNA sequence obtained from whole-genome shotgun sequencing (WGSS) of 34 strains of S. cerevisiae. Contrary to expectation, we find significant rDNA sequence variation exists within individual genomes. Many of the detected polymorphisms are not fully resolved. For this type of sequence variation, we introduce the term partial single nucleotide polymorphism, or pSNP. Comparative analysis of the complete data set reveals that different S. cerevisiae genomes possess different patterns of rDNA polymorphism, with much of the variation located within the rapidly evolving nontranscribed intergenic spacer (IGS) region. Furthermore, we find that strains known to have either structured or mosaic/hybrid genomes can be distinguished from one another based on rDNA pSNP number, indicating that pSNP dynamics may provide a reliable new measure of genome origin and stability.

Molecular evidence for the existence of natural hybrids in the genus Zygosaccharomyces.

26S rDNA D1/D2 sequencing was used to characterise a number of food-associated Zygosaccharomyces rouxii strains held at the National Collection of Yeast Cultures. In the course of this study, four strains (NCYC 1682, NCYC 3042, NCYC 3060 and NCYC 3061) were identified which appeared, based on their D1/D2 sequences, to belong to a novel Zygosaccharomyces species. However, subsequent sequence analysis showed that NCYC 1682, NCYC 3060 and NCYC 3061 possess two highly divergent copies of the nuclear-encoded ADE2, HIS3 and SOD2 genes, indicating these three strains are in fact hybrids. NCYC 3042, however, does appear to represent a novel species which may be hypothesized to have crossed with Z. rouxii and given rise to hybrid strains. Additional approaches to define precise taxonomic status and mechanisms of hybrid genome formation amongst yeast species are discussed.

Determining a unique defining DNA sequence for yeast species using hashing techniques.

MOTIVATION: Yeasts are often still identified with physiological growth tests, which are both time consuming and unsuitable for detection of a mixture of organisms. Hence, there is a need for molecular methods to identify yeast species. RESULTS: A hashing technique has been developed to search for unique DNA sequences in 702 26S rRNA genes. A unique DNA sequence has been found for almost every yeast species described to date. The locations of the unique defining sequences are in accordance with the variability map of large subunit ribosomal RNA and provide detail of the evolution of the D1/D2 region. This approach will be applicable to the rapid identification of unique sequences in other DNA sequence sets. AVAILABILITY: Freely available upon request from the authors. SUPPLEMENTARY INFORMATION: Results are available at http://www.sys.uea.ac.uk/~jjw/project/paper