Using in vivo transcriptomics and RNA enrichment to identify genes involved in virulence of Candida glabrata

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

Transcriptomic meta-analysis to identify potential antifungal targets in Candida albicans

BACKGROUND

Candida albicans is a fungal pathogen causing human infections. Here we investigated differential gene expression patterns and functional enrichment in C. albicans strains grown under different conditions.

METHODS

A systematic GEO database search identified 239 "Candida albicans" datasets, of which 14 were selected after rigorous criteria application. Retrieval of raw sequencing data from the ENA database was accompanied by essential metadata extraction from dataset descriptions and original articles. Pre-processing via the tailored nf-core pipeline for C. albicans involved alignment, gene/transcript quantification, and diverse quality control measures. Quality assessment via PCA and DESeq2 identified significant genes (FDR < = 0.05, log2-fold change > = 1 or <= -1), while topGO conducted GO term enrichment analysis. Exclusions were made based on data quality and strain relevance, resulting in the selection of seven datasets from the SC5314 strain background for in-depth investigation.

RESULTS

The meta-analysis of seven selected studies unveiled a substantial number of genes exhibiting significant up-regulation (24,689) and down-regulation (18,074). These differentially expressed genes were further categorized into 2,497 significantly up-regulated and 2,573 significantly down-regulated Gene Ontology (GO) IDs. GO term enrichment analysis clustered these terms into distinct groups, providing insights into the functional implications. Three target gene lists were compiled based on previous studies, focusing on central metabolism, ion homeostasis, and pathogenicity. Frequency analysis revealed genes with higher occurrence within the identified GO clusters, suggesting their potential as antifungal targets. Notably, the genes TPS2, TPS1, RIM21, PRA1, SAP4, and SAP6 exhibited higher frequencies within the clusters. Through frequency analysis within the GO clusters, several key genes emerged as potential targets for antifungal therapies. These include RSP5, GLC7, SOD2, SOD5, SOD1, SOD6, SOD4, SOD3, and RIM101 which exhibited higher occurrence within the identified clusters.

CONCLUSION

This comprehensive study significantly advances our understanding of the dynamic nature of gene expression in C. albicans. The identification of genes with enhanced potential as antifungal drug targets underpins their value for future interventions. The highlighted genes, including TPS2, TPS1, RIM21, PRA1, SAP4, SAP6, RSP5, GLC7, SOD2, SOD5, SOD1, SOD6, SOD4, SOD3, and RIM101, hold promise for the development of targeted antifungal therapies.

Bibliometric analysis and thematic review of Candida pathogenesis: Fundamental omics to applications as potential antifungal drugs and vaccines

Invasive candidiasis caused by the pathogenic Candida yeast species has resulted in elevating global mortality. The pathogenicity of Candida spp. is not only originated from its primary invasive yeast-to-hyphal transition; virulence factors (transcription factors, adhesins, invasins, and enzymes), biofilm, antifungal drug resistance, stress tolerance, and metabolic adaptation have also contributed to a greater clinical burden. However, the current research theme in fungal pathogenicity could hardly be delineated with the increasing research output. Therefore, our study analysed the research trends in Candida pathogenesis over the past 37 years via a bibliometric approach against the Scopus and Web of Science databases. Based on the 3993 unique documents retrieved, significant international collaborations among researchers were observed, especially between Germany (Bernhard Hube) and the UK (Julian Naglik), whose focuses are on Candida proteinases, adhesins, and candidalysin. The prominent researchers (Neils Gow, Alistair Brown, and Frank Odds) at the University of Exeter and the University of Aberdeen (second top performing affiliation) UK contribute significantly to the mechanisms of Candida adaptation, tolerance, and stress response. However, the science mapping of co-citation analysis performed herein could not identify a hub representative of subsequent work since the clusters were semi-redundant. The co-word analysis that was otherwise adopted, revealed three research clusters; the cluster-based thematic analyses indicated the severeness of Candida biofilm and antifungal resistance as well as the elevating trend on molecular mechanism elucidation for drug screening and repurposing. Importantly, the in vivo pathogen adaptation and interactions with hosts are crucial for potential vaccine development.

A Putative NADPH Oxidase Gene in Unicellular Pathogenic Candida glabrata Is Required for Fungal ROS Production and Oxidative Stress Response

Most previous studies on fungal NADPH oxidases (Nox) focused on multicellular fungi and highlighted the important roles of Nox-derived reactive oxygen species (ROS) in cellular differentiation and signaling communication. However, there are few reports about Nox in unicellular fungi. A novel NOX ortholog, CAGL0K05863g (named CgNOX1), in Candida glabrata was investigated in this study. Deletion of CgNOX1 led to a decrease in both intracellular and extracellular ROS production. In addition, the Cgnox1∆ mutant exhibited hypersensitivity to hydrogen peroxide and menadione. Also, the wild-type strain showed higher levels of both CgNOX1 mRNA expression and ROS production under oxidative stress. Moreover, the absence of CgNOX1 resulted in impaired ferric reductase activity. Although there was no effect on in vitro biofilm formation, the CgNOX1 mutant did not produce hepatic apoptosis, which might be mediated by fungal Nox-derived ROS during co-incubation. Together, these results indicated that the novel NOX gene plays important roles in unicellular pathogenic C. glabrata and its interaction with host cells.

Pseudomonas aeruginosa responds to altered membrane phospholipid composition by adjusting the production of two-component systems, proteases and iron uptake proteins

