PerSVade: personalized structural variant detection in any species of interest

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.

Corrected and republished from: "Extracellular Vesicle Formation in Cryptococcus deuterogattii Impacts Fungal Virulence"

Small molecules are components of fungal extracellular vesicles (EVs), but their biological roles are only superficially known. NOP16 is a eukaryotic gene that is required for the activity of benzimidazoles against Cryptococcus deuterogattii. In this study, during the phenotypic characterization of C. deuterogattii mutants expected to lack NOP16 expression, we observed a reduced EV production. Whole-genome sequencing, RNA-Seq, and cellular proteomics revealed that, contrary to our initial findings, these mutants expressed Nop16 but exhibited altered expression of 14 genes potentially involved in sugar transport. Based on this observation, we designated these mutant strains as Past1 and Past2, representing potentially altered sugar transport. Analysis of the small molecule composition of EVs produced by wild-type cells and the Past1 and Past2 mutant strains revealed not only a reduced number of EVs but also an altered small molecule composition. In a Galleria mellonella model of infection, the Past1 and Past2 mutant strains were hypovirulent. The hypovirulent phenotype was reverted when EVs produced by wild-type cells, but not mutant EVs, were co-injected with the mutant cells in G. mellonella. These results connect EV biogenesis, cargo, and cryptococcal virulence.

Calling Structural Variants with Confidence from Short-Read Data in Wild Bird Populations

Comprehensive characterization of structural variation in natural populations has only become feasible in the last decade. To investigate the population genomic nature of structural variation, reproducible and high-confidence structural variation callsets are first required. We created a population-scale reference of the genome-wide landscape of structural variation across 33 Nordic house sparrows (Passer domesticus). To produce a consensus callset across all samples using short-read data, we compare heuristic-based quality filtering and visual curation (Samplot/PlotCritic and Samplot-ML) approaches. We demonstrate that curation of structural variants is important for reducing putative false positives and that the time invested in this step outweighs the potential costs of analyzing short-read-discovered structural variation data sets that include many potential false positives. We find that even a lenient manual curation strategy (e.g. applied by a single curator) can reduce the proportion of putative false positives by up to 80%, thus enriching the proportion of high-confidence variants. Crucially, in applying a lenient manual curation strategy with a single curator, nearly all (>99%) variants rejected as putative false positives were also classified as such by a more stringent curation strategy using three additional curators. Furthermore, variants rejected by manual curation failed to reflect the expected population structure from SNPs, whereas variants passing curation did. Combining heuristic-based quality filtering with rapid manual curation of structural variants in short-read data can therefore become a time- and cost-effective first step for functional and population genomic studies requiring high-confidence structural variation callsets.

Genomics of the expanding pine pathogen Lecanosticta acicola reveals patterns of ongoing genetic admixture

Lecanosticta acicola is the causal agent for brown spot needle blight that affects pine trees across the northern hemisphere. Based on marker genes and microsatellite data, two distinct lineages have been identified that were introduced into Europe on two separate occasions. Despite their overall distinct geographic distribution, they have been found to coexist in regions of northern Spain and France. Here, we present the first genome-wide study of Lecanosticta acicola, including assembly of the reference genome and a population genomics analysis of 70 natural isolates from northern Spain. We show that most of the isolates belong to the southern lineage but show signs of introgression with northern lineage isolates, indicating mating between the two lineages. We also identify phenotypic differences between the two lineages based on the activity profiles of 20 enzymes, with introgressed strains being more phenotypically similar to members of the southern lineage. In conclusion, we show undergoing genetic admixture between the two main lineages of L. acicola in a region of recent expansion.

IMPORTANCE

Lecanosticta acicola is a fungal pathogen causing severe defoliation, growth reduction, and even death in more than 70 conifer species. Despite the increasing incidence of this species, little is known about its population dynamics. Two divergent lineages have been described that have now been found together in regions of France and Spain, but it is unknown how these mixed populations evolve. Here we present the first reference genome for this important plant pathogenic fungi and use it to study the population genomics of 70 isolates from an affected forest in the north of Spain. We find signs of introgression between the two main lineages, indicating that active mating is occurring in this region which could propitiate the appearance of novel traits in this species. We also study the phenotypic differences across this population based on enzymatic activities on 20 compounds.

Recent gene selection and drug resistance underscore clinical adaptation across Candida species

