YMAP: a pipeline for visualization of copy number variation and loss of heterozygosity in eukaryotic pathogens

aCNViewer: Comprehensive genome-wide visualization of absolute copy number and copy neutral variations

Motivation

Copy number variations (CNV) include net gains or losses of part or whole chromosomal regions. They differ from copy neutral loss of heterozygosity (cn-LOH) events which do not induce any net change in the copy number and are often associated with uniparental disomy. These phenomena have long been reported to be associated with diseases and particularly in cancer. Losses/gains of genomic regions are often correlated with lower/higher gene expression. On the other hand, loss of heterozygosity (LOH) and cn-LOH are common events in cancer and may be associated with the loss of a functional tumor suppressor gene. Therefore, identifying recurrent CNV and cn-LOH events can be important as they may highlight common biological components and give insights into the development or mechanisms of a disease. However, no currently available tools allow a comprehensive whole-genome visualization of recurrent CNVs and cn-LOH in groups of samples providing absolute quantification of the aberrations leading to the loss of potentially important information.

Results

To overcome these limitations, we developed aCNViewer (Absolute CNV Viewer), a visualization tool for absolute CNVs and cn-LOH across a group of samples. aCNViewer proposes three graphical representations: dendrograms, bi-dimensional heatmaps showing chromosomal regions sharing similar abnormality patterns, and quantitative stacked histograms facilitating the identification of recurrent absolute CNVs and cn-LOH. We illustrated aCNViewer using publically available hepatocellular carcinomas (HCCs) Affymetrix SNP Array data (Fig 1A). Regions 1q and 8q present a similar percentage of total gains but significantly different copy number gain categories (p-value of 0.0103 with a Fisher exact test), validated by another cohort of HCCs (p-value of 5.6e-7) (Fig 2B).

Availability and implementation

aCNViewer is implemented in python and R and is available with a GNU GPLv3 license on GitHub https://github.com/FJD-CEPH/aCNViewer and Docker https://hub.docker.com/r/fjdceph/acnviewer/.

Contact

aCNViewer@cephb.fr

Ploidy dynamics and evolvability in fungi

Rapid responses to acute stresses are essential for stress survival and are critical to the ability of fungal pathogens to adapt to new environments or hosts. The rapid emergence of drug resistance is used as a model for how fungi adapt and survive stress conditions that inhibit the growth of progenitor cells. Aneuploidy and loss of heterozygosity (LOH), which are large-scale genome shifts involving whole chromosomes or chromosome arms, occur at higher frequency than point mutations and have the potential to mediate stress survival. Furthermore, the stress of exposure to an antifungal drug can induce elevated levels of LOH and can promote the formation of aneuploids. This occurs via mitotic defects that first produce tetraploid progeny with extra spindles, followed by chromosome mis-segregation. Thus, drug exposure induces elevated levels of aneuploidy, which can alter the copy number of genes that improve survival in a given stress or drug. Selection then acts to increase the proportion of adaptive aneuploids in the population. Because aneuploidy is a common property of many pathogenic fungi, including those posing emerging threats to plants, animals and humans, we propose that aneuploid formation and LOH often accompanying it contribute to the rapid generation of diversity that can facilitate the emergence of fungal pathogens to new environmental niches and/or new hosts, as well as promote antifungal drug resistance that makes emerging fungal infections ever more difficult to contain.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Loss of C-5 Sterol Desaturase Activity Results in Increased Resistance to Azole and Echinocandin Antifungals in a Clinical Isolate of Candida parapsilosis

Among emerging non-albicans Candida species, Candida parapsilosis is of particular concern as a cause of nosocomial bloodstream infections in neonatal and intensive care unit patients. While fluconazole and echinocandins are considered effective treatments for such infections, recent reports of fluconazole and echinocandin resistance in C. parapsilosis indicate a growing problem. The present study describes a novel mechanism of antifungal resistance in this organism affecting susceptibility to azole and echinocandin antifungals in a clinical isolate obtained from a patient with prosthetic valve endocarditis. Transcriptome analysis indicated differential expression of several genes in the resistant isolate, including upregulation of ergosterol biosynthesis pathway genes ERG2, ERG5, ERG6, ERG11, ERG24, ERG25, and UPC2 Whole-genome sequencing revealed that the resistant isolate possessed an ERG3 mutation resulting in a G111R amino acid substitution. Sterol profiles indicated a reduction in sterol desaturase activity as a result of this mutation. Replacement of both mutant alleles in the resistant isolate with the susceptible isolate's allele restored wild-type susceptibility to all azoles and echinocandins tested. Disruption of ERG3 in the susceptible and resistant isolates resulted in a loss of sterol desaturase activity, high-level azole resistance, and an echinocandin-intermediate to -resistant phenotype. While disruption of ERG3 in C. albicans resulted in azole resistance, echinocandin MICs, while elevated, remained within the susceptible range. This work demonstrates that the G111R substitution in Erg3 is wholly responsible for the altered azole and echinocandin susceptibilities observed in this C. parapsilosis isolate and is the first report of an ERG3 mutation influencing susceptibility to the echinocandins.

