Analysis of clinical Candida parapsilosis isolates reveals copy number variation in key fluconazole resistance genes.

We used whole-genome sequencing to analyze a collection of 35 fluconazole-resistant and 7 susceptible Candida parapsilosis isolates together with coverage analysis and GWAS techniques to identify new mechanisms of fluconazole resistance. Phylogenetic analysis shows that although the collection is diverse, two persistent clinical lineages were identified. We identified copy number variation (CNV) of two genes, ERG11 and CDR1B, in resistant isolates. Two strains have a CNV at the ERG11 locus; the entire ORF is amplified in one, and only the promoter region is amplified in the other. We show that the annotated telomeric gene CDR1B is actually an artifactual in silico fusion of two highly similar neighboring CDR genes due to an assembly error in the C. parapsilosis CDC317 reference genome. We report highly variable copy numbers of the CDR1B region across the collection. Several strains have increased the expansion of the two genes into a tandem array of new chimeric genes. Other strains have experienced a deletion between the two genes creating a single gene with a reciprocal chimerism. We find translocations, duplications, and gene conversion across the CDR gene family in the C. parapsilosis species complex, showing that it is a highly dynamic family.

WHO elements - A new category of selfish genetic elements at the borderline between homing elements and transposable elements.

Homing genetic elements are a form of selfish DNA that inserts into a specific target site in the genome and spreads through the population by a process of biased inheritance. Two well-known types of homing element, called inteins and homing introns, were discovered decades ago. In this review we describe WHO elements, a newly discovered type of homing element that constitutes a distinct third category but is rare, having been found only in a few yeast species so far. WHO elements are inferred to spread using the same molecular homing mechanism as inteins and introns: they encode a site-specific endonuclease that cleaves the genome at the target site, making a DNA break that is subsequently repaired by copying the element. For most WHO elements, the target site is in the glycolytic gene FBA1. WHO elements differ from inteins and homing introns in two fundamental ways: they do not interrupt their host gene (FBA1), and they occur in clusters. The clusters were formed by successive integrations of different WHO elements into the FBA1 locus, the result of an 'arms race' between the endonuclease and its target site. We also describe one family of WHO elements (WHO10) that is no longer specifically associated with the FBA1 locus and instead appears to have become transposable, inserting at random genomic sites in Torulaspora globosa with up to 26 copies per strain. The WHO family of elements is therefore at the borderline between homing genetic elements and transposable elements.

Genome Analysis of a Newly Discovered Yeast Species, Hanseniaspora menglaensis.

Annual surveys of Irish soil samples identified three isolates, CBS 16921 (UCD88), CBS 18246 (UCD443), and CBS 18247 (UCD483), of an apiculate yeast species within the Hanseniaspora genus. The internal transcribed spacer (ITS) and D1/D2 region of the large subunit (LSU) rRNA sequences showed that these are isolates of the recently described species Hanseniaspora menglaensis, first isolated from Southwest China. No genome sequence for H. menglaensis is currently available. The genome sequences of the three Irish isolates were determined using short-read (Illumina) sequencing, and the sequence of one isolate (CBS 16921) was assembled to chromosome level using long-read sequencing (Oxford Nanopore Technologies). Phylogenomic analysis shows that H. menglaensis belongs to the fast-evolving lineage (FEL) of Hanseniaspora. Only one MAT idiomorph (encoding MATα1) was identified in all three sequenced H. menglaensis isolates, consistent with one mating type of a heterothallic species. Genome comparisons showed that there has been a rearrangement near MATα of FEL species compared to isolates from the slowly evolving lineage (SEL).

Draft genome sequence of the yeast Schwanniomyces capriottii strain UCD805, isolated from soil in Ireland.

Schwanniomyces capriottii is a member of the Debaryomycetaceae family in the order Saccharomycetales. Here, we present the genome sequence of S. capriottii UCD805, which was isolated from soil in Dublin, Ireland. This genome is 12.2 Mb and was assembled into 14 scaffolds plus a mitochondrial genome scaffold.

Genome sequences of two isolates of the yeast Candida zeylanoides: UCD849 from soil in Ireland, and AWD from an African wild dog.

We report genome sequences of two new isolates of the budding yeast Candida zeylanoides. Strain UCD849 from soil in Ireland was assembled into eight complete chromosomes. Strain AWD from an African Wild Dog in a US zoo was sequenced to draft level. The assemblies are 10.6 Mb and 99.57% identical.

