Candidate genes for productivity identified by genome-wide association study with indicators of class in the Russian meat merino sheep breed

Genome-wide association studies allow identification of loci and polymorphisms associated with the formation of relevant phenotypes. When conducting a full genome analysis of sheep, particularly promising is the study of individuals with outstanding productivity indicators – exhibition animals, representatives of the super-elite class. The aim of this study was to identify new candidate genes for economically valuable traits based on the search for single nucleotide polymorphisms (SNPs) associated with belonging to different evaluation classes in rams of the Russian meat merino breed. Animal genotyping was performed using Ovine Infinium HD BeadChip 600K DNA, association search was performed using PLINK v. 1.07 software. Highly reliable associations were found between animals belonging to different evaluation classes and the frequency of occurrence of individual SNPs on chromosomes 2, 6, 10, 13, and 20. Most of the substitutions with high association reliability are concentrated on chromosome 10 in the region 10: 30859297–31873769. To search for candidate genes, 15 polymorphisms with the highest association reliability were selected (–log₁₀(р) > 9). Determining the location of the analyzed SNPs relative to the latest annotation Oar_rambouillet_v1.0 allowed to identify 11 candidate genes presumably associated with the formation of a complex of phenotypic traits of animals in the exhibition group: RXFP2, ALOX5AP, MEDAG, OPN5, PRDM5, PTPRT, TRNAS-GGA, EEF1A1, FRY, ZBTB21-like, and B3GLCT-like. The listed genes encode proteins involved in the control of the cell cycle and DNA replication, regulation of cell proliferation and apoptosis, lipid and carbohydrate metabolism, the development of the inflammatory process and the work of circadian rhythms. Thus, the candidate genes under consideration can influence the formation of exterior features and productive qualities of sheep. However, further research is needed to confirm the influence of genes and determine the exact mechanisms for implementing this influence on the phenotype.

1. Abdoli R., Mirhoseini S.Z., Ghavi Hossein-Zadeh N., Zamani P., Moradi M.H., Ferdosi M.H., Gondro C. Genome-wide association study of first lambing age and lambing interval in sheep. Small Rumin. Res. 2019;178:43-45. DOI 10.1016/j.smallrumres.2019.07.014.

2. Amerhanov H.A., Egorov M.V., Selionova M.I., Shumaenko S.N., Efimova N.I. A new breed of sheep: Russian meat Merino. Sel’skokhozyajstvennyj Zhurnal = Agricultural Journal. 2018; 1(11):50-56. DOI 10.25930/0372-3054-2018-1-11-65-71. (in Russian)

3. Bermingham M.L., Bishop S.C., Woolliams J.A., Pong-Wong R., Allen A.R., McBride S.H., Ryder J.J., Wright D.M., Skuce R.A., McDowell S.W., Glass E.J. Genome-wide association study identifies novel loci associated with resistance to bovine tuberculosis. Heredity. 2014;112(5):543-551. DOI 10.1038/hdy.2013.137.

4. Boudon S., Henry-Berger J., Cassar-Malek I. Aggregation of omic data and secretome prediction enable the discovery of candidate plasma biomarkers for beef tenderness. In Int. J. Mol. Sci. 2020; 21(2):54-63. DOI 10.3390/ijms21020664.

5. Buhr E.D., Vemaraju S., Diaz N., Lang R.A., Van Gelder R.N. Neuropsin (OPN5) Mediates Local Light-Dependent Induction of Circadian Clock Genes and Circadian Photoentrainment in Exposed Murine Skin. Curr. Biol. 2019;29(20):3478-3487. DOI 10.1016/j.cub.2019.08.063.

6. Burkitt Wright E.M.M., Spencer H.L., Daly S.B., Manson F.D.C., Zeef L.A.H., Urquhart J., Zoppi N., Bonshek R., Tosounidis I., Mohan M., Madden C., Dodds A., Chandler K. E., Banka S., Au L., Clayton-Smith J., Khan N., Biesecker L.G., Wilson M., Black G.C.M. Mutations in PRDM5 in Brittle Cornea Syndrome Identify a Pathway Regulating Extracellular Matrix Development and Maintenance. Am. J. Hum. Genet. 2011;88(6):767-777. DOI 10.1016/j.ajhg.2011.05.007.

7. Dakhlan A., Moghaddar N., Gondro C., Werf J.H.J. Gene by birth type interaction in merino lamb. Proc. Assoc. Advmt. Anim. Breed. Genet. 2018;22:45-48.

