TGx-DDI, a Transcriptomic Biomarker for Genotoxicity Hazard Assessment of Pharmaceuticals and Environmental Chemicals.

Genotoxicity testing is an essential component of the safety assessment paradigm required by regulatory agencies world-wide for analysis of drug candidates, and environmental and industrial chemicals. Current genotoxicity testing batteries feature a high incidence of irrelevant positive findings-particularly for in vitro chromosomal damage (CD) assays. The risk management of compounds with positive in vitro findings is a major challenge and requires complex, time consuming, and costly follow-up strategies including animal testing. Thus, regulators are urgently in need of new testing approaches to meet legislated mandates. Using machine learning, we identified a set of transcripts that responds predictably to DNA-damage in human cells that we refer to as the TGx-DDI biomarker, which was originally referred to as TGx-28.65. We proposed to use this biomarker in conjunction with current genotoxicity testing batteries to differentiate compounds with irrelevant "false" positive findings in the in vitro CD assays from true DNA damaging agents (i.e., for de-risking agents that are clastogenic in vitro but not in vivo). We validated the performance of the TGx-DDI biomarker to identify true DNA damaging agents, assessed intra- and inter- laboratory reproducibility, and cross-platform performance. Recently, to augment the application of this biomarker, we developed a high-throughput cell-based genotoxicity testing system using the NanoString nCounter® technology. Here, we review the status of TGx-DDI development, its integration in the genotoxicity testing paradigm, and progress to date in its qualification at the US Food and Drug Administration (FDA) as a drug development tool. If successfully validated and implemented, the TGx-DDI biomarker assay is expected to significantly augment the current strategy for the assessment of genotoxic hazards for drugs and chemicals.

Development and validation of a high-throughput transcriptomic biomarker to address 21st century genetic toxicology needs.

Interpretation of positive genotoxicity findings using the current in vitro testing battery is a major challenge to industry and regulatory agencies. These tests, especially mammalian cell assays, have high sensitivity but suffer from low specificity, leading to high rates of irrelevant positive findings (i.e., positive results in vitro that are not relevant to human cancer hazard). We developed an in vitro transcriptomic biomarker-based approach that provides biological relevance to positive genotoxicity assay data, particularly for in vitro chromosome damage assays, and propose its application for assessing the relevance of the in vitro positive results to carcinogenic hazard. The transcriptomic biomarker TGx-DDI (previously known as TGx-28.65) readily distinguishes DNA damage-inducing (DDI) agents from non-DDI agents. In this study, we demonstrated the ability of the biomarker to classify 45 test agents across a broad set of chemical classes as DDI or non-DDI. Furthermore, we assessed the biomarker's utility in derisking known irrelevant positive agents and evaluated its performance across analytical platforms. We correctly classified 90% (9 of 10) of chemicals with irrelevant positive findings in in vitro chromosome damage assays as negative. We developed a standardized experimental and analytical protocol for our transcriptomics biomarker, as well as an enhanced application of TGx-DDI for high-throughput cell-based genotoxicity testing using nCounter technology. This biomarker can be integrated in genetic hazard assessment as a follow-up to positive chromosome damage findings. In addition, we propose how it might be used in chemical screening and assessment. This approach offers an opportunity to significantly improve risk assessment and reduce cost.

A predictive toxicogenomics signature to classify genotoxic versus non-genotoxic chemicals in human TK6 cells.

Genotoxicity testing is a critical component of chemical assessment. The use of integrated approaches in genetic toxicology, including the incorporation of gene expression data to determine the DNA damage response pathways involved in response, is becoming more common. In companion papers previously published in Environmental and Molecular Mutagenesis, Li et al. (2015) [6] developed a dose optimization protocol that was based on evaluating expression changes in several well-characterized stress-response genes using quantitative real-time PCR in human lymphoblastoid TK6 cells in culture. This optimization approach was applied to the analysis of TK6 cells exposed to one of 14 genotoxic or 14 non-genotoxic agents, with sampling 4 h post-exposure. Microarray-based transcriptomic analyses were then used to develop a classifier for genotoxicity using the nearest shrunken centroids method. A panel of 65 genes was identified that could accurately classify toxicants as genotoxic or non-genotoxic. In Buick et al. (2015) [1], the utility of the biomarker for chemicals that require metabolic activation was evaluated. In this study, TK6 cells were exposed to increasing doses of four chemicals (two genotoxic that require metabolic activation and two non-genotoxic chemicals) in the presence of rat liver S9 to demonstrate that S9 does not impair the ability to classify genotoxicity using this genomic biomarker in TK6cells.

