Environment Tunes Propagation of Cell-to-Cell Variation in the Human Macrophage Gene Network.

Cell-to-cell variation in gene expression and the propagation of such variation (PoV or "noise propagation") from one gene to another in the gene network, as reflected by gene-gene correlation across single cells, are commonly observed in single-cell transcriptomic studies and can shape the phenotypic diversity of cell populations. While gene network "rewiring" is known to accompany cellular adaptation to different environments, how PoV changes between environments and its underlying regulatory mechanisms are less understood. Here, we systematically explored context-dependent PoV among genes in human macrophages, utilizing different cytokines as natural perturbations of multiple molecular parameters that may influence PoV. Our single-cell, epigenomic, computational, and stochastic simulation analyses reveal that environmental adaptation can tune PoV to potentially shape cellular heterogeneity by changing parameters such as the degree of phosphorylation and transcription factor-chromatin interactions. This quantitative tuning of PoV may be a widespread, yet underexplored, property of cellular adaptation to distinct environments.

Empowered genome community: leveraging a bioinformatics platform as a citizen-scientist collaboration tool.

There is on-going effort in the biomedical research community to leverage Next Generation Sequencing (NGS) technology to identify genetic variants that affect our health. The main challenge facing researchers is getting enough samples from individuals either sick or healthy - to be able to reliably identify the few variants that are causal for a phenotype among all other variants typically seen among individuals. At the same time, more and more individuals are having their genome sequenced either out of curiosity or to identify the cause of an illness. These individuals may benefit from of a way to view and understand their data. QIAGEN's Ingenuity Variant Analysis is an online application that allows users with and without extensive bioinformatics training to incorporate information from published experiments, genetic databases, and a variety of statistical models to identify variants, from a long list of candidates, that are most likely causal for a phenotype as well as annotate variants with what is already known about them in the literature and databases. Ingenuity Variant Analysis is also an information sharing platform where users may exchange samples and analyses. The Empowered Genome Community (EGC) is a new program in which QIAGEN is making this on-line tool freely available to any individual who wishes to analyze their own genetic sequence. EGC members are then able to make their data available to other Ingenuity Variant Analysis users to be used in research. Here we present and describe the Empowered Genome Community in detail. We also present a preliminary, proof-of-concept study that utilizes the 200 genomes currently available through the EGC. The goal of this program is to allow individuals to access and understand their own data as well as facilitate citizen-scientist collaborations that can drive research forward and spur quality scientific dialogue in the general public.

Predictive computational modeling of the mucosal immune responses during Helicobacter pylori infection.

T helper (Th) cells play a major role in the immune response and pathology at the gastric mucosa during Helicobacter pylori infection. There is a limited mechanistic understanding regarding the contributions of CD4+ T cell subsets to gastritis development during H. pylori colonization. We used two computational approaches: ordinary differential equation (ODE)-based and agent-based modeling (ABM) to study the mechanisms underlying cellular immune responses to H. pylori and how CD4+ T cell subsets influenced initiation, progression and outcome of disease. To calibrate the model, in vivo experimentation was performed by infecting C57BL/6 mice intragastrically with H. pylori and assaying immune cell subsets in the stomach and gastric lymph nodes (GLN) on days 0, 7, 14, 30 and 60 post-infection. Our computational model reproduced the dynamics of effector and regulatory pathways in the gastric lamina propria (LP) in silico. Simulation results show the induction of a Th17 response and a dominant Th1 response, together with a regulatory response characterized by high levels of mucosal Treg) cells. We also investigated the potential role of peroxisome proliferator-activated receptor γ (PPARγ) activation on the modulation of host responses to H. pylori by using loss-of-function approaches. Specifically, in silico results showed a predominance of Th1 and Th17 cells in the stomach of the cell-specific PPARγ knockout system when compared to the wild-type simulation. Spatio-temporal, object-oriented ABM approaches suggested similar dynamics in induction of host responses showing analogous T cell distributions to ODE modeling and facilitated tracking lesion formation. In addition, sensitivity analysis predicted a crucial contribution of Th1 and Th17 effector responses as mediators of histopathological changes in the gastric mucosa during chronic stages of infection, which were experimentally validated in mice. These integrated immunoinformatics approaches characterized the induction of mucosal effector and regulatory pathways controlled by PPARγ during H. pylori infection affecting disease outcomes.

Systems modeling of molecular mechanisms controlling cytokine-driven CD4+ T cell differentiation and phenotype plasticity.

