Transcriptome-wide association study identifies candidate genomic regions for lactation persistency in Murrah buffaloes

PANKAJ  CHAYAL; VIKAS  VOHRA; SUPRIYA  CHHOTARAY; AKSHATA  PATIL; HEMLATA  VALMIKI; DAMINI  SHARMA

doi:10.56093/ijans.v96i5.178691

Authors

PANKAJ CHAYAL ICAR-National Dairy Research Institute, Karnal 132 001 Haryana India
VIKAS VOHRA ICAR-National Dairy Research Institute, Karnal 132 001 Haryana India
SUPRIYA CHHOTARAY ICAR- Central Institute for Research on Buffaloes, Hisar-125 001, Haryana, India
AKSHATA PATIL ICAR-National Dairy Research Institute, Karnal 132 001 Haryana India
HEMLATA VALMIKI ICAR-National Dairy Research Institute, Karnal 132 001 Haryana India
DAMINI SHARMA ICAR- Central Institute for Research on Buffaloes, Hisar-125 001, Haryana, India

DOI:

https://doi.org/10.56093/ijans.v96i5.178691

Keywords:

Lactation persistency, Murrah, Post-peak decline rate, TWAS

Abstract

Lactation persistency is an important economic trait of dairy animals. Understanding the genetic architecture underlying lactation persistency remains a challenge, as conventional genome-wide association studies (GWAS) often lack functional interpretability. The present study aimed to identify genes associated with lactation persistency in Murrah buffaloes using a transcriptome-wide association study (TWAS) approach. Test-day milk yield records from 120 Murrah buffaloes up to 245 days post-calving were used to model lactation curves using Wood’s incomplete gamma function. The post-peak decline rate was considered as a measure of lactation persistency to perform TWAS. TWAS was conducted using an already predicted gene expression panel comprising 26,956 genes. No gene-level association surpassed the Bonferroni-corrected genome-wide significance threshold or the suggestive significance threshold. Although the limited sample size reduced power for detection and no statistically significant association was observed, functional annotation of the top-ranked genes based on their lowest p-values suggested their involvement in biological processes related to mammary epithelial cell survival, proliferation, and differentiation, which are critical for maintaining milk production after peak lactation. The present exploratory analysis highlighted potential genes and pathways that may contribute to regulation of lactation persistency in buffaloes. Further validations could substantiate these genes as candidates for genetic selection in Murrah buffaloes.

Downloads

Download data is not yet available.

References

Abdalla M and Abdalla M. 2022. A general framework for predicting the transcriptomic consequences of non-coding variation and small molecules. PLOS Computational Biology 18(4): p.e1010028. DOI: https://doi.org/10.1371/journal.pcbi.1010028

Buckley CE and St Johnston D. 2022. Apical–basal polarity and the control of epithelial form and function. Nature Reviews Molecular Cell Biology 23(8): 559–577. DOI: https://doi.org/10.1038/s41580-022-00465-y

Cano-Gamez E and Trynka G. 2020. From GWAS to function: using functional genomics to identify the mechanisms underlying complex diseases. Frontiers in Genetics 11: 424. DOI: https://doi.org/10.3389/fgene.2020.00424

Cao C, Ding B, Li Q, Kwok D, Wu J and Long Q. 2021. Power analysis of transcriptome-wide association study: Implications for practical protocol choice. PLOS Genetics 17(2): p.e1009405. DOI: https://doi.org/10.1371/journal.pgen.1009405

Capuco AV, Ellis SE, Hale SA, Long E, Erdman RA, Zhao X and Paape MJ. 2003. Lactation persistency: insights from mammary cell proliferation studies. Journal of Animal Science 81(15): 18–31. DOI: https://doi.org/10.2527/2003.81suppl_318x

Chhotaray S, Vohra V, Uttam V, Santhosh A, Saxena P, Gahlyan RK and Gowane G. 2023. TWAS revealed significant causal loci for milk production and its composition in Murrah buffaloes. Scientific Reports 13(1): 22401. DOI: https://doi.org/10.1038/s41598-023-49767-x

Cole JB, Null DJ and VanRaden PM. 2009. Best prediction of yields for long lactations. Journal of Dairy Science 92(4): 1796–1810. DOI: https://doi.org/10.3168/jds.2007-0976

Deng TX, Ma XY, Duan A, Lu XR and Abdel-Shafy H. 2024. Genome-wide copy number variant analysis reveals candidate genes associated with milk production traits in water buffalo (Bubalus bubalis). Journal of Dairy Science 107(9): 7022– 7037. DOI: https://doi.org/10.3168/jds.2023-24614

Deng T, Liang A, Liang S, Ma X, Lu X, Duan A, Pang C, Hua G, Liu S, Campanile G and Salzano A. 2019. Integrative analysis of transcriptome and GWAS data to identify the hub genes associated with milk yield trait in buffalo. Frontiers in Genetics 10: 36. DOI: https://doi.org/10.3389/fgene.2019.00036

