Flow Cytometry Analysis of Fungal Ploidy

Spindle architecture constrains karyotype evolution

The eukaryotic cell division machinery must rapidly and reproducibly duplicate and partition the cell's chromosomes in a carefully coordinated process. However, chromosome numbers vary dramatically between genomes, even on short evolutionary timescales. We sought to understand how the mitotic machinery senses and responds to karyotypic changes by using a series of budding yeast strains in which the native chromosomes have been successively fused. Using a combination of cell biological profiling, genetic engineering and experimental evolution, we show that chromosome fusions are well tolerated up until a critical point. Cells with fewer than five centromeres lack the necessary number of kinetochore-microtubule attachments needed to counter outward forces in the metaphase spindle, triggering the spindle assembly checkpoint and prolonging metaphase. Our findings demonstrate that spindle architecture is a constraining factor for karyotype evolution.

Chromosome-level genome assembly of an auxotrophic strain of the pathogenic yeast Candida parapsilosis

We report the genome sequence of the pathogenic yeast Candida parapsilosis strain SR23 (CBS 7157) used in a number of experimental studies. The nuclear genome assembly consists of eight chromosome-sized contigs with a total size of 13.04 Mbp (N50 2.09 Mbp) and a G+C content of 38.7%.

Erg251 has complex and pleiotropic effects on sterol composition, azole susceptibility, filamentation, and stress response phenotypes

Ergosterol is essential for fungal cell membrane integrity and growth, and numerous antifungal drugs target ergosterol. Inactivation or modification of ergosterol biosynthetic genes can lead to changes in antifungal drug susceptibility, filamentation and stress response. Here, we found that the ergosterol biosynthesis gene ERG251 is a hotspot for point mutations during adaptation to antifungal drug stress within two distinct genetic backgrounds of Candida albicans. Heterozygous point mutations led to single allele dysfunction of ERG251 and resulted in azole tolerance in both genetic backgrounds. This is the first known example of point mutations causing azole tolerance in C. albicans. Importantly, single allele dysfunction of ERG251 in combination with recurrent chromosome aneuploidies resulted in bona fide azole resistance. Homozygous deletions of ERG251 caused increased fitness in low concentrations of fluconazole and decreased fitness in rich medium, especially at low initial cell density. Homozygous deletions of ERG251 resulted in accumulation of ergosterol intermediates consistent with the fitness defect in rich medium. Dysfunction of ERG251, together with FLC exposure, resulted in decreased accumulation of the toxic sterol (14-ɑ-methylergosta-8,24(28)-dien-3β,6α-diol) and increased accumulation of non-toxic alternative sterols. The altered sterol composition of the ERG251 mutants had pleiotropic effects on transcription, filamentation, and stress responses including cell membrane, osmotic and oxidative stress. Interestingly, while dysfunction of ERG251 resulted in azole tolerance, it also led to transcriptional upregulation of ZRT2, a membrane-bound Zinc transporter, in the presence of FLC, and overexpression of ZRT2 is sufficient to increase azole tolerance in wild-type C. albicans. Finally, in a murine model of systemic infection, homozygous deletion of ERG251 resulted in decreased virulence while the heterozygous deletion mutants maintain their pathogenicity. Overall, this study demonstrates that single allele dysfunction of ERG251 is a recurrent and effective mechanism of acquired azole tolerance. We propose that altered sterol composition resulting from ERG251 dysfunction mediates azole tolerance as well as pleiotropic effects on stress response, filamentation and virulence.

Chromosome-level genome assembly of the yeast Lodderomyces beijingensis reveals the genetic nature of metabolic adaptations and identifies subtelomeres as hotspots for amplification of mating type loci

Lodderomyces beijingensis is an ascosporic ascomycetous yeast. In contrast to related species Lodderomyces elongisporus, which is a recently emerging human pathogen, L. beijingensis is associated with insects. To provide an insight into its genetic makeup, we investigated the genome of its type strain, CBS 14171. We demonstrate that this yeast is diploid and describe the high contiguity nuclear genome assembly consisting of eight chromosome-sized contigs with a total size of about 15.1 Mbp. We find that the genome sequence contains multiple copies of the mating type loci and codes for essential components of the mating pheromone response pathway, however, the missing orthologs of several genes involved in the meiotic program raise questions about the mode of sexual reproduction. We also show that L. beijingensis genome codes for the 3-oxoadipate pathway enzymes, which allow the assimilation of protocatechuate. In contrast, the GAL gene cluster underwent a decay resulting in an inability of L. beijingensis to utilize galactose. Moreover, we find that the 56.5 kbp long mitochondrial DNA is structurally similar to known linear mitochondrial genomes terminating on both sides with covalently closed single-stranded hairpins. Finally, we discovered a new double-stranded RNA mycovirus from the Totiviridae family and characterized its genome sequence.