Integration of metabolic activation with a predictive toxicogenomics signature to classify genotoxic versus nongenotoxic chemicals in human TK6 cells.

The use of integrated approaches in genetic toxicology, including the incorporation of gene expression data to determine the molecular pathways involved in the response, is becoming more common. In a companion article, a genomic biomarker was developed in human TK6 cells to classify chemicals as genotoxic or nongenotoxic. Because TK6 cells are not metabolically competent, we set out to broaden the utility of the biomarker for use with chemicals requiring metabolic activation. Specifically, chemical exposures were conducted in the presence of rat liver S9. The ability of the biomarker to classify genotoxic (benzo[a]pyrene, BaP; aflatoxin B1, AFB1) and nongenotoxic (dexamethasone, DEX; phenobarbital, PB) agents correctly was evaluated. Cells were exposed to increasing chemical concentrations for 4 hr and collected 0 hr, 4 hr, and 20 hr postexposure. Relative survival, apoptosis, and micronucleus frequency were measured at 24 hr. Transcriptome profiles were measured with Agilent microarrays. Statistical modeling and bioinformatics tools were applied to classify each chemical using the genomic biomarker. BaP and AFB1 were correctly classified as genotoxic at the mid- and high concentrations at all three time points, whereas DEX was correctly classified as nongenotoxic at all concentrations and time points. The high concentration of PB was misclassified at 24 hr, suggesting that cytotoxicity at later time points may cause misclassification. The data suggest that the use of S9 does not impair the ability of the biomarker to classify genotoxicity in TK6 cells. Finally, we demonstrate that the biomarker is also able to accurately classify genotoxicity using a publicly available dataset derived from human HepaRG cells.

Development of a toxicogenomics signature for genotoxicity using a dose-optimization and informatics strategy in human cells.

The development of in vitro molecular biomarkers to accurately predict toxicological effects has become a priority to advance testing strategies for human health risk assessment. The application of in vitro transcriptomic biomarkers promises increased throughput as well as a reduction in animal use. However, the existing protocols for predictive transcriptional signatures do not establish appropriate guidelines for dose selection or account for the fact that toxic agents may have pleiotropic effects. Therefore, comparison of transcriptome profiles across agents and studies has been difficult. Here we present a dataset of transcriptional profiles for TK6 cells exposed to a battery of well-characterized genotoxic and nongenotoxic chemicals. The experimental conditions applied a new dose optimization protocol that was based on evaluating expression changes in several well-characterized stress-response genes using quantitative real-time PCR in preliminary dose-finding studies. The subsequent microarray-based transcriptomic analyses at the optimized dose revealed responses to the test chemicals that were typically complex, often exhibiting substantial overlap in the transcriptional responses between a variety of the agents making analysis challenging. Using the nearest shrunken centroids method we identified a panel of 65 genes that could accurately classify toxicants as genotoxic or nongenotoxic. To validate the 65-gene panel as a genomic biomarker of genotoxicity, the gene expression profiles of an additional three well-characterized model agents were analyzed and a case study demonstrating the practical application of this genomic biomarker-based approach in risk assessment was performed to demonstrate its utility in genotoxicity risk assessment.

Comparison of toxicogenomics and traditional approaches to inform mode of action and points of departure in human health risk assessment of benzo[a]pyrene in drinking water.

Toxicogenomics is proposed to be a useful tool in human health risk assessment. However, a systematic comparison of traditional risk assessment approaches with those applying toxicogenomics has never been done. We conducted a case study to evaluate the utility of toxicogenomics in the risk assessment of benzo[a]pyrene (BaP), a well-studied carcinogen, for drinking water exposures. Our study was intended to compare methodologies, not to evaluate drinking water safety. We compared traditional (RA1), genomics-informed (RA2) and genomics-only (RA3) approaches. RA2 and RA3 applied toxicogenomics data from human cell cultures and mice exposed to BaP to determine if these data could provide insight into BaP's mode of action (MOA) and derive tissue-specific points of departure (POD). Our global gene expression analysis supported that BaP is genotoxic in mice and allowed the development of a detailed MOA. Toxicogenomics analysis in human lymphoblastoid TK6 cells demonstrated a high degree of consistency in perturbed pathways with animal tissues. Quantitatively, the PODs for traditional and transcriptional approaches were similar (liver 1.2 vs. 1.0 mg/kg-bw/day; lungs 0.8 vs. 3.7 mg/kg-bw/day; forestomach 0.5 vs. 7.4 mg/kg-bw/day). RA3, which applied toxicogenomics in the absence of apical toxicology data, demonstrates that this approach provides useful information in data-poor situations. Overall, our study supports the use of toxicogenomics as a relatively fast and cost-effective tool for hazard identification, preliminary evaluation of potential carcinogens, and carcinogenic potency, in addition to identifying current limitations and practical questions for future work.

