Novel Approaches for Systems Biology of Metabolism-Oriented Pathogen-Human Interactions: A Mini-Review

Pathogenic microorganisms exploit host metabolism for sustained survival by rewiring its metabolic interactions. Therefore, several metabolic changes are induced in both pathogen and host cells in the course of infection. A systems-based approach to elucidate those changes includes the integrative use of genome-scale metabolic networks and molecular omics data, with the overall goal of better characterizing infection mechanisms for novel treatment strategies. This review focuses on novel aspects of metabolism-oriented systems-based investigation of pathogen-human interactions. The reviewed approaches are the generation of dual-omics data for the characterization of metabolic signatures of pathogen-host interactions, the reconstruction of pathogen-host integrated genome-scale metabolic networks, which has a high potential to be applied to pathogen-gut microbiota interactions, and the structure-based analysis of enzymes playing role in those interactions. The integrative use of those approaches will pave the way for the identification of novel biomarkers and drug targets for the prediction and prevention of infectious diseases.

Introduction

One major focus of systems biology in medical research is the integrated use of high-throughput omics data and molecular interaction networks, with the overall goal of elucidating molecular mechanisms of diseases to identify potential biomarkers and personalized treatment strategies. A wealth of research has been conducted to apply this approach on infectious diseases by focusing on cell metabolism. This included collecting transcriptomic, proteomic, or metabolomic data from pathogens or infected hosts on one side, and developing interaction models focusing on the enzymes and reactions of crucial metabolic pathways on the other side (Bumann, 2009).

Rather than investigating host and pathogen metabolisms separately, novel systems biology based approaches consider the pathogen-host interaction (PHI) system as a whole to elucidate the mechanisms of infection. Here, omics data is collected from both organisms in a single experiment, termed dual-omics data, focusing on all the molecules within the PHI system (Westermann et al., 2017). Partly possible thanks to the developments in Next Generation Sequencing technologies, this enables the generation of valuable omics data in shorter times.

Reconstruction of microorganism-specific metabolic networks based on genome sequences is a powerful approach when combined with the computational tools that enable the prediction of phenotypic indicators such as rates of metabolic reactions, activities of metabolic pathways, effect of gene deletion on growth. When applied to pathogenic organisms, it is possible to identify potential biomarkers and drug targets with this approach (Kim et al., 2012; Dunphy and Papin, 2018). There are several studies that report the reconstruction of genome-scale metabolic networks (GMNs) for pathogens, which are reviewed elsewhere (Cesur et al., 2018). The generation of dual-omics data makes it possible to extend the genome-scale metabolic network modeling of pathogens to PHI systems, by simultaneously considering pathogen and host metabolisms.

One crucial contribution of PHI-based integrated metabolic approach is to identify host non-homologous metabolic pathways of pathogenic resistant microorganisms to propose new drug targets. Resistant pathogens use their unique metabolic pathways to overcome the harsh conditions produced by the administered antibiotics. Identification of those pathways and related enzymes using protein homology and structure is important to predict novel drug targets for infectious diseases.

This mini-review focuses on the recently emerging PHI-based approach to metabolism in human infection, covering studies using dual-omics approach, pathogen-host integrated GMNs and structure and homology based analysis of enzymes. The reviewed approaches are illustrated in Figure 1.

FIGURE 1

Figure 1. Pathogen-host metabolic interactions. (A) When a host cell is intracellularly infected by a pathogen, several metabolite exchanges are possible between the pathogen and the host. (B) In response to infection several metabolic pathways will get activated/inactivated in both pathogen and host in addition to the activation of some metabolic exchanges between the two cells. Dual-omics studies provide experimental tools to take a snapshot of the PHI systems, enabling prediction of active/inactive metabolic reactions in a specific infection time/condition. Constraint-based analysis approaches, on the other hand, provide computational tools to map omic data on integrated pathogen-host metabolic networks to identify active parts. (C) Such metabolic-network based frameworks can also be developed to predict metabolic interactions between pathogens and gut microbes in gut. (D) The integrated metabolic networks can be processed for the prediction of selective drug targets and corresponding drugs.

Metabolism-Oriented Dual Omics of Pathogen-Host Interactions

