Novel Drivers of Virulence in Clostridioides difficile Identified via Context-Specific Metabolic Network Analysis

Dual RNA-seq study of the dynamics of coding and non-coding RNA expression during Clostridioides difficile infection in a mouse model

Clostridioides difficile is the leading cause of healthcare-associated diarrhea in industrialized countries. Many questions remain to be answered about the mechanisms governing its interaction with the host during infection. Non-coding RNAs (ncRNAs) contribute to shape virulence in many pathogens and modulate host responses; however, their role in C. difficile infection (CDI) has not been explored. To better understand the dynamics of ncRNA expression contributing to C. difficile infectious cycle and host response, we used a dual RNA-seq approach in a conventional murine model. From the pathogen side, this transcriptomic analysis revealed the upregulation of virulence factors, metabolism, and sporulation genes, as well as the identification of 61 ncRNAs differentially expressed during infection that correlated with the analysis of available raw RNA-seq data sets from two independent studies. From these data, we identified 118 potential new transcripts in C. difficile, including 106 new ncRNA genes. From the host side, we observed the induction of several pro-inflammatory pathways, and among the 185 differentially expressed ncRNAs, the overexpression of microRNAs (miRNAs) previously associated to inflammatory responses or unknown long ncRNAs and miRNAs. A particular host gene expression profile could be associated to the symptomatic infection. In accordance, the metatranscriptomic analysis revealed specific microbiota changes accompanying CDI and specific species associated with symptomatic infection in mice. This first adaptation of in vivo dual RNA-seq to C. difficile contributes to unravelling the regulatory networks involved in C. difficile infectious cycle and host response and provides valuable resources for further studies of RNA-based mechanisms during CDI.IMPORTANCEClostridioides difficile is a major cause of nosocomial infections associated with antibiotic therapy classified as an urgent antibiotic resistance threat. This pathogen interacts with host and gut microbial communities during infection, but the mechanisms of these interactions remain largely to be uncovered. Noncoding RNAs contribute to bacterial virulence and host responses, but their expression has not been explored during C. difficile infection. We took advantage of the conventional mouse model of C. difficile infection to look simultaneously to the dynamics of gene expression in pathogen, its host, and gut microbiota composition, providing valuable resources for future studies. We identified a number of ncRNAs that could mediate the adaptation of C. difficile inside the host and the crosstalk with the host immune response. Promising inflammation markers and potential therapeutic targets emerged from this work open new directions for RNA-based and microbiota-modulatory strategies to improve the efficiency of C. difficile infection treatments.

Elucidating human gut microbiota interactions that robustly inhibit diverse Clostridioides difficile strains across different nutrient landscapes

The human gut pathogen Clostridioides difficile displays substantial inter-strain genetic variability and confronts a changeable nutrient landscape in the gut. We examined how human gut microbiota inter-species interactions influence the growth and toxin production of various C. difficile strains across different nutrient environments. Negative interactions influencing C. difficile growth are prevalent in an environment containing a single highly accessible resource and sparse in an environment containing C. difficile-preferred carbohydrates. C. difficile toxin production displays significant community-context dependent variation and does not trend with growth-mediated inter-species interactions. C. difficile strains exhibit differences in interactions with Clostridium scindens and the ability to compete for proline. Further, C. difficile shows substantial differences in transcriptional profiles in co-culture with C. scindens or Clostridium hiranonis. C. difficile exhibits massive alterations in metabolism and other cellular processes in co-culture with C. hiranonis, reflecting their similar metabolic niches. C. hiranonis uniquely inhibits the growth and toxin production of diverse C. difficile strains across different nutrient environments and robustly ameliorates disease severity in mice. In sum, understanding the impact of C. difficile strain variability and nutrient environments on inter-species interactions could help improve the effectiveness of anti-C. difficile strategies.

Emerging tools and best practices for studying gut microbial community metabolism