Spindle architecture constrains karyotype in budding yeast

The eukaryotic cell division machinery must rapidly and reproducibly duplicate and partition the cell's chromosomes in a carefully coordinated process. However, chromosome number varies dramatically between genomes, even on short evolutionary timescales. We sought to understand how the mitotic machinery senses and responds to karyotypic changes by using a series of budding yeast strains in which the native chromosomes have been successively fused. Using a combination of cell biological profiling, genetic engineering, and experimental evolution, we show that chromosome fusions are well tolerated up until a critical point. Cells with fewer than five centromeres lack the necessary number of kinetochore-microtubule attachments needed to counter outward forces in the metaphase spindle, triggering the spindle assembly checkpoint and prolonging metaphase. Our findings demonstrate that spindle architecture is a constraining factor for karyotype evolution.

Whole-genome duplication in the Multicellularity Long Term Evolution Experiment

Whole-genome duplication (WGD) is widespread across eukaryotes and can promote adaptive evolution1-4. However, given the instability of newly-formed polyploid genomes5-7, understanding how WGDs arise in a population, persist, and underpin adaptations remains a challenge. Using our ongoing Multicellularity Long Term Evolution Experiment (MuLTEE)8, we show that diploid snowflake yeast (Saccharomyces cerevisiae) under selection for larger multicellular size rapidly undergo spontaneous WGD. From its origin within the first 50 days of the experiment, tetraploids persist for the next 950 days (nearly 5,000 generations, the current leading edge of our experiment) in ten replicate populations, despite being genomically unstable. Using synthetic reconstruction, biophysical modeling, and counter-selection experiments, we found that tetraploidy evolved because it confers immediate fitness benefits in this environment, by producing larger, longer cells that yield larger clusters. The same selective benefit also maintained tetraploidy over long evolutionary timescales, inhibiting the reversion to diploidy that is typically seen in laboratory evolution experiments. Once established, tetraploidy facilitated novel genetic routes for adaptation, playing a key role in the evolution of macroscopic multicellular size via the origin of evolutionarily conserved aneuploidy. These results provide unique empirical insights into the evolutionary dynamics and impacts of WGD, showing how it can initially arise due to its immediate adaptive benefits, be maintained by selection, and fuel long-term innovations by creating additional dimensions of heritable genetic variation.

Unlocking the genome of the non-sourdough Kazachstania humilis MAW1: insights into inhibitory factors and phenotypic properties

BACKGROUND

Ascomycetous budding yeasts are ubiquitous environmental microorganisms important in food production and medicine. Due to recent intensive genomic research, the taxonomy of yeast is becoming more organized based on the identification of monophyletic taxa. This includes genera important to humans, such as Kazachstania. Until now, Kazachstania humilis (previously Candida humilis) was regarded as a sourdough-specific yeast. In addition, any antibacterial activity has not been associated with this species.

RESULTS

Previously, we isolated a yeast strain that impaired bio-hydrogen production in a dark fermentation bioreactor and inhibited the growth of Gram-positive and Gram-negative bacteria. Here, using next generation sequencing technologies, we sequenced the genome of this strain named K. humilis MAW1. This is the first genome of a K. humilis isolate not originating from a fermented food. We used novel phylogenetic approach employing the 18 S-ITS-D1-D2 region to show the placement of the K. humilis MAW1 among other members of the Kazachstania genus. This strain was examined by global phenotypic profiling, including carbon sources utilized and the influence of stress conditions on growth. Using the well-recognized bacterial model Escherichia coli AB1157, we show that K. humilis MAW1 cultivated in an acidic medium inhibits bacterial growth by the disturbance of cell division, manifested by filament formation. To gain a greater understanding of the inhibitory effect of K. humilis MAW1, we selected 23 yeast proteins with recognized toxic activity against bacteria and used them for Blast searches of the K. humilis MAW1 genome assembly. The resulting panel of genes present in the K. humilis MAW1 genome included those encoding the 1,3-β-glucan glycosidase and the 1,3-β-glucan synthesis inhibitor that might disturb the bacterial cell envelope structures.