It is possible to perform simultaneous omics analysis of the pathogen and host together, generating dual omics data focusing on both organisms (Durmuş et al., 2015; Westermann et al., 2017). First examples of metabolism-oriented analysis of dual transcriptomic data have recently been available, which aim to gain insights specifically on the metabolism of the pathogen and the host during infection. These examples are reviewed below.

Zimmermann et al. (2017) collected metabolomics and dual RNA-Seq data of the macrophages infected by Mycobacterium tuberculosis. Integrating the omics data with the combined genome-scale metabolic network model of macrophage and M. tuberculosis provided insights on the metabolic robustness and resistance of the bacteria to metabolic interventions. Calculation of metabolic fluxes using the combined GMN constrained by the dual RNA-Seq data generated predictions of co-utilization of 33 different carbon sources. The results enlightened the substrates directly used by the pathogen as biomass precursors and the ones further metabolized for energy or building blocks. On the other hand, pathogen/host joint metabolomics and dual transcriptomics data were investigated together to shed light on the metabolic changes during Toxoplasma gondii infection of human (Olson et al., 2018). Paired analysis of joint metabolome and dual transcriptome data uncovered the manipulation of the host metabolome by Toxoplasma and identified sedoheptulose biphosphatase driven ribose synthesis from glucose as a novel metabolic capability of the parasite. The identified metabolic enzyme was proposed as a potential drug target since it was not present in human. In another study, Tucey et al. (2018) investigated the crosstalk between glucose metabolism of immune cells and that of pathogenic fungus Candida albicans, focusing on differentially expressed pathogen and host genes during infection. From the analysis of dual RNA-seq data, the glucose competition by C. albicans was found to be responsible for the death of lots of infected macrophages. The results provided evidence for a key role of host glucose homeostasis in vivo during animal infection, and it was proposed that a glucose-rich diet improved host outcomes in Candida infection. In a recent study, ulcer-associated pathogen Haemophilus ducreyi was investigated in detail by collecting dual transcriptomic and host metabolomic data from infected human tissue (Griesenauer et al., 2019). The results suggest the consumption of ascorbic acid and adaptation of anaerobic metabolism as survival mechanisms by the pathogen in glucose-poor abscess environment.

There are a number of other dual-transcriptomic analyses of PHI systems in literature without a specific focus on metabolism, but they also briefly report associated metabolic alterations. These studies are reviewed in Table 1.

TABLE 1

Table 1. Key metabolic findings from dual transcriptome analysis of pathogen-mammalian host systems, given in chronological order.

Genome-Scale Metabolic Networks of Pathogen-Host Interactions

Since pathogens tightly interact with the host cells during infection, they cause alterations in the level of available metabolites in the host cell, leading to rewiring of host metabolic pathways. Therefore, focusing solely on pathogen metabolism in infection research has limitations. The genome-scale reconstruction of pathogen-host integrated metabolic networks is therefore crucial for a holistic analysis of metabolic changes in infection. Guidelines specific to the reconstruction of such integrated metabolic networks were published (Jamshidi and Raghunathan, 2015; Raghunathan and Jamshidi, 2018). Contrary to its utmost importance, studies with such integrated metabolic networks are currently limited in number in literature. The details of the available studies are reviewed in this section.

The first integrated genome-scale metabolic network of pathogen-host interactions was reported for the Plasmodium falciparum infection in erythrocytes, where infection-specific gene expression data of the parasite was incorporated into flux prediction algorithm (Huthmacher et al., 2010). The authors reconstructed a metabolic network of human erythrocytes with 349 reactions and a network of 998 reactions controlled by 579 genes for the malaria pathogen. In the integrated network simulations, P. falciparum was forced to consume some host metabolites to better represent infection characteristics. Compared to the use of only the metabolic network of the pathogen, the simulation of the integrated metabolic network better predicted the metabolites exchanged between pathogen and the host when integrated with transcriptomic data. The authors later applied the same approach for the hepatocyte infection of P. falciparum, the first site of infection in human body for this pathogen (Bazzani et al., 2012). The GMN of human hepatocyte, Hepatocyte1, includes 2,539 reactions controlled by 704 genes. The network was integrated with an updated network of P. falciparum with 1,394 reactions and 579 genes. By leaving all pathogen-host metabolite exchange rates unconstrained, they performed gene deletion and reduced fitness simulations. The integrated analysis enabled the prediction of 24 enzymes as selective drug targets, which are essential in the pathogen but non-essential in hepatocytes.