The human gut microbiome vastly extends the set of metabolic reactions catalysed by our own cells, with far-reaching consequences for host health and disease. However, our knowledge of gut microbial metabolism relies on a handful of model organisms, limiting our ability to interpret and predict the metabolism of complex microbial communities. In this Perspective, we discuss emerging tools for analysing and modelling the metabolism of gut microorganisms and for linking microorganisms, pathways and metabolites at the ecosystem level, highlighting promising best practices for researchers. Continued progress in this area will also require infrastructure development to facilitate cross-disciplinary synthesis of scientific findings. Collectively, these efforts can enable a broader and deeper understanding of the workings of the gut ecosystem and open new possibilities for microbiome manipulation and therapy.

Clostridioides difficile-mucus interactions encompass shifts in gene expression, metabolism, and biofilm formation

In a healthy colon, the stratified mucus layer serves as a crucial innate immune barrier to protect the epithelium from microbes. Mucins are complex glycoproteins that serve as a nutrient source for resident microflora and can be exploited by pathogens. We aimed to understand how the intestinal pathogen, Clostridioides difficile, independently uses or manipulates mucus to its benefit, without contributions from members of the microbiota. Using a 2-D primary human intestinal epithelial cell model to generate physiologic mucus, we assessed C. difficile-mucus interactions through growth assays, RNA-Seq, biophysical characterization of mucus, and contextualized metabolic modeling. We found that host-derived mucus promotes C. difficile growth both in vitro and in an infection model. RNA-Seq revealed significant upregulation of genes related to central metabolism in response to mucus, including genes involved in sugar uptake, the Wood-Ljungdahl pathway, and the glycine cleavage system. In addition, we identified differential expression of genes related to sensing and transcriptional control. Analysis of mutants with deletions in highly upregulated genes reflected the complexity of C. difficile-mucus interactions, with potential interplay between sensing and growth. Mucus also stimulated biofilm formation in vitro, which may in turn alter the viscoelastic properties of mucus. Context-specific metabolic modeling confirmed differential metabolism and the predicted importance of enzymes related to serine and glycine catabolism with mucus. Subsequent growth experiments supported these findings, indicating mucus is an important source of serine. Our results better define responses of C. difficile to human gastrointestinal mucus and highlight flexibility in metabolism that may influence pathogenesis.

IMPORTANCE

Clostridioides difficile results in upward of 250,000 infections and 12,000 deaths annually in the United States. Community-acquired infections continue to rise, and recurrent disease is common, emphasizing a vital need to understand C. difficile pathogenesis. C. difficile undoubtedly interacts with colonic mucus, but the extent to which the pathogen can independently respond to and take advantage of this niche has not been explored extensively. Moreover, the metabolic complexity of C. difficile remains poorly understood but likely impacts its capacity to grow and persist in the host. Here, we demonstrate that C. difficile uses native colonic mucus for growth, indicating C. difficile possesses mechanisms to exploit the mucosal niche. Furthermore, mucus induces metabolic shifts and biofilm formation in C. difficile, which has potential ramifications for intestinal colonization. Overall, our work is crucial to better understand the dynamics of C. difficile-mucus interactions in the context of the human gut.

Elucidating human gut microbiota interactions that robustly inhibit diverse Clostridioides difficile strains across different nutrient landscapes

The human gut pathogen Clostridioides difficile displays extreme genetic variability and confronts a changeable nutrient landscape in the gut. We mapped gut microbiota inter-species interactions impacting the growth and toxin production of diverse C. difficile strains in different nutrient environments. Although negative interactions impacting C. difficile are prevalent in environments promoting resource competition, they are sparse in an environment containing C. difficile-preferred carbohydrates. C. difficile strains display differences in interactions with Clostridium scindens and the ability to compete for proline. C. difficile toxin production displays substantial community-context dependent variation and does not trend with growth-mediated inter-species interactions. C. difficile shows substantial differences in transcriptional profiles in the presence of the closely related species C. hiranonis or C. scindens. In co-culture with C. hiranonis, C. difficile exhibits massive alterations in metabolism and other cellular processes, consistent with their high metabolic overlap. Further, Clostridium hiranonis inhibits the growth and toxin production of diverse C. difficile strains across different nutrient environments and ameliorates the disease severity of a C. difficile challenge in a murine model. In sum, strain-level variability and nutrient environments are major variables shaping gut microbiota interactions with C. difficile.

