3D Hydrogel Culture System Recapitulates Key Tuberculosis Phenotypes and Demonstrates Pyrazinamide Efficacy

The mortality caused by tuberculosis (TB) infections is a global concern, and there is a need to improve understanding of the disease. Current in vitro infection models to study the disease have limitations such as short investigation durations and divergent transcriptional signatures. This study aims to overcome these limitations by developing a 3D collagen culture system that mimics the biomechanical and extracellular matrix (ECM) of lung microenvironment (collagen fibers, stiffness comparable to in vivo conditions) as the infection primarily manifests in the lungs. The system incorporates Mycobacterium tuberculosis (Mtb) infected human THP-1 or primary monocytes/macrophages. Dual RNA sequencing reveals higher mammalian gene expression similarity with patient samples than 2D macrophage infections. Similarly, bacterial gene expression more accurately recapitulates in vivo gene expression patterns compared to bacteria in 2D infection models. Key phenotypes observed in humans, such as foamy macrophages and mycobacterial cords, are reproduced in the model. This biomaterial system overcomes challenges associated with traditional platforms by modulating immune cells and closely mimicking in vivo infection conditions, including showing efficacy with clinically relevant concentrations of anti-TB drug pyrazinamide, not seen in any other in vitro infection model, making it reliable and readily adoptable for tuberculosis studies and drug screening.

Contribution of Microorganisms with the Clade II Nitrous Oxide Reductase to Suppression of Surface Emissions of Nitrous Oxide

The sources and sinks of nitrous oxide, as control emissions to the atmosphere, are generally poorly constrained for most environmental systems. Initial depth-resolved analysis of nitrous oxide flux from observation wells and the proximal surface within a nitrate contaminated aquifer system revealed high subsurface production but little escape from the surface. To better understand the environmental controls of production and emission at this site, we used a combination of isotopic, geochemical, and molecular analyses to show that chemodenitrification and bacterial denitrification are major sources of nitrous oxide in this subsurface, where low DO, low pH, and high nitrate are correlated with significant nitrous oxide production. Depth-resolved metagenomes showed that consumption of nitrous oxide near the surface was correlated with an enrichment of Clade II nitrous oxide reducers, consistent with a growing appreciation of their importance in controlling release of nitrous oxide to the atmosphere. Our work also provides evidence for the reduction of nitrous oxide at a pH of 4, well below the generally accepted limit of pH 5.

An atlas of causal and mechanistic drivers of interpatient heterogeneity in glioma

Grade IV glioma, formerly known as glioblastoma multiforme (GBM) is the most aggressive and lethal type of brain tumor, and its treatment remains challenging in part due to extensive interpatient heterogeneity in disease driving mechanisms and lack of prognostic and predictive biomarkers. Using mechanistic inference of node-edge relationship (MINER), we have analyzed multiomics profiles from 516 patients and constructed an atlas of causal and mechanistic drivers of interpatient heterogeneity in GBM (gbmMINER). The atlas has delineated how 30 driver mutations act in a combinatorial scheme to causally influence a network of regulators (306 transcription factors and 73 miRNAs) of 179 transcriptional "programs", influencing disease progression in patients across 23 disease states. Through extensive testing on independent patient cohorts, we share evidence that a machine learning model trained on activity profiles of programs within gbmMINER significantly augments risk stratification, identifying patients who are super-responders to standard of care and those that would benefit from 2 nd line treatments. In addition to providing mechanistic hypotheses regarding disease prognosis, the activity of programs containing targets of 2 nd line treatments accurately predicted efficacy of 28 drugs in killing glioma stem-like cells from 43 patients. Our findings demonstrate that interpatient heterogeneity manifests from differential activities of transcriptional programs, providing actionable strategies for mechanistically characterizing GBM from a systems perspective and developing better prognostic and predictive biomarkers for personalized medicine.

Individualized dynamic risk assessment for multiple myeloma

Background

Individualized treatment decisions for patients with multiple myeloma (MM) requires accurate risk stratification that takes into account patient-specific consequences of genetic abnormalities and tumor microenvironment on disease outcome and therapy responsiveness.

Methods

Previously, SYstems Genetic Network AnaLysis (SYGNAL) of multi-omics tumor profiles from 881 MM patients generated the mmSYGNAL network, which uncovered different causal and mechanistic drivers of genetic programs associated with disease progression across MM subtypes. Here, we have trained a machine learning (ML) algorithm on activities of mmSYGNAL programs within individual patient tumor samples to develop a risk classification scheme for MM that significantly outperformed cytogenetics, International Staging System, and multi-gene biomarker panels in predicting risk of PFS across four independent patient cohorts.

Results

We demonstrate that, unlike other tests, mmSYGNAL can accurately predict disease progression risk at primary diagnosis, pre- and post-transplant and even after multiple relapses, making it useful for individualized dynamic risk assessment throughout the disease trajectory.

Conclusion

mmSYGNAL provides improved individualized risk stratification that accounts for a patient's distinct set of genetic abnormalities and can monitor risk longitudinally as each patient's disease characteristics change.

Gene regulatory network topology governs resistance and treatment escape in glioma stem-like cells

Poor prognosis and drug resistance in glioblastoma (GBM) can result from cellular heterogeneity and treatment-induced shifts in phenotypic states of tumor cells, including dedifferentiation into glioma stem-like cells (GSCs). This rare tumorigenic cell subpopulation resists temozolomide, undergoes proneural-to-mesenchymal transition (PMT) to evade therapy, and drives recurrence. Through inference of transcriptional regulatory networks (TRNs) of patient-derived GSCs (PD-GSCs) at single-cell resolution, we demonstrate how the topology of transcription factor interaction networks drives distinct trajectories of cell state transitions in PD-GSCs resistant or susceptible to cytotoxic drug treatment. By experimentally testing predictions based on TRN simulations, we show that drug treatment drives surviving PD-GSCs along a trajectory of intermediate states, exposing vulnerability to potentiated killing by siRNA or a second drug targeting treatment-induced transcriptional programs governing non-genetic cell plasticity. Our findings demonstrate an approach to uncover TRN topology and use it to rationally predict combinatorial treatments that disrupts acquired resistance in GBM.

Personalized Clostridioides difficile engraftment risk prediction and probiotic therapy assessment in the human gut

