The D1/D2 domain of the large-subunit rDNA of the yeast species is unusually polymorphic

High intragenomic, intergenomic, and phenotypic diversity in pulcherrimin-producing Metschnikowia yeasts indicates a special mode of genome evolution

In molecular systematics, the delimitation of yeast species is based on the notion that the barcode differences are smaller within species than between them. The most widely used barcodes are segments of the chromosomal repeats coding for ribosomal RNAs that are homogenised in yeasts. The analysis of these segments of the type strains of ten species recently merged in Metschnikowia pulcherrima and 37 new isolates demonstrated that this is not the case in this species. The intragenomic diversity significantly exceeded the threshold gaps used to differentiate related yeast species. Large segments of the D1/D2 domains were not diverse within the genomes and could therefore be used to determine the taxonomic affiliation of the isolates. The genome structures of the isolates were compared by RAPD and the RFLP of the mitochondrial DNA. Both patterns were highly heterogeneous. The sequence analysis of the PUL4 gene (a member of the PUL gene cluster involved in pulcherrimin production) revealed very high intragenomic differences, suggesting that the genomes may be chimerised. Three phenotypic traits related to the antimicrobial antagonism characteristic of the species were also highly diverse and prone to reversible segregation resembling epigenetic processes (silencing and reactivation of regulators) rather than mutations and back-mutations. These features make M. pulcherrima unique among yeasts and indicate that it evolves in a non-standard way.

Galactomyces candidus diversity in the complex mycobiota of cow-milk bryndza cheese comprising antagonistic and sensitive strains

Galactomyces candidus (orthographic variant: Galactomyces candidum) is a heterogeneous species of Saccharomycetales that comprises dimorphic yeasts described previously under various names (e.g. Geotrichum, Dipodascus). Its strains are common components of the cheese surface mycobiota. This study identified genetically and physiologically heterogeneous G. candidus strains in the complex mycobiota of artisanal cow-milk bryndza samples. The traditional Slovak bryndza is a cheese produced from ewe's milk in cooler mountainous regions and from cow's milk in warmer low-land regions. The taxonomic analysis of the culturable yeasts of the latter version carried out in this study revealed considerable differences from the yeast biota previously described for ovine bryndza. However, the conventional D1/D2- and ITS-based barcode analyses could not assign unanimously all isolates to species because of the intragenomic barcode diversity in certain groups and the discordance between the D1/D2 and ITS results in other groups. The identified species and groups of isolates had different abilities to utilise the carbon and energy sources (lactose, lactate, lipids and proteins) available in milk and ripening cheese. The G. candidus strains did not metabolise lactose and lactate, hydrolysed milk proteins with diverse, usually moderate efficiency and only could grow on certain amino acids as only energy sources. Their preferred substrate was lipid. Under aerobic conditions, its hyphae penetrated the lipid droplets and degraded their content from inside by developing a dense internal mycelium. Sporulation and different MLST (multilocus sequence typing) patterns indicated that the Galactomyces strains could sexually interact and their genomes could recombine. The Galactomyces and Kluyveromyces isolates had antagonistic effects against other members of the mycobiota.

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.

Characterization of a new Blastobotrys navarrensis strain indicates that it is not a later synonym of Blastobotrys proliferans

The species Blastobotrys navarrensis Sesma and Ramirez was delineated based on the description of the single strain CBS 139.77T. Based on its phenotypic similarities to Blastobotrys proliferans, B. navarrensis CBS 139.77T was later considered a synonym of B. proliferans. In the present study, we isolated the yeast strain IST 508 (=PYCC 8784=CBS 16671) from the soil surrounding an olive tree in Ferreira do Alentejo, Portugal. The phylogenetic analysis of D1/D2 domain and ITS sequences from strain IST 508 indicates that is closely related to B. navarrensis and B. proliferans. Although strain IST 508 differs from B. navarrensis CBS 139.77T by 14 substitutions and 20 indels (6.6 % divergence) in the ITS sequence, no divergence was detected at the level of D1/D2 domain, mitochondrial small subunit rDNA, and cytochrome oxidase II sequences. On the other hand, strains IST 508 and CBS 139.77 differ from B. proliferans NRRL Y-17577T by eight substitutions (1.4 % divergence) in the D1/D2 domain sequence, by 16 substitutions (2.7 % divergence) in the cytochrome oxidase II sequence, and by 16 substitutions (3.7 % divergence) in the mitochondrial small subunit rDNA sequence. Due to the high number of variable phenotypic tests in B. proliferans and B. navarrensis, strains from the two species are difficult to distinguish. Contrasting with what is described for other Blastobotrys species, no differences were detected at the level of micromorphology between the two species. Nevertheless, based on the molecular differences between the two strains, CBS 139.77 and IST 508, and B. proliferans NRRL Y-17577T and their phylogenetic analysis, strains CBS 139.77 and IST 508 are from B. navarrensis and this species should be considered as an independent species and not a later synonym of B. proliferans. We propose an emended description of B. navarrensis.

Yeasts with Fermentative Potential Associated with Fruits of Camu-Camu (Myrciaria dubia, Kunth) from North of Brazilian Amazon

