Khushboo Borah
Ye Xu
Johnjoe McFadden- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, United Kingdom
Tuberculosis (TB) is a devastating infectious disease that kills over a million people every year. There is an increasing burden of multi drug resistance (MDR) and extensively drug resistance (XDR) TB. New and improved therapies are urgently needed to overcome the limitations of current treatment. The causative agent, Mycobacterium tuberculosis (Mtb) is one of the most successful pathogens that can manipulate host cell environment for adaptation, evading immune defences, virulence, and pathogenesis of TB infection. Host-pathogen interaction is important to establish infection and it involves a complex set of processes. Metabolic cross talk between the host and pathogen is a facet of TB infection and has been an important topic of research where there is growing interest in developing therapies and drugs that target these interactions and metabolism of the pathogen in the host. Mtb scavenges multiple nutrient sources from the host and has adapted its metabolism to survive in the intracellular niche. Advancements in systems-based omic technologies have been successful to unravel host-pathogen interactions in TB. In this review we discuss the application and usefulness of omics in TB research that provides promising interventions for developing anti-TB therapies.
Introduction
Tuberculosis (TB) remains a global pandemic and the biggest infectious killer despite being a preventable and curable disease (1). This communicable disease is caused by the pathogen Mycobacterium tuberculosis (Mtb) and it kills over a million people every year. In 2016, WHO estimated around 1.7 billion people were reported to have latent TB infection (1, 2). Drug resistance in TB has been an escalating problem worldwide with nearly half a million people developing rifampicin resistance in 2019. The pathogenic mycobacterium becomes resistant to the first line anti-TB drugs including isoniazid, rifampicin, fluoroquinolone and second line drugs causing multidrug resistant (MDR) and extensively drug resistant (XDR) TB, which are a serious global threat with MDR-TB accounting for one third of the deaths worldwide due to antimicrobial resistance (1, 3). MDR and XDR-TB cases are developed mainly due to poor drug regimen, treatment misuse and poor patient compliance due to the lengthy anti-TB treatments and the associated side-effects of the medications (3). This rapidly growing problem of drug resistance needs urgent attention. We need to develop new therapies to overcome the limitations of current anti-TB drug treatments. We also need improved and efficient diagnostic platforms for detection and surveillance for disease control and management of MDR and XDR-TB. COVID-19 pandemic posed significant challenges and disruptions to TB control and cure by restricting the health services, infrastructure, workforce and research efforts which were diverted away from TB and other diseases (4, 5). These disruptions from COVID-19 pandemic are predicted to increase global TB deaths by 20% over the next 5 years (4). We need to harness global research efforts to develop new and effective treatments for TB and to mitigate the limitations on treatment and cure imposed by the COVID-19 pandemic for avoiding additional morbidity and mortality.
The causative agent of TB has evolved to persist in humans. Mtb interacts with alveolar macrophages which are the primary cellular site of infection (6, 7). Macrophages in the lungs of an infected individual phagocytize TB bacilli which were inhaled as aerosols from another infected individual (6). Inside macrophages Mtb are contained in phagosomes which exhibit bactericidal activity through maturation and acidification (8, 9). Macrophages produce reactive oxygen and nitrogen species to eliminate intraphagosomal Mtb (6, 10). Macrophages initiate a pro-inflammatory response with other immune cells being recruited to the site of infection and forming a granulomatous lesion to arrest bacterial replication and stop dissemination. However, Mtb has evolved adaptation mechanisms to avoid the hostile host environment and immune defences including intraphagosomal survival and escape into the cytoplasm (9, 11). Mtb adapts its life and metabolism to survive in the hostile phagosome environments including hypoxia, low pH, antimicrobial peptides, reactive oxygen and nitrogen species and nutrient limitations (12–14). Metabolism of the TB pathogen is key for its survival in the human host. Decades of research using in vivo and in vitro TB models has enabled identification of the metabolic characteristics of Mtb within the intracellular niche. It is now well-established that Mtb uses multiple host metabolites such as lipids and amino acids as nutrient sources during infection (6, 13–15). Therefore, identifying vulnerable metabolic targets of Mtb provide avenues for the development of new antimicrobials and it has been an accelerating research area. Scientific advances in the development and application of systems-wide tools and platforms to screen genome, transcriptome, proteome, metabolome, metabolic modelling and fluxome have revolutionised research in untangling the interactions between human host and Mtb (Figure 1). In this review, we discuss advancements in TB research achieved using systems-based omics.
Figure 1 Overview of systems-based omics used for exploring Mtb and host metabolic interactions. Created with BioRender.com.
Genomics, Transcriptomics, and Proteomics
The discovery of next-generation sequencing (NGS) facilitates the genomic and transcriptomic research of mycobacteria. The genomic or transcriptomic analysis focuses on the variation of genomic features, such as gene expression level, DNA sequencing or regulatory elements annotation, which enables the identification of essential genetic or regulatory targets under certain conditions (16–18). Resistance mutation of mycobacterium to drugs can be detected by genomic analysis of clinical isolates (19). The massive online database generated from genomics analysis of resistance mutation enables the generation of genetic interaction networks to different antibiotics (18, 20). RNAseq and methylome analysis of different clades of Mtb complex detected differential gene expression involved in host interaction and metabolism, which is further linked to the varied phenotypes and host susceptibility between the Mtb complex (21). RNAseq has been used to identify the relevant virulent genes in clinical Mtb strains such as the MKR mutant (17). Genes involved in cholesterol degradation as well as ESX-1 secretion system were up-regulated in the MDK strain but not in Mtb H37Rv (17). Drugs targeting the cellular pathway could decrease the intracellular survival of the MDK strain, which further proves the application of transcriptional analysis of Mtb strains could facilitate the identification of drug targets. Some studies combining both genomics and transcriptomic analysis among the same set of mycobacterial isolates identify specific single nucleotide polymorphisms (SNPs) and varied gene expression levels, as well as the link to phenotypes of the bacterial population (22, 23). Genomics or transcriptomics analyses are adopted to investigate the influence of deletion or over-expression mutants on the gene expression of the whole population to identify relevant cellular pathways (24, 25). Genomics and transcriptomics have provided insights into TB pathogenesis and host-pathogen interactions and facilitated identification of biomarkers of TB and vaccine development (26–29). RNA seq analyses by Kaufmann et al., demonstrated changes in the transcriptome profiles of hematopoietic stem cells (HSCs) and multipotent progenitors in BCG intravenously vaccinated mice vs. control mice. These epigenetic changes in the HSCs enhanced myelopoiesis and generated BCG trained macrophages that provided enhanced protection against Mtb (27). Esterhuyse et al. employed DNA methylome, transcriptome and proteome analyses of monocytes and granulocytes to demonstrate differences in DNA methylation profiles which correlated to changes in immune cell activation, transcription and inflammation between active TB and latent TB infection (LTBI) patients (29). A recent study using transcriptomic analyses compared the lung biosignatures between TB, lung adenocarcinoma (LUAD) and sarcoidosis patients and identified MK167 which is a mediator of Mtb-promoted tumor cell proliferation, migration, and invasion to be overexpressed in both TB and LUAD patients (26).
Proteomic analyses are used to identify and quantify proteins, which are functional biochemical entity of an organism, and these analyses provide direct functional information complementary to the genomics and transcriptomics analyses (30). Development and application of proteomic platforms has been researched for over a decade. In TB research, systems-level proteomic analysis has been applied with the motive to develop diagnostic markers, vaccines and therapeutics (31). Proteomic methods involve gel electrophoresis separation of proteins (30, 32) and liquid chromatography-mass spectrometry (LC-MS/MS) analyses of protein samples which can be conducted, broadly in two ways. A widely used approach, also known as bottom-up proteomics, involves digestion of proteins into peptides, followed by their identification and mapping into the proteins using their mass-to-charge ratios and further fragmentation of the ions for quantitative analyses (30, 33). The second approach, called top-down approach involves fragmentation of the total intact protein into molecular ions followed by protein identification and quantitation (31). Proteomic studies in TB have provided information on the TB pathogen’s protein expression profiles in in vitro and within the host. Discovery proteomic analyses by Rosenkrands et al. (34) identified 82 novel proteins out of 49 extracellular culture filtrate and 118 cellular lysate proteins in Mtb. Comparative proteomic analyses in clinical isolates and virulent and avirulent mycobacterial vaccine strains and Mtb complexes (MTBC) have revealed strain-specific characteristics (35, 36). Jungblut et al. (35) identified six proteins including L-alanine dehydrogenase (Rv2780), isopropyl malate synthase (Rv3710), nicotinate-nucleotide pyrophosphatase (Rv1596), MPT64 (Rv1980c) and two conserved hypotheticals (Rv2449c and Rv0036c) with missing counterparts in M. bovis BCG (Chicago strain). Chicago and Copenhagen BCG strains exhibited highly similar proteomes with only three identified variants while Mtb H37Rv and Erdman strains had 18 variants between them (35). Clinical strains JAL and BND showed distinct variations in the Esx and mce1 operon proteins, which contributed to their virulence and drug resistance, as compared with H37Ra and H37Rv strains (37). A tandem mass tag (TMT) labelled proteomics identified differential protein expression and phosphorylation in PE/PPE/PE-PGRS between H37Rv and H37Ra which contributes to the virulence of H37Ra in human cells (38). Differences in expression of nitrate metabolism proteins between Mtb H37Rv and drug-susceptible Beijing and multidrug-resistant Beijing strains were identified to be important in the pathogens’ adaptation to stressful intracellular environments (39). The non-replicative persistent (NRP) state of Mtb was investigated using isotope coded affinity tag-based (ICAT) proteomics which uses isotopic labels for quantitation of proteins (40). Cho et al. (40) cultivated Mtb in an oxygen depleted fermentor to achieve early and late NRP states of cultures and measured relative expression of proteins. This proteomic analysis revealed different expression profiles in the two NRP states which were associated with energy metabolism and degradation (40). In guinea pig models of TB, shot gun tandem mass spectrometry (MS/MS) proteomics at early (30 days-post infection) and chronic (90 days-post infection) infection stages identified over 500 Mtb proteins in infected lung tissues (41). This research identified heterogeneity in two protein classes, belonging to the cell wall processes and respiration and metabolism between the two stages highlighting these processes as necessary adaptations in persistent Mtb (41).
Proteomic approaches have been used to investigate host-pathogen interactions and immunological responses in TB, and in identification of diagnostic markers. The cell wall of Mtb is important for its virulence the cell wall lipids and proteins are antigenic and elicit host immunomodulations. Lipoproteins, T and B cell antigens and proteins associated with small and macromolecule metabolism were identified in Mtb’s cell wall which likely facilitated the transport of nutrients between the cytosol and extracellular milieu and cell wall re-modelling (42). Antibody screening and immunological responses using systems level proteomics has attracted great attention in recent years. TMT based quantitative proteomics showed differentially expressed proteins in THP-1 macrophages infected with Mtb H37Rv and H37Ra strains, and include proteins involved in apoptosis, blood coagulation and oxidative phosphorylation providing evidences to strain-specific host macrophage responses (43). T-cell IFN- Ɣ immunological response elicited by Mtb antigens in culture filtrate and cellular extracts were screened in splenocytes from Mtb-infected mice which allowed identification of 17 novel T-cell antigens (44). Proteomic analyses by Penn et al. using affinity tag purification mass spectrometry mapped 187 Mtb-human protein-protein interactions and identified two factors, Mtb’s secreted protein LpqN and the human ubiquitin ligase CBL involved in Mtb’s pathogenesis (45). Robust quantification of proteins in complex samples representing the intracellular milieu is critical to not only advance the TB biology but also to develop therapeutic and diagnostic interventions. Schubert et al. (46) developed Mtb proteome library, a public resource of definitive MS assays using single reaction monitoring (SRM) technique to quantify proteins in complex biological samples. Label-free quantitation of plasma proteins in patients with pulmonary (PTB) and (LTBI) provided a diagnostic model which showed alpha-1-antichymotrypsin(ACT), alpha-1-acid glycoprotein 1 (AGP1), and E-cadherin (CDH1) with >80% sensitivity, specificity and accuracy in distinguishing LTBI from PTB (47). Another study using TMT-based quantitative proteomics of human serum demonstrated different signatures of inflammatory proteins and apolipoprotein A and serotransferrin proteins (involved in lipid transport and iron metabolism) in LTBI and active TB cohorts (48). A multidimensional and stable isotope labelled q3D LC-MS quantitative plasma proteomics enabled discovery of novel protein biomarkers and 5-protein signature comprising of complement factor H related 5 (CFHR5), interleukin enhancer binding factor 2 (ILF2), leucine-rich alpha-2 glycoprotein (LRG1), LPS-binding protein (LBP), serum amyloid A (SAA), plasma C-reactive protein (CRP) and E3 ubiquitin-protein ligase listerin (LTN) to improve diagnostic accuracy of TB (49). Immunological proteome screening using integrative host and pathogen biochemical datasets allowed identification of antibodies in human serum of TB patients. A large-scale screening used protein microarray platforms which included 4099 Mtb proteins; these microarrays were probed with serum from over 500 individuals with and without active TB. Antibodies were screened against the entire Mtb proteome to identify immunoproteome for active TB which were extracellular proteins comprising 0.5% of the entire proteome (50) demonstrating well-established application of proteomics in identification of host-pathogen interactions for serodiagnostic markers.
Metabolomics and Isotopic Tracing Studies
Metabolomics is the one of the newest omic technologies used to discover, identify, and quantify metabolites which are biochemical entities produced and consumed in metabolic reactions which in turn drives metabolism and energy production in biological systems. Metabolomics allow big data integration across different omics for comprehensive determination of the consequences of all metabolites on cellular function and physiology (51–53). Metabolomics combined with isotope labelling strategies have identified nutrient sources for Mtb in in vivo environments and measured their uptake and assimilation by Mtb metabolic network; this is not feasible with any other omic platforms. Metabolomic analysis involves quenching of metabolites from living cells and cellular extracts and polar/nonpolar metabolites recovery, followed by analytical identification and quantification. Gas chromatography-mass spectrometry (GC-MS), nuclear magnetic resonance (NMR) and LC-MS are three most widely used analytical platforms for metabolomics. Depending on the purpose of the study, metabolites are analysed either by untargeted or unbiased methods/platforms, mainly used for qualitative analysis or by targeted analyses involving SRM or multiple reaction monitoring (MRM). Metabolism is at the heart of survival, growth, and virulence of Mtb, and other pathogens inside the human host cells. Metabolism of immune cells is equally important in driving a variety of cellular and immunological responses during infection and therefore, in recent years, immunometabolism has emerged as a rapidly growing area with the motive to develop new diagnostics and therapeutics for TB. The advancements in analytical tools and platforms have also allowed metabolomics and metabolic investigations at single cell and subcellular (cellular microcompartment) levels (54, 55). In this section, we have reviewed the latest developments in metabolomic technologies and their application to achieving breakthroughs in characterisation of the metabolic cross talk between the host and Mtb.
Metabolomics has been extensively used to identify metabolic pathways and to decipher metabolism of Mtb in in vitro and in vivo animal models. Activity-based metabolomic profiling which uses cellular metabolome as the platform using protein and time-dependent production and consumption of small molecules, assigned functional and structural characteristics of Rv1692 as a glycerol 3-phosphate phosphatase in glycerophospholipid metabolism (56). Rv322c has been recently assigned a nitrogen metabolic role in aspartate biosynthesis and this enzyme is essential for Mtb’s survival in murine macrophages and in mice (57). Lineage-specific metabolic differences, predominantly in amino acid, glycolysis, and tricarboxylic acid (TCA) cycle metabolism between six strains of MTBC complex were identified in exometabolomes of culture supernatants by untargeted time-of-flight (TOF) MS analysis. This study involved integration of metabolomics data into a constraint-based model which predicted the SNPs and genetic basis of the metabolic variations (58).
Complementary to other omic platforms, metabolomics has been applied to successfully identify Mtb and host metabolites in human biological fluids. High resolution orbitrap MS identified lipid metabolites such as phospha-tidylglycerol (16:0_18:1), lysophosphatidylinositol (18:0) and acylphosphati-dylinositolmannoside (Ac1PIM1) elevated in plasma of patients with active pulmonary TB (59). Another study applied orbitrap MS/MS and measured host responses to active TB and LTBI and identified elevated tryptophan catabolism to kynurenine mediated by indoleamine 2,3-dioxygenase-1 (IDO-1) (60). Metabolic biosignatures measured in human urine samples identified pre and post therapy clinical responses to TB antimicrobials which included 23 molecular features changing during the treatment (61). Drug phenotypes and their mechanism of Mtb's growth inhibition has been investigated using metabolomics and isotopic labelling analyses. Bedaquiline (BDQ), a newly approved anti-TB drug exhibits Mtb killing through inhibition of ATP synthase. The mode-of-action of BDQ is complex and has been investigated by multiple research demonstrating that BDQ induces Mtb’s glycolytic and gluconeogenic metabolic dependencies and that pyruvate phosphate dikinase (PPDK) is a vulnerable metabolic target (62). Knoll et al. (63) measured metabolite markers involved in TCA cycle metabolism, cell wall and DNA synthesis with GCxGC-TOF-MS in Mtb grown with and without ciprofloxacin, (a quinolone potent against drug resistance TB) to identify mechanism of action and Mtb’s adaptation to antimicrobials (63). Comparison of drug resistant and drug susceptible metabolic profiles showed proline and isoleucine levels significantly reduced in drug resistant strains (64).
Host and pathogen interaction in infectious diseases has been investigated in detail by multiple studies across decades using metabolomics combined with isotopic tracing studies. Metabolite profiling and 13C tracer studies have re-assigned identity and functionality of Mtb enzymes and elaborated its central carbon metabolic (CCM) pathways such as Rv1248c as 2-hydroxy-3-oxoadipate synthase catalysing C-C bond formation between α-ketoglutarate and glyoxylate (65), the role of phosphoenolpyruvate carboxykinase (PEPCK) in gluconeogenesis and its essentiality for Mtb’s survival in mice models (66) and the bifunctional role of Rv0812 coupling nucleic acid and cell wall biosynthesis (67). Untargeted metabolite profiling and isotopic analysis of Mtb grown on 13C-labelled carbon substrates dextrose, acetate and glycerol demonstrated that Mtb catabolizes multiple carbon substrates simultaneously to achieve monophasic growth providing a snapshot of metabolic adaptations in Mtb to host intracellular niche (68). GC-MS based isotopologue profiling using 13C carbon isotopic labels of Mtb infected-human THP-1 macrophages measured utilization of glucose by Mtb and essentiality of the anaplerotic node enzymes pyruvate carboxylase (PCA), PEPCK, malic enzyme (MEZ), and PPDK in gluconeogenesis, propionate and cholesterol detoxification and lipid synthesis (15, 69). Similar isotopologue analyses have been successfully applied across various host-pathogen models to measure intracellular metabolism of these pathogens, such as 13C glucose profiles of M. leprae growing in Schwann cells (70) and Trypanosoma brucei in bloodstream form (71) and 13C glycerol profiles in legionella pneumophila replicating in macrophages (72). 13C-isotopic analysis of Mtb in vitro cultures in a chemostat setup at controlled growth rate demonstrated efficient isotope incorporation and co-catabolism of 13C-labelled acetate and cholesterol, which are intracellular nutrient sources for Mtb (73). Recent advances in MS technologies have enabled the use of several stable isotopes including both single and mixed labelled species such as 15N and 2H. 2H cholesterol labelling of Mtb in vitro cultures confirmed that Mtb utilized cholesterol to synthesize amino acids (73). 15N-asparagine was used to measure the uptake and assimilation of nitrogen in murine macrophages using NANO-Secondary ion mass spectrometry (SIMS) (74). Agapova et al. (75) performed an LC-MS and 15N-isotopologue analysis to measure ammonium and amino acid utilization in in vitro Mtb and demonstrated that amino acids were preferred over ammonium and were utilized at similar rates, and that alanine dehydrogenase ald was important for alanine utilisation as a nitrogen source (75). 15N-isotopic labelling and GC-MS isotopologue analyses of Mtb infected-THP-1 macrophages identified nitrogen assimilation of amino acids in intracellular Mtb replicating human cells and further identified phosphoserine transaminase serC as a potential drug target (13).
Constraint-Based Modelling
Metabolism in a biological system involves numerous complex biochemical processes. Each biochemical reaction involves consumption and production of metabolites or chemical species. Metabolic flux through a reaction defines the rate at which a substrate is utilized and is the accurate measure of the activity of a pathway in the network. Metabolomics provides a measure of the pool sizes of the metabolites that can indicate active pathways but cannot measure the fluxes through the pathways. Therefore, measuring metabolic fluxes is the way to confirm network activity and function, and define the metabolic phenotype of a system.
Mathematical modelling of Mtb is a recent innovation in TB research and is popular for prediction of genes, metabolic and drug phenotypes of the pathogen and host-pathogen interactions. These models include a network of biochemical reactions including enzymatic, spontaneous or transport reactions which are reconstructed from annotated genome and literatures and implements computational methods such as flux balance analysis (FBA), flux variability analysis (FVA) and gene essentially predictions (76). Constraint based in silico modelling comprises of stoichiometric, mass, charge and energy balanced reactions and is used for prediction of steady-state phenotypes. The models include an objective function, typically a biomass equation. Flux distributions through the network are computationally optimized to maximize biomass function and growth rate to predict phenotypes. The accuracy of these predictions needs to be further validated using experimental approaches. These genome scale models are powerful as these platforms can be used to integrate transcriptomic, proteomic, and metabolic data to provide a holistic view of an organism’s metabolic and biochemical network in various environmental and nutritional conditions, stress and disease which is otherwise very challenging to derive from experimental and previously discussed omic platforms alone. Such metabolic modelling has been applied to study metabolic flexibility of Mtb and to predict mutant phenotypes (76–78). GSMN-TB, the first genome scale model of Mtb predicted the requirement of the enzyme isocitrate lyase during slow growth of Mtb in a continuous culture (78). Bonde et al. developed a differential producibility analysis (DPA) algorithm to extract metabolite production profiles through integration of micro array gene expression data into GSMN-TB model (79). DPA using GSMN-ML, the first genome scale model of M. leprae, was applied to interrogate RNA-seq data to derive in vivo metabolite production of leprae bacillus and nutritional status (80). Metabolic variations within the Mtb complex were identified through comparisons of substrate utilisation rates, gene essentiality data and growth rates calculated from GSMN-TB, GSMN-MB (M. bovis) and GSMN-BCG (M. bovis BCG) models (81). Lofthouse et al. tested a range of carbon and nitrogen sources to compare substrate utilisation between the models; the experimental validations identified several discrepancies with the predictions, highlighting the need to verify in silico predictions and to continually update and curate the models with updated genetic and enzymatic knowledge. Recent integration of metabolomic and genomic data into the genome scale models of Mtb complex predicted genetic mutations and their associated metabolic vulnerabilities in 18 clinical strains and providing an alternative approach to the genomic approaches such as GWAS to relate genetic variations to clinical phenotypes, and additionally, predicting related metabolic variations (53). Genome scale modelling has been used to identify drug targets through gene essentiality predictions (77). Rienksma et al. modelled metabolic state of Mtb during infection using condition-specific objective function and predicted nutrient uptake and gene essentiality to identify vulnerable drug targets in Mtb (76, 77). Till date, there are 16 genome scale models that vary in network topology and biomass equations, and these inconsistencies amongst the existing models poses limitations in their use and applications. It is therefore important that these models are updated and are standardized for predicting physiologically relevant phenotypes. Kavvas et al. (82) unified sMtb and iOSDD, two of the newest reconstructions and generated iEK1011, a standardized model that explored metabolic phenotypes of Mtb under various environmental conditions (82). Recently López-Agudelo et al. (83) conducted useful comparisons between eight Mtb models and identified the two best sMtb and iEK1011 models for metabolic state predictions (83). Although constraint-based genome scale modelling provides useful predictions on metabolic network operation and gene essentialities, these predictions are calculated based on the maximization or minimization of the assumed objective function in the model and therefore may not be an exact representation of the metabolic phenotype (84). The predictions need to be validated through measurement of metabolically active fluxes supported by experimental data.
Metabolic Flux Analysis and Fluxomics
Metabolic Flux Analysis (MFA), which is an amalgamation of computational and experimental techniques, provides precise quantification of metabolic fluxes in the network. The general workflow of steady state MFA is elaborated in Figure 2. The experimental part of MFA involves cultivation of microbial or eukaryotic cells in a growth medium containing isotopically labelled substrates. Majority of the MFA studies are conducted on steady state cultures where the cells are grown at a controlled growth rate using a bioreactor/chemostat setup (73, 85, 86) or assuming a pseudo-steady state during exponential growth mainly for eukaryotic cells (87, 88). Metabolic steady state of the system is confirmed through measurement of substrate uptake, CO2, and biomass production rates (73). Metabolically steady state cultures are next checked for isotopic steady state through addition of isotopic substrates and their incorporation into the biomass at various time points of cultivation. Cells at isotopic steady state are harvested for extraction of metabolites and follow up metabolomic and mass isotopomer distribution (MID) analyses using mass spectrometry (GC-MS, LC-MS) or NMR. The computational part of MFA involves construction of a metabolic model of the system under study; the model includes atomic transitions for the central metabolic network reactions, free and measured fluxes and a biomass reaction constructed from experimentally derived measurements of macromolecular compositions required to constrain the model. The MIDs are included in the model as the experimental data set, and the model is simulated using linear programming to iteratively fit the experimental MIDs and derive intracellular fluxes. Several computational tools such as 13C-FLUX2 (86), Isotopomer Network Compartmental Analysis (INCA) (73), WUFlux (89) are available for MFA calculations and for statistical evaluations of measured fluxes.
Figure 2 General workflow of steady state 13C-Metabolic Flux Analysis (MFA). The application is shown for measuring metabolic fluxes of Mtb in vitro in chemostat and during intracellular growth in macrophages. The methodology includes growth of the biological system in isotopically labelled media, followed by achievement of metabolic and isotopic steady state for chemostat set up or assuming pseudo steady state from macromolecular measurements at various time points of isotopic labelling. Metabolomics and mass isotopomer distributions (MIDs) are measured using GC-MS, LC-MS or NMR. A minimum of 106-107 cells is required for robust MID measurements. The computational part of 13C-MFA includes construction of a metabolic model consisting of atomic transitions for central metabolic reactions. The model is constrained with extracellular flux and biomass measurements. Measured MIDs are incorporated into the model and the model is cimulated with MFA computational platform to derive best-fit metabolic fluxes that defines the metabolic phenotype of the system under study. Created with BioRender.com.
Steady state 13C-MFA, which uses 13C-isotopically labelled substrates has been used to quantify intracellular central carbon metabolic fluxes of Mtb in vitro and in human THP-1 macrophages. 13C-MFA of Mtb in vitro cultures cultivated at slow and fast growth rate using 13C-glycerol as the labelled carbon substrate identified a novel pathway for pyruvate dissimilation through glyoxylate shunt and anaplerotic reactions and demonstrated the ability of Mtb to fix CO2 for synthesizing biomass (86). Recently, the intracellular carbon fluxes that support co-catabolism of multiple carbon substrates by Mtb were identified using 13C-MFA (73). This study demonstrated that flux partitioning between the TCA cycle and glyoxylate shunt with a reversible methyl citrate cycle enables Mtb to co-catabolise cholesterol and acetate. 13C-MFA derived metabolic flux profiles of BDQ-treated Mtb identified the operation of a bifurcated TCA cycle and requirement of the anaplerotic node and methylcitrate cycle, thereby elaborating the metabolic changes in BDQ-treated Mtb (62). Application of classical 13C-MFA to measure metabolic fluxes of Mtb inside human cells is challenging as the system is at non-steady state. To overcome this challenge Beste et al. (15) developed 13C-Flux Spectral Analysis (FSA), a computational tool that probed the uptake and utilisation of carbon sources including a mixture of amino acids, C1 and C2 compounds by Mtb from human THP-1 macrophages (15). In addition to carbon, Mtb utilises several nitrogen sources during infection. This was recently demonstrated by Borah et al. (13) through the application of 15N-isotopically labelled amino acids to Mtb-THP-1 macrophage and development of 15N-Flux Spectral Ratio Analysis (FSRA) to identify the nitrogen sources acquired by Mtb from the host cell and their intracellular assimilation (13). This research demonstrated glutamine as the primary nitrogen donor for intracellular Mtb and that serine biosynthesis is essential for survival of Mtb in human macrophages.
Future Perspectives and Conclusions
Till date, several omic technologies have been developed and applied to investigate host-pathogen interactions in TB and in other diseases. Each of the five omics discussed in this review uses unique tools and methodology and specializes in the measurement of particular type of biomolecules and their functions. Genomic technologies such as GWAS has led to the discovery of new biological mechanisms in TB. Transcriptomic approaches are powerful tools to derive host and pathogen molecular signatures and for identification of biomarkers for diagnosis and prognosis. However, both these approaches have their limitations that needs to be considered while designing future studies. There are inconsistencies in genomic studies due to the failure in replicating genetic associations across various studies and experimental settings; previous transcriptomic studies have highlighted the need to conduct comparisons across diverse disease cohorts and between whole blood, tissues, and cells, thus limiting the precise identification of genetic signatures between active and latent TB (90, 91). Proteomics, the complementary approach to genomics and transcriptomics which has been used to identify host-Mtb interactions and pathogen physiology, also suffers from limitations such as low instrument sensitivity (92). Although metabolomics provides the advanced omic platform for detection of small molecules, the precision of detection and measurement are affected by extraction method and conditions. There are variations in metabolomic analyses performed across research labs due to the variations in instrument, sample preparations, data integration and statistical analysis which currently limits the reproducibility, sensitivity and specificity required for clinical applications (93). The most recent fluxomic technology measures the metabolic phenotype of a biological system and provides a systems-wide identification of active pathways. However, there exists limitations in application of fluxomics to TB disease cohorts, firstly due to the unavailability of a human TB metabolic and mathematical model of disease and secondly, fluxomics is currently limited to measuring mainly central carbon fluxes in host and pathogen. Isotopic labelling studies in Mtb are mainly focussed on single isotope species. Co-labelling studies with different isotopes such as 13C and 15N are limited mainly because of the low sensitivities of MS instruments to precisely distinguish between different isotopic species. As evident from several studies, there are multiple nutrients that the TB pathogen utilises inside the human host. To gain a complete picture of the metabolic characteristics of Mtb pathogen in humans, there is a need to further develop the fluxomic technologies to model multiple substrate utilisation and carbon and nitrogen co-metabolic reactions. There is a need to develop metabolomics and MS sensitivities to precisely distinguish and measure isotopic species-specific MIDs obtained from labelling experiments conducted with multiple isotopic substrates. Finally, the flux computational platform needs to be developed to accurately model fluxes for the utilisation of multiple atoms through the metabolic network.
In conclusion, systems-based omic technologies have advanced the frontiers of host and pathogen biology and their interactions in TB disease. Over decades, there has been several successful attempts in developing and applying various omic platforms for the identification of diagnostic markers and therapeutic interventions. However, the current limitations in these omic approaches limit the investigation of host-pathogen molecular interactions in clinical settings and to reproducibly identify and validate diagnostic and therapeutic markers. The integration of various omic approaches, data sharing, cross lab validations and application to studying multiple disease cohorts are plausible ways forward to harness the calibre of various omic approaches for discovering clinically relevant biological mechanisms in TB.
Author Contributions
KB, conceptualisation, writing, reviewing and editing, and funding acquisition. JM, reviewing and editing and funding acquisition. YX, writing and editing. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Biotechnology and Biological Sciences Research Council (BBSRC) grant, BB/V010611/1 United Kingdom. The funders have no role in preparation of this work.
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. World Health Organization. Global Tuberculosis Report 2020. Geneva: World Health Organization (2020).
2. Cohen A, Mathiasen VD, Schön T, Wejse C. The Global Prevalence of Latent Tuberculosis: A Systematic Review and Meta-Analysis. Eur Respir J (2019) 54(3):1900655. doi: 10.1183/13993003.00655-2019
3. Iacobino A, Fattorini L, Giannoni F. Drug-Resistant Tuberculosis 2020: Where We Stand. Appl Sci (Switzerland) (2020) 10(8):1342–60. doi: 10.3390/app10062153
4. Hogan AB, Jewell BL, Sherrard-Smith E, Vesga JF, Watson OJ, Whittaker C, et al. Potential Impact of the COVID-19 Pandemic on HIV, Tuberculosis, and Malaria in Low-Income and Middle-Income Countries: A Modelling Study. Lancet Global Health (2020) 8:e1132–41. doi: 10.1016/S2214-109X(20)30288-6
5. Khurana AK, Aggarwal D. The (in)Significance of TB and COVID-19 Co-Infection. Eur Respir J (2020) 56(2):2002105. doi: 10.1183/13993003.02105-2020
6. Russell DG, VanderVen BC, Lee W, Abramovitch RB, Kim MJ, Homolka S, et al. Mycobacterium Tuberculosis Wears What It Eats. Cell Host Microbe (2010) 8:68–76. doi: 10.1016/j.chom.2010.06.002
7. Cohen SB, Gern BH, Delahaye JL, Adams KN, Plumlee CR, Winkler JK, et al. Alveolar Macrophages Provide an Early Mycobacterium Tuberculosis Niche and Initiate Dissemination. Cell Host Microbe (2018) 24:439–46.e4. doi: 10.1016/j.chom.2018.08.001
8. Prashar A, Schnettger L, Bernard EM, Gutierrez MG. Rab GTPases in Immunity and Inflammation. Front Cell Infect Microbiol (2017) 7:435. doi: 10.3389/fcimb.2017.00435
9. Queval CJ, Brosch R, Simeone R. The Macrophage: A Disputed Fortress in the Battle Against Mycobacterium Tuberculosis. Front Microbiol (2017) 8:2284. doi: 10.3389/fmicb.2017.02284
10. Matta SK, Kumar D. Hypoxia and Classical Activation Limits Mycobacterium Tuberculosis Survival by Akt-Dependent Glycolytic Shift in Macrophages. Cell Death Discov (2016) 2:16022. doi: 10.1038/cddiscovery.2016.22
11. Quigley J, Hughitt VK, Velikovsky CA, Mariuzza RA, El-Sayed NM, Briken V. The Cell Wall Lipid PDIM Contributes to Phagosomal Escape and Host Cell Exit of Mycobacterium Tuberculosis. mBio (2017) 8(2):e00148–17. doi: 10.1128/mBio.00148-17
12. Gengenbacher M, Kaufmann SHE. Mycobacterium Tuberculosis: Success Through Dormancy. FEMS Microbiol Rev (2012) 36:514–32. doi: 10.1111/j.1574-6976.2012.00331.x
13. Borah K, Beyß M, Theorell A, Wu H, Basu P, Mendum TA, et al. Intracellular Mycobacterium Tuberculosis Exploits Multiple Host Nitrogen Sources During Growth in Human Macrophages. Cell Rep (2019) 29:3580–91.e4. doi: 10.1016/j.celrep.2019.11.037
14. Gouzy A, Poquet Y, Neyrolles O. Nitrogen Metabolism in Mycobacterium Tuberculosis Physiology and Virulence. Nat Rev Microbiol (2014) 12:729–37. doi: 10.1038/nrmicro3349
15. Beste DJV, Nöh K, Niedenführ S, Mendum TA, Hawkins ND, Ward JL, et al. 13C-Flux Spectral Analysis of Host-Pathogen Metabolism Reveals a Mixed Diet for Intracellular Mycobacterium Tuberculosis. Chem Biol (2013) 20:1012–21. doi: 10.1016/j.chembiol.2013.06.012
16. Flores-Valdez MA, Aceves-Sánchez M de J, Peterson EJR, Baliga N, Bravo-Madrigal J, de la Cruz-Villegas MÁ, et al. Transcriptional Portrait of M. Bovis BCG During Biofilm Production Shows Genes Differentially Expressed During Intercellular Aggregation and Substrate Attachment. Sci Rep (2020) 10:12578. doi: 10.1038/s41598-020-69152-2
17. Aiewsakun P, Prombutara P, Siregar TAP, Laopanupong T, Kanjanasirirat P, Khumpanied T, et al. Transcriptional Response to the Host Cell Environment of a Multidrug-Resistant Mycobacterium Tuberculosis Clonal Outbreak Beijing Strain Reveals Its Pathogenic Features. Sci Rep (2021) 11(1):3199. doi: 10.1038/s41598-021-82905-x
18. Kavvas ES, Catoiu E, Mih N, Yurkovich JT, Seif Y, Dillon N, et al. Machine Learning and Structural Analysis of Mycobacterium Tuberculosis Pan-Genome Identifies Genetic Signatures of Antibiotic Resistance. Nat Commun (2018) 9(1):4306. doi: 10.1038/s41467-018-06634-y
19. Farhat MR, Freschi L, Calderon R, Ioerger T, Snyder M, Meehan CJ, et al. GWAS for Quantitative Resistance Phenotypes in Mycobacterium Tuberculosis Reveals Resistance Genes and Regulatory Regions. Nat Commun (2019) 10(1):2128. doi: 10.1038/s41467-019-10110-6
20. Coll F, Phelan J, Hill-Cawthorne GA, Nair MB, Mallard K, Ali S, et al. Genome-Wide Analysis of Multi- and Extensively Drug-Resistant Mycobacterium Tuberculosis. Nat Genet (2018) 50:307–16. doi: 10.1038/s41588-017-0029-0
21. Chiner-Oms Á, Berney M, Boinett C, González-Candelas F, Young DB, Gagneux S, et al. Genome-Wide Mutational Biases Fuel Transcriptional Diversity in the Mycobacterium Tuberculosis Complex. Nat Commun (2019) 10(1):3994. doi: 10.1038/s41467-019-11948-6
22. Gomez-Gonzalez PJ, Andreu N, Phelan JE, de Sessions PF, Glynn JR, Crampin AC, et al. An Integrated Whole Genome Analysis of Mycobacterium Tuberculosis Reveals Insights Into Relationship Between Its Genome, Transcriptome and Methylome. Sci Rep (2019) 9(1):5204. doi: 10.1038/s41598-019-41692-2
23. Zimmermann M, Kogadeeva M, Gengenbacher M, McEwen G, Mollenkopf H-J, Zamboni N, et al. Integration of Metabolomics and Transcriptomics Reveals a Complex Diet of Mycobacterium Tuberculosis During Early Macrophage Infection. mSystems (2017) 2(4):e00057–17. doi: 10.1128/msystems.00057-17
24. Hasenoehrl EJ, Rae Sajorda D, Berney-Meyer L, Johnson S, Tufariello JAM, Fuhrer T, et al. Derailing the Aspartate Pathway of Mycobacterium Tuberculosis to Eradicate Persistent Infection. Nat Commun (2019) 10(1):4215. doi: 10.1038/s41467-019-12224-3
25. Brzostek A, Płociński P, Minias A, Ciszewska A, Gąsior F, Pawełczyk J, et al. Dissecting the RecA-(In)dependent Response to Mitomycin C in Mycobacterium Tuberculosis Using Transcriptional Profiling and Proteomics Analyses. Cells (2021) 10(5):1168. doi: 10.3390/cells10051168
26. Chai Q, Lu Z, Liu Z, Zhong Y, Zhang F, Qiu C, et al. Lung Gene Expression Signatures Suggest Pathogenic Links and Molecular Markers for Pulmonary Tuberculosis, Adenocarcinoma and Sarcoidosis. Commun Biol (2020) 3:604. doi: 10.1038/s42003-020-01318-0
27. Kaufmann E, Sanz J, Dunn JL, Khan N, Mendonça LE, Pacis A, et al. BCG Educates Hematopoietic Stem Cells to Generate Protective Innate Immunity Against Tuberculosis. Cell (2018) 172:176–90.e19. doi: 10.1016/j.cell.2017.12.031
28. Altman MC, Rinchai D, Baldwin N, Toufiq M, Whalen E, Garand M, et al. Development of a Fixed Module Repertoire for the Analysis and Interpretation of Blood Transcriptome Data. Nat Commun (2021) 12(1):4385. doi: 10.1038/s41467-021-24584-w
29. Esterhuyse MM, Weiner J, Caron E, Loxton AG, Iannaccone M, Wagman C, et al. Epigenetics and Proteomics Join Transcriptomics in the Quest for Tuberculosis Biomarkers. mBio (2015) 6(5):e01187–15. doi: 10.1128/mBio.01187-15
30. Banaei-Esfahani A, Nicod C, Aebersold R, Collins BC. Systems Proteomics Approaches to Study Bacterial Pathogens: Application to Mycobacterium Tuberculosis. Curr Opin Microbiol (2017) 39:64–72. doi: 10.1016/j.mib.2017.09.013
31. Jungblut PR, Müller EC, Mattow J, Kaufmann SHE. Proteomics Reveals Open Reading Frames in Mycobacterium Tuberculosis H37Rv Not Predicted by Genomics. Infect Immun (2001) 69:5905–7. doi: 10.1128/IAI.69.9.5905-5907.2001
32. Young YB, Suk AK, Kim JS, Euh HJ, Bai GH, Cho SN, et al. Antigens Secreted From Mycobacterium Tuberculosis: Identification by Proteomics Approach and Test for Diagnostic Marker. Proteomics (2004) 4:3299–307. doi: 10.1002/pmic.200400980
33. Gillet LC, Navarro P, Tate S, Röst H, Selevsek N, Reiter L, et al. Targeted Data Extraction of the MS/MS Spectra Generated by Data-Independent Acquisition: A New Concept for Consistent and Accurate Proteome Analysis. Mol Cell Proteomics (2012) 11(6):O111.016717. doi: 10.1074/mcp.O111.016717
34. Rosenkrands I, King A, Weldingh K, Moniatte M, Moertz E, Andersen P. Towards the Proteome of Mycobacterium Tuberculosis. Electrophoresis 21(17):3740–56. doi: 10.1002/1522-2683(200011)21:17<3740::AID-ELPS3740>3.0
35. Jungblut PR, Schaible UE, Mollenkopf H-J, Zimny-Arndt U, Raupach B, Mattow J, et al. Comparative Proteome Analysis of Mycobacterium Tuberculosis and Mycobacterium Bovis BCG Strains: Towards Functional Genomics of Microbial Pathogens. Mol Microbiol 33(6):1103–17. doi: 10.1046/j.1365-2958.1999.01549.x
36. Betts JC, Dodson P, Quan s, Lewis AP, Thomas PJ, Duncan K, et al. Comparison of the Proteome of Mycobacterium Tuberculosis Strain H37Rv With Clinical Isolate CDC 1551. Microbiology (2000) 146(Pt 12):3205–16. doi: 10.1099/00221287-146-12-3205
37. Jhingan GD, Kumari S, Jamwal SV, Kalam H, Arora D, Jain N, et al. Comparative Proteomic Analyses of Avirulent, Virulent, and Clinical Strains of Mycobacterium Tuberculosis Identify Strain-Specific Patterns. J Biol Chem (2016) 291:14257–73. doi: 10.1074/jbc.M115.666123
38. Verma R, Pinto SM, Patil AH, Advani J, Subba P, Kumar M, et al. Quantitative Proteomic and Phosphoproteomic Analysis of H37Ra and H37Rv Strains of Mycobacterium Tuberculosis. J Proteome Res (2017) 16:1632–45. doi: 10.1021/acs.jproteome.6b00983
39. Saiboonjan B, Roytrakul S, Sangka A, Lulitanond V, Faksri K, Namwat W. Proteomic Analysis of Drug-Susceptible and Multidrug-Resistant Nonreplicating Beijing Strains of Mycobacterium Tuberculosis Cultured In Vitro. Biochem Biophys Rep (2021) 26:100960. doi: 10.1016/j.bbrep.2021.100960
40. Cho SH, Goodlett D, Franzblau S. ICAT-Based Comparative Proteomic Analysis of Non-Replicating Persistent Mycobacterium Tuberculosis. Tuberculosis (2006) 86:445–60. doi: 10.1016/j.tube.2005.10.002
41. Kruh NA, Troudt J, Izzo A, Prenni J, Dobos KM. Portrait of a Pathogen: The Mycobacterium Tuberculosis Proteome In Vivo. PLoS One (2010) 5(11):e13938. doi: 10.1371/journal.pone.0013938
42. Wolfe LM, Mahaffey SB, Kruh NA, Dobos KM. Proteomic Definition of the Cell Wall of Mycobacterium Tuberculosis. J Proteome Res (2010) 9:5816–26. doi: 10.1021/pr1005873
43. Li H, Wei S, Fang Y, Li M, Li X, Li Z, et al. Quantitative Proteomic Analysis of Host Responses Triggered by Mycobacterium Tuberculosis Infection in Human Macrophage Cells. Acta Biochim Biophys Sin (2017) 49:835–44. doi: 10.1093/abbs/gmx080
44. Covert BA, Spencer JS, Orme IM, Belisle JT. The Application of Proteomics in Defining the T Cell Antigens of Mycobacterium Tuberculosis. Proteomics (2001) 1(4):574–86. doi: 10.1002/1615-9861(200104)1:4<574::AID-PROT574>3.0.CO;2-8
45. Penn BH, Netter Z, Johnson JR, von Dollen J, Jang GM, Johnson T, et al. An Mtb-Human Protein-Protein Interaction Map Identifies a Switch Between Host Antiviral and Antibacterial Responses. Mol Cell (2018) 71:637–48.e5. doi: 10.1016/j.molcel.2018.07.010
46. Schubert OT, Mouritsen J, Ludwig C, Röst HL, Rosenberger G, Arthur PK, et al. The Mtb Proteome Library: A Resource of Assays to Quantify the Complete Proteome of Mycobacterium Tuberculosis. Cell Host Microbe (2013) 13:602–12. doi: 10.1016/j.chom.2013.04.008
47. Sun H, Pan L, Jia H, Zhang Z, Gao M, Huang M, et al. Label-Free Quantitative Proteomics Identifies Novel Plasma Biomarkers for Distinguishing Pulmonary Tuberculosis and Latent Infection. Front Microbiol (2018) 9:1267. doi: 10.3389/fmicb.2018.01267
48. Mateos J, Estévez O, González-Fernández Á, Anibarro L, Pallarés Á, Reljic R, et al. Serum Proteomics of Active Tuberculosis Patients and Contacts Reveals Unique Processes Activated During Mycobacterium Tuberculosis Infection. Sci Rep (2020) 10:3844. doi: 10.1038/s41598-020-60753-5
49. Garay-Baquero DJ, White CH, Walker NF, Tebruegge M, Schiff HF, Ugarte-Gil C, et al. Comprehensive Plasma Proteomic Profiling Reveals Biomarkers for Active Tuberculosis. JCI Insight (2020) 5(18):e137427. doi: 10.1172/JCI.INSIGHT.137427
50. Kunnath-Velayudhan S, Salamon H, Wang HY, Davidow AL, Molina DM, Huynh VT, et al. Dynamic Antibody Responses to the Mycobacterium Tuberculosis Proteome. Proc Natl Acad Sci USA (2010) 107:14703–8. doi: 10.1073/pnas.1009080107
51. Diray-Arce J, Conti MG, Petrova B, Kanarek N, Angelidou A, Levy O. Integrative Metabolomics to Identify Molecular Signatures of Responses to Vaccines and Infections. Metabolites (2020) 10:1–18. doi: 10.3390/metabo10120492
52. Mirsaeidi M, Banoei MM, Winston BW, Schraufnagel DE. Metabolomics: Applications and Promise in Mycobacterial Disease. Ann Am Thorac Soc (2015) 12:1278–87. doi: 10.1513/AnnalsATS.201505-279PS
53. Johnson CH, Ivanisevic J, Siuzdak G. Metabolomics: Beyond Biomarkers and Towards Mechanisms. Nat Rev Mol Cell Biol (2016) 17:451–9. doi: 10.1038/nrm.2016.25
54. Borah K, Rickman OJ, Voutsina N, Ampong I, Gao D, Baple EL, et al. A Quantitative LC-MS/MS Method for Analysis of Mitochondrial -Specific Oxysterol Metabolism. Redox Biol (2020) 36:101595. doi: 10.1016/j.redox.2020.101595
55. Taylor MJ, Lukowski JK, Anderton CR. Spatially Resolved Mass Spectrometry at the Single Cell: Recent Innovations in Proteomics and Metabolomics. J Am Soc Mass Spectrom (2021) 32:872–94. doi: 10.1021/jasms.0c00439
56. Larrouy-Maumus G, Biswas T, Hunt DM, Kelly G, Tsodikov OV, de Carvalho LPS. Discovery of a Glycerol 3-Phosphate Phosphatase Reveals Glycerophospholipid Polar Head Recycling in Mycobacterium Tuberculosis. Proc Natl Acad Sci USA (2013) 110:11320–5. doi: 10.1073/pnas.1221597110
57. Jansen RS, Mandyoli L, Hughes R, Wakabayashi S, Pinkham JT, Selbach B, et al. Aspartate Aminotransferase Rv3722c Governs Aspartate-Dependent Nitrogen Metabolism in Mycobacterium Tuberculosis. Nat Commun (2020) 11(1):1960. doi: 10.1038/s41467-020-15876-8
58. Øyås O, Borrell S, Trauner A, Zimmermann M, Feldmann J, Liphardt T, et al. Model-Based Integration of Genomics and Metabolomics Reveals SNP Functionality in Mycobacterium Tuberculosis. Proc Natl Acad Sci USA (2020) 117:(15):8494–502. doi: 10.1073/pnas.1915551117
59. Collins JM, Walker DI, Jones DP, Tukvadze N, Liu KH, Tran VT, et al. High-Resolution Plasma Metabolomics Analysis to Detect Mycobacterium Tuberculosisassociated Metabolites That Distinguish Active Pulmonary Tuberculosis in Humans. PLoS One (2018) 13(10):e0205398. doi: 10.1371/journal.pone.0205398
60. Collins JM, Siddiqa A, Jones DP, Liu K, Kempker RR, Nizam A, et al. Tryptophan Catabolism Reflects Disease Activity in Human Tuberculosis. JCI Insight (2020) 5(10):e137131. doi: 10.1172/JCI.INSIGHT.137131
61. Mahapatra S, Hess AM, Johnson JL, Eisenach KD, DeGroote MA, Gitta P, et al. A Metabolic Biosignature of Early Response to Anti-Tuberculosis Treatment. BMC Infect Dis (2014) 14:53. doi: 10.1186/1471-2334-14-53
62. Mackenzie JS, Lamprecht DA, Asmal R, Adamson JH, Borah K, Beste DJV, et al. Bedaquiline Reprograms Central Metabolism to Reveal Glycolytic Vulnerability in Mycobacterium Tuberculosis. Nat Commun (2020) 11(1):6092. doi: 10.1038/s41467-020-19959-4
63. Knoll KE, Lindeque Z, Adeniji AA, Oosthuizen CB, Lall N, Loots DT. Elucidating the Antimycobacterial Mechanism of Action of Decoquinate Derivative Rmb041 Using Metabolomics. Antibiotics (2021) 10(6):693. doi: 10.3390/antibiotics10060693
64. Rêgo AM, Alves da Silva D, Ferreira NV, de Pina LC, Evaristo JAM, Caprini Evaristo GP, et al. Metabolic Profiles of Multidrug Resistant and Extensively Drug Resistant Mycobacterium Tuberculosis Unveiled by Metabolomics. Tuberculosis (2021) 126. doi: 10.1016/j.tube.2020.102043
65. de Carvalho LPS, Zhao H, Dickinson CE, Arango NM, Lima CD, Fischer SM, et al. Activity-Based Metabolomic Profiling of Enzymatic Function: Identification of Rv1248c as a Mycobacterial 2-Hydroxy-3-Oxoadipate Synthase. Chem Biol (2010) 17:323–32. doi: 10.1016/j.chembiol.2010.03.009
66. Marrero J, Rhee KY, Schnappinger D, Pethe K, Ehrt S. Gluconeogenic Carbon Flow of Tricarboxylic Acid Cycle Intermediates Is Critical for Mycobacterium Tuberculosis to Establish and Maintain Infection. Proc Natl Acad Sci USA (2010) 107:9819–24. doi: 10.1073/pnas.1000715107
67. Black KA, Duan L, Mandyoli L, Selbach BP, Xu W, Ehrt S, et al. Metabolic Bifunctionality of Rv0812 Couples Folate and Peptidoglycan Biosynthesis in Mycobacterium Tuberculosis. J Exp Med (2021) 218(7):e20191957. doi: 10.1084/jem.20191957
68. de Carvalho LPS, Fischer SM, Marrero J, Nathan C, Ehrt S, Rhee KY. Metabolomics of Mycobacterium Tuberculosis Reveals Compartmentalized Co-Catabolism of Carbon Substrates. Chem Biol (2010) 17:1122–31. doi: 10.1016/j.chembiol.2010.08.009
69. Basu P, Sandhu N, Bhatt A, Singh A, Balhana R, Gobe I, et al. The Anaplerotic Node Is Essential for the Intracellular Survival of Mycobacterium Tuberculosis. J Biol Chem (2018) 293:5695–704. doi: 10.1074/jbc.RA118.001839
70. Borah K, Girardi K do C de V, Mendum TA, Lery LMS, Beste DJV, Lara FA, et al. Intracellular Mycobacterium Leprae Utilizes Host Glucose as a Carbon Source in Schwann Cells. mBio (2019) 10(6):e02351–19. doi: 10.1128/mBio.02351-19
71. Creek DJ, Mazet M, Achcar F, Anderson J, Kim DH, Kamour R, et al. Probing the Metabolic Network in Bloodstream-Form Trypanosoma Brucei Using Untargeted Metabolomics With Stable Isotope Labelled Glucose. PLoS Pathog (2015) 11:1–25. doi: 10.1371/journal.ppat.1004689
72. Häuslein I, Manske C, Goebel W, Eisenreich W, Hilbi H. Pathway Analysis Using 13C-Glycerol and Other Carbon Tracers Reveals a Bipartite Metabolism of Legionella Pneumophila. Mol Microbiol (2016) 100:229–46. doi: 10.1111/mmi.13313
73. Borah K, Mendum TA, Hawkins ND, Ward JL, Beale MH, Larrouy-Maumus G, et al. Metabolic Fluxes for Nutritional Flexibility of Mycobacterium Tuberculosis. Mol Syst Biol (2021) 17(5):e10280. doi: 10.15252/msb.202110280
74. Gouzy A, Larrouy-Maumus G, Bottai D, Levillain F, Dumas A, Wallach JB, et al. Mycobacterium Tuberculosis Exploits Asparagine to Assimilate Nitrogen and Resist Acid Stress During Infection. PLoS Pathog (2014) 10(2):e1003928. doi: 10.1371/journal.ppat.1003928
75. Agapova A, Serafini A, Petridis M, Hunt DM, Garza-Garcia A, Sohaskey CD, et al. Flexible Nitrogen Utilisation by the Metabolic Generalist Pathogen Mycobacterium Tuberculosis. Elife (2019) 8:e41129. doi: 10.7554/eLife.41129
76. Rienksma RA, Schaap PJ, dos Santos VAPM, Suarez-Diez M. Modeling the Metabolic State of Mycobacterium Tuberculosis Upon Infection. Front Cell Infect Microbiol (2018) 8:264. doi: 10.3389/fcimb.2018.00264
77. Rienksma RA, Suarez-Diez M, Spina L, Schaap PJ, Martins Dos Santos VAP. Systems-Level Modeling of Mycobacterial Metabolism for the Identification of New (Multi-)Drug Targets. Semin Immunol (2014) 26:610–22. doi: 10.1016/j.smim.2014.09.013
78. Beste DJV, Hooper T, Stewart G, Bonde B, Avignone-Rossa C, Bushell ME, et al. GSMN-TB: A Web-Based Genome-Scale Network Model of Mycobacterium Tuberculosis Metabolism. Genome Biol (2007) 8(5):R89. doi: 10.1186/gb-2007-8-5-r89
79. Bonde BK, Beste DJ, Laing E, Kierzek AM, McFadden J. Differential Producibility Analysis (DPA) of Transcriptomic Data With Metabolic Networks: Deconstructing the Metabolic Response of M. tuberculosis. PLoS Comput Biol (2011) 7(6):e1002060. doi: 10.1371/journal.pcbi.1002060
80. Borah K, Kearney JL, Banerjee R, Vats P, Wu H, Dahale S, et al. GSMN-ML-A Genome Scale Metabolic Network Reconstruction of the Obligate Human Pathogen Mycobacterium Leprae. PLoS Negl Trop Dis (2020) 14:1–20. doi: 10.1371/journal.pntd.0007871
81. Lofthouse EK, Wheeler PR, Beste DJV, Khatri BL, Wu H, Mendum TA, et al. Systems-Based Approaches to Probing Metabolic Variation Within the Mycobacterium Tuberculosis Complex. PLoS One (2013) 8(9):e75913. doi: 10.1371/journal.pone.0075913
82. Kavvas ES, Seif Y, Yurkovich JT, Norsigian C, Poudel S, Greenwald WW, et al. Updated and Standardized Genome-Scale Reconstruction of Mycobacterium Tuberculosis H37Rv, Iek1011, Simulates Flux States Indicative of Physiological Conditions. BMC Syst Biol (2018) 12(1):25. doi: 10.1186/s12918-018-0557-y
83. López-Agudelo VA, Mendum TA, Laing E, Wu HH, Baena A, Barrera LF, et al. A Systematic Evaluation of Mycobacterium Tuberculosis Genome-Scale Metabolic Networks. PLoS Comput Biol (2020) 16(6):e1007533. doi: 10.1371/journal.pcbi.1007533
84. Wiechert W, Moo M, Ren Petersen S, de Graaf AA. A Universal Framework for 13 C Metabolic Flux Analysis (2001). Available at: www.idealibrary.com.
85. Graf M, Haas T, Teleki A, Feith A, Cerff M, Wiechert W, et al. Revisiting the Growth Modulon of Corynebacterium Glutamicum Under Glucose Limited Chemostat Conditions. Front Bioeng Biotechnol (2020) 8:584614. doi: 10.3389/fbioe.2020.584614
86. Beste DJV, Bonde B, Hawkins N, Ward JL, Beale MH, Noack S, et al. 13c Metabolic Flux Analysis Identifies an Unusual Route for Pyruvate Dissimilation in Mycobacteria Which Requires Isocitrate Lyase and Carbon Dioxide Fixation. PLoS Pathog (2011) 7(7):e1002091. doi: 10.1371/journal.ppat.1002091
87. Antoniewicz MR. A Guide to 13C Metabolic Flux Analysis for the Cancer Biologist. Exp Mol Med (2018) 50(4):1–13. doi: 10.1038/s12276-018-0060-y
88. Kruger NJ, Ratcliffe RG. Insights Into Plant Metabolic Networks From Steady-State Metabolic Flux Analysis. Biochimie (2009) 91:697–702. doi: 10.1016/j.biochi.2009.01.004
89. He L, Wu SG, Zhang M, Chen Y, Tang YJ. WUFlux: An Open-Source Platform for 13C Metabolic Flux Analysis of Bacterial Metabolism. BMC Bioinformatics (2016) 17:444. doi: 10.1186/s12859-016-1314-0
90. Stein CM, Sausville L, Wejse C, Sobota RS, Zetola NM, Hill PC, et al. Genomics of Human Pulmonary Tuberculosis: From Genes to Pathways. Curr Genet Med Rep (2017) 5:149–66. doi: 10.1007/s40142-017-0130-9
91. Burel JG, Babor M, Pomaznoy M, Lindestam Arlehamn CS, Khan N, Sette A, et al. Host Transcriptomics as a Tool to Identify Diagnostic and Mechanistic Immune Signatures of Tuberculosis. Front Immunol (2019) 10:221. doi: 10.3389/fimmu.2019.00221
92. Goff A, Cantillon D, Wildner LM, Waddell SJ. Multi-Omics Technologies Applied to Tuberculosis Drug Discovery. Appl Sci (Switzerland) (2020) 10(13):4629. doi: 10.3390/app10134629
Keywords: Mycobacterium tuberculosis, tuberculosis, systems biology, macrophage, omic technology
Citation: Borah K, Xu Y and McFadden J (2021) Dissecting Host-Pathogen Interactions in TB Using Systems-Based Omic Approaches. Front. Immunol. 12:762315. doi: 10.3389/fimmu.2021.762315
Received: 21 August 2021; Accepted: 18 October 2021;
Published: 02 November 2021.
Edited by:
Sangeeta Tiwari, The University of Texas at El Paso, United StatesReviewed by:
Shanmuga Priyaa Madhukaran, University of Texas Southwestern Medical Center, United StatesQiyao Chai, Chinese Academy of Sciences (CAS), China
Copyright © 2021 Borah, Xu and McFadden. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Khushboo Borah, k.borah@surrey.ac.uk; Johnjoe McFadden, j.mcfadden@surrey.ac.uk