Clostridioides difficile colonizes up to 30-40% of community-dwelling adults without causing disease. C. difficile infections (CDIs) are the leading cause of antibiotic-associated diarrhea in the U.S. and typically develop in individuals following disruption of the gut microbiota due to antibiotic or chemotherapy treatments. Further treatment of CDI with antibiotics is not always effective and can lead to antibiotic resistance and recurrent infections (rCDI). The most effective treatment for rCDI is the reestablishment of an intact microbiota via fecal microbiota transplants (FMTs). However, the success of FMTs has been difficult to generalize because the microbial interactions that prevent engraftment and facilitate the successful clearance of C. difficile are still only partially understood. Here we show how microbial community-scale metabolic models (MCMMs) accurately predicted known instances of C. difficile colonization susceptibility or resistance in vitro and in vivo. MCMMs provide detailed mechanistic insights into the ecological interactions that govern C. difficile engraftment, like cross-feeding or competition involving metabolites like succinate, trehalose, and ornithine, which differ from person to person. Indeed, three distinct C. difficile metabolic niches emerge from our MCMMs, two associated with positive growth rates and one that represents non-growth, which are consistently observed across 15,204 individuals from five independent cohorts. Finally, we show how MCMMs can predict personalized engraftment and C. difficile growth suppression for a probiotic cocktail (VE303) designed to replace FMTs for the treatment rCDI. Overall, this powerful modeling approach predicts personalized C. difficile engraftment risk and can be leveraged to assess probiotic treatment efficacy. MCMMs could be extended to understand the mechanistic underpinnings of personalized engraftment of other opportunistic bacterial pathogens, beneficial probiotic organisms, or more complex microbial consortia.

Genetic manipulation of Patescibacteria provides mechanistic insights into microbial dark matter and the epibiotic lifestyle

Patescibacteria, also known as the candidate phyla radiation (CPR), are a diverse group of bacteria that constitute a disproportionately large fraction of microbial dark matter. Its few cultivated members, belonging mostly to Saccharibacteria, grow as epibionts on host Actinobacteria. Due to a lack of suitable tools, the genetic basis of this lifestyle and other unique features of Patescibacteira remain unexplored. Here, we show that Saccharibacteria exhibit natural competence, and we exploit this property for their genetic manipulation. Imaging of fluorescent protein-labeled Saccharibacteria provides high spatiotemporal resolution of phenomena accompanying epibiotic growth, and a transposon-insertion sequencing (Tn-seq) genome-wide screen reveals the contribution of enigmatic Saccharibacterial genes to growth on their hosts. Finally, we leverage metagenomic data to provide cutting-edge protein structure-based bioinformatic resources that support the strain Southlakia epibionticum and its corresponding host, Actinomyces israelii, as a model system for unlocking the molecular underpinnings of the epibiotic lifestyle.

A comprehensive spectral assay library to quantify the Halobacterium salinarum NRC-1 proteome by DIA/SWATH-MS

Data-Independent Acquisition (DIA) is a mass spectrometry-based method to reliably identify and reproducibly quantify large fractions of a target proteome. The peptide-centric data analysis strategy employed in DIA requires a priori generated spectral assay libraries. Such assay libraries allow to extract quantitative data in a targeted approach and have been generated for human, mouse, zebrafish, E. coli and few other organisms. However, a spectral assay library for the extreme halophilic archaeon Halobacterium salinarum NRC-1, a model organism that contributed to several notable discoveries, is not publicly available yet. Here, we report a comprehensive spectral assay library to measure 2,563 of 2,646 annotated H. salinarum NRC-1 proteins. We demonstrate the utility of this library by measuring global protein abundances over time under standard growth conditions. The H. salinarum NRC-1 library includes 21,074 distinct peptides representing 97% of the predicted proteome and provides a new, valuable resource to confidently measure and quantify any protein of this archaeon. Data and spectral assay libraries are available via ProteomeXchange (PXD042770, PXD042774) and SWATHAtlas (SAL00312-SAL00319).

Growth phase estimation for abundant bacterial populations sampled longitudinally from human stool metagenomes

Longitudinal sampling of the stool has yielded important insights into the ecological dynamics of the human gut microbiome. However, human stool samples are available approximately once per day, while commensal population doubling times are likely on the order of minutes-to-hours. Despite this mismatch in timescales, much of the prior work on human gut microbiome time series modeling has assumed that day-to-day fluctuations in taxon abundances are related to population growth or death rates, which is likely not the case. Here, we propose an alternative model of the human gut as a stationary system, where population dynamics occur internally and the bacterial population sizes measured in a bolus of stool represent a steady-state endpoint of these dynamics. We formalize this idea as stochastic logistic growth. We show how this model provides a path toward estimating the growth phases of gut bacterial populations in situ. We validate our model predictions using an in vitro Escherichia coli growth experiment. Finally, we show how this method can be applied to densely-sampled human stool metagenomic time series data. We discuss how these growth phase estimates may be used to better inform metabolic modeling in flow-through ecosystems, like animal guts or industrial bioreactors.

MtrA modulates Mycobacterium tuberculosis cell division in host microenvironments to mediate intrinsic resistance and drug tolerance

The success of Mycobacterium tuberculosis (Mtb) is largely attributed to its ability to physiologically adapt and withstand diverse localized stresses within host microenvironments. Here, we present a data-driven model (EGRIN 2.0) that captures the dynamic interplay of environmental cues and genome-encoded regulatory programs in Mtb. Analysis of EGRIN 2.0 shows how modulation of the MtrAB two-component signaling system tunes Mtb growth in response to related host microenvironmental cues. Disruption of MtrAB by tunable CRISPR interference confirms that the signaling system regulates multiple peptidoglycan hydrolases, among other targets, that are important for cell division. Further, MtrA decreases the effectiveness of antibiotics by mechanisms of both intrinsic resistance and drug tolerance. Together, the model-enabled dissection of complex MtrA regulation highlights its importance as a drug target and illustrates how EGRIN 2.0 facilitates discovery and mechanistic characterization of Mtb adaptation to specific host microenvironments within the host.

Discovery of a glutathione utilization pathway in Francisella that shows functional divergence between environmental and pathogenic species

Glutathione (GSH) is an abundant metabolite within eukaryotic cells that can act as a signal, a nutrient source, or serve in a redox capacity for intracellular bacterial pathogens. For Francisella, GSH is thought to be a critical in vivo source of cysteine; however, the cellular pathways permitting GSH utilization by Francisella differ between strains and have remained poorly understood. Using genetic screening, we discovered a unique pathway for GSH utilization in Francisella. Whereas prior work suggested GSH catabolism initiates in the periplasm, the pathway we define consists of a major facilitator superfamily (MFS) member that transports intact GSH and a previously unrecognized bacterial cytoplasmic enzyme that catalyzes the first step of GSH degradation. Interestingly, we find that the transporter gene for this pathway is pseudogenized in pathogenic Francisella, explaining phenotypic discrepancies in GSH utilization among Francisella spp. and revealing a critical role for GSH in the environmental niche of these bacteria.

Genetic manipulation of candidate phyla radiation bacteria provides functional insights into microbial dark matter