Considering the high biotechnological potential of yeasts associated to edible fruits, a screening for these microorganisms, capable of alcoholic fermentation, was performed in ripe fruits of camu-camu (Myrciaria dubia, Kunth). The fruits were collected from north of Brazilian Amazon, in the floodplain of the Cauamé River. Yeasts were isolated, and fermentation capability was evaluated using Durham tubes. Quantitative assays were performed to calculate ethanol yield (g g^-1), specific growth rate (h^-1), and ethanol productivity (g L^-1·h^-1). Taxonomic identification was performed by ribosomal gene nucleotide sequence analysis by alignment using BLASTN. A total of fifteen yeast colonies were isolated, and three of them presented promising ability to ferment glucose to ethanol. These isolates were identified as Candida orthopsilosis, Pichia kudriavzevii, and Meyerozyma caribbica. When cultured in broth containing 180 g·L^-1 of glucose, M. caribbica CC003 reached 91.7 percent of the maximum theoretical ethanol concentration (84.4 g·L^-1), presenting an ethanol yield and productivity of 0.4688 g·g^-1 and 0.781 g·L^-1·h^-1, respectively. These results indicate a promising potential of this isolate for bioprocess applications. This paper is a rare report of C. orthopsilosis with endophytic habit because most of the references indicate it as a human pathogen. Besides this, M. caribbica is a promising fermenter for alcoholic beverages due to its osmotolerance and high ethanol yield. This is the first paper reporting endophytic yeasts associated with fruits of Myrciaria dubia.

The evolving species concepts used for yeasts: from phenotypes and genomes to speciation networks

Here we review how evolving species concepts have been applied to understand yeast diversity. Initially, a phenotypic species concept was utilized taking into consideration morphological aspects of colonies and cells, and growth profiles. Later the biological species concept was added, which applied data from mating experiments. Biophysical measurements of DNA similarity between isolates were an early measure that became more broadly applied with the advent of sequencing technology, leading to a sequence-based species concept using comparisons of parts of the ribosomal DNA. At present phylogenetic species concepts that employ sequence data of rDNA and other genes are universally applied in fungal taxonomy, including yeasts, because various studies revealed a relatively good correlation between the biological species concept and sequence divergence. The application of genome information is becoming increasingly common, and we strongly recommend the use of complete, rather than draft genomes to improve our understanding of species and their genome and genetic dynamics. Complete genomes allow in-depth comparisons on the evolvability of genomes and, consequently, of the species to which they belong. Hybridization seems a relatively common phenomenon and has been observed in all major fungal lineages that contain yeasts. Note that hybrids may greatly differ in their post-hybridization development. Future in-depth studies, initially using some model species or complexes may shift the traditional species concept as isolated clusters of genetically compatible isolates to a cohesive speciation network in which such clusters are interconnected by genetic processes, such as hybridization.

Clavispora santaluciae f.a., sp. nov., a novel ascomycetous yeast species isolated from grapes

During a study of yeast diversity in Azorean vineyards, four strains were isolated which were found to represent a novel yeast species based on the sequences of the internal transcribed spacer (ITS) region (ITS1-5.8S-ITS2) and of the D1/D2 domain of the large subunit (LSU) rRNA gene, together with their physiological characteristics. An additional strain isolated from Drosophila suzukii in Italy had identical D1/D2 sequences and very similar ITS regions (five nucleotide substitutions) to the Azorean strains. Phylogenetic analysis using sequences of the ITS region and D1/D2 domain showed that the five strains are closely related to Clavispora lusitaniae, although with 56 nucleotide differences in the D2 domain. Intraspecies variation revealed between two and five nucleotide differences, considering the five strains of Clavispora santaluciae. Some phenotypic discrepancies support the separation of the new species from their closely related ones, such as the inability to grow at temperatures above 35 °C, to produce acetic acid and the capacity to assimilate starch. Neither conjugations nor ascospore formation were observed in any of the strains. The name Clavispora santaluciae f.a., sp. nov., is proposed to accommodate the above noted five strains (holotype, CBS 16465^T; MycoBank no., MB 835794).

Candida metrosideri pro tempore sp. nov. and Candida ohialehuae pro tempore sp. nov., two antifungal-resistant yeasts associated with Metrosideros polymorpha flowers in Hawaii

Flowers produce an array of nutrient-rich exudates in which microbes can thrive, making them hotspots for microbial abundance and diversity. During a diversity study of yeasts inhabiting the flowers of Metrosideros polymorpha (Myrtaceae) in the Hawai’i Volcanoes National Park (HI, USA), five isolates were found to represent two novel species. Morphological and physiological characterization, and sequence analysis of the small subunit ribosomal RNA (rRNA) genes, the D1/D2 domains of the large subunit rRNA genes, the internal transcribed spacer (ITS) regions, and the genes encoding the largest and second largest subunits of the RNA polymerase II (RPB1 and RPB2, respectively), classified both species in the family Metschnikowiaceae, and we propose the names Candida metrosideri pro tempore sp. nov. (JK22^T = CBS 16091 = MUCL 57821) and Candida ohialehuae pro tempore sp. nov. (JK58.2^T = CBS 16092 = MUCL 57822) for such new taxa. Both novel Candida species form a well-supported subclade in the Metschnikowiaceae containing species associated with insects, flowers, and a few species of clinical importance. The ascosporic state of the novel species was not observed. The two novel yeast species showed elevated minimum inhibitory concentrations to the antifungal drug amphotericin B (>4 μg/mL). The ecology and phylogenetic relationships of C. metrosideri and C. ohialehuae are also discussed.

Isolation of a new Papiliotrema laurentii strain that displays capacity to achieve high lipid content from xylose