Metabolic model-based ecological modeling for probiotic design

The microbial community composition in the human gut has a profound effect on human health. This observation has lead to extensive use of microbiome therapies, including over-the-counter 'probiotic' treatments intended to alter the composition of the microbiome. Despite so much promise and commercial interest, the factors that contribute to the success or failure of microbiome-targeted treatments remain unclear. We investigate the biotic interactions that lead to successful engraftment of a novel bacterial strain introduced to the microbiome as in probiotic treatments. We use pairwise genome-scale metabolic modeling with a generalized resource allocation constraint to build a network of interactions between taxa that appear in an experimental engraftment study. We create induced sub-graphs using the taxa present in individual samples and assess the likelihood of invader engraftment based on network structure. To do so, we use a generalized Lotka-Volterra model, which we show has strong ability to predict if a particular invader or probiotic will successfully engraft into an individual's microbiome. Furthermore, we show that the mechanistic nature of the model is useful for revealing which microbe-microbe interactions potentially drive engraftment.

Clostridioides difficile-mucus interactions encompass shifts in gene expression, metabolism, and biofilm formation

In a healthy colon, the stratified mucus layer serves as a crucial innate immune barrier to protect the epithelium from microbes. Mucins are complex glycoproteins that serve as a nutrient source for resident microflora and can be exploited by pathogens. We aimed to understand how the intestinal pathogen, Clostridioides diffiicile, independently uses or manipulates mucus to its benefit, without contributions from members of the microbiota. Using a 2-D primary human intestinal epithelial cell model to generate physiologic mucus, we assessed C. difficile-mucus interactions through growth assays, RNA-Seq, biophysical characterization of mucus, and contextualized metabolic modeling. We found that host-derived mucus promotes C. difficile growth both in vitro and in an infection model. RNA-Seq revealed significant upregulation of genes related to central metabolism in response to mucus, including genes involved in sugar uptake, the Wood-Ljungdahl pathway, and the glycine cleavage system. In addition, we identified differential expression of genes related to sensing and transcriptional control. Analysis of mutants with deletions in highly upregulated genes reflected the complexity of C. difficile-mucus interactions, with potential interplay between sensing and growth. Mucus also stimulated biofilm formation in vitro, which may in turn alter viscoelastic properties of mucus. Context-specific metabolic modeling confirmed differential metabolism and predicted importance of enzymes related to serine and glycine catabolism with mucus. Subsequent growth experiments supported these findings, indicating mucus is an important source of serine. Our results better define responses of C. difficile to human gastrointestinal mucus and highlight a flexibility in metabolism that may influence pathogenesis.

Applications of genome-scale metabolic models to investigate microbial metabolic adaptations in response to genetic or environmental perturbations

Environmental perturbations are encountered by microorganisms regularly and will require metabolic adaptations to ensure an organism can survive in the newly presenting conditions. In order to study the mechanisms of metabolic adaptation in such conditions, various experimental and computational approaches have been used. Genome-scale metabolic models (GEMs) are one of the most powerful approaches to study metabolism, providing a platform to study the systems level adaptations of an organism to different environments which could otherwise be infeasible experimentally. In this review, we are describing the application of GEMs in understanding how microbes reprogram their metabolic system as a result of environmental variation. In particular, we provide the details of metabolic model reconstruction approaches, various algorithms and tools for model simulation, consequences of genetic perturbations, integration of '-omics' datasets for creating context-specific models and their application in studying metabolic adaptation due to the change in environmental conditions.