Differentiation of CD4+ T cells into effector or regulatory phenotypes is tightly controlled by the cytokine milieu, complex intracellular signaling networks and numerous transcriptional regulators. We combined experimental approaches and computational modeling to investigate the mechanisms controlling differentiation and plasticity of CD4+ T cells in the gut of mice. Our computational model encompasses the major intracellular pathways involved in CD4+ T cell differentiation into T helper 1 (Th1), Th2, Th17 and induced regulatory T cells (iTreg). Our modeling efforts predicted a critical role for peroxisome proliferator-activated receptor gamma (PPARγ) in modulating plasticity between Th17 and iTreg cells. PPARγ regulates differentiation, activation and cytokine production, thereby controlling the induction of effector and regulatory responses, and is a promising therapeutic target for dysregulated immune responses and inflammation. Our modeling efforts predict that following PPARγ activation, Th17 cells undergo phenotype switch and become iTreg cells. This prediction was validated by results of adoptive transfer studies showing an increase of colonic iTreg and a decrease of Th17 cells in the gut mucosa of mice with colitis following pharmacological activation of PPARγ. Deletion of PPARγ in CD4+ T cells impaired mucosal iTreg and enhanced colitogenic Th17 responses in mice with CD4+ T cell-induced colitis. Thus, for the first time we provide novel molecular evidence in vivo demonstrating that PPARγ in addition to regulating CD4+ T cell differentiation also plays a major role controlling Th17 and iTreg plasticity in the gut mucosa.

ENteric Immunity SImulator: a tool for in silico study of gastroenteric infections.

Clinical symptoms of microbial infection of the gastrointestinal (GI) tract are often exacerbated by inflammation induced pathology. Identifying novel avenues for treating and preventing such pathologies is necessary and complicated by the complexity of interacting immune pathways in the gut, where effector and inflammatory immune cells are regulated by anti-inflammatory or regulatory cells. Here we present new advances in the development of the ENteric Immunity SImulator (ENISI), a simulator of GI immune mechanisms in response to resident commensal bacteria as well as invading pathogens and the effect on the development of intestinal lesions. ENISI is a tool for identifying potential treatment strategies that reduce inflammation-induced damage and, at the same time, ensure pathogen removal by allowing one to test plausibility of in vitro observed behavior as explanations for observations in vivo, propose behaviors not yet tested in vitro that could explain these tissue-level observations, and conduct low-cost, preliminary experiments of proposed interventions/treatments. An example of such application is shown in which we simulate dysentery resulting from Brachyispira hyodysenteriae infection and identify aspects of the host immune pathways that lead to continued inflammation-induced tissue damage even after pathogen elimination.

Model of colonic inflammation: immune modulatory mechanisms in inflammatory bowel disease.

Inflammatory bowel disease (IBD) is an immunoinflammatory illness of the gut initiated by an immune response to bacteria in the microflora. The resulting immunopathogenesis leads to lesions in epithelial lining of the colon through which bacteria may infiltrate the tissue causing recurring bouts of diarrhea, rectal bleeding, and malnutrition. In healthy individuals such immunopathogenesis is avoided by the presence of regulatory cells that inhibit the inflammatory pathway. Highly relevant to the search for treatment strategies is the identification of components of the inflammatory pathway that allow regulatory mechanisms to be overridden and immunopathogenesis to proceed. In vitro techniques have identified cellular interactions involved in inflammation-regulation crosstalk. However, tracing immunological mechanisms discovered at the cellular level confidently back to an in vivo context of multiple, simultaneous interactions has met limited success. To explore the impact of specific interactions, we have constructed a system of 29 ordinary differential equations representing different phenotypes of T-cells, macrophages, dendritic cells, and epithelial cells as they move and interact with bacteria in the lumen, lamina propria, and lymphoid tissue of the colon. Simulations revealed the positive inflammatory feedback loop formed by inflammatory M1 macrophage activation of T-cells as a driving force underlying the immunopathology of IBD. Furthermore, strategies that remove M1 from the site of infection, by either (i) increasing its potential to switch to a regulatory M2 phenotype or (ii) increasing the rate of reversion (for M1 and M2 alike) to a resting state, cease immunopathogenesis even as bacteria are eliminated by other inflammatory cells. Based on these results, we identify macrophages and their mechanisms of plasticity as key targets for mucosal inflammation intervention strategies. In addition, we propose that the primary mechanism behind the association of PPARgamma mutation with IBD is its ability to mediate the M1 to M2 switch.

An analysis of enzyme kinetics data for mitochondrial DNA strand termination by nucleoside reverse transcription inhibitors.

Nucleoside analogs used in antiretroviral treatment have been associated with mitochondrial toxicity. The polymerase-gamma hypothesis states that this toxicity stems from the analogs' inhibition of the mitochondrial DNA polymerase (polymerase-gamma) leading to mitochondrial DNA (mtDNA) depletion. We have constructed a computational model of the interaction of polymerase-gamma with activated nucleoside and nucleotide analog drugs, based on experimentally measured reaction rates and base excision rates, together with the mtDNA genome size, the human mtDNA sequence, and mitochondrial dNTP concentrations. The model predicts an approximately 1000-fold difference in the activated drug concentration required for a 50% probability of mtDNA strand termination between the activated di-deoxy analogs d4T, ddC, and ddI (activated to ddA) and the activated forms of the analogs 3TC, TDF, AZT, FTC, and ABC. These predictions are supported by experimental and clinical data showing significantly greater mtDNA depletion in cell culture and patient samples caused by the di-deoxy analog drugs. For zidovudine (AZT) we calculated a very low mtDNA replication termination probability, in contrast to its reported mitochondrial toxicity in vitro and clinically. Therefore AZT mitochondrial toxicity is likely due to a mechanism that does not involve strand termination of mtDNA replication.