CONCLUSIONS

We characterized a non-sourdough-derived strain of K. humilis, including its genome sequence and physiological aspects. The MAW1, together with other K. humilis strains, shows the new organization of the mating-type locus. The revealed here pH-dependent ability to inhibit bacterial growth has not been previously recognized in this species. Our study contributes to the building of genome sequence-based classification systems; better understanding of K.humilis as a cell factory in fermentation processes and exploring bacteria-yeast interactions in microbial communities.

Host-Pathogen Interactions and Correlated Factors That Are Affected in Replicative-Aged Cryptococcus neoformans

Cryptococcus neoformans is a facultative intracellular fungal pathogen. Ten-generation-old (10GEN) C. neoformans cells are more resistant to phagocytosis and killing by macrophages than younger daughter cells. However, mechanisms that mediate this resistance and intracellular parasitism are poorly understood. Here, we identified important factors for the intracellular survival of 10GEN C. neoformans, such as urease activity, capsule synthesis, and DNA content using flow cytometry and fluorescent microscopy techniques. The real-time visualization of time-lapse imaging was applied to determine the phagosomal acidity, membrane permeability, and vomocytosis (non-lytic exocytosis) rate in J774 macrophages that phagocytosed C. neoformans of different generational ages. Our results showed that old C. neoformans exhibited higher urease activity and enhanced Golgi activity. In addition, old C. neoformans were more likely to be arrested in the G2 phase, resulting in the occasional formation of aberrant trimera-like cells. To finish, the advanced generational age of the yeast cells slightly reduced vomocytosis events within host cells, which might be associated with increased phagolysosome pH and membrane permeability. Altogether, our results suggest that old C. neoformans prevail within acidic phagolysosomes and can manipulate the phagosome pH. These strategies may be used by old C. neoformans to resist phagosomal killing and drive cryptococcosis pathogenesis. The comprehension of these essential host-pathogen interactions could further shed light on mechanisms that bring new insights for novel antifungal therapeutic design.

On the origin of bird's nest fungi: Phylogenomic analyses of fungi in the Nidulariaceae (Agaricales, Basidiomycota)

Nidulariaceae, also known as bird's nest fungi, is an understudied group of mushroom-forming fungi. The common name is derived from their nest-like morphology. Bird's nest fungi are ubiquitous wood decomposers or saprobes on dung. Recent studies showed that species in the Nidulariaceae form a monophyletic group with five sub-clades. However, phylogenetic relationships among genera and placement of Nidulariaceae are still unclear. We present phylogenomic analyses of bird's nest fungi and related Agaricales fungi to gain insight into the evolution of Nidulariaceae. A species tree with 17 newly generated genomes of bird's nest fungi and representatives from all major clades of Agaricales was constructed using 1044 single-copy genes to explore the intergeneric relationships and pinpoint the placement of Nidulariaceae within Agaricales. We corroborated the hypothesis that bird's nest fungi are sister to Squamanitaceae, which includes mushroom-shaped fungi with a stipe and pileus that are saprobes and mycoparasites. Lastly, stochastic character mapping of discrete traits on phylogenies (SIMMAP) suggests that the ancestor of bird's nest fungi likely possessed an evanescent, globose peridium without strings attaching to the spore packets (funiculi). This analysis suggests that the funiculus was gained twice and that the persistent, cupulate peridium form was gained at least four times and lost once. However, alternative coding schemes and datasets with a wider array of Agaricales produced conflicting results during ancestral state reconstruction, indicating that there is some uncertainty in the number of peridium transitions and that taxon sampling may significantly alter ancestral state reconstructions. Overall, our results suggest that several key morphological characters of Nidulariaceae have been subject to homoplasy.

Chromosome-level Subgenome-aware de novo Assembly of Saccharomyces bayanus Provides Insight into Genome Divergence after Hybridization

