Investigation of de novo unique differentially expressed genes related to evolution in exercise response during domestication in Thoroughbred race horses.
Abstract: Previous studies of horse RNA-seq were performed by mapping sequence reads to the reference genome during transcriptome analysis. However in this study, we focused on two main ideas. First, differentially expressed genes (DEGs) were identified by de novo-based analysis (DBA) in RNA-seq data from six Thoroughbreds before and after exercise, here-after referred to as "de novo unique differentially expressed genes" (DUDEG). Second, by integrating both conventional DEGs and genes identified as being selected for during domestication of Thoroughbred and Jeju pony from whole genome re-sequencing (WGS) data, we give a new concept to the definition of DEG. We identified 1,034 and 567 DUDEGs in skeletal muscle and blood, respectively. DUDEGs in skeletal muscle were significantly related to exercise-induced stress biological process gene ontology (BP-GO) terms: 'immune system process'; 'response to stimulus'; and, 'death' and a KEGG pathways: 'JAK-STAT signaling pathway'; 'MAPK signaling pathway'; 'regulation of actin cytoskeleton'; and, 'p53 signaling pathway'. In addition, we found TIMELESS, EIF4A3 and ZNF592 in blood and CHMP4C and FOXO3 in skeletal muscle, to be in common between DUDEGs and selected genes identified by evolutionary statistics such as FST and Cross Population Extended Haplotype Homozygosity (XP-EHH). Moreover, in Thoroughbreds, three out of five genes (CHMP4C, EIF4A3 and FOXO3) related to exercise response showed relatively low nucleotide diversity compared to the Jeju pony. DUDEGs are not only conceptually new DEGs that cannot be attained from reference-based analysis (RBA) but also supports previous RBA results related to exercise in Thoroughbred. In summary, three exercise related genes which were selected for during domestication in the evolutionary history of Thoroughbred were identified as conceptually new DEGs in this study.
The Equine Research Bank provides access to a large database of publicly available scientific literature. Inclusion in the Research Bank does not imply endorsement of study methods or findings by Mad Barn.
This research summary has been generated with artificial intelligence and may contain errors and omissions. Refer to the original study to confirm details provided. Submit correction.
The research article discusses a study aimed at identifying unique genes that are differentially expressed during exercise in Thoroughbred race horses and could possibly be related to their evolution during domestication.
Introduction to the Study
The study focuses on two primary aspects. Firstly, it identifies differentially expressed genes (DEGs) using an independent, de novo-based analysis (DBA) in RNA-seq data from six Thoroughbreds before and after exercise. These are called “de novo unique differentially expressed genes” (DUDEG).
Secondly, the study combines the traditional DEGs and genes identified as being favored during domestication of Thoroughbred and Jeju pony using whole genome re-sequencing (WGS) data.
This approach brings a new perspective to how DEGs are defined.
Key Findings of the Study
Researchers identified 1,034 DUDEGs in skeletal muscle and 567 DUDEGs in blood.
The DUDEGs in skeletal muscle were significantly linked with certain exercise-induced stress biological process gene ontology (BP-GO) terms like ‘immune system process’, ‘response to stimulus’, and ‘death’ along with KEGG pathways like ‘JAK-STAT signaling pathway’, ‘MAPK signaling pathway’, ‘regulation of actin cytoskeleton’, and ‘p53 signaling pathway’.
Several genes, such as TIMELESS, EIF4A3, and ZNF592 in blood, and CHMP4C and FOXO3 in skeletal muscle, were found common between DUDEGs and selected genes identified by evolutionary statistics such as FST and Cross Population Extended Haplotype Homozygosity (XP-EHH).
In Thoroughbreds, three out of five exercise-related genes (CHMP4C, EIF4A3, and FOXO3) demonstrated comparatively lower nucleotide diversity than in the Jeju pony.
Study Conclusion
DUDEGs are conceptually new DEGs that cannot be obtained from reference-based analysis (RBA), but support previous RBA results related to exercise in Thoroughbred.
In summary, the study highlighted three exercise-related genes which were selected for during the domestication in the evolutionary history of Thoroughbred.
Cite This Article
APA
Park W, Kim J, Kim HJ, Choi J, Park JW, Cho HW, Kim BW, Park MH, Shin TS, Cho SK, Park JK, Kim H, Hwang JY, Lee CK, Lee HK, Cho S, Cho BW.
(2014).
Investigation of de novo unique differentially expressed genes related to evolution in exercise response during domestication in Thoroughbred race horses.
PLoS One, 9(3), e91418.
https://doi.org/10.1371/journal.pone.0091418
Department of Agricultural Biotechnology and Research Institute for Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea.
Kim, Jaemin
Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, Republic of Korea.
Kim, Hyeon Jeong
C&K genomics, Seoul National University, Seoul, Republic of Korea.
Choi, JaeYoung
Department of Animal Science, College of Life Sciences, Pusan National University, Miryang, Republic of Korea.
Park, Jeong-Woong
Department of Animal Science, College of Life Sciences, Pusan National University, Miryang, Republic of Korea.
Cho, Hyun-Woo
Department of Animal Science, College of Life Sciences, Pusan National University, Miryang, Republic of Korea.
Kim, Byeong-Woo
Department of Animal Science, College of Life Sciences, Pusan National University, Miryang, Republic of Korea.
Park, Myung Hum
TNT Research, Anyang, Republic of Korea.
Shin, Teak-Soon
Department of Animal Science, College of Life Sciences, Pusan National University, Miryang, Republic of Korea.
Cho, Seong-Keun
Department of Animal Science, College of Life Sciences, Pusan National University, Miryang, Republic of Korea.
Park, Jun-Kyu
Leaders in Industry-university Cooperation, Pusan National University, Miryang, Republic of Korea.
Kim, Heebal
Department of Agricultural Biotechnology and Research Institute for Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea; C&K genomics, Seoul National University, Seoul, Republic of Korea.
Hwang, Jae Yeon
Department of Agricultural Biotechnology and Research Institute for Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea.
Lee, Chang-Kyu
Department of Agricultural Biotechnology and Research Institute for Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea.
Lee, Hak-Kyo
Genomic Informatics Center, Hankyong National University, Anseong, Republic of Korea.
Cho, Seoae
C&K genomics, Seoul National University, Seoul, Republic of Korea.
Cho, Byung-Wook
Department of Animal Science, College of Life Sciences, Pusan National University, Miryang, Republic of Korea.
MeSH Terms
Animals
Gene Expression Regulation
Genome
Genotype
Horses / genetics
Horses / metabolism
Physical Conditioning, Animal
Transcriptome
Conflict of Interest Statement
MHP has the following financial competing interest: Paid employment in TNT Research Company Limited, Dongan-gu, Anyang-si, Republic of Korea. The rest of the authors have declared that no competing interests exist. This does not alter the authors' adherence to all the PLoS ONE policies on sharing data and materials. HK has the following financial competing interest: Paid employment in C&K Genomics company INC, C&K genomics Main Bldg. #514, SNU Research Park, Seoul National University Mt.4-2. SC has the following financial competing interest: CEO of C&K Genomics company INC, C&K genomics Main Bldg. #514, SNU Research Park, Seoul National University Mt.4-2.
References
This article includes 69 references
Weatherby J. An Introduction to a General Stud Book. .
Poole D. Current concepts of oxygen transport during exercise. Equine and Comparative Exercise Physiology 1: 5–22.
Wright SI, Bi IV, Schroeder SG, Yamasaki M, Doebley JF, McMullen MD, Gaut BS. The effects of artificial selection on the maize genome.. Science 2005 May 27;308(5726):1310-4.
Chen G, Li R, Shi L, Qi J, Hu P, Luo J, Liu M, Shi T. Revealing the missing expressed genes beyond the human reference genome by RNA-Seq.. BMC Genomics 2011 Dec 2;12:590.
Kinsella M, Harismendy O, Nakano M, Frazer KA, Bafna V. Sensitive gene fusion detection using ambiguously mapping RNA-Seq read pairs.. Bioinformatics 2011 Apr 15;27(8):1068-75.
Birzele F, Schaub J, Rust W, Clemens C, Baum P, Kaufmann H, Weith A, Schulz TW, Hildebrandt T. Into the unknown: expression profiling without genome sequence information in CHO by next generation sequencing.. Nucleic Acids Res 2010 Jul;38(12):3999-4010.
Burset M, Seledtsov IA, Solovyev VV. Analysis of canonical and non-canonical splice sites in mammalian genomes.. Nucleic Acids Res 2000 Nov 1;28(21):4364-75.
Robertson G, Schein J, Chiu R, Corbett R, Field M, Jackman SD, Mungall K, Lee S, Okada HM, Qian JQ, Griffith M, Raymond A, Thiessen N, Cezard T, Butterfield YS, Newsome R, Chan SK, She R, Varhol R, Kamoh B, Prabhu AL, Tam A, Zhao Y, Moore RA, Hirst M, Marra MA, Jones SJ, Hoodless PA, Birol I. De novo assembly and analysis of RNA-seq data.. Nat Methods 2010 Nov;7(11):909-12.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-Seq data without a reference genome.. Nat Biotechnol 2011 May 15;29(7):644-52.
Birol I, Jackman SD, Nielsen CB, Qian JQ, Varhol R, Stazyk G, Morin RD, Zhao Y, Hirst M, Schein JE, Horsman DE, Connors JM, Gascoyne RD, Marra MA, Jones SJ. De novo transcriptome assembly with ABySS.. Bioinformatics 2009 Nov 1;25(21):2872-7.
Schulz MH, Zerbino DR, Vingron M, Birney E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels.. Bioinformatics 2012 Apr 15;28(8):1086-92.
Martin J, Bruno VM, Fang Z, Meng X, Blow M, Zhang T, Sherlock G, Snyder M, Wang Z. Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads.. BMC Genomics 2010 Nov 24;11:663.
Surget-Groba Y, Montoya-Burgos JI. Optimization of de novo transcriptome assembly from next-generation sequencing data.. Genome Res 2010 Oct;20(10):1432-40.
Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J. SOAP2: an improved ultrafast tool for short read alignment.. Bioinformatics 2009 Aug 1;25(15):1966-7.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.. Bioinformatics 2010 Jan 1;26(1):139-40.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The Sequence Alignment/Map format and SAMtools.. Bioinformatics 2009 Aug 15;25(16):2078-9.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data.. Genome Res 2010 Sep;20(9):1297-303.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, McKenna A, Fennell TJ, Kernytsky AM, Sivachenko AY, Cibulskis K, Gabriel SB, Altshuler D, Daly MJ. A framework for variation discovery and genotyping using next-generation DNA sequencing data.. Nat Genet 2011 May;43(5):491-8.
Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, McVean G, Durbin R. The variant call format and VCFtools.. Bioinformatics 2011 Aug 1;27(15):2156-8.
Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.. Nat Protoc 2009;4(1):44-57.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method.. Methods 2001 Dec;25(4):402-8.
Browning BL, Yu Z. Simultaneous genotype calling and haplotype phasing improves genotype accuracy and reduces false-positive associations for genome-wide association studies.. Am J Hum Genet 2009 Dec;85(6):847-61.
Reynolds J, Weir BS, Cockerham CC. Estimation of the coancestry coefficient: basis for a short-term genetic distance.. Genetics 1983 Nov;105(3):767-79.
Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows.. Mol Ecol Resour 2010 May;10(3):564-7.
Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, Xie X, Byrne EH, McCarroll SA, Gaudet R, Schaffner SF, Lander ES, Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, Pasternak S, Wheeler DA, Willis TD, Yu F, Yang H, Zeng C, Gao Y, Hu H, Hu W, Li C, Lin W, Liu S, Pan H, Tang X, Wang J, Wang W, Yu J, Zhang B, Zhang Q, Zhao H, Zhao H, Zhou J, Gabriel SB, Barry R, Blumenstiel B, Camargo A, Defelice M, Faggart M, Goyette M, Gupta S, Moore J, Nguyen H, Onofrio RC, Parkin M, Roy J, Stahl E, Winchester E, Ziaugra L, Altshuler D, Shen Y, Yao Z, Huang W, Chu X, He Y, Jin L, Liu Y, Shen Y, Sun W, Wang H, Wang Y, Wang Y, Xiong X, Xu L, Waye MM, Tsui SK, Xue H, Wong JT, Galver LM, Fan JB, Gunderson K, Murray SS, Oliphant AR, Chee MS, Montpetit A, Chagnon F, Ferretti V, Leboeuf M, Olivier JF, Phillips MS, Roumy S, Sallée C, Verner A, Hudson TJ, Kwok PY, Cai D, Koboldt DC, Miller RD, Pawlikowska L, Taillon-Miller P, Xiao M, Tsui LC, Mak W, Song YQ, Tam PK, Nakamura Y, Kawaguchi T, Kitamoto T, Morizono T, Nagashima A, Ohnishi Y, Sekine A, Tanaka T, Tsunoda T, Deloukas P, Bird CP, Delgado M, Dermitzakis ET, Gwilliam R, Hunt S, Morrison J, Powell D, Stranger BE, Whittaker P, Bentley DR, Daly MJ, de Bakker PI, Barrett J, Chretien YR, Maller J, McCarroll S, Patterson N, Pe'er I, Price A, Purcell S, Richter DJ, Sabeti P, Saxena R, Schaffner SF, Sham PC, Varilly P, Altshuler D, Stein LD, Krishnan L, Smith AV, Tello-Ruiz MK, Thorisson GA, Chakravarti A, Chen PE, Cutler DJ, Kashuk CS, Lin S, Abecasis GR, Guan W, Li Y, Munro HM, Qin ZS, Thomas DJ, McVean G, Auton A, Bottolo L, Cardin N, Eyheramendy S, Freeman C, Marchini J, Myers S, Spencer C, Stephens M, Donnelly P, Cardon LR, Clarke G, Evans DM, Morris AP, Weir BS, Tsunoda T, Johnson TA, Mullikin JC, Sherry ST, Feolo M, Skol A, Zhang H, Zeng C, Zhao H, Matsuda I, Fukushima Y, Macer DR, Suda E, Rotimi CN, Adebamowo CA, Ajayi I, Aniagwu T, Marshall PA, Nkwodimmah C, Royal CD, Leppert MF, Dixon M, Peiffer A, Qiu R, Kent A, Kato K, Niikawa N, Adewole IF, Knoppers BM, Foster MW, Clayton EW, Watkin J, Gibbs RA, Belmont JW, Muzny D, Nazareth L, Sodergren E, Weinstock GM, Wheeler DA, Yakub I, Gabriel SB, Onofrio RC, Richter DJ, Ziaugra L, Birren BW, Daly MJ, Altshuler D, Wilson RK, Fulton LL, Rogers J, Burton J, Carter NP, Clee CM, Griffiths M, Jones MC, McLay K, Plumb RW, Ross MT, Sims SK, Willey DL, Chen Z, Han H, Kang L, Godbout M, Wallenburg JC, L'Archevêque P, Bellemare G, Saeki K, Wang H, An D, Fu H, Li Q, Wang Z, Wang R, Holden AL, Brooks LD, McEwen JE, Guyer MS, Wang VO, Peterson JL, Shi M, Spiegel J, Sung LM, Zacharia LF, Collins FS, Kennedy K, Jamieson R, Stewart J. Genome-wide detection and characterization of positive selection in human populations.. Nature 2007 Oct 18;449(7164):913-8.
Peng Z, Cheng Y, Tan BC, Kang L, Tian Z, Zhu Y, Zhang W, Liang Y, Hu X, Tan X, Guo J, Dong Z, Liang Y, Bao L, Wang J. Comprehensive analysis of RNA-Seq data reveals extensive RNA editing in a human transcriptome.. Nat Biotechnol 2012 Feb 12;30(3):253-60.
Garber M, Grabherr MG, Guttman M, Trapnell C. Computational methods for transcriptome annotation and quantification using RNA-seq.. Nat Methods 2011 Jun;8(6):469-77.
Bray NJ, Buckland PR, Owen MJ, O'Donovan MC. Cis-acting variation in the expression of a high proportion of genes in human brain.. Hum Genet 2003 Jul;113(2):149-53.
Gregg C, Zhang J, Weissbourd B, Luo S, Schroth GP, Haig D, Dulac C. High-resolution analysis of parent-of-origin allelic expression in the mouse brain.. Science 2010 Aug 6;329(5992):643-8.
Xia Z, Xu H, Zhai J, Li D, Luo H, He C, Huang X. RNA-Seq analysis and de novo transcriptome assembly of Hevea brasiliensis.. Plant Mol Biol 2011 Oct;77(3):299-308.
Kingston SG, Hoffman-Goetz L. Effect of environmental enrichment and housing density on immune system reactivity to acute exercise stress.. Physiol Behav 1996 Jul;60(1):145-50.
Stefansson S, Yepes M, Gorlatova N, Day DE, Moore EG, Zabaleta A, McMahon GA, Lawrence DA. Mutants of plasminogen activator inhibitor-1 designed to inhibit neutrophil elastase and cathepsin G are more effective in vivo than their endogenous inhibitors.. J Biol Chem 2004 Jul 16;279(29):29981-7.
Dousset E, Avela J, Ishikawa M, Kallio J, Kuitunen S, Kyröláinen H, Linnamo V, Komi PV. Bimodal recovery pattern in human skeletal muscle induced by exhaustive stretch-shortening cycle exercise.. Med Sci Sports Exerc 2007 Mar;39(3):453-60.
Sun L, Ma K, Wang H, Xiao F, Gao Y, Zhang W, Wang K, Gao X, Ip N, Wu Z. JAK1-STAT1-STAT3, a key pathway promoting proliferation and preventing premature differentiation of myoblasts.. J Cell Biol 2007 Oct 8;179(1):129-38.
O'Shea JJ, Pesu M, Borie DC, Changelian PS. A new modality for immunosuppression: targeting the JAK/STAT pathway.. Nat Rev Drug Discov 2004 Jul;3(7):555-64.
Long YC, Widegren U, Zierath JR. Exercise-induced mitogen-activated protein kinase signalling in skeletal muscle.. Proc Nutr Soc 2004 May;63(2):227-32.
Cecchetti L, Tolley ND, Michetti N, Bury L, Weyrich AS, Gresele P. Megakaryocytes differentially sort mRNAs for matrix metalloproteinases and their inhibitors into platelets: a mechanism for regulating synthetic events.. Blood 2011 Aug 18;118(7):1903-11.
Giorgi C, Yeo GW, Stone ME, Katz DB, Burge C, Turrigiano G, Moore MJ. The EJC factor eIF4AIII modulates synaptic strength and neuronal protein expression.. Cell 2007 Jul 13;130(1):179-91.
Sandri M, Lin J, Handschin C, Yang W, Arany ZP, Lecker SH, Goldberg AL, Spiegelman BM. PGC-1alpha protects skeletal muscle from atrophy by suppressing FoxO3 action and atrophy-specific gene transcription.. Proc Natl Acad Sci U S A 2006 Oct 31;103(44):16260-5.
Zhao J, Brault JJ, Schild A, Cao P, Sandri M, Schiaffino S, Lecker SH, Goldberg AL. FoxO3 coordinately activates protein degradation by the autophagic/lysosomal and proteasomal pathways in atrophying muscle cells.. Cell Metab 2007 Dec;6(6):472-83.
Jamart C, Benoit N, Raymackers JM, Kim HJ, Kim CK, Francaux M. Autophagy-related and autophagy-regulatory genes are induced in human muscle after ultraendurance exercise.. Eur J Appl Physiol 2012 Aug;112(8):3173-7.
Arora D, Srikanth K, Lee J, Lee D, Park N, Wy S, Kim H, Park JE, Chai HH, Lim D, Cho IC, Kim J, Park W. Integration of multi-omics approaches for functional characterization of muscle related selective sweep genes in Nanchukmacdon. Sci Rep 2021 Mar 30;11(1):7219.
Bicskei B, Taggart JB, Bron JE, Glover KA. Transcriptomic comparison of communally reared wild, domesticated and hybrid Atlantic salmon fry under stress and control conditions. BMC Genet 2020 May 29;21(1):57.
Mattucci F, Galaverni M, Lyons LA, Alves PC, Randi E, Velli E, Pagani L, Caniglia R. Genomic approaches to identify hybrids and estimate admixture times in European wildcat populations. Sci Rep 2019 Aug 12;9(1):11612.
Bicskei B, Taggart JB, Glover KA, Bron JE. Comparing the transcriptomes of embryos from domesticated and wild Atlantic salmon (Salmo salar L.) stocks and examining factors that influence heritability of gene expression. Genet Sel Evol 2016 Mar 17;48:20.
Kim H, Caetano-Anolles K, Seo M, Kwon YJ, Cho S, Seo K, Kim H. Prediction of Genes Related to Positive Selection Using Whole-Genome Resequencing in Three Commercial Pig Breeds. Genomics Inform 2015 Dec;13(4):137-45.