Genome sequence data of the strongly antagonistic yeast Pichia kluyveri isolate APC 11.10 B as a foundation for analysing biocontrol mechanisms.

Pichia kluyveri strain APC 11.10 B was isolated from apple bark in Switzerland and exhibited strong antagonistic activity against plant pathogenic fungi in vitro (e.g., Botrytis, Fusarium or Monilinia isolates). In order to identify the mechanisms underlying this antagonism, we have sequenced the genome of this isolate by long- and short-read sequencing technologies. The sequence data were de novo assembled into nine scaffolds and a fully resolved circularized mitogenome. The total genome size was 10.9 Mbp and 7451 potential open reading frames (ORFs) and 202 tRNA genes were predicted. In comparison to two P. kluyveri genomes deposited at the NCBI (of strains X31-10 and CBA6002), the APC 11.10 B strain seemed to represent a hybrid because backmapping of sequencing reads resulted in a high rate of heterozygous and structural variants in the nuclear genome (this was not observed for the mitochondrial genome). The P. kluyveri (APC 11.10 B) draft genome represents a first step and resource for genome mining, comparative and functional genomics (e.g., identifying the biocontrol mode of action), and evolutionary studies. Since the genus Pichia comprises many biotechnologically relevant yeasts, the genome data may be used in a variety of fields and disciplines.

Analysis of clinical Candida parapsilosis isolates reveals copy number variation in key fluconazole resistance genes.

We used whole-genome sequencing to analyse a collection of 35 fluconazole resistant and 7 susceptible Candida parapsilosis isolates together with coverage analysis and GWAS techniques to identify new mechanisms of fluconazole resistance. Phylogenetic analysis shows that although the collection is diverse, two probable outbreak groups were identified. We identified copy number variation of two genes, ERG11 and CDR1B, in resistant isolates. Two strains have a CNV at the ERG11 locus; the entire ORF is amplified in one, and only the promoter region is amplified in the other. We show the annotated telomeric gene CDR1B is actually an artefactual in silico fusion of two highly similar neighbouring CDR genes due to an assembly error in the C. parapsilosis CDC317 reference genome. We report highly variable copy numbers of the CDR1B region across the collection. Several strains have increased expansion of the two genes into a tandem array of new chimeric genes. Other strains have experienced a deletion between the two genes creating a single gene with a reciprocal chimerism. We find translocations, duplications, and gene conversion across the CDR gene family in the C. parapsilosis species complex, showing that it is a highly dynamic family.

Identification of European isolates of the lager yeast parent Saccharomyces eubayanus.

Lager brewing first occurred in Bavaria in the 15th century, associated with restrictions of brewing to colder months. The lager yeast, Saccharomyces pastorianus, is cold tolerant. It is a hybrid between Saccharomyces cerevisiae and Saccharomyces eubayanus, and has been found only in industrial settings. Natural isolates of S. eubayanus were first discovered in Patagonia 11 years ago. They have since been isolated from China, Tibet, New Zealand, and North America, but not from Europe. Here, we describe the first European strains UCD646 and UCD650, isolated from a wooded area on a university campus in Dublin, Ireland. We generated complete chromosome level assemblies of both genomes using long- and short-read sequencing. The UCD isolates belong to the Holarctic clade. Genome analysis shows that isolates similar to the Irish strains contributed to the S. eubayanus component of S. pastorianus, but isolates from Tibet made a larger contribution.

Identification of genetic variants of the industrial yeast Komagataella phaffii (Pichia pastoris) that contribute to increased yields of secreted heterologous proteins.

The yeast Komagataella phaffii (formerly called Pichia pastoris) is used widely as a host for secretion of heterologous proteins, but only a few isolates of this species exist and all the commonly used expression systems are derived from a single genetic background, CBS7435 (NRRL Y-11430). We hypothesized that other genetic backgrounds could harbor variants that affect yields of secreted proteins. We crossed CBS7435 with 2 other K. phaffii isolates and mapped quantitative trait loci (QTLs) for secretion of a heterologous protein, ß-glucosidase, by sequencing individual segregant genomes. A major QTL mapped to a frameshift mutation in the mannosyltransferase gene HOC1, which gives CBS7435 a weaker cell wall and higher protein secretion than the other isolates. Inactivation of HOC1 in the other isolates doubled ß-glucosidase secretion. A second QTL mapped to an amino acid substitution in IRA1 that tripled ß-glucosidase secretion in 1-week batch cultures but reduced cell viability, and its effects are specific to this heterologous protein. Our results demonstrate that QTL analysis is a powerful method for dissecting the basis of biotechnological traits in nonconventional yeasts, and a route to improving their industrial performance.