Understanding how microbial pathogens adapt to treatments, humans and clinical environments is key to infer mechanisms of virulence, transmission and drug resistance. This may help improve therapies and diagnostics for infections with a poor prognosis, such as those caused by fungal pathogens, including Candida. Here we analysed genomic variants across approximately 2,000 isolates from six Candida species (C. glabrata, C. auris, C. albicans, C. tropicalis, C. parapsilosis and C. orthopsilosis) and identified genes under recent selection, suggesting a highly complex clinical adaptation. These involve species-specific and convergently affected adaptive mechanisms, such as adhesion. Using convergence-based genome-wide association studies we identified known drivers of drug resistance alongside potentially novel players. Finally, our analyses reveal an important role of structural variants and suggest an unexpected involvement of (para)sexual recombination in the spread of resistance. Our results provide insights on how opportunistic pathogens adapt to human-related environments and unearth candidate genes that deserve future attention.

Overlooked Candida glabrata petites are echinocandin tolerant, induce host inflammatory responses, and display poor in vivo fitness

IMPORTANCE

Candida glabrata is a major fungal pathogen, which is able to lose mitochondria and form small and slow-growing colonies, called "petite." This attenuated growth rate has created controversies and questioned the clinical importance of petiteness. Herein, we have employed multiple omics technologies and in vivo mouse models to critically assess the clinical importance of petite phenotype. Our WGS identifies multiple genes potentially underpinning petite phenotype. Interestingly, petite C. glabrata cells engulfed by macrophages are dormant and, therefore, are not killed by the frontline antifungal drugs. Interestingly, macrophages infected with petite cells mount distinct transcriptomic responses. Consistent with our ex vivo observations, mitochondrial-proficient parental strains outcompete petites during systemic and gut colonization. Retrospective examination of C. glabrata isolates identified petite prevalence a rare entity, which can significantly vary from country to country. Collectively, our study overcomes the existing controversies and provides novel insights regarding the clinical relevance of petite C. glabrata isolates.

Proteogenomics in Cancer: Then and Now

For years, the paths of sequencing technologies and mass spectrometry have occurred in isolation, with each developing its own unique culture and expertise. These two technologies are crucial for inspecting complementary aspects of the molecular phenotype across the central dogma. Integrative multiomics strives to bridge the analysis gap among different fields to complete more comprehensive mechanisms of life events and diseases. Proteogenomics is one integrated multiomics field. Here in this review, we mainly summarize and discuss three aspects: workflow of proteogenomics, proteogenomics applications in cancer research, and the SWOT (Strengths, Weaknesses, Opportunities, Threats) analysis of proteogenomics in cancer research. In conclusion, proteogenomics has a promising future as it clarifies the functional consequences of many unannotated genomic abnormalities or noncanonical variants and identifies driver genes and novel therapeutic targets across cancers, which would substantially accelerate the development of precision oncology.

Overlooked Candida glabrata petites are echinocandin tolerant, induce host inflammatory responses, and display poor in vivo fitness

Small colony variants (SCVs) are relatively common among some bacterial species and are associated with poor prognosis and recalcitrant infections. Similarly, Candida glabrata - a major intracellular fungal pathogen - produces small and slow-growing respiratory-deficient colonies, termed "petite." Despite reports of clinical petite C . glabrata strains, our understanding of petite behavior in the host remains obscure. Moreover, controversies exist regarding in-host petite fitness and its clinical relevance. Herein, we employed whole-genome sequencing (WGS), dual-RNAseq, and extensive ex vivo and in vivo studies to fill this knowledge gap. WGS identified multiple petite-specific mutations in nuclear and mitochondrially-encoded genes. Consistent with dual-RNAseq data, petite C . glabrata cells did not replicate inside host macrophages and were outcompeted by their non-petite parents in macrophages and in gut colonization and systemic infection mouse models. The intracellular petites showed hallmarks of drug tolerance and were relatively insensitive to the fungicidal activity of echinocandin drugs. Petite-infected macrophages exhibited a pro-inflammatory and type I IFN-skewed transcriptional program. Interrogation of international C . glabrata blood isolates ( n =1000) showed that petite prevalence varies by country, albeit at an overall low prevalence (0-3.5%). Collectively, our study sheds new light on the genetic basis, drug susceptibility, clinical prevalence, and host-pathogen responses of a clinically overlooked phenotype in a major fungal pathogen.