Adaptive Mistranslation Accelerates the Evolution of Fluconazole Resistance and Induces Major Genomic and Gene Expression Alterations in Candida albicans

Regulated erroneous protein translation (adaptive mistranslation) increases proteome diversity and produces advantageous phenotypic variability in the human pathogen Candida albicans. It also increases fitness in the presence of fluconazole, but the underlying molecular mechanism is not understood. To address this question, we evolved hypermistranslating and wild-type strains in the absence and presence of fluconazole and compared their fluconazole tolerance and resistance trajectories during evolution. The data show that mistranslation increases tolerance and accelerates the acquisition of resistance to fluconazole. Genome sequencing, array-based comparative genome analysis, and gene expression profiling revealed that during the course of evolution in fluconazole, the range of mutational and gene deregulation differences was distinctively different and broader in the hypermistranslating strain, including multiple chromosome duplications, partial chromosome deletions, and polyploidy. Especially, the increased accumulation of loss-of-heterozygosity events, aneuploidy, translational and cell surface modifications, and differences in drug efflux seem to mediate more rapid drug resistance acquisition under mistranslation. Our observations support a pivotal role for adaptive mistranslation in the evolution of drug resistance in C. albicans. IMPORTANCE Infectious diseases caused by drug-resistant fungi are an increasing threat to public health because of the high mortality rates and high costs associated with treatment. Thus, understanding of the molecular mechanisms of drug resistance is of crucial interest for the medical community. Here we investigated the role of regulated protein mistranslation, a characteristic mechanism used by C. albicans to diversify its proteome, in the evolution of fluconazole resistance. Such codon ambiguity is usually considered highly deleterious, yet recent studies found that mistranslation can boost adaptation in stressful environments. Our data reveal that CUG ambiguity diversifies the genome in multiple ways and that the full spectrum of drug resistance mechanisms in C. albicans goes beyond the traditional pathways that either regulate drug efflux or alter the interactions of drugs with their targets. The present work opens new avenues to understand the molecular and genetic basis of microbial drug resistance.

Tolerance to Caspofungin in Candida albicans Is Associated with at Least Three Distinctive Mechanisms That Govern Expression of FKS Genes and Cell Wall Remodeling

Expanding echinocandin use to prevent or treat invasive fungal infections has led to an increase in the number of breakthrough infections due to resistant Candida species. Although it is uncommon, echinocandin resistance is well documented for Candida albicans, which is among the most prevalent bloodstream organisms. A better understanding is needed to assess the cellular factors that promote tolerance and predispose infecting cells to clinical breakthrough. We previously showed that some mutants that were adapted to growth in the presence of toxic sorbose due to loss of one chromosome 5 (Ch5) also became more tolerant to caspofungin. We found here, following direct selection of mutants on caspofungin, that tolerance can be conferred by at least three mechanisms: (i) monosomy of Ch5, (ii) combined monosomy of the left arm and trisomy of the right arm of Ch5, and (iii) an aneuploidy-independent mechanism. Tolerant mutants possessed cell walls with elevated chitin and showed downregulation of genes involved in cell wall biosynthesis, namely, FKS, located outside Ch5, and CHT2, located on Ch5, irrespective of Ch5 ploidy. Also irrespective of Ch5 ploidy, the CNB1 and MID1 genes on Ch5, which are involved in the calcineurin signaling pathway, were expressed at the diploid level. Thus, multiple mechanisms can affect the relative expression of the aforementioned genes, controlling them in similar ways. Although breakthrough mutations in two specific regions of FKS1 have previously been associated with caspofungin resistance, we found mechanisms of caspofungin tolerance that are independent of FKS1 and thus represent an earlier event in resistance development.

Ploidy tug-of-war: Evolutionary and genetic environments influence the rate of ploidy drive in a human fungal pathogen