A widespread inversion polymorphism conserved among Saccharomyces species is caused by recurrent homogenization of a sporulation gene family.

Saccharomyces genomes are highly collinear and show relatively little structural variation, both within and between species of this yeast genus. We investigated the only common inversion polymorphism known in S. cerevisiae, which affects a 24-kb 'flip/flop' region containing 15 genes near the centromere of chromosome XIV. The region exists in two orientations, called reference (REF) and inverted (INV). Meiotic recombination in this region is suppressed in crosses between REF and INV orientation strains such as the BY x RM cross. We find that the inversion polymorphism is at least 17 million years old because it is conserved across the genus Saccharomyces. However, the REF and INV isomers are not ancient alleles but are continually being re-created by re-inversion of the region within each species. Inversion occurs due to continual homogenization of two almost identical 4-kb sequences that form an inverted repeat (IR) at the ends of the flip/flop region. The IR consists of two pairs of genes that are specifically and strongly expressed during the late stages of sporulation. We show that one of these gene pairs, YNL018C/YNL034W, codes for a protein that is essential for spore formation. YNL018C and YNL034W are the founder members of a gene family, Centroid, whose members in other Saccharomycetaceae species evolve fast, duplicate frequently, and are preferentially located close to centromeres. We tested the hypothesis that Centroid genes are a meiotic drive system, but found no support for this idea.

Draft Genome Sequence of the Yeast Blastobotrys aristata Strain UCD613, Isolated from Soil in Ireland.

Blastobotrys aristata is a member of the Trichomonascaceae family in the order Saccharomycetales. Here, we present the genome sequence of B. aristata UCD613, which was isolated from soil in Dublin, Ireland. This genome is 13.3 Mb and was assembled into 4 chromosome-size scaffolds of >2.2 Mb in size plus a mitochondrial genome scaffold.

Draft Genome Sequence of the Yeast Torulaspora quercuum Strain UCD657, Isolated from Soil in Ireland.

Torulaspora quercuum is an ascomycete yeast first isolated in 2009. Here, we present the genome sequence of T. quercuum isolate UCD657, which was isolated from soil in Ireland. This genome is 10.4 Mb and was assembled into 8 chromosome-sized scaffolds of >1 Mb in size, plus a mitochondrial genome scaffold.

Systematic Analysis of Copy Number Variations in the Pathogenic Yeast Candida parapsilosis Identifies a Gene Amplification in RTA3 That is Associated with Drug Resistance.

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 (CPAR2_104610) and ARR3 (CPAR2_601050), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3, which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3, 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700, coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes. IMPORTANCE Copy number variations (CNVs) are an important source of genomic diversity that have been associated with drug resistance. We identify two unusual CNVs in the human fungal pathogen Candida parapsilosis. Both target a single gene (RTA3 or ARR3), and they have occurred multiple times in multiple isolates. The copy number of RTA3, a putative floppase that controls the inward translocation of lipids in the cell membrane, correlates with resistance to miltefosine, a derivative of phosphatidylcholine (PC) that was originally developed as an anticancer drug. In 2021, miltefosine was designated an orphan drug by the United States Food and Drug Administration for the treatment of invasive candidiasis. Importantly, we find that resistance to miltefosine is also caused by mutations in flippases, which control the outward movement of lipids, and that many C. parapsilosis isolates are prone to easily acquiring an increased resistance to miltefosine.

Mating-Type Switching in Budding Yeasts, from Flip/Flop Inversion to Cassette Mechanisms.

Mating-type switching is a natural but unusual genetic control process that regulates cell identity in ascomycete yeasts. It involves physically replacing one small piece of genomic DNA by another, resulting in replacement of the master regulatory genes in the mating pathway and hence a switch of cell type and mating behavior. In this review, we concentrate on recent progress that has been made on understanding the origins and evolution of mating-type switching systems in budding yeasts (subphylum Saccharomycotina). Because of the unusual nature and the complexity of the mechanism in Saccharomyces cerevisiae, mating-type switching was assumed until recently to have originated only once or twice during yeast evolution. However, comparative genomics analysis now shows that switching mechanisms arose many times independently-at least 11 times in budding yeasts and once in fission yeasts-a dramatic example of convergent evolution. Most of these lineages switch mating types by a flip/flop mechanism that inverts a section of a chromosome and is simpler than the well-characterized 3-locus cassette mechanism (MAT/HML/HMR) used by S. cerevisiae. Mating-type switching (secondary homothallism) is one of the two possible mechanisms by which a yeast species can become self-fertile. The other mechanism (primary homothallism) has also emerged independently in multiple evolutionary lineages of budding yeasts, indicating that homothallism has been favored strongly by natural selection. Recent work shows that HO endonuclease, which makes the double-strand DNA break that initiates switching at the S. cerevisiae MAT locus, evolved from an unusual mobile genetic element that originally targeted a glycolytic gene, FBA1.