Interspecies hybridization is prevalent in various eukaryotic lineages and plays important roles in phenotypic diversification, adaption, and speciation. To better understand the changes that occurred in the different subgenomes of a hybrid species and how they facilitated adaptation, we completed chromosome-level de novo assemblies of all 16 pairs chromosomes for a recently formed hybrid yeast, Saccharomyces bayanus strain CBS380 (IFO11022), using Nanopore MinION long-read sequencing. Characterization of S. bayanus subgenomes and comparative analysis with the genomes of its parent species, S. uvarum and S. eubayanus, provide several new insights into understanding genome evolution after a relatively recent hybridization. For instance, multiple recombination events between the two subgenomes have been observed in each chromosome, followed by loss of heterozygosity (LOH) in most chromosomes in nine chromosome pairs. In addition to maintaining nearly all gene content and synteny from its parental genomes, S. bayanus has acquired many genes from other yeast species, primarily through the introgression of S. cerevisiae, such as those involved in the maltose metabolism. In addition, the patterns of recombination and LOH suggest an allotetraploid origin of S. bayanus. The gene acquisition and rapid LOH in the hybrid genome probably facilitated its adaption to maltose brewing environments and mitigated the maladaptive effect of hybridization.

Erg251 has complex and pleiotropic effects on azole susceptibility, filamentation, and stress response phenotypes

Ergosterol is essential for fungal cell membrane integrity and growth, and numerous antifungal drugs target ergosterol. Inactivation or modification of ergosterol biosynthetic genes can lead to changes in antifungal drug susceptibility, filamentation and stress response. Here, we found that the ergosterol biosynthesis gene ERG251 is a hotspot for point mutations during adaptation to antifungal drug stress within two distinct genetic backgrounds of Candida albicans. Heterozygous point mutations led to single allele dysfunction of ERG251 and resulted in azole tolerance in both genetic backgrounds. This is the first known example of point mutations causing azole tolerance in C. albicans. Importantly, single allele dysfunction of ERG251 in combination with recurrent chromosome aneuploidies resulted in bona fide azole resistance. Homozygous deletions of ERG251 caused increased fitness in low concentrations of fluconazole and decreased fitness in rich medium, especially at low initial cell density. Dysfunction of ERG251 resulted in transcriptional upregulation of the alternate sterol biosynthesis pathway and ZRT2, a Zinc transporter. Notably, we determined that overexpression of ZRT2 is sufficient to increase azole tolerance in C. albicans. Our combined transcriptional and phenotypic analyses revealed the pleiotropic effects of ERG251 on stress responses including cell wall, osmotic and oxidative stress. Interestingly, while loss of either allele of ERG251 resulted in similar antifungal drug responses, we observed functional divergence in filamentation regulation between the two alleles of ERG251 (ERG251-A and ERG251-B) with ERG251-A exhibiting a dominant role in the SC5314 genetic background. Finally, in a murine model of systemic infection, homozygous deletion of ERG251 resulted in decreased virulence while the heterozygous deletion mutants maintain their pathogenicity. Overall, this study provides extensive genetic, transcriptional and phenotypic analysis for the effects of ERG251 on drug susceptibility, fitness, filamentation and stress responses.

Determining mating type and ploidy in Rhodotorula toruloides and its effect on growth on sugars from lignocellulosic biomass

Rhodotorula toruloides is being developed for the use in industrial biotechnology processes because of its favorable physiology. This includes its ability to produce and store large amounts of lipids in the form of intracellular lipid bodies. Nineteen strains were characterized for mating type, ploidy, robustness for growth, and accumulation of lipids on inhibitory switchgrass hydrolysate (SGH). Mating type was determined using a novel polymerase chain reaction (PCR)-based assay, which was validated using the classical microscopic test. Three of the strains were heterozygous for mating type (A1/A2). Ploidy analysis revealed a complex pattern. Two strains were triploid, eight haploid, and eight either diploid or aneuploid. Two of the A1/A2 strains were compared to their parents for growth on 75%v/v concentrated SGH. The A1/A2 strains were much more robust than the parental strains, which either did not grow or had extended lag times. The entire set was evaluated in 60%v/v SGH batch cultures for growth kinetics and biomass and lipid production. Lipid titers were 2.33-9.40 g/L with a median of 6.12 g/L, excluding the two strains that did not grow. Lipid yields were 0.032-0.131 (g/g) and lipid contents were 13.5-53.7% (g/g). Four strains had significantly higher lipid yields and contents. One of these strains, which had among the highest lipid yield in this study (0.131 ± 0.007 g/g), has not been previously described in the literature.

SUMMARY

The yeast Rhodotorula toruloides was used to produce oil using sugars extracted from a bioenergy grass.