A systems approach to predict oncometabolites via context-specific genome-scale metabolic networks.

Altered metabolism in cancer cells has been viewed as a passive response required for a malignant transformation. However, this view has changed through the recently described metabolic oncogenic factors: mutated isocitrate dehydrogenases (IDH), succinate dehydrogenase (SDH), and fumarate hydratase (FH) that produce oncometabolites that competitively inhibit epigenetic regulation. In this study, we demonstrate in silico predictions of oncometabolites that have the potential to dysregulate epigenetic controls in nine types of cancer by incorporating massive scale genetic mutation information (collected from more than 1,700 cancer genomes), expression profiling data, and deploying Recon 2 to reconstruct context-specific genome-scale metabolic models. Our analysis predicted 15 compounds and 24 substructures of potential oncometabolites that could result from the loss-of-function and gain-of-function mutations of metabolic enzymes, respectively. These results suggest a substantial potential for discovering unidentified oncometabolites in various forms of cancers.

Genome-scale models of metabolism and gene expression extend and refine growth phenotype prediction.

Growth is a fundamental process of life. Growth requirements are well-characterized experimentally for many microbes; however, we lack a unified model for cellular growth. Such a model must be predictive of events at the molecular scale and capable of explaining the high-level behavior of the cell as a whole. Here, we construct an ME-Model for Escherichia coli--a genome-scale model that seamlessly integrates metabolic and gene product expression pathways. The model computes ~80% of the functional proteome (by mass), which is used by the cell to support growth under a given condition. Metabolism and gene expression are interdependent processes that affect and constrain each other. We formalize these constraints and apply the principle of growth optimization to enable the accurate prediction of multi-scale phenotypes, ranging from coarse-grained (growth rate, nutrient uptake, by-product secretion) to fine-grained (metabolic fluxes, gene expression levels). Our results unify many existing principles developed to describe bacterial growth.

COBRApy: COnstraints-Based Reconstruction and Analysis for Python.

BACKGROUND: COnstraint-Based Reconstruction and Analysis (COBRA) methods are widely used for genome-scale modeling of metabolic networks in both prokaryotes and eukaryotes. Due to the successes with metabolism, there is an increasing effort to apply COBRA methods to reconstruct and analyze integrated models of cellular processes. The COBRA Toolbox for MATLAB is a leading software package for genome-scale analysis of metabolism; however, it was not designed to elegantly capture the complexity inherent in integrated biological networks and lacks an integration framework for the multiomics data used in systems biology. The openCOBRA Project is a community effort to promote constraints-based research through the distribution of freely available software. RESULTS: Here, we describe COBRA for Python (COBRApy), a Python package that provides support for basic COBRA methods. COBRApy is designed in an object-oriented fashion that facilitates the representation of the complex biological processes of metabolism and gene expression. COBRApy does not require MATLAB to function; however, it includes an interface to the COBRA Toolbox for MATLAB to facilitate use of legacy codes. For improved performance, COBRApy includes parallel processing support for computationally intensive processes. CONCLUSION: COBRApy is an object-oriented framework designed to meet the computational challenges associated with the next generation of stoichiometric constraint-based models and high-density omics data sets. AVAILABILITY: http://opencobra.sourceforge.net/

GIM3E: condition-specific models of cellular metabolism developed from metabolomics and expression data.