The study of bacteria has yielded fundamental insights into cellular biology and physiology, biotechnological advances and many therapeutics. Yet due to a lack of suitable tools, the significant portion of bacterial diversity held within the candidate phyla radiation (CPR) remains inaccessible to such pursuits. Here we show that CPR bacteria belonging to the phylum Saccharibacteria exhibit natural competence. We exploit this property to develop methods for their genetic manipulation, including the insertion of heterologous sequences and the construction of targeted gene deletions. Imaging of fluorescent protein-labeled Saccharibacteria provides high spatiotemporal resolution of phenomena accompanying epibiotic growth and a transposon insertion sequencing genome-wide screen reveals the contribution of enigmatic Saccharibacterial genes to growth on their Actinobacteria hosts. Finally, we leverage metagenomic data to provide cutting-edge protein structure-based bioinformatic resources that support the strain Southlakia epibionticum and its corresponding host, Actinomyces israelii , as a model system for unlocking the molecular underpinnings of the epibiotic lifestyle.

A single-cell based precision medicine approach using glioblastoma patient-specific models

Glioblastoma (GBM) is a heterogeneous tumor made up of cell states that evolve over time. Here, we modeled tumor evolutionary trajectories during standard-of-care treatment using multi-omic single-cell analysis of a primary tumor sample, corresponding mouse xenografts subjected to standard of care therapy, and recurrent tumor at autopsy. We mined the multi-omic data with single-cell SYstems Genetics Network AnaLysis (scSYGNAL) to identify a network of 52 regulators that mediate treatment-induced shifts in xenograft tumor-cell states that were also reflected in recurrence. By integrating scSYGNAL-derived regulatory network information with transcription factor accessibility deviations derived from single-cell ATAC-seq data, we developed consensus networks that modulate cell state transitions across subpopulations of primary and recurrent tumor cells. Finally, by matching targeted therapies to active regulatory networks underlying tumor evolutionary trajectories, we provide a framework for applying single-cell-based precision medicine approaches to an individual patient in a concurrent, adjuvant, or recurrent setting.

A non-canonical role of caspase-1 in regulating bacterial physiology and antimicrobial resistance

Caspase-1 is a key effector molecule involved in inflammasome activation and has a well-established role in restricting the growth of intracellular pathogens like Salmonella by triggering a form of cell death called pyroptosis. Here we reveal a non-canonical, cell death independent role for caspase-1 in controlling the transcriptional state and drug resistance of intracellular Salmonella. Using Pathogen-sequencing, a method for sensitive transcriptional profiling of miniscule numbers of intracellular bacteria from infected macrophages, we show that that caspase-1 regulates key processes involved in survival and antimicrobial susceptibility of intracellular Salmonella. Host caspase-1 increased susceptibility of Salmonella to endogenous cationic antimicrobial peptides, as well as to a cationic polypeptide antibiotic used as a last-line drug in Gram-negative bacterial infections. These effects of caspase-1 were independent of its enzymatic activity but dependent on its ability to repress activation of a two-component signal transduction system in intracellular bacteria. These effects were also independent of caspase-11. Our data suggest a “backup” role for caspase-1 in dampening antimicrobial resistance of intracellular Salmonella which evade initial innate immune detection and restriction by caspase-1. These findings also reveal the role of a key innate immune effector in altering pathogen physiology, broadening our view of host-pathogen crosstalk with possible implications for targeting caspase-1 in host-directed therapy to combat antimicrobial resistance.

Supported by grant from NIH (R01 AI155685-A1 ) and AAI Careers in Immunology Fellowship (2018-2019)

Mathematical models to study the biology of pathogens and the infectious diseases they cause

Mathematical models have many applications in infectious diseases: epidemiologists use them to forecast outbreaks and design containment strategies; systems biologists use them to study complex processes sustaining pathogens, from the metabolic networks empowering microbial cells to ecological networks in the microbiome that protects its host. Here, we (1) review important models relevant to infectious diseases, (2) draw parallels among models ranging widely in scale. We end by discussing a minimal set of information for a model to promote its use by others and to enable predictions that help us better fight pathogens and the diseases they cause.

Transcriptome signature of cell viability predicts drug response and drug interaction in Mycobacterium tuberculosis

There is an urgent need for new drug regimens to rapidly cure tuberculosis. Here, we report the development of drug response assayer (DRonA) and "MLSynergy," algorithms to perform rapid drug response assays and predict response of Mycobacterium tuberculosis (Mtb) to drug combinations. Using a transcriptome signature for cell viability, DRonA detects Mtb killing by diverse mechanisms in broth culture, macrophage infection, and patient sputum, providing an efficient and more sensitive alternative to time- and resource-intensive bacteriologic assays. Further, MLSynergy builds on DRonA to predict synergistic and antagonistic multidrug combinations using transcriptomes of Mtb treated with single drugs. Together, DRonA and MLSynergy represent a generalizable framework for rapid monitoring of drug effects in host-relevant contexts and accelerate the discovery of efficacious high-order drug combinations.

Quantitative prediction of conditional vulnerabilities in regulatory and metabolic networks using PRIME

The ability of Mycobacterium tuberculosis (Mtb) to adopt heterogeneous physiological states underlies its success in evading the immune system and tolerating antibiotic killing. Drug tolerant phenotypes are a major reason why the tuberculosis (TB) mortality rate is so high, with over 1.8 million deaths annually. To develop new TB therapeutics that better treat the infection (faster and more completely), a systems-level approach is needed to reveal the complexity of network-based adaptations of Mtb. Here, we report a new predictive model called PRIME (Phenotype of Regulatory influences Integrated with Metabolism and Environment) to uncover environment-specific vulnerabilities within the regulatory and metabolic networks of Mtb. Through extensive performance evaluations using genome-wide fitness screens, we demonstrate that PRIME makes mechanistically accurate predictions of context-specific vulnerabilities within the integrated regulatory and metabolic networks of Mtb, accurately rank-ordering targets for potentiating treatment with frontline drugs.

In vivo commensal control of Clostridioides difficile virulence

Leveraging systems biology approaches, we illustrate how metabolically distinct species of Clostridia protect against or worsen Clostridioides difficile infection in mice by modulating the pathogen's colonization, growth, and virulence to impact host survival. Gnotobiotic mice colonized with the amino acid fermenter Paraclostridium bifermentans survive infection with reduced disease severity, while mice colonized with the butyrate-producer, Clostridium sardiniense, succumb more rapidly. Systematic in vivo analyses revealed how each commensal alters the gut-nutrient environment to modulate the pathogen's metabolism, gene regulatory networks, and toxin production. Oral administration of P. bifermentans rescues conventional, clindamycin-treated mice from lethal C. difficile infection in a manner similar to that of monocolonized animals, thereby supporting the therapeutic potential of this commensal species. Our findings lay the foundation for mechanistically informed therapies to counter C. difficile disease using systems biology approaches to define host-commensal-pathogen interactions in vivo.

Predictive regulatory and metabolic network models for systems analysis of Clostridioides difficile