Antifungal Drug Concentration Impacts the Spectrum of Adaptive Mutations in Candida albicans

Invasive fungal infections are a leading global cause of human mortality. Only three major classes of antifungal drugs are widely used, and resistance to all three classes can arise rapidly. The most widely prescribed antifungal drug, fluconazole, disseminates rapidly and reaches a wide range of concentrations throughout the body. The impact of drug concentration on the spectrum and effect of mutations acquired during adaptation is not known for any fungal pathogen, and how the specific level of a given stress influences the distribution of beneficial mutations has been poorly explored in general. We evolved 144 lineages from three genetically distinct clinical isolates of Candida albicans to four concentrations of fluconazole (0, 1, 8, and 64 μg/ml) and performed comprehensive phenotypic and genomic comparisons of ancestral and evolved populations. Adaptation to different fluconazole concentrations resulted in distinct adaptive trajectories. In general, lineages evolved to drug concentrations close to their MIC50 (the level of drug that reduces growth by 50% in the ancestor) tended to rapidly evolve an increased MIC50 and acquired distinct segmental aneuploidies and copy number variations. By contrast, lineages evolved to drug concentrations above their ancestral MIC50 tended to acquire a different suite of mutational changes and increased in drug tolerance (the ability of a subpopulation of cells to grow slowly above their MIC50). This is the first evidence that different concentrations of drug can select for different genotypic and phenotypic outcomes in vitro and may explain observed in vivo drug response variation.

Copy Number Variation and Allele Ratio Analysis in Candida albicans Using Whole Genome Sequencing Data

Whole genome sequencing of human fungal pathogens has revolutionized the speed and accuracy in which sequence variants that cause antifungal resistance can be identified. Genome rearrangements resulting in copy number variation (CNV) are a significant source of acquired antifungal drug resistance across diverse fungal species. Some CNVs are transient in nature, while other CNVs are stable and well tolerated even in the absence of antifungal drugs. By visualizing whole genome sequencing read depth as a function of genomic location, CNVs and CNV breakpoints (genomic positions where the copy number changes occur relative to the rest of the genome) are rapidly identified. A similar analysis can be used to visualize allele ratio changes that occur across the genomes of heterozygous fungal species, both in the presence and absence of CNVs. This protocol walks through the bioinformatic analysis of CNVs and allele ratios utilizing free, open-source visualization tools. We provide code to use with an example dataset (matched antifungal drug-sensitive and drug-resistant Candida albicans isolates) and notes on how to expand this protocol to other fungal genomes.

Reduction of Ethyl Carbamate in an Alcoholic Beverage by CRISPR/Cas9-Based Genome Editing of the Wild Yeast

Ethyl carbamate (EC) is a naturally occurring substance in alcoholic beverages from the reaction of ethanol with urea during fermentation and storage. EC can cause dizziness and vomiting when consumed in small quantities and develop kidney cancer when consumed in excess. Thus, the reduction of EC formation in alcoholic beverages is important for food safety and human health. One of the strategies for reducing EC contents in alcoholic beverages is developing a new yeast starter strain to enable less formation of EC during fermentation. In this study, we isolated a polyploid wild-type yeast Saccharomyces cerevisiae strain from the Nuruk (Korean traditional grain-based inoculum of wild yeast and mold) and developed a starter culture by genome engineering to reduce EC contents in alcoholic beverages. We deleted multiple copies of the target genes involved in the EC formation in the S. cerevisiae by a CRISPR/Cas9-based genome editing tool. First, the CAR1 gene encoding for the arginase enzyme responsible for the formation of urea was completely deleted in the genome of S. cerevisiae. Additionally, the GZF3 gene encoding the transcription factor controlling expression levels of several genes (DUR1, 2, and DUR3) related to urea absorption and degradation was deleted in S. cerevisiae to further reduce the EC formation. The effects of gene deletion were validated by RT-qPCR to confirm changes in transcriptional levels of the EC-related genes. The resulting strain of S. cerevisiae carrying a double deletion of CAR1 and GZF3 genes successfully reduced the EC contents in the fermentation medium without significant changes in alcohol contents and fermentation profiles when compared to the wild-type strain. Finally, we brewed the Korean traditional rice wine Makgeolli using the double deletion strain of S. cerevisiae dCAR1&GZF3, resulting in a significant reduction of the EC content in Makgeolli up to 41.6% when compared to the wild-type strain. This study successfully demonstrated the development of a starter culture to reduce the EC formation in an alcoholic beverage by CRISPR/Cas9 genome editing of the wild yeast.