Variation in baseline ploidy is seen throughout the tree of life, yet the factors that determine why one ploidy level is maintained over another remain poorly understood. Experimental evolution studies using asexual fungal microbes with manipulated ploidy levels intriguingly reveals a propensity to return to the historical baseline ploidy, a phenomenon that we term "ploidy drive." We evolved haploid, diploid, and polyploid strains of the human fungal pathogen Candida albicans under three different nutrient limitation environments to test whether these conditions, hypothesized to select for low ploidy levels, could counteract ploidy drive. Strains generally maintained or acquired smaller genome sizes (measured as total nuclear DNA through flow cytometry) in minimal medium and under phosphorus depletion compared to in a complete medium, while mostly maintained or acquired increased genome sizes under nitrogen depletion. Improvements in fitness often ran counter to changes in genome size; in a number of scenarios lines that maintained their original genome size often increased in fitness more than lines that converged toward diploidy (the baseline ploidy of C. albicans). Combined, this work demonstrates a role for both the environment and genotype in determination of the rate of ploidy drive, and highlights questions that remain about the force(s) that cause genome size variation.

Haplotyping a Non-meiotic Diploid Fungal Pathogen Using Induced Aneuploidies and SNP/CGH Microarray Analysis

The generation of haplotype information has recently become very attractive due to its utility for identifying mutations associated with human disease and for the development of personalized medicine. Haplotype information also is crucial for studying recombination mechanisms and genetic diversity, and for analyzing allele-specific gene expression. Classic haplotyping methods require the analysis of hundreds of meiotic progeny. To facilitate haplotyping in the non-meiotic human fungal pathogen Candida albicans, we exploited trisomic heterozygous chromosomes generated via the UAU1 selection strategy. Using this system, we obtained phasing information from allelic biases, detected by SNP/CGH microarray analysis. This strategy has the potential to be applicable to other diploid, asexual Candida species that are important causes of human disease.

Analysis of Repair Mechanisms following an Induced Double-Strand Break Uncovers Recessive Deleterious Alleles in the Candida albicans Diploid Genome

The diploid genome of the yeast Candida albicans is highly plastic, exhibiting frequent loss-of-heterozygosity (LOH) events. To provide a deeper understanding of the mechanisms leading to LOH, we investigated the repair of a unique DNA double-strand break (DSB) in the laboratory C. albicans SC5314 strain using the I-SceI meganuclease. Upon I-SceI induction, we detected a strong increase in the frequency of LOH events at an I-SceI target locus positioned on chromosome 4 (Chr4), including events spreading from this locus to the proximal telomere. Characterization of the repair events by single nucleotide polymorphism (SNP) typing and whole-genome sequencing revealed a predominance of gene conversions, but we also observed mitotic crossover or break-induced replication events, as well as combinations of independent events. Importantly, progeny that had undergone homozygosis of part or all of Chr4 haplotype B (Chr4B) were inviable. Mining of genome sequencing data for 155 C. albicans isolates allowed the identification of a recessive lethal allele in the GPI16 gene on Chr4B unique to C. albicans strain SC5314 which is responsible for this inviability. Additional recessive lethal or deleterious alleles were identified in the genomes of strain SC5314 and two clinical isolates. Our results demonstrate that recessive lethal alleles in the genomes of C. albicans isolates prevent the occurrence of specific extended LOH events. While these and other recessive lethal and deleterious alleles are likely to accumulate in C. albicans due to clonal reproduction, their occurrence may in turn promote the maintenance of corresponding nondeleterious alleles and, consequently, heterozygosity in the C. albicans species.

IMPORTANCE

Recessive lethal alleles impose significant constraints on the biology of diploid organisms. Using a combination of an I-SceI meganuclease-mediated DNA DSB, a fluorescence-activated cell sorter (FACS)-optimized reporter of LOH, and a compendium of 155 genome sequences, we were able to unmask and identify recessive lethal and deleterious alleles in isolates of Candida albicans, a diploid yeast and the major fungal pathogen of humans. Accumulation of recessive deleterious mutations upon clonal reproduction of C. albicans could contribute to the maintenance of heterozygosity despite the high frequency of LOH events in this species.

Polyploid titan cells produce haploid and aneuploid progeny to promote stress adaptation

