doi: 10.15389/agrobiology.2020.2.257eng
UDC: 636.2:636.082:577.21
Supported financially by the Ministry of Science and Higher Education of the Russian Federation (theme No. AAAA-A18-118021590134-3) and the Ministry of Agriculture of the Russian Federation (theme No. АААА-А19-119031590044-3)
GENOMIC VARIABILITY ASSESS FOR BREEDING TRAITS IN HOLSTEINIZATED RUSSIAN BLACK-AND-WHITE CATTLE USING GWAS ANALYSIS AND ROH PATTERNS
A.A. Sermyagin1, O.A. Bykova2, O.G. Loretts2, O.V. Kostyunina1, N.A. Zinovieva1
1Ernst Federal Science Center for Animal Husbandry, 60, pos. Dubrovitsy, Podolsk, Moscow Province, 142132 Russia, e-mail alex_sermyagin85@mail.ru (✉ corresponding author), kostolan@mail.ru, n_zinovieva@mail.ru;
2Ural State Agricultural University, 42, ul. Karl Liebknecht, Еkaterinburg, Sverdlovsk Province, 620075 Russia, e-mail olbyk75@mail.ru, rector.urgau@yandex.ru
ORCID:
Sermyagin A.A. orcid.org/0000-0002-1799-6014
Kostyunina O.V. orcid.org/0000-0001-8206-3221
Bykova O.A. orcid.org/0000-0002-0753-1539
Zinovieva N.A. orcid.org/0000-0003-4017-6863
Loretts O.G. orcid.org/0000-0002-9945-5691
Received January 31, 2020
For using the pedigree information recorded in dairy herds rational, the creation of reference groups with a high reliability the estimated breeding value of animals is required. This is a necessary requirement not only for the genetic assessment procedure, but also for the introduction of genomic selection methods. However, without high-quality phenotyping cannot be this achieved. Therefore, it seems relevant to conduct model studies using the highly productive herd of Holsteinizated Russian Black-and-White cattle in the Urals as an example to study the genomic variability of phenotypic traits in animals’ different generations using data from genome-wide analysis (GWAS) and runs of homozygosity (ROH), as well as to replenish the Russian reference population. The novelty of the work is to assess the dynamics of genomic inbreeding variability (FROH) in animal generations (“dam of mother”—“mother”-“daughter”) and its comparison with direct genomic value (DGV), as well as the search for new mechanisms to confirm the scientific hypothesis of using a limited experimental dataset in GWAS. The object of this research was 76 cows and heifers of Holsteinizated Russian Black-and-White breed, as well as 9 Holstein sires genotyped for 139 thousand SNPs by the Bovine GGP HD platform (Illumina/Neogen, USA). To calculate the DGVs of studied animals the Russian reference bulls and cows group that include 591 individuals was used. GWAS and ROH analyzes were performed based on 110448 SNPs. The reliability (p-values) of genome-wide associations with direct cows’ phenotypes ranged from 2.31×10-5 to 1.08×10-7. Quantitative traits loci on autosomes BTA1, BTA2, BTA5, BTA7, BTA8, BTA10, BTA11, BTA12, BTA14, BTA16, BTA20, BTA21, and BTA26 were found. For milk yield a region on BTA14 (1.44-1.59 Mb) with the genes ZNF16, ARHGAP39, ZNF7 associated with an increased fat milk yield was detected. For the number of inseminations found SNPs included into the genes (ARHGAP31) or located close to the genes (SERPINA5) and associated with the growth intensity to a mature state, as well as ovarian function in animals. The characteristic of ROHs depending on the length of their fragments in the genome is given. Conservative homozygous regions on BTA12, BTA14, BTA26, and BTA29 and the most significant genes entering them were identified, which are potentially associated with selection pressure in the studied population mainly by milk production traits, reproduction, and udder type measurement parameters. The value of FROH significantly (р < 0.05-0.001) increased in the offspring—parent generations: by +0.012 or 1.2 % for mothers, and +0.029 for daughters. The highest values of FROH = 0.135 were noted for bulls that were signed as fathers of cows (generation of “mother”) and heifers (generation of “daughter”). Each subsequent generation of individuals showed an average increase in DGV for milk yield by +94.2 kg, fat milk yield by +4.4 kg and protein milk yield by +3.0 kg, reflecting clear the strategy to improve milk production traits in the herd for obtaining new cattle genotypes. Thus, the possibilities for assessing the variability of direct cattle phenotypes, as a model for studying genomic variability in a single herd, based on the search for associations and loci in the genome under selection pressure are shown.
Keywords: cattle, GWAS, ROH, reference population, genomic inbreeding, genomic evaluation, milk production, fertility.
REFERENCES
- Zinov'eva N., Strekozov N., Yanchukov I., Ermilov A., Eskin G. Zhivotnovodstvo Rossii, 2015, 3: 27-30 (in Russ.).
- Kudinov A.A., Maslennikova E.S., Plemyashov K.V. Zootekhniya, 2019, 1: 2-6 (in Russ.).
- Sermyagin A.A., Gladyr E.A., Plemyashov K.V., Kudinov A.A., Dotsev A.V., Deniskova T.E., Zinovieva N.A. Genome-wide association studies for milk production traits in Russian population of Holstein and Black-and-White cattle. Proceedings of the Scientific-Practical Conference «Research and Development — 2016». K. Anisimov et al. (eds.). Springer, Cham, 2018: 591-599 CrossRef
- Qanbari S., Pimentel E.C.G., Tetens J., Thaller G., Lichtner P., Sharifi A.R., Simianer H. A genome-wide scan for signatures of recent selection in Holstein cattle. Animal Genetics, 2010, 41(4): 377-389 CrossRef
- Wang X., Wurmser C., Pausch H., Jung S., Reinhardt F., Tetens J., Thaller G., Fries R. Identification and dissection of four major QTL affecting milk fat content in the German Holstein-Friesian population. PLoS ONE, 2012, 7(7): e40711 CrossRef
- Maurice-Van Eijndhoven M.H.T., Bovenhuis H., Veerkamp R.F., Calus M.P.L. Overlap in genomic variation associated with milk fat composition in Holstein Friesian and Dutch native dual-purpose breeds. Journal of Dairy Science, 2015, 98(9): 6510-6521 CrossRef
- Chamberlain A.J., Hayes B.J., Savin K., Bolormaa S., McPartlan H.C., Bowman P.J., Van Der Jagt C., MacEachern S., Goddard M.E. Validation of single nucleotide polymorphisms associated with milk production traits in dairy cattle. Journal of Dairy Science, 2012, 95(2): 864-875 CrossRef
- Fang M., Fu W., Jiang D., Zhang Q., Sun D., Ding X., Liu J. A multiple-SNP approach for genome-wide association study of milk production traits in Chinese Holstein cattle. PLoS ONE, 2014, 9(8): e99544 CrossRef
- Xiang R., MacLeod I.M., Bolormaa S., Goddard M.E. Genome-wide comparative analyses of correlated and uncorrelated phenotypes identify major pleiotropic variants in dairy cattle. Scientific Reports, 2017, 7: 9248 CrossRef
- Zhang F., Wang Y., Mukiibi R., Chen L,, Vinsky M., Plastow G., Basarab J., Stothard P., Li G. Genetic architecture of quantitative traits in beef cattle revealed by genome wide association studies of imputed whole genome sequence variants: I: feed efficiency and component traits. BMC Genomics, 2020, 21: 36 CrossRef
- Reyer H., Shirali M., Ponsuksili S., Murani E., Varley P.F., Jensen J., Wimmers K. Exploring the genetics of feed efficiency and feeding behaviour traits in a pig line highly selected for performance characteristics. Molecular Genetics and Genomics: MGG, 2017, 292(5): 1001-1011 CrossRef
- Mel'nikova E.E., Bardukov N.V., Fornara M.S., Kostyunina O.V., Sermyagin A.A., Brem G., Zinov'eva N.A. The study of effect of genotypes for DNA marker on reproductive qualities of sows of Large White and Landrace breeds. Sel'skokhozyaistvennaya biologiya [Agricultural Biology], 2019, 54(2): 227-238 CrossRef
- Belous A.A., Sermyagin A.A., Kostyunina O.V., Brem G., Zinov'eva N.A. Study of genetic architecture of feed conversion rate in duroc young boars (Sus scrofa) based on the genome-wide SNP analysis. Sel'skokhozyaistvennaya biologiya [Agricultural Biology], 2019, 54(4): 705-712 CrossRef
- Pryce J.E., Johnston J., Hayes B.J., Sahana G., Weigel K.A., McParland S., Spurlock D., Krattenmacher N., Spelman R.J., Wall E., Calus M.P.L. Imputation of genotypes from low density (50,000 markers) to high density (700,000 markers) of cows from research herds in Europe, North America, and Australasia using 2 reference populations. Journal of Dairy Science, 2014, 97(3): 1799-1811 CrossRef
- Ferenčaković M., Sölkner J., Curik I. Estimating autozygosity from high-throughput information: effects of SNP density and genotyping errors. Genetics Selection Evolution, 2013, 45: 42 CrossRef
- Ignat'eva L.P., Belous A.A., Nedashkovskii I.S., Kostyunina O.V., Sermyagin A.A., Zinov'eva N.A. Sibirskii vestnik sel'skokhozyaistvennoi nauki, 2019, 49(6): 79-87 CrossRef (in Russ.).
- Bakoev C.Yu., Getmantseva L.V. Dostizheniya nauki i tekhniki APK, 2019, 33(11): 63-68 CrossRef (in Russ.).
- Smaragdov M.G., Kudinov A.A. Dostizheniya nauki i tekhniki APK, 2019, 33(6): 51-53 CrossRef (in Russ.).
- Sermyagin A.A., Zinov'eva N.A. Dostizheniya nauki i tekhniki APK, 2019, 33(12): 77-82 CrossRef (in Russ.).
- Chang C.C., Chow C.C., Tellier L.CAM., Vattikuti S., Purcell S.M., Lee J.J. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience, 2015, 4: 7 CrossRef
- Zhang L., Orloff M.S., Reber S., Li S., Zhao Y., Eng C. aTOH: Extended Approach for Identifying Tracts of Homozygosity. PLoS ONE, 2013, 8(3): e57772 CrossRef
- Hu Z.-L., Park C.A., Reecy J.M. Building a livestock genetic and genomic information knowledgebase through integrative developments of Animal QTLdb and CorrDB. Nucleic Acids Research, 2019, 47: D701-710 CrossRef
- VanRaden P.M. Efficient methods to compute genomic predictions. Journal of Dairy Science, 2008, 91(11): 4414-4423 CrossRef
- Jiang J., Ma L., Prakapenka D., VanRaden P.M., Cole J.B., Da Y. A large-scale genome-wide association study in U.S. Holstein cattle. Front. Genet., 2019, 10: 412 CrossRef
- Sun J., Shan L., Zhang C., Chen H. Haplotype combination of the bovine PCSK1 gene sequence variants and association with growth traits in Jiaxian cattle. Journal of Genetics, 2014, 93(3): e123-e129. Rezhim dostupa: http://www.ias.ac.in/jgenet/OnlineResources/93/e123.pdf. Data obrashcheniya: 27.01.2020 goda.
- Friedrich J., Brand B., Ponsuksili S., Graunke K.L., Langbein J., Knaust J., Kühn C., Schwerin M. Detection of genetic variants affecting cattle behaviour and their impact on milk production: a genome-wide association study. Anim. Genet., 2016, 47(1): 12-18 CrossRef
- Karisa B.K., Thomson J., Wang Z., Stothard P., Moore S.S., Plastow G.S. Candidate genes and single nucleotide polymorphisms associated with variation in residual feed intake in beef cattle. Journal of Animal Science, 2013, 91(8): 3502-3513 CrossRef
- Liu Y., Zan L., Xin Y., Tian W., Li L., Wang H. ZBTB38 gene polymorphism associated with body measurement traits in native Chinese cattle breeds. Gene, 2013, 513(2): 272-277 CrossRef
- Valente T.S., Baldi F., Sant'Anna A.C., Albuquerque L.G., Paranhos da Costa M.J. Genome-wide association study between single nucleotide polymorphisms and flight speed in Nellore cattle. PLoS ONE, 2016, 11(6): e0156956 CrossRef
- Lima A.O.D., Oliveira P.S.N., Tizioto P.C., Somavilla A.L., Diniz W.J.S., Silva J.V.D., Andrade S.C.S., Boschiero C., Cesar A.S.M., Souza M.M., Rocha M.I.P., Afonso J., Buss C.E., Mudadu M.A., Mourao G.B., Coutinho L.L., Regitano L.C.A. PRUNE2 gene has a potential effect on residual feed intake in Nellore cattle. Journal of Animal Science, 2016,94(E-Suppl. 5): 152-153 CrossRef
- Cole J.B., Wiggans G.R., Ma L., Sonstegard T.S., Lawlor T.J., Crooker B.A., Van Tassell C.P., Yang J., Wang S., Matukumalli L.K., Da Y. Genome-wide association analysis of thirty one production, health, reproduction and body conformation traits in contemporary U.S. Holstein cows. BMC Genomics, 2011, 12: 408 CrossRef
- Szewczuk M. Polymorphism in exon 2 encoding the putative ligand-binding pocket of the bovine insulin-like growth factor 1 receptor affects milk traits in four different cattle breeds. Journal of Animal Breeding and Genetics, 2017, 134(1): 34-42 CrossRef
- Leal-Gutiérrez J.D., Elzo M.A., Johnson D.D., Hamblen H., Mateescu R.G. Genome wide association and gene enrichment analysis reveal membrane anchoring and structural proteins associated with meat quality in beef. BMC Genomics, 2019, 20(1): 151 CrossRef
- Hayashi K., Ushizawa K., Hosoe M., Takahashi T. Differential gene expression of serine protease inhibitors in bovine ovarian follicle: possible involvement in follicular growth and atresia. Reproductive Biology and Endocrinology, 2011, 9: 72 CrossRef
- Zhuang Z., Xu L., Yang J., Gao H., Zhang L., Gao X., Li J., Zhu B. Weighted single step genome-wide association study for growth traits in Chinese Simmental beef cattle. Genes, 2020, 11: 189 CrossRef
- Zhang K., Rajput S.K., Wang S., Folger J.K., Knott J.G., Smith G.W. CHD1 regulates deposition of histone variant H3.3 during bovine early embryonic development. Biology of Reproduction, 2016, 94(6): 140, 1-8 CrossRef
- Pryce J.E., Hayes B.J., Bolormaa S., Goddard M.E. Polymorphic regions affecting human height also control stature in cattle. Genetics, 2011, 187(3): 981-984 CrossRef
- Twomey A.J., Berry D.P., Evans R.D., Doherty M.L., Graham D.A., Purfield D.C. Genome-wide association study of endo-parasite phenotypes using imputed whole-genome sequence data in dairy and beef cattle. Genetics Selection Evolution, 2019, 51: 15 CrossRef
- Wang X.G., Ju Z.H., Hou M.H., Jiang Q., Yang C.H., Zhang Y., Sun Y., Li R.L., Wang C.F., Zhong J.F., Huang J.M. Deciphering transcriptome and complex alternative splicing transcripts in mammary gland tissues from cows naturally infected with Staphylococcus aureus mastitis. PLoS ONE, 2016, 11(12): e0167666 CrossRef
- Edea Z., Jeoung J.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-Australasian Journal of Animal Sciences (AJAS), 2018, 31(3): 327-334 CrossRef
- Utsunomiya Y.T., Milanesi M., Utsunomiya A.T.H., Torrecilha R.B.P., Kim E.S., Costa M.S., Aguiar T.S., Schroeder S., do Carmo A.S., Carvalheiro R., Neves H.H.R., Padula R.C.M., Sussai T.S., Zavarez L.B., Cipriano R.S., Caminhas M.M.T., Hambrecht G., Colli L., Eufemi E., Ajmone-Marsan P., Cesana D., Sannazaro M., Buora M., Morgante M., Liu G., Bickhart D., Van Tassell C.P., Sölkner J., Sonstegard T.S., Garcia J.F. A PLAG1 mutation contributed to stature recovery in modern cattle. Scientific Reports, 2017, 7(1): 17140 CrossRef
- Li C., Sun D., Zhang S., Liu L., Alim M.A., Zhang Q. A post-GWAS confirming the SCD gene associated with milk medium- and long-chain unsaturated fatty acids in Chinese Holstein population. Animal Genetics, 2016, 47(4): 483-490 CrossRef
- Hering D.M., Olenski K., Kaminski S. Genome-wide association study for poor sperm motility in Holstein-Friesian bulls. Animal Reproduction Science, 2014, 146(3-4): 89-97 CrossRef
- Santos D.J.A., Cole J.B., Null D.J., Byrem T.M., Ma L. Genetic and nongenetic profiling of milk pregnancy-associated glycoproteins in Holstein cattle. Journal of Dairy Science, 2018, 101(11): 9987-10000 CrossRef
- Ibeagha-Awemu E.M., Akwanji K.A., Beaudoin F., Zhao X. Associations between variants of FADS genes and omega-3 and omega-6 milk fatty acids of Canadian Holstein cows. BMC Genetics, 2014, 15: 25 CrossRef