A dissemination-prone morphotype enhances extrapulmonary organ entry by Cryptococcus neoformans

Environmental pathogens move from ecological niches to mammalian hosts, requiring adaptation to dramatically different environments. Microbes that disseminate farther, including the fungal meningitis pathogen Cryptococcus neoformans, require additional adaptation to diverse tissues. We demonstrate that the formation of a small C. neoformans morphotype-called "seed" cells due to their colonizing ability-is critical for extrapulmonary organ entry. Seed cells exhibit changes in fungal cell size and surface expression that result in an enhanced macrophage update. Seed cell formation is triggered by environmental factors, including C. neoformans' environmental niche, and pigeon guano with phosphate plays a central role. Seed cells show the enhanced expression of phosphate acquisition genes, and mutants unable to acquire phosphate fail to adopt the seed cell morphotype. Additionally, phosphate can be released by tissue damage, potentially establishing a feed-forward loop of seed cell formation and dissemination. Thus, C. neoformans' size variation represent inducible morphotypes that change host interactions to facilitate microbe spread.

Candida auris on Apples: Diversity and Clinical Significance

Candida auris is a multidrug-resistant nosocomial fungal pathogen. While the marine environment was recently identified as a natural niche for C. auris, the environment(s) that might have contributed to the development and spread of antifungal resistance in C. auris remains a mystery. Because stored fruits are often treated with fungicides to prevent postharvest spoilage, we hypothesized that stored fruits could serve as a possible selective force for and a transmission reservoir of antifungal-resistant isolates of pathogenic yeasts, including C. auris. To test this hypothesis, we screened fruits to study the diversity of pathogenic yeasts and their antifungal susceptibility profiles. Among the 62 screened apples, the surfaces of 8 were positive for C. auris, and all were stored apples. Whole-genome sequencing (WGS) showed that C. auris strains from apples were genetically diverse and exhibited broad phylogenetic distribution among the subclades within clade I. Interestingly, strains from apples had closely related strains from other sources in India, including from patients, hospitals, and marine environments, and from clinical strains from other parts of the world. A broad range of fungicides, including dimethyl inhibitors (DMIs), were detected in stored apples, and all C. auris isolates exhibited reduced sensitivity to DMIs. Interestingly, C. auris was not isolated from freshly picked apples. Together, the results suggest a potentially complex ecology for C. auris with agriculture fungicide application on stored fruits as a significant selective force for drug resistance in clinics.

Aneuploidy Underlies Tolerance and Cross-Tolerance to Drugs in Candida parapsilosis

Candida species are the most common human fungal pathogens worldwide. Although C. albicans remains the predominant cause of candidiasis, infections caused by non-albicans Candida species, including C. parapsilosis, are increasing. In C. albicans, genome plasticity has been shown to be a prevalent strategy of adaptation to stresses. However, the role of aneuploidy in C. parapsilosis is largely unknown. In this study, we found that six different aneuploid karyotypes conferred adaptation to the endoplasmic reticulum stress inducer tunicamycin (TUN) in C. parapsilosis. Interestingly, a specific aneuploidy including trisomy of chromosome 6 (Chr6x3) also enabled cross-tolerance to aureobasidin A (AbA), a sphingolipid biosynthesis inhibitor. Consistent with this, selection on AbA identified adaptors with three different aneuploid karyotypes, including Chr6x3, which also enabled cross-tolerance to both AbA and TUN. Therefore, as in other Candida species, recurrent aneuploid karyotypes enable the adaptation of C. parapsilosis to specific stresses, and specific aneuploidies enable cross-adaptation to different stresses. IMPORTANCE Candida parapsilosis is an emerging human fungal pathogen, especially prevalent in neonates. Aneuploidy, having uneven numbers of chromosomes, is a well-known mechanism for adapting to stress in Candida albicans, the most common human fungal pathogen. In this study, we exposed C. parapsilosis to two very different drugs and selected for rare cells that grew in one of the drugs. We found that the majority of isolates that grew in the drugs had acquired an extra copy of one of several aneuploid chromosomes and that specific aneuploid chromosomes appeared in several independent cell clones. Importantly, an extra copy of chromosome 6 was detected following selection in either one of the drugs, and this extra chromosome conferred the ability to grow in both drugs, a property called cross-adaptation, or cross-tolerance. Thus, this study highlights the genome plasticity of C. parapsilosis and the ability of an extra copy of a single chromosome to promote cell growth in the presence of more than one drug.

