1. Academic Validation
  2. Global Analysis of Transcription Start Sites in the New Ovine Reference Genome ( Oar rambouillet v1.0)

Global Analysis of Transcription Start Sites in the New Ovine Reference Genome ( Oar rambouillet v1.0)

  • Front Genet. 2020 Oct 23;11:580580. doi: 10.3389/fgene.2020.580580.
Mazdak Salavati 1 2 Alex Caulton 3 4 Richard Clark 5 Iveta Gazova 1 6 Timothy P L Smith 7 Kim C Worley 8 Noelle E Cockett 9 Alan L Archibald 1 Shannon M Clarke 3 Brenda M Murdoch 10 Emily L Clark 1 2
Affiliations

Affiliations

  • 1 The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Edinburgh, United Kingdom.
  • 2 Centre for Tropical Livestock Genetics and Health (CTLGH), Roslin Institute, University of Edinburgh, Midlothian, United Kingdom.
  • 3 AgResearch, Invermay Agricultural Centre, Mosgiel, New Zealand.
  • 4 Genetics Otago, Department of Biochemistry, University of Otago, Dunedin, New Zealand.
  • 5 Genetics Core, Edinburgh Clinical Research Facility, The University of Edinburgh, Edinburgh, United Kingdom.
  • 6 MRC Human Genetics Unit, The University of Edinburgh, Edinburgh, United Kingdom.
  • 7 USDA, Agricultural Research Service, U.S. Meat Animal Research Center, Clay Center, NE, United States.
  • 8 Baylor College of Medicine, Houston, TX, United States.
  • 9 Department of Animal, Dairy and Veterinary Sciences, Utah State University, Logan, UT, United States.
  • 10 Department of Animal, Veterinary and Food Sciences, University of Idaho, Moscow, ID, United States.
Abstract

The overall aim of the Ovine FAANG project is to provide a comprehensive annotation of the new highly contiguous sheep reference genome sequence (Oar rambouillet v1.0). Mapping of transcription start sites (TSS) is a key first step in understanding transcript regulation and diversity. Using 56 tissue samples collected from the reference ewe Benz2616, we have performed a global analysis of TSS and TSS-Enhancer clusters using Cap Analysis Gene Expression (CAGE) sequencing. CAGE measures RNA expression by 5' cap-trapping and has been specifically designed to allow the characterization of TSS within promoters to single-nucleotide resolution. We have adapted an analysis pipeline that uses TagDust2 for clean-up and trimming, Bowtie2 for mapping, CAGEfightR for clustering, and the Integrative Genomics Viewer (IGV) for visualization. Mapping of CAGE tags indicated that the expression levels of CAGE tag clusters varied across tissues. Expression profiles across tissues were validated using corresponding polyA+ mRNA-Seq data from the same samples. After removal of CAGE tags with <10 read counts, 39.3% of TSS overlapped with 5' ends of 31,113 transcripts that had been previously annotated by NCBI (out of a total of 56,308 from the NCBI annotation). For 25,195 of the transcripts, previously annotated by NCBI, no TSS meeting stringent criteria were identified. A further 14.7% of TSS mapped to within 50 bp of annotated promoter regions. Intersecting these predicted TSS regions with annotated promoter regions (±50 bp) revealed 46% of the predicted TSS were "novel" and previously un-annotated. Using whole-genome bisulfite sequencing data from the same tissues, we were able to determine that a proportion of these "novel" TSS were hypo-methylated (32.2%) indicating that they are likely to be reproducible rather than "noise". This global analysis of TSS in sheep will significantly enhance the annotation of gene models in the new ovine reference assembly. Our analyses provide one of the highest resolution annotations of transcript regulation and diversity in a livestock species to date.

Keywords

CAGE; FAANG; TSS; WGBS; enhancer; ovine; promoter; transcriptome.

Figures