MOTIVATION: Genome-scale metabolic models have been used extensively to investigate alterations in cellular metabolism. The accuracy of these models to represent cellular metabolism in specific conditions has been improved by constraining the model with omics data sources. However, few practical methods for integrating metabolomics data with other omics data sources into genome-scale models of metabolism have been developed. RESULTS: GIM(3)E (Gene Inactivation Moderated by Metabolism, Metabolomics and Expression) is an algorithm that enables the development of condition-specific models based on an objective function, transcriptomics and cellular metabolomics data. GIM(3)E establishes metabolite use requirements with metabolomics data, uses model-paired transcriptomics data to find experimentally supported solutions and provides calculations of the turnover (production/consumption) flux of metabolites. GIM(3)E was used to investigate the effects of integrating additional omics datasets to create increasingly constrained solution spaces of Salmonella Typhimurium metabolism during growth in both rich and virulence media. This integration proved to be informative and resulted in a requirement of additional active reactions (12 in each case) or metabolites (26 or 29, respectively). The addition of constraints from transcriptomics also impacted the allowed solution space, and the cellular metabolites with turnover fluxes that were necessarily altered by the change in conditions increased from 118 to 271 of 1397. AVAILABILITY: GIM(3)E has been implemented in Python and requires a COBRApy 0.2.x. The algorithm and sample data described here are freely available at: http://opencobra.sourceforge.net/ CONTACTS: brianjamesschmidt@gmail.com

Identifying radiation exposure biomarkers from mouse blood transcriptome.

Ionising radiation is a pleiotropic stress agent that may induce a variety of adverse effects. Molecular biomarker approaches possess promise to assess radiation exposure, however, the pleiotropic nature of ionising radiation induced transcriptional responses and the historically poor inter-laboratory performance of omics-derived biomarkers serve as barriers to identification of unequivocal biomarker sets. Here, we present a whole-genome survey of the murine transcriptomic response to physiologically relevant radiation doses, 2 Gy and 8 Gy. We used this dataset with the Random Forest algorithm to correctly classify independently generated data and to identify putative metabolite biomarkers for radiation exposure.

Salmonella modulates metabolism during growth under conditions that induce expression of virulence genes.

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a facultative pathogen that uses complex mechanisms to invade and proliferate within mammalian host cells. To investigate possible contributions of metabolic processes to virulence in S. Typhimurium grown under conditions known to induce expression of virulence genes, we used a metabolomics-driven systems biology approach coupled with genome-scale modeling. First, we identified distinct metabolite profiles associated with bacteria grown in either rich or virulence-inducing media and report the most comprehensive coverage of the S. Typhimurium metabolome to date. Second, we applied an omics-informed genome-scale modeling analysis of the functional consequences of adaptive alterations in S. Typhimurium metabolism during growth under our conditions. Modeling efforts highlighted a decreased cellular capability to both produce and utilize intracellular amino acids during stationary phase culture in virulence conditions, despite significant abundance increases for these molecules as observed by our metabolomics measurements. Furthermore, analyses of omics data in the context of the metabolic model indicated rewiring of the metabolic network to support pathways associated with virulence. For example, cellular concentrations of polyamines were perturbed, as well as the predicted capacity for secretion and uptake.

Analysis of omics data with genome-scale models of metabolism.

Over the past decade a massive amount of research has been dedicated to generating omics data to gain insight into a variety of biological phenomena, including cancer, obesity, biofuel production, and infection. Although most of these omics data are available publicly, there is a growing concern that much of these data sit in databases without being used or fully analyzed. Statistical inference methods have been widely applied to gain insight into which genes may influence the activities of others in a given omics data set, however, they do not provide information on the underlying mechanisms or whether the interactions are direct or distal. Biochemically, genetically, and genomically consistent knowledge bases are increasingly being used to extract deeper biological knowledge and understanding from these data sets than possible by inferential methods. This improvement is largely due to knowledge bases providing a validated biological context for interpreting the data.

Exposure to ionizing radiation induced persistent gene expression changes in mouse mammary gland.