In this work, we isolated and selected oleaginous yeasts from rock field soils from two National Parks in Brazil (Caparaó and Serra dos Órgãos) with the potential to accumulate oil from xylose, the main pentose sugar found in lignocellulosic biomass. From the 126 isolates, two were selected based on their lipid contents. They were taxonomically identified as Papiliotrema laurentii (UFV-1 and UFV-2). Of the two, P. laurentii UFV-1 was selected as the best lipid producer. Under unoptimized conditions, lipid production by P. laurentii UFV-1 was higher in glucose than in xylose. To improve its lipid production from xylose, we applied response surface methodology (RSM) with a face-centered central composite design (CCF). We evaluated the effects of agitation rate, initial cell biomass (OD₆₀₀), carbon/nitrogen ratio (C/N ratio) and pH on lipid production. P. laurentii UFV-1 recorded the highest lipid content, 63.5% (w/w) of the cell dry mass, under the following conditions: C/N ratio = 100:1, pH value = 7.0, initial OD₆₀₀ = 0.8 and agitation = 300 rpm. Under these optimized conditions, biomass, lipid titer and volumetric lipid productivity were 9.31 g/L, 5.90 g/L and 0.082 g/L.h, respectively. Additionally, we determined the fatty acid composition of P. laurentii UFV-1 as follows: C14:0 (0.5%), C16:0 (28.4-29.4%), C16:1 (0.2%), C18:0 (9.5-11%), C18:1 (58.6-60.5%), and C20:0 (0.7-0.8%). Based on this composition, the predicted properties of biodiesel showed that P. laurentii UFV-1 oil is suitable for use as feedstock in biodiesel production.

Towards yeast taxogenomics: lessons from novel species descriptions based on complete genome sequences

In recent years, ‘multi-omic’ sciences have affected all aspects of fundamental and applied biological research. Yeast taxonomists, though somewhat timidly, have begun to incorporate complete genomic sequences into the description of novel taxa, taking advantage of these powerful data to calculate more reliable genetic distances, construct more robust phylogenies, correlate genotype with phenotype and even reveal cryptic sexual behaviors. However, the use of genomic data in formal yeast species descriptions is far from widespread. The present review examines published examples of genome-based species descriptions of yeasts, highlights relevant bioinformatic approaches, provides recommendations for new users and discusses some of the challenges facing the genome-based systematics of yeasts.

Metschnikowia pulcherrima and Related Pulcherrimin-Producing Yeasts: Fuzzy Species Boundaries and Complex Antimicrobial Antagonism

Yeasts affiliated with the Metschnikowia pulcherrima clade (subclade) of the large ascomycetous genus Metschnikowia frequently turn out to produce the characteristic maroon-red pulcherrimin when tested for pigment production and prove to exert antagonistic effects on many types of microorganisms. The determination of the exact taxonomic position of the strains is hampered by the shortage of distinctive morphological and physiological properties of the species of the clade and the lack of rDNA barcode gaps. The rDNA repeats of the type strains of the species are not homogenized and are assumed to evolve by a birth-and-death mechanism combined with reticulation. The taxonomic division is further hampered by the incomplete biological (reproductive) isolation of the species: certain type strains can be hybridized and genome sequencing revealed chimeric genome structures in certain strains that might have evolved from interspecies hybrids (alloploid genome duplication). Various mechanisms have been proposed for the antimicrobial antagonism. One is related to pulcherrimin production. The diffusible precursor of pulcherrimin, the pulcherriminic acid is secreted by the cells into the environment where it forms the insoluble pulcherrimin with the ferric ions. The lack of free iron caused by the immobilization of ferric ions inhibits the growth of many microorganisms. Recent results of research into the complexity of the taxonomic division of the pulcherrimin-producing Metschnikowia yeasts and the mechanism(s) underlying their antimicrobial antagonism are discussed in this review.

Hyphopichia lachancei, f.a., sp. nov., a yeast species from diverse origins

Three strains originating from insect frass in South Africa, yellow foxglove in Hungary and soil in France, were characterised phenotypically and by sequencing of the D1/D2 domain of the large subunit and the ITS1-5.8S-ITS2 (ITS)-region of the rRNA gene. The strains have identical D1/D2 domain sequences and only one strain shows a 1 bp indel in a 9 bp homopolymer A/T repeat within the ITS-region. Based on sequence analysis Hyphopichia burtonii is the closest related species. The investigated strains differ from the type strain of H. burtonii by 1.9% (9 substitutions and an indel) in the D1/D2 domain and by 23 substitutions and 21-22 indels in the ITS-region. Since the sequence variability is very low among the three strains and the sequence divergence with the closely related H. burtonii exceeds the level generally encountered between species we propose the new species Hyphopichia lachancei f.a., sp. nov. to accommodate the three novel strains. From H. burtonii the new species can be distinguished phenotypically by its inability to ferment cellobiose and by the formation of endospores (Holotype: CBS 5999^T; Isotype: NCAIM Y.02228^T; MycoBank no.: MB833616).

Isolation of xylose-assimilating yeasts and optimization of xylitol production by a new Meyerozyma guilliermondii strain

Production of xylitol from lignocellulosic biomass is of interest to modern biorefineries, because this biomass should be processed into a spectrum of chemicals (bio-based products) and not only energy. The isolation of new yeast strains capable of efficiently converting xylose into xylitol and withstanding inhibitors released from biomass hydrolysis can contribute to making its production feasible in biorefineries. Forty-three out of 128 yeast strains isolated from the gut of Passalidae beetles were capable of assimilating xylose as the sole carbon source. Meyerozyma guilliermondii UFV-1 was selected due to its ability to grow and ferment D-xylose in a synthetic medium. This yeast assimilated the broad range of sugars present in lignocellulosic biomass hydrolysates, such as xylose, raffinose, cellobiose, rhamnose, arabinose, and glucose. Its optimum growth conditions were pH 8.0 and a temperature of 30 °C. In concentrations of 0.07 mol/L acetic acid, 0.05 mol/L 5-hydroximethylfurfural, and 0.04 mol/L furfural, M. guilliermondii UFV-1 did not grow. Maximum xylitol production in aerobiosis and hypoxia were 51.88 and 27.73 g/L, respectively. Under aerobic condition, xylose concentration and agitation rate were the factors which were statistically significant, while only the agitation rate was significant in hypoxia. We fitted a response surface (RS) that estimated the best agitation rate (113.33 rpm) and xylose concentration (90 g/L) for maximum xylitol production in aerobiosis. Therefore, M. guilliermondii UFV-1 displays potential for being used for xylitol production in biorefineries.