Another early example of integrated metabolic network approach is for the pathogen M. tuberculosis and its infection of alveolar macrophages (Bordbar et al., 2010). To this aim, the authors used a genome scale reconstruction of the pathogen with 661 genes and 1,027 reactions. The host metabolic network was derived from the generic human metabolic reconstruction Recon1 with the help of macrophage-specific gene expression data, leading to 3,394 reactions and 1,410 genes. This corresponds to the removal of 86 genes from the Recon1 model. The exchange of metabolites between the two networks was allowed based on literature information. The authors also updated the composition of biomass reaction in the pathogen network to better represent in vivo infection conditions. They showed that gene deletion simulations with the integrated network leads to improved gene essentiality predictions when compared to the use of the pathogen metabolic network alone. The pathogen-host metabolic network was also integrated with the transcriptome data from infected macrophages. This shed light on the effect of three different types of M. tuberculosis infection on the metabolisms of the pathogen and the host. Later, the same integrated metabolic network was integrated with dual RNA-seq data, as mentioned in the previous section (Zimmermann et al., 2017). The results predicted a multi-nutrient strategy for the pathogen during early infection. Additional gene essentiality simulations indicated that the number of essential genes decreases during infection. Recently, others investigated the same pathogenic infection of macrophages by using improved metabolic networks (Rienksma et al., 2019). The M. tuberculosis metabolic network used included 1,192 reactions and 915 genes. As the host GMN, a more recent reconstruction of generic human metabolism, Recon 2.2, was used, with 7,785 reactions and 1,675 associated genes. Dual RNA-seq data from macrophage-like THP-1 cells infected with a close relative, M. bovis BCG, was integrated with the reconstructed pathogen-host metabolic network. First, condition-specific biomass reactions were created by the help of transcriptome data for both cell types. Later, the response of host and pathogen metabolic networks to 11 metabolically active anti-tuberculosis drugs was individually simulated by gradually increasing the flux through the affected reactions. This enabled the identification of reroutings in the pathogen metabolism.

Genome-Scale Metabolic Networks for the Interaction of Gut Microbiota and Infectious Pathogens

Microbial ecology has witnessed tremendous progress over the last decade empowered by new sequencing technologies. These innovations in DNA sequencing resulted in a paradigm shift concerning our understanding of pathogenic processes and gave birth to a new scientific concept, namely the “pathobiome” (Vayssier-Taussat et al., 2014). In the pathobiome concept, the microbial community containing a pathogenic agent exerts a major influence on the persistence, transmission, and evolution of the pathogen. For example, reduction of indigenous Lachnospiraceae in the human gut following antibiotic treatment creates a niche for Clostridium difficile infections (Fleming-Davies et al., 2017; Jenior et al., 2017). Sexually transmitted infections such as Chlamydia trachomatis and human papilloma virus have been associated with reduction of Lactobacillus species that results in high-diversity vaginal microbiota (Sewankambo et al., 1997). Interestingly, metabolic products or functions of the microbiota can be critical for the survival of invading species. The microbiota regulates the bile acid pool whereas specific bile acids are required for the maximum germination of C. difficile. Short chain fatty acids levels, which can limit the growth and colonization of pathogens by disrupting the intracellular pH homeostasis, as in the case of Salmonella enterica, are also influenced by the composition of the microbiota (McHan and Shotts, 1993; Horswill et al., 2001; Jacobson et al., 2018). It is therefore mandatory to identify those species that contribute, either as partners or antagonists, to pathogen's survival, extinction or dispersal by applying concepts and approaches of systems biology. Developing systems biology frameworks would allow us to go beyond the identification of mere statistical co-occurrence patterns and instead develop computational models that describe the spatio-temporal dynamics of the pathogen and the functional interactions with its biotic environment.