We present predictive models for comprehensive systems analysis of Clostridioides difficile, the etiology of pseudomembranous colitis. By leveraging 151 published transcriptomes, we generated an EGRIN model that organizes 90% of C. difficile genes into a transcriptional regulatory network of 297 co-regulated modules, implicating genes in sporulation, carbohydrate transport, and metabolism. By advancing a metabolic model through addition and curation of metabolic reactions including nutrient uptake, we discovered 14 amino acids, diverse carbohydrates, and 10 metabolic genes as essential for C. difficile growth in the intestinal environment. Finally, we developed a PRIME model to uncover how EGRIN-inferred combinatorial gene regulation by transcription factors, such as CcpA and CodY, modulates essential metabolic processes to enable C. difficile growth relative to commensal colonization. The C. difficile interactive web portal provides access to these model resources to support collaborative systems-level studies of context-specific virulence mechanisms in C. difficile.

Autophagy as a Mechanism for Adaptive Prediction-Mediated Emergence of Drug Resistance

Drug resistance is a major problem in treatment of microbial infections and cancers. There is growing evidence that a transient drug tolerant state may precede and potentiate the emergence of drug resistance. Therefore, understanding the mechanisms leading to tolerance is critical for combating drug resistance and for the development of effective therapeutic strategy. Through laboratory evolution of yeast, we recently demonstrated that adaptive prediction (AP), a strategy employed by organisms to anticipate and prepare for a future stressful environment, can emerge within 100 generations by linking the response triggered by a neutral cue (caffeine) to a mechanism of protection against a lethal agent (5-fluoroorotic acid, 5-FOA). Here, we demonstrate that mutations selected across multiple laboratory-evolved lines had linked the neutral cue response to core genes of autophagy. Across these evolved lines, conditional activation of autophagy through AP conferred tolerance, and potentiated subsequent selection of mutations in genes specific to overcoming the toxicity of 5-FOA. These results offer a new perspective on how extensive genome-wide genetic interactions of autophagy could have facilitated the emergence of AP over short evolutionary timescales to potentiate selection of 5-FOA resistance-conferring mutations.

Contrasting heat stress response patterns of coral holobionts across the Red Sea suggest distinct mechanisms of thermal tolerance

Corals from the northern Red Sea, in particular the Gulf of Aqaba (GoA), have exceptionally high bleaching thresholds approaching >5℃ above their maximum monthly mean (MMM) temperatures. These elevated thresholds are thought to be due to historical selection, as corals passed through the warmer Southern Red Sea during recolonization from the Arabian Sea. To test this hypothesis, we determined thermal tolerance thresholds of GoA versus central Red Sea (CRS) Stylophora pistillata corals using multi-temperature acute thermal stress assays to determine thermal thresholds. Relative thermal thresholds of GoA and CRS corals were indeed similar and exceptionally high (~7℃ above MMM). However, absolute thermal thresholds of CRS corals were on average 3℃ above those of GoA corals. To explore the molecular underpinnings, we determined gene expression and microbiome response of the coral holobiont. Transcriptomic responses differed markedly, with a strong response to the thermal stress in GoA corals and their symbiotic algae versus a remarkably muted response in CRS colonies. Concomitant to this, coral and algal genes showed temperature-induced expression in GoA corals, while exhibiting fixed high expression (front-loading) in CRS corals. Bacterial community composition of GoA corals changed dramatically under heat stress, whereas CRS corals displayed stable assemblages. We interpret the response of GoA corals as that of a resilient population approaching a tipping point in contrast to a pattern of consistently elevated thermal resistance in CRS corals that cannot further attune. Such response differences suggest distinct thermal tolerance mechanisms that may affect the response of coral populations to ocean warming.

Synergistic epistasis enhances the co-operativity of mutualistic interspecies interactions

Early evolution of mutualism is characterized by big and predictable adaptive changes, including the specialization of interacting partners, such as through deleterious mutations in genes not required for metabolic cross-feeding. We sought to investigate whether these early mutations improve cooperativity by manifesting in synergistic epistasis between genomes of the mutually interacting species. Specifically, we have characterized evolutionary trajectories of syntrophic interactions of Desulfovibrio vulgaris (Dv) with Methanococcus maripaludis (Mm) by longitudinally monitoring mutations accumulated over 1000 generations of nine independently evolved communities with analysis of the genotypic structure of one community down to the single-cell level. We discovered extensive parallelism across communities despite considerable variance in their evolutionary trajectories and the perseverance within many evolution lines of a rare lineage of Dv that retained sulfate-respiration (SR+) capability, which is not required for metabolic cross-feeding. An in-depth investigation revealed that synergistic epistasis across pairings of Dv and Mm genotypes had enhanced cooperativity within SR- and SR+ assemblages, enabling their coexistence within the same community. Thus, our findings demonstrate that cooperativity of a mutualism can improve through synergistic epistasis between genomes of the interacting species, enabling the coexistence of mutualistic assemblages of generalists and their specialized variants.

A Systems Approach to Brain Tumor Treatment

Brain tumors are among the most lethal tumors. Glioblastoma, the most frequent primary brain tumor in adults, has a median survival time of approximately 15 months after diagnosis or a five-year survival rate of 10%; the recurrence rate is nearly 90%. Unfortunately, this prognosis has not improved for several decades. The lack of progress in the treatment of brain tumors has been attributed to their high rate of primary therapy resistance. Challenges such as pronounced inter-patient variability, intratumoral heterogeneity, and drug delivery across the blood-brain barrier hinder progress. A comprehensive, multiscale understanding of the disease, from the molecular to the whole tumor level, is needed to address the intratumor heterogeneity resulting from the coexistence of a diversity of neoplastic and non-neoplastic cell types in the tumor tissue. By contrast, inter-patient variability must be addressed by subtyping brain tumors to stratify patients and identify the best-matched drug(s) and therapies for a particular patient or cohort of patients. Accomplishing these diverse tasks will require a new framework, one involving a systems perspective in assessing the immense complexity of brain tumors. This would in turn entail a shift in how clinical medicine interfaces with the rapidly advancing high-throughput (HTP) technologies that have enabled the omics-scale profiling of molecular features of brain tumors from the single-cell to the tissue level. However, several gaps must be closed before such a framework can fulfill the promise of precision and personalized medicine for brain tumors. Ultimately, the goal is to integrate seamlessly multiscale systems analyses of patient tumors and clinical medicine. Accomplishing this goal would facilitate the rational design of therapeutic strategies matched to the characteristics of patients and their tumors. Here, we discuss some of the technologies, methodologies, and computational tools that will facilitate the realization of this vision to practice.

Intricate Genetic Programs Controlling Dormancy in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) displays the remarkable ability to transition in and out of dormancy, a hallmark of the pathogen’s capacity to evade the immune system and exploit susceptible individuals. Uncovering the gene regulatory programs that underlie the phenotypic shifts in MTB during disease latency and reactivation has posed a challenge. We develop an experimental system to precisely control dissolved oxygen levels in MTB cultures in order to capture the transcriptional events that unfold as MTB transitions into and out of hypoxia-induced dormancy. Using a comprehensive genome-wide transcription factor binding map and insights from network topology analysis, we identify regulatory circuits that deterministically drive sequential transitions across six transcriptionally and functionally distinct states encompassing more than three-fifths of the MTB genome. The architecture of the genetic programs explains the transcriptional dynamics underlying synchronous entry of cells into a dormant state that is primed to infect the host upon encountering favorable conditions.