Genome Sequence of the Yeast Clavispora lusitaniae Type Strain CBS 6936

Clavispora lusitaniae, an environmental saprophytic yeast belonging to the CTG clade of Candida, can behave occasionally as an opportunistic pathogen in humans. We report here the genome sequence of the type strain CBS 6936. Comparison with sequences of strain ATCC 42720 indicates conservation of chromosomal structure but significant nucleotide divergence.

Contrasting Patterns of rDNA Homogenization within the Zygosaccharomyces rouxii Species Complex

Arrays of repetitive ribosomal DNA (rDNA) sequences are generally expected to evolve as a coherent family, where repeats within such a family are more similar to each other than to orthologs in related species. The continuous homogenization of repeats within individual genomes is a recombination process termed concerted evolution. Here, we investigated the extent and the direction of concerted evolution in 43 yeast strains of the Zygosaccharomyces rouxii species complex (Z. rouxii, Z. sapae, Z. mellis), by analyzing two portions of the 35S rDNA cistron, namely the D1/D2 domains at the 5’ end of the 26S rRNA gene and the segment including the internal transcribed spacers (ITS) 1 and 2 (ITS regions). We demonstrate that intra-genomic rDNA sequence variation is unusually frequent in this clade and that rDNA arrays in single genomes consist of an intermixing of Z. rouxii, Z. sapae and Z. mellis-like sequences, putatively evolved by reticulate evolutionary events that involved repeated hybridization between lineages. The levels and distribution of sequence polymorphisms vary across rDNA repeats in different individuals, reflecting four patterns of rDNA evolution: I) rDNA repeats that are homogeneous within a genome but are chimeras derived from two parental lineages via recombination: Z. rouxii in the ITS region and Z. sapae in the D1/D2 region; II) intra-genomic rDNA repeats that retain polymorphisms only in ITS regions; III) rDNA repeats that vary only in their D1/D2 domains; IV) heterogeneous rDNA arrays that have both polymorphic ITS and D1/D2 regions. We argue that an ongoing process of homogenization following allodiplodization or incomplete lineage sorting gave rise to divergent evolutionary trajectories in different strains, depending upon temporal, structural and functional constraints. We discuss the consequences of these findings for Zygosaccharomyces species delineation and, more in general, for yeast barcoding.

Intragenomic polymorphism and intergenomic recombination in the ribosomal RNA genes of strains belonging to a yeast species Pichia membranifaciens

A concerted evolution model has been proposed to explain the observed lack of sequence variation among the multiple ribosomal RNA (rRNA) gene copies in many different eukaryotic species. Recent studies on the level of intragenomic variations in the rRNA gene repeats of fungi resulted in controversial conclusions. In this study, we clearly showed that significant polymorphisms of the internal transcribed spacers (ITS1 and ITS2) of ribosomal DNA (rDNA) exist within the genome of a strain of the yeast species Pichia membranifaciens. More interestingly, we showed that the intragenomic ITS sequence polymorphisms were formed by intergenomic rDNA recombination among different P. membranifaciens strains with significantly different ITS sequences. Intergenomic rDNA recombination was also responsible for the diversification of rDNA sequences in different strains of the species. After the events bring together different rDNA types in individual genomes of the P. membranifaciens strains compared, rDNA sequence heterogeneity has remained in the genome of one but eliminated by homogenisation in the genomes of other strains. Our findings show new clue for further investigation on the mechanism of concerted evolution of rRNA genes in eukaryotes.

Advances in yeast systematics and phylogeny and their use as predictors of biotechnologically important metabolic pathways

Detection, identification and classification of yeasts have undergone a major transformation in the last decade and a half following application of gene sequence analyses and genome comparisons. Development of a database (barcode) of easily determined DNA sequences from domains 1 and 2 (D1/D2) of the nuclear large subunit rRNA gene and from ITS now permits many laboratories to identify species quickly and accurately, thus replacing the laborious and often inaccurate phenotypic tests previously used. Phylogenetic analysis of gene sequences is leading to a major revision of yeast systematics that will result in redefinition of nearly all genera. This new understanding of species relationships has prompted a change of rules for naming and classifying yeasts and other fungi, and these new rules are presented in the recently implemented International Code of Nomenclature for algae, fungi, and plants (Melbourne Code). The use of molecular methods for species identification and the impact of Code changes on classification will be discussed, as will use of phylogeny for prediction of biotechnological applications.

International Society of Human and Animal Mycology (ISHAM)-ITS reference DNA barcoding database--the quality controlled standard tool for routine identification of human and animal pathogenic fungi

Human and animal fungal pathogens are a growing threat worldwide leading to emerging infections and creating new risks for established ones. There is a growing need for a rapid and accurate identification of pathogens to enable early diagnosis and targeted antifungal therapy. Morphological and biochemical identification methods are time-consuming and require trained experts. Alternatively, molecular methods, such as DNA barcoding, a powerful and easy tool for rapid monophasic identification, offer a practical approach for species identification and less demanding in terms of taxonomical expertise. However, its wide-spread use is still limited by a lack of quality-controlled reference databases and the evolving recognition and definition of new fungal species/complexes. An international consortium of medical mycology laboratories was formed aiming to establish a quality controlled ITS database under the umbrella of the ISHAM working group on "DNA barcoding of human and animal pathogenic fungi." A new database, containing 2800 ITS sequences representing 421 fungal species, providing the medical community with a freely accessible tool at http://www.isham.org/ and http://its.mycologylab.org/ to rapidly and reliably identify most agents of mycoses, was established. The generated sequences included in the new database were used to evaluate the variation and overall utility of the ITS region for the identification of pathogenic fungi at intra-and interspecies level. The average intraspecies variation ranged from 0 to 2.25%. This highlighted selected pathogenic fungal species, such as the dermatophytes and emerging yeast, for which additional molecular methods/genetic markers are required for their reliable identification from clinical and veterinary specimens.

