An integrated approach in gene-expression landscape profiling to identify housekeeping and tissue-specific genes in cattle
Peng Li A , Yun Zhu A , Xiaolong Kang A , Xingang Dan A , Yun Ma A and Yuangang Shi A BA College of Agriculture, Ningxia University, Helan Mountain West Road 489, 750021, Yinchuan, Ningxia Hui Autonomous Region, China.
B Corresponding author. Email: shi_yg@nxu.edu.cn
Animal Production Science 61(16) 1643-1651 https://doi.org/10.1071/AN20689
Submitted: 25 December 2020 Accepted: 7 July 2021 Published: 13 September 2021
Journal Compilation © CSIRO 2021 Open Access CC BY-NC-ND
Abstract
Context: High-throughput transcriptome sequencing (RNA-Seq) has been widely applied in cattle studies. Public databases such as the National Center for Biotechnology Information (NCBI) contain large collections of gene expression data from various cattle tissues that can be used in gene expression analysis research
Aims: This study was conducted to investigate patterns of transcriptome variation across tissues of cattle through large-scale identification of housekeeping genes (i.e. those crucial to maintaining basic cellular activity) and tissue-specific genes in cattle tissues.
Methods: Using data available in the NCBI Sequence Read Archive database, we analysed 1377 transcriptome data sequences from 60 bovine tissue types, identified tissue-specific and housekeeping genes, and set up a web-based bovine gene expression analysis tool.
Key results: We found 101 genes widely expressed in almost all tissue and screened out five housekeeping genes: RPL35A, eIF4A2, GAPDH, IPO5 and PAK2. Focusing on 12 major organs, we found 861 genes specifically expressing in these tissues. Furthermore, 187 significantly differentially expressed genes were found among six types of muscle tissues. All expression data were made available at our new website http://cattleExp.org, which can be freely accessed for future gene expression analyses.
Conclusions: The housekeeping genes and tissue-specific genes identified will provide more information for researchers studying gene expression in cattle.
Implications: The web-based cattle gene expression analysis tool will make it easy for researchers to access large public datasets. Users can easily access all publicly available RNA data and upload their own RNA-Seq data.
Key words: beef cattle, bioinformatics, differentially expressed genes, gene regulation, genomics, muscle tissue, public database, reverse genetics, RNA-Seq.
References
Ardlie KG, DeLuca DS, Segrè AV, Sullivan TJ, Young TR, Gelfand ET, et al (2015) The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science 348, 648–660.| The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans.Crossref | GoogleScholarGoogle Scholar |
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G (2000) Gene Ontology: tool for the unification of biology. Nature Genetics 25, 25
| Gene Ontology: tool for the unification of biology.Crossref | GoogleScholarGoogle Scholar | 10802651PubMed |
Barber RD, Harmer DW, Coleman RA, Clark BJ (2005) GAPDH as a housekeeping gene: analysis of GAPDH mRNA expression in a panel of 72 human tissues. Physiological Genomics 21, 389–395.
| GAPDH as a housekeeping gene: analysis of GAPDH mRNA expression in a panel of 72 human tissues.Crossref | GoogleScholarGoogle Scholar | 15769908PubMed |
Chen S, Zhou Y, Chen Y, Gu J (2018) Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890.
| Fastp: an ultra-fast all-in-one FASTQ preprocessor.Crossref | GoogleScholarGoogle Scholar | 30423086PubMed |
Curina A, Termanini A, Barozzi I, Prosperini E, Simonatto M, Polletti S, Silvola A, Soldi M, Austenaa L, Bonaldi T, Ghisletti S, Natoli G (2017) High constitutive activity of a broad panel of housekeeping and tissue-specific cis-regulatory elements depends on a subset of ETS proteins. Genes & Development 31, 399–412.
| High constitutive activity of a broad panel of housekeeping and tissue-specific cis-regulatory elements depends on a subset of ETS proteins.Crossref | GoogleScholarGoogle Scholar |
Edlinger L, Berger-Becvar A, Menzl I, Hoermann G, Greiner G, Grundschober E, Bago-Horvath Z, Al-Zoughbi W, Hoefler G, Brostjan C, Gille L, Moriggl R, Spittler A, Sexl V, Hoelbl-Kovacic A (2017) Expansion of BCR/ABL1+ cells requires PAK2 but not PAK1. British Journal of Haematology 179, 229–241.
| Expansion of BCR/ABL1+ cells requires PAK2 but not PAK1.Crossref | GoogleScholarGoogle Scholar | 28707321PubMed |
Eisenberg E, Levanon EY (2013) Human housekeeping genes, revisited. Trends in Genetics 29, 569–574.
| Human housekeeping genes, revisited.Crossref | GoogleScholarGoogle Scholar | 23810203PubMed |
Fang L, Sahana G, Su G, Yu Y, Zhang S, Lund MS, Sørensen P (2017) Integrating sequence-based GWAS and RNA-seq provides novel insights into the genetic basis of mastitis and milk production in dairy cattle. Scientific Reports 7, 45560
| Integrating sequence-based GWAS and RNA-seq provides novel insights into the genetic basis of mastitis and milk production in dairy cattle.Crossref | GoogleScholarGoogle Scholar | 28358110PubMed |
Kanehisa M, Goto S (2000) KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Research 28, 27–30.
| KEGG: Kyoto encyclopedia of genes and genomes.Crossref | GoogleScholarGoogle Scholar | 10592173PubMed |
Khansefid M, Millen CA, Chen Y, Pryce JE, Chamberlain AJ, Vander Jagt CJ, Gondro C, Goddard ME (2017) Gene expression analysis of blood, liver, and muscle in cattle divergently selected for high and low residual feed intake. Journal of Animal Science 95, 4764–4775.
| Gene expression analysis of blood, liver, and muscle in cattle divergently selected for high and low residual feed intake.Crossref | GoogleScholarGoogle Scholar | 29293712PubMed |
Lee JS, Mo Y, Gan H, Burgess RJ, Baker DJ, van Deursen JM, Zhang Z (2019) Pak2 kinase promotes cellular senescence and organismal aging. Proceedings of the National Academy of Sciences of the United States of America 116, 13311–13319.
| Pak2 kinase promotes cellular senescence and organismal aging.Crossref | GoogleScholarGoogle Scholar | 31209047PubMed |
Mukiibi R, Vinsky M, Keogh K, Fitzsimmons C, Stothard P, Waters SM, Li C (2019) Liver transcriptome profiling of beef steers with divergent growth rate, feed intake, or metabolic body weight phenotypes. Journal of Animal Science 97, 4386–4404.
| Liver transcriptome profiling of beef steers with divergent growth rate, feed intake, or metabolic body weight phenotypes.Crossref | GoogleScholarGoogle Scholar | 31583405PubMed |
Ramírez-González RH, Borrill P, Lang D, Harrington SA, Brinton J, Venturini L, Davey M, Jacobs J, Van Ex F, Pasha A, Khedikar Y, Robinson SJ, Cory AT, Florio T, Concia L, Juery C, Schoonbeek H, Steuernagel B, Xiang D, Ridout CJ, Chalhoub B, Mayer KFX, Benhamed M, Latrasse D, Bendahmane A, Wulff BBH, Appels R, Tiwari V, Datla R, Choulet F, Pozniak CJ, Provart NJ, Sharpe AG, Paux E, Spannagl M, Bräutigam A, Uauy C (2018) The transcriptional landscape of polyploid wheat. Science 361, eaar6089
| The transcriptional landscape of polyploid wheat.Crossref | GoogleScholarGoogle Scholar | 30115783PubMed |
Robinson MD, Oshlack A (2010) A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biology 11, R25
| A scaling normalization method for differential expression analysis of RNA-seq data.Crossref | GoogleScholarGoogle Scholar | 20196867PubMed |
Ruan H, Li Y, Wang X, Sun B, Fang W, Jiang S, Liang C (2020) CRYAB inhibits migration and invasion of bladder cancer cells through the PI3K/AKT and ERK pathways. Japanese Journal of Clinical Oncology 50, 254–260.
| CRYAB inhibits migration and invasion of bladder cancer cells through the PI3K/AKT and ERK pathways.Crossref | GoogleScholarGoogle Scholar | 31829429PubMed |
Rubie C, Kempf K, Hans J, Su T, Tilton B, Georg T, Brittner B, Ludwig B, Schilling M (2005) Housekeeping gene variability in normal and cancerous colorectal, pancreatic, esophageal, gastric and hepatic tissues. Molecular and Cellular Probes 19, 101–109.
| Housekeeping gene variability in normal and cancerous colorectal, pancreatic, esophageal, gastric and hepatic tissues.Crossref | GoogleScholarGoogle Scholar | 15680211PubMed |
Sultan M, Schulz MH, Richard H, Magen A, Klingenhoff A, Scherf M, Seifert M, Borodina T, Soldatov A, Parkhomchuk D, Schmidt D, O’Keeffe S, Haas S, Vingron M, Lehrach H, Yaspo ML (2008) A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science 321, 956–960.
| A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome.Crossref | GoogleScholarGoogle Scholar | 18599741PubMed |
Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennan G, Grisar T, Igout A, Heinen E (1999) Housekeeping genes as internal standards: use and. Journal of Biotechnology 75, 291–295.
| Housekeeping genes as internal standards: use and.Crossref | GoogleScholarGoogle Scholar | 10617337PubMed |
Vandesompele J, Preter KD, Patty NF, Poppe B, Speleman F (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology 3, research0034.1
| Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes.Crossref | GoogleScholarGoogle Scholar | 12184808PubMed |
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB (2008a) Alternative isoform regulation in human tissue transcriptomes. Nature 456, 470–476.
| Alternative isoform regulation in human tissue transcriptomes.Crossref | GoogleScholarGoogle Scholar | 18978772PubMed |
Wang X, Sun Q, McGrath SD, Mardis ER, Soloway PD, Clark AG (2008b) Transcriptome-wide identification of novel imprinted genes in neonatal mouse brain. PLoS One 3, e3839
| Transcriptome-wide identification of novel imprinted genes in neonatal mouse brain.Crossref | GoogleScholarGoogle Scholar | 19066630PubMed |
Weijers SR, De Jonge J, Van Zanten O, Benedetti L, Langeveld J, Menkveld HW, Van Nieuwenhuijzen AF (2012) KALLISTO: cost effective and integrated optimization of the urban wastewater system Eindhoven. Water Practice and Technology 7, wpt2012036
| KALLISTO: cost effective and integrated optimization of the urban wastewater system Eindhoven.Crossref | GoogleScholarGoogle Scholar |
Wu J, Mao X, Cai T, Luo J, Wei L (2006) KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Research 34, W720–W724.
| KOBAS server: a web-based platform for automated annotation and pathway identification.Crossref | GoogleScholarGoogle Scholar | 16845106PubMed |
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li C-Y, Wei L (2011) KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Research 39, W316–W322.
| KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases.Crossref | GoogleScholarGoogle Scholar | 21715386PubMed |
Zhu B, Xue F, Li G, Zhang C (2020) CRYAB promotes osteogenic differentiation of human bone marrow stem cells via stabilizing β-catenin and promoting the Wnt signalling. Cell Proliferation 53, 1–12.
| CRYAB promotes osteogenic differentiation of human bone marrow stem cells via stabilizing β-catenin and promoting the Wnt signalling.Crossref | GoogleScholarGoogle Scholar |