Regional centromeres in the yeast Candida lusitaniae lack pericentromeric heterochromatin

Development of a Shuttle Vector That Transforms at High Frequency for the Emerging Human Fungal Pathogen: Candida auris

Routine molecular manipulation of any organism is inefficient and difficult without the existence of a plasmid. Although transformation is possible in C. auris, no plasmids are available that can serve as cloning or shuttle vectors. C. auris centromeres have been well characterized but have not been explored further as molecular tools. We tested C. auris centromeric sequences to identify which, if any, could be used to create a plasmid that was stably maintained after transformation. We cloned all seven C. auris centromeric sequences and tested them for transformation frequency and stability. Transformation frequency varied significantly; however, one was found to transform at a very high frequency. A 1.7 Kb subclone of this sequence was used to construct a shuttle vector. The vector was stable with selection and maintained at ~1 copy per cell but could be easily lost when selection was removed, which suggested that the properties of the centromeric sequence were more Autonomously Replicating Sequence (ARS)-like than centromere-like when part of a plasmid. Rescue of this plasmid from transformed C. auris cells into E. coli revealed that it remained intact after the initial C. auris transformation, even when carrying large inserts. The plasmid was found to be able to transform all four clades of C. auris, with varying frequencies. This plasmid is an important new reagent in the C. auris molecular toolbox, which will enhance the investigation of this human fungal pathogen.

Telomere-to-Telomere Genome Assembly of Tibetan Medicinal Mushroom Ganoderma leucocontextum and the First Copia Centromeric Retrotransposon in Macro-Fungi Genome

A complete telomere-to-telomere (T2T) genome has been a longstanding goal in the field of genomic research. By integrating high-coverage and precise long-read sequencing data using multiple assembly strategies, we present here the first T2T gap-free genome assembly of Ganoderma leucocontextum strain GL72, a Tibetan medicinal mushroom. The T2T genome, with a size of 46.69 Mb, consists 13 complete nuclear chromosomes and typical telomeric repeats (CCCTAA)n were detected at both ends of 13 chromosomes. The high mapping rate, uniform genome coverage, a complete BUSCOs of 99.7%, and base accuracy exceeding 99.999% indicate that this assembly represents the highest level of completeness and quality. Regions characterized by distinct structural attributes, including highest Hi-C interaction intensity, high repeat content, decreased gene density, low GC content, and minimal or no transcription levels across all chromosomes may represent potential centromeres. Sequence analysis revealed the first Copia centromeric retrotransposon in macro-fungi genome. Phylogenomic analysis identified that G. leucocontextum and G. tsugae diverged from the other Ganoderma species approximately 9.8-17.9 MYA. The prediction of secondary metabolic clusters confirmed the capability of this fungus to produce a substantial quantity of metabolites. This T2T gap-free genome will contribute to the genomic 'dark matter' elucidation and server as a great reference for genetics, genomics, and evolutionary studies of G. leucocontextum.

Sirtuins in Epigenetic Silencing and Control of Gene Expression in Model and Pathogenic Fungi

Fungi, including yeasts, molds, and mushrooms, proliferate on decaying matter and then adopt quiescent forms once nutrients are depleted. This review explores how fungi use sirtuin deacetylases to sense and respond appropriately to changing nutrients. Because sirtuins are NAD⁺-dependent deacetylases, their activity is sensitive to intracellular NAD⁺ availability. This allows them to transmit information about a cell's metabolic state on to the biological processes they influence. Fungal sirtuins are primarily known to deacetylate histones, repressing transcription and modulating genome stability. Their target genes include those involved in NAD⁺ homeostasis, metabolism, sporulation, secondary metabolite production, and virulence traits of pathogenic fungi. By targeting different genes over evolutionary time, sirtuins serve as rewiring points that allow organisms to evolve novel responses to low NAD⁺ stress by bringing relevant biological processes under the control of sirtuins. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Peering Into Candida albicans Pir Protein Function and Comparative Genomics of the Pir Family

The fungal cell wall, comprised primarily of protein and polymeric carbohydrate, maintains cell structure, provides protection from the environment, and is an important antifungal drug target. Pir proteins (proteins with internal repeats) are linked to cell wall β-1,3-glucan and are best studied in Saccharomyces cerevisiae. Sequential deletion of S. cerevisiae PIR genes produces strains with increasingly notable cell wall damage. However, a true null mutant lacking all five S. cerevisiae PIR genes was never constructed. Because only two PIR genes (PIR1, PIR32) were annotated in the Candida albicans genome, the initial goal of this work was to construct a true Δpir/Δpir null strain in this species. Unexpectedly, the phenotype of the null strain was almost indistinguishable from its parent, leading to the search for other proteins with Pir function. Bioinformatic approaches revealed nine additional C. albicans proteins that share a conserved Pir functional motif (minimally DGQ). Examination of the protein sequences revealed another conserved motif (QFQFD) toward the C-terminal end of each protein. Sequence similarities and presence of the conserved motif(s) were used to identify a set of 75 proteins across 16 fungal species that are proposed here as Pir proteins. The Pir family is greatly expanded in C. albicans and C. dubliniensis compared to other species and the orthologs are known to have specialized function during chlamydospore formation. Predicted Pir structures showed a conserved core of antiparallel beta-sheets and sometimes-extensive loops that contain amino acids with the potential to form linkages to cell wall components. Pir phylogeny demonstrated emergence of specific ortholog groups among the fungal species. Variation in gene expression patterns was noted among the ortholog groups during growth in rich medium. PIR allelic variation was quite limited despite the presence of a repeated sequence in many loci. Results presented here demonstrate that the Pir family is larger than previously recognized and lead to new hypotheses to test to better understand Pir proteins and their role in the fungal cell wall.

