A Method for Inferring Gene Regulatory Networks from Time

Total Page:16

File Type:pdf, Size:1020Kb

A Method for Inferring Gene Regulatory Networks from Time www.nature.com/scientificreports OPEN BGRMI: A method for inferring gene regulatory networks from time-course gene expression data Received: 08 September 2016 Accepted: 24 October 2016 and its application in breast cancer Published: 23 November 2016 research Luis F. Iglesias-Martinez1, Walter Kolch1,2,3 & Tapesh Santra1 Reconstructing gene regulatory networks (GRNs) from gene expression data is a challenging problem. Existing GRN reconstruction algorithms can be broadly divided into model-free and model–based methods. Typically, model-free methods have high accuracy but are computation intensive whereas model-based methods are fast but less accurate. We propose Bayesian Gene Regulation Model Inference (BGRMI), a model-based method for inferring GRNs from time-course gene expression data. BGRMI uses a Bayesian framework to calculate the probability of different models of GRNs and a heuristic search strategy to scan the model space efficiently. Using benchmark datasets, we show that BGRMI has higher/comparable accuracy at a fraction of the computational cost of competing algorithms. Additionally, it can incorporate prior knowledge of potential gene regulation mechanisms and TF hetero-dimerization processes in the GRN reconstruction process. We incorporated existing ChIP-seq data and known protein interactions between TFs in BGRMI as sources of prior knowledge to reconstruct transcription regulatory networks of proliferating and differentiating breast cancer (BC) cells from time-course gene expression data. The reconstructed networks revealed key driver genes of proliferation and differentiation in BC cells. Some of these genes were not previously studied in the context of BC, but may have clinical relevance in BC treatment. Cellular functions depend on the precise regulation of thousands of genes which are activated or silenced by transcription factors (TFs)1. The networks representing the interactions between TFs and their target genes are typically known as GRNs and can be reconstructed from temporal measurements of gene expressions2–5. The GRN reconstruction methods can be classified into two main categories; model-based and model-free methods. Model-based methods aim to capture the regulatory interactions by fitting mathematical models of gene regula- tion to observed gene expression data2,4,5. On the other hand, model-free approaches use information-theoretic criteria to infer the structure of the network3,6. Though the performances of GRN reconstruction methods depend on several aspects such as data type, network properties of the GRN etc7., generally, model-based methods tend to be faster but have lower predictive performance than model-free methods 20156. However, model-free methods are often not scalable enough to reconstruct genome-wide GRNs5,6 in reasonable time. Typically, model-based methods formulate the expression of a gene as a function of its regulators, evaluate competing models containing different sets of regulators and chose those which closely predict the target gene expression4,5,8. Although vast majority of model based methods assume that the expressions of a gene and its regulators are linearly depend- ent5,9–12, these methods use different model search algorithms e.g. Least Absolute Shrinkage and Selection Operator (LASSO), Dantzig Selector, elastic net, Markov Chain Monte Carlo and Heuristic search5,12–19. Some of these methods such as LASSO and elastic net choose the best model, others such as MCMC or Heuristic search based Bayesian Model Averaging (BMA) methods5,12,13,15,18 select multiple models that provide close fits to the 1Systems Biology Ireland, University College Dublin, Belfield, Dublin 4, Republic of Ireland.2 Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Belfield, Dublin 4, Ireland.3 School of Medicine and Medical Science, University College Dublin, Belfield, Dublin 4, Ireland. Correspondence and requests for materials should be addressed to T.S. (email: [email protected]) SCIENTIFIC REPORTS | 6:37140 | DOI: 10.1038/srep37140 1 www.nature.com/scientificreports/ data and use these to estimate an average model along with its confidence interval. It is also possible to incorpo- rate different types of existing data in BMA5,12,18 to increase the accuracy of the reconstructed GRN. We developed BGRMI, a model-based method that relies on the principles of BMA for inferring GRNs from time course gene expression data. BGRMI uses discretized ordinary differential equation (DODE) based mathe- matical models to formulate the interactions between each gene and its regulators. It formulates the rate of change in a gene’s expression as a function of the expressions of its regulators, takes basal expression and self-regulation into account and therefore provides a more realistic model of gene regulation than many existing methods. These models are then used in a Bayesian framework to evaluate how likely a set of TFs is to regulate a certain gene. We developed a greedy heuristic search algorithm to explore different combinations of TFs and find the most likely TF combinations for each gene. The proposed algorithm is faster and more scalable than many existing methods. The average of some of the most likely models was then used to represent the regulatory model of the gene. We compared the accuracy of BGRMI against other methods using in-silico and in-vivo benchmark datasets. BGRMI consistently out-performed most of the other competing methods in our benchmarking study. We then showed how additional data sources, e.g. ChIP-seq and the protein-protein interaction (PPI) between TFs can be incorpo- rated as prior knowledge in the core BGRMI formulation. Finally, we applied BGRMI to study the transcriptional mechanisms that lead to proliferation and differentiation in BC cells by combining ChIP-seq, PPI and time course gene expression profiles. Our study uncovered previously unknown transcriptional mechanisms that drive phe- notypic changes in BC cells. Method A brief overview of the BGRMI algorithm is as follows. We first developed a mathematical model of TF-mediated gene regulations. The model can predict the temporal changes in the expressions of target genes using the tempo- ral expression patterns of TFs as input. This model is then used to evaluate different combinations of TFs to find those that closely predict target gene expressions. This is done by iteratively exploring different TF combinations, calculating the posterior probability of each of these combinations to predict target gene expression, and selecting those with high probabilities. The average of selected gene regulation models for each gene is used as its regulation model. Below we describe each step of our algorithm in detail. The mathematical model of gene regulation. We used a discretized form of ODEs to formulate the dependence of a target gene on its regulators as shown below. [(mRNA tm)(−−RN At ∆t)] mRNA ()tt−∆ ∆=ii=+αεβTi + mRNAi ()t ii i ()t ∆t TFi()tt−∆ (1) Here, mRNAi(t) is the expression of gene i at time t, αi is the basal gene expression rate, βi is a vector that con- tains the coefficients of self regulation (by means of degradation, auto-activation/inhibition) and the regulation by a set of TFs (TFi), TFi(t − Δt ) are the expressions of the TFs that regulate gene i at time (t − Δt ). εi(t) is the model fitting error caused by the measurement noise in expression data. Since measurement noise is random εi(t) is a 2 2 12,18 random variable and typically has Gaussian distribution with zero mean and variance σ , i.e. εi(t)~N(0, σ ) . The error variance (σ2) depends on many factors such as biological variability and measurement noise, and is typically unknown. The posterior probability of a gene regulation model. We used Bayesian statistics to calculate this probability. There are two main components of Bayesian formulations; (a) the prior probability (p(Mk)) which represents how well a model (Mk) is supported by prior knowledge, and (b) the likelihood function (p(mR- 20 NAi|Mk)) which evaluates how well a model explains experimental data. By Bayes’ rule , the posterior probability (p(Mk|mRNAi)) is proportional to the product of these two entities and represents how well a model (Mk) is sup- ported by prior knowledge and experimental data combined. In the absence of prior knowledge, we assumed that sparse regulatory models (i.e. those involving fewer TFs) are a priori more likely than dense models (i.e. those involving a large number of TFs). This assumption was for- 5 −2.66 mulated by assigning the following prior distributions over regulatory models (Mk) : p(Mk) = L , where L is the number of regulators in the model. However, there is a wealth of publicly available information about GRNs of several organisms such as yeast, E. coli, humans etc. This information can be used to formulate more informative priors for reconstructing GRNs of these organisms. We shall discuss the formulation of priors for human GRNs using ChIP-seq data12 in a later section where we describe the implementation of BGRMI on human transcrip- tomic data. The likelihood of a gene regulation model (Mk) is the probability that the observed expression pattern 5,12 (mRNAi) of gene i, can be predicted by the model (Mk), and has the following form : T 22TimRNA ()tt−∆ PN(mRNA TF ,,,ασββ) = , σ iii i ∏ i −∆ tt=∆ TFi()tt (2) 2 The likelihood function in Eq. 2 depends on the model parameters (αi, βi σ ), whose values are typically unknown. Therefore, we evaluated the average of the likelihood (Eq. 2) over all possible values of the model parameters. The average likelihood is also called the marginal likelihood P( (mRNAi|Mk)). To analytically calculate the marginal likelihood we assigned conjugate prior distributions to each of the unknown variables. These distri- butions represent our prior knowledge of how likely a parameter is to have a certain value. Following Fernandez 13 13,21 et al. we assigned uninformative Jeffrey’s prior distribution to the basal expression rates αi, p(αi) = 1, which SCIENTIFIC REPORTS | 6:37140 | DOI: 10.1038/srep37140 2 www.nature.com/scientificreports/ implies that αi is equally likely to have any real value.
Recommended publications
  • Screening and Identification of Key Biomarkers in Clear Cell Renal Cell Carcinoma Based on Bioinformatics Analysis
    bioRxiv preprint doi: https://doi.org/10.1101/2020.12.21.423889; this version posted December 23, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Screening and identification of key biomarkers in clear cell renal cell carcinoma based on bioinformatics analysis Basavaraj Vastrad1, Chanabasayya Vastrad*2 , Iranna Kotturshetti 1. Department of Biochemistry, Basaveshwar College of Pharmacy, Gadag, Karnataka 582103, India. 2. Biostatistics and Bioinformatics, Chanabasava Nilaya, Bharthinagar, Dharwad 580001, Karanataka, India. 3. Department of Ayurveda, Rajiv Gandhi Education Society`s Ayurvedic Medical College, Ron, Karnataka 562209, India. * Chanabasayya Vastrad [email protected] Ph: +919480073398 Chanabasava Nilaya, Bharthinagar, Dharwad 580001 , Karanataka, India bioRxiv preprint doi: https://doi.org/10.1101/2020.12.21.423889; this version posted December 23, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Abstract Clear cell renal cell carcinoma (ccRCC) is one of the most common types of malignancy of the urinary system. The pathogenesis and effective diagnosis of ccRCC have become popular topics for research in the previous decade. In the current study, an integrated bioinformatics analysis was performed to identify core genes associated in ccRCC. An expression dataset (GSE105261) was downloaded from the Gene Expression Omnibus database, and included 26 ccRCC and 9 normal kideny samples. Assessment of the microarray dataset led to the recognition of differentially expressed genes (DEGs), which was subsequently used for pathway and gene ontology (GO) enrichment analysis.
    [Show full text]
  • Figure S1. Representative Report Generated by the Ion Torrent System Server for Each of the KCC71 Panel Analysis and Pcafusion Analysis
    Figure S1. Representative report generated by the Ion Torrent system server for each of the KCC71 panel analysis and PCaFusion analysis. (A) Details of the run summary report followed by the alignment summary report for the KCC71 panel analysis sequencing. (B) Details of the run summary report for the PCaFusion panel analysis. A Figure S1. Continued. Representative report generated by the Ion Torrent system server for each of the KCC71 panel analysis and PCaFusion analysis. (A) Details of the run summary report followed by the alignment summary report for the KCC71 panel analysis sequencing. (B) Details of the run summary report for the PCaFusion panel analysis. B Figure S2. Comparative analysis of the variant frequency found by the KCC71 panel and calculated from publicly available cBioPortal datasets. For each of the 71 genes in the KCC71 panel, the frequency of variants was calculated as the variant number found in the examined cases. Datasets marked with different colors and sample numbers of prostate cancer are presented in the upper right. *Significantly high in the present study. Figure S3. Seven subnetworks extracted from each of seven public prostate cancer gene networks in TCNG (Table SVI). Blue dots represent genes that include initial seed genes (parent nodes), and parent‑child and child‑grandchild genes in the network. Graphical representation of node‑to‑node associations and subnetwork structures that differed among and were unique to each of the seven subnetworks. TCNG, The Cancer Network Galaxy. Figure S4. REVIGO tree map showing the predicted biological processes of prostate cancer in the Japanese. Each rectangle represents a biological function in terms of a Gene Ontology (GO) term, with the size adjusted to represent the P‑value of the GO term in the underlying GO term database.
    [Show full text]
  • An Amino-Terminal Domain of Mxil Mediates Anti-Myc Oncogenic Activity and Interacts with a Homolog of the Yeast Transcriptional Repressor SIN3
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Elsevier - Publisher Connector Cell, Vol, 80, 777-786, March 10, 1995, Copyright © 1995 by Cell Press An Amino-Terminal Domain of Mxil Mediates Anti-Myc Oncogenic Activity and Interacts with a Homolog of the Yeast Transcriptional Repressor SIN3 Nicole Schreiber-Agus,*t Lynda Chin,*tt Ken Chen,t et al., 1990), and a carboxy-terminal a-helical domain re- Richard Torres, t Govinda Rao,§ Peter Guida,t quired for dimerization with another basic region-helix- Arthur h Skoultchi,§ and Ronald A. DePinhot Ioop-helix-leucine zipper (bHLH-LZ) protein, Max (Black- "rDepartments of Microbiology and Immunology wood and Eisenman, 1991; Prendergast et al., 1991). and of Medicine Many of the biochemical and biological activities of Myc §Department of Cell Biology appear to be highly dependent upon its association with ~Division of Dermatology Max (Blackwood and Eisenman, 1991 ; Prendergast et al., Albert Einstein College of Medicine 1991; Kretzner et al., 1992; Amati et al., 1993a, 1993b). Bronx, New York 10461 In addition to its key role as an obligate partner in transacti- vation-competent Myc-Max complexes, Max may also re- press Myc-responsive genes through the formation of Summary transactivation-inert complexes that are capable of bind- ing the Myc-Max recognition sequence (Blackwood et al., Documented interactions among members of the Myc 1992; Kato et al., 1992; Kretzner et al., 1992; Makela et superfamily support a yin-yang model for the regula- al., 1992; Mukherjee et al., 1992; Prendergast et al., 1992; tion of Myc-responsive genes in which t ransactivation- Ayer et al., 1993; Zervos et al., 1993).
    [Show full text]
  • Transcriptional Control of Tissue-Resident Memory T Cell Generation
    Transcriptional control of tissue-resident memory T cell generation Filip Cvetkovski Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the Graduate School of Arts and Sciences COLUMBIA UNIVERSITY 2019 © 2019 Filip Cvetkovski All rights reserved ABSTRACT Transcriptional control of tissue-resident memory T cell generation Filip Cvetkovski Tissue-resident memory T cells (TRM) are a non-circulating subset of memory that are maintained at sites of pathogen entry and mediate optimal protection against reinfection. Lung TRM can be generated in response to respiratory infection or vaccination, however, the molecular pathways involved in CD4+TRM establishment have not been defined. Here, we performed transcriptional profiling of influenza-specific lung CD4+TRM following influenza infection to identify pathways implicated in CD4+TRM generation and homeostasis. Lung CD4+TRM displayed a unique transcriptional profile distinct from spleen memory, including up-regulation of a gene network induced by the transcription factor IRF4, a known regulator of effector T cell differentiation. In addition, the gene expression profile of lung CD4+TRM was enriched in gene sets previously described in tissue-resident regulatory T cells. Up-regulation of immunomodulatory molecules such as CTLA-4, PD-1, and ICOS, suggested a potential regulatory role for CD4+TRM in tissues. Using loss-of-function genetic experiments in mice, we demonstrate that IRF4 is required for the generation of lung-localized pathogen-specific effector CD4+T cells during acute influenza infection. Influenza-specific IRF4−/− T cells failed to fully express CD44, and maintained high levels of CD62L compared to wild type, suggesting a defect in complete differentiation into lung-tropic effector T cells.
    [Show full text]
  • Dissertation
    Dissertation Submitted to the Combined Faculties for the Natural Sciences and for Mathematics of the Ruperto-Carola University of Heidelberg, Germany for the Degree of Doctor of Natural Sciences Presented by Ann-Cathrin Hofer (M.Sc.) Born in Heidelberg, Germany Oral Examination: 12th of September 2016 Regulatory T cells protect the neonatal liver and secure the hepatic circadian rhythm Referees 1st Referee: Prof. Dr. Peter Angel 2nd Referee: Dr. Markus Feuerer This dissertation was performed and written during the period from November 2012 to May 2016 in the German Cancer Research Center (DKFZ) under the supervision of Prof. Dr. Peter Angel and direct supervision of Dr. Markus Feuerer. The dissertation was submitted to the Combined Faculties for the Natural Sciences and for Mathematics of the Ruperto-Carola University of Heidelberg, Germany in June 2016. German Cancer Research Center (DKFZ) Immune Tolerance (D100) Im Neuenheimer Feld 280 69120 Heidelberg, Germany I II Confirmation Hereby, I confirm that I have written this thesis independently, using only the results of my investigation unless otherwise stated. Furthermore, I declare that I have not submitted this thesis for a degree to any other academic or similar institution. Parts of this dissertation have been submitted for publishing: Regulatory T cells protect the liver early in life and safeguard the hepatic circadian rhythm Ann-Cathrin Hofer, Thomas Hielscher, David M. Richards, Michael Delacher, Ulrike Träger, Sophia Föhr, Artyom Vlasov, Marvin Wäsch, Marieke Esser, Annette Kopp-Schneider, Achim Breiling, Frank Lyko, Ursula Klingmüller, Peter Angel, Jakub Abramson, Jeroen Krijgsveld & Markus Feuerer Parts of the experiments in this dissertation were performed in collaboration with other research groups as follows: CG methylation analysis with the 454 pyrosequencing technology: Division of Epigenetics, DKFZ, Heidelberg Dr.
    [Show full text]
  • The Id-Protein Family in Developmental and Cancer-Associated Pathways Cornelia Roschger and Chiara Cabrele*
    Roschger and Cabrele Cell Communication and Signaling (2017) 15:7 DOI 10.1186/s12964-016-0161-y REVIEW Open Access The Id-protein family in developmental and cancer-associated pathways Cornelia Roschger and Chiara Cabrele* Abstract Inhibitors of DNA binding and cell differentiation (Id) proteins are members of the large family of the helix-loop- helix (HLH) transcription factors, but they lack any DNA-binding motif. During development, the Id proteins play a key role in the regulation of cell-cycle progression and cell differentiation by modulating different cell-cycle regulators both by direct and indirect mechanisms. Several Id-protein interacting partners have been identified thus far, which belong to structurally and functionally unrelated families, including, among others, the class I and II bHLH transcription factors, the retinoblastoma protein and related pocket proteins, the paired-box transcription factors, and the S5a subunit of the 26 S proteasome. Although the HLH domain of the Id proteins is involved in most of their protein-protein interaction events, additional motifs located in their N-terminal and C-terminal regions are required for the recognition of diverse protein partners. The ability of the Id proteins to interact with structurally different proteins is likely to arise from their conformational flexibility: indeed, these proteins contain intrinsically disordered regions that, in the case of the HLH region, undergo folding upon self- or heteroassociation. Besides their crucial role for cell-fate determination and cell-cycle progression during development, other important cellular events have been related to the Id-protein expression in a number of pathologies. Dysregulated Id-protein expression has been associated with tumor growth, vascularization, invasiveness, metastasis, chemoresistance and stemness, as well as with various developmental defects and diseases.
    [Show full text]
  • Hypoxia and Aging Eui-Ju Yeo1
    Yeo Experimental & Molecular Medicine (2019) 51:67 https://doi.org/10.1038/s12276-019-0233-3 Experimental & Molecular Medicine REVIEW ARTICLE Open Access Hypoxia and aging Eui-Ju Yeo1 Abstract Eukaryotic cells require sufficient oxygen (O2) for biological activity and survival. When the oxygen demand exceeds its supply, the oxygen levels in local tissues or the whole body decrease (termed hypoxia), leading to a metabolic crisis, threatening physiological functions and viability. Therefore, eukaryotes have developed an efficient and rapid oxygen sensing system: hypoxia-inducible factors (HIFs). The hypoxic responses are controlled by HIFs, which induce the expression of several adaptive genes to increase the oxygen supply and support anaerobic ATP generation in eukaryotic cells. Hypoxia also contributes to a functional decline during the aging process. In this review, we focus on the molecular mechanisms regulating HIF-1α and aging-associated signaling proteins, such as sirtuins, AMP-activated protein kinase, mechanistic target of rapamycin complex 1, UNC-51-like kinase 1, and nuclear factor κB, and their roles in aging and aging-related diseases. In addition, the effects of prenatal hypoxia and obstructive sleep apnea (OSA)- induced intermittent hypoxia have been reviewed due to their involvement in the progression and severity of many diseases, including cancer and other aging-related diseases. The pathophysiological consequences and clinical manifestations of prenatal hypoxia and OSA-induced chronic intermittent hypoxia are discussed in detail. Introduction chemoreceptors stimulates the neurotransmitter release 3 1234567890():,; 1234567890():,; 1234567890():,; 1234567890():,; Oxygen (O2) plays critical roles in aerobic respiration pathway and modulates the activity of a neutral endo- and metabolism as the final electron acceptor of the peptidase, neprilysin (NEP), which modifies the cellular mitochondrial electron transport chain, which is respon- response to hypoxia by hydrolyzing substance P4.
    [Show full text]
  • 1.1.4.1. Tumour Suppressor Genes
    To my Mother and Father Learn from yesterday, live for today, hope for tomorrow. The important thing is not to stop questioning. Albert Einstein Inhibition of Tumourigenicity of Small Cell Lung Cancer by Simultaneous Suppression of Id1 and Id3 Expression Danqing Chen ABSTRACT Inhibitor of DNA binding (Id) proteins are a group of transcription factors belonging to the basic helix-loop-helix (bHLH) family and play a wide range of roles in differentiation, proliferation and cell cycle progression. Id proteins act as negative dominant regulators of other bHLH factors by making dimers to these factors to prevent them from binding to E-box of DNA and, hence, to inhibit transcription of target genes. In this work, we first established SCLC cell line N417-derived sublines expressing reduced levels of Id1 and Id3 by transfection of a single vector constructed to co-express two shRNAs simultaneously. Then we investigated the effect of either singly or jointly suppressed Id1 or Id3 on tumourigenicity of SCLC cells in vitro and in vivo. The molecular mechanisms involved in the functional roles of Id1 and Id3 were also assessed. Id1-suppressed cells and Id1 and Id3 double knockdown cells produced significant reductions in proliferation rate by more than 1.4- and 3.9-fold respectively when compared with the control. Soft agar assay showed the number of colonies produced by Id1-suppressed cells and Id1 and Id3 double knockdown cells were reduced by more than 13.7- and 233-fold respectively compared with the control. The suppression effect was also observed in the invasion assay which showed that Id1-suppressed cells and Id1 and Id3 double knockdown cells produced more than 1.7- and 4.6- fold reduction respectively in relative invasiveness.
    [Show full text]
  • Ltrs Activated by Epstein-Barr Virus-Induced Transformation of B Cells Alter the Transcriptome
    Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press LTRs activated by Epstein-Barr virus-induced transformation of B cells alter the transcriptome Authors: Amy Leung1,*, Candi Trac1, Hiroyuki Kato1, Kevin R Costello1,2, Zhaoxia Chen1, Rama Natarajan1,2, Dustin E Schones1,2,* 1Department of Diabetes Complications and Metabolism, Beckman Research Institute, City of Hope, Duarte, CA, USA 2Irell & Manella Graduate School of Biological Sciences, City of Hope, Duarte, CA, USA *Corresponding Authors: Amy Leung ([email protected]), Dustin E. Schones ([email protected]) Amy Leung, Ph.D. Department of Diabetes Complications and Metabolism City of Hope Duarte, CA, 91016 Email: [email protected] Phone: 626-218-2278 Dustin Schones, Ph.D. Department of Diabetes Complications and Metabolism City of Hope Duarte, CA, 91016 Email: [email protected] Phone: 626-218-9319 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press Running title: Activated LTRs alter the transcriptome of B cells (50 characters) Keywords: Retrotransposons, LTRs, H3K4me3, EBV, cryptic promoters Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press ABSTRACT Endogenous retroviruses (ERVs) are ancient viral elements that have accumulated in the genome through retrotransposition events. Although they have lost their ability to transpose, many of the long terminal repeats (LTRs) that originally flanked full-length ERVs maintain the ability to regulate transcription. While these elements are typically repressed in somatic cells, they can function as transcriptional enhancers and promoters when this repression is lost. Epstein-Barr virus (EBV), which transforms primary B cells into continuously proliferating cells, is a tumor virus associated with lymphomas.
    [Show full text]
  • Sox2 Promotes Malignancy in Glioblastoma by Regulating
    Volume 16 Number 3 March 2014 pp. 193–206.e25 193 www.neoplasia.com Artem D. Berezovsky*, Laila M. Poisson†, Sox2 Promotes Malignancy in ‡ ‡ David Cherba , Craig P. Webb , Glioblastoma by Regulating Andrea D. Transou*, Nancy W. Lemke*, Xin Hong*, Laura A. Hasselbach*, Susan M. Irtenkauf*, Plasticity and Astrocytic Tom Mikkelsen*,§ and Ana C. deCarvalho* Differentiation1,2 *Department of Neurosurgery, Henry Ford Hospital, Detroit, MI; †Department of Public Health Sciences, Henry Ford Hospital, Detroit, MI; ‡Program of Translational Medicine, Van Andel Research Institute, Grand Rapids, MI; §Department of Neurology, Henry Ford Hospital, Detroit, MI Abstract The high-mobility group–box transcription factor sex-determining region Y–box 2 (Sox2) is essential for the maintenance of stem cells from early development to adult tissues. Sox2 can reprogram differentiated cells into pluripotent cells in concert with other factors and is overexpressed in various cancers. In glioblastoma (GBM), Sox2 is a marker of cancer stemlike cells (CSCs) in neurosphere cultures and is associated with the proneural molecular subtype. Here, we report that Sox2 expression pattern in GBM tumors and patient-derived mouse xenografts is not restricted to a small percentage of cells and is coexpressed with various lineage markers, suggesting that its expression extends beyond CSCs to encompass more differentiated neoplastic cells across molecular subtypes. Employing a CSC derived from a patient with GBM and isogenic differentiated cell model, we show that Sox2 knockdown in the differentiated state abolished dedifferentiation and acquisition of CSC phenotype. Furthermore, Sox2 deficiency specifically impaired the astrocytic component of a biphasic gliosarcoma xenograft model while allowing the formation of tumors with sarcomatous phenotype.
    [Show full text]
  • Mir-155: a Novel Target in Allergic Asthma
    Int. J. Mol. Sci. 2016, 17, 1773; doi: 10.3390/ijms17101773 S1 of S2 Supplementary Materials: miR-155: A Novel Target in Allergic Asthma Hong Zhou, Junyao Li, Peng Gao, Qi Wang and Jie Zhang Table S1. miR-155 targets identified by TargetScan, miRTarBase and both bioinformatic analyses. Names Total Elements LCORL CHD8 TOMM20 SLC33A1 CHD9 ZNF248 IRF2BP2 DNAJB1 C10orf12 PALLD CARD11 GNAS ZBTB38 RAPH1 ETNK2 MSH6 ARL5B CCDC41 MMP16 RHEB TOMM34 MEF2A RICTOR RAB11FIP2 FAM135A ZBTB18 TMEM33 TCF12 KRAS TM6SF1 DHX40 PICALM MYO10 TCF4 FUBP1 ATP6V1C1 SERTAD2 SH3PXD2A UBQLN2 YWHAZ AGO4 CHAF1A ZNF236 MORC3 MEIS1 WWC1 TAB2 NAA50 PRKAR1A CSNK1G2 PHC2 HBP1 SPRED1 ADAM10 KANSL1 MIDN ZNF644 NFAT5 IL17RB STRN3 MAP3K10 ZSWIM6 DMTF1 ITK PDE3A ZIC3 PELI1 CSNK1A1 ARID2 GSK3B SPIN1 TSPAN14 PTAR1 FOXK1 WEE1 PKN2 TPD52 CARHSP1 MYBL1 WBP1L SAP30L VEZF1 EEF2 FLT1 PHF17 RCOR1 SMAD2 CBFB RORA HIVEP2 CHD7 RAP1B TargetScan 190 SPI1 PEA15 FGF7 RREB1 CBL MYLK S1PR1 TMEM136 PIK3CA NKX3-1 CTLA4 miRTarBase RAB3B SMAD1 ANKFY1 FOS SKIV2L2 SMARCA4 TP53INP1 TSHZ3 PSMG1 FGF2 SKI CPEB4 JARID2 MSI2 SWSAP1 LRRC40 ETS1 COPS3 IKBKE SOCS1 TRIM32 LRRC59 CDC73 RAB5C CAB39 LNX2 NSA2 CDC37 MBNL3 MAFB INPP5D E2F2 PKIA RAB30 CEP41 DET1 UBTD2 C3orf18 BACH1 RAPGEF2 CREBRF SHANK2 PAXBP1 BAG5 KBTBD2 KIF3A HHIP EHD1 HERC4 PALD1 HNRNPA3 N4BP1 PIK3R1 PTPRJ NOVA1 GPM6B CKAP5 TAPT1 CLDN1 SIRT1 SEPT11 COLGALT1 HMGCS1 TLE4 TERF1 ZNF703 FOXO3 KCTD3 APC INADL BCAT1 WNK1 CEBPB TRPS1 CSF1R KDM3A MYO1D RNF123 TADA2B AAK1 RBAK USP8 RCN2 SMAD5 PDE12 ZNF652 MYB SMNDC1 RPTOR PLCE1 KIF26B TNIK RTKN2 ZPLD1 ARRB2
    [Show full text]
  • Single Cell Transcriptomics of Regulatory T Cells Reveals Trajectories of Tissue Adaptation ​ ​ ​ ​ ​ ​
    bioRxiv preprint doi: https://doi.org/10.1101/217489; this version posted November 22, 2017. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. Single cell transcriptomics of regulatory T cells reveals trajectories of tissue adaptation ​ ​ ​ ​ ​ ​ Authors: 1,2,a 1,a 3,4 7 Ricardo J Miragaia ,​ Tomás Gomes ,​ Agnieszka Chomka ,​ Laura Jardine ,​ Angela ​ ​ ​ ​ 8 3,4,5 1,9 1 3,4,6 Riedel ,​ Ahmed N. Hegazy ,​ Ida Lindeman ,​ Guy Emerton ,​ Thomas Krausgruber ,​ ​ ​ ​ ​ ​ 8 7 3,4 1,10,11,* Jacqueline Shields ,​ Muzlifah Haniffa ,​ Fiona Powrie ,​ Sarah A. Teichmann ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ Affiliations: 1 Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, UK ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ 2 Centre of Biological Engineering, University of Minho, Braga, Portugal ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ 3 Kennedy Institute of Rheumatology, Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences, University of Oxford, Oxford, UK ​ ​ ​ ​ ​ ​ ​ ​ 4 Translational Gastroenterology Unit, Experimental Medicine Division Nuffield Department of Clinical Medicine, University of Oxford, John Radcliffe Hospital, Oxford, UK ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ 5 Current affiliation: Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health,
    [Show full text]