Inferring microbial interactions with their environment from genomic and metagenomic data

Microbial communities assemble through a complex set of interactions between microbes and their environment, and the resulting metabolic impact on the host ecosystem can be profound. Microbial activity is known to impact human health, plant growth, water quality, and soil carbon storage which has lead to the development of many approaches and products meant to manipulate the microbiome. In order to understand, predict, and improve microbial community engineering, genome-scale modeling techniques have been developed to translate genomic data into inferred microbial dynamics. However, these techniques rely heavily on simulation to draw conclusions which may vary with unknown parameters or initial conditions, rather than more robust qualitative analysis. To better understand microbial community dynamics using genome-scale modeling, we provide a tool to investigate the network of interactions between microbes and environmental metabolites over time. Using our previously developed algorithm for simulating microbial communities from genome-scale metabolic models (GSMs), we infer the set of microbe-metabolite interactions within a microbial community in a particular environment. Because these interactions depend on the available environmental metabolites, we refer to the networks that we infer as metabolically contextualized, and so name our tool MetConSIN: Metabolically Contextualized Species Interaction Networks.

Transcriptome-guided metabolic network analysis reveals rearrangements of carbon flux distribution in Neisseria gonorrhoeae during neutrophil co-culture

The ability of bacterial pathogens to metabolically adapt to the environmental conditions of their hosts is critical to both colonization and invasive disease. Infection with Neisseria gonorrhoeae (the gonococcus, Gc) is characterized by the influx of neutrophils [polymorphonuclear leukocytes (PMNs)], which fail to clear the bacteria and make antimicrobial products that can exacerbate tissue damage. The inability of the human host to clear Gc infection is particularly concerning in light of the emergence of strains that are resistant to all clinically recommended antibiotics. Bacterial metabolism represents a promising target for the development of new therapeutics against Gc. Here, we generated a curated genome-scale metabolic network reconstruction (GENRE) of Gc strain FA1090. This GENRE links genetic information to metabolic phenotypes and predicts Gc biomass synthesis and energy consumption. We validated this model with published data and in new results reported here. Contextualization of this model using the transcriptional profile of Gc exposed to PMNs revealed substantial rearrangements of Gc central metabolism and induction of Gc nutrient acquisition strategies for alternate carbon source use. These features enhanced the growth of Gc in the presence of neutrophils. From these results, we conclude that the metabolic interplay between Gc and PMNs helps define infection outcomes. The use of transcriptional profiling and metabolic modeling to reveal new mechanisms by which Gc persists in the presence of PMNs uncovers unique aspects of metabolism in this fastidious bacterium, which could be targeted to block infection and thereby reduce the burden of gonorrhea in the human population. IMPORTANCE The World Health Organization designated Gc as a high-priority pathogen for research and development of new antimicrobials. Bacterial metabolism is a promising target for new antimicrobials, as metabolic enzymes are widely conserved among bacterial strains and are critical for nutrient acquisition and survival within the human host. Here we used genome-scale metabolic modeling to characterize the core metabolic pathways of this fastidious bacterium and to uncover the pathways used by Gc during culture with primary human immune cells. These analyses revealed that Gc relies on different metabolic pathways during co-culture with human neutrophils than in rich media. Conditionally essential genes emerging from these analyses were validated experimentally. These results show that metabolic adaptation in the context of innate immunity is important to Gc pathogenesis. Identifying the metabolic pathways used by Gc during infection can highlight new therapeutic targets for drug-resistant gonorrhea.

Systems-ecology designed bacterial consortium protects from severe Clostridioides difficile infection