Genetic diversity of Clavispora lusitaniae isolated from Agave fourcroydes Lem, as revealed by DNA fingerprinting

This study characterized Clavispora lusitaniae strains isolated from different stages of the processing and early fermentation of a henequen (Agave fourcroydes) spirit produced in Yucatan, Mexico using a molecular technique. Sixteen strains identified based on morphological features, obtained from different substrates, were typed molecularly. Nine different versions of the divergent D1/D2 domain of the large-subunit ribosomal DNA sequence were identified among the C. lusitaniae strains. The greatest degree of polymorphism was found in the 90-bp structural motif of the D2 domain. The MSP-PCR technique was able to differentiate 100% of the isolates. This study provides significant insight into the genetic diversity of the mycobiota present during the henequen fermentation process, especially that of C. lusitaniae, for which only a few studies in plants have been published. The applied MSP-PCR markers were very efficient in revealing olymorphisms between isolates of this species.

Use of gene sequence analyses and genome comparisons for yeast systematics

Detection, identification and classification of yeasts have undergone a major transformation in the past decade and a half following application of gene sequence analyses and genome comparisons. Development of a database (barcode) of easily determined gene sequences from domains 1 and 2 (D1/D2) of large subunit rRNA and from the internal transcribed spacer (ITS) now permits many laboratories to identify species accurately and this has led to a doubling in the number of known species of yeasts over the past decade. Phylogenetic analysis of gene sequences has resulted in major revision of yeast systematics, resulting in redefinition of nearly all genera. Future work calls for application of genomics to refine our understanding of the species concept and to provide a better understanding of the boundaries of genera and higher levels of classification. This increased understanding of phylogeny is expected to allow prediction of the genetic potential of various clades and species for biotechnological applications and adaptation to environmental changes.

Metschnikowia laotica f.a., sp. nov., a dimorphic, pigment-producing yeast species isolated from fruit

Eight strains with identical sequences of the D1/D2 domains of the large subunit rRNA genes were isolated from fallen fruits in two distant localities in Laos. These strains represent a novel dimorphic budding yeast species producing invasive pseudohyphae and a brown pigment when growing on media containing quinic acid as the sole carbon source or tryptophan as the sole nitrogen source. Phylogenetic analysis of the sequences of the D1/D2 domains, the internal transcribed spacer (ITS) regions and the 18S rRNA genes placed the novel species in the Metschnikowia clade close to Candida torresii, Metschnikowia drosophilae and Candida danieliae. The taxonomic name Metschnikowia laotica f.a., sp. nov., reflecting the geographical origin of the isolates, is proposed for the novel species. The type strain is 11-524(T) ( = CBS 12961(T) = NCAIM Y.02124(T) = CCY 64-4-1(T)). The Mycobank number is MB 807383.

Yarrowia porcina sp. nov. and Yarrowia bubula f.a. sp. nov., two yeast species from meat and river sediment

Eleven yeast strains representing two hitherto undescribed species were isolated from different kinds of meat samples in Hungary and one from the sediment of a tropical freshwater river in Southeastern Brazil. The analysis of the sequences of their large subunit rRNA gene D1/D2 domain and the internal transcribed spacer (ITS) regions placed the two new species in the Yarrowia clade. Some of the seven strains representing the first new species can mate and give rise to asci and form ascospores embedded in capsular material, which qualifies it as the third teleomorph species of the Yarrowia clade. The name Yarrowia porcina sp. nov. (type strain: NCAIM Y.02100(T) = CBS 12935(T) = NRRL Y-63669(T), allotype strain UFMG-RD131(A) = CBS 12932(A)) is proposed for this new yeast species, which, based on physiological characters, is indistinguishable from Yarrowia lipolytica and some other species of the genus. Considerable intraspecific variability was detected among the sequences of the large subunit rRNA gene D1/D2 domains of the seven strains. The variability among the D1/D2 sequences exceeded the divergence observed among the ITS sequences and in some cases more than 1 % substitution among the D1/D2 sequences was detected. The conspecificity of these strains was supported by the low (0-3 substitutions) sequence divergence among their ITS sequences, the result of a parsimony network analysis utilizing the concatenated ITS and D1/D2 sequences and also by the fingerprint patterns generated by microsatellite primed PCR. No ascospore formation was observed in the group of the other five strains representing the second new species. These strains shared identical D1/D2 and ITS sequences. Yarrowia bubula f.a., sp. nov. (type strain: NCAIM Y.01998(T) = CBS 12934(T) = NRRL Y-63668(T)) is proposed to accommodate these strains.

Comparative analyses of classical phenotypic method and ribosomal RNA gene sequencing for identification of medically relevant Candida species

As the distribution of Candida species and their susceptibility to antifungal agents have changed, a new means of accurately and rapidly identifying these species is necessary for the successful early resolution of infection and the subsequent reduction of morbidity and mortality. The current work aimed to evaluate ribosomal RNA gene sequencing for the identification of medically relevant Candida species in comparison with a standard phenotypic method. Eighteen reference strains (RSs), 69 phenotypically identified isolates and 20 inconclusively identified isolates were examined. Internal transcribed spaces (ITSs) and D1/D2 of the 26S ribosomal RNA gene regions were used as targets for sequencing. Additionally, the sequences of the ITS regions were used to establish evolutionary relationships. The sequencing of the ITS regions was successful for 88% (94/107) of the RS and isolates, whereas 100% of the remaining 12% (13/107) of the samples were successfully analysed by sequencing the D1/D2 region. Similarly, genotypic analysis identified all of the RS and isolates, including the 20 isolates that were not phenotypically identified. Phenotypic analysis, however, misidentified 10% (7/69) of the isolates. Phylogenetic analysis allowed the confirmation of the relationships between evolutionarily close species. Currently, the use of genotypic methods is necessary for the correct identification of Candida species.

