DOCSLIB.ORG

Screening for Implicated Genes in Colorectal Cancer Using Whole‑Genome Gene Expression Profiling

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text

Load more

8260 MOLECULAR MEDICINE REPORTS 17: 8260-8268, 2018 Screening for implicated genes in colorectal cancer using whole‑genome gene expression profiling LONG-CI SUN1,2 and HAI-XIN QIAN1 1Department of General Surgery, The First Affiliated Hospital of Soochow University, Suzhou, Jiangsu 215006; 2Department of Gastrointestinal Surgery, Renji Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai 200127, P.R. China Received October 24, 2017; Accepted March 2, 2018 DOI: 10.3892/mmr.2018.8862 Abstract. To identify biologically relevant genes associated WD repeat domains may represent novel target molecules for with the pathogenesis of colorectal cancer (CRC), genome future treatments of CRC. wide expression profiles of 17 pairs of CRC tumor and adjacent tissues, previously published in a DNA microarray study, were Introduction analyzed. Cytoscape, String tools and DAVID tools were used to investigate the biological pathways encoded by the genes Colorectal cancer (CRC) is one of the most commonly diag- identified as being either upregulated or downregulated in nosed cancers worldwide (1). It is the second most prevalent CRC, to determine protein-protein interactions and to identify cancer type among women, following breast cancer, and the potential hub genes associated with CRC. As a result, a total of third most prevalent cancer type among men (2,3). In recent 3,264 genes were identified as being differentially expressed decades, diagnoses and treatment of CRC have significantly in CRC and adjacent tissues, including 1,594 downregulated improved. However, the disease remains a major public health and 1,670 upregulated genes. Furthermore, 306 genes were problem worldwide. In 2016, the American Cancer Society revealed to be clustered in a complex interaction network, estimated 70,820 new male cases and 63,670 new female and the top 20 hub genes in this network were determined by cases of CRC would be diagnosed, and that there would be application of the Matthews Correlation Coefficient algorithm. 49,190 CRC associated deaths in the USA (4). The high rate In addition, the patterns of the expression levels of the 20 hub of mortality of the disease is in part due to limitations in the genes were investigated using reverse transcription‑quantitative currently available therapies for CRC treatment, which is the polymerase chain reaction. Gene Ontology analysis revealed result of an incomplete understanding of the biological mecha- that four of the 20 hub genes encoded small subunit processome nisms underlying the disease. components (UTP3 small subunit processome component; It has been well established that genes have an important UTP14 small subunit processome component; UTP 18 role in tumor pathogenesis, and that the transcriptional small subunit processome component; and UTP20 small activation or inhibition of certain genes serves an important role subunit processome component) and a further four encoded in the development and progression of the majority of human WD repeat domains (WD repeat-containing protein 3, WD tumors (5). The analysis of genome‑wide gene expression repeat domain 12, WD repeat-containing protein 43 and WD levels via DNA microarray experiments is a recently developed repeat‑containing protein 75). In conclusion, the present DNA systematic approach that is used to gain comprehensive microarray study identified genes involved in the pathogenesis information regarding tumor associated gene transcription of CRC. Furthermore, it was revealed that hub genes identified profiles (6). This information can lead to the identification of from among the total identified upregulated and downregulated prognostic biomarkers (7,8), enable the discrimination between genes in CRC encoding subunit processome components and various histological subtypes of tumors (9,10). Interactions between genes can affect the pathogenesis of CRC and tumor progression. The aim of the present study was to identify biologically relevant genes that may be associated with the pathogenesis of CRC, via analysis of genome-wide expression profiles in CRC. Following the identification of Correspondence to: Professor Hai-Xin Qian, Department differentially expressed genes (DEGs), the present study further of General Surgery, The First Affiliated Hospital of Soochow University, 899 Pinghai Road, Suzhou, Jiangsu 215006, P.R. China aimed to determine the hub genes involved in the regulatory E‑mail: [email protected] pathways of CRC progression, by Gene Ontology (GO) analysis, and to subsequently derive a gene interaction network. Key words: colorectal cancer, microarray, hub gene, expression Identification of hub genes in CRC and the development of a gene profiling interaction network may reveal the biological characteristics of CRC tumor development and enable the development of novel molecular‑targeted therapies for the treatment of CRC. SUN and QIAN: SCREENING FOR IMPLICATED GENES IN COLORECTAL CANCER 8261 Materials and methods Hospital Affiliated to Shanghai Jiaotong University School of Medicine (Shanghai, China) between November 2017 and Data collection. The gene expression profiling dataset from December 2017. The inclusion criteria for patients were as previously published relevant study was retrieved from follows: i) colorectal adenocarcinoma confirmed by patho- the NCBI GEO database (accession number GSE32323; logical examination; ii) the patients must have not received a www.ncbi.nlm.nih.gov/geo) (11). Data regarding a total of radical resection prior and must not have been treated with 17 paired samples of CRC tumors and adjacent tissues were preoperative chemoradiotherapy. Characteristics of patients available from the GSE32323 dataset. included in the present study are presented in Table I. Written informed consent was obtained from all patients included in Data processing. Raw RNA expression profile datasets the present study. The Ethics Committee of Renji Hospital were preprocessed using R 3.4.1 statistical software Affiliated to Shanghai Jiaotong University School of Medicine (https://www.r‑project.org/) together with Affy package (12). approved the present study and protocols. CRC tissue samples In accordance with the R‑package, the gcrma algorithm (13) (diameter of 0.5 cm) were obtained via tumor dissection. was used to adjust for background intensities in the Affymtrix Paracancerous tissue (distal tissue 3 cm away from the tumor) array data by including optical noise and non‑specific was also obtained from patients via dissection. Paracancerous binding. The background adjusted probe intensities were then tissues and cancerous tissues were collected as normal tissues converted into expression measures using the normaliza- and tumor tissues, respectively. RNA was extracted from these tion and summarization methods included in the multiarray tissues by TRIzol (Invitrogen; Thermo Fisher Scientific, Inc., average algorithm (12). The k‑Nearest Neighbor algorithm (14) Waltham, MA, USA) for future experiments. was subsequently used to generate the missing values. Reverse transcription‑quantitative polymerase chain reaction Microarray and statistical analysis. The R‑package and (RT‑qPCR). In order to further investigate the results of limma (www.bioconductor.org.uk) (15) were used to inves- the microarray analysis, the 20 hub genes identified to be tigate gene expression microarray data. Limma uses linear associated with CRC were subjected to RT‑qPCR analysis. A models, as well as empirical Bayesian methods, to investigate total of 1 µg of RNA from tumor or normal tissue samples was the differential expression levels in the GSE32323 dataset. reverse transcribed to cDNA using SuperScript™ II reverse Fold‑change (FC) values were calculated using the ratio of transcriptase (Invitrogen; Thermo Fisher Scientific, Inc.; the geometric means of the gene expression levels between heat mixture to 65˚C for 5 min, incubate at 42˚C for 50 min, the two treatment groups (CRC tissue vs. adjacent tissue). A inactivate the reaction by heating at 70˚C for 15 min. RT‑qPCR Wilcoxon Rank‑Sum Test was performed to determine the was performed using cDNA, iQ™ SYBR Green Supermix statistical significance of differences between the groups (Invitrogen; Thermo Fisher Scientific, Inc.) and primers (P<0.05 was considered statistically significant). Genes were (2.5 µM). The reaction was performed using the Bio‑Rad selected for if the log2‑FC values were >1 (equivalent to a RT‑qPCR detection system (Bio‑Rad Laboratories, Inc., FC value of >2.0 or <0.5), and P‑values determined by the Hercules, CA, USA) under the following PCR conditions: Initial Wilcoxon Rank‑Sum Test were P<0.0002. Hierarchical clus- denaturation at 95˚C for 5 min; 35 cycles at 95˚C for 30 sec; tering of the selected genes was performed with R software 60˚C for 30 sec; 72˚C for 30 sec; and final extension at 72˚C using the Euclidean distance and complete linkage method. for 5 min. Expression levels of these genes were normalized Regarding clustering, the expression data were standardized to the house‑keeping gene, β-actin via the 2-ΔΔCq method (19). via conversion to z scores (mean=0, variance=1) for each probe. Primer sequences used were as follows: WD repeat-containing Following this, protein‑protein interactions among the selected protein 3 (WDR3) forward, 5 ' ‑ A A T CCA GCG GGT GAC genes were analyzed using the reference data contained in TAA‑3' and reverse, 5 ' ‑ A C A GGC TGA GGA GTA GGC‑3'; WD the String database of known protein‑protein interactions, repeat‑containing protein 43 (WDR43)
Recommended publications
  • The DEAD-Box RNA Helicase-Like Utp25 Is an SSU Processome Component

    The DEAD-Box RNA Helicase-Like Utp25 Is an SSU Processome Component

    Downloaded from rnajournal.cshlp.org on September 25, 2021 - Published by Cold Spring Harbor Laboratory Press The DEAD-box RNA helicase-like Utp25 is an SSU processome component J. MICHAEL CHARETTE1,2 and SUSAN J. BASERGA1,2,3 1Department of Molecular Biophysics & Biochemistry, Yale University School of Medicine, New Haven, Connecticut 06520, USA 2Department of Therapeutic Radiology, Yale University School of Medicine, New Haven, Connecticut 06520, USA 3Department of Genetics, Yale University School of Medicine, New Haven, Connecticut 06520, USA ABSTRACT The SSU processome is a large ribonucleoprotein complex consisting of the U3 snoRNA and at least 43 proteins. A database search, initiated in an effort to discover additional SSU processome components, identified the uncharacterized, conserved and essential yeast nucleolar protein YIL091C/UTP25 as one such candidate. The C-terminal DUF1253 motif, a domain of unknown function, displays limited sequence similarity to DEAD-box RNA helicases. In the absence of the conserved DEAD-box sequence, motif Ia is the only clearly identifiable helicase element. Since the yeast homolog is nucleolar and interacts with components of the SSU processome, we examined its role in pre-rRNA processing. Genetic depletion of Utp25 resulted in slowed growth. Northern analysis of pre-rRNA revealed an 18S rRNA maturation defect at sites A0,A1, and A2. Coimmunoprecipitation confirmed association with U3 snoRNA and with Mpp10, and with components of the t-Utp/UtpA, UtpB, and U3 snoRNP subcomplexes. Mutation of the conserved motif Ia residues resulted in no discernable temperature- sensitive or cold-sensitive growth defects, implying that this motif is dispensable for Utp25 function.
  • Analysis of Gene Expression Data for Gene Ontology

    Analysis of Gene Expression Data for Gene Ontology

    ANALYSIS OF GENE EXPRESSION DATA FOR GENE ONTOLOGY BASED PROTEIN FUNCTION PREDICTION A Thesis Presented to The Graduate Faculty of The University of Akron In Partial Fulfillment of the Requirements for the Degree Master of Science Robert Daniel Macholan May 2011 ANALYSIS OF GENE EXPRESSION DATA FOR GENE ONTOLOGY BASED PROTEIN FUNCTION PREDICTION Robert Daniel Macholan Thesis Approved: Accepted: _______________________________ _______________________________ Advisor Department Chair Dr. Zhong-Hui Duan Dr. Chien-Chung Chan _______________________________ _______________________________ Committee Member Dean of the College Dr. Chien-Chung Chan Dr. Chand K. Midha _______________________________ _______________________________ Committee Member Dean of the Graduate School Dr. Yingcai Xiao Dr. George R. Newkome _______________________________ Date ii ABSTRACT A tremendous increase in genomic data has encouraged biologists to turn to bioinformatics in order to assist in its interpretation and processing. One of the present challenges that need to be overcome in order to understand this data more completely is the development of a reliable method to accurately predict the function of a protein from its genomic information. This study focuses on developing an effective algorithm for protein function prediction. The algorithm is based on proteins that have similar expression patterns. The similarity of the expression data is determined using a novel measure, the slope matrix. The slope matrix introduces a normalized method for the comparison of expression levels throughout a proteome. The algorithm is tested using real microarray gene expression data. Their functions are characterized using gene ontology annotations. The results of the case study indicate the protein function prediction algorithm developed is comparable to the prediction algorithms that are based on the annotations of homologous proteins.
  • A Computational Approach for Defining a Signature of Β-Cell Golgi Stress in Diabetes Mellitus

    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.
  • 4-6 Weeks Old Female C57BL/6 Mice Obtained from Jackson Labs Were Used for Cell Isolation

    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).
  • Proteomics Provides Insights Into the Inhibition of Chinese Hamster V79

    Proteomics Provides Insights Into the Inhibition of Chinese Hamster V79

    www.nature.com/scientificreports OPEN Proteomics provides insights into the inhibition of Chinese hamster V79 cell proliferation in the deep underground environment Jifeng Liu1,2, Tengfei Ma1,2, Mingzhong Gao3, Yilin Liu4, Jun Liu1, Shichao Wang2, Yike Xie2, Ling Wang2, Juan Cheng2, Shixi Liu1*, Jian Zou1,2*, Jiang Wu2, Weimin Li2 & Heping Xie2,3,5 As resources in the shallow depths of the earth exhausted, people will spend extended periods of time in the deep underground space. However, little is known about the deep underground environment afecting the health of organisms. Hence, we established both deep underground laboratory (DUGL) and above ground laboratory (AGL) to investigate the efect of environmental factors on organisms. Six environmental parameters were monitored in the DUGL and AGL. Growth curves were recorded and tandem mass tag (TMT) proteomics analysis were performed to explore the proliferative ability and diferentially abundant proteins (DAPs) in V79 cells (a cell line widely used in biological study in DUGLs) cultured in the DUGL and AGL. Parallel Reaction Monitoring was conducted to verify the TMT results. γ ray dose rate showed the most detectable diference between the two laboratories, whereby γ ray dose rate was signifcantly lower in the DUGL compared to the AGL. V79 cell proliferation was slower in the DUGL. Quantitative proteomics detected 980 DAPs (absolute fold change ≥ 1.2, p < 0.05) between V79 cells cultured in the DUGL and AGL. Of these, 576 proteins were up-regulated and 404 proteins were down-regulated in V79 cells cultured in the DUGL. KEGG pathway analysis revealed that seven pathways (e.g.
  • Supp Material.Pdf

    Supp Material.Pdf

    Simon et al. Supplementary information: Table of contents p.1 Supplementary material and methods p.2-4 • PoIy(I)-poly(C) Treatment • Flow Cytometry and Immunohistochemistry • Western Blotting • Quantitative RT-PCR • Fluorescence In Situ Hybridization • RNA-Seq • Exome capture • Sequencing Supplementary Figures and Tables Suppl. items Description pages Figure 1 Inactivation of Ezh2 affects normal thymocyte development 5 Figure 2 Ezh2 mouse leukemias express cell surface T cell receptor 6 Figure 3 Expression of EZH2 and Hox genes in T-ALL 7 Figure 4 Additional mutation et deletion of chromatin modifiers in T-ALL 8 Figure 5 PRC2 expression and activity in human lymphoproliferative disease 9 Figure 6 PRC2 regulatory network (String analysis) 10 Table 1 Primers and probes for detection of PRC2 genes 11 Table 2 Patient and T-ALL characteristics 12 Table 3 Statistics of RNA and DNA sequencing 13 Table 4 Mutations found in human T-ALLs (see Fig. 3D and Suppl. Fig. 4) 14 Table 5 SNP populations in analyzed human T-ALL samples 15 Table 6 List of altered genes in T-ALL for DAVID analysis 20 Table 7 List of David functional clusters 31 Table 8 List of acquired SNP tested in normal non leukemic DNA 32 1 Simon et al. Supplementary Material and Methods PoIy(I)-poly(C) Treatment. pIpC (GE Healthcare Lifesciences) was dissolved in endotoxin-free D-PBS (Gibco) at a concentration of 2 mg/ml. Mice received four consecutive injections of 150 μg pIpC every other day. The day of the last pIpC injection was designated as day 0 of experiment.
  • Gene Expression Profiling Analysis Contributes to Understanding the Association Between Non-Syndromic Cleft Lip and Palate, and Cancer

    Gene Expression Profiling Analysis Contributes to Understanding the Association Between Non-Syndromic Cleft Lip and Palate, and Cancer

    2110 MOLECULAR MEDICINE REPORTS 13: 2110-2116, 2016 Gene expression profiling analysis contributes to understanding the association between non-syndromic cleft lip and palate, and cancer HONGYI WANG, TAO QIU, JIE SHI, JIULONG LIANG, YANG WANG, LIANGLIANG QUAN, YU ZHANG, QIAN ZHANG and KAI TAO Department of Plastic Surgery, General Hospital of Shenyang Military Area Command, PLA, Shenyang, Liaoning 110016, P.R. China Received March 10, 2015; Accepted December 18, 2015 DOI: 10.3892/mmr.2016.4802 Abstract. The present study aimed to investigate the for NSCL/P were implicated predominantly in the TGF-β molecular mechanisms underlying non-syndromic cleft lip, signaling pathway, the cell cycle and in viral carcinogenesis. with or without cleft palate (NSCL/P), and the association The TP53, CDK1, SMAD3, PIK3R1 and CASP3 genes were between this disease and cancer. The GSE42589 data set found to be associated, not only with NSCL/P, but also with was downloaded from the Gene Expression Omnibus data- cancer. These results may contribute to a better understanding base, and contained seven dental pulp stem cell samples of the molecular mechanisms of NSCL/P. from children with NSCL/P in the exfoliation period, and six controls. Differentially expressed genes (DEGs) were Introduction screened using the RankProd method, and their potential functions were revealed by pathway enrichment analysis and Non-syndromic cleft lip, with or without cleft palate (NSCL/P) construction of a pathway interaction network. Subsequently, is one of the most common types of congenital defect and cancer genes were obtained from six cancer databases, and affects 3.4-22.9/10,000 individuals worldwide (1).
  • Utpa and Utpb Chaperone Nascent Pre-Ribosomal RNA and U3 Snorna to Initiate Eukaryotic Ribosome Assembly

    Utpa and Utpb Chaperone Nascent Pre-Ribosomal RNA and U3 Snorna to Initiate Eukaryotic Ribosome Assembly

    ARTICLE Received 6 Apr 2016 | Accepted 27 May 2016 | Published 29 Jun 2016 DOI: 10.1038/ncomms12090 OPEN UtpA and UtpB chaperone nascent pre-ribosomal RNA and U3 snoRNA to initiate eukaryotic ribosome assembly Mirjam Hunziker1,*, Jonas Barandun1,*, Elisabeth Petfalski2, Dongyan Tan3, Cle´mentine Delan-Forino2, Kelly R. Molloy4, Kelly H. Kim5, Hywel Dunn-Davies2, Yi Shi4, Malik Chaker-Margot1,6, Brian T. Chait4, Thomas Walz5, David Tollervey2 & Sebastian Klinge1 Early eukaryotic ribosome biogenesis involves large multi-protein complexes, which co-transcriptionally associate with pre-ribosomal RNA to form the small subunit processome. The precise mechanisms by which two of the largest multi-protein complexes—UtpA and UtpB—interact with nascent pre-ribosomal RNA are poorly understood. Here, we combined biochemical and structural biology approaches with ensembles of RNA–protein cross-linking data to elucidate the essential functions of both complexes. We show that UtpA contains a large composite RNA-binding site and captures the 50 end of pre-ribosomal RNA. UtpB forms an extended structure that binds early pre-ribosomal intermediates in close proximity to architectural sites such as an RNA duplex formed by the 50 ETS and U3 snoRNA as well as the 30 boundary of the 18S rRNA. Both complexes therefore act as vital RNA chaperones to initiate eukaryotic ribosome assembly. 1 Laboratory of Protein and Nucleic Acid Chemistry, The Rockefeller University, New York, New York 10065, USA. 2 Wellcome Trust Centre for Cell Biology, University of Edinburgh, Michael Swann Building, Max Born Crescent, Edinburgh EH9 3BF, UK. 3 Department of Cell Biology, Harvard Medical School, Boston, Massachusetts 02115, USA.
  • A Chromosome Level Genome of Astyanax Mexicanus Surface Fish for Comparing Population

    A Chromosome Level Genome of Astyanax Mexicanus Surface Fish for Comparing Population

    bioRxiv preprint doi: https://doi.org/10.1101/2020.07.06.189654; this version posted July 6, 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. 1 Title 2 A chromosome level genome of Astyanax mexicanus surface fish for comparing population- 3 specific genetic differences contributing to trait evolution. 4 5 Authors 6 Wesley C. Warren1, Tyler E. Boggs2, Richard Borowsky3, Brian M. Carlson4, Estephany 7 Ferrufino5, Joshua B. Gross2, LaDeana Hillier6, Zhilian Hu7, Alex C. Keene8, Alexander Kenzior9, 8 Johanna E. Kowalko5, Chad Tomlinson10, Milinn Kremitzki10, Madeleine E. Lemieux11, Tina 9 Graves-Lindsay10, Suzanne E. McGaugh12, Jeff T. Miller12, Mathilda Mommersteeg7, Rachel L. 10 Moran12, Robert Peuß9, Edward Rice1, Misty R. Riddle13, Itzel Sifuentes-Romero5, Bethany A. 11 Stanhope5,8, Clifford J. Tabin13, Sunishka Thakur5, Yamamoto Yoshiyuki14, Nicolas Rohner9,15 12 13 Authors for correspondence: Wesley C. Warren ([email protected]), Nicolas Rohner 14 ([email protected]) 15 16 Affiliation 17 1Department of Animal Sciences, Department of Surgery, Institute for Data Science and 18 Informatics, University of Missouri, Bond Life Sciences Center, Columbia, MO 19 2 Department of Biological Sciences, University of Cincinnati, Cincinnati, OH 20 3 Department of Biology, New York University, New York, NY 21 4 Department of Biology, The College of Wooster, Wooster, OH 22 5 Harriet L. Wilkes Honors College, Florida Atlantic University, Jupiter FL 23 6 Department of Genome Sciences, University of Washington, Seattle, WA 1 bioRxiv preprint doi: https://doi.org/10.1101/2020.07.06.189654; this version posted July 6, 2020.
  • Evolutionary Fate of Retroposed Gene Copies in the Human Genome

    Evolutionary Fate of Retroposed Gene Copies in the Human Genome

    Evolutionary fate of retroposed gene copies in the human genome Nicolas Vinckenbosch*, Isabelle Dupanloup*†, and Henrik Kaessmann*‡ *Center for Integrative Genomics, University of Lausanne, Ge´nopode, 1015 Lausanne, Switzerland; and †Computational and Molecular Population Genetics Laboratory, Zoological Institute, University of Bern, 3012 Bern, Switzerland Communicated by Wen-Hsiung Li, University of Chicago, Chicago, IL, December 30, 2005 (received for review December 14, 2005) Given that retroposed copies of genes are presumed to lack the and rodent genomes (7–12). In addition, three recent studies regulatory elements required for their expression, retroposition using EST data (13, 14) and tiling-microarray data from chro- has long been considered a mechanism without functional rele- mosome 22 (15) indicated that retrocopy transcription may be vance. However, through an in silico assay for transcriptional widespread, although these surveys were limited, and potential activity, we identify here >1,000 transcribed retrocopies in the functional implications were not addressed. human genome, of which at least Ϸ120 have evolved into bona To further explore the functional significance of retroposition fide genes. Among these, Ϸ50 retrogenes have evolved functions in the human genome, we systematically screened for signatures in testes, more than half of which were recruited as functional of selection related to retrocopy transcription. Our results autosomal counterparts of X-linked genes during spermatogene- suggest that retrocopy transcription is not rare and that the sis. Generally, retrogenes emerge ‘‘out of the testis,’’ because they pattern of transcription of human retrocopies has been pro- are often initially transcribed in testis and later evolve stronger and foundly shaped by natural selection, acting both for and against sometimes more diverse spatial expression patterns.
  • A Unique Retrogene Within a Gene That Has Acquired Multiple Promoters and a Specific Function in Spermatogenesis

    A Unique Retrogene Within a Gene That Has Acquired Multiple Promoters and a Specific Function in Spermatogenesis

    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector Developmental Biology 304 (2007) 848–859 www.elsevier.com/locate/ydbio Genomes & Developmental Control Utp14b: A unique retrogene within a gene that has acquired multiple promoters and a specific function in spermatogenesis Ming Zhao a,1, Jan Rohozinski b,1, Manju Sharma c, Jun Ju a, Robert E. Braun c, ⁎ Colin E. Bishop b,d, Marvin L. Meistrich a, a Department of Experimental Radiation Oncology, University of Texas M.D. Anderson Cancer Center, Box 066, 1515 Holcombe Blvd, Houston, TX 77030, USA b Department of Obstetrics and Gynecology, Baylor College of Medicine, 1709 Dryden Road, Houston, TX 77030, USA c Department of Genome Sciences, University of Washington School of Medicine, Box 357730, 1705 N.E. Pacific Street, Seattle, WA 98195, USA d Department of Molecular and Human Genetics, Baylor College of Medicine, One Baylor Plaza, Houston, TX 77030, USA Received for publication 28 September 2006; revised 9 December 2006; accepted 3 January 2007 Available online 9 January 2007 Abstract The mouse retrogene Utp14b is essential for male fertility, and a mutation in its sequence results in the sterile juvenile spermatogonial depletion (jsd) phenotype. It is a retrotransposed copy of the Utp14a gene, which is located on the X chromosome, and is inserted within an intron of the autosomal acyl-CoA synthetase long-chain family member 3 (Acsl3) gene. To elucidate the roles of the Utp14 genes in normal spermatogenic cell development as a basis for understanding the defects that result in the jsd phenotype, we analyzed the various mRNAs produced from the Utp14b retrogene and their expression in different cell types.
  • WO 2012/174282 A2 20 December 2012 (20.12.2012) P O P C T

    WO 2012/174282 A2 20 December 2012 (20.12.2012) P O P C T

    (12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) (19) World Intellectual Property Organization International Bureau (10) International Publication Number (43) International Publication Date WO 2012/174282 A2 20 December 2012 (20.12.2012) P O P C T (51) International Patent Classification: David [US/US]; 13539 N . 95th Way, Scottsdale, AZ C12Q 1/68 (2006.01) 85260 (US). (21) International Application Number: (74) Agent: AKHAVAN, Ramin; Caris Science, Inc., 6655 N . PCT/US20 12/0425 19 Macarthur Blvd., Irving, TX 75039 (US). (22) International Filing Date: (81) Designated States (unless otherwise indicated, for every 14 June 2012 (14.06.2012) kind of national protection available): AE, AG, AL, AM, AO, AT, AU, AZ, BA, BB, BG, BH, BR, BW, BY, BZ, English (25) Filing Language: CA, CH, CL, CN, CO, CR, CU, CZ, DE, DK, DM, DO, Publication Language: English DZ, EC, EE, EG, ES, FI, GB, GD, GE, GH, GM, GT, HN, HR, HU, ID, IL, IN, IS, JP, KE, KG, KM, KN, KP, KR, (30) Priority Data: KZ, LA, LC, LK, LR, LS, LT, LU, LY, MA, MD, ME, 61/497,895 16 June 201 1 (16.06.201 1) US MG, MK, MN, MW, MX, MY, MZ, NA, NG, NI, NO, NZ, 61/499,138 20 June 201 1 (20.06.201 1) US OM, PE, PG, PH, PL, PT, QA, RO, RS, RU, RW, SC, SD, 61/501,680 27 June 201 1 (27.06.201 1) u s SE, SG, SK, SL, SM, ST, SV, SY, TH, TJ, TM, TN, TR, 61/506,019 8 July 201 1(08.07.201 1) u s TT, TZ, UA, UG, US, UZ, VC, VN, ZA, ZM, ZW.