Genome plasticity in Candida albicans is driven by long repeat sequences

Positive selection in cilia-related genes may facilitate deep-sea adaptation of Thermocollonia jamsteci

Deep-sea hydrothermal vents are characterized by high hydrostatic pressure, hypoxia, darkness and toxic substances. However, how organisms adapt to such extreme marine ecosystems remain poorly understood. We hypothesize that adaptive evolution plays an essential role in generating novelty for evolutionary adaptation to the deep-sea environment because adaptive evolution has been found to be critical for species origin and evolution. In this project, the chromosome-level genome of the deep-sea hydrothermal vent gastropod T. jamsteci was constructed for the first time to examine molecular mechanisms of its adaptation to the deep-sea environment. The genome size was large (2.54 Gb), ranking at the top of all species in the Vetigastropoda subclass, driven primarily by the bursts of transposable elements (TEs). The transposition of TEs may also trigger chromosomal changes including both inter-chromosomal fusions and intra-chromosomal activities involving chromosome inversions, rearrangements and fusions, as revealed by comparing the genomes of T. jamsteci and its closely related shallow-sea species Gibbula magus. Innovative changes including the expansion of the ABC transporter gene family that may facilitate detoxification, duplication of genes related to endocytosis, immunity, apoptosis, and anti-apoptotic domains that may help T. jamsteci fight against microbial pathogens, were identified. Furthermore, comparative analysis identified positive selection signals in a large number of genes including the hypoxia up-regulated protein 1, which is a chaperone that may promote adaptation of the T. jamsteci to hypoxic deepsea environments, hox2, Rx2, Pax6 and cilia-related genes BBS1, BBS2, BBS9 and RFX4. Notably, because of the critical importance of cilia and IFT in development, positive selection in cilia-related genes may play a critical role in facilitating T. jamsteci to adapt to the high-pressure deep-sea ecosystem. Results from this study thus revealed important molecular clues that may facilitate further research on the adaptation of molluscs to deep-sea hydrothermal vents.

Evolution and strain diversity advance exploration of Candida albicans biology

Fungi were some of the earliest organismal systems used to explore mutational processes and its phenotypic consequences on members of a species. Yeasts that cause significant human disease were quickly incorporated into these investigations to define the genetic and phenotypic drivers of virulence. Among Candida species, Candida albicans has emerged as a model for studying genomic processes of evolution because of its clinical relevance, relatively small genome, and ability to tolerate complex chromosomal changes. Here, we describe major recent findings that used evolution of strains from defined genetic backgrounds to delineate mutational and adaptative processes and include how nascent exploration into naturally occurring variation is contributing to these conceptual frameworks. Ultimately, efforts to discern adaptive mechanisms used by C. albicans will continue to divulge new biology and can better inform treatment regimens for the increasing prevalence of fungal disease.

Long-term stability of acquired drug resistance and resistance associated mutations in the fungal pathogen Nakaseomyces glabratus (Candida glabrata)

The limited number of available antifungal drugs and the increasing number of fungal isolates that show drug or multidrug resistance pose a serious medical threat. Several yeast pathogens, such as Nakaseomyces glabratus (Candida glabrata), show a remarkable ability to develop drug resistance during treatment through the acquisition of genetic mutations. However, how stable this resistance and the underlying mutations are in non-selective conditions remains poorly characterized. The stability of acquired drug resistance has fundamental implications for our understanding of the appearance and spread of drug-resistant outbreaks and for defining efficient strategies to combat them. Here, we used an in vitro evolution approach to assess the stability under optimal growth conditions of resistance phenotypes and resistance-associated mutations that were previously acquired under exposure to antifungals. Our results reveal a remarkable stability of the resistant phenotype and the underlying mutations in a significant number of evolved populations, which conserved their phenotype for at least two months in the absence of drug-selective pressure. We observed a higher stability of anidulafungin resistance over fluconazole resistance, and of resistance-conferring point mutations as compared with aneuploidies. In addition, we detected accumulation of novel mutations in previously altered resistance-associated genes in non-selective conditions, which suggest a possible compensatory role. We conclude that acquired resistance, particularly to anidulafungin, is a long-lasting phenotype, which has important implications for the persistence and propagation of drug-resistant clinical outbreaks.

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

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

The genomic landscape of transposable elements in yeast hybrids is shaped by structural variation and genotype-specific modulation of transposition rate

Transposable elements (TEs) are major contributors to structural genomic variation by creating interspersed duplications of themselves. In return, structural variants (SVs) can affect the genomic distribution of TE copies and shape their load. One long-standing hypothesis states that hybridization could trigger TE mobilization and thus increase TE load in hybrids. We previously tested this hypothesis (Hénault et al., 2020) by performing a large-scale evolution experiment by mutation accumulation (MA) on multiple hybrid genotypes within and between wild populations of the yeasts Saccharomyces paradoxus and Saccharomyces cerevisiae. Using aggregate measures of TE load with short-read sequencing, we found no evidence for TE load increase in hybrid MA lines. Here, we resolve the genomes of the hybrid MA lines with long-read phasing and assembly to precisely characterize the role of SVs in shaping the TE landscape. Highly contiguous phased assemblies of 127 MA lines revealed that SV types like polyploidy, aneuploidy, and loss of heterozygosity have large impacts on the TE load. We characterized 18 de novo TE insertions, indicating that transposition only has a minor role in shaping the TE landscape in MA lines. Because the scarcity of TE mobilization in MA lines provided insufficient resolution to confidently dissect transposition rate variation in hybrids, we adapted an in vivo assay to measure transposition rates in various S. paradoxus hybrid backgrounds. We found that transposition rates are not increased by hybridization, but are modulated by many genotype-specific factors including initial TE load, TE sequence variants, and mitochondrial DNA inheritance. Our results show the multiple scales at which TE load is shaped in hybrid genomes, being highly impacted by SV dynamics and finely modulated by genotype-specific variation in transposition rates.

Tandem gene duplications contributed to high-level azole resistance in a rapidly expanding Candida tropicalis population