Metschnikowia Species Share a Pool of Diverse rRNA Genes Differing in Regions That Determine Hairpin-Loop Structures and Evolve by Reticulation

Modern taxonomy of yeasts is mainly based on phylogenetic analysis of conserved DNA and protein sequences. By far the most frequently used sequences are those of the repeats of the chromosomal rDNA array. It is generally accepted that the rDNA repeats of a genome have identical sequences due to the phenomenon of sequence homogenisation and can thus be used for identification and barcoding of species. Here we show that the rDNA arrays of the type strains of Metschnikowia andauensis and M. fructicola are not homogenised. Both have arrays consisting of diverse repeats that differ from each other in the D1/D2 domains by up to 18 and 25 substitutions. The variable sites are concentrated in two regions that correspond to back-folding stretches of hairpin loops in the predicted secondary structure of the RNA molecules. The substitutions do not alter significantly the overall hairpin-loop structure due to wobble base pairing at sites of C-T transitions and compensatory mutations in the complementary strand of the hairpin stem. The phylogenetic and network analyses of the cloned sequences revealed that the repeats had not evolved in a vertical tree-like way but reticulation might have shaped the rDNA arrays of both strains. The neighbour-net analysis of all cloned sequences of the type strains and the database sequences of different strains further showed that these species share a continuous pool of diverse repeats that appear to evolve by reticulate evolution.

Yamadazyma siamensis sp. nov., Yamadazyma phyllophila sp. nov. and Yamadazyma paraphyllophila sp. nov., three novel yeast species isolated from phylloplane in Thailand and Taiwan

Four strains representing three novel anamorphic yeast species were isolated from the external surface of sugarcane leaves (DMKU-RK254(T)), corn leaves (DMKU-RK548(T)), bean leaves (K129) in Thailand and hengchun pencilwood leaves (TrB1-1(T)) in Taiwan. On the basis of morphological, biochemical, physiological and chemotaxonomic characteristics, the sequence analysis of the D1/D2 region of the large subunit (LSU) rRNA gene, the internal transcribed spacer (ITS) region, the actin gene (ACT1) and the elongation factor 2 gene (EF2), the four strains were determined to represent novel Yamadazyma species although formation of ascospores was not observed. Strain DMKU-RK254(T) was determined to be related to Candida diddensiae, Candida naeodendra and Candida kanchanaburiensis but with 1.8, 1.8 and 2.0 % nucleotide substitutions in the D1/D2 region of the LSU rRNA gene, respectively. It was assigned to Yamadazyma siamensis sp. nov. (type strain DMKU-RK254(T) = BCC 50730(T) = NBRC 108901(T) = CBS 12573(T)). The sequences of the D1/D2 region of the LSU rRNA gene, the ITS region, ACT1 gene and EF2 gene of two strains (DMKU-RK548(T) and K129) were identical but differed from that of strain TrB1-1(T) by 0.6, 1.0, 3.3 and 5.9 % nucleotide substitutions, respectively. Therefore, the two strains (DMKU-RK548(T) and K129) and strain TrB1-1(T) were assigned to be two separate species. The closest species in terms of pairwise sequences similarity of the D1/D2 region to the two novel species was Yamadazyma philogaea but with 1.1-1.7 % nucleotide substitutions. The two strains (DMKU-RK548(T) and K129) were assigned to Yamadazyma phyllophila sp. nov. (type strain DMKU-RK548(T) = BCC 50736(T) = NBRC 108906(T) = CBS 12572(T)) and the strain TrB1-1(T) was named Yamadazyma paraphyllophila sp. nov. (type strain TrB1-1(T) = BCRC 23030(T) = CCTCC AY 204005(T) = CBS 9928(T)).

Candida kuoi sp. nov., an anamorphic species of the Starmerella yeast clade that synthesizes sophorolipids

A novel strain of anamorphic yeast, designated strain NRRL Y-27208T, was isolated from concentrated grape juice in Cape Province, South Africa. Analysis of nuclear large subunit rRNA gene sequences from the D1/D2 domains separated the novel isolate from strains of Starmerella bombicola and Starmerella meliponinorum, as well as from species of the genus Candida that are members of the Starmerella clade. Compared to previously described species, strain NRRL Y-27208T is most closely related to S. bombicola but can be separated from this species by its ability to grow on d-ribose and erythritol. Strain NRRL Y-27208T produced sophorolipids that have an open chain structure similar to Candida batistae, Candida riodocensis and Candida stellata, which is in contrast to the closed chain sophorolipids produced by S. bombicola and Candida apicola. The analyses showed that NRRL Y-27208T ( = CBS 7267T) represents a novel species distinct from previously described species, for which the name Candida kuoi sp. nov. is proposed.

Species reassignment of Geotrichum bryndzae, Geotrichum phurueaensis, Geotrichum silvicola and Geotrichum vulgare based on phylogenetic analyses and mating compatibility

The present classification of Galactomyces and its anamorph, Geotrichum, is based on various studies that used morphology, ecology, biochemistry, DNA-DNA reassociation comparisons and gene sequencing. In this study, the identities of strains of the Centraalbureau voor Schimmelcultures yeast culture collection, as well as seven strains from South Africa, were examined by analyses of the nucleotide divergence in the internal transcribed spacer regions of the nuclear rRNA gene (nrRNA) operon, the D1/D2 domains of the 26S rRNA gene and partial actin gene sequences as well as compatibility studies. The South African strains were assigned to species in the genus Galactomyces. The phylogenetic analyses and mating studies revealed that Geotrichum silvicola and Geotrichum bryndzae are synonyms of Galactomyces candidus and that Geotrichum vulgare is a synonym of Galactomyces pseudocandidus.