Mycobacterium tuberculosis (MTB) persists within the host by counteracting disparate stressors including hypoxia. Peterson et al. report a transcriptional program that coordinates sequential state transitions to drive MTB in and out of hypoxia-induced dormancy. Among varied properties, this program encodes advanced preparedness to infect the host in favorable conditions.

History and Future Perspectives on the Discipline of Quantitative Systems Pharmacology Modeling and Its Applications

Mathematical biology and pharmacology models have a long and rich history in the fields of medicine and physiology, impacting our understanding of disease mechanisms and the development of novel therapeutics. With an increased focus on the pharmacology application of system models and the advances in data science spanning mechanistic and empirical approaches, there is a significant opportunity and promise to leverage these advancements to enhance the development and application of the systems pharmacology field. In this paper, we will review milestones in the evolution of mathematical biology and pharmacology models, highlight some of the gaps and challenges in developing and applying systems pharmacology models, and provide a vision for an integrated strategy that leverages advances in adjacent fields to overcome these challenges.

Mechanism Across Scales: A Holistic Modeling Framework Integrating Laboratory and Field Studies for Microbial Ecology

Over the last century, leaps in technology for imaging, sampling, detection, high-throughput sequencing, and -omics analyses have revolutionized microbial ecology to enable rapid acquisition of extensive datasets for microbial communities across the ever-increasing temporal and spatial scales. The present challenge is capitalizing on our enhanced abilities of observation and integrating diverse data types from different scales, resolutions, and disciplines to reach a causal and mechanistic understanding of how microbial communities transform and respond to perturbations in the environment. This type of causal and mechanistic understanding will make predictions of microbial community behavior more robust and actionable in addressing microbially mediated global problems. To discern drivers of microbial community assembly and function, we recognize the need for a conceptual, quantitative framework that connects measurements of genomic potential, the environment, and ecological and physical forces to rates of microbial growth at specific locations. We describe the Framework for Integrated, Conceptual, and Systematic Microbial Ecology (FICSME), an experimental design framework for conducting process-focused microbial ecology studies that incorporates biological, chemical, and physical drivers of a microbial system into a conceptual model. Through iterative cycles that advance our understanding of the coupling across scales and processes, we can reliably predict how perturbations to microbial systems impact ecosystem-scale processes or vice versa. We describe an approach and potential applications for using the FICSME to elucidate the mechanisms of globally important ecological and physical processes, toward attaining the goal of predicting the structure and function of microbial communities in chemically complex natural environments.

Lettuce (Lactuca sativa) productivity influenced by microbial inocula under nitrogen-limited conditions in aquaponics

The demand for food will outpace productivity of conventional agriculture due to projected growth of the human population, concomitant with shrinkage of arable land, increasing scarcity of freshwater, and a rapidly changing climate. While aquaponics has potential to sustainably supplement food production with minimal environmental impact, there is a need to better characterize the complex interplay between the various components (fish, plant, microbiome) of these systems to optimize scale up and productivity. Here, we investigated how the commonly-implemented practice of continued microbial community transfer from pre-existing systems might promote or impede productivity of aquaponics. Specifically, we monitored plant growth phenotypes, water chemistry, and microbiome composition of rhizospheres, biofilters, and fish feces over 61-days of lettuce (Lactuca sativa var. crispa) growth in nitrogen-limited aquaponic systems inoculated with bacteria that were either commercially sourced or originating from a pre-existing aquaponic system. Lettuce above- and below-ground growth were significantly reduced across replicates treated with a pre-existing aquaponic system inoculum when compared to replicates treated with a commercial inoculum. Reduced productivity was associated with enrichment in specific bacterial genera in plant roots, including Pseudomonas, following inoculum transfer from pre-existing systems. Increased productivity was associated with enrichment of nitrogen-fixing Rahnella in roots of plants treated with the commercial inoculum. Thus, we show that inoculation from a pre-existing system, rather than from a commercial inoculum, is associated with lower yields. Further work will be necessary to test the putative mechanisms involved.

A comprehensive spectral assay library to quantify the Escherichia coli proteome by DIA/SWATH-MS

Data-Independent Acquisition (DIA) is a method to improve consistent identification and precise quantitation of peptides and proteins by mass spectrometry (MS). The targeted data analysis strategy in DIA relies on spectral assay libraries that are generally derived from a priori measurements of peptides for each species. Although Escherichia coli (E. coli) is among the best studied model organisms, so far there is no spectral assay library for the bacterium publicly available. Here, we generated a spectral assay library for 4,014 of the 4,389 annotated E. coli proteins using one- and two-dimensional fractionated samples, and ion mobility separation enabling deep proteome coverage. We demonstrate the utility of this high-quality library with robustness in quantitation of the E. coli proteome and with rapid-chromatography to enhance throughput by targeted DIA-MS. The spectral assay library supports the detection and quantification of 91.5% of all E. coli proteins at high-confidence with 56,182 proteotypic peptides, making it a valuable resource for the scientific community. Data and spectral libraries are available via ProteomeXchange (PXD020761, PXD020785) and SWATHAtlas (SAL00222-28).

Crippling life support for SARS-CoV-2 and other viruses through synthetic lethality

With the rapid global spread of SARS-CoV-2, we have become acutely aware of the inadequacies of our ability to respond to viral epidemics. Although disrupting the viral life cycle is critical for limiting viral spread and disease, it has proven challenging to develop targeted and selective therapeutics. Synthetic lethality offers a promising but largely unexploited strategy against infectious viral disease; as viruses infect cells, they abnormally alter the cell state, unwittingly exposing new vulnerabilities in the infected cell. Therefore, we propose that effective therapies can be developed to selectively target the virally reconfigured host cell networks that accompany altered cellular states to cripple the host cell that has been converted into a virus factory, thus disrupting the viral life cycle.

The State of Systems Genetics in 2017

Cell Systems invited 16 experts to share their views on the field of systems genetics. In questions repeated in the headings, we asked them to define systems genetics, highlight its relevance to researchers outside the field, discuss what makes a strong systems genetics paper, and paint a picture of where the field is heading in the coming years. Their responses, ordered by the journal but otherwise unedited, make it clear that deciphering genotype to phenotype relationships is a central challenge of systems genetics and will require understanding how networks and higher-order properties of biological systems underlie complex traits. In addition, our experts illuminate the applications and relevance of systems genetics to human disease, the gut microbiome, development of tools that connect the global research community, sustainability, drug discovery, patient-specific disease and network models, and personalized treatments. Finally, a table of suggested reading provides a sample of influential work in the field.