GMNs have only recently been used for gaining a mechanistic insight into the interactions between communities and their host, as well as, between ecosystem members. For the former, one of the first attempts was based on the “supra-organism” concept. By defining the metagenome metabolic network of 124 healthy, obese and with inflammatory bowel disease individuals, Jacobsen et al. (2013), could reveal clusters of high, medium and low metabolic potential. The abundance of specific genera was shown to be a key factor for the metabolic potential of the gut community. Recently, Garza et al. (2018) combined GMNs of individual species to predict the metabolic status of the whole community. Using an optimization framework that takes as input metagenomic species abundance levels and their GMNs, the authors could produce meaningful predictions for the metabolic environment. This approach provides a starting point for assessing the risk of a microbial infection based on the full composition of the environment. Recent pioneering studies providing a large number of high-quality GMNs for gut bacteria enable in silico studies of gut metabolic function and interaction (Magnúsdóttir et al., 2017; Machado et al., 2018). These resources have exemplified the study of gut microbes and their respective pairwise interactions but have not been used so far to study the interaction with pathogenic species. From these >1,000 human gut bacteria with GMNs available, we could create a “Potential Therapeutic Microbial Pool (PTMP)” that consists of those microbes that (i) are not pathogens or previously associated with diseases, and (ii) are more stable against perturbations in healthy individuals. Subsequently, and based on computed pairwise growth interactions with gut pathogens, such as C. difficile, E. coli, Enteroccocus spp., and Campylobacter spp., among others, we can identify microbes from the aforementioned PTMP that have a competitive relationship with the pathogen. Such computational frameworks could enable the design of microbiome-based prophylactic and therapeutic interventions for vulnerable-to-infections individuals.

Structure-Based Approaches to Analyze Enzymes in Pathogen-Host Interactions

PHI-based integrated metabolic approach can be applied to discover enzymes as potential drug targets. Computational subtractive genomics is one such approach to shortlist or prioritize enzymes in unique metabolic pathways of infectious pathogens as drug targets. The enzymes are filtered based on the non-host protein homology as well as essentiality for the survival of the pathogen. The structure-based methods such as homology modeling, molecular docking and molecular dynamics simulations are used to classify the candidate enzymes in terms of structure and function, as reviewed below.

Hossain et al. (2018) presented a very interesting application of the subtractive genomics method on a bacterial species Campylobacter jejuni RM1221 (CjR). The bacteria CjR is now resistant to many antibiotics. Therefore, there is a need to predict new drug targets based on unique and new pathways of the resistant CjR. In the first phase, the authors used the complete protein set of the resistant CjR to prioritize the essential proteins of CjR taking role in unique metabolic pathways. These unique pathways were pathogen specific only, therefore, absent in the host organism. In the next step, the authors filtered out the human homolog proteins. Finally, using this subtractive genomics strategy, they identified 38 eligible therapeutic targets, about 10 of which were metabolic enzymes. They studied corresponding 3D structural interactions with FDA-approved drugs. Additionally, they identified peptidoglycan biosynthesis as a pathway based drug target, to be explored by experimental methods. To help future research, a comprehensive database was also made, and it has 3D structure and other related data of the listed drug targets.

Sulfur metabolism pathway is one of the essential pathways of microorganisms to survive as a counter strategy to overcome the oxidative defense of a human host during the long latent phase of their life cycle, particularly, in anaerobic pathogens such as Entamoeba histolytica, Leishmania donovani, Trichomonas vaginalis, and Salmonella typhimurium. Since the Sulfur metabolism pathway is absent in human and all other mammals, the enzymes taking role in this pathway were best considered as drug targets and subjected to the structure based methods to predict potential strong binders (Mazumder and Gourinath, 2016). The Sulfur metabolism in microorganisms is through the cysteine biosynthetic pathway, which helps the microorganisms in their survival and pathogenicity. Cysteine biosynthesis pathway involves two major enzymes and these are: (i) serine acetyl transferase (SAT) and (ii) O-acetylserine sulfhydrylase (OASS). Among these two enzymes, the availability of the crystal structure of OASS made it best target to focus for the structure based computational studies. In addition to its absence in human host, the OASS is highly conserved between above-mentioned pathogens. In the study of Mazumder and Gourinath (2016) the authors have shown that the intermediate state of the enzyme could be considered as the best option to design inhibitor against OASS instead of closed and open conformation of this enzyme.