Membrane protein and phospholipid (PL) composition changes in response to environmental cues and during infections. To achieve these, bacteria use adaptation mechanisms involving covalent modification and remodelling of the acyl chain length of PLs. However, little is known about bacterial pathways regulated by PLs. Here, we investigated proteomic changes in the biofilm of P. aeruginosa phospholipase mutant (∆plaF) with altered membrane PL composition. The results revealed profound alterations in the abundance of many biofilm-related two-component systems (TCSs), including accumulation of PprAB, a key regulator of the transition to biofilm. Furthermore, a unique phosphorylation pattern of transcriptional regulators, transporters and metabolic enzymes, as well as differential production of several proteases, in ∆plaF, indicate that PlaF-mediated virulence adaptation involves complex transcriptional and posttranscriptional response. Moreover, proteomics and biochemical assays revealed the depletion of pyoverdine-mediated iron uptake pathway proteins in ∆plaF, while proteins from alternative iron-uptake systems were accumulated. These suggest that PlaF may function as a switch between different iron-acquisition pathways. The observation that PL-acyl chain modifying and PL synthesis enzymes were overproduced in ∆plaF reveals the interconnection of degradation, synthesis and modification of PLs for proper membrane homeostasis. Although the precise mechanism by which PlaF simultaneously affects multiple pathways remains to be elucidated, we suggest that alteration of PL composition in ∆plaF plays a role for the global adaptive response in P. aeruginosa mediated by TCSs and proteases. Our study revealed the global regulation of virulence and biofilm by PlaF and suggests that targeting this enzyme may have therapeutic potential.

The Candida glabrata Parent Strain Trap: How Phenotypic Diversity Affects Metabolic Fitness and Host Interactions

Reference strains improve reproducibility by standardizing observations and methodology, which has ultimately led to important insights into fungal pathogenesis. However, recent investigations have highlighted significant genotypic and phenotypic heterogeneity across isolates that influence genetic circuitry and virulence within a species. Candida glabrata is the second leading cause of candidiasis, a life-threatening infection, and undergoes extensive karyotype and phenotypic changes in response to stress. Much of the work conducted on this pathogen has focused on two sequenced strains, CBS138 (ATCC 2001) and BG2. Few studies have compared these strains in detail, but key differences include mating type and altered patterns of expression of EPA adhesins. In fact, most C. glabrata isolates and BG2 are MATa, while CBS138 is MATα. However, it is not known if other phenotypic differences between these strains play a role in our understanding of C. glabrata pathogenesis. Thus, we set out to characterize metabolic, cell wall, and host-interaction attributes for CBS138 and BG2. We found that BG2 utilized a broader range of nitrogen sources and had reduced cell wall size and carbohydrate exposure than CBS138, which we hypothesized results in differences in innate immune interactions and virulence. We observed that, although both strains were phagocytosed to a similar extent, BG2 replicated to higher numbers in macrophages and was more virulent during Galleria mellonella infection than CBS138 in a dose-dependent manner. Interestingly, deletion of SNF3, a major nutrient sensor, did not affect virulence in G. mellonella for BG2, but significantly enhanced larval killing in the CBS138 background compared to the parent strain. Understanding these fundamental differences in metabolism and host interactions will allow more robust conclusions to be drawn in future studies of C. glabrata pathogenesis. IMPORTANCE Reference strains provide essential insights into the mechanisms underlying virulence in fungal pathogens. However, recent studies in Candida albicans and other species have revealed significant genotypic and phenotypic diversity within clinical isolates that are challenging paradigms regarding key virulence factors and their regulation. Candida glabrata is the second leading cause of candidiasis, and many studies use BG2 or CBS138 for their investigations. Therefore, we aimed to characterize important virulence-related phenotypes for both strains that might alter conclusions about C. glabrata pathogenesis. Our study provides context for metabolic and cell wall changes and how these may influence host interaction phenotypes. Understanding these differences is necessary to support robust conclusions about how virulence factors may function in these and other very different strain backgrounds.

Multiplexed target enrichment of coding and non-coding transcriptomes enables studying Candida spp. infections from human derived samples

The study of transcriptomic interactions between host and pathogens in in vivo conditions is challenged by the low relative amounts of the pathogen RNA. Yeast opportunistic pathogens of the genus Candida can cause life-threatening systemic infections in immunocompromised patients, and are of growing medical concern. Four phylogenetically diverse species account for over 90% of Candida infections, and their specific interactions with various human tissues are still poorly understood. To enable in vivo transcriptomic analysis in these species, we designed and validated pan-Candida target capture probes to enrich protein-coding and non-coding transcriptomes. The probe-based enrichment approach outperformed enrichment based on differential lysis of host cells, and showed similar enrichment performance as an existing capture design, yet achieving better fidelity of expression levels, enabling species multiplexing and capturing of lncRNAs. In addition, we show that our probe-based enrichment strategy allows robust genotype-based identification of the infecting strain present in the sample.

A novel Candida glabrata doxycycline-inducible system for in vitro/in vivo use

Candida glabrata is an important pathogen causing superficial to invasive disease in human. Conditional expression systems are helpful in addressing the function of genes and especially when they can be applied to in vivo studies. Tetracycline-dependent regulation systems have been used in diverse fungi to turn-on (Tet-on) or turn-off (Tet-off) gene expression either in vitro but also in vivo in animal models. Up to now, only a Tet-off expression has been constructed for gene expression in C. glabrata. Here, we report a Tet-on gene expression system which can be used in vitro and in vivo in any C. glabrata genetic background. This system was used in a mice model of systemic infection to demonstrate that the general amino acid permease Gap1 is important for C. glabrata virulence.