Modeling Gene Regulation from Paired Expression and Chromatin Accessibility Data

Total Page:16

File Type:pdf, Size:1020Kb

Modeling Gene Regulation from Paired Expression and Chromatin Accessibility Data Modeling gene regulation from paired expression and PNAS PLUS chromatin accessibility data Zhana Durena,b,c, Xi Chenb, Rui Jiangd,1, Yong Wanga,c,1, and Wing Hung Wongb,1 aAcademy of Mathematics and Systems Science, National Center for Mathematics and Interdisciplinary Sciences, Chinese Academy of Sciences, Beijing 100080, China; bDepartment of Statistics, Department of Biomedical Data Science, Bio-X Program, Stanford University, Stanford, CA 94305; cSchool of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China; and dMinistry of Education Key Laboratory of Bioinformatics, Bioinformatics Division and Center for Synthetic & Systems Biology, Tsinghua National Laboratory for Information Science and Technology, Department of Automation, Tsinghua University, Beijing 100084, China Contributed by Wing Hung Wong, May 8, 2017 (sent for review March 20, 2017; reviewed by Christina Kendziorski and Sheng Zhong) The rapid increase of genome-wide datasets on gene expression, gene expression data, accessibility data are available for a diverse set chromatin states, and transcription factor (TF) binding locations offers of cellular contexts (Fig. 1, blue boxes). In fact, we expect the an exciting opportunity to interpret the information encoded in amount of matched expression and accessibility data (i.e., measured genomes and epigenomes. This task can be challenging as it requires on the same sample) will increase very rapidly in the near future. joint modeling of context-specific activation of cis-regulatory ele- The purpose of the present work is to show that, by using ments (REs) and the effects on transcription of associated regulatory matched expression and accessibility data across diverse cellular factors. To meet this challenge, we propose a statistical approach contexts, it is possible to recover a significant portion of the in- based on paired expression and chromatin accessibility (PECA) data formation in the missing data on binding location and chromatin across diverse cellular contexts. In our approach, we model (i)the state and to achieve accurate inference of the gene regulatory localization to REs of chromatin regulators (CRs) based on their in- relations. In our approach, key events in the regulatory process, teraction with sequence-specific TFs, (ii) the activation of REs due to such as recruitment of chromatin remodeling factors to a regu- CRs that are localized to them, and (iii) the effect of TFs bound to latory element, activation of regulatory elements, etc., are activated REs on the transcription of target genes (TGs). The transcrip- regarded as latent unobserved variables in a statistical model that tional regulatory network inferred by PECA provides a detailed view describes the relations among these variables and the gene ex- trans cis of how -and -regulatory elements work together to affect pression variables, conditional on accessibility data on the regula- BIOPHYSICS AND gene expression in a context-specific manner. We illustrate the fea- tory elements. By fitting this model to expression and accessibility COMPUTATIONAL BIOLOGY sibility of this approach by analyzing paired expression and accessi- data across a large number of cellular contexts, we can infer many bility data from the mouse Encyclopedia of DNA Elements (ENCODE) details of the gene regulatory system helpful in the interpretation of and explore various applications of the resulting model. new data or the generation of new hypotheses. We end this Introduction with comments on related works. gene regulation | transcription factor | regulatory element | Several methods have recently been proposed to detect tran- chromatin regulator | chromatin activity scription factor (TF) binding sites by “footprinting” in which the presence of a bound TF is reflected by the shape of the DNase- ver since the emergence of high-throughput gene expression seq (or ATAC-seq) profile around its binding site (6–8). These Eexperiments (1), computational biologists have been inter- works focus on the effect of TF binding on the frequency of ested in the inference of gene regulatory relationships from gene cleavage near the site and do not attempt to model gene regulatory expression data across diverse cellular contexts corresponding to diverse cell types and experimental conditions (Fig. 1, red boxes). Significance However, progress has been hindered by the fact that gene expres- sion measurements provide little information on underlying regula- tory mechanisms such as transcription factor binding and chromatin Chromatin plays a critical role in the regulation of gene ex- modification. To fill this gap, chromatin immunoprecipitation- pression. Interactions among chromatin regulators, sequence- based methods (2, 3) have been developed for the genome-wide specific transcription factors, and cis-regulatory sequence ele- mapping of transcriptional regulator binding locations and the ments are the main driving forces shaping context-specific detection of epigenetic marks characteristic of specific chromatin chromatin structure and gene expression. However, because states. For example, by performing thousands of ChIP-seq exper- of the large number of such interactions, direct data on them iments, the Encyclopedia of DNA Elements (ENCODE) consor- are often missing in most cellular contexts. The purpose of the tium has generated such data for many chromatin marks and present work is to show that, by modeling matched expression transcriptional regulators on a small number of cell lines (Fig. 1, and accessibility data across diverse cellular contexts, it is green boxes). However, because a large number of transcriptional possible to recover a significant portion of the information in regulators and chromatin marks have to be analyzed one by one, it the missing data on binding locations and chromatin states and is unlikely that such comprehensive data will become available for to achieve accurate inference of gene regulatory relations. many other cell lines. For most cellular contexts, the desired data will remain missing in the foreseeable future (Fig. 1, gray boxes). Author contributions: R.J., Y.W., and W.H.W. designed research; Z.D. and X.C. performed – research; Z.D., R.J., Y.W., and W.H.W. analyzed data; and Z.D., R.J., Y.W., and W.H.W. On the other hand, it is known that many of the protein DNA wrote the paper. interactions important for gene regulation occur in regulatory Reviewers: C.K., University of Wisconsin; and S.Z., University of California, San Diego. elements (REs) such as enhancers and insulators, which com- The authors declare no conflict of interest. pose only a small portion of the noncoding sequences in a ge- nome. The REs active in gene regulation in a given cellular state Freely available online through the PNAS open access option. tend to have an open chromatin structure so that they are ac- Data deposition: The sequence data reported in this paper have been deposited in the Gene Expression Omnibus (GEO) database, https://www.ncbi.nlm.nih.gov/geo (accession cessible for binding by relevant transcriptional regulators. This no. GSE98479). suggests that many of the relevant regulatory relations may be 1To whom correspondence may be addressed. Email: [email protected], ruijiang@ revealed by analyzing the accessible REs. Fortunately, genome-wide tsinghua.edu.cn, or [email protected]. measurement of chromatin accessibility is now straightforward by This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. recentmethodssuchasDNase-seq(4)orATAC-seq(5).Similarto 1073/pnas.1704553114/-/DCSupplemental. www.pnas.org/cgi/doi/10.1073/pnas.1704553114 PNAS Early Edition | 1of10 Downloaded by guest on September 25, 2021 ChIP-seq Biological RNA-seq DNase-seq Ctcf Zmiz1 Bhlhe40 Cebpb Chd1 Chd2 Ep300 Fli1 Gabpa Gata1 Gata2 Hcfc1 Jund Mafk Maz Myb Myod1 Myog Pax5 Rad21 Rdbp Sin3a Srf Tbp Tcf12 Ubtf Usf1 Usf2 Zc3h11a Zkscan1 Znf384 Fosl1 Kat2a Max Polr2a Smc3 system Cell types ENCODE sample ID E2f4 Ets1 Jun Mxi1 Myc Rcor1 Rest Tal1 Tcf3 Muscular SkMuscle SkmuscleC57bl6MAdult8wks Circulatory G1E-ER4 G1eer4S129ME0Diffd24h G1 E G1eS12 9ME0 Cerebrum CerebrumC57bl6MAdult8wks Nervous Cerebellum CerebellumC57bl6MAdult8wks WholeBrain WbrainC57bl6ME18half Respiratory Lung LungC57bl6MAdult8wks NIH-3T3 Nih3t3NihsMImmortal LgIntestine LgintC57bl6MAdult8wks Liver liver129dlcrME14half Digestive Liver LiverC57bl6MAdult8wks Liver LiverC57bl6ME14half Excretory Kidney KidneyC57bl6MAdult8wks Endocrine FatPad FatC57bl6MAdult8wks GenitalFatPad GfatC57bl6MAdult8wks 416B 416bC57bl6MAdult8wks A20 A20BalbcannMAdult8wks B-cell(CD19+) Bcellcd19pC57bl6MAdult8wks B-cell(CD43-) Bcellcd43nC57bl6MAdult8wks MEL MelC57bl6MAdult8wks Spleen SpleenC57bl6MAdult8wks Lymphatic Thymus ThymusC57bl6MAdult8wks T-Naïve TnaiveC57bl6MAdult8wks Fig. 1. Genome-wide data for gene-regulatory inference. Each row represents a particular cellular context under which multiple types of genome-wide data may be available. In this paper we illustrate our method by analyzing data from 25 contexts studied in the mouse ENCODE project covering a variety of mouse cell types and developmental stages. Expression data (from RNA-seq, red boxes) and chromatin accessibility data (from DNase-seq or ATAC-seq, blue boxes) are available for each context, but most of the location data (from ChIP-seq data, green boxes) for transcriptional regulators are missing. We expect that the number of contexts (i.e., number of rows) with expression and accessibility data will increase rapidly in
Recommended publications
  • AP2A2 Antibody Cat
    AP2A2 Antibody Cat. No.: 63-513 AP2A2 Antibody Formalin-fixed and paraffin-embedded human Flow cytometric analysis of HepG2 cells using AP2A2 hepatocarcinoma with AP2A2 Antibody , which was Antibody (bottom histogram) compared to a negative peroxidase-conjugated to the secondary antibody, control cell (top histogram). FITC-conjugated goat-anti- followed by DAB staining. rabbit secondary antibodies were used for the analysis. Specifications HOST SPECIES: Rabbit SPECIES REACTIVITY: Human This AP2A2 antibody is generated from rabbits immunized with a KLH conjugated IMMUNOGEN: synthetic peptide between 610-637 amino acids from the Central region of human AP2A2. TESTED APPLICATIONS: Flow, IHC-P, WB For WB starting dilution is: 1:1000 APPLICATIONS: For IHC-P starting dilution is: 1:10~50 For FACS starting dilution is: 1:10~50 September 25, 2021 1 https://www.prosci-inc.com/ap2a2-antibody-63-513.html PREDICTED MOLECULAR 104 kDa WEIGHT: Properties This antibody is purified through a protein A column, followed by peptide affinity PURIFICATION: purification. CLONALITY: Polyclonal ISOTYPE: Rabbit Ig CONJUGATE: Unconjugated PHYSICAL STATE: Liquid BUFFER: Supplied in PBS with 0.09% (W/V) sodium azide. CONCENTRATION: batch dependent Store at 4˚C for three months and -20˚C, stable for up to one year. As with all antibodies STORAGE CONDITIONS: care should be taken to avoid repeated freeze thaw cycles. Antibodies should not be exposed to prolonged high temperatures. Additional Info OFFICIAL SYMBOL: AP2A2 AP-2 complex subunit alpha-2, 100 kDa coated vesicle
    [Show full text]
  • In Silico Prediction of High-Resolution Hi-C Interaction Matrices
    ARTICLE https://doi.org/10.1038/s41467-019-13423-8 OPEN In silico prediction of high-resolution Hi-C interaction matrices Shilu Zhang1, Deborah Chasman 1, Sara Knaack1 & Sushmita Roy1,2* The three-dimensional (3D) organization of the genome plays an important role in gene regulation bringing distal sequence elements in 3D proximity to genes hundreds of kilobases away. Hi-C is a powerful genome-wide technique to study 3D genome organization. Owing to 1234567890():,; experimental costs, high resolution Hi-C datasets are limited to a few cell lines. Computa- tional prediction of Hi-C counts can offer a scalable and inexpensive approach to examine 3D genome organization across multiple cellular contexts. Here we present HiC-Reg, an approach to predict contact counts from one-dimensional regulatory signals. HiC-Reg pre- dictions identify topologically associating domains and significant interactions that are enri- ched for CCCTC-binding factor (CTCF) bidirectional motifs and interactions identified from complementary sources. CTCF and chromatin marks, especially repressive and elongation marks, are most important for HiC-Reg’s predictive performance. Taken together, HiC-Reg provides a powerful framework to generate high-resolution profiles of contact counts that can be used to study individual locus level interactions and higher-order organizational units of the genome. 1 Wisconsin Institute for Discovery, 330 North Orchard Street, Madison, WI 53715, USA. 2 Department of Biostatistics and Medical Informatics, University of Wisconsin-Madison, Madison, WI 53715, USA. *email: [email protected] NATURE COMMUNICATIONS | (2019) 10:5449 | https://doi.org/10.1038/s41467-019-13423-8 | www.nature.com/naturecommunications 1 ARTICLE NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-019-13423-8 he three-dimensional (3D) organization of the genome has Results Temerged as an important component of the gene regulation HiC-Reg for predicting contact count using Random Forests.
    [Show full text]
  • Evolution, Dynamic Expression Changes and Regulatory Characteristics of Gene Families Involved in the Glycerophosphate Pathway O
    www.nature.com/scientificreports OPEN Evolution, dynamic expression changes and regulatory characteristics of gene families Received: 2 October 2018 Accepted: 14 August 2019 involved in the glycerophosphate Published: xx xx xxxx pathway of triglyceride synthesis in chicken (Gallus gallus) Liyu Yang1, Ziming Liu1, Kepeng Ou2, Taian Wang1, Zhuanjian Li1,3,4, Yadong Tian1,3,4, Yanbin Wang1,3,4, Xiangtao Kang1,3,4, Hong Li1,3,4 & Xiaojun Liu1,3,4 It is well documented that four gene families, including the glycerophosphate acyltransferases (GPATs), acylglycerophosphate acyltransferases (AGPATs), lipid phosphate phosphohydrolases (LPINs) and diacylglycerol acyltransferases (DGATs), are involved in the glycerophosphate pathway of de novo triglyceride (TG) biosynthesis in mammals. However, no systematic analysis has been conducted to characterize the gene families in poultry. In this study, the sequences of gene family members in the glycerophosphate pathway were obtained by screening the public databases. The phylogenetic tree, gene structures and conserved motifs of the corresponding proteins were evaluated. Dynamic expression changes of the genes at diferent developmental stages were analyzed by qRT-PCR. The regulatory characteristics of the genes were analyzed by in vivo experiments. The results showed that the GPAT, AGPAT and LPIN gene families have 2, 7 and 2 members, respectively, and they were classifed into 2, 4 and 2 cluster respectively based on phylogenetic analysis. All of the genes except AGPAT1 were extensively expressed in various tissues. Estrogen induction upregulated the expression of GPAM and AGPAT2, downregulated the expression of AGPAT3, AGPAT9, LPIN1 and LPIN2, and had no efect on the expression of the other genes. These fndings provide a valuable resource for further investigation of lipid metabolism in liver of chicken.
    [Show full text]
  • Global-Scale Analysis of the Dynamic Transcriptional Adaptations Within Skeletal Muscle During Hypertrophic Growth
    University of Kentucky UKnowledge Theses and Dissertations--Physiology Physiology 2015 GLOBAL-SCALE ANALYSIS OF THE DYNAMIC TRANSCRIPTIONAL ADAPTATIONS WITHIN SKELETAL MUSCLE DURING HYPERTROPHIC GROWTH Tyler Kirby University of Kentucky, [email protected] Right click to open a feedback form in a new tab to let us know how this document benefits ou.y Recommended Citation Kirby, Tyler, "GLOBAL-SCALE ANALYSIS OF THE DYNAMIC TRANSCRIPTIONAL ADAPTATIONS WITHIN SKELETAL MUSCLE DURING HYPERTROPHIC GROWTH" (2015). Theses and Dissertations--Physiology. 22. https://uknowledge.uky.edu/physiology_etds/22 This Doctoral Dissertation is brought to you for free and open access by the Physiology at UKnowledge. It has been accepted for inclusion in Theses and Dissertations--Physiology by an authorized administrator of UKnowledge. For more information, please contact [email protected]. STUDENT AGREEMENT: I represent that my thesis or dissertation and abstract are my original work. Proper attribution has been given to all outside sources. I understand that I am solely responsible for obtaining any needed copyright permissions. I have obtained needed written permission statement(s) from the owner(s) of each third-party copyrighted matter to be included in my work, allowing electronic distribution (if such use is not permitted by the fair use doctrine) which will be submitted to UKnowledge as Additional File. I hereby grant to The University of Kentucky and its agents the irrevocable, non-exclusive, and royalty-free license to archive and make accessible my work in whole or in part in all forms of media, now or hereafter known. I agree that the document mentioned above may be made available immediately for worldwide access unless an embargo applies.
    [Show full text]
  • Gene Regulation and Speciation in House Mice
    Downloaded from genome.cshlp.org on September 26, 2021 - Published by Cold Spring Harbor Laboratory Press Research Gene regulation and speciation in house mice Katya L. Mack,1 Polly Campbell,2 and Michael W. Nachman1 1Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, Berkeley, California 94720-3160, USA; 2Department of Integrative Biology, Oklahoma State University, Stillwater, Oklahoma 74078, USA One approach to understanding the process of speciation is to characterize the genetic architecture of post-zygotic isolation. As gene regulation requires interactions between loci, negative epistatic interactions between divergent regulatory elements might underlie hybrid incompatibilities and contribute to reproductive isolation. Here, we take advantage of a cross between house mouse subspecies, where hybrid dysfunction is largely unidirectional, to test several key predictions about regulatory divergence and reproductive isolation. Regulatory divergence between Mus musculus musculus and M. m. domesticus was charac- terized by studying allele-specific expression in fertile hybrid males using mRNA-sequencing of whole testes. We found ex- tensive regulatory divergence between M. m. musculus and M. m. domesticus, largely attributable to cis-regulatory changes. When both cis and trans changes occurred, they were observed in opposition much more often than expected under a neutral model, providing strong evidence of widespread compensatory evolution. We also found evidence for lineage-specific positive se- lection on a subset of genes related to transcriptional regulation. Comparisons of fertile and sterile hybrid males identified a set of genes that were uniquely misexpressed in sterile individuals. Lastly, we discovered a nonrandom association between these genes and genes showing evidence of compensatory evolution, consistent with the idea that regulatory interactions might contribute to Dobzhansky-Muller incompatibilities and be important in speciation.
    [Show full text]
  • Seq2pathway Vignette
    seq2pathway Vignette Bin Wang, Xinan Holly Yang, Arjun Kinstlick May 19, 2021 Contents 1 Abstract 1 2 Package Installation 2 3 runseq2pathway 2 4 Two main functions 3 4.1 seq2gene . .3 4.1.1 seq2gene flowchart . .3 4.1.2 runseq2gene inputs/parameters . .5 4.1.3 runseq2gene outputs . .8 4.2 gene2pathway . 10 4.2.1 gene2pathway flowchart . 11 4.2.2 gene2pathway test inputs/parameters . 11 4.2.3 gene2pathway test outputs . 12 5 Examples 13 5.1 ChIP-seq data analysis . 13 5.1.1 Map ChIP-seq enriched peaks to genes using runseq2gene .................... 13 5.1.2 Discover enriched GO terms using gene2pathway_test with gene scores . 15 5.1.3 Discover enriched GO terms using Fisher's Exact test without gene scores . 17 5.1.4 Add description for genes . 20 5.2 RNA-seq data analysis . 20 6 R environment session 23 1 Abstract Seq2pathway is a novel computational tool to analyze functional gene-sets (including signaling pathways) using variable next-generation sequencing data[1]. Integral to this tool are the \seq2gene" and \gene2pathway" components in series that infer a quantitative pathway-level profile for each sample. The seq2gene function assigns phenotype-associated significance of genomic regions to gene-level scores, where the significance could be p-values of SNPs or point mutations, protein-binding affinity, or transcriptional expression level. The seq2gene function has the feasibility to assign non-exon regions to a range of neighboring genes besides the nearest one, thus facilitating the study of functional non-coding elements[2]. Then the gene2pathway summarizes gene-level measurements to pathway-level scores, comparing the quantity of significance for gene members within a pathway with those outside a pathway.
    [Show full text]
  • A Computational Approach for Defining a Signature of Β-Cell Golgi Stress in Diabetes Mellitus
    Page 1 of 781 Diabetes A Computational Approach for Defining a Signature of β-Cell Golgi Stress in Diabetes Mellitus Robert N. Bone1,6,7, Olufunmilola Oyebamiji2, Sayali Talware2, Sharmila Selvaraj2, Preethi Krishnan3,6, Farooq Syed1,6,7, Huanmei Wu2, Carmella Evans-Molina 1,3,4,5,6,7,8* Departments of 1Pediatrics, 3Medicine, 4Anatomy, Cell Biology & Physiology, 5Biochemistry & Molecular Biology, the 6Center for Diabetes & Metabolic Diseases, and the 7Herman B. Wells Center for Pediatric Research, Indiana University School of Medicine, Indianapolis, IN 46202; 2Department of BioHealth Informatics, Indiana University-Purdue University Indianapolis, Indianapolis, IN, 46202; 8Roudebush VA Medical Center, Indianapolis, IN 46202. *Corresponding Author(s): Carmella Evans-Molina, MD, PhD ([email protected]) Indiana University School of Medicine, 635 Barnhill Drive, MS 2031A, Indianapolis, IN 46202, Telephone: (317) 274-4145, Fax (317) 274-4107 Running Title: Golgi Stress Response in Diabetes Word Count: 4358 Number of Figures: 6 Keywords: Golgi apparatus stress, Islets, β cell, Type 1 diabetes, Type 2 diabetes 1 Diabetes Publish Ahead of Print, published online August 20, 2020 Diabetes Page 2 of 781 ABSTRACT The Golgi apparatus (GA) is an important site of insulin processing and granule maturation, but whether GA organelle dysfunction and GA stress are present in the diabetic β-cell has not been tested. We utilized an informatics-based approach to develop a transcriptional signature of β-cell GA stress using existing RNA sequencing and microarray datasets generated using human islets from donors with diabetes and islets where type 1(T1D) and type 2 diabetes (T2D) had been modeled ex vivo. To narrow our results to GA-specific genes, we applied a filter set of 1,030 genes accepted as GA associated.
    [Show full text]
  • Supplementary Materials
    1 Supplementary Materials: Supplemental Figure 1. Gene expression profiles of kidneys in the Fcgr2b-/- and Fcgr2b-/-. Stinggt/gt mice. (A) A heat map of microarray data show the genes that significantly changed up to 2 fold compared between Fcgr2b-/- and Fcgr2b-/-. Stinggt/gt mice (N=4 mice per group; p<0.05). Data show in log2 (sample/wild-type). 2 Supplemental Figure 2. Sting signaling is essential for immuno-phenotypes of the Fcgr2b-/-lupus mice. (A-C) Flow cytometry analysis of splenocytes isolated from wild-type, Fcgr2b-/- and Fcgr2b-/-. Stinggt/gt mice at the age of 6-7 months (N= 13-14 per group). Data shown in the percentage of (A) CD4+ ICOS+ cells, (B) B220+ I-Ab+ cells and (C) CD138+ cells. Data show as mean ± SEM (*p < 0.05, **p<0.01 and ***p<0.001). 3 Supplemental Figure 3. Phenotypes of Sting activated dendritic cells. (A) Representative of western blot analysis from immunoprecipitation with Sting of Fcgr2b-/- mice (N= 4). The band was shown in STING protein of activated BMDC with DMXAA at 0, 3 and 6 hr. and phosphorylation of STING at Ser357. (B) Mass spectra of phosphorylation of STING at Ser357 of activated BMDC from Fcgr2b-/- mice after stimulated with DMXAA for 3 hour and followed by immunoprecipitation with STING. (C) Sting-activated BMDC were co-cultured with LYN inhibitor PP2 and analyzed by flow cytometry, which showed the mean fluorescence intensity (MFI) of IAb expressing DC (N = 3 mice per group). 4 Supplemental Table 1. Lists of up and down of regulated proteins Accession No.
    [Show full text]
  • 4-6 Weeks Old Female C57BL/6 Mice Obtained from Jackson Labs Were Used for Cell Isolation
    Methods Mice: 4-6 weeks old female C57BL/6 mice obtained from Jackson labs were used for cell isolation. Female Foxp3-IRES-GFP reporter mice (1), backcrossed to B6/C57 background for 10 generations, were used for the isolation of naïve CD4 and naïve CD8 cells for the RNAseq experiments. The mice were housed in pathogen-free animal facility in the La Jolla Institute for Allergy and Immunology and were used according to protocols approved by the Institutional Animal Care and use Committee. Preparation of cells: Subsets of thymocytes were isolated by cell sorting as previously described (2), after cell surface staining using CD4 (GK1.5), CD8 (53-6.7), CD3ε (145- 2C11), CD24 (M1/69) (all from Biolegend). DP cells: CD4+CD8 int/hi; CD4 SP cells: CD4CD3 hi, CD24 int/lo; CD8 SP cells: CD8 int/hi CD4 CD3 hi, CD24 int/lo (Fig S2). Peripheral subsets were isolated after pooling spleen and lymph nodes. T cells were enriched by negative isolation using Dynabeads (Dynabeads untouched mouse T cells, 11413D, Invitrogen). After surface staining for CD4 (GK1.5), CD8 (53-6.7), CD62L (MEL-14), CD25 (PC61) and CD44 (IM7), naïve CD4+CD62L hiCD25-CD44lo and naïve CD8+CD62L hiCD25-CD44lo were obtained by sorting (BD FACS Aria). Additionally, for the RNAseq experiments, CD4 and CD8 naïve cells were isolated by sorting T cells from the Foxp3- IRES-GFP mice: CD4+CD62LhiCD25–CD44lo GFP(FOXP3)– and CD8+CD62LhiCD25– CD44lo GFP(FOXP3)– (antibodies were from Biolegend). In some cases, naïve CD4 cells were cultured in vitro under Th1 or Th2 polarizing conditions (3, 4).
    [Show full text]
  • Mechanisms of Synaptic Plasticity Mediated by Clathrin Adaptor-Protein Complexes 1 and 2 in Mice
    Mechanisms of synaptic plasticity mediated by Clathrin Adaptor-protein complexes 1 and 2 in mice Dissertation for the award of the degree “Doctor rerum naturalium” at the Georg-August-University Göttingen within the doctoral program “Molecular Biology of Cells” of the Georg-August University School of Science (GAUSS) Submitted by Ratnakar Mishra Born in Birpur, Bihar, India Göttingen, Germany 2019 1 Members of the Thesis Committee Prof. Dr. Peter Schu Institute for Cellular Biochemistry, (Supervisor and first referee) University Medical Center Göttingen, Germany Dr. Hans Dieter Schmitt Neurobiology, Max Planck Institute (Second referee) for Biophysical Chemistry, Göttingen, Germany Prof. Dr. med. Thomas A. Bayer Division of Molecular Psychiatry, University Medical Center, Göttingen, Germany Additional Members of the Examination Board Prof. Dr. Silvio O. Rizzoli Department of Neuro-and Sensory Physiology, University Medical Center Göttingen, Germany Dr. Roland Dosch Institute of Developmental Biochemistry, University Medical Center Göttingen, Germany Prof. Dr. med. Martin Oppermann Institute of Cellular and Molecular Immunology, University Medical Center, Göttingen, Germany Date of oral examination: 14th may 2019 2 Table of Contents List of abbreviations ................................................................................. 5 Abstract ................................................................................................... 7 Chapter 1: Introduction ............................................................................
    [Show full text]
  • Supplemental Information
    Supplemental information Dissection of the genomic structure of the miR-183/96/182 gene. Previously, we showed that the miR-183/96/182 cluster is an intergenic miRNA cluster, located in a ~60-kb interval between the genes encoding nuclear respiratory factor-1 (Nrf1) and ubiquitin-conjugating enzyme E2H (Ube2h) on mouse chr6qA3.3 (1). To start to uncover the genomic structure of the miR- 183/96/182 gene, we first studied genomic features around miR-183/96/182 in the UCSC genome browser (http://genome.UCSC.edu/), and identified two CpG islands 3.4-6.5 kb 5’ of pre-miR-183, the most 5’ miRNA of the cluster (Fig. 1A; Fig. S1 and Seq. S1). A cDNA clone, AK044220, located at 3.2-4.6 kb 5’ to pre-miR-183, encompasses the second CpG island (Fig. 1A; Fig. S1). We hypothesized that this cDNA clone was derived from 5’ exon(s) of the primary transcript of the miR-183/96/182 gene, as CpG islands are often associated with promoters (2). Supporting this hypothesis, multiple expressed sequences detected by gene-trap clones, including clone D016D06 (3, 4), were co-localized with the cDNA clone AK044220 (Fig. 1A; Fig. S1). Clone D016D06, deposited by the German GeneTrap Consortium (GGTC) (http://tikus.gsf.de) (3, 4), was derived from insertion of a retroviral construct, rFlpROSAβgeo in 129S2 ES cells (Fig. 1A and C). The rFlpROSAβgeo construct carries a promoterless reporter gene, the β−geo cassette - an in-frame fusion of the β-galactosidase and neomycin resistance (Neor) gene (5), with a splicing acceptor (SA) immediately upstream, and a polyA signal downstream of the β−geo cassette (Fig.
    [Show full text]
  • Investigation of the Underlying Hub Genes and Molexular Pathogensis in Gastric Cancer by Integrated Bioinformatic Analyses
    bioRxiv preprint doi: https://doi.org/10.1101/2020.12.20.423656; this version posted December 22, 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. Investigation of the underlying hub genes and molexular pathogensis in gastric cancer by integrated bioinformatic analyses Basavaraj Vastrad1, Chanabasayya Vastrad*2 1. Department of Biochemistry, Basaveshwar College of Pharmacy, Gadag, Karnataka 582103, India. 2. Biostatistics and Bioinformatics, Chanabasava Nilaya, Bharthinagar, Dharwad 580001, Karanataka, India. * Chanabasayya Vastrad [email protected] Ph: +919480073398 Chanabasava Nilaya, Bharthinagar, Dharwad 580001 , Karanataka, India bioRxiv preprint doi: https://doi.org/10.1101/2020.12.20.423656; this version posted December 22, 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 The high mortality rate of gastric cancer (GC) is in part due to the absence of initial disclosure of its biomarkers. The recognition of important genes associated in GC is therefore recommended to advance clinical prognosis, diagnosis and and treatment outcomes. The current investigation used the microarray dataset GSE113255 RNA seq data from the Gene Expression Omnibus database to diagnose differentially expressed genes (DEGs). Pathway and gene ontology enrichment analyses were performed, and a proteinprotein interaction network, modules, target genes - miRNA regulatory network and target genes - TF regulatory network were constructed and analyzed. Finally, validation of hub genes was performed. The 1008 DEGs identified consisted of 505 up regulated genes and 503 down regulated genes.
    [Show full text]