Do DN, Bissonnette N, Lacasse P, Miglior F, Sargolzaei M, Zhao X and Ibeagha-Awemu EM. 2017. Genome-wide association analysis and pathways enrichment for lactation persistency in Canadian Holstein cattle. Journal of Dairy Science 100(3): 1955–1970. DOI: https://doi.org/10.3168/jds.2016-11910

Farhadian M, Rafat SA, Panahi B and Mayack C. 2021. Weighted gene co-expression network analysis identifies modules and functionally enriched pathways in the lactation process. Scientific Reports 11(1): 2367. DOI: https://doi.org/10.1038/s41598-021-81888-z

Fata JE, Mori H, Ewald AJ, Zhang H, Yao E, Werb Z and Bissell MJ. 2007. The MAPKERK-1, 2 pathway integrates distinct and antagonistic signals from TGFα and FGF7 in morphogenesis of mouse mammary epithelium. Developmental Biology 306(1): 193–207. DOI: https://doi.org/10.1016/j.ydbio.2007.03.013

Hu Z, Lu Y, Cai J, Liu J and Wang D. 2024. Cross-oxygen gradients transcriptomic comparison revealed the central role of MAPK and hippo in hypoxia-mediated mammary proliferation inhibition. Antioxidants 13(3): 288. DOI: https://doi.org/10.3390/antiox13030288

Huang DW, Sherman BT and Lempicki RA. 2009. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature Protocols 4(1): 44–57. DOI: https://doi.org/10.1038/nprot.2008.211

Itoh M, Nelson CM, Myers CA and Bissell MJ. 2007. Rap1 integrates tissue polarity, lumen formation, and tumorigenic potential in human breast epithelial cells. Cancer Research 67(10): 4759–4766. DOI: https://doi.org/10.1158/0008-5472.CAN-06-4246

Laodim T, Koonawootrittriron S, Elzo MA, Suwanasopee T, Jattawa D and Sarakul M. 2023. Genetic factors influencing milk and fat yields in tropically adapted dairy cattle: insights from quantitative trait loci analysis and gene associations. Animal Bioscience 37(4): 576. DOI: https://doi.org/10.5713/ab.23.0246

Li S, Qin C, Chen Y, Wei D, Tan Z and Meng J. 2021. Implications of cell division cycle associated 4 on the Wilm’s tumor cells viability via AKT/mTOR signaling pathway. Renal Failure 43(1): 1470–1478. DOI: https://doi.org/10.1080/0886022X.2021.1994994

Lund LR, Bjørn SF, Sternlicht MD, Nielsen BS, Solberg H, Usher PA, Østerby R, Christensen IJ, Stephens RW, Bugge TH and Danø K. 2000. Lactational competence and involution of the mouse mammary gland require plasminogen. Development 127(20): 4481–4492. DOI: https://doi.org/10.1242/dev.127.20.4481

Lund LR, Rømer J, Thomasset N, Solberg H, Pyke C, Bissell MJ, Danø K and Werb Z. 1996. Two distinct phases of apoptosis in mammary gland involution: proteinase-independent and- dependent pathways. Development 122(1): 181–193. DOI: https://doi.org/10.1242/dev.122.1.181

Macciotta NPP, Vicario D and Cappio-Borlino A. 2005. Detection of different shapes of lactation curve for milk yield in dairy cattle by empirical mathematical models. Journal of Dairy Science 88(3): 1178-1191. DOI: https://doi.org/10.3168/jds.S0022-0302(05)72784-3

Meurette O, Stylianou S, Rock R, Collu GM, Gilmore AP and Brennan K. 2009. Notch activation induces Akt signaling via an autocrine loop to prevent apoptosis in breast epithelial cells. Cancer Research 69(12): 5015–5022. DOI: https://doi.org/10.1158/0008-5472.CAN-08-3478

Mu T, Hu H, Ma Y, Wen H, Yang C, Feng X, Wen W, Zhang J and Gu Y. 2022. Identifying key genes in milk fat metabolism by weighted gene co-expression network analysis. Scientific Reports 12(1): 6836. DOI: https://doi.org/10.1038/s41598-022-10435-1

Muir BL, Fatehi J and Schaeffer LR. 2004. Genetic relationships between persistency and reproductive performance in first- lactation Canadian Holsteins. Journal of Dairy Science 87(9): 3029–3037. DOI: https://doi.org/10.3168/jds.S0022-0302(04)73435-9

Das PK, Sejian V, Mukherjee J and Banerjee D. eds. 2023. Textbook of veterinary physiology: 513–568. Singapore: Springer. DOI: https://doi.org/10.1007/978-981-19-9410-4_22