Invasive diseases caused by the globally distributed commensal yeast Candida tropicalis are associated with mortality rates of greater than 50%. Notable increases of azole resistance have been observed in this species, particularly within Asia-Pacific regions. Here, we carried out a genetic population study on 1571 global C. tropicalis isolates using multilocus sequence typing (MLST). In addition, whole-genome sequencing (WGS) analysis was conducted on 629 of these strains, comprising 448 clinical invasive strains obtained in this study and 181 genomes sourced from public databases. We found that MLST clade 4 is the predominant azole-resistant clone. WGS analyses demonstrated that dramatically increasing rates of azole resistance are associated with a rapid expansion of cluster AZR, a sublineage of clade 4. Cluster AZR isolates exhibited a distinct high-level azole resistance, which was induced by tandem duplications of the ERG11A395T gene allele. Ty3/gypsy-like retrotransposons were found to be highly enriched in this population. The alarming expansion of C. tropicalis cluster AZR population underscores the urgent need for strategies against growing threats of antifungal resistance.

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

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

The Microevolution of Antifungal Drug Resistance in Pathogenic Fungi

The mortality rates of invasive fungal infections remain high because of the limited number of antifungal drugs available and antifungal drug resistance, which can rapidly evolve during treatment. Mutations in key resistance genes such as ERG11 were postulated to be the predominant cause of antifungal drug resistance in the clinic. However, recent advances in whole genome sequencing have revealed that there are multiple mechanisms leading to the microevolution of resistance. In many fungal species, resistance can emerge through ERG11-independent mechanisms and through the accumulation of mutations in many genes to generate a polygenic resistance phenotype. In addition, genome sequencing has revealed that full or partial aneuploidy commonly occurs in clinical or microevolved in vitro isolates to confer antifungal resistance. This review will provide an overview of the mutations known to be selected during the adaptive microevolution of antifungal drug resistance and focus on how recent advances in genome sequencing technology have enhanced our understanding of this process.

Recurrent copy number variations in the human fungal pathogen Candida parapsilosis