Genome sequence data of the antagonistic soil-borne yeast Cyberlindnera sargentensis (SHA 17.2).

Cyberlindnera sargentensis strain SHA 17.2, isolated from a Swiss soil sample, exhibited strong antagonistic activity against several plant pathogenic fungi in vitro and was highly competitive against other yeasts in soil. As a basis for identifying the mechanisms underlying its strong antagonistic activity, we have sequenced the genome of C. sargentensis (SHA 17.2) by long- and short read sequencing, de novo assembled them into seven contigs/chromosomes and a mitogenome (total genome size 11.4 Mbp), and annotated 5455 genes. This high-quality genome is the reference for transcriptome and proteome analyses aiming at elucidating the mode of action of C. sargentensis against fungal plant pathogens. It will thus serve as a resource for identifying potential biocontrol genes and performing comparative genomics analyses of yeast genomes.

Genomic diversity, chromosomal rearrangements, and interspecies hybridization in the Ogataea polymorpha species complex.

The methylotrophic yeast Ogataea polymorpha has long been a useful system for recombinant protein production, as well as a model system for methanol metabolism, peroxisome biogenesis, thermotolerance, and nitrate assimilation. It has more recently become an important model for the evolution of mating-type switching. Here, we present a population genomics analysis of 47 isolates within the O. polymorpha species complex, including representatives of the species O. polymorpha, Ogataea parapolymorpha, Ogataea haglerorum, and Ogataea angusta. We found low levels of nucleotide sequence diversity within the O. polymorpha species complex and identified chromosomal rearrangements both within and between species. In addition, we found that one isolate is an interspecies hybrid between O. polymorpha and O. parapolymorpha and present evidence for loss of heterozygosity following hybridization.

Draft Genome Sequence of the Yeast Ogataea degrootiae Strain UCD465, Isolated from Soil in Ireland.

Ogataea degrootiae is an ascomycete yeast that was first isolated in the Netherlands in 2017. It is a member of the Pichiaceae clade. Here, we present the genome sequence of O. degrootiae UCD465, which was isolated from soil in Ireland. This genome is 14.6 Mb and haploid.

Draft Genome Sequence of a Diploid and Hybrid Candida Strain, Candida sanyaensis UCD423, Isolated from Compost in Ireland.

Candida sanyaensis is a CUG-Ser1 clade yeast that is associated with soil. Assembly of short-read and long-read data shows that C. sanyaensis has a diploid and hybrid genome, with approximately 97% identity between the haplotypes. The haploid genome size is approximately 15.4 Mb.

Giant GAL gene clusters for the melibiose-galactose pathway in Torulaspora.

In many yeast species, the three genes at the centre of the galactose catabolism pathway, GAL1, GAL10 and GAL7, are neighbours in the genome and form a metabolic gene cluster. We report here that some yeast strains in the genus Torulaspora have much larger GAL clusters that include genes for melibiase (MEL1), galactose permease (GAL2), glucose transporter (HGT1), phosphoglucomutase (PGM1) and the transcription factor GAL4, in addition to GAL1, GAL10, and GAL7. Together, these eight genes encode almost all the steps in the pathway for catabolism of extracellular melibiose (a disaccharide of galactose and glucose). We show that a progenitor 5-gene cluster containing GAL 7-1-10-4-2 was likely present in the common ancestor of Torulaspora and Zygotorulaspora. It added PGM1 and MEL1 in the ancestor of most Torulaspora species. It underwent further expansion in the T. pretoriensis clade, involving the fusion of three progenitor clusters in tandem and the gain of HGT1. These giant GAL clusters are highly polymorphic in structure, and subject to horizontal transfers, pseudogenization and gene losses. We identify recent horizontal transfers of complete GAL clusters from T. franciscae into one strain of T. delbrueckii, and from a relative of T. maleeae into one strain of T. globosa. The variability and dynamic evolution of GAL clusters in Torulaspora indicates that there is strong natural selection on the GAL pathway in this genus.