Identification of a newly isolated erythritol-producing yeast and cloning of its erythrose [corrected] reductase genes

A new erythritol-producing yeast (strain BH010) was isolated in this study. Analysis of the D1/D2 domain of the 26S rDNA sequence, the ITS/5.8S rDNA sequence [corrected] and the 18S rDNA sequence allowed the taxonomic position of strain BH010 to be discussed, [corrected] and it was identified and named Moniliella sp. BH010. Physiological characteristics were described. Scanning electron micrography clearly indicated that the cells were cylindrical to elliptical with an average size of 5 × 10 μm when growing in liquid medium [corrected] and that pseudohyphae and blastoconidia were observed when cultivated in agar plates. The erythrose [corrected] reductase genes were cloned, sequenced, and analyzed. BLAST analysis and multiple sequence alignment demonstrated that erythrose [corrected] reductase genes of Moniliella sp. BH010 shared very high homology with that of Trichosporonoides megachiliensis SNG-42 except for the presence of introns. The deduced amino acid sequences showed high homology to the aldo-keto reductase superfamily.

Flowers as islands: spatial distribution of nectar-inhabiting microfungi among plants of Mimulus aurantiacus, a hummingbird-pollinated shrub

Microfungi that inhabit floral nectar offer unique opportunities for the study of microbial distribution and the role that dispersal limitation may play in generating distribution patterns. Flowers are well-replicated habitat islands, among which the microbes disperse via pollinators. This metapopulation system allows for investigation of microbial distribution at multiple spatial scales. We examined the distribution of the yeast, Metschnikowia reukaufii, and other fungal species found in the floral nectar of the sticky monkey flower, Mimulus aurantiacus, a hummingbird-pollinated shrub, at a California site. We found that the frequency of nectar-inhabiting microfungi on a given host plant was not significantly correlated with light availability, nectar volume, or the percent cover of M. aurantiacus around the plant, but was significantly correlated with the location of the host plant and loosely correlated with the density of flowers on the plant. These results suggest that dispersal limitation caused by spatially nonrandom foraging by pollinators may be a primary factor driving the observed distribution pattern.

The airway microbiota in cystic fibrosis: a complex fungal and bacterial community--implications for therapeutic management

Background

Given the polymicrobial nature of pulmonary infections in patients with cystic fibrosis (CF), it is essential to enhance our knowledge on the composition of the microbial community to improve patient management. In this study, we developed a pyrosequencing approach to extensively explore the diversity and dynamics of fungal and prokaryotic populations in CF lower airways.

Methodology and Principal Findings

Fungi and bacteria diversity in eight sputum samples collected from four adult CF patients was investigated using conventional microbiological culturing and high-throughput pyrosequencing approach targeting the ITS2 locus and the 16S rDNA gene. The unveiled microbial community structure was compared to the clinical profile of the CF patients. Pyrosequencing confirmed recently reported bacterial diversity and observed complex fungal communities, in which more than 60% of the species or genera were not detected by cultures. Strikingly, the diversity and species richness of fungal and bacterial communities was significantly lower in patients with decreased lung function and poor clinical status. Values of Chao1 richness estimator were statistically correlated with values of the Shwachman-Kulczycki score, body mass index, forced vital capacity, and forced expiratory volume in 1 s (p = 0.046, 0.047, 0.004, and 0.001, respectively for fungal Chao1 indices, and p = 0.010, 0.047, 0.002, and 0.0003, respectively for bacterial Chao1 values). Phylogenetic analysis showed high molecular diversities at the sub-species level for the main fungal and bacterial taxa identified in the present study. Anaerobes were isolated with Pseudomonas aeruginosa, which was more likely to be observed in association with Candida albicans than with Aspergillus fumigatus.

Conclusions

In light of the recent concept of CF lung microbiota, we viewed the microbial community as a unique pathogenic entity. We thus interpreted our results to highlight the potential interactions between microorganisms and the role of fungi in the context of improving survival in CF.

First case report of bloodstream infection due to a Candida species closely related to the novel species Candida pseudorugosa

Candida pseudorugosa is a novel species closely related to Candida rugosa for which only one case has been reported. We report the first case of a bloodstream infection in humans caused by a Candida sp. closely related to C. pseudorugosa. We contribute evidence to show this organism as a potential human pathogen that may be misidentified by conventional methods, also pointing out its lower sensitivity to azoles and other antifungal agents.

NEW METHOD FOR YEAST IDENTIFICATION USING BURROWS–WHEELER TRANSFORM

The explosive growth in biological data in recent years has led to the development of new methods to identify DNA sequences. Many algorithms have recently been developed that search DNA sequences looking for unique DNA sequences. This paper considers the application of the Burrows–Wheeler transform (BWT) to the problem of unique DNA sequence identification. The BWT transforms a block of data into a format that is extremely well suited for compression. This paper presents a time-efficient algorithm to search for unique DNA sequences in a set of genes. This algorithm is applicable to the identification of yeast species and other DNA sequence sets.

Ribosomal DNA polymorphisms in the yeast Geotrichum candidum

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

Five novel Wickerhamomyces- and Metschnikowia-related yeast species, Wickerhamomyces chaumierensis sp. nov., Candida pseudoflosculorum sp. nov., Candida danieliae sp. nov., Candida robnettiae sp. nov. and Candida eppingiae sp. nov., isolated from plants