Candida parapsilosis is an opportunistic fungal pathogen with increasing incidence in hospital settings worldwide; however, we lack a comprehensive understanding of the mechanisms promoting its virulence and drug resistance. Bergin et al. systematically quantify the frequency and effect of copy number variation (CNV) across 170 diverse clinical and environmental isolates of C. parapsilosis (Bergin SA, Zhao F, Ryan AP, Müller CA, Nieduszynski CA, Zhai B, Rolling T, Hohl TM, Morio F, Scully J, Wolfe KH, Butler G, 2022, mBio, https://doi.org/10.1128/mbio.01777-22). Using a combination of both short- and long-read whole genome sequencing techniques, they determine the structure and copy number of two CNVs that arose recurrently throughout the evolution of these isolates. Each CNV predominantly amplifies one coding sequence (ARR3 or RTA3); however, the amplitude and recombination breakpoints are variable across the isolates. Amplification of RTA3 correlates with drug resistance and deletion causes drug susceptibility. This study highlights the need for further research into the mechanisms and dynamics of CNV formation and the impact of these CNVs on virulence and drug resistance across diverse fungal pathogens.

The First Genome Survey of the Snail Provanna glabra Inhabiting Deep-Sea Hydrothermal Vents

The snail P. glabra is an endemic species in deep-sea chemosynthetic ecosystems of the Northwest Pacific Ocean. To obtain more genetic information on this species and provide the basis for subsequent whole-genome map construction, a genome survey was performed on this snail from the hydrothermal vent of Okinawa Trough. The genomic size of P. glabra was estimated to be 1.44 Gb, with a heterozygosity of 1.91% and a repeated sequence content of 69.80%. Based on the sequencing data, a draft genome of 1.32 Gb was assembled. Transposal elements (TEs) accounted for 40.17% of the entire genome, with DNA transposons taking the highest proportion. It was found that most TEs were inserted in the genome recently. In the simple sequence repeats, the dinucleotide motif was the most enriched microsatellite type, accounting for 53% of microsatellites. A complete mitochondrial genome of P. glabra with a total length of 16,268 bp was assembled from the sequencing data. After comparison with the published mitochondrial genome of Provanna sp. from a methane seep, 331 potential single nucleotide polymorphism (SNP) sites were identified in protein-coding genes (PCGs). Except for the cox1 gene, nad2, nad4, nad5, and cob genes are expected to be candidate markers for population genetic and phylogenetic studies of P. glabra and other deep-sea snails. Compared with shallow-water species, three mitochondrial genes of deep-sea gastropods exhibited a higher evolutionary rate, indicating strong selection operating on mitochondria of deep-sea species. This study provides insights into the genome characteristics of P. glabra and supplies genomic resources for further studies on the adaptive evolution of the snail in extreme deep-sea chemosynthetic environments.

Genetic adaptations of sea anemone to hydrothermal environment

Hydrothermal vent habitats are characterized by high hydrostatic pressure, darkness, and the continuous release of toxic metal ions into the surrounding environment where sea anemones and other invertebrates thrive. Nevertheless, the understanding of metazoan metal ion tolerances and environmental adaptations remains limited. We assembled a chromosome-level genome for the vent sea anemone, Alvinactis idsseensis sp. nov. Comparative genomic analyses revealed gene family expansions and gene innovations in A. idsseensis sp. nov. as a response to high concentrations of metal ions. Impressively, the metal tolerance proteins MTPs is a unique evolutionary response to the high concentrations of Fe2+ and Mn2+ present in the environments of these anemones. We also found genes associated with high concentrations of polyunsaturated fatty acids that may respond to high hydrostatic pressure and found sensory and circadian rhythm-regulated genes that were essential for adaptations to darkness. Overall, our results provide insights into metazoan adaptation to metal ions, high pressure, and darkness in hydrothermal vents.

The Dynamic Fungal Genome: Polyploidy, Aneuploidy and Copy Number Variation in Response to Stress

Fungal species have dynamic genomes and often exhibit genomic plasticity in response to stress. This genome plasticity often comes with phenotypic consequences that affect fitness and resistance to stress. Fungal pathogens exhibit genome plasticity in both clinical and agricultural settings and often during adaptation to antifungal drugs, posing significant challenges to human health. Therefore, it is important to understand the rates, mechanisms, and impact of large genomic changes. This review addresses the prevalence of polyploidy, aneuploidy, and copy number variation across diverse fungal species, with special attention to prominent fungal pathogens and model species. We also explore the relationship between environmental stress and rates of genomic changes and highlight the mechanisms underlying genotypic and phenotypic changes. A comprehensive understanding of these dynamic fungal genomes is needed to identify novel solutions for the increase in antifungal drug resistance.

Confronting antifungal resistance, tolerance, and persistence: Advances in drug target discovery and delivery systems

Human pathogenic fungi pose a serious threat to human health and safety. Unfortunately, the limited number of antifungal options is exacerbated by the continuous emergence of drug-resistant variants, leading to frequent drug treatment failures. Recent studies have also highlighted the clinical importance of other modes of fungal survival of antifungal treatment, including drug tolerance and persistence, pointing to the complexity of the fungal response to antifungal drugs. A lack of understanding of the fungal drug response has hampered the identification of new targets, the development of alternative antifungal strategies and the design of appropriate delivery systems. In this review we summarize recent advances in the study of antifungal resistance, tolerance and persistence, with an emphasis on promising drug targets and drug delivery systems that may yield important insights into the development of new or improved antifungal therapies against fungal infections.

Pattern and variation in simple sequence repeat (SSR) at different genomic regions and its implications to maize evolution and breeding

BACKGROUND

Repetitive DNA sequences accounts for over 80% of maize genome. Although simple sequence repeats (SSRs) account for only 0.03% of the genome, they have been widely used in maize genetic research and breeding as highly informative codominant DNA markers. The genome-wide distribution and polymorphism of SSRs are not well studied due to the lack of high-quality genome DNA sequence data.

RESULTS

In this study, using data from high-quality de novo-sequenced maize genomes of five representative maize inbred lines, we revealed that SSRs were more densely present in telomeric region than centromeric region, and were more abundant in genic sequences than intergenic sequences. On genic sequences, tri- and hexanucleotide motifs were more abundant in CDS sequence and some mono- and dinucleotide motifs were more abundant in UTR sequences. Median length and chromosomal density of SSRs were both narrowly range-bound, with median length of 14-18 bp and genome-wide average density of 3355.77 bp/Mbp. LTR-RTs of < 0.4 Mya had higher SSR density (4498-4992 bp/Mbp). The genome-specific and motif-specific SSR polymorphism were studied. Their potential breeding applications were discussed.

CONCLUSIONS

We found that the median length of SSR sequences of different SSR motifs was nearly constant. SSR density in genic regions was much higher than intergenic regions. In addition, SSR density at LTR-RTs of different evolutionary ages varied in a narrow range. The SSRs and their LTR-RT carriers evolved at an equal rate. All these observations indicated that SSR length and density were under control of yet unknown evolutionary forces. The chromosome region-specific and motif-specific SSR polymorphisms we observed supported the notion that SSR polymorphism was invaluable genome resource for developing highly informative genome and gene markers in maize genetic research and molecular breeding.

Development and applications of a CRISPR activation system for facile genetic overexpression in Candida albicans

For the fungal pathogen Candida albicans, genetic overexpression readily occurs via a diversity of genomic alterations, such as aneuploidy and gain-of-function mutations, with important consequences for host adaptation, virulence, and evolution of antifungal drug resistance. Given the important role of overexpression on C. albicans biology, it is critical to develop and harness tools that enable the analysis of genes expressed at high levels in the fungal cell. Here, we describe the development, optimization, and application of a novel, single-plasmid-based CRISPR activation (CRISPRa) platform for targeted genetic overexpression in C. albicans, which employs a guide RNA to target an activator complex to the promoter region of a gene of interest, thus driving transcriptional expression of that gene. Using this system, we demonstrate the ability of CRISPRa to drive high levels of gene expression in C. albicans, and we assess optimal guide RNA targeting for robust and constitutive overexpression. We further demonstrate the specificity of the system via RNA sequencing. We highlight the application of CRISPR activation to overexpress genes involved in pathogenesis and drug susceptibility, and contribute toward the identification of novel phenotypes. Consequently, this tool will facilitate a broad range of applications for the study of C. albicans genetic overexpression.

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

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.

Stress combined with loss of the Candida albicans SUMO protease Ulp2 triggers selection of aneuploidy via a two-step process

A delicate balance between genome stability and instability ensures genome integrity while generating genetic diversity, a critical step for evolution. Indeed, while excessive genome instability is harmful, moderated genome instability can drive adaptation to novel environments by maximising genetic variation. Candida albicans, a human fungal pathogen that colonises different parts of the human body, adapts rapidly and frequently to different hostile host microenvironments. In this organism, the ability to generate large-scale genomic variation is a key adaptative mechanism triggering dangerous infections even in the presence of antifungal drugs. Understanding how fitter novel karyotypes are selected is key to determining how C. albicans and other microbial pathogens establish infections. Here, we identified the SUMO protease Ulp2 as a regulator of C. albicans genome integrity through genetic screening. Deletion of ULP2 leads to increased genome instability, enhanced genome variation and reduced fitness in the absence of additional stress. The combined stress caused by the lack of ULP2 and antifungal drug treatment leads to the selection of adaptive segmental aneuploidies that partially rescue the fitness defects of ulp2Δ/Δ cells. Short and long-read genomic sequencing demonstrates that these novel genotypes are selected via a two-step process leading to the formation of novel chromosomal fragments with breakpoints at microhomology regions and DNA repeats.

Single nucleotide polymorphisms and chromosomal copy number variation may impact the Sporothrix brasiliensis antifungal susceptibility and sporotrichosis clinical outcomes

Feline-transmitted sporotrichosis has garnered attention due to the recent high incidence and the lack of efficient control in the epicenter of the epidemic, Rio de Janeiro, Brazil. Sporothrix brasiliensis is the major pathogen involved in feline-to-human sporotrichosis in Brazil and displays more virulent genotypes than the closely related species S. schenckii. Over the last two decades, several reports of antifungal-resistant strains have emerged. Sequencing and comparison analysis of the outbreak strains allowed us to observe that the azole non-wild-type S. brasiliensis strain CFP 1054 had significant chromosomal variations compared to wild-type strains. One of these variants includes a region of 231 Kb containing 75 duplicated genes, which were overrepresented for lipid and isoprenoid metabolism. We also identified an additional strain (CFP 1055) that was resistant to itraconazole and amphotericin B, which had a single nucleotide polymorphism in the tac1 gene. The patients infected with these two strains showed protracted clinical course and sequelae. Even though our sample size is modest, these results suggest the possibility of identifying specific point mutations and large chromosomal duplications potentially associated with antifungal resistance and clinical outcomes of sporotrichosis.

Transposable element expansion and low-level piRNA silencing in grasshoppers may cause genome gigantism

Background

Transposable elements (TEs) have been likened to parasites in the genome that reproduce and move ceaselessly in the host, continuously enlarging the host genome. However, the Piwi-interacting RNA (piRNA) pathway defends animal genomes against the harmful consequences of TE invasion by imposing small-RNA-mediated silencing. Here we compare the TE activity of two grasshopper species with different genome sizes in Acrididae (Locusta migratoria manilensis♀1C = 6.60 pg, Angaracris rhodopa♀1C = 16.36 pg) to ascertain the influence of piRNAs.

Results

We discovered that repetitive sequences accounted for 74.56% of the genome in A. rhodopa, more than 56.83% in L. migratoria, and the large-genome grasshopper contained a higher TEs proportions. The comparative analysis revealed that 41 TEs (copy number > 500) were shared in both species. The two species exhibited distinct “landscapes” of TE divergence. The TEs outbreaks in the small-genome grasshopper occurred at more ancient times, while the large-genome grasshopper maintains active transposition events in the recent past. Evolutionary history studies on TEs suggest that TEs may be subject to different dynamics and resistances in these two species. We found that TE transcript abundance was higher in the large-genome grasshopper and the TE-derived piRNAs abundance was lower than in the small-genome grasshopper. In addition, we found that the piRNA methylase HENMT, which is underexpressed in the large-genome grasshopper, impedes the piRNA silencing to a lower level.

Conclusions

Our study revealed that the abundance of piRNAs is lower in the gigantic genome grasshopper than in the small genome grasshopper. In addition, the key gene HENMT in the piRNA biogenesis pathway (Ping-Pong cycle) in the gigantic genome grasshopper is underexpressed. We hypothesize that low-level piRNA silencing unbalances the original positive correlation between TEs and piRNAs, and triggers TEs to proliferate out of control, which may be one of the reasons for the gigantism of grasshopper genomes.

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

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

Acquisition of cross-azole tolerance and aneuploidy in Candida albicans strains evolved to posaconazole

A number of in vitro studies have examined the acquisition of drug resistance to the triazole fluconazole, a first-line treatment for many Candida infections. Much less is known about posaconazole, a newer triazole. We conducted the first in vitro experimental evolution of replicates from 8 diverse strains of Candida albicans in a high level of the fungistatic drug posaconazole. Approximately half of the 132 evolved replicates survived 50 generations of evolution, biased toward some of the strain backgrounds. We found that although increases in drug resistance were rare, increases in drug tolerance (the slow growth of a subpopulation of cells in a level of drug above the resistance level) were common across strains. We also found that adaptation to posaconazole resulted in widespread cross-tolerance to other azole drugs. Widespread aneuploidy was observed in evolved replicates from some strain backgrounds. Trisomy of at least one of chromosomes 3, 6, and R was identified in 11 of 12 whole-genome sequenced evolved SC5314 replicates. These findings document rampant evolved cross-tolerance among triazoles and highlight that increases in drug tolerance can evolve independently of drug resistance in a diversity of C. albicans strain backgrounds.

Experimental and in-host evolution of triazole resistance in human pathogenic fungi

The leading fungal pathogens causing systemic infections in humans are Candida spp., Aspergillus fumigatus, and Cryptococcus neoformans. The major class of antifungals used to treat such infections are the triazoles, which target the cytochrome P450 lanosterol 14-α-demethylase, encoded by the ERG11 (yeasts)/cyp51A (molds) genes, catalyzing a key step in the ergosterol biosynthetic pathway. Triazole resistance in clinical fungi is a rising concern worldwide, causing increasing mortality in immunocompromised patients. This review describes the use of serial clinical isolates and in-vitro evolution toward understanding the mechanisms of triazole resistance. We outline, compare, and discuss how these approaches have helped identify the evolutionary pathways taken by pathogenic fungi to acquire triazole resistance. While they all share a core mechanism (mutation and overexpression of ERG11/cyp51A and efflux transporters), their timing and mechanism differs: Candida and Cryptococcus spp. exhibit resistance-conferring aneuploidies and copy number variants not seen in A. fumigatus. Candida spp. have a proclivity to develop resistance by undergoing mutations in transcription factors (TAC1, MRR1, PDR5) that increase the expression of efflux transporters. A. fumigatus is especially prone to accumulate resistance mutations in cyp51A early during the evolution of resistance. Recently, examination of serial clinical isolates and experimental lab-evolved triazole-resistant strains using modern omics and gene editing tools has begun to realize the full potential of these approaches. As a result, triazole-resistance mechanisms can now be analyzed at increasingly finer resolutions. This newfound knowledge will be instrumental in formulating new molecular approaches to fight the rapidly emerging epidemic of antifungal resistant fungi.

Collaboration between Antagonistic Cell Type Regulators Governs Natural Variation in the Candida albicans Biofilm and Hyphal Gene Expression Network

Candida albicans is among the most significant human fungal pathogens. However, the vast majority of C. albicans studies have focused on a single clinical isolate and its marked derivatives. We investigated natural variation among clinical C. albicans isolates in gene regulatory control of biofilm formation, a process crucial to virulence. The transcription factor Efg1 is required for biofilm-associated gene expression and biofilm formation. Previously, we found extensive variation in Efg1-responsive gene expression among 5 diverse clinical isolates. However, chromatin immunoprecipitation sequencing analysis showed that Efg1 binding to genomic loci was uniform among the isolates. Functional dissection of strain differences identified three transcription factors, Brg1, Tec1, and Wor1, for which small changes in expression levels reshaped the Efg1 regulatory network. Brg1 and Tec1 are known biofilm activators, and their role in Efg1 network variation may be expected. However, Wor1 is a known repressor of EFG1 expression and an inhibitor of biofilm formation. In contrast, we found that a modest increase in WOR1 RNA levels, reflecting the expression differences between C. albicans strains, could augment biofilm formation and expression of biofilm-related genes. The analysis of natural variation here reveals a novel function for a well-characterized gene and illustrates that strain diversity offers a unique resource for elucidation of network interactions.

The contribution of DNA repair pathways to genome editing and evolution in filamentous pathogens

DNA double-strand breaks require repair or risk corrupting the language of life. To ensure genome integrity and viability, multiple DNA double-strand break repair pathways function in eukaryotes. Two such repair pathways, canonical non-homologous end joining and homologous recombination, have been extensively studied, while other pathways such as microhomology-mediated end joint and single-strand annealing, once thought to serve as back-ups, now appear to play a fundamental role in DNA repair. Here, we review the molecular details and hierarchy of these four DNA repair pathways, and where possible, a comparison for what is known between animal and fungal models. We address the factors contributing to break repair pathway choice, and aim to explore our understanding and knowledge gaps regarding mechanisms and regulation in filamentous pathogens. We additionally discuss how DNA double-strand break repair pathways influence genome engineering results, including unexpected mutation outcomes. Finally, we review the concept of biased genome evolution in filamentous pathogens, and provide a model, termed Biased Variation, that links DNA double-strand break repair pathways with properties of genome evolution. Despite our extensive knowledge for this universal process, there remain many unanswered questions, for which the answers may improve genome engineering and our understanding of genome evolution.

Multiple Stochastic Parameters Influence Genome Dynamics in a Heterozygous Diploid Eukaryotic Model

The heterozygous diploid genome of Candida albicans displays frequent genomic rearrangements, in particular loss-of-heterozygosity (LOH) events, which can be seen on all eight chromosomes and affect both laboratory and clinical strains. LOHs, which are often the consequence of DNA damage repair, can be observed upon stresses reminiscent of the host environment, and result in homozygous regions of various sizes depending on the molecular mechanisms at their origins. Recent studies have shed light on the biological importance of these frequent and ubiquitous LOH events in C. albicans. In diploid Saccharomyces cerevisiae, LOH facilitates the passage of recessive beneficial mutations through Haldane’s sieve, allowing rapid evolutionary adaptation. This also appears to be true in C. albicans, where the full potential of an adaptive mutation is often only observed upon LOH, as illustrated in the case of antifungal resistance and niche adaptation. To understand the genome-wide dynamics of LOH events in C. albicans, we constructed a collection of 15 strains, each one carrying a LOH reporter system on a different chromosome arm. This system involves the insertion of two fluorescent marker genes in a neutral genomic region on both homologs, allowing spontaneous LOH events to be detected by monitoring the loss of one of the fluorescent markers using flow cytometry. Using this collection, we observed significant LOH frequency differences between genomic loci in standard laboratory growth conditions; however, we further demonstrated that comparable heterogeneity was also observed for a given genomic locus between independent strains. Additionally, upon exposure to stress, three outcomes could be observed in C. albicans, where individual strains displayed increases, decreases, or no effect of stress in terms of LOH frequency. Our results argue against a general stress response triggering overall genome instability. Indeed, we showed that the heterogeneity of LOH frequency in C. albicans is present at various levels, inter-strain, intra-strain, and inter-chromosomes, suggesting that LOH events may occur stochastically within a cell, though the genetic background potentially impacts genome stability in terms of LOH throughout the genome in both basal and stress conditions. This heterogeneity in terms of genome stability may serve as an important adaptive strategy for the predominantly clonal human opportunistic pathogen C. albicans, by quickly generating a wide spectrum of genetic variation combinations potentially permitting subsistence in a rapidly evolving environment.

Using genomics to understand the mechanisms of virulence and drug resistance in fungal pathogens

Fungal pathogens pose an increasingly worrying threat to human health, food security and ecosystem diversity. To tackle fungal infections and improve current diagnostic and therapeutic tools it is necessary to understand virulence and antifungal drug resistance mechanisms in diverse species. Recent advances in genomics approaches have provided a suitable framework to understand these phenotypes, which ultimately depend on genetically encoded determinants. In this work, we review how the study of genome sequences has been key to ascertain the bases of virulence and drug resistance traits. We focus on the contribution of comparative genomics, population genomics and directed evolution studies. In addition, we discuss how different types of genomic mutations (small or structural variants) contribute to intraspecific differences in virulence or drug resistance. Finally, we review current challenges in the field and anticipate future directions to solve them. In summary, this work provides a short overview of how genomics can be used to understand virulence and drug resistance in fungal pathogens.

Genome of a giant isopod, Bathynomus jamesi, provides insights into body size evolution and adaptation to deep-sea environment

Background

The deep-sea may be regarded as a hostile living environment, due to low temperature, high hydrostatic pressure, and limited food and light. Isopods, a species-rich group of crustaceans, are widely distributed across different environments including the deep sea and as such are a useful model for studying adaptation, migration, and speciation. Similar to other deep-sea organisms, giant isopods have larger body size than their shallow water relatives and have large stomachs and fat bodies presumably to store organic reserves. In order to shed light on the genetic basis of these large crustaceans adapting to the oligotrophic environment of deep-sea, the high-quality genome of a deep-sea giant isopod Bathynomus jamesi was sequenced and assembled.

Results

B. jamesi has a large genome of 5.89 Gb, representing the largest sequenced crustacean genome to date. Its large genome size is mainly attributable to the remarkable proliferation of transposable elements (84%), which may enable high genome plasticity for adaptive evolution. Unlike its relatives with small body size, B. jamesi has expanded gene families related to pathways of thyroid and insulin hormone signaling that potentially contribute to its large body size. Transcriptomic analysis showed that some expanded gene families related to glycolysis and vesicular transport were specifically expressed in its digestive organs. In addition, comparative genomics and gene expression analyses in six tissues suggested that B. jamesi has inefficient lipid degradation, low basal metabolic rate, and bulk food storage, suggesting giant isopods adopt a more efficient mechanism of nutrient absorption, storage, and utilization to provide sustained energy supply for their large body size.

Conclusions

Taken together, the giant isopod genome may provide a valuable resource for understanding body size evolution and adaptation mechanisms of macrobenthic organisms to deep-sea environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01302-6.

Genome plasticity in Candida albicans: A cutting-edge strategy for evolution, adaptation, and survival

Candida albicans is the most implicated fungal species that grows as a commensal or opportunistic pathogen in the human host. It is associated with many life-threatening infections, especially in immunocompromised persons. The genome of Candida albicans is very flexible and can withstand a wide assortment of variations in a continuously changing environment. Thus, genome plasticity is central to its adaptation and has long been of considerable interest. C. albicans has a diploid heterozygous genome that is highly dynamic and can display variation from small to large scale chromosomal rearrangement and aneuploidy, which have implications in drug resistance, virulence, and pathogenicity. This review presents an up-to-date overview of recent genomic studies involving C. albicans. It discusses the accumulating evidence that shows how mitotic recombination events, ploidy dynamics, aneuploidy, and loss of heterozygosity (LOH) influence evolution, adaptation, and survival in C. albicans. Understanding the factors that affect the genome is crucial for a proper understanding of species and rapid development and adjustment of therapeutic strategies to mitigate their spread.

Emerging roles of SWI/SNF remodelers in fungal pathogens

Fungal pathogens constantly sense and respond to the environment they inhabit, and this interaction is vital for their survival inside hosts and exhibiting pathogenic traits. Since such responses often entail specific patterns of gene expression, regulators of chromatin structure contribute to the fitness and virulence of the pathogens by modulating DNA accessibility to the transcriptional machinery. Recent studies in several human and plant fungal pathogens have uncovered the SWI/SNF group of chromatin remodelers as an important determinant of pathogenic traits and provided insights into their mechanism of function. Here, we review these studies and highlight the differential functions of these remodeling complexes and their subunits in regulating fungal fitness and pathogenicity. As an extension of our previous study, we also show that loss of specific RSC subunits can predispose the human fungal pathogen Candida albicans cells to filamentous growth in a context-dependent manner. Finally, we consider the potential of targeting the fungal SWI/SNF remodeling complexes for antifungal interventions.

Growth of Aspergillus fumigatus in Biofilms in Comparison to Candida albicans

Biofilm formation during infections with the opportunistic pathogen Aspergillus fumigatus can be very problematic in clinical settings, since it provides the fungal cells with a protective environment. Resistance against drug treatments, immune recognition as well as adaptation to the host environment allows fungal survival in the host. The exact molecular mechanisms behind most processes in the formation of biofilms are unclear. In general, the formation of biofilms can be categorized roughly in a few stages; adhesion, conidial germination and development of hyphae, biofilm maturation and cell dispersion. Fungi in biofilms can adapt to the in-host environment. These adaptations can occur on a level of phenotypic plasticity via gene regulation. However, also more substantial genetic changes of the genome can result in increased resistance and adaptation in the host, enhancing the survival chances of fungi in biofilms. Most research has focused on the development of biofilms. However, to tackle developing microbial resistance and adaptation in biofilms, more insight in mechanisms behind genetic adaptations is required to predict which defense mechanisms can be expected. This can be helpful in the development of novel and more targeted antifungal treatments to combat fungal infections.

Chromatin profiling reveals heterogeneity in clinical isolates of the human pathogen Aspergillus fumigatus

Invasive Pulmonary Aspergillosis, which is caused by the filamentous fungus Aspergillus fumigatus, is a life-threatening infection for immunosuppressed patients. Chromatin structure regulation is important for genome stability maintenance and has the potential to drive genome rearrangements and affect virulence and pathogenesis of pathogens. Here, we performed the first A. fumigatus global chromatin profiling of two histone modifications, H3K4me3 and H3K9me3, focusing on the two most investigated A. fumigatus clinical isolates, Af293 and CEA17. In eukaryotes, H3K4me3 is associated with active transcription, while H3K9me3 often marks silent genes, DNA repeats, and transposons. We found that H3K4me3 deposition is similar between the two isolates, while H3K9me3 is more variable and does not always represent transcriptional silencing. Our work uncovered striking differences in the number, locations, and expression of transposable elements between Af293 and CEA17, and the differences are correlated with H3K9me3 modifications and higher genomic variations among strains of Af293 background. Moreover, we further showed that the Af293 strains from different laboratories actually differ in their genome contents and found a frequently lost region in chromosome VIII. For one such Af293 variant, we identified the chromosomal changes and demonstrated their impacts on its secondary metabolites production, growth and virulence. Overall, our findings not only emphasize the influence of genome heterogeneity on A. fumigatus fitness, but also caution about unnoticed chromosomal variations among common laboratory strains.

Parasexuality of Candida Species

While most fungi have the ability to reproduce sexually, multiple independent lineages have lost meiosis and developed parasexual cycles in its place. Emergence of parasexual cycles is particularly prominent in medically relevant fungi from the CUG paraphyletic group of Candida species. Since the discovery of parasex in C. albicans roughly two decades ago, it has served as the model for Candida species. Importantly, parasex in C. albicans retains hallmarks of meiosis including genetic recombination and chromosome segregation, making it a potential driver of genetic diversity. Furthermore, key meiotic genes play similar roles in C. albicans parasex and highlights parallels between these processes. Yet, the evolutionary role of parasex in Candida adaptation and the extent of resulting genotypic and phenotypic diversity remain as key knowledge gaps in this facultative reproductive program. Here, we present our current understanding of parasex, the mechanisms governing its regulation, and its relevance to Candida biology.

Dynamic genome plasticity during unisexual reproduction in the human fungal pathogen Cryptococcus deneoformans

Genome copy number variation occurs during each mitotic and meiotic cycle and it is crucial for organisms to maintain their natural ploidy. Defects in ploidy transitions can lead to chromosome instability, which is a hallmark of cancer. Ploidy in the haploid human fungal pathogen Cryptococcus neoformans is exquisitely orchestrated and ranges from haploid to polyploid during sexual development and under various environmental and host conditions. However, the mechanisms controlling these ploidy transitions are largely unknown. During C. deneoformans (formerly C. neoformans var. neoformans, serotype D) unisexual reproduction, ploidy increases prior to the onset of meiosis, can be independent from cell-cell fusion and nuclear fusion, and likely occurs through an endoreplication pathway. To elucidate the molecular mechanisms underlying this ploidy transition, we identified twenty cell cycle-regulating genes encoding cyclins, cyclin-dependent kinases (CDK), and CDK regulators. We characterized four cyclin genes and two CDK regulator genes that were differentially expressed during unisexual reproduction and contributed to diploidization. To detect ploidy transition events, we generated a ploidy reporter, called NURAT, which can detect copy number increases via double selection for nourseothricin-resistant, uracil-prototrophic cells. Utilizing this ploidy reporter, we showed that ploidy transition from haploid to diploid can be detected during the early phases of unisexual reproduction. Interestingly, selection for the NURAT reporter revealed several instances of segmental aneuploidy of multiple chromosomes, which conferred azole resistance in some isolates. These findings provide further evidence of ploidy plasticity in fungi with significant biological and public health implications.

Ploidy is an intrinsic fundamental feature of all eukaryotic organisms, and ploidy variation and maintenance are critical to the organism survival and evolution. Fungi exhibit exquisite plasticity in ploidy variation in adaptation to various environmental stresses. For example, the haploid opportunistic human fungal pathogen C. deneoformans can generate diploid blastospores during unisexual reproduction and also forms polyploid titan cells during host infection; however, the mechanisms underlying these ploidy transitions are largely unknown. In this study, we elucidated the genetic regulatory circuitry governing ploidy duplication during C. deneoformans unisexual reproduction through the identification and characterization of cell cycle regulators that are differentially expressed during unisexual reproduction. We showed that four cyclin and two cyclin-dependent kinase regulator genes function in concert to orchestrate ploidy transition during unisexual reproduction. To trace and track ploidy transition events, we also generated a ploidy reporter and revealed the formation of segmental aneuploidy in addition to diploidization, illustrating the diverse mechanisms of genome plasticity in C. deneoformans.

The Genome of the CTG(Ser1) Yeast Scheffersomyces stipitis Is Plastic

Microorganisms need to adapt to environmental changes, and genome plasticity can lead to rapid adaptation to hostile environments by increasing genetic diversity. Here, we investigate genome plasticity in the CTG(Ser1) yeast Scheffersomyces stipitis, an organism with an enormous potential for second-generation biofuel production. We demonstrate that S. stipitis has an intrinsically plastic genome and that different S. stipitis isolates have genomes with distinct chromosome organizations. Real-time evolution experiments show that S. stipitis genome plasticity is common and rapid since extensive genomic changes with fitness benefits are detected following in vitro evolution experiments. Hybrid MinION Nanopore and Illumina genome sequencing identify retrotransposons as major drivers of genome diversity. Indeed, the number and position of retrotransposons are different in different S. stipitis isolates, and retrotransposon-rich regions of the genome are sites of chromosome rearrangements. Our findings provide important insights into the adaptation strategies of the CTG(Ser1) yeast clade and have critical implications in the development of second-generation biofuels. These data highlight that genome plasticity is an essential factor for developing sustainable S. stipitis platforms for second-generation biofuels production. IMPORTANCE Genomes contain genes encoding the information needed to build the organism and allow it to grow and develop. Genomes are described as stable structures where genes have specific positions within a chromosome. Changes in gene dosage and position are viewed as harmful. However, it is becoming increasingly clear that genome plasticity can benefit microbial organisms that need to adapt rapidly to environmental changes. Mechanisms of genome plasticity are still poorly understood. This study focuses on Scheffersomyces stipitis, a yeast that holds great potential for second-generation biofuel production generated from forestry and agriculture waste. We demonstrate that S. stipitis chromosomes are easily reshuffled and that chromosome reshuffling is linked to adaptation to hostile environments. Genome sequencing demonstrates that mobile genetic elements, called transposons, mediate S. stipitis genome reshuffling. These data highlight that understanding genome plasticity is important for developing sustainable S. stipitis platforms for second-generation biofuels production.

Integrative multi-omics profiling reveals cAMP-independent mechanisms regulating hyphal morphogenesis in Candida albicans

Microbial pathogens grow in a wide range of different morphologies that provide distinct advantages for virulence. In the fungal pathogen Candida albicans, adenylyl cyclase (Cyr1) is thought to be a master regulator of the switch to invasive hyphal morphogenesis and biofilm formation. However, faster growing cyr1Δ/Δ pseudorevertant (PR) mutants were identified that form hyphae in the absence of cAMP. Isolation of additional PR mutants revealed that their improved growth was due to loss of one copy of BCY1, the negative regulatory subunit of protein kinase A (PKA) from the left arm of chromosome 2. Furthermore, hyphal morphogenesis was improved in some of PR mutants by multigenic haploinsufficiency resulting from loss of large regions of the left arm of chromosome 2, including global transcriptional regulators. Interestingly, hyphal-associated genes were also induced in a manner that was independent of cAMP. This indicates that basal protein kinase A activity is an important prerequisite to induce hyphae, but activation of adenylyl cyclase is not needed. Instead, phosphoproteomic analysis indicated that the Cdc28 cyclin-dependent kinase and the casein kinase 1 family member Yck2 play key roles in promoting polarized growth. In addition, integrating transcriptomic and proteomic data reveals hyphal stimuli induce increased production of key transcription factors that contribute to polarized morphogenesis.

Genome and transcriptome of a pathogenic yeast, Candida nivariensis

We present a highly contiguous genome and transcriptome of the pathogenic yeast, Candida nivariensis. We sequenced both the DNA and RNA of this species using both the Oxford Nanopore Technologies and Illumina platforms. We assembled the genome into an 11.8 Mb draft composed of 16 contigs with an N50 of 886 Kb, including a circular mitochondrial sequence of 28 Kb. Using direct RNA nanopore sequencing and Illumina cDNA sequencing, we constructed an annotation of our new assembly, supplemented by lifting over genes from Saccharomyces cerevisiae and Candida glabrata.

Resistance to Antifungal Drugs

Pathogenic fungi have several mechanisms of resistance to antifungal drugs, driven by the genetic plasticity and versatility of their homeostatic responses to stressful environmental cues. We critically review the molecular mechanisms of resistance and cellular adaptations of pathogenic fungi in response to antifungals and discuss the factors contributing to such resistance. We offer suggestions for the translational and clinical research agenda of this rapidly evolving and medically important field. A better understanding of antifungal resistance should assist in developing better detection tools and inform optimal strategies for preventing and treating refractory mycoses in the future.

Simple sequence repeats drive genome plasticity and promote adaptive evolution in penaeid shrimp

Simple sequence repeats (SSRs) are rare (approximately 1%) in most genomes and are generally considered to have no function. However, penaeid shrimp genomes have a high proportion of SSRs (>23%), raising the question of whether these SSRs play important functional and evolutionary roles in these SSR-rich species. Here, we show that SSRs drive genome plasticity and adaptive evolution in two penaeid shrimp species, Fenneropenaeus chinensis and Litopenaeus vannamei. Assembly and comparison of genomes of these two shrimp species at the chromosome-level revealed that transposable elements serve as carriers for SSR expansion, which is still occurring. The remarkable genome plasticity identified herein might have been shaped by significant SSR expansions. SSRs were also found to regulate gene expression by multi-omics analyses, and be responsible for driving adaptive evolution, such as the variable osmoregulatory capacities of these shrimp under low-salinity stress. These data provide strong evidence that SSRs are an important driver of the adaptive evolution in penaeid shrimp.

Comparative genomics of white and opaque cell states supports an epigenetic mechanism of phenotypic switching in Candida albicans

Several Candida species can undergo a heritable and reversible transition from a 'white' state to a mating proficient 'opaque' state. This ability relies on highly interconnected transcriptional networks that control cell-type-specific gene expression programs over multiple generations. Candida albicans, the most prominent pathogenic Candida species, provides a well-studied paradigm for the white-opaque transition. In this species, a network of at least eight transcriptional regulators controls the balance between white and opaque states that have distinct morphologies, transcriptional profiles, and physiological properties. Given the reversible nature and the high frequency of white-opaque transitions, it is widely assumed that this switch is governed by epigenetic mechanisms that occur independently of any changes in DNA sequence. However, a direct genomic comparison between white and opaque cells has yet to be performed. Here, we present a whole-genome comparative analysis of C. albicans white and opaque cells. This analysis revealed rare genetic changes between cell states, none of which are linked to white-opaque switching. This result is consistent with epigenetic mechanisms controlling cell state differentiation in C. albicans and provides direct evidence against a role for genetic variation in mediating the switch.

Genome compartmentalization predates species divergence in the plant pathogen genus Zymoseptoria

Background

Antagonistic co-evolution can drive rapid adaptation in pathogens and shape genome architecture. Comparative genome analyses of several fungal pathogens revealed highly variable genomes, for many species characterized by specific repeat-rich genome compartments with exceptionally high sequence variability. Dynamic genome structure may enable fast adaptation to host genetics. The wheat pathogen Zymoseptoria tritici with its highly variable genome, has emerged as a model organism to study genome evolution of plant pathogens. Here, we compared genomes of Z. tritici isolates and of sister species infecting wild grasses to address the evolution of genome composition and structure.

Results

Using long-read technology, we sequenced and assembled genomes of Z. ardabiliae, Z. brevis, Z. pseudotritici and Z. passerinii, together with two isolates of Z. tritici. We report a high extent of genome collinearity among Zymoseptoria species and high conservation of genomic, transcriptomic and epigenomic signatures of compartmentalization. We identify high gene content variability both within and between species. In addition, such variability is mainly limited to the accessory chromosomes and accessory compartments. Despite strong host specificity and non-overlapping host-range between species, predicted effectors are mainly shared among Zymoseptoria species, yet exhibiting a high level of presence-absence polymorphism within Z. tritici. Using in planta transcriptomic data from Z. tritici, we suggest different roles for the shared orthologs and for the accessory genes during infection of their hosts.

Conclusion

Despite previous reports of high genomic plasticity in Z. tritici, we describe here a high level of conservation in genomic, epigenomic and transcriptomic composition and structure across the genus Zymoseptoria. The compartmentalized genome allows the maintenance of a functional core genome co-occurring with a highly variable accessory genome.

Expandable and reversible copy number amplification drives rapid adaptation to antifungal drugs

Previously, we identified long repeat sequences that are frequently associated with genome rearrangements, including copy number variation (CNV), in many diverse isolates of the human fungal pathogen Candida albicans (Todd et al., 2019). Here, we describe the rapid acquisition of novel, high copy number CNVs during adaptation to azole antifungal drugs. Single-cell karyotype analysis indicates that these CNVs appear to arise via a dicentric chromosome intermediate and breakage-fusion-bridge cycles that are repaired using multiple distinct long inverted repeat sequences. Subsequent removal of the antifungal drug can lead to a dramatic loss of the CNV and reversion to the progenitor genotype and drug susceptibility phenotype. These findings support a novel mechanism for the rapid acquisition of antifungal drug resistance and provide genomic evidence for the heterogeneity frequently observed in clinical settings.

Candida and Candidiasis-Opportunism Versus Pathogenicity: A Review of the Virulence Traits

One of the most important questions in microbiology nowadays, is how apparently harmless, commensal yeasts like Candida spp. can cause a rising number of infections. The occurrence of the disease requires firstly the attachment to the host cells, followed by the invasion of the tissue. The adaptability translates into a rapid ability to respond to stress factors, to take up nutrients or to multiply under different conditions. By forming complex intracellular networks such as biofilms, Candida spp. become not only more refractive to antifungal therapies but also more prone to cause disease. The inter-microbial interactions can enhance the virulence of a strain. In vivo, the fungal cells face a multitude of challenges and, as a result, they develop complex strategies serving one ultimate goal: survival. This review presents the virulence factors of the most important Candida spp., contributing to a better understanding of the onset of candidiasis and raising awareness of the highly complex interspecies interactions that can change the outcome of the disease.

Identification and Characterization of a Neutral Locus for Knock-in Purposes in C. parapsilosis

Invasive fungal infections caused by Candida species affect approximately 700,000 people worldwide resulting in 300,000 deaths annually. Besides Candida albicans, other members of the genus have gained relevance in the last two decades, including C. parapsilosis whose incidence is particularly high amongst low birth weight neonates. To investigate the virulence properties of this pathogen several techniques have been developed for generating knock-out mutants, however, no target locus for knock-in approaches have been published so far. Here we report CpNEUT5L (N5L), an intergenic locus in C. parapsilosis, and introduce an integrative Gateway^TM and a classical ligation based replacement plasmid to target it with. As a proof of principle, we fluorescently tagged laboratory and prototroph strains and established that this locus is also suitable for reintegration purposes. We concluded that GFP-expressing constructs integrated into this region provide strong, homogenous fluorescent signals while alteration of this locus affects neither the growth of the mutants in liquid or on solid media, even in the presence of different stressors, nor their basic virulence properties. Hence, our findings demonstrate that N5L is a highly effective neutral locus for knock-in approaches in C. parapsilosis.

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.

Drug resistance and tolerance in fungi

Systemic fungal infections pose a serious clinical problem. Treatment options are limited, and antifungal drug resistance is increasing. In addition, a substantial proportion of patients do not respond to therapy despite being infected with fungi that are susceptible to the drug. The discordance between overall treatment outcome and low levels of clinical resistance may be attributable to antifungal drug tolerance. In this Review, we define and distinguish resistance and tolerance and discuss the current understanding of the molecular, genetic and physiological mechanisms that contribute to those phenomena. Distinguishing tolerance from resistance might provide important insights into the reasons for treatment failure in some settings.