8. Dapas B., Pozzato G., Zorzet S., Capolla S., Paolo M., ScaggianteB., Coan M., Guerra C., Gnan C., Gattei V., Zanconati F., Grassi G. Effects of eEF1A1 targeting by aptamer/siRNA in chronic lymphocytic leukaemia cells. Int. J. Pharm. 2020;57(4):48-59. DOI 10.1016/j.ijpharm.2019.118895.

9. Dominik S., Henshall J.M., Hayes B.J. A single nucleotide polymorphism on chromosome 10 is highly predictive for the polled phenotype in Australian Merino sheep. Anim. Genet. 2012; 43(4):468-470. DOI 10.1111/j.1365-2052.2011.02271.x.

10. Duijvesteijn N., Bolormaa S., Daetwyler H.D., Van Der Werf J.H.J. Genomic prediction of the polled and horned phenotypes in Merino sheep. Genet. Sel. Evol. 2018;50(1):1-11. DOI 10.1186/s12711-018-0398-6.

11. Edea Z., Jeoung Y.H., Shin S.S., Ku J., Seo S., Kim I.H., Kim S.W., Kim K.S. Genome–wide association study of carcass weight in commercial Hanwoo cattle. Asian-Australas. J. Anim. Sci. 2018; 31(3): 327-334. DOI 10.5713/ajas.17.0276.

12. García-Gámez E., Gutiérrez-Gil B., Sahana G., Sánchez J.P., Bayón Y., Arranz J.J. GWA Analysis for Milk Production Traits in Dairy Sheep and Genetic Support for a QTN Influencing Milk Protein Percentage in the LALBA Gene. PLoS ONE. 2012; 7(10):1-9. DOI 10.1371/journal.pone.0047782.

13. Georges M., Charlier C., Hayes B. Harnessing genomic information for livestock improvement. Nat. Rev. Genet. 2019;20(3): 135-156. DOI 10.1038/s41576-018-0082-2.

14. Gudmundsdottir O.O. Genome-wide association study of muscle traits in Icelandic sheep. Agricultural University of Iceland, 2015.

15. Johnston S.E., Gratten J., Berenos C., Pilkington J.G., CluttonBrock T.H., Pemberton J.M., Slate J. Life history trade-offs at a single locus maintain sexually selected genetic variation. Nature. 2013;502(7469):93-95. DOI 10.1038/nature12489.

16. Johnston S.E., McEwan J.C., Pickering N.K., Kijas J.W., Beraldi D., Pilkington J.G., Pemberton J. M., Slate J. Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population. Mol. Ecol. 2011;20(12):2555-2566. DOI 10.1111/j.1365-294X.2011.05076.x.

17. Kamiński S., Hering D.M., Oleński K., Lecewicz M., Kordan W. Genome-wide association study for sperm membrane integrity in frozen-thawed semen of Holstein-Friesian bulls. Anim. Reprod. Sci. 2016;170:135-140. DOI 10.1016/j.anireprosci.2016.05.002.

18. Kijas J.W., Serrano M., Mcculloch R., Li Y., Salces Ortiz J., Calvo J.H., Pérez-Guzmán M.D. Genomewide association for a dominant pigmentation gene in sheep. J. Anim. Breed. Genet. 2013;130(6):468-475. DOI 10.1111/jbg.12048.

19. Lee J.R. Protein tyrosine phosphatase PTPRT as a regulator of synaptic formation and neuronal development. BMB Rep. 2015; 48(5):249-255. DOI 10.5483/BMBRep.2015.48.5.037.

20. Liu Y., Chen X., Gong Z., Zhang H., Fei F., Tang X., Wang J., Xu P., Zarbl H., Ren X. Fry Is Required for Mammary Gland Development During Pregnant Periods and Affects the Morphology and Growth of Breast Cancer Cells. Front. Oncol. 2019;4(2):1-12. DOI 10.3389/fonc.2019.01279.

21. Mastrangelo S., Moioli B., Ahbara A., Latairish S., Portolano B., Pilla F., Ciani E. Genome-wide scan of fat-tail sheep identifies signals of selection for fat deposition and adaptation. Anim. Prod. Sci. 2019;59(5):835-842. DOI 10.1071/AN17753.

22. Moioli B., Pilla F., Ciani E. Signatures of selection identify loci associated with fat tail in sheep. J. Anim. Sci. 2015;93(10):4660- 4669. DOI 10.2527/jas.2015-9389.