On the basis of nucleotide divergences in the D1/D2 domain of the 26S rRNA gene and the internal transcribed spacers (ITS) domain of the rRNA gene, five novel yeast species, Wickerhamomyces chaumierensis sp. nov. (CBS 8565T  = JCM 17246T), Candida pseudoflosculorum sp. nov. (CBS 8584T  = JCM 17242T), Candida danieliae sp. nov. (CBS 8533T  = JCM 17247T), Candida robnettiae sp. nov. (CBS 8580T  = JCM 17243T) and Candida eppingiae sp. nov. (CBS 8586T  = JCM 17241T), isolated from plants in Thailand and Guyana, are proposed in this study.

Ribosomal DNA sequence polymorphism and the delineation of two ascosporic yeast species: Metschnikowia agaves and Starmerella bombicola

The relationship between mating success and sequence divergence in the internal transcribed spacer (ITS)/5.8S-D1/D2 rDNA region was examined in isolates tentatively assigned to Metschnikowia agaves and Starmerella bombicola. Both species are haplontic and heterothallic, such that the formation of mature asci can be used as a measure of genetic compatibility. Parsimony haplotype network analysis and mating success confirmed that all known isolates of M. agaves are conspecific. The previously reported D1/D2 polymorphism of five substitutions was not corroborated; the maximum divergence observed between any two strains was three substitutions, four with ITS. Of 39 putative S. bombicola strains, 36 formed an ITS-D1/D2 haplotype network using the 95% criterion. Thirty-five strains could mate with one or more compatible partner. The excluded strains did not mate. Mature asci arose from crosses between individuals differing by as many as five, but not six or seven substitutions in the D1/D2 domain. All strains capable of mating formed mature asci with at least one partner and all network members could be linked to another member by three or fewer substitutions. These results support the use of sequence divergence as a criterion for species delineation, but caution against describing poorly sampled species solely on the basis of that criterion.

Intragenomic variation in the ITS rDNA region obscures phylogenetic relationships and inflates estimates of operational taxonomic units in genus Laetiporus

Regions of rDNA are commonly used to infer phylogenetic relationships among fungal species and as DNA barcodes for identification. These regions occur in large tandem arrays, and concerted evolution is believed to reduce intragenomic variation among copies within these arrays, although some variation still might exist. Phylogenetic studies typically use consensus sequencing, which effectively conceals most intragenomic variation, but cloned sequences containing intragenomic variation are becoming prevalent in DNA databases. To understand effects of using cloned rDNA sequences in phylogenetic analyses we amplified and cloned the ITS region from pure cultures of six Laetiporus species and one Wolfiporia species (Basidiomycota, Polyporales). An average of 66 clones were selected randomly and sequenced from 21 cultures, producing a total of 1399 interpretable sequences. Significant variation (≥ 5% variation in sequence similarity) was observed among ITS copies within six cultures from three species clades (L. cincinnatus, L. sp. clade J, and Wolfiporia dilatohypha) and phylogenetic analyses with the cloned sequences produced different trees relative to analyses with consensus sequences. Cloned sequences from L. cincinnatus fell into more than one species clade and numerous cloned L. cincinnatus sequences fell into entirely new clades, which if analyzed on their own most likely would be recognized as "undescribed" or "novel" taxa. The use of a 95% cut off for defining operational taxonomic units (OTUs) produced seven Laetiporus OTUs with consensus ITS sequences and 20 OTUs with cloned ITS sequences. The use of cloned rDNA sequences might be problematic in fungal phylogenetic analyses, as well as in fungal bar-coding initiatives and efforts to detect fungal pathogens in environmental samples.

Species richness of yeast communities in floral nectar of southern Spanish plants

Floral nectar of insect-pollinated plants often contains dense yeast populations, yet little quantitative information exists on patterns and magnitude of species richness of nectar-dwelling yeasts in natural plant communities. This study evaluates yeast species richness at both the plant community and plant species levels in a montane forest area in southern Spain, and also explores possible correlations between the incidence of different yeast species in nectar and their reported tolerance to high sugar concentrations, and between yeast diversity and pollinator composition. Yeast species occurring in a total of 128 field-collected nectar samples from 24 plant species were identified by sequencing the D1/D2 domain of the large subunit rDNA, and rarefaction-based analyses were used to estimate yeast species richness at the plant community and plant species levels, using nectar drops as elemental sampling units. Individual nectar samples were generally characterized by very low species richness (1.2 yeast species/sample, on average), with the ascomycetous Metschnikowia reukaufii and Metschnikowia gruessii accounting altogether for 84.7% of the 216 isolates identified. Other yeasts recorded included species in the genera Aureobasidium, Rhodotorula, Cryptococcus, Sporobolomyces, and Lecythophora. The shapes and slopes of observed richness accumulation curves were quite similar for the nectar drop and plant species approaches, but the two approaches yielded different expected richness estimates. Expected richness was higher for plant species-based than for nectar drop-based analyses, showing that the coverage of nectar yeast species occurring in the region would be improved by sampling additional host plant species. A significant correlation was found between incidence of yeast species in nectar and their reported ability to grow in a medium containing 50% glucose. Neither diversity nor incidence of yeasts was correlated with pollinator composition across plant species.

Lachancea mirantina sp. nov., an ascomycetous yeast isolated from the cachaca fermentation process

In the present work, a novel ascomycete species, Lachancea mirantina sp. nov., isolated from the fermentation process that produces cachaça, a Brazilian spirit, is proposed. Nucleotide sequence analysis of the 26S D1/D2 rDNA locus showed that L. mirantina sp. nov. was genetically related to Lachancea cidri and Lachancea fermentati, although some physiological traits showed remarkable differences. Analysis of the D1/D2 large-subunit rDNA molecular marker showed a clear distinction among all three species, confirming that L. mirantina sp. nov. belongs to a separate taxonomic species in the Lachancea clade. The type strain of Lachancea mirantina sp. nov. is URM 5925(T) ( = CLIB 1160(T)  = CBS 11717(T)).