BACKGROUND: Breast tissue is among the most sensitive tissues to the carcinogenic actions of ionizing radiation and epidemiological studies have linked radiation exposure to breast cancer. Currently, molecular understanding of radiation carcinogenesis in mammary gland is hindered due to the scarcity of in vivo long-term follow up data. We undertook this study to delineate radiation-induced persistent alterations in gene expression in mouse mammary glands 2-month after radiation exposure. METHODS: Six to eight week old female C57BL/6J mice were exposed to 2 Gy of whole body γ radiation and mammary glands were surgically removed 2-month after radiation. RNA was isolated and microarray hybridization performed for gene expression analysis. Ingenuity Pathway Analysis (IPA) was used for biological interpretation of microarray data. Real time quantitative PCR was performed on selected genes to confirm the microarray data. RESULTS: Compared to untreated controls, the mRNA levels of a total of 737 genes were significantly (p<0.05) perturbed above 2-fold of control. More genes (493 genes; 67%) were upregulated than the number of downregulated genes (244 genes; 33%). Functional analysis of the upregulated genes mapped to cell proliferation and cancer related canonical pathways such as 'ERK/MAPK signaling', 'CDK5 signaling', and '14-3-3-mediated signaling'. We also observed upregulation of breast cancer related canonical pathways such as 'breast cancer regulation by Stathmin1', and 'HER-2 signaling in breast cancer' in IPA. Interestingly, the downregulated genes mapped to fewer canonical pathways involved in cell proliferation. We also observed that a number of genes with tumor suppressor function (GPRC5A, ELF1, NAB2, Sema4D, ACPP, MAP2, RUNX1) persistently remained downregulated in response to radiation exposure. Results from qRT-PCR on five selected differentially expressed genes confirmed microarray data. The PCR data on PPP4c, ELF1, MAPK12, PLCG1, and E2F6 showed similar trend in up and downregulation as has been observed with the microarray. CONCLUSIONS: Exposure to a clinically relevant radiation dose led to long-term activation of mammary gland genes involved in proliferative and metabolic pathways, which are known to have roles in carcinogenesis. When considered along with downregulation of a number of tumor suppressor genes, our study has implications for breast cancer initiation and progression after therapeutic radiation exposure.

In silico method for modelling metabolism and gene product expression at genome scale.

Transcription and translation use raw materials and energy generated metabolically to create the macromolecular machinery responsible for all cellular functions, including metabolism. A biochemically accurate model of molecular biology and metabolism will facilitate comprehensive and quantitative computations of an organism's molecular constitution as a function of genetic and environmental parameters. Here we formulate a model of metabolism and macromolecular expression. Prototyping it using the simple microorganism Thermotoga maritima, we show our model accurately simulates variations in cellular composition and gene expression. Moreover, through in silico comparative transcriptomics, the model allows the discovery of new regulons and improving the genome and transcription unit annotations. Our method presents a framework for investigating molecular biology and cellular physiology in silico and may allow quantitative interpretation of multi-omics data sets in the context of an integrated biochemical description of an organism.

Model-driven multi-omic data analysis elucidates metabolic immunomodulators of macrophage activation.

Macrophages are central players in immune response, manifesting divergent phenotypes to control inflammation and innate immunity through release of cytokines and other signaling factors. Recently, the focus on metabolism has been reemphasized as critical signaling and regulatory pathways of human pathophysiology, ranging from cancer to aging, often converge on metabolic responses. Here, we used genome-scale modeling and multi-omics (transcriptomics, proteomics, and metabolomics) analysis to assess metabolic features that are critical for macrophage activation. We constructed a genome-scale metabolic network for the RAW 264.7 cell line to determine metabolic modulators of activation. Metabolites well-known to be associated with immunoactivation (glucose and arginine) and immunosuppression (tryptophan and vitamin D3) were among the most critical effectors. Intracellular metabolic mechanisms were assessed, identifying a suppressive role for de-novo nucleotide synthesis. Finally, underlying metabolic mechanisms of macrophage activation are identified by analyzing multi-omic data obtained from LPS-stimulated RAW cells in the context of our flux-based predictions. Our study demonstrates metabolism's role in regulating activation may be greater than previously anticipated and elucidates underlying connections between activation and metabolic effectors.

Comparison of mouse urinary metabolic profiles after exposure to the inflammatory stressors γ radiation and lipopolysaccharide.