The pathogens and host relationship is governed by the Protein-Protein Interactions (PPIs) and several of those interactions can involve enzymes. With methicillin resistant Staphylococcus aureus (MRSA), which is the main source of nosocomial infection worldwide, Uddin et al. (2017) developed a network of PPI by using the Interolog method. The interolog method is based on the hypothesis that if two proteins interact in one species then there is a likely chance that a set of the ortholog of interacting proteins will also interact in the cross species. Hence, in the study of Uddin et al. (2017) the pathogen-host PPIs between MRSA and Humans were predicted as possible drug targets by considering various filters. A distant homolog approach was employed to find the MRSA homologs in human host using PSI-BLAST. The Database of Interacting Protein (DIP) was also used to find the direct homologs of the pathogen proteins in human host. As a result, the common pathogen-host PPIs were predicted and out of which the most repeated entry i.e., MRSA Histone Deactylase (HDAC) was proposed as a new drug target, against which one can propose new chemical compounds as potential inhibitors. The structure of the HDAC enzyme was modeled using the SWISS-MODEL and a small-scale chemical compounds library was docked on the binding site of the MRSA HDAC. This computational screening resulted in compounds that showed theoretical strong binding with the HDAC according to their estimated energetics. The predicted strong binders were proposed as potential inhibitors against MRSA HDAC.

Conclusion

With the birth of dual-omics technologies, we have entered a new era where interactions between pathogens and their host can be directly predicted to unravel molecular mechanisms of infectious diseases. Although integrated pathogen-host GMNs are currently scarce, we believe that the increase in the number of dual omics datasets, especially dual RNA-seq studies, will pave the way for the reconstruction of more and more integrated GMNs (Figure 1). Mapping dual-omics data from different states of infection on those models will lead to infection-specific prediction of metabolic phenotypes. Taking into account the crucial role of pathobiomes in infection, the balance between pathogens and microbiota in human gut is another dimension to analyze infectious diseases by using GMNs. We will witness reconstruction of integrated microbiota-pathogen GMNs or even multi models covering microbiota-pathogen-host GMNs to enable deeper interpretation of metagenomics datasets. The combined use of here-reviewed approaches for metabolism-oriented PHIs will certainly contribute positively to the drug research via identifying novel candidate enzymes as targets.

Author Contributions

TÇ conceived the study. TÇ, GP, RU, and SD wrote the manuscript.

Funding

This work was financially supported by TUBITAK, The Scientific and Technological Research Council of Turkey (Project Code: 316S005), by Deutsche Forschungsgemeinschaft (DFG) CRC/Transregio 124 Pathogenic fungi and their human host: Networks of interaction, subproject B5 and INF, and by PSF, The Pakistan Science Foundation [Project Code: PSF-TUBITAK/S-HEJ (04)].

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Baddal, B., Muzzi, A., Censini, S., Calogero, R. A., Torricelli, G., Guidotti, S., et al. (2015). Dual RNA-seq of nontypeable Haemophilus influenzae and host cell transcriptomes reveals novel insights into host-pathogen cross talk. mBio 6, e01765–e01715. doi: 10.1128/mBio.01765-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Bazzani, S., Hoppe, A., and Holzhütter, H. G. (2012). Network-based assessment of the selectivity of metabolic drug targets in Plasmodium falciparum with respect to human liver metabolism. BMC Syst. Biol. 6:118. doi: 10.1186/1752-0509-6-118

PubMed Abstract | CrossRef Full Text | Google Scholar

Bordbar, A., Lewis, N. E., Schellenberger, J., Palsson, B. Ø., and Jamshidi, N. (2010). Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions. Mol. Syst. Biol. 6:422. doi: 10.1038/msb.2010.68

PubMed Abstract | CrossRef Full Text | Google Scholar

Bumann, D. (2009). System-level analysis of Salmonella metabolism during infection. Curr. Opin. Microbiol. 12, 559–567. doi: 10.1016/j.mib.2009.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Cesur, M. F., Abdik, E., Güven-Gülhan, Ü., Durmuş, S., and Çakir, T. (2018). “Computational systems biology of metabolism in infection” in Metabolic Interaction in Infection, Experientia Supplementum, vol. 109, eds R. Silvestre and E. Torrado (Cham: Springer), 235–282. doi: 10.1007/978-3-319-74932-7_6

CrossRef Full Text | Google Scholar