Efficient and crucial quality control of HAP1 cell ploidy status

The near-haploid human cell line HAP1 recently became a popular subject for CRISPR/Cas9 editing, since only one allele requires modification. Through the gene-editing service at Horizon Discovery, there are at present more than 7500 edited cell lines available and the number continuously increases. The haploid nature of HAP1 is unstable as cultures become diploid with time. Here, we demonstrated some fundamental differences between haploid and diploid HAP1 cells, hence underlining the need for taking control over ploidy status in HAP1 cultures prior to phenotyping. Consequently, we optimized a procedure to determine the ploidy of HAP1 by flow cytometry in order to obtain diploid cultures and avoid ploidy status as an interfering variable in experiments. Furthermore, in order to facilitate this quality control, we validated a size-based cell sorting procedure to obtain the diploid culture more rapidly. Hence, we provide here two streamlined protocols for quality controlling the ploidy of HAP1 cells and document their validity and necessity.

This article has an associated First Person interview with the co-first authors of the paper.

Summary: Sharing an effective procedure to quality control the near-haploid HAP1 cells for standardized comparison to CRISPR/Cas9 modified versions and demonstrating the need for controlling the spontaneous diploidization of HAP1 cultures.

A Target Capture-Based Method to Estimate Ploidy From Herbarium Specimens

Whole genome duplication (WGD) events are common in many plant lineages, but the ploidy status and possible occurrence of intraspecific ploidy variation are unknown for most species. Standard methods for ploidy determination are chromosome counting and flow cytometry approaches. While flow cytometry approaches typically use fresh tissue, an increasing number of studies have shown that recently dried specimens can be used to yield ploidy data. Recent studies have started to explore whether high-throughput sequencing (HTS) data can be used to assess ploidy levels by analyzing allelic frequencies from single copy nuclear genes. Here, we compare different approaches using a range of yam (Dioscorea) tissues of varying ages, drying methods and quality, including herbarium tissue. Our aims were to: (1) explore the limits of flow cytometry in estimating ploidy level from dried samples, including herbarium vouchers collected between 1831 and 2011, and (2) optimize a HTS-based method to estimate ploidy by considering allelic frequencies from nuclear genes obtained using a target-capture method. We show that, although flow cytometry can be used to estimate ploidy levels from herbarium specimens collected up to fifteen years ago, success rate is low (5.9%). We validated our HTS-based estimates of ploidy using 260 genes by benchmarking with dried samples of species of known ploidy (Dioscorea alata, D. communis, and D. sylvatica). Subsequently, we successfully applied the method to the 85 herbarium samples analyzed with flow cytometry, and successfully provided results for 91.7% of them, comprising species across the phylogenetic tree of Dioscorea. We also explored the limits of using this HTS-based approach for identifying high ploidy levels in herbarium material and the effects of heterozygosity and sequence coverage. Overall, we demonstrated that ploidy diversity within and between species may be ascertained from historical collections, allowing the determination of polyploidization events from samples collected up to two centuries ago. This approach has the potential to provide insights into the drivers and dynamics of ploidy level changes during plant evolution and crop domestication.

Genome plasticity in Candida albicans is driven by long repeat sequences

Genome rearrangements resulting in copy number variation (CNV) and loss of heterozygosity (LOH) are frequently observed during the somatic evolution of cancer and promote rapid adaptation of fungi to novel environments. In the human fungal pathogen Candida albicans, CNV and LOH confer increased virulence and antifungal drug resistance, yet the mechanisms driving these rearrangements are not completely understood. Here, we unveil an extensive array of long repeat sequences (65-6499 bp) that are associated with CNV, LOH, and chromosomal inversions. Many of these long repeat sequences are uncharacterized and encompass one or more coding sequences that are actively transcribed. Repeats associated with genome rearrangements are predominantly inverted and separated by up to ~1.6 Mb, an extraordinary distance for homology-based DNA repair/recombination in yeast. These repeat sequences are a significant source of genome plasticity across diverse strain backgrounds including clinical, environmental, and experimentally evolved isolates, and represent previously uncharacterized variation in the reference genome.