Cryptococcus neoformans is a major life-threatening fungal pathogen. In response to the stress of the host environment, C. neoformans produces large polyploid titan cells. Titan cell production enhances the virulence of C. neoformans, yet whether the polyploid aspect of titan cells is specifically influential remains unknown. We show that titan cells were more likely to survive and produce offspring under multiple stress conditions than typical cells and that even their normally sized daughters maintained an advantage over typical cells in continued exposure to stress. Although polyploid titan cells generated haploid daughter cell progeny upon in vitro replication under nutrient-replete conditions, titan cells treated with the antifungal drug fluconazole produced fluconazole-resistant diploid and aneuploid daughter cells. Interestingly, a single titan mother cell was capable of generating multiple types of aneuploid daughter cells. The increased survival and genomic diversity of titan cell progeny promote rapid adaptation to new or high-stress conditions.

IMPORTANCE

The ability to adapt to stress is a key element for survival of pathogenic microbes in the host and thus plays an important role in pathogenesis. Here we investigated the predominantly haploid human fungal pathogen Cryptococcus neoformans, which is capable of ploidy and cell size increases during infection through production of titan cells. The enlarged polyploid titan cells are then able to rapidly undergo ploidy reduction to generate progeny with reduced ploidy and/or aneuploidy. Under stressful conditions, titan cell progeny have a growth and survival advantage over typical cell progeny. Understanding how titan cells enhance the rate of cryptococcal adaptation under stress conditions may assist in the development of novel drugs aimed at blocking ploidy transitions.

Budding off: bringing functional genomics to Candida albicans

Candida species are the most prevalent human fungal pathogens, with Candida albicans being the most clinically relevant species. Candida albicans resides as a commensal of the human gastrointestinal tract but is a frequent cause of opportunistic mucosal and systemic infections. Investigation of C. albicans virulence has traditionally relied on candidate gene approaches, but recent advances in functional genomics have now facilitated global, unbiased studies of gene function. Such studies include comparative genomics (both between and within Candida species), analysis of total RNA expression, and regulation and delineation of protein-DNA interactions. Additionally, large collections of mutant strains have begun to aid systematic screening of clinically relevant phenotypes. Here, we will highlight the development of functional genomics in C. albicans and discuss the use of these approaches to addressing both commensalism and pathogenesis in this species.

Shift and adapt: the costs and benefits of karyotype variations

Variation is the spice of life or, in the case of evolution, variation is the necessary material on which selection can act to enable adaptation. Karyotypic variation in ploidy (the number of homologous chromosome sets) and aneuploidy (imbalance in the number of chromosomes) are fundamentally different than other types of genomic variants. Karyotypic variation emerges through different molecular mechanisms than other mutational events, and unlike mutations that alter the genome at the base pair level, rapid reversion to the wild type chromosome number is often possible. Although karyotypic variation has long been noted and discussed by biologists, interest in the importance of karyotypic variants in evolutionary processes has spiked in recent years, and much remains to be discovered about how karyotypic variants are produced and subsequently selected.

Parasexual Ploidy Reduction Drives Population Heterogeneity Through Random and Transient Aneuploidy in Candida albicans

The opportunistic pathogen Candida albicans has a large repertoire of mechanisms to generate genetic and phenotypic diversity despite the lack of meiosis in its life cycle. Its parasexual cycle enables shifts in ploidy, which in turn facilitate recombination, aneuploidy, and homozygosis of whole chromosomes to fuel rapid adaptation. Here we show that the tetraploid state potentiates ploidy variation and drives population heterogeneity. In tetraploids, the rate of losing a single heterozygous marker [loss of heterozygosity (LOH)] is elevated ∼30-fold higher than the rate in diploid cells. Furthermore, isolates recovered after selection for LOH of one, two, or three markers were highly aneuploid, with a broad range of karyotypes including strains with a combination of di-, tri-, and tetrasomic chromosomes. We followed the ploidy trajectories for these tetraploid- and aneuploid-derived isolates, using a combination of flow cytometry and double-digestion restriction-site-associated DNA analyzed with next-generation sequencing. Isolates derived from either tetraploid or aneuploid isolates predominately resolved to a stable euploid state. The majority of isolates reduced to the conventional diploid state; however, stable triploid and tetraploid states were observed in ∼30% of the isolates. Notably, aneuploid isolates were more transient than tetraploid isolates, resolving to a euploid state within a few passages. Furthermore, the likelihood that a particular isolate will resolve to the same ploidy state in replicate evolution experiments is only ∼50%, supporting the idea that the chromosome loss process of the parasexual cycle is random and does not follow trajectories involving specific combinations of chromosomes. Together, our results indicate that tetraploid progenitors can produce populations of progeny cells with a high degree of genomic diversity, from altered ploidy to homozygosis, providing an excellent source of genetic variation upon which selection can act.