Importance

Candida glabrata is a major fungal pathogen, which is able to lose mitochondria and form small and slow-growing colonies, called "petite". This attenuated growth rate has created controversies and questioned the clinical importance of petiteness. Herein, we have employed multiple omicstechnologies and in vivo mouse models to critically assess the clinical importance of petite phenotype. Our WGS identifies multiple genes potentially underpinning petite phenotype. Interestingly, petite C. glabrata cells engulfed by macrophages are dormant and therefore are not killed by the frontline antifungal drugs. Interestingly, macrophages infected with petite cells mount distinct transcriptomic responses. Consistent with our ex-vivo observations, mitochondrial-proficient parental strains outcompete petites during systemic and gut colonization. Retrospective examination of C. glabrata isolates identified petite prevalence a rare entity, can significantly vary from country to country. Collectively, our study overcomes the existing controversies and provides novel insights regarding the clinical relevance of petite C. glabrata isolates.

Whole-genome sequencing confirms a persistent candidaemia clonal outbreak due to multidrug-resistant Candida parapsilosis

OBJECTIVES

Although perceived as a rare clinical entity, recent studies have noted the emergence of MDR C. parapsilosis (MDR-Cp) isolates from single patients (resistant to both azole and echinocandins). We previously reported a case series of MDR-Cp isolates carrying a novel FKS1R658G mutation. Herein, we identified an echinocandin-naive patient infected with MDR-Cp a few months after the previously described isolates. WGS and CRISPR-Cas9 editing were used to explore the origin of the new MDR-Cp isolates, and to determine if the novel mutation confers echinocandin resistance.

METHODS

WGS was applied to assess the clonality of these isolates and CRISPR-Cas9 editing and a Galleria mellonella model were used to examine whether FKS1R658G confers echinocandin resistance.

RESULTS

Fluconazole treatment failed, and the patient was successfully treated with liposomal amphotericin B (LAMB). WGS proved that all historical and novel MDR-Cp strains were clonal and distant from the fluconazole-resistant outbreak cluster in the same hospital. CRISPR-Cas9 editing and G. mellonella virulence assays confirmed that FKS1R658G confers echinocandin resistance in vitro and in vivo. Interestingly, the FKS1R658G mutant showed a very modest fitness cost compared with the parental WT strain, consistent with the persistence of the MDR-Cp cluster in our hospital.

CONCLUSIONS

Our study showcases the emergence of MDR-Cp isolates as a novel threat in clinical settings, which undermines the efficacy of the two most widely used antifungal drugs against candidiasis, leaving only LAMB as a last resort. Additionally, surveillance studies and WGS are warranted to effectively establish infection control and antifungal stewardship strategies.

The challenges of the genome-based identification of antifungal resistance in the clinical routine

The increasing number of chronic and life-threatening infections caused by antimicrobial resistant fungal isolates is of critical concern. Low DNA sequencing cost may facilitate the identification of the genomic profile leading to resistance, the resistome, to rationally optimize the design of antifungal therapies. However, compared to bacteria, initiatives for resistome detection in eukaryotic pathogens are underdeveloped. Firstly, reported mutations in antifungal targets leading to reduced susceptibility must be extensively collected from the literature to generate comprehensive databases. This information should be complemented with specific laboratory screenings to detect the highest number possible of relevant genetic changes in primary targets and associations between resistance and other genomic markers. Strikingly, some drug resistant strains experience high-level genetic changes such as ploidy variation as much as duplications and reorganizations of specific chromosomes. Such variations involve allelic dominance, gene dosage increments and target expression regime effects that should be explicitly parameterized in antifungal resistome prediction algorithms. Clinical data indicate that predictors need to consider the precise pathogen species and drug levels of detail, instead of just genus and drug class. The concomitant needs for mutation accuracy and assembly quality assurance suggest hybrid sequencing approaches involving third-generation methods will be utilized. Moreover, fatal fast infections, like fungemia and meningitis, will further require both sequencing and analysis facilities are available in-house. Altogether, the complex nature of antifungal resistance demands extensive sequencing, data acquisition and processing, bioinformatic analysis pipelines, and standard protocols to be accomplished prior to genome-based protocols are applied in the clinical setting.