Damron, F. H., Oglesby-Sherrouse, A. G., Wilks, A., and Barbier, M. (2016). Dual-seq transcriptomics reveals the battle for iron during Pseudomonas aeruginosa acute murine pneumonia. Sci. Rep. 6:39172. doi: 10.1038/srep39172

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunphy, L. J., and Papin, J. A. (2018). Biomedical applications of genome-scale metabolic network reconstructions of human pathogens. Curr. Opin. Biotechnol. 51, 70–79. doi: 10.1016/j.copbio.2017.11.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Durmuş, S., Çakir, T., Özgür, A., and Guthke, R. (2015). A review on computational systems biology of pathogen–host interactions. Front. Microbiol. 6:235. doi: 10.3389/fmicb.2015.00235

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandes, M. C., Dillon, L. A., Belew, A. T., Bravo, H. C., Mosser, D. M., and El-Sayed, N. M. (2016). Dual transcriptome profiling of Leishmania-infected human macrophages reveals distinct reprogramming signatures. mBio 7, e00027–e00016. doi: 10.1128/mBio.00027-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Fleming-Davies, A., Jabbar, S., Robertson, S. L., Asih, T. S. N., Lanzas, C., Lenhart, S., et al. (2017). “Mathematical modeling of the effects of nutrient competition and bile acid metabolism by the gut microbiota on colonization resistance against Clostridium difficile,” in Women in Mathematical Biology, eds A. Layton and L. Miller (New York, NY: Springer), 137–161. doi: 10.1007/978-3-319-60304-9_8

CrossRef Full Text | Google Scholar

Garza, D. R., van Verk, M. C., Huynen, M. A., and Dutilh, B. E. (2018). Towards predicting the environmental metabolome from metagenomics with a mechanistic model. Nat. Microbiol. 3, 456–460. doi: 10.1038/s41564-018-0124-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Griesenauer, B., Tran, T. M., Fortney, K. R., Janowicz, D. M., Johnson, P., Gao, H., et al. (2019). Determination of an interaction network between an extracellular bacterial pathogen and the human host. mBio 10, e01193–e01119. doi: 10.1128/mBio.01193-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Horswill, A. R., Dudding, A. R., and Escalante-Semerena, J. C. (2001). Studies of propionate toxicity in Salmonella enterica identify 2-methylcitrate as a potent inhibitor of cell growth. J. Biol. Chem. 276, 19094–19101. doi: 10.1074/jbc.M100244200

PubMed Abstract | CrossRef Full Text | Google Scholar

Hossain, M. U., Omar, T. M., Alam, I., Das, K. C., Mohiuddin, A. K. M., Keya, C. A., et al. (2018). Pathway based therapeutic targets identification and development of an interactive database CampyNIBase of Campylobacter jejuni RM1221 through non-redundant protein dataset. PLoS ONE 13:e0198170. doi: 10.1371/journal.pone.0198170

PubMed Abstract | CrossRef Full Text | Google Scholar

Humphrys, M. S., Creasy, T., Sun, Y., Shetty, A. C., Chibucos, M. C., Drabek, E. F., et al. (2013). Simultaneous transcriptional profiling of bacteria and their host cells. PLoS ONE 8:e80597. doi: 10.1371/journal.pone.0080597

PubMed Abstract | CrossRef Full Text | Google Scholar