Fecal Microbiota Transplant (FMT) is an emerging therapy that has had remarkable success in treatment and prevention of recurrent Clostridioides difficile infection (rCDI). FMT has recently been associated with adverse outcomes such as inadvertent transfer of antimicrobial resistance, necessitating development of more targeted bacteriotherapies. To address this challenge, we developed a novel systems biology pipeline to identify candidate probiotic strains that would be predicted to interrupt C. difficile pathogenesis. Utilizing metagenomic characterization of human FMT donor samples, we identified those metabolic pathways most associated with successful FMTs and reconstructed the metabolism of encoding species to simulate interactions with C. difficile . This analysis resulted in predictions of high levels of cross-feeding for amino acids in species most associated with FMT success. Guided by these in silico models, we assembled consortia of bacteria with increased amino acid cross-feeding which were then validated in vitro . We subsequently tested the consortia in a murine model of CDI, demonstrating total protection from severe CDI through decreased toxin levels, recovered gut microbiota, and increased intestinal eosinophils. These results support the novel framework that amino acid cross-feeding is likely a critical mechanism in the initial resolution of CDI by FMT. Importantly, we conclude that our predictive platform based on predicted and testable metabolic interactions between the microbiota and C. difficile led to a rationally designed biotherapeutic framework that may be extended to other enteric infections.

Understanding Antimicrobial Resistance Using Genome-Scale Metabolic Modeling

The urgent necessity to fight antimicrobial resistance is universally recognized. In the search of new targets and strategies to face this global challenge, a promising approach resides in the study of the cellular response to antimicrobial exposure and on the impact of global cellular reprogramming on antimicrobial drugs' efficacy. The metabolic state of microbial cells has been shown to undergo several antimicrobial-induced modifications and, at the same time, to be a good predictor of the outcome of an antimicrobial treatment. Metabolism is a promising reservoir of potential drug targets/adjuvants that has not been fully exploited to date. One of the main problems in unraveling the metabolic response of cells to the environment resides in the complexity of such metabolic networks. To solve this problem, modeling approaches have been developed, and they are progressively gaining in popularity due to the huge availability of genomic information and the ease at which a genome sequence can be converted into models to run basic phenotype predictions. Here, we review the use of computational modeling to study the relationship between microbial metabolism and antimicrobials and the recent advances in the application of genome-scale metabolic modeling to the study of microbial responses to antimicrobial exposure.

Network analysis of toxin production in Clostridioides difficile identifies key metabolic dependencies

Clostridioides difficile pathogenesis is mediated through its two toxin proteins, TcdA and TcdB, which induce intestinal epithelial cell death and inflammation. It is possible to alter C. difficile toxin production by changing various metabolite concentrations within the extracellular environment. However, it is unknown which intracellular metabolic pathways are involved and how they regulate toxin production. To investigate the response of intracellular metabolic pathways to diverse nutritional environments and toxin production states, we use previously published genome-scale metabolic models of C. difficile strains CD630 and CDR20291 (iCdG709 and iCdR703). We integrated publicly available transcriptomic data with the models using the RIPTiDe algorithm to create 16 unique contextualized C. difficile models representing a range of nutritional environments and toxin states. We used Random Forest with flux sampling and shadow pricing analyses to identify metabolic patterns correlated with toxin states and environment. Specifically, we found that arginine and ornithine uptake is particularly active in low toxin states. Additionally, uptake of arginine and ornithine is highly dependent on intracellular fatty acid and large polymer metabolite pools. We also applied the metabolic transformation algorithm (MTA) to identify model perturbations that shift metabolism from a high toxin state to a low toxin state. This analysis expands our understanding of toxin production in C. difficile and identifies metabolic dependencies that could be leveraged to mitigate disease severity.

Enterococci enhance Clostridioides difficile pathogenesis