Metabolomics on easily accessible biofluids has the potential to provide rapid identification and distinction between stressors and inflammatory states. In the event of a radiological event, individuals with underlying medical conditions could present with similar symptoms to radiation poisoning, prominently nausea, diarrhea, vomiting and fever. Metabolomics of radiation exposure in mice has provided valuable biomarkers, and in this study we aimed to identify biomarkers of lipopolysaccharide (LPS) exposure to compare and contrast with ionizing radiation. LPS treatment leads to a severe inflammatory response and a cytokine storm, events similar to radiation exposure, and LPS exposure can recapitulate many of the responses seen in sepsis. Urine from control mice, LPS-treated mice, and mice irradiated with 3, 8 and 15 Gy of γ rays was analyzed by LCMS, and markers were extracted using SIMCA-P(+) and Random Forests. Markers were validated through tandem mass spectrometry against pure chemicals. Five metabolites, cytosine, cortisol, adenine, O-propanoylcarnitine and isethionic acid, showed increased excretion at 24 h after LPS treatment (P < 0.0001, 0.0393, 0.0393, <0.0001 and 0.0004, respectively). Of these, cytosine, adenine and O-propanoylcarnitine showed specificity to LPS treatment when compared to radiation. On the other hand, increased excretion of cortisol after LPS and radiation treatments indicated a rapid systemic response to inflammatory agents. Isethionic acid excretion, however, showed elevated levels not only after LPS treatment but also after a very high dose of radiation (15 Gy), while additional metabolites showed responsiveness to radiation but not LPS. Metabolomics therefore has the potential to distinguish between different inflammatory responses based on differential ion signatures. It can also provide quick and reliable assessment of medical conditions in a mass casualty radiological scenario and aid in effective triaging.

An experimentally-supported genome-scale metabolic network reconstruction for Yersinia pestis CO92.

BACKGROUND: Yersinia pestis is a gram-negative bacterium that causes plague, a disease linked historically to the Black Death in Europe during the Middle Ages and to several outbreaks during the modern era. Metabolism in Y. pestis displays remarkable flexibility and robustness, allowing the bacterium to proliferate in both warm-blooded mammalian hosts and cold-blooded insect vectors such as fleas. RESULTS: Here we report a genome-scale reconstruction and mathematical model of metabolism for Y. pestis CO92 and supporting experimental growth and metabolite measurements. The model contains 815 genes, 678 proteins, 963 unique metabolites and 1678 reactions, accurately simulates growth on a range of carbon sources both qualitatively and quantitatively, and identifies gaps in several key biosynthetic pathways and suggests how those gaps might be filled. Furthermore, our model presents hypotheses to explain certain known nutritional requirements characteristic of this strain. CONCLUSIONS: Y. pestis continues to be a dangerous threat to human health during modern times. The Y. pestis genome-scale metabolic reconstruction presented here, which has been benchmarked against experimental data and correctly reproduces known phenotypes, provides an in silico platform with which to investigate the metabolism of this important human pathogen.

Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0.

Over the past decade, a growing community of researchers has emerged around the use of constraint-based reconstruction and analysis (COBRA) methods to simulate, analyze and predict a variety of metabolic phenotypes using genome-scale models. The COBRA Toolbox, a MATLAB package for implementing COBRA methods, was presented earlier. Here we present a substantial update of this in silico toolbox. Version 2.0 of the COBRA Toolbox expands the scope of computations by including in silico analysis methods developed since its original release. New functions include (i) network gap filling, (ii) (13)C analysis, (iii) metabolic engineering, (iv) omics-guided analysis and (v) visualization. As with the first version, the COBRA Toolbox reads and writes systems biology markup language-formatted models. In version 2.0, we improved performance, usability and the level of documentation. A suite of test scripts can now be used to learn the core functionality of the toolbox and validate results. This toolbox lowers the barrier of entry to use powerful COBRA methods.

Technologies and approaches to elucidate and model the virulence program of salmonella.

Salmonella is a primary cause of enteric diseases in a variety of animals. During its evolution into a pathogenic bacterium, Salmonella acquired an elaborate regulatory network that responds to multiple environmental stimuli within host animals and integrates them resulting in fine regulation of the virulence program. The coordinated action by this regulatory network involves numerous virulence regulators, necessitating genome-wide profiling analysis to assess and combine efforts from multiple regulons. In this review we discuss recent high-throughput analytic approaches used to understand the regulatory network of Salmonella that controls virulence processes. Application of high-throughput analyses have generated large amounts of data and necessitated the development of computational approaches for data integration. Therefore, we also cover computer-aided network analyses to infer regulatory networks, and demonstrate how genome-scale data can be used to construct regulatory and metabolic systems models of Salmonella pathogenesis. Genes that are coordinately controlled by multiple virulence regulators under infectious conditions are more likely to be important for pathogenesis. Thus, reconstructing the global regulatory network during infection or, at the very least, under conditions that mimic the host cellular environment not only provides a bird's eye view of Salmonella survival strategy in response to hostile host environments but also serves as an efficient means to identify novel virulence factors that are essential for Salmonella to accomplish systemic infection in the host.