Huthmacher, C., Hoppe, A., Bulik, S., and Holzhütter, H. G. (2010). Antimalarial drug targets in Plasmodium falciparum predicted by stage-specific metabolic network analysis. BMC Syst. Biol. 4:120. doi: 10.1186/1752-0509-4-120

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacobsen, U. P., Nielsen, H. B., Hildebrand, F., Raes, J., Sicheritz-Ponten, T., Kouskoumvekaki, I., et al. (2013). The chemical interactome space between the human host and the genetically defined gut metabotypes. ISME J. 7, 730–742. doi: 10.1038/ismej.2012.141

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacobson, A., Lam, L., Rajendram, M., Tamburini, F., Honeycutt, C., Pham, T., et al. (2018). A gut commensal-produced metabolite mediates colonization resistance to Salmonella infection. Cell Host Microbe 24, 296–307. doi: 10.1016/j.chom.2018.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacquet, R., LaBauve, A. E., Akoolo, L., Patel, S., Alqarzaee, A. A., Lung, T. W. F., et al. (2019). Dual gene expression analysis identifies factors associated with Staphylococcus aureus virulence in diabetic mice. Infect. Immun. 87, e00163–e00119. doi: 10.1128/IAI.00163-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Jamshidi, N., and Raghunathan, A. (2015). Cell scale host-pathogen modeling: another branch in the evolution of constraint-based methods. Front. Microbiol. 6:1032. doi: 10.3389/fmicb.2015.01032

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenior, M. L., Leslie, J. L., Young, V. B., and Schloss, P. D. (2017). Clostridium difficile colonizes alternative nutrient niches during infection across distinct murine gut microbiomes. mSystems 2, e00063–e00017. doi: 10.1128/mSystems.00063-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiedrowski, M. R., Gaston, J. R., Kocak, B. R., Coburn, S. L., Lee, S., Pilewski, J. M., et al. (2018). Staphylococcus aureus Biofilm growth on cystic fibrosis airway epithelial cells is enhanced during respiratory syncytial virus coinfection. mSphere 3, e00341–e00318. doi: 10.1128/mSphere.00341-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, H. U., Sohn, S. B., and Lee, S. Y. (2012). Metabolic network modeling and simulation for drug targeting and discovery. Biotechnol. J. 7, 330–342. doi: 10.1002/biot.201100159

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, P., Xu, Z., Sun, X., Yin, Y., Fan, Y., Zhao, J., et al. (2017). Transcript profiling of the immunological interactions between Actinobacillus pleuropneumoniae serotype 7 and the host by dual RNA-seq. BMC Microbiol. 17:193. doi: 10.1186/s12866-017-1105-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Machado, D., Andrejev, S., Tramontano, M., and Patil, K. R. (2018). Fast automated reconstruction of genome-scale metabolic models for microbial species and communities. Nucleic Acids Res. 46, 7542–7553. doi: 10.1093/nar/gky537

PubMed Abstract | CrossRef Full Text | Google Scholar

Magnúsdóttir, S., Heinken, A., Kutt, L., Ravcheev, D. A., Bauer, E., Noronha, A., et al. (2017). Generation of genome-scale metabolic reconstructions for 773 members of the human gut microbiota. Nat. Biotechnol. 35, 81–89. doi: 10.1038/nbt.3703

PubMed Abstract | CrossRef Full Text | Google Scholar

Mazumder, M., and Gourinath, S. (2016). Structure-based design of inhibitors of the crucial cysteine biosynthetic pathway enzyme O-acetyl serine sulfhydrylase. Curr. Top. Med. Chem. 16, 948–959. doi: 10.2174/1568026615666150825142422

PubMed Abstract | CrossRef Full Text | Google Scholar

McHan, F., and Shotts, E. B. (1993). Effect of short-chain fatty acids on the growth of Salmonella typhimurium in an in vitro system. Avian Dis. 37, 396–398. doi: 10.2307/1591664

PubMed Abstract | CrossRef Full Text | Google Scholar

Minhas, V., Aprianto, R., McAllister, L. J., Wang, H., David, S. C., McLean, K. T., et al. (2019). In vivo dual RNA-seq analysis reveals the basis for differential tissue tropism of clinical isolates of Streptococcus pneumoniae. bioRxiv 862755. doi: 10.1101/862755

CrossRef Full Text | Google Scholar

Muñoz, J. F., Delorey, T., Ford, C. B., Li, B. Y., Thompson, D. A., Rao, R. P., et al. (2019). Coordinated host-pathogen transcriptional dynamics revealed using sorted subpopulations and single macrophages infected with Candida albicans. Nat. Commun. 10:1607. doi: 10.1038/s41467-019-09599-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Niemiec, M. J., Grumaz, C., Ermert, D., Desel, C., Shankar, M., Lopes, J. P., et al. (2017). Dual transcriptome of the immediate neutrophil and Candida albicans interplay. BMC Genomics 18:696. doi: 10.1186/s12864-017-4097-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Nuss, A. M., Beckstette, M., Pimenova, M., Schmühl, C., Opitz, W., Pisano, F., et al. (2017). Tissue dual RNA-seq allows fast discovery of infection-specific functions and riboregulators shaping host–pathogen transcriptomes. Proc. Nat. Acad. Sci. U.S.A. 114, E791–E800. doi: 10.1073/pnas.1613405114