Nagpal S, Meng X, Epstein MP, Tsoi LC, Patrick M, Gibson G, De Jager PL, Bennett DA, Wingo AP, Wingo TS and Yang J. 2019. TIGAR: an improved Bayesian tool for transcriptomic data imputation enhances gene mapping of complex traits. The American Journal of Human Genetics 105(2): 258–266. DOI: https://doi.org/10.1016/j.ajhg.2019.05.018

Nayeri S, Sargolzaei M, Abo-Ismail MK, Miller S, Schenkel F, Moore SS and Stothard P. 2017. Genome-wide association study for lactation persistency, female fertility, longevity, and lifetime profit index traits in Holstein dairy cattle. Journal of Dairy Science 100(2): 1246–1258. DOI: https://doi.org/10.3168/jds.2016-11770

Pedrosa VB, Schenkel FS, Chen SY, Oliveira HR, Casey TM, Melka MG and Brito LF. 2021. Genomewide association analyses of lactation persistency and milk production traits in Holstein cattle based on imputed whole-genome sequence data. Genes 12(11): 1830. DOI: https://doi.org/10.3390/genes12111830

Politi K, Feirt N and Kitajewski J. 2004. Notch in mammary gland development and breast cancer. In Seminars in Cancer Biology 14(5): 341–347. DOI: https://doi.org/10.1016/j.semcancer.2004.04.013

Sehested J, Gaillard C, Lehmann JO, Maciel GM, Vestergaard M, Weisbjerg MR, Mogensen L, Larsen LB, Poulsen NA and Kristensen T. 2019. Extended lactation in dairy cattle. Animal 13(S1): 65–74. DOI: https://doi.org/10.1017/S1751731119000806

Stefanon B, Colitti M, Gabai G, Knight CH and Wilde CJ. 2002. Mammary apoptosis and lactation persistency in dairy animals. Journal of Dairy Research 69(1): 37–52. DOI: https://doi.org/10.1017/S0022029901005246

Turner SD. 2014. qqman: an R package for visualizing GWAS results using QQ and manhattan plots. Biorxiv: 005165. DOI: https://doi.org/10.1101/005165

Visscher PM, Wray NR, Zhang Q, Sklar P, McCarthy MI, Brown MA and Yang J. 2017. 10 years of GWAS discovery: biology, function, and translation. The American Journal of Human Genetics 101(1): 5–22. DOI: https://doi.org/10.1016/j.ajhg.2017.06.005

Vohra V, Chhotaray S, Gowane G, Alex R, Mukherjee A, Verma A and Deb SM. 2021. Genome-wide association studies in Indian Buffalo revealed genomic regions for lactation and fertility. Frontiers in Genetics 12: p.696109. DOI: https://doi.org/10.3389/fgene.2021.696109

Wilson JM, Prokop JW, Lorimer E, Ntantie E and Williams CL. 2016. Differences in the phosphorylation-dependent regulation of prenylation of Rap1A and Rap1B. Journal of Molecular Biology 428(24): 4929–4945. DOI: https://doi.org/10.1016/j.jmb.2016.10.016

Wood PDP. 1969. Factors affecting the shape of the lactation curve in cattle. Animal Science 11(3): 307–316. DOI: https://doi.org/10.1017/S0003356100026945

Xu W, Wang Y, Zhou J, Zhu X, Zhang S, Yuan W, Liu X, Shi Y, Cao L, Zeng Q and Jiang Z. 2016. Cardiac Specific Overexpression of hHole Attenuates Isoproterenol–Induced Hypertrophic Remodeling through Inhibition of Extracellular Signal-Regulated Kinases (ERKs) Signalling. Current Molecular Medicine 16(5): 515–523. DOI: https://doi.org/10.2174/1566524016666160523143704

Yalcin-Ozuysal Ö, Fiche M, Guitierrez M, Wagner KU, Raffoul W and Brisken C. 2010. Antagonistic roles of Notch and p63 in controlling mammary epithelial cell fates. Cell Death and Differentiation 17(10): 1600–1612. DOI: https://doi.org/10.1038/cdd.2010.37

Yang B, Cai Y, Zhu P, Jiang Z, Ao J, Zhang Q, Yuan W, Peng Z, Chen J, Wen Y and Chen Y. 2022. Transmembrane protein 121 as a novel inhibitor of cervical cancer metastasis. Experimental and Therapeutic Medicine 24(3): 572. DOI: https://doi.org/10.3892/etm.2022.11509

Zeng Y, Gao T, Huang W, Yang Y, Qiu R, Hou Y, Yu W, Leng S, Feng D, Liu W and Teng X. 2019. MicroRNA-455-3p mediates GATA3 tumor suppression in mammary epithelial cells by inhibiting TGF-β signaling. Journal of Biological Chemistry 294(43): 15808–15825. DOI: https://doi.org/10.1074/jbc.RA119.010800