23. Purcell S., Neale B., Todd-Brown K., Thomas L., Ferreira M.A.R., Bender D., Maller J., Sklar P., Bakker P.I.W., Daly M.J., Sham P.C. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007; 81(3):559-575. DOI 10.1086/519795.

24. Rossi U.A., Hasenauer F.C., Caffaro M.E., Neumann R., Salatin A., Poli M.A., Rossetti C.A. A haplotype at intron 8 of PTPRT gene is associated with resistance to Brucella infection in Argentinian creole goats. Vet. Microbiol. 2017;2(3):133 137. DOI 10.1016/j.vetmic.2017.06.001.

25. Rovadoscki G.A., Pertile S.F.N., Alvarenga A.B., Cesar A.S.M., Pértille F., Petrini J., Franzo V., Soares W.V.B., Morota G., Spangler M.L., Pinto L.F.B., Carvalho G.G.P., Lanna D.P.D., Coutinho L.L., Mourão G.B. Estimates of genomic heritability and genome-wide association study for fatty acids profile in Santa Inеs sheep. BMC Genom. 2018;19(1):1-14. DOI 10.1186/ s12864-018-4777-8.

26. Scott D.J., Rosengren K.J., Bathgate R.A.D. The different ligandbinding modes of relaxin family peptide receptors RXFP1 and RXFP2. Mol. Endocrinol. 2012;26(11):1896-1906. DOI 10.1210/me.2012-1188.

27. Selionova M.I., Shumaenko S.N., Efimova N.I. Surov A.I. Bobrishov S.S. Target indicators and characteristics of the Russian Meat Merino breed: Proceedings of the Research Institute for Sheep and Goat Farming. Sel’skokhozyajstvennyj Zhurnal = Agricultural Journal. 2017;2(10):10-16. (in Russian)

28. Wang J., Kudoh J., Takayanagi A., Shimizu N. Novel human BTB/ POZ domain-containing zinc finger protein ZNF295 is directly associated with ZFP161. Biochem. Biophys. Res. Commun. 2005;327(2):615-627. DOI 10.1016/j.bbrc.2004.12.048.

29. Wang Z., Zhang H., Yang H., Wang S., Rong E., Pei W., Li H., Wang N. Genome-wide association study for wool production traits in a Chinese merino sheep population. PLoS ONE. 2014; 9(9):3-10. DOI 10.1371/journal.pone.0107101.

30. Weh E., Takeuchi H., Muheisen S., Haltiwanger R.S., Semina E.V. Functional characterization of zebrafish orthologs of the human Beta 3-Glucosyltransferase B3GLCT gene mutated in Peters Plus Syndrome. PLoS ONE. 2017;12(9):e0184903. DOI 10.1371/journal.pone.0184903.

31. Wiedemar N., Drögemüller C.A 1.8-kb insertion in the 3′-UTR of RXFP2 is associated with polledness in sheep. Anim. Genet. 2015;46(4):457-461. DOI 10.1111/age.12309.

32. Xu S.S., Gao L., Xie X.L., Ren Y.L., Shen Z.Q., Wang F., Shen M., Eypórsdóttir E., Hallsson J.H., Kiseleva T., Kantanen J., Li M.H. Genome-wide association analyses highlight the potential for dif¬ferent genetic mechanisms for litter size among sheep breeds. Front. Genet. 2018;6(2):1-14. DOI 10.3389/fgene.2018.00118.

33. Yan J., Dukkipati V., Blair H.T., Biggs P.J., Hamie J.C., Aw G. A genome-wide scan of positive selection signature using Ovine Infinium HD SNP BeadChip in two Romney lines, selected for resistance or resilience to nematodes. Anim. Genet. 2017;4: 87-94.

34. Zhang H., Chen X., Sairam M.R. Novel genes of visceral adiposity: Identification of mouse and human Mesenteric Estrogen-Dependent Adipose (MEDA)-4 gene and its adipogenic function. Endocrinology. 2012;153(6):2665-2676. DOI 10.1210/en.2011- 2008.

35. Zhang T., Gao H., Sahana G., Zan Y., Fan H., Liu J., Shi L., Liu J., Du L., Wang L., Zhao F. Genome-wide association studies revealed candidate genes for tail fat deposition and body size in the Hulun Buir sheep. J. Anim. Breed. Genet. 2019;6(1):1-9. DOI 10.1111/jbg.12402.