PubMed Abstract | CrossRef Full Text | Google Scholar

Olson, W. J., Stevenson, D., Amador-Noguez, D., and Knoll, L. J. (2018). Dual metabolomic profiling uncovers Toxoplasma manipulation of the host metabolome and the discovery of a novel parasite metabolic capability. bioRxiv 463075. doi: 10.1101/463075

CrossRef Full Text | Google Scholar

Petrucelli, M. F., Peronni, K., Sanches, P. R., Komoto, T. T., Matsuda, J. B., da Silva, W. A., et al. (2018). Dual RNA-Seq analysis of Trichophyton rubrum and HaCat Keratinocyte co-culture highlights important genes for fungal-host interaction. Genes 9:362. doi: 10.3390/genes9070362

PubMed Abstract | CrossRef Full Text | Google Scholar

Pittman, K. J., Aliota, M. T., and Knoll, L. J. (2014). Dual transcriptional profiling of mice and Toxoplasma gondii during acute and chronic infection. BMC Genomics 15:806. doi: 10.1186/1471-2164-15-806

PubMed Abstract | CrossRef Full Text | Google Scholar

Raghunathan, A., and Jamshidi, N. (2018). “Integrated host-pathogen metabolic reconstructions” in Metabolic Network Reconstruction and Modeling, ed M. Fondi (New York, NY: Humana Press), 197–217. doi: 10.1007/978-1-4939-7528-0_9

PubMed Abstract | CrossRef Full Text | Google Scholar

Rienksma, R. A., Schaap, P. J., dos Santos, V. A. M., and Suarez-Diez, M. (2019). Modeling host-pathogen interaction to elucidate the metabolic drug response of intracellular Mycobacterium tuberculosis. Front. Cell. Infect. Microbiol. 9:144. doi: 10.3389/fcimb.2019.00144

PubMed Abstract | CrossRef Full Text | Google Scholar

Sewankambo, N., Gray, R. H., Wawer, M. J., Paxton, L., McNaim, D., Wabwire-Mangen, F., et al. (1997). HIV-1 infection associated with abnormal vaginal flora morphology and bacterial vaginosis. Lancet 350, 546–550. doi: 10.1016/S0140-6736(97)01063-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Tucey, T. M., Verma, J., Harrison, P. F., Snelgrove, S. L., Lo, T. L., Scherer, A. K., et al. (2018). Glucose homeostasis is important for immune cell viability during Candida challenge and host survival of systemic fungal infection. Cell Metab. 27, 988–1006. doi: 10.1016/j.cmet.2018.03.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Uddin, R., Tariq, S. S., Azam, S. S., Wadood, A., and Moin, S. T. (2017). Identification of histone deacetylase (HDAC) as a drug target against MRSA via interolog method of protein-protein interaction prediction. Eur. J. Pharm. Sci. 106, 198–211. doi: 10.1016/j.ejps.2017.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Vayssier-Taussat, M., Albina, E., Citti, C., Cosson, J. F., Jacques, M. A., Lebrun, M. H., et al. (2014). Shifting the paradigm from pathogens to pathobiome: new concepts in the light of meta-omics. Front. Cell. Infect. Microbiol. 4:29. doi: 10.3389/fcimb.2014.00029

PubMed Abstract | CrossRef Full Text | Google Scholar

Westermann, A. J., Barquist, L., and Vogel, J. (2017). Resolving host–pathogen interactions by dual RNA-seq. PLoS Pathog. 13:e1006033. doi: 10.1371/journal.ppat.1006033

PubMed Abstract | CrossRef Full Text | Google Scholar

Zimmermann, M., Kogadeeva, M., Gengenbacher, M., McEwen, G., Mollenkopf, H. J., Zamboni, N., et al. (2017). Integration of metabolomics and transcriptomics reveals a complex diet of Mycobacterium tuberculosis during early macrophage infection. mSystems 2, e00057–e00017. doi: 10.1128/mSystems.00057-17

PubMed Abstract | CrossRef Full Text | Google Scholar