Enteric pathogens are exposed to a dynamic polymicrobial environment in the gastrointestinal tract¹. This microbial community has been shown to be important during infection, but there are few examples illustrating how microbial interactions can influence the virulence of invading pathogens². Here we show that expansion of a group of antibiotic-resistant, opportunistic pathogens in the gut-the enterococci-enhances the fitness and pathogenesis of Clostridioides difficile. Through a parallel process of nutrient restriction and cross-feeding, enterococci shape the metabolic environment in the gut and reprogramme C. difficile metabolism. Enterococci provide fermentable amino acids, including leucine and ornithine, which increase C. difficile fitness in the antibiotic-perturbed gut. Parallel depletion of arginine by enterococci through arginine catabolism provides a metabolic cue for C. difficile that facilitates increased virulence. We find evidence of microbial interaction between these two pathogenic organisms in multiple mouse models of infection and patients infected with C. difficile. These findings provide mechanistic insights into the role of pathogenic microbiota in the susceptibility to and the severity of C. difficile infection.

d-Proline Reductase Underlies Proline-Dependent Growth of Clostridioides difficile

Clostridioides difficile is a nosocomial pathogen that colonizes the gut and causes diarrhea, colitis, and severe inflammation. Recently, C. difficile has been shown to use toxin-mediated inflammation to promote host collagen degradation, which releases several amino acids into the environment. Amino acids act as electron donors and acceptors in Stickland metabolism, an anaerobic process involving redox reactions between pairs of amino acids. Proline, glycine, and hydroxyproline are the three main constituents of collagen and are assumed to act as electron acceptors, but their exact effects on the growth and physiology of C. difficile are still unclear. Using three standard culture media (supplemented brain heart infusion [BHIS], tryptone-yeast [TY], and C. difficile minimal medium [CDMM]) supplemented with proline, glycine, or hydroxyproline, we grew C. difficile strains R20291, JIR8094, and a panel of mutants unable to express the Stickland selenoenzymes d-proline reductase and glycine reductase. In the wild-type strains, growth yields in rich media (BHIS and TY) were higher with proline and hydroxyproline but not glycine; moreover, proline-stimulated growth yields required the activity of d-proline reductase, whereas hydroxyproline-stimulated growth yields were independent of its activity. While assumed to be a proline auxotroph, C. difficile could surprisingly grow in a defined medium (CDMM) without proline but only if d-proline reductase was absent. We believe the mere presence of this enzyme ultimately determines the organism's strict dependence on proline and likely defines the bioenergetic priorities for thriving in the host. Finally, we demonstrated that addition of proline and hydroxyproline to the culture medium could reduce toxin production but not in cells lacking selenoproteins. IMPORTANCE Stickland metabolism is a core facet of C. difficile physiology that likely plays a major role in host colonization. Here, we carefully delineate the effects of each amino acid on the growth of C. difficile with respect to the selenoenzymes d-proline reductase and glycine reductase. Moreover, we report that d-proline reductase forces C. difficile to strictly depend on proline for growth. Finally, we provide evidence that proline and hydroxyproline suppress toxin production and that selenoproteins are involved in this mechanism. Our findings highlight the significance of selenium-dependent Stickland reactions and may provide insight on what occurs during host infection, especially as it relates to the decision to colonize based on proline as a nutrient.

Computational approaches to understanding Clostridioides difficile metabolism and virulence

The progress of infection by Clostridioides difficile is strongly influenced by metabolic cues it encounters as it colonizes the gastrointestinal tract. Both colonization and regulation of virulence have a multi-factorial interaction between host, microbiome, and gene expression cascades. While these connections with metabolism have been understood for some time, many mechanisms of control have remained difficult to directly assay due to high metabolic variability among C. difficile isolates and difficult genetic systems. Computational systems offer a means to interrogate structure of complex or noisy datasets and generate useful, tractable hypotheses to be tested in the laboratory. Recently, in silico techniques have provided powerful insights into metabolic elements of C. difficile infection ranging from virulence regulation to interactions with the gut microbiota. In this review, we introduce and provide context to the methods of computational modeling that have been applied to C. difficile metabolism and virulence thus far. The techniques discussed here have laid the foundation for future multi-scale efforts aimed at understanding the complex interplay of metabolic activity between pathogen, host, and surrounding microbial community in the regulation of C. difficile pathogenesis.