The Curious Case of Nonrepetitive Centromeric DNA Sequences in Candida auris and Related Species

2009 saw the first description of Candida auris, a yeast pathogen of humans. C. auris has since grown into a global problem in intensive care settings, where it causes systemic infections in patients with underlying health issues. Recent whole-genome sequencing has discerned five C. auris clades with distinct phenotypic features which display genomic divergence on a DNA sequence and a chromosome structure level. In the absence of sexual reproduction in C. auris, the mechanism(s) behind the rapid genomic evolution of this emerging killer yeast has remained obscure. Yet, one important bit of information about chromosome organization was missing, the identification of the centromeres. In a recent study, Sanyal and coworkers (A. Narayanan, R. N. Vadnala, P. Ganguly, P. Selvakumar, et al., mBio 12:e00905-21, 2021, https://doi.org/10.1128/mBio.00905-21) filled this knowledge gap by mapping the centromeres in C. auris and its close relatives. This represents a major advance in the chromosome biology of the Candida/Clavispora clade.

Stable Positions of Epigenetically Inherited Centromeres in the Emerging Fungal Pathogen Candida auris and Its Relatives

Candida auris is an emerging fungal pathogen that is thermotolerant and often resistant to standard antifungal treatments. To trace its evolutionary history, the Sanyal lab conducted a comparative genomic study focusing on the positions of centromeres in C. auris and eight other species from the Clavispora/Candida clade of yeasts (A. Narayanan et al., mBio 12:e00905-12, 2021). These researchers discovered that these species possess small regional centromeres that are highly stable, having remained in the same syntenic positions for over 100 million years. This stability is remarkable, given the lack of a conserved sequence underlying the centromeres and the relative ease with which other yeasts form neocentromeres. Thus, this work provides an opportunity to investigate the molecular mechanism of centromere inheritance in a genetically tractable and medically important yeast.

Clade-specific chromosomal rearrangements and loss of subtelomeric adhesins in Candida auris

Candida auris is an emerging fungal pathogen of rising concern due to global spread, the ability to cause healthcare-associated outbreaks, and antifungal resistance. Genomic analyses revealed that early contemporaneously detected cases of C. auris were geographically stratified into four major clades. While Clades I, III, and IV are responsible for ongoing outbreaks of invasive and multidrug-resistant infections, Clade II, also termed the East Asian clade, consists primarily of cases of ear infection, is often susceptible to all antifungal drugs, and has not been associated with outbreaks. Here, we generate chromosome-level assemblies of twelve isolates representing the phylogenetic breadth of these four clades and the only isolate described to date from Clade V. This Clade V genome is highly syntenic with those of Clades I, III, and IV, although the sequence is highly divergent from the other clades. Clade II genomes appear highly rearranged, with translocations occurring near GC-poor regions, and large subtelomeric deletions in most chromosomes, resulting in a substantially different karyotype. Rearrangements and deletion lengths vary across Clade II isolates, including two from a single patient, supporting ongoing genome instability. Deleted subtelomeric regions are enriched in Hyr/Iff-like cell-surface proteins, novel candidate cell wall proteins, and an ALS-like adhesin. Cell wall proteins from these families and other drug-related genes show clade-specific signatures of selection in Clades I, III, and IV. Subtelomeric dynamics and the conservation of cell surface proteins in the clades responsible for global outbreaks causing invasive infections suggest an explanation for the different phenotypes observed between clades.

Functional and Comparative Analysis of Centromeres Reveals Clade-Specific Genome Rearrangements in Candida auris and a Chromosome Number Change in Related Species

The thermotolerant multidrug-resistant ascomycete Candida auris rapidly emerged since 2009 causing systemic infections worldwide and simultaneously evolved in different geographical zones. The molecular events that orchestrated this sudden emergence of the killer fungus remain mostly elusive. Here, we identify centromeres in C. auris and related species, using a combined approach of chromatin immunoprecipitation and comparative genomic analyses. We find that C. auris and multiple other species in the Clavispora/Candida clade shared a conserved small regional GC-poor centromere landscape lacking pericentromeres or repeats. Further, a centromere inactivation event led to karyotypic alterations in this species complex. Interspecies genome analysis identified several structural chromosomal changes around centromeres. In addition, centromeres are found to be rapidly evolving loci among the different geographical clades of the same species of C. auris Finally, we reveal an evolutionary trajectory of the unique karyotype associated with clade 2 that consists of the drug-susceptible isolates of C. aurisIMPORTANCECandida auris, the killer fungus, emerged as different geographical clades, exhibiting multidrug resistance and high karyotype plasticity. Chromosomal rearrangements are known to play key roles in the emergence of new species, virulence, and drug resistance in pathogenic fungi. Centromeres, the genomic loci where microtubules attach to separate the sister chromatids during cell division, are known to be hot spots of breaks and downstream rearrangements. We identified the centromeres in C. auris and related species to study their involvement in the evolution and karyotype diversity reported in C. auris We report conserved centromere features in 10 related species and trace the events that occurred at the centromeres during evolution. We reveal a centromere inactivation-mediated chromosome number change in these closely related species. We also observe that one of the geographical clades, the East Asian clade, evolved along a unique trajectory, compared to the other clades and related species.

Implications of the Evolutionary Trajectory of Centromeres in the Fungal Kingdom

Chromosome segregation during the cell cycle is an evolutionarily conserved, fundamental biological process. Dynamic interaction between spindle microtubules and the kinetochore complex that assembles on centromere DNA is required for faithful chromosome segregation. The first artificial minichromosome was constructed by cloning the centromere DNA of the budding yeast Saccharomyces cerevisiae. Since then, centromeres have been identified in >60 fungal species. The DNA sequence and organization of the sequence elements are highly diverse across these fungal centromeres. In this article, we provide a comprehensive view of the evolution of fungal centromeres. Studies of this process facilitated the identification of factors influencing centromere specification, maintenance, and propagation through many generations. Additionally, we discuss the unique features and plasticity of centromeric chromatin and the involvement of centromeres in karyotype evolution. Finally, we discuss the implications of recurrent loss of RNA interference (RNAi) and/or heterochromatin components on the trajectory of the evolution of fungal centromeres and propose the centromere structure of the last common ancestor of three major fungal phyla-Ascomycota, Basidiomycota, and Mucoromycota.

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.

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.

Evolution of Distinct Responses to Low NAD+ Stress by Rewiring the Sir2 Deacetylase Network in Yeasts

Evolutionary adaptation increases the fitness of a species in its environment. It can occur through rewiring of gene regulatory networks, such that an organism responds appropriately to environmental changes. We investigated whether sirtuin deacetylases, which repress transcription and require NAD⁺ for activity, serve as transcriptional rewiring points that facilitate the evolution of potentially adaptive traits. If so, bringing genes under the control of sirtuins could enable organisms to mount appropriate responses to stresses that decrease NAD⁺ levels. To explore how the genomic targets of sirtuins shift over evolutionary time, we compared two yeast species, Saccharomyces cerevisiae and Kluyveromyces lactis, that display differences in cellular metabolism and life cycle timing in response to nutrient availability. We identified sirtuin-regulated genes through a combination of chromatin immunoprecipitation and RNA expression. In both species, regulated genes were associated with NAD⁺ homeostasis, mating, and sporulation, but the specific genes differed. In addition, regulated genes in K. lactis were associated with other processes, including utilization of nonglucose carbon sources, detoxification of arsenic, and production of the siderophore pulcherrimin. Consistent with the species-restricted regulation of these genes, sirtuin deletion affected relevant phenotypes in K. lactis but not S. cerevisiae Finally, sirtuin-regulated gene sets were depleted for broadly conserved genes, consistent with sirtuins regulating processes restricted to a few species. Taken together, these results are consistent with the notion that sirtuins serve as rewiring points that allow species to evolve distinct responses to low NAD⁺ stress.

What makes a centromere?

Centromeres are the eukaryotic chromosomal sites at which the kinetochore forms and attaches to spindle microtubules to orchestrate chromosomal segregation in mitosis and meiosis. Although centromeres are essential for cell division, their sequences are not conserved and evolve rapidly. Centromeres vary dramatically in size and organization. Here we categorize their diversity and explore the evolutionary forces shaping them. Nearly all centromeres favor AT-rich DNA that is gene-free and transcribed at a very low level. Repair of frequent centromere-proximal breaks probably contributes to their rapid sequence evolution. Point centromeres are only ~125 bp and are specified by common protein-binding motifs, whereas short regional centromeres are 1-5 kb, typically have unique sequences, and may have pericentromeric repeats adapted to facilitate centromere clustering. Transposon-rich centromeres are often ~100-300 kb and are favored by RNAi machinery that silences transposons, by suppression of meiotic crossovers at centromeres, and by the ability of some transposons to target centromeres. Megabase-length satellite centromeres arise in plants and animals with asymmetric female meiosis that creates centromere competition, and favors satellite monomers one or two nucleosomes in length that position and stabilize centromeric nucleosomes. Holocentromeres encompass the length of a chromosome and may differ dramatically between mitosis and meiosis. We propose a model in which low level transcription of centromeres facilitates the formation of non-B DNA that specifies centromeres and promotes loading of centromeric nucleosomes.

Loss of centromere function drives karyotype evolution in closely related Malassezia species

Genomic rearrangements associated with speciation often result in variation in chromosome number among closely related species. Malassezia species show variable karyotypes ranging between six and nine chromosomes. Here, we experimentally identified all eight centromeres in M. sympodialis as 3-5-kb long kinetochore-bound regions that span an AT-rich core and are depleted of the canonical histone H3. Centromeres of similar sequence features were identified as CENP-A-rich regions in Malassezia furfur, which has seven chromosomes, and histone H3 depleted regions in Malassezia slooffiae and Malassezia globosa with nine chromosomes each. Analysis of synteny conservation across centromeres with newly generated chromosome-level genome assemblies suggests two distinct mechanisms of chromosome number reduction from an inferred nine-chromosome ancestral state: (a) chromosome breakage followed by loss of centromere DNA and (b) centromere inactivation accompanied by changes in DNA sequence following chromosome-chromosome fusion. We propose that AT-rich centromeres drive karyotype diversity in the Malassezia species complex through breakage and inactivation.

Comparative Genomics for the Elucidation of Multidrug Resistance in Candida lusitaniae

Antifungal resistance is an inevitable phenomenon when fungal pathogens are exposed to antifungal drugs. These drugs can be grouped in four distinct classes (azoles, candins, polyenes, and pyrimidine analogs) and are used in different clinical settings. Failures in therapy implicate the sequential or combined use of these different drug classes, which can result in some cases in the development of multidrug resistance (MDR). MDR is particularly challenging in the clinic since it drastically reduces possible treatment alternatives. In this study, we report the rapid development of MDR in Candida lusitaniae in a patient, which became resistant to all known antifungal agents used until now in medicine. To understand how MDR developed in C. lusitaniae, whole-genome sequencing followed by comparative genome analysis was undertaken in sequential MDR isolates. This helped to detect all specific mutations linked to drug resistance and explained the different MDR patterns exhibited by the clinical isolates.

Multidrug resistance (MDR) has emerged in hospitals due to the use of several agents administered in combination or sequentially to the same individual. We reported earlier MDR in Candida lusitaniae during therapy with amphotericin B (AmB), azoles, and candins. Here, we used comparative genomic approaches between the initial susceptible isolate and 4 other isolates with different MDR profiles. From a total of 18 nonsynonymous single nucleotide polymorphisms (NSS) in genome comparisons with the initial isolate, six could be associated with MDR. One of the single nucleotide polymorphisms (SNPs) occurred in a putative transcriptional activator (MRR1) resulting in a V668G substitution in isolates resistant to azoles and 5-fluorocytosine (5-FC). We demonstrated by genome editing that MRR1 acted by upregulation of MFS7 (a multidrug transporter) in the presence of the V668G substitution. MFS7 itself mediated not only azole resistance but also 5-FC resistance, which represents a novel resistance mechanism for this drug class. Three other distinct NSS occurred in FKS1 (a glucan synthase gene that is targeted by candins) in three candin-resistant isolates. Last, two other NSS in ERG3 and ERG4 (ergosterol biosynthesis) resulting in nonsense mutations were revealed in AmB-resistant isolates, one of which accumulated the two ERG NSS. AmB-resistant isolates lacked ergosterol and exhibited sterol profiles, consistent with ERG3 and ERG4 defects. In conclusion, this genome analysis combined with genetics and metabolomics helped decipher the resistance profiles identified in this clinical case. MDR isolates accumulated six different mutations conferring resistance to all antifungal agents used in medicine. This case study illustrates the capacity of C. lusitaniae to rapidly adapt under drug pressure within the host.

Early Diverging Fungus Mucor circinelloides Lacks Centromeric Histone CENP-A and Displays a Mosaic of Point and Regional Centromeres

Centromeres are rapidly evolving across eukaryotes, despite performing a conserved function to ensure high-fidelity chromosome segregation. CENP-A chromatin is a hallmark of a functional centromere in most organisms. Due to its critical role in kinetochore architecture, the loss of CENP-A is tolerated in only a few organisms, many of which possess holocentric chromosomes. Here, we characterize the consequence of the loss of CENP-A in the fungal kingdom. Mucor circinelloides, an opportunistic human pathogen, lacks CENP-A along with the evolutionarily conserved CENP-C but assembles a monocentric chromosome with a localized kinetochore complex throughout the cell cycle. Mis12 and Dsn1, two conserved kinetochore proteins, were found to co-localize to a short region, one in each of nine large scaffolds, composed of an ∼200-bp AT-rich sequence followed by a centromere-specific conserved motif that echoes the structure of budding yeast point centromeres. Resembling fungal regional centromeres, these core centromere regions are embedded in large genomic expanses devoid of genes yet marked by Grem-LINE1s, a novel retrotransposable element silenced by the Dicer-dependent RNAi pathway. Our results suggest that these hybrid features of point and regional centromeres arose from the absence of CENP-A, thus defining novel mosaic centromeres in this early-diverging fungus.

Cis- and Trans-chromosomal Interactions Define Pericentric Boundaries in the Absence of Conventional Heterochromatin

The diploid budding yeast Candida albicans harbors unique CENPA-rich 3- to 5-kb regions that form the centromere (CEN) core on each of its eight chromosomes. The epigenetic nature of these CENs does not permit the stabilization of a functional kinetochore on an exogenously introduced CEN plasmid. The flexible nature of such centromeric chromatin is exemplified by the reversible silencing of a transgene upon its integration into the CENPA-bound region. The lack of a conventional heterochromatin machinery and the absence of defined boundaries of CENPA chromatin makes the process of CEN specification in this organism elusive. Additionally, upon native CEN deletion, C. albicans can efficiently activate neocentromeres proximal to the native CEN locus, hinting at the importance of CEN-proximal regions. In this study, we examine this CEN-proximity effect and identify factors for CEN specification in C. albicans We exploit a counterselection assay to isolate cells that can silence a transgene when integrated into the CEN-flanking regions. We show that the frequency of reversible silencing of the transgene decreases from the central core of CEN7 to its peripheral regions. Using publicly available C. albicans high-throughput chromosome conformation capture data, we identify a 25-kb region centering on the CENPA-bound core that acts as CEN-flanking compact chromatin (CFCC). Cis- and trans-chromosomal interactions associated with the CFCC spatially segregates it from bulk chromatin. We further show that neocentromere activation on chromosome 7 occurs within this specified region. Hence, this study identifies a specialized CEN-proximal domain that specifies and restricts the centromeric activity to a unique region.

Cellular Dynamics and Genomic Identity of Centromeres in Cereal Blast Fungus

Magnaporthe oryzae is an important fungal pathogen that causes a loss of 10% to 30% of the annual rice crop due to the devastating blast disease. In most organisms, kinetochores are clustered together or arranged at the metaphase plate to facilitate synchronized anaphase separation of sister chromatids in mitosis. In this study, we showed that the initially clustered kinetochores separate and position randomly prior to anaphase in M. oryzae. Centromeres in M. oryzae occupy large genomic regions and form on AT-rich DNA without any common sequence motifs. Overall, this study identified atypical kinetochore dynamics and mapped functional centromeres in M. oryzae to define the roles of centromeric and pericentric boundaries in kinetochore assembly on epigenetically specified centromere loci. This study should pave the way for further understanding of the contribution of heterochromatin in genome stability and virulence of the blast fungus and its related species of high economic importance.

Precise kinetochore-microtubule interactions ensure faithful chromosome segregation in eukaryotes. Centromeres, identified as scaffolding sites for kinetochore assembly, are among the most rapidly evolving chromosomal loci in terms of the DNA sequence and length and organization of intrinsic elements. Neither the centromere structure nor the kinetochore dynamics is well studied in plant-pathogenic fungi. Here, we sought to understand the process of chromosome segregation in the rice blast fungus Magnaporthe oryzae. High-resolution imaging of green fluorescent protein (GFP)-tagged inner kinetochore proteins CenpA and CenpC revealed unusual albeit transient declustering of centromeres just before anaphase separation of chromosomes in M. oryzae. Strikingly, the declustered centromeres positioned randomly at the spindle midzone without an apparent metaphase plate per se. Using CenpA chromatin immunoprecipitation followed by deep sequencing, all seven centromeres in M. oryzae were found to be regional, spanning 57-kb to 109-kb transcriptionally poor regions. Highly AT-rich and heavily methylated DNA sequences were the only common defining features of all the centromeres in rice blast. Lack of centromere-specific DNA sequence motifs or repetitive elements suggests an epigenetic specification of centromere function in M. oryzae. PacBio genome assemblies and synteny analyses facilitated comparison of the centromeric/pericentromeric regions in distinct isolates of rice blast and wheat blast and in Magnaporthiopsis poae. Overall, this study revealed unusual centromere dynamics and precisely identified the centromere loci in the top model fungal pathogens that belong to Magnaporthales and cause severe losses in the global production of food crops and turf grasses.

Transcriptional silencing of centromere repeats by heterochromatin safeguards chromosome integrity

The centromere region of chromosomes consists of repetitive DNA sequences, and is, therefore, one of the fragile sites of chromosomes in many eukaryotes. In the core region, the histone H3 variant CENP-A forms centromere-specific nucleosomes that are required for kinetochore formation. In the pericentromeric region, histone H3 is methylated at lysine 9 (H3K9) and heterochromatin is formed. The transcription of pericentromeric repeats by RNA polymerase II is strictly repressed by heterochromatin. However, the role of the transcriptional silencing of the pericentromeric repeats remains largely unclear. Here, we focus on the chromosomal rearrangements that occur at the repetitive centromeres, and highlight our recent studies showing that transcriptional silencing by heterochromatin suppresses gross chromosomal rearrangements (GCRs) at centromeres in fission yeast. Inactivation of the Clr4 methyltransferase, which is essential for the H3K9 methylation, increased GCRs with breakpoints located in centromeric repeats. However, mutations in RNA polymerase II or the transcription factor Tfs1/TFIIS, which promotes restart of RNA polymerase II following its backtracking, reduced the GCRs that occur in the absence of Clr4, demonstrating that heterochromatin suppresses GCRs by repressing the Tfs1-dependent transcription. We also discuss how the transcriptional restart gives rise to chromosomal rearrangements at centromeres.

Heterochromatin suppresses gross chromosomal rearrangements at centromeres by repressing Tfs1/TFIIS-dependent transcription

Heterochromatin, characterized by histone H3 lysine 9 (H3K9) methylation, assembles on repetitive regions including centromeres. Although centromeric heterochromatin is important for correct segregation of chromosomes, its exact role in maintaining centromere integrity remains elusive. Here, we found in fission yeast that heterochromatin suppresses gross chromosomal rearrangements (GCRs) at centromeres. Mutations in Clr4/Suv39 methyltransferase increased the formation of isochromosomes, whose breakpoints were located in centromere repeats. H3K9A and H3K9R mutations also increased GCRs, suggesting that Clr4 suppresses centromeric GCRs via H3K9 methylation. HP1 homologs Swi6 and Chp2 and the RNAi component Chp1 were the chromodomain proteins essential for full suppression of GCRs. Remarkably, mutations in RNA polymerase II (RNAPII) or Tfs1/TFIIS, the transcription factor that facilitates restart of RNAPII after backtracking, specifically bypassed the requirement of Clr4 for suppressing GCRs. These results demonstrate that heterochromatin suppresses GCRs by repressing Tfs1-dependent transcription of centromere repeats.

Identification of an Exceptionally Long Intron in the HAC1 Gene of Candida parapsilosis

The unfolded protein response (UPR) responds to the build-up of misfolded proteins in the endoplasmic reticulum. The UPR has wide-ranging functions from fungal pathogenesis to applications in biotechnology. The UPR is regulated through the splicing of an unconventional intron in the HAC1 gene. This intron has been described in many fungal species and is of variable length. Until now it was believed that some members of the CTG-Ser1 clade such as C. parapsilosis did not contain an intron in HAC1, suggesting that the UPR was regulated in a different manner. Here we demonstrate that HAC1 plays an important role in regulating the UPR in C. parapsilosis. We also identified an unusually long intron (626 bp) in C. parapsilosisHAC1. Further analysis showed that HAC1 orthologs in several species in the CTG-Ser1 clade contain long introns.

The unfolded protein response (UPR) in the endoplasmic reticulum (ER) is well conserved in eukaryotes from metazoa to yeast. The transcription factor HAC1 is a major regulator of the UPR in many eukaryotes. Deleting HAC1 in the yeast Candida parapsilosis rendered cells more sensitive to DTT, a known inducer of the UPR. The deletion strain was also sensitive to Congo red, calcofluor white, and the antifungal drug ketoconazole, indicating that HAC1 has a role in cell wall maintenance. Transcriptomic analysis revealed that treatment of the wild type with DTT resulted in the increased expression of 368 genes. Comparison with mutant cells treated with DTT reveals that expression of 137 of these genes requires HAC1. Enriched GO term analysis includes response to ER stress, cell wall biogenesis and glycosylation. Orthologs of many of these are associated with UPR in Saccharomyces cerevisiae and Candida albicans. Unconventional splicing of an intron from HAC1 mRNA is required to produce a functional transcription factor. The spliced intron varies in length from 19 bases in C. albicans to 379 bases in Candida glabrata, but has not been previously identified in Candida parapsilosis and related species. We used RNA-seq data and in silico analysis to identify the HAC1 intron in 12 species in the CTG-Ser1 clade. We show that the intron has undergone major contractions and expansions in this clade, reaching up to 848 bases. Exposure to DTT induced splicing of the long intron in C. parapsilosisHAC1, inducing the UPR.

IMPORTANCE The unfolded protein response (UPR) responds to the build-up of misfolded proteins in the endoplasmic reticulum. The UPR has wide-ranging functions from fungal pathogenesis to applications in biotechnology. The UPR is regulated through the splicing of an unconventional intron in the HAC1 gene. This intron has been described in many fungal species and is of variable length. Until now it was believed that some members of the CTG-Ser1 clade such as C. parapsilosis did not contain an intron in HAC1, suggesting that the UPR was regulated in a different manner. Here we demonstrate that HAC1 plays an important role in regulating the UPR in C. parapsilosis. We also identified an unusually long intron (626 bp) in C. parapsilosisHAC1. Further analysis showed that HAC1 orthologs in several species in the CTG-Ser1 clade contain long introns.

RNAi is a critical determinant of centromere evolution in closely related fungi

The “centromere paradox” refers to rapidly evolving and highly diverse centromere DNA sequences even in closely related eukaryotes. However, factors contributing to this rapid divergence are largely unknown. Here, we identified large regional, LTR retrotransposon-rich centromeres in a group of human fungal pathogens belonging to the Cryptococcus species complex. We provide evidence that loss-of-functional RNAi machinery and possibly cytosine DNA methylation trigger instability of the genome by activation of centromeric retrotransposons presumably suppressed by RNAi. We propose that RNAi, together with cytosine DNA methylation, serves as a critical determinant that maintains repetitive transposon-rich centromere structures. This study explores the direct link between RNAi and centromere structure evolution.

The centromere DNA locus on a eukaryotic chromosome facilitates faithful chromosome segregation. Despite performing such a conserved function, centromere DNA sequence as well as the organization of sequence elements is rapidly evolving in all forms of eukaryotes. The driving force that facilitates centromere evolution remains an enigma. Here, we studied the evolution of centromeres in closely related species in the fungal phylum of Basidiomycota. Using ChIP-seq analysis of conserved inner kinetochore proteins, we identified centromeres in three closely related Cryptococcus species: two of which are RNAi-proficient, while the other lost functional RNAi. We find that the centromeres in the RNAi-deficient species are significantly shorter than those of the two RNAi-proficient species. While centromeres are LTR retrotransposon-rich in all cases, the RNAi-deficient species lost all full-length retroelements from its centromeres. In addition, centromeres in RNAi-proficient species are associated with a significantly higher level of cytosine DNA modifications compared with those of RNAi-deficient species. Furthermore, when an RNAi-proficient Cryptococcus species and its RNAi-deficient mutants were passaged under similar conditions, the centromere length was found to be occasionally shortened in RNAi mutants. In silico analysis of predicted centromeres in a group of closely related Ustilago species, also belonging to the Basidiomycota, were found to have undergone a similar transition in the centromere length in an RNAi-dependent fashion. Based on the correlation found in two independent basidiomycetous species complexes, we present evidence suggesting that the loss of RNAi and cytosine DNA methylation triggered transposon attrition, which resulted in shortening of centromere length during evolution.

Evolving Centromeres and Kinetochores

The genetic material, contained on chromosomes, is often described as the "blueprint for life." During nuclear division, the chromosomes are pulled into each of the two daughter nuclei by the coordination of spindle microtubules, kinetochores, centromeres, and chromatin. These four functional units must link the chromosomes to the microtubules, signal to the cell when the attachment is made so that division can proceed, and withstand the force generated by pulling the chromosomes to either daughter cell. To perform each of these functions, kinetochores are large protein complexes, approximately 5MDa in size, and they contain at least 45 unique proteins. Many of the central components in the kinetochore are well conserved, yielding a common core of proteins forming consistent structures. However, many of the peripheral subcomplexes vary between different taxonomic groups, including changes in primary sequence and gain or loss of whole proteins. It is still unclear how significant these changes are, and answers to this question may provide insights into adaptation to specific lifestyles or progression of disease that involve chromosome instability.

Proteogenomics produces comprehensive and highly accurate protein-coding gene annotation in a complete genome assembly of Malassezia sympodialis

Complete and accurate genome assembly and annotation is a crucial foundation for comparative and functional genomics. Despite this, few complete eukaryotic genomes are available, and genome annotation remains a major challenge. Here, we present a complete genome assembly of the skin commensal yeast Malassezia sympodialis and demonstrate how proteogenomics can substantially improve gene annotation. Through long-read DNA sequencing, we obtained a gap-free genome assembly for M. sympodialis (ATCC 42132), comprising eight nuclear and one mitochondrial chromosome. We also sequenced and assembled four M. sympodialis clinical isolates, and showed their value for understanding Malassezia reproduction by confirming four alternative allele combinations at the two mating-type loci. Importantly, we demonstrated how proteomics data could be readily integrated with transcriptomics data in standard annotation tools. This increased the number of annotated protein-coding genes by 14% (from 3612 to 4113), compared to using transcriptomics evidence alone. Manual curation further increased the number of protein-coding genes by 9% (to 4493). All of these genes have RNA-seq evidence and 87% were confirmed by proteomics. The M. sympodialis genome assembly and annotation presented here is at a quality yet achieved only for a few eukaryotic organisms, and constitutes an important reference for future host-microbe interaction studies.

Diatom centromeres suggest a mechanism for nuclear DNA acquisition

Centromeres are essential for cell division and growth in all eukaryotes, and knowledge of their sequence and structure guides the development of artificial chromosomes for functional cellular biology studies. Centromeric proteins are conserved among eukaryotes; however, centromeric DNA sequences are highly variable. We combined forward and reverse genetic approaches with chromatin immunoprecipitation to identify centromeres of the model diatom Phaeodactylum tricornutum We observed 25 unique centromere sequences typically occurring once per chromosome, a finding that helps to resolve nuclear genome organization and indicates monocentric regional centromeres. Diatom centromere sequences contain low-GC content regions but lack repeats or other conserved sequence features. Native and foreign sequences with similar GC content to P. tricornutum centromeres can maintain episomes and recruit the diatom centromeric histone protein CENH3, suggesting nonnative sequences can also function as diatom centromeres. Thus, simple sequence requirements may enable DNA from foreign sources to persist in the nucleus as extrachromosomal episomes, revealing a potential mechanism for organellar and foreign DNA acquisition.

Genome Diversity and Evolution in the Budding Yeasts (Saccharomycotina)

Considerable progress in our understanding of yeast genomes and their evolution has been made over the last decade with the sequencing, analysis, and comparisons of numerous species, strains, or isolates of diverse origins. The role played by yeasts in natural environments as well as in artificial manufactures, combined with the importance of some species as model experimental systems sustained this effort. At the same time, their enormous evolutionary diversity (there are yeast species in every subphylum of Dikarya) sparked curiosity but necessitated further efforts to obtain appropriate reference genomes. Today, yeast genomes have been very informative about basic mechanisms of evolution, speciation, hybridization, domestication, as well as about the molecular machineries underlying them. They are also irreplaceable to investigate in detail the complex relationship between genotypes and phenotypes with both theoretical and practical implications. This review examines these questions at two distinct levels offered by the broad evolutionary range of yeasts: inside the best-studied Saccharomyces species complex, and across the entire and diversified subphylum of Saccharomycotina. While obviously revealing evolutionary histories at different scales, data converge to a remarkably coherent picture in which one can estimate the relative importance of intrinsic genome dynamics, including gene birth and loss, vs. horizontal genetic accidents in the making of populations. The facility with which novel yeast genomes can now be studied, combined with the already numerous available reference genomes, offer privileged perspectives to further examine these fundamental biological questions using yeasts both as eukaryotic models and as fungi of practical importance.

Sporadic Gene Loss After Duplication Is Associated with Functional Divergence of Sirtuin Deacetylases Among Candida Yeast Species

Gene duplication promotes the diversification of protein functions in several ways. Ancestral functions can be partitioned between the paralogs, or a new function can arise in one paralog. These processes are generally viewed as unidirectional. However, paralogous proteins often retain related functions and can substitute for one another. Moreover, in the event of gene loss, the remaining paralog might regain ancestral functions that had been shed. To explore this possibility, we focused on the sirtuin deacetylase SIR2 and its homolog HST1 in the CTG clade of yeasts. HST1 has been consistently retained throughout the clade, whereas SIR2 is only present in a subset of species. These NAD⁺-dependent deacetylases generate condensed chromatin that represses transcription and stabilizes tandemly repeated sequences. By analyzing phylogenetic trees and gene order, we found that a single duplication of the SIR2/HST1 gene occurred, likely prior to the emergence of the CTG clade. This ancient duplication was followed by at least two independent losses of SIR2. Functional characterization of Sir2 and Hst1 in three species revealed that these proteins have not maintained consistent functions since the duplication. In particular, the rDNA locus is deacetylated by Sir2 in Candida albicans, by Hst1 in C. lusitaniae, and by neither paralog in C. parapsilosis. In addition, the subtelomeres in C. albicans are deacetylated by Sir2 rather than by Hst1, which is orthologous to the sirtuin associated with Saccharomyces cerevisiae subtelomeres. These differences in function support the model that sirtuin deacetylases can regain ancestral functions to compensate for gene loss.

Centromeres of the Yeast Komagataella phaffii (Pichia pastoris) Have a Simple Inverted-Repeat Structure

Centromere organization has evolved dramatically in one clade of fungi, the Saccharomycotina. These yeasts have lost the ability to make normal eukaryotic heterochromatin with histone H3K9 methylation, which is a major component of pericentromeric regions in other eukaryotes. Following this loss, several different types of centromere emerged, including two types of sequence-defined (“point”) centromeres, and the epigenetically defined “small regional” centromeres of Candida albicans. Here we report that centromeres of the methylotrophic yeast Komagataella phaffii (formerly called Pichia pastoris) are structurally defined. Each of its four centromeres consists of a 2-kb inverted repeat (IR) flanking a 1-kb central core (mid) region. The four centromeres are unrelated in sequence. CenH3 (Cse4) binds strongly to the cores, with a decreasing gradient along the IRs. This mode of organization resembles Schizosaccharomyces pombe centromeres but is much more compact and lacks the extensive flanking heterochromatic otr repeats. Different isolates of K. phaffii show polymorphism for the orientation of the mid regions, due to recombination in the IRs. CEN4 is located within a 138-kb region that changes orientation during mating-type switching, but switching does not induce recombination of centromeric IRs. Our results demonstrate that evolutionary transitions in centromere organization have occurred in multiple yeast clades.

Histone Deacetylases and Their Inhibition in Candida Species

Fungi are generally benign members of the human mucosal flora or live as saprophytes in the environment. However, they can become pathogenic, leading to invasive and life threatening infections in vulnerable patients. These invasive fungal infections are regarded as a major public health problem on a similar scale to tuberculosis or malaria. Current treatment for these infections is based on only four available drug classes. This limited therapeutic arsenal and the emergence of drug-resistant strains are a matter of concern due to the growing number of patients to be treated, and new therapeutic strategies are urgently needed. Adaptation of fungi to drug pressure involves transcriptional regulation, in which chromatin dynamics and histone modifications play a major role. Histone deacetylases (HDACs) remove acetyl groups from histones and actively participate in controlling stress responses. HDAC inhibition has been shown to limit fungal development, virulence, biofilm formation, and dissemination in the infected host, while also improving the efficacy of existing antifungal drugs toward Candida spp. In this article, we review the functional roles of HDACs and the biological effects of HDAC inhibitors on Candida spp., highlighting the correlations between their pathogenic effects in vitro and in vivo. We focus on how HDAC inhibitors could be used to treat invasive candidiasis while also reviewing recent developments in their clinical evaluation.

The Chromatin of Candida albicans Pericentromeres Bears Features of Both Euchromatin and Heterochromatin

Centromeres, sites of kinetochore assembly, are important for chromosome stability and integrity. Most eukaryotes have regional centromeres epigenetically specified by the presence of the histone H3 variant CENP-A. CENP-A chromatin is often surrounded by pericentromeric regions packaged into transcriptionally silent heterochromatin. Candida albicans, the most common human fungal pathogen, possesses small regional centromeres assembled into CENP-A chromatin. The chromatin state of C. albicans pericentromeric regions is unknown. Here, for the first time, we address this question. We find that C. albicans pericentromeres are assembled into an intermediate chromatin state bearing features of both euchromatin and heterochromatin. Pericentromeric chromatin is associated with nucleosomes that are highly acetylated, as found in euchromatic regions of the genome; and hypomethylated on H3K4, as found in heterochromatin. This intermediate chromatin state is inhibitory to transcription and partially represses expression of proximal genes and inserted marker genes. Our analysis identifies a new chromatin state associated with pericentromeric regions.

Repeat-Associated Fission Yeast-Like Regional Centromeres in the Ascomycetous Budding Yeast Candida tropicalis

The centromere, on which kinetochore proteins assemble, ensures precise chromosome segregation. Centromeres are largely specified by the histone H3 variant CENP-A (also known as Cse4 in yeasts). Structurally, centromere DNA sequences are highly diverse in nature. However, the evolutionary consequence of these structural diversities on de novo CENP-A chromatin formation remains elusive. Here, we report the identification of centromeres, as the binding sites of four evolutionarily conserved kinetochore proteins, in the human pathogenic budding yeast Candida tropicalis. Each of the seven centromeres comprises a 2 to 5 kb non-repetitive mid core flanked by 2 to 5 kb inverted repeats. The repeat-associated centromeres of C. tropicalis all share a high degree of sequence conservation with each other and are strikingly diverged from the unique and mostly non-repetitive centromeres of related Candida species--Candida albicans, Candida dubliniensis, and Candida lusitaniae. Using a plasmid-based assay, we further demonstrate that pericentric inverted repeats and the underlying DNA sequence provide a structural determinant in CENP-A recruitment in C. tropicalis, as opposed to epigenetically regulated CENP-A loading at centromeres in C. albicans. Thus, the centromere structure and its influence on de novo CENP-A recruitment has been significantly rewired in closely related Candida species. Strikingly, the centromere structural properties along with role of pericentric repeats in de novo CENP-A loading in C. tropicalis are more reminiscent to those of the distantly related fission yeast Schizosaccharomyces pombe. Taken together, we demonstrate, for the first time, fission yeast-like repeat-associated centromeres in an ascomycetous budding yeast.