EcoFABs: advancing microbiome science through standardized fabricated ecosystems

Microbiomes play critical roles in ecosystems and human health, yet in most cases scientists lack standardized and reproducible model microbial communities. The development of fabricated microbial ecosystems, which we term EcoFABs, will provide such model systems for microbiome studies.

Path-seq identifies an essential mycolate remodeling program for mycobacterial host adaptation

The success of Mycobacterium tuberculosis (MTB) stems from its ability to remain hidden from the immune system within macrophages. Here, we report a new technology (Path‐seq) to sequence miniscule amounts of MTB transcripts within up to million‐fold excess host RNA. Using Path‐seq and regulatory network analyses, we have discovered a novel transcriptional program for in vivo mycobacterial cell wall remodeling when the pathogen infects alveolar macrophages in mice. We have discovered that MadR transcriptionally modulates two mycolic acid desaturases desA1/desA2 to initially promote cell wall remodeling upon in vitro macrophage infection and, subsequently, reduces mycolate biosynthesis upon entering dormancy. We demonstrate that disrupting MadR program is lethal to diverse mycobacteria making this evolutionarily conserved regulator a prime antitubercular target for both early and late stages of infection.

Ocean acidification conditions increase resilience of marine diatoms

The fate of diatoms in future acidified oceans could have dramatic implications on marine ecosystems, because they account for ~40% of marine primary production. Here, we quantify resilience of Thalassiosira pseudonana in mid-20th century (300 ppm CO2) and future (1000 ppm CO2) conditions that cause ocean acidification, using a stress test that probes its ability to recover from incrementally higher amount of low-dose ultraviolet A (UVA) and B (UVB) radiation and re-initiate growth in day-night cycles, limited by nitrogen. While all cultures eventually collapse, those growing at 300 ppm CO2 succumb sooner. The underlying mechanism for collapse appears to be a system failure resulting from "loss of relational resilience," that is, inability to adopt physiological states matched to N-availability and phase of the diurnal cycle. Importantly, under elevated CO2 conditions diatoms sustain relational resilience over a longer timeframe, demonstrating increased resilience to future acidified ocean conditions. This stress test framework can be extended to evaluate and predict how various climate change associated stressors may impact microbial community resilience.

The BCGΔBCG1419c Vaccine Candidate Reduces Lung Pathology, IL-6, TNF-α, and IL-10 During Chronic TB Infection

Mycobacterium tuberculosis (M. tuberculosis), the causative agent of human tuberculosis (TB), is estimated to be harbored by up to 2 billion people in a latent TB infection (LTBI) state. The only TB vaccine approved for use in humans, BCG, does not confer protection against establishment of or reactivation from LTBI, so new vaccine candidates are needed to specifically address this need. Following the hypothesis that mycobacterial biofilms resemble aspects of LTBI, we modified BCG by deleting the BCG1419c gene to create the BCGΔBCG1419c vaccine strain. In this study, we compared cytokine profiles, bacterial burden, and lung lesions after immunization with BCG or BCGΔBCG1419c before and after 6 months of aerosol infection with M. tuberculosis H37Rv in the resistant C57BL/6 mouse model. Our results show that in infected mice, BCGΔBCG1419c significantly reduced lung lesions and IL-6 in comparison to the unmodified BCG strain, and was the only vaccine that decreased production of TNF-α and IL-10 compared to non-vaccinated mice, while vaccination with BCG or BCGΔBCG1419c significantly reduced IFN-γ production. Moreover, transcriptome profiling of BCGΔBCG1419c suggests that compared to BCG, it has decreased expression of genes involved in mycolic acids (MAs) metabolism, and antigenic chaperones, which might be involved in reduced pathology compared to BCG-vaccinated mice.

Transcriptional regulatory networks underlying gene expression changes in Huntington's disease

Transcriptional changes occur presymptomatically and throughout Huntington's disease (HD), motivating the study of transcriptional regulatory networks (TRNs) in HD We reconstructed a genome-scale model for the target genes of 718 transcription factors (TFs) in the mouse striatum by integrating a model of genomic binding sites with transcriptome profiling of striatal tissue from HD mouse models. We identified 48 differentially expressed TF-target gene modules associated with age- and CAG repeat length-dependent gene expression changes in Htt CAG knock-in mouse striatum and replicated many of these associations in independent transcriptomic and proteomic datasets. Thirteen of 48 of these predicted TF-target gene modules were also differentially expressed in striatal tissue from human disease. We experimentally validated a specific model prediction that SMAD3 regulates HD-related gene expression changes using chromatin immunoprecipitation and deep sequencing (ChIP-seq) of mouse striatum. We found CAG repeat length-dependent changes in the genomic occupancy of SMAD3 and confirmed our model's prediction that many SMAD3 target genes are downregulated early in HD.

Adaptive Prediction Emerges Over Short Evolutionary Time Scales

Adaptive prediction is a capability of diverse organisms, including microbes, to sense a cue and prepare in advance to deal with a future environmental challenge. Here, we investigated the timeframe over which adaptive prediction emerges when an organism encounters an environment with novel structure. We subjected yeast to laboratory evolution in a novel environment with repetitive, coupled exposures to a neutral chemical cue (caffeine), followed by a sublethal dose of a toxin (5-FOA), with an interspersed requirement for uracil prototrophy to counter-select mutants that gained constitutive 5-FOA resistance. We demonstrate the remarkable ability of yeast to internalize a novel environmental pattern within 50-150 generations by adaptively predicting 5-FOA stress upon sensing caffeine. We also demonstrate how novel environmental structure can be internalized by coupling two unrelated response networks, such as the response to caffeine and signaling-mediated conditional peroxisomal localization of proteins.

Robustness of a model microbial community emerges from population structure among single cells of a clonal population

Microbial populations can withstand, overcome and persist in the face of environmental fluctuation. Previously, we demonstrated how conditional gene regulation in a fluctuating environment drives dilution of condition-specific transcripts, causing a population of Desulfovibrio vulgaris Hildenborough (DvH) to collapse after repeatedly transitioning from sulfate respiration to syntrophic conditions with the methanogen Methanococcus maripaludis. Failure of the DvH to successfully transition contributed to the collapse of this model community. We investigated the mechanistic basis for loss of robustness by examining whether conditional gene regulation altered heterogeneity in gene expression across individual DvH cells. We discovered that robustness of a microbial population across environmental transitions was attributable to the retention of cells in two states that exhibited different condition-specific gene expression patterns. In our experiments, a population with disrupted conditional regulation successfully alternated between cell states. Meanwhile, a population with intact conditional regulation successfully switched between cell states initially, but collapsed after repeated transitions, possibly due to the high energy requirements of regulation. These results demonstrate that the survival of this entire model microbial community is dependent on the regulatory system's influence on the distribution of distinct cell states among individual cells within a clonal population.

Combining inferred regulatory and reconstructed metabolic networks enhances phenotype prediction in yeast

Gene regulatory and metabolic network models have been used successfully in many organisms, but inherent differences between them make networks difficult to integrate. Probabilistic Regulation Of Metabolism (PROM) provides a partial solution, but it does not incorporate network inference and underperforms in eukaryotes. We present an Integrated Deduced And Metabolism (IDREAM) method that combines statistically inferred Environment and Gene Regulatory Influence Network (EGRIN) models with the PROM framework to create enhanced metabolic-regulatory network models. We used IDREAM to predict phenotypes and genetic interactions between transcription factors and genes encoding metabolic activities in the eukaryote, Saccharomyces cerevisiae. IDREAM models contain many fewer interactions than PROM and yet produce significantly more accurate growth predictions. IDREAM consistently outperformed PROM using any of three popular yeast metabolic models and across three experimental growth conditions. Importantly, IDREAM's enhanced accuracy makes it possible to identify subtle synthetic growth defects. With experimental validation, these novel genetic interactions involving the pyruvate dehydrogenase complex suggested a new role for fatty acid-responsive factor Oaf1 in regulating acetyl-CoA production in glucose grown cells.

Mechanism for microbial population collapse in a fluctuating resource environment

Managing trade-offs through gene regulation is believed to confer resilience to a microbial community in a fluctuating resource environment. To investigate this hypothesis, we imposed a fluctuating environment that required the sulfate-reducer Desulfovibrio vulgaris to undergo repeated ecologically relevant shifts between retaining metabolic independence (active capacity for sulfate respiration) and becoming metabolically specialized to a mutualistic association with the hydrogen-consuming Methanococcus maripaludis Strikingly, the microbial community became progressively less proficient at restoring the environmentally relevant physiological state after each perturbation and most cultures collapsed within 3-7 shifts. Counterintuitively, the collapse phenomenon was prevented by a single regulatory mutation. We have characterized the mechanism for collapse by conducting RNA-seq analysis, proteomics, microcalorimetry, and single-cell transcriptome analysis. We demonstrate that the collapse was caused by conditional gene regulation, which drove precipitous decline in intracellular abundance of essential transcripts and proteins, imposing greater energetic burden of regulation to restore function in a fluctuating environment.

The Transcription Factor Nfatc2 Regulates β-Cell Proliferation and Genes Associated with Type 2 Diabetes in Mouse and Human Islets

Human genome-wide association studies (GWAS) have shown that genetic variation at >130 gene loci is associated with type 2 diabetes (T2D). We asked if the expression of the candidate T2D-associated genes within these loci is regulated by a common locus in pancreatic islets. Using an obese F2 mouse intercross segregating for T2D, we show that the expression of ~40% of the T2D-associated genes is linked to a broad region on mouse chromosome (Chr) 2. As all but 9 of these genes are not physically located on Chr 2, linkage to Chr 2 suggests a genomic factor(s) located on Chr 2 regulates their expression in trans. The transcription factor Nfatc2 is physically located on Chr 2 and its expression demonstrates cis linkage; i.e., its expression maps to itself. When conditioned on the expression of Nfatc2, linkage for the T2D-associated genes was greatly diminished, supporting Nfatc2 as a driver of their expression. Plasma insulin also showed linkage to the same broad region on Chr 2. Overexpression of a constitutively active (ca) form of Nfatc2 induced β-cell proliferation in mouse and human islets, and transcriptionally regulated more than half of the T2D-associated genes. Overexpression of either ca-Nfatc2 or ca-Nfatc1 in mouse islets enhanced insulin secretion, whereas only ca-Nfatc2 was able to promote β-cell proliferation, suggesting distinct molecular pathways mediating insulin secretion vs. β-cell proliferation are regulated by NFAT. Our results suggest that many of the T2D-associated genes are downstream transcriptional targets of NFAT, and may act coordinately in a pathway through which NFAT regulates β-cell proliferation in both mouse and human islets.

Genome-wide association studies (GWAS) and linkage studies provide a powerful way to establish a causal connection between a gene locus and a physiological or pathophysiological phenotype. We wondered if candidate genes associated with type 2 diabetes in human populations, in addition to being causal for the disease, could also be intermediate traits in a pathway leading to disease. In addition, we wished to know if there were any regulatory loci that could coordinately drive the expression of these genes in pancreatic islets and thus complete a pathway; i.e. Driver → GWAS candidate expression → type 2 diabetes. Using data from a mouse intercross between a diabetes-susceptible and a diabetes-resistant mouse strain, we found that the expression of ~40% of >130 candidate GWAS genes genetically mapped to a hot spot on mouse chromosome 2. Using a variety of statistical methods, we identified the transcription factor Nfatc2 as the candidate driver. Follow-up experiments showed that overexpression of Nfatc2 does indeed affect the expression of the GWAS genes and regulates β-cell proliferation and insulin secretion. The work shows that in addition to being causal, GWAS candidate genes can be intermediate traits in a pathway leading to disease. Model organisms can be used to explore these novel causal pathways.

A genome wide dosage suppressor network reveals genomic robustness

Genomic robustness is the extent to which an organism has evolved to withstand the effects of deleterious mutations. We explored the extent of genomic robustness in budding yeast by genome wide dosage suppressor analysis of 53 conditional lethal mutations in cell division cycle and RNA synthesis related genes, revealing 660 suppressor interactions of which 642 are novel. This collection has several distinctive features, including high co-occurrence of mutant-suppressor pairs within protein modules, highly correlated functions between the pairs and higher diversity of functions among the co-suppressors than previously observed. Dosage suppression of essential genes encoding RNA polymerase subunits and chromosome cohesion complex suggests a surprising degree of functional plasticity of macromolecular complexes, and the existence of numerous degenerate pathways for circumventing the effects of potentially lethal mutations. These results imply that organisms and cancer are likely able to exploit the genomic robustness properties, due the persistence of cryptic gene and pathway functions, to generate variation and adapt to selective pressures.

Causal Mechanistic Regulatory Network for Glioblastoma Deciphered Using Systems Genetics Network Analysis

We developed the transcription factor (TF)-target gene database and the Systems Genetics Network Analysis (SYGNAL) pipeline to decipher transcriptional regulatory networks from multi-omic and clinical patient data, and we applied these tools to 422 patients with glioblastoma multiforme (GBM). The resulting gbmSYGNAL network predicted 112 somatically mutated genes or pathways that act through 74 TFs and 37 microRNAs (miRNAs) (67 not previously associated with GBM) to dysregulate 237 distinct co-regulated gene modules associated with patient survival or oncogenic processes. The regulatory predictions were associated to cancer phenotypes using CRISPR-Cas9 and small RNA perturbation studies and also demonstrated GBM specificity. Two pairwise combinations (ETV6-NFKB1 and romidepsin-miR-486-3p) predicted by the gbmSYGNAL network had synergistic anti-proliferative effects. Finally, the network revealed that mutations in NF1 and PIK3CA modulate IRF1-mediated regulation of MHC class I antigen processing and presentation genes to increase tumor lymphocyte infiltration and worsen prognosis. Importantly, SYGNAL is widely applicable for integrating genomic and transcriptomic measurements from other human cohorts.

Sense overlapping transcripts in IS1341-type transposase genes are functional non-coding RNAs in archaea

The existence of sense overlapping transcripts that share regulatory and coding information in the same genomic sequence shows an additional level of prokaryotic gene expression complexity. Here we report the discovery of ncRNAs associated with IS1341-type transposase (tnpB) genes, at the 3'-end of such elements, with examples in archaea and bacteria. Focusing on the model haloarchaeon Halobacterium salinarum NRC-1, we show the existence of sense overlapping transcripts (sotRNAs) for all its IS1341-type transposases. Publicly available transcriptome compendium show condition-dependent differential regulation between sotRNAs and their cognate genes. These sotRNAs allowed us to find a UUCA tetraloop motif that is present in other archaea (ncRNA family HgcC) and in a H. salinarum intergenic ncRNA derived from a palindrome associated transposable elements (PATE). Overexpression of one sotRNA and the PATE-derived RNA harboring the tetraloop motif improved H. salinarum growth, indicating that these ncRNAs are functional.

Genome-wide CRISPR-Cas9 Screens Reveal Loss of Redundancy between PKMYT1 and WEE1 in Glioblastoma Stem-like Cells

To identify therapeutic targets for glioblastoma (GBM), we performed genome-wide CRISPR-Cas9 knockout (KO) screens in patient-derived GBM stem-like cells (GSCs) and human neural stem/progenitors (NSCs), non-neoplastic stem cell controls, for genes required for their in vitro growth. Surprisingly, the vast majority GSC-lethal hits were found outside of molecular networks commonly altered in GBM and GSCs (e.g., oncogenic drivers). In vitro and in vivo validation of GSC-specific targets revealed several strong hits, including the wee1-like kinase, PKMYT1/Myt1. Mechanistic studies demonstrated that PKMYT1 acts redundantly with WEE1 to inhibit cyclin B-CDK1 activity via CDK1-Y15 phosphorylation and to promote timely completion of mitosis in NSCs. However, in GSCs, this redundancy is lost, most likely as a result of oncogenic signaling, causing GBM-specific lethality.

A refined genome-scale reconstruction of Chlamydomonas metabolism provides a platform for systems-level analyses

Microalgae have reemerged as organisms of prime biotechnological interest due to their ability to synthesize a suite of valuable chemicals. To harness the capabilities of these organisms, we need a comprehensive systems-level understanding of their metabolism, which can be fundamentally achieved through large-scale mechanistic models of metabolism. In this study, we present a revised and significantly improved genome-scale metabolic model for the widely-studied microalga, Chlamydomonas reinhardtii. The model, iCre1355, represents a major advance over previous models, both in content and predictive power. iCre1355 encompasses a broad range of metabolic functions encoded across the nuclear, chloroplast and mitochondrial genomes accounting for 1355 genes (1460 transcripts), 2394 and 1133 metabolites. We found improved performance over the previous metabolic model based on comparisons of predictive accuracy across 306 phenotypes (from 81 mutants), lipid yield analysis and growth rates derived from chemostat-grown cells (under three conditions). Measurement of macronutrient uptake revealed carbon and phosphate to be good predictors of growth rate, while nitrogen consumption appeared to be in excess. We analyzed high-resolution time series transcriptomics data using iCre1355 to uncover dynamic pathway-level changes that occur in response to nitrogen starvation and changes in light intensity. This approach enabled accurate prediction of growth rates, the cessation of growth and accumulation of triacylglycerols during nitrogen starvation, and the temporal response of different growth-associated pathways to increased light intensity. Thus, iCre1355 represents an experimentally validated genome-scale reconstruction of C. reinhardtii metabolism that should serve as a useful resource for studying the metabolic processes of this and related microalgae.

cMonkey2: Automated, systematic, integrated detection of co-regulated gene modules for any organism

The cMonkey integrated biclustering algorithm identifies conditionally co-regulated modules of genes (biclusters). cMonkey integrates various orthogonal pieces of information which support evidence of gene co-regulation, and optimizes biclusters to be supported simultaneously by one or more of these prior constraints. The algorithm served as the cornerstone for constructing the first global, predictive Environmental Gene Regulatory Influence Network (EGRIN) model for a free-living cell, and has now been applied to many more organisms. However, due to its computational inefficiencies, long run-time and complexity of various input data types, cMonkey was not readily usable by the wider community. To address these primary concerns, we have significantly updated the cMonkey algorithm and refactored its implementation, improving its usability and extendibility. These improvements provide a fully functioning and user-friendly platform for building co-regulated gene modules and the tools necessary for their exploration and interpretation. We show, via three separate analyses of data for E. coli, M. tuberculosis and H. sapiens, that the updated algorithm and inclusion of novel scoring functions for new data types (e.g. ChIP-seq and transcription factor over-expression [TFOE]) improve discovery of biologically informative co-regulated modules. The complete cMonkey2 software package, including source code, is available at https://github.com/baliga-lab/cmonkey2.

Mapping and manipulating the Mycobacterium tuberculosis transcriptome using a transcription factor overexpression-derived regulatory network

BACKGROUND

Mycobacterium tuberculosis senses and responds to the shifting and hostile landscape of the host. To characterize the underlying intertwined gene regulatory network governed by approximately 200 transcription factors of M. tuberculosis, we have assayed the global transcriptional consequences of overexpressing each transcription factor from an inducible promoter.

RESULTS

We cloned and overexpressed 206 transcription factors in M. tuberculosis to identify the regulatory signature of each. We identified 9,335 regulatory consequences of overexpressing each of 183 transcription factors, providing evidence of regulation for 70% of the M. tuberculosis genome. These transcriptional signatures agree well with previously described M. tuberculosis regulons. The number of genes differentially regulated by transcription factor overexpression varied from hundreds of genes to none, with the majority of expression changes repressing basal transcription. Exploring the global transcriptional maps of transcription factor overexpressing (TFOE) strains, we predicted and validated the phenotype of a regulator that reduces susceptibility to a first line anti-tubercular drug, isoniazid. We also combined the TFOE data with an existing model of M. tuberculosis metabolism to predict the growth rates of individual TFOE strains with high fidelity.

CONCLUSION

This work has led to a systems-level framework describing the transcriptome of a devastating bacterial pathogen, characterized the transcriptional influence of nearly all individual transcription factors in M. tuberculosis, and demonstrated the utility of this resource. These results will stimulate additional systems-level and hypothesis-driven efforts to understand M. tuberculosis adaptations that promote disease.