Weighted single-step GWAS and RNA sequencing reveals key candidate genes associated with physiological indicators of heat stress in Holstein cattle
Tóm tắt
The study of molecular processes regulating heat stress response in dairy cattle is paramount for developing mitigation strategies to improve heat tolerance and animal welfare. Therefore, we aimed to identify quantitative trait loci (QTL) regions associated with three physiological indicators of heat stress response in Holstein cattle, including rectal temperature (RT), respiration rate score (RS), and drooling score (DS). We estimated genetic parameters for all three traits. Subsequently, a weighted single-step genome-wide association study (WssGWAS) was performed based on 3200 genotypes, 151,486 phenotypic records, and 38,101 animals in the pedigree file. The candidate genes located within the identified QTL regions were further investigated through RNA sequencing (RNA-seq) analyses of blood samples for four cows collected in April (non-heat stress group) and four cows collected in July (heat stress group). The heritability estimates for RT, RS, and DS were 0.06, 0.04, and 0.03, respectively. Fourteen, 19, and 20 genomic regions explained 2.94%, 3.74%, and 4.01% of the total additive genetic variance of RT, RS, and DS, respectively. Most of these genomic regions are located in the Bos taurus autosome (BTA) BTA3, BTA6, BTA8, BTA12, BTA14, BTA21, and BTA24. No genomic regions overlapped between the three indicators of heat stress, indicating the polygenic nature of heat tolerance and the complementary mechanisms involved in heat stress response. For the RNA-seq analyses, 2627 genes were significantly upregulated and 369 downregulated in the heat stress group in comparison to the control group. When integrating the WssGWAS, RNA-seq results, and existing literature, the key candidate genes associated with physiological indicators of heat stress in Holstein cattle are: PMAIP1, SBK1, TMEM33, GATB, CHORDC1, RTN4IP1, and BTBD7. Physiological indicators of heat stress are heritable and can be improved through direct selection. Fifty-three QTL regions associated with heat stress indicators confirm the polygenic nature and complex genetic determinism of heat tolerance in dairy cattle. The identified candidate genes will contribute for optimizing genomic evaluation models by assigning higher weights to genetic markers located in these regions as well as to the design of SNP panels containing polymorphisms located within these candidate genes.
Tài liệu tham khảo
Schär C, Vidale PL, Lüthi D, Frei C, Häberli C, Liniger MA, et al. The role of increasing temperature variability in European summer heatwaves. Nature. 2004;427(6972):332–6. https://doi.org/10.1038/nature02300.
Das R, Sailo L, Verma N, Bharti P, Saikia J, Imtiwati KR. Impact of heat stress on health and performance of dairy animals: a review. Vet World. 2016;9(3):260–8. https://doi.org/10.14202/vetworld.2016.260-268.
Polsky L, von Keyserlingk M. Invited review: effects of heat stress on dairy cattle welfare. J Dairy Sci. 2017;100(11):8645–57. https://doi.org/10.3168/jds.2017-12651.
Sejian V, Bhatta R, Gaughan JB, Dunshea FR, Lacetera N. Review: adaptation of animals to heat stress. Animal. 2018;12(s2):s431–44. https://doi.org/10.1017/S1751731118001945.
Ju XH, Xu HJ, Yong YH, An LL, Jiao PR, Liao M. Heat stress upregulation of toll-like receptors 2/4 and acute inflammatory cytokines in peripheral blood mononuclear cell (PBMC) of Bama miniature pigs: An in vivo and in vitro study. Animal. 2014;8(9):1462–8. https://doi.org/10.1017/S1751731114001268.
Jin Y, Hu Y, Han D, Wang M. Chronic heat stress weakened the innate immunity and increased the virulence of highly pathogenic avian influenza virus H5N1 in mice. J Biomed Biotechnol. 2011;2011:367846. https://doi.org/10.1155/2011/367846.
Bagath M, Krishnan G, Devaraj C, Rashamol VP, Pragna P, Lees AM, et al. The impact of heat stress on the immune system in dairy cattle: a review. Res Vet Sci. 2019;126:94–102. https://doi.org/10.1016/j.rvsc.2019.08.011.
Sammad A, Wang YJ, Umer S, Lirong H, Khan I, Khan A, et al. Nutritional physiology and biochemistry of dairy cattle under the influence of heat stress: consequences and opportunities. Animals. 2020;10(5):793. https://doi.org/10.3390/ani10050793.
Parkhurst AM. Model for understanding thermal hysteresis during heat stress: a matter of direction. Int J Biometeorol. 2010;54(6):637–45. https://doi.org/10.1007/s00484-009-0299-z.
Luo H, Brito LF, Li X, Su G, Dou J, Xu W, et al. Genetic parameters for rectal temperature, respiration rate, and drooling score in Holstein cattle and their relationships with various fertility, production, body conformation, and health traits. J Dairy Sci. 2021;104(4):4390–403. https://doi.org/10.3168/jds.2020-19192.
Chen Z, Brito LF, Luo H, Shi R, Chang Y, Liu L, et al. Genetic and genomic analyses of service sire effect on female reproductive traits in Holstein cattle. Front Genet. 2021;12:713575. https://doi.org/10.3389/fgene.2021.713575.
Luo H, Li X, Hu L, Xu W, Chu Q, Liu A, et al. Genomic analyses and biological validation of candidate genes for rectal temperature as an indicator of heat stress in Holstein cattle. J Dairy Sci. 2021;104(4):4441–51. https://doi.org/10.3168/jds.2020-18725.
Wang H, Misztal I, Aguilar I, Legarra A, Muir W. Genome-wide association mapping including phenotypes from relatives without genotypes. Genet Res. 2012;94:73–83. https://doi.org/10.1017/S0016672312000274.
Freitas PHF, Oliveira HR, Silva FF, Fleming A, Miglior F, Schenkel FS, et al. Genomic analyses for predicted Milk fatty acid composition throughout lactation in north American Holstein cattle. J Dairy Sci. 2020;103(7):6318–31. https://doi.org/10.3168/jds.2019-17628.
Oliveira HR, Lourenco DAL, Masuda Y, Misztal I, Tsuruta S, Jamrozik J, et al. Single-step genome-wide association for longitudinal traits of Canadian Ayrshire, Holstein, and Jersey dairy cattle. J Dairy Sci. 2019;102(11):9995–10011. https://doi.org/10.3168/jds.2019-16821.
Guarini AR, Lourenco DAL, Brito LF, Sargolzaei M, Baes CF, Miglior F, et al. Genetics and genomics of reproductive disorders in Canadian Holstein cattle. J Dairy Sci. 2019;102(2):1341–53. https://doi.org/10.3168/jds.2018-15038.
Shi R, Brito LF, Liu A, Luo H, Chen Z, Liu L, et al. Genotype-by-environment interaction in Holstein heifer fertility traits using single-step genomic reaction norm models. BMC Genomics. 2021;22:193. https://doi.org/10.1186/s12864-021-07496-3.
Zhao B, Luo H, Huang X, Wei C, Di J, Tian Y, et al. Integration of a single-step genome-wide association study with a multi-tissue transcriptome analysis provides novel insights into the genetic basis of wool and weight traits in sheep. Genet Sel Evol. 2021;53:56. https://doi.org/10.1186/s12711-021-00649-8.
Ozsolak F, Milos PM. RNA sequencing: advances, challenges and opportunities. Nat Rev Genet. 2011;12(2):87–98. https://doi.org/10.1038/nrg2934.
Browning BL, Zhou Y, Browning SR. A one-penny imputed genome from next-generation reference panels. Am J Hum Genet. 2018;103(3):338–48. https://doi.org/10.1016/j.ajhg.2018.07.015.
Aguilar I, Misztal I, Johnson DL, Legarra A, Tsuruta S, Lawlor TJ. Hot topic: a unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. J Dairy Sci. 2010;93(2):743–52. https://doi.org/10.3168/jds.2009-2730.
VanRaden PM. Efficient methods to compute genomic predictions. J Dairy Sci. 2008;91(11):4414–23. https://doi.org/10.3168/jds.2007-0980.
Misztal I, Tsuruta S, Lourenco DAL, Masuda Y, Aguilar I, Legarra A, et al. Manual for BLUPF90 family programs: University of Georgia; 2018. http://nce.ads.uga.edu/wiki/doku.php?id=documentation.
Wang H, Misztal I, Aguilar I, Legarra A, Fernando RL, Vitezica Z, et al. Genome-wide association mapping including phenotypes from relatives without genotypes in a single-step (ssGWAS) for 6-week body weight in broiler chickens. Front Genet. 2014;5:134. https://doi.org/10.3389/fgene.2014.00134.
Team RC. R: A Language and environment for statistical computing. R Foundation for statistical computing, Vienna, Austria. 2021. URL https://www.r-project.org.
Durinck S, Spellman PT, Birney E, Huber W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protoc. 2009;4(8):1184–91. https://doi.org/10.1038/nprot.2009.97.
Yu G, Wang L, Han Y, He Q. ClusterProfiler: An R package for comparing biological themes among gene clusters. OMICS: J Integr Biol. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Liu G, Liao Y, Sun B, Guo Y, Deng M, Li Y, et al. Effects of chronic heat stress on mRNA and miRNA expressions in dairy cows. Gene. 2020;742:144550. https://doi.org/10.1016/j.gene.2020.144550.
Yue S, Wang Z, Wang L, Peng Q, Xue B. Transcriptome functional analysis of mammary gland of cows in heat stress and thermoneutral condition. Animals. 2020;10(6):1015.
Rio DC, Ares MJ, Hannon GJ, Nilsen TW. Purification of RNA using TRIzol (TRI reagent). Cold Spring Harb Protoc. 2010;2010(6):t5439. https://doi.org/10.1101/pdb.prot5439.
Chen S, Zhou Y, Chen Y, Gu J. Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90. https://doi.org/10.1093/bioinformatics/bty560.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21. https://doi.org/10.1093/bioinformatics/bts635.
Picard. URL http://broadinstitute.github.io/picard/. Accessed 9 Dec 2020.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75. https://doi.org/10.1086/519795.
McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40(10):4288–97. https://doi.org/10.1093/nar/gks042.
Raivo Kolde. Pheatmap: pretty heatmaps. 2019. URL https://CRAN.R-project.org/package=pheatmap.
Dikmen S, Cole JB, Null DJ, Hansen PJ. Heritability of rectal temperature and genetic correlations with production and reproduction traits in dairy cattle. J Dairy Sci. 2012;95(6):3401–5. https://doi.org/10.3168/jds.2011-4306.
Riley DG, Chase CC, Coleman SW, Olson TA. Genetic assessment of rectal temperature and coat score in Brahman, Angus, and Romosinuano crossbred and straightbred cows and calves under subtropical summer conditions. Livest Sci. 2012;148(1–2):109–18. https://doi.org/10.1016/j.livsci.2012.05.017.
Muranty H, Troggio M, Sadok IB, Rifaï MA, Auwerkerken A, Banchi E, et al. Accuracy and responses of genomic selection on key traits in apple breeding. Hortic Res-England. 2015;2:15060. https://doi.org/10.1038/hortres.2015.60.
Song H, Zhang J, Zhang Q, Ding X. Using different single-step strategies to improve the efficiency of genomic prediction on body measurement traits in pig. Front Genet. 2019;9:730. https://doi.org/10.3389/fgene.2018.00730.
Wiggans GR, Cooper TA, VanRaden PM, Cole JB. Technical note: adjustment of traditional cow evaluations to improve accuracy of genomic predictions. J Dairy Sci. 2011;94(12):6188–93. https://doi.org/10.3168/jds.2011-4481.
Mancin E, Lourenco D, Bermann M, Mantovani R, Misztal I. Accounting for population structure and phenotypes from relatives in association mapping for farm animals: a simulation study. Front Genet. 2021;12:642065. https://doi.org/10.3389/fgene.2021.642065.
Marques DBD, Bastiaansen JWM, Broekhuijse MLWJ, Lopes MS, Knol EF, Harlizius B, et al. Weighted single-step GWAS and gene network analysis reveal new candidate genes for semen traits in pigs. Genet Sel Evol. 2018;50:40. https://doi.org/10.1186/s12711-018-0412-z.
Lourenco DAL, Fragomeni BO, Bradford HL, Menezes IR, Ferraz JBS, Aguilar I, et al. Implications of SNP weighting on single-step genomic predictions for different reference population sizes. J Anim Breed Genet. 2017;134(6):463–71. https://doi.org/10.1111/jbg.12288.
Cáceres P, Barría A, Christensen KA, Bassini LN, Correa K, Garcia B, et al. Genome-scale comparative analysis for host resistance against sea lice between Atlantic Salmon and Rainbow trout. Sci Rep-UK. 2021;11:13231. https://doi.org/10.1038/s41598-021-92425-3.
Dikmen S, Cole JB, Null DJ, Hansen PJ. Genome-wide association mapping for identification of quantitative trait loci for rectal temperature during heat stress in Holstein cattle. PLoS One. 2013;8(7):e69202. https://doi.org/10.1371/journal.pone.0069202.
Srikanth K, Kwon A, Lee E, Chung H. Characterization of genes and pathways that respond to heat stress in Holstein calves through transcriptome analysis. Cell Stress Chaperones. 2017;22(1):29–42. https://doi.org/10.1007/s12192-016-0739-8.
Sengar GS, Deb R, Singh U, Junghare V, Hazra S, Raja TV, et al. Identification of differentially expressed microRNAs in Sahiwal (Bos Indicus) breed of cattle during thermal stress. Cell Stress Chaperones. 2018;23(5):1019–32. https://doi.org/10.1007/s12192-018-0911-4.
Khan A, Dou J, Wang Y, Jiang X, Khan MZ, Luo H, et al. Evaluation of heat stress effects on cellular and transcriptional adaptation of bovine granulosa cells. J Anim Sci Biotechno. 2020;11:25. https://doi.org/10.1186/s40104-019-0408-8.
Li Q, Yang C, Du J, Zhang B, He Y, Hu Q, et al. Characterization of miRNA profiles in the mammary tissue of dairy cattle in response to heat stress. BMC Genomics. 2018;19:975. https://doi.org/10.1186/s12864-018-5298-1.
Wen H, Luo H, Yang M, Augustino SMA, Wang D, Mi S, et al. Genetic parameters and weighted single-step genome-wide association study for supernumerary teats in Holstein cattle. J Dairy Sci. 2021;104(11):11867–77. https://doi.org/10.3168/jds.2020-19943.
Zhao BR, Fu XF, Tian KC, Huang XX, Jiang D, Yan B, et al. Identification of SNPs and expression patterns of FZD3 gene and its effect on wool traits in Chinese merino sheep (Xinjiang type). J Integr Agric. 2019;18(10):2351–60. https://doi.org/10.1016/S2095-3119(19)62735-8.
Dikmen S, Wang XZ, Ortega MS, Cole JB, Null DJ, Hansen PJ. Single nucleotide polymorphisms associated with thermoregulation in lactating dairy cows exposed to heat stress. J Anim Breed Genet. 2015;132(6):409–19. https://doi.org/10.1111/jbg.12176.
Zhao B, Luo H, He J, Huang X, Chen S, Fu X, et al. Comprehensive transcriptome and Methylome analysis delineates the biological basis of hair follicle development and wool-related traits in merino sheep. BMC Biol. 2021;19:197. https://doi.org/10.1186/s12915-021-01127-9.
Elgendy M, Sheridan C, Brumatti G, Martin SJ. Oncogenic ras-induced expression of Noxa and Beclin-1 promotes autophagic cell death and limits clonogenic survival. Mol Cell. 2011;42(1):23–35. https://doi.org/10.1016/j.molcel.2011.02.009.
Lowman XH, McDonnell MA, Kosloske A, Odumade OA, Jenness C, Karim CB, et al. The proapoptotic function of Noxa in human leukemia cells is regulated by the kinase Cdk5 and by glucose. Mol Cell. 2010;40(5):823–33. https://doi.org/10.1016/j.molcel.2010.11.035.
Janus P, Toma-Jonik A, Vydra N, Mrowiec K, Korfanty J, Chadalski M, et al. Pro-death signaling of cytoprotective heat shock factor 1: upregulation of NOXA leading to apoptosis in heat-sensitive cells. Cell Death Differ. 2020;27(7):2280–92. https://doi.org/10.1038/s41418-020-0501-8.
Wang P, Guo J, Wang F, Shi T, Ma D. Human SBK1 is dysregulated in multiple cancers and promotes survival of ovary cancer SK-OV-3 cells. Mol Biol Rep. 2011;38(5):3551–9. https://doi.org/10.1007/s11033-010-0465-8.
Sakabe I, Hu R, Jin L, Clarke R, Kasid UN. TMEM33: a new stress-inducible endoplasmic reticulum transmembrane protein and modulator of the unfolded protein response signaling. Breast Cancer Res TR. 2015;153(2):285–97. https://doi.org/10.1007/s10549-015-3536-7.
Madsen ML, Nettleton D, Thacker EL, Edwards R, Minion FC. Transcriptional profiling of mycoplasma hyopneumoniae during heat shock using microarrays. Infect Immun. 2006;74(1):160–6. https://doi.org/10.1128/IAI.74.1.160-166.2006.
Friederich MW, Timal S, Powell CA, Dallabona C, Kurolap A, Palacios-Zambrano S, et al. Pathogenic variants in glutamyl-tRNAGln Amidotransferase subunits cause a lethal mitochondrial cardiomyopathy disorder. Nat Commun. 2018;9:4065. https://doi.org/10.1038/s41467-018-06250-w.
Liang W, Bi Y, Wang H, Dong S, Li K, Li J. Gene expression profiling of clostridium botulinum under heat shock stress. Biomed Res Int. 2013;2013:760904. https://doi.org/10.1155/2013/760904.
Dou J, Khan A, Khan MZ, Mi S, Wang Y, Yu Y, et al. Heat stress impairs the physiological responses and regulates genes coding for extracellular exosomal proteins in rat. Genes-Basel. 2020;11(3):306. https://doi.org/10.3390/genes11030306.
Signorelli M, Ebrahimpoor M, Veth O, Hettne K, Verwey N, García‐Rodríguez R, et al. Peripheral blood transcriptome profiling enables monitoring disease progression in dystrophic mice and patients. EMBO Mol Med. 2021;13(4):e13328. https://doi.org/10.15252/emmm.202013328.
Toni M, Angiulli E, Miccoli G, Cioni C, Alleva E, Frabetti F, et al. Environmental temperature variation affects brain protein expression and cognitive abilities in adult zebrafish (danio Rerio): a proteomic and behavioural study. J Proteomics. 2019;204:103396. https://doi.org/10.1016/j.jprot.2019.103396.
Rouillard AD, Gundersen GW, Fernandez NF, Wang Z, Monteiro CD, McDermott MG, et al. The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins. Database. 2016;2016:baw100. https://doi.org/10.1093/database/baw100.
Sakai T. Development and regeneration of salivary gland toward for clinical application. Oral Sci Int. 2016;13(1):7–14. https://doi.org/10.1016/S1348-8643(15)00040-3.
Tao YM, Huang JL, Zeng S, Zhang S, Fan XG, Wang ZM, et al. BTB/POZ domain-containing protein 7: epithelial-mesenchymal transition promoter and prognostic biomarker of hepatocellular carcinoma. Hepatology (Baltimore, Md). 2013;57(6):2326–37. https://doi.org/10.1002/hep.26268.
Abbas Z, Sammad A, Hu L, Fang H, Xu Q, Wang Y. Glucose metabolism and dynamics of facilitative glucose transporters (GLUTs) under the influence of heat stress in dairy cattle. Metabolites. 2020;10(8):312. https://doi.org/10.3390/metabo10080312.
Collier RJ, Baumgard LH, Zimbelman RB, Xiao Y. Heat stress: physiology of acclimation and adaptation. Anim Front. 2019;9:12–9.