Identification of a novel Candida metapsilosis isolate reveals multiple hybridization events.

Candida metapsilosis is a member of the Candida parapsilosis species complex, a group of opportunistic human pathogens. Of all the members of this complex, C. metapsilosis is the least virulent, and accounts for a small proportion of invasive Candida infections. Previous studies established that all C. metapsilosis isolates are hybrids, originating from a single hybridization event between two lineages, parent A and parent B. Here, we use MinION and Illumina sequencing to characterize a C. metapsilosis isolate that originated from a separate hybridization. One of the parents of the new isolate is very closely related to parent A. However, the other parent (parent C) is not the same as parent B. Unlike C. metapsilosis AB isolates, the C. metapsilosis AC isolate has not undergone introgression at the mating type-like locus. In addition, the A and C haplotypes are not fully collinear. The C. metapsilosis AC isolate has undergone loss of heterozygosity with a preference for haplotype A, indicating that this isolate is in the early stages of genome stabilization.

Population genomics of the pathogenic yeast Candida tropicalis identifies hybrid isolates in environmental samples.

Candida tropicalis is a human pathogen that primarily infects the immunocompromised. Whereas the genome of one isolate, C. tropicalis MYA-3404, was originally sequenced in 2009, there have been no large-scale, multi-isolate studies of the genetic and phenotypic diversity of this species. Here, we used whole genome sequencing and phenotyping to characterize 77 isolates of C. tropicalis from clinical and environmental sources from a variety of locations. We show that most C. tropicalis isolates are diploids with approximately 2-6 heterozygous variants per kilobase. The genomes are relatively stable, with few aneuploidies. However, we identified one highly homozygous isolate and six isolates of C. tropicalis with much higher heterozygosity levels ranging from 36-49 heterozygous variants per kilobase. Our analyses show that the heterozygous isolates represent two different hybrid lineages, where the hybrids share one parent (A) with most other C. tropicalis isolates, but the second parent (B or C) differs by at least 4% at the genome level. Four of the sequenced isolates descend from an AB hybridization, and two from an AC hybridization. The hybrids are MTLa/α heterozygotes. Hybridization, or mating, between different parents is therefore common in the evolutionary history of C. tropicalis. The new hybrids were predominantly found in environmental niches, including from soil. Hybridization is therefore unlikely to be associated with virulence. In addition, we used genotype-phenotype correlation and CRISPR-Cas9 editing to identify a genome variant that results in the inability of one isolate to utilize certain branched-chain amino acids as a sole nitrogen source.

Polymorphic centromere locations in the pathogenic yeast Candida parapsilosis.

Centromeres pose an evolutionary paradox: strongly conserved in function but rapidly changing in sequence and structure. However, in the absence of damage, centromere locations are usually conserved within a species. We report here that isolates of the pathogenic yeast species Candida parapsilosis show within-species polymorphism for the location of centromeres on two of its eight chromosomes. Its old centromeres have an inverted-repeat (IR) structure, whereas its new centromeres have no obvious structural features but are located within 30 kb of the old site. Centromeres can therefore move naturally from one chromosomal site to another, apparently spontaneously and in the absence of any significant changes in DNA sequence. Our observations are consistent with a model in which all centromeres are genetically determined, such as by the presence of short or long IRs or by the ability to form cruciforms. We also find that centromeres have been hotspots for genomic rearrangements in the C. parapsilosis clade.

Spatial inter-centromeric interactions facilitated the emergence of evolutionary new centromeres.

Centromeres of Candida albicans form on unique and different DNA sequences but a closely related species, Candida tropicalis, possesses homogenized inverted repeat (HIR)-associated centromeres. To investigate the mechanism of centromere type transition, we improved the fragmented genome assembly and constructed a chromosome-level genome assembly of C. tropicalis by employing PacBio sequencing, chromosome conformation capture sequencing (3C-seq), chromoblot, and genetic analysis of engineered aneuploid strains. Further, we analyzed the 3D genome organization using 3C-seq data, which revealed spatial proximity among the centromeres as well as telomeres of seven chromosomes in C. tropicalis. Intriguingly, we observed evidence of inter-centromeric translocations in the common ancestor of C. albicans and C. tropicalis. Identification of putative centromeres in closely related Candida sojae, Candida viswanathii and Candida parapsilosis indicates loss of ancestral HIR-associated centromeres and establishment of evolutionary new centromeres (ENCs) in C. albicans. We propose that spatial proximity of the homologous centromere DNA sequences facilitated karyotype rearrangements and centromere type transitions in human pathogenic yeasts of the CUG-Ser1 clade.

Draft Genome Sequence of the Yeast Kazachstania telluris CBS 16338 Isolated from Forest Soil in Ireland.

Draft genomes of several Kazachstania species, a large group of ascomycetous budding yeasts, have been recently published. However, there is none yet available for pathogenic species from the K. telluris clade, including K. telluris sensu stricto, an opportunistic yeast that has been isolated from various niches and human clinical samples. Here we provide the first draft genome sequence of the strain K. telluris CBS 16338, that we isolated from forest soil in Ireland.

Draft Genome Sequence of a Red Basidiomycete Yeast, Symmetrospora coprosmae Strain UCD350, Isolated from Soil in Ireland.

Symmetrospora coprosmae is a red yeast from the subphylum Pucciniomycotina in the phylum Basidiomycota. Here, we present the first genome sequence of S. coprosmae strain UCD350, from an isolate collected from soil in Ireland. The genome size is 20.2 Mb.

Draft Genome Sequences of Two Isolates of the Yeast Kazachstania servazzii Recovered from Soil in Ireland.

We sequenced two isolates of Kazachstania servazzii, UCD13 and UCD335, from soil in Ireland. Heterozygosity in these diploid genomes differs 19-fold between the two strains. Most currently available K. servazzii genome sequences come from Korean kimchi isolates, so our data will facilitate analysis of diversity in this species.

Draft Genome Sequence of the Birch Tree Fungal Pathogen Taphrina betulina UCD315.

Taphrina betulina is the ascomycete yeast that causes the formation of witches' brooms in birch trees. Here, we report the first draft genome sequence of T. betulina, from strain UCD315, isolated from soil in Ireland. The genome is haploid and 12.5 Mb long.

Coverage-Versus-Length Plots, a Simple Quality Control Step for de Novo Yeast Genome Sequence Assemblies.

Illumina sequencing has revolutionized yeast genomics, with prices for commercial draft genome sequencing now below $200. The popular SPAdes assembler makes it simple to generate a de novo genome assembly for any yeast species. However, whereas making genome assemblies has become routine, understanding what they contain is still challenging. Here, we show how graphing the information that SPAdes provides about the length and coverage of each scaffold can be used to investigate the nature of an assembly, and to diagnose possible problems. Scaffolds derived from mitochondrial DNA, ribosomal DNA, and yeast plasmids can be identified by their high coverage. Contaminating data, such as cross-contamination from other samples in a multiplex sequencing run, can be identified by its low coverage. Scaffolds derived from the bacteriophage PhiX174 and Lambda DNAs that are frequently used as molecular standards in Illumina protocols can also be detected. Assemblies of yeast genomes with high heterozygosity, such as interspecies hybrids, often contain two types of scaffold: regions of the genome where the two alleles assembled into two separate scaffolds and each has a coverage level C, and regions where the two alleles co-assembled (collapsed) into a single scaffold that has a coverage level 2C Visualizing the data with Coverage-vs.-Length (CVL) plots, which can be done using Microsoft Excel or Google Sheets, provides a simple method to understand the structure of a genome assembly and detect aberrant scaffolds or contigs. We provide a Python script that allows assemblies to be filtered to remove contaminants identified in CVL plots.

Draft Genome Sequences of Two Natural Isolates of the Yeast Barnettozyma californica from Ireland, UCD09 and UCD89.

No genome sequence of a species from Barnettozyma, a yeast genus in the family Phaffomycetaceae, is currently available. We isolated two B. californica strains from soils in Ireland and generated draft sequences of their 11.7-Mb genomes. Single nucleotide polymorphism (SNP) analysis showed 20,490 differences between the strains and suggests that B. californica is haploid.

Draft Genome Sequence of the Yeast Nadsonia starkeyi-henricii UCD142, Isolated from Forest Soil in Ireland.

We report a draft genome sequence of a strain of the nonfermentative yeast Nadsonia starkeyi-henricii, isolated from soil in a forest in Ireland. Comparison to Nadsonia fulvescens shows few rearrangements and a level of divergence similar to that of Saccharomyces cerevisiae versus Saccharomyces paradoxus Its mitochondrial genome lacks NAD genes.

Draft Genome Sequence of a Highly Heterozygous Yeast Strain from the Metschnikowia pulcherrima Subclade, UCD127.

Metschnikowia strain UCD127 was isolated from soil in Ireland and sequenced. It is a highly heterozygous diploid strain with 385,000 single nucleotide polymorphisms (SNPs). Its ribosomal DNA has the highest similarity to that of M. chrysoperlae, but its ACT1 and TEF1 loci and mitochondrial genome show affinity to those of M. fructicola, whose genome is significantly larger.

Genome analysis of the yeast Diutina catenulata, a member of the Debaryomycetaceae/Metschnikowiaceae (CTG-Ser) clade.

Diutina catenulata (Candida catenulata) is an ascomycetous yeast that has been isolated from humans, animals and environmental sources. The species is a contaminant of dairy products, and has been linked to superficial and invasive infections in both humans and animals. Previous phylogenetic analyses have assigned the species to the Saccharomycetales, but failed to identify its specific clade. Here, we report the genome sequence of an environmental isolate of D. catenulata. Examination of the tRNA repertoire and coding potential of this species shows that it translates the CUG codon as serine and not leucine. In addition, two phylogenetic analyses using 204 ubiquitous gene family alignments and 3,826 single-copy genes both confirm the placement of the species in the Debaryomycetaceae/Metschnikowiaceae, or CTG-Ser clade. The sequenced isolate contains an MTLα idiomorph. However, unlike most MTL loci in related species, poly (A) polymerase (PAP) is not adjacent to MTLα1.