The SNP of WBP1 is associated with heifer reproductive performance in the Korean native cattle Hanwoo

ANIMAL
Jiyeon Jeong1Seung-Hwan Lee1Inchul Choi1*

Abstract

It is well documented that intensive selection in dairy cattle for economic value such as increased milk yield led to a decline in reproductive performance. Recent studies using genome-wide association studies (GWASs) discovered candidate genes involved in the lower fertility including embryo development and conception rates. However, the information, which showed a lower reproductive performance, is limited to dairy cattle, especially Holstein, and the candidate genes were not examined in the Korean native cattle Hanwoo which has been intensively selected and bred for meat in the last few decades. We selected the candidate genes WBP1 and PARM1 reported to be associated with cow and/or heifer conception in dairy cattle and analyzed the genotype because those genes have non-synonymous single nucleotide polymorphisms (SNPs). To determine the single base change, we used the high resolution melting (HRM) assay which is rapid and cost-effective for a small number of genes. We found that most heifers with higher conception (1: service per conception) have the AA genotype coding Threonine rather than Proline in the WBP1 gene. We did not detect an association for a SNP in PARM1 in our analysis. In conclusion, the genetic variation of WBP1 can be used as a selective marker gene to improve reproductive performance, and HRM assay can be used to identify common SNP genotypes rapidly and cost effectively.

Keyword



Introduction

It is well documented that intensive selection for increased milk yield led to decline in the reproductive efficiency of dairy cattle, and genome-wide association study (GWAS) identified single nucleotide polymorphisms (SNPs) in candidate genes associated with reproductive performance including bull fertility, fertilization, preimplantation embryo development (Penagaricano and Khatib, 2012; Penagaricano et al., 2012; Cochran et al., 2013a; Cochran et al., 2013b; Penagaricano et al., 2013). Particularly, Cochran et al. (2013a) reported SNPs associated with daughter pregnancy rate, heifer conception rate, cow conception rate, productive life, milk yield in Holstein cattle. However, GWAS studies of reproduction traits were mainly performed in dairy cattle because of economic impact of reproductive performance and breeding strategies ( De Vries, 2006; Galvão et al., 2013; Azizian et al., 2016) rather than in beef cattle. Recently, we reported that Korean native cattle, Hanwoo showed lower reproductive performance, for example calving interval was increased 4.3% and pregnancy rate was decreased about up-to 2.8% year-on-year (Cho et al., 2016; Choi and Cho, 2016).

Here, we postulated that the decline in the reproductive performance in Hanwoo is due to intensive breeding programs using artificial insemination (AI) for economic purpose (Lim et al., 2016). In this context, we carried out in silico studies to select candidate genes associated with heifer pregnancy rate in Korean native cattle based on previous findings in dairy cattle (Cochran et al., 2013a; Ortega et al., 2017) and selected genes (WBP1 and PARM1) containing non-synonymous (missense) SNPs because the amino acid substitutions may have more functional impact on reproductive trait. In this study, we employed a new approach based on high resolution melting (HRM) that is a rapid, reliable and robust method (Garritano et al., 2009; Kysel’ová et al., 2012).

Materials and Methods

Genomic DNA isolation and HRM reaction

To develop the HRM method, we analysed a total of 29 Korean native cattle, Hanwoo from a university farm (Chungyang, Korea). DNA was extracted from whole blood samples using standard protocols for DNA Blood & Tissue kit (Qiagen, Germany). High-resolution melting analysis primer WBP 1 F (5′- ccagcctatgaggatgtggt-3′), WBP1 R (5′- acgtcttccacgtttgttcc-3′), PARM 1 F (5′- cacactcacctcccctcaag -3′), PARM1 R (5′- agatggaacctcaggtggtg -3′) were designed from the bovine genomic sequence and dbSNP (WBP1 rs134282828; PARM1 rs111027720) using the program Beacon Designer (PREMIER Biosoft International, CA, USA). Reactions were run in a volume of 25 μL using Qiagen Type-it® HRM™ PCR mix (Qiagen, Germany). Each reaction contained final concentrations of 1x HRM Type-it™ mix (Qiagen, UK) and 0.7 μM forward and reverse primer (for each target). Each sample was tested in duplicate. Reactions were run on a Qiagen Rotorgene-Q 5-plex machine (Qiagen, Germany) with temperature cycling parameters for the amplification stage being on hold at 95℃ for 5 minutes, followed by 40 cycles of 95℃ for 10 seconds, 55℃ for 30 seconds, and 72℃ for 20 seconds. For the HRM stage, fluorescence recordings were made over the range of 65 - 95℃ by increments of 0.1℃ for 2 seconds. A normalised graph was generated using normalisation regions of 73 - 75℃ and 86 - 88℃.

Results and Discussion

We carried out HRM assay for genotype of WBP1 and PARM1. SNP of WBP is missense variant in which ACT is changed into CCT, consequently amino acid residue Threonine is changed into Proline. Interestingly, in the group of heifers that are pregnant after one AI service, a high proportion of heifers have AA while lower conception or failed heifers have hetero GA or GG (Table 1). However, we did not find relationship between PARM1 and conception rates of heifers. A recent study demonstrated depletion of WBP1 decreased the blastocyst development with a decrease number of trophectoderm cells in the bovine blastocyst (Ortega et al., 2017), suggesting that changes of protein folding pattern by an amino acid substitution (Thr-Pro) may affect protein-protein interaction during embryo development.

Table 1. Genotype of WBP1 and PARM1 using high resolution melting (HRM)

http://dam.zipot.com:8080/sites/kjoas/images/N0030460103_image/Table_KJOAS_46_01_03_T1.jpg

S/C, service per conception.

Conclusion

In line with a loss of function study of WBP1 in dairy cattle, our HRM analysis demonstrated that WBP1 is closely related with conception and maintenance of pregnancy in Korean native beef cattle as well. In summary, our results indicated that genetic variation of WBP1 can be used as a selective marker gene for cattle breeding programs to improve reproductive performance and minimize cost of AI. In addition, HRM assay is able to clearly distinguish a single base variance and can be used for common SNP genotypes rapidly and cost effectively when the targets are selected.

Authors Information

Inchul Choi, https://orcid.org/0000-0001-5011-2658

Jiyeon Jeong, https://orcid.org/0000-0003-0755-0497

Acknowledgements

This work was supported by research fund of Chungnam National University (CNU-2016 to I Choi).

References

1  Azizian S, Shadparvar AA, Ghavi Hossein-Zadeh N, Joezy-Shekalgorabi S. 2016. Effect of increasing accuracy of genomic evaluations on economic efficiency of dairy cattle breeding programmes. Italian Journal of Animal Science 15:379-385.  

2  Cho J, Do C, Choi I. 2016. Reproductive performance of Korean native cattle (Hanwoo) focusing on calving interval and parity. Journal of Embryo Transfer 31:273-279. [In Korean]  

3  Choi I, Cho J. 2016. Reproduction and marketing plans for improving profitability of Korean native cattle (Hanwoo) farm. Journal of Embryo Transfer 31:267-272. [In Korean]  

4  Cochran SD, Cole JB, Null DJ, Hansen PJ. 2013a. Discovery of single nucleotide polymorphisms in candidate genes associated with fertility and production traits in Holstein cattle. Bio Med Central genetics 14:49.  

5  Cochran SD, Cole JB, Null DJ, Hansen PJ. 2013b. Single nucleotide polymorphisms in candidate genes associated with fertilizing ability of sperm and subsequent embryonic development in cattle. Biology of reproduction 89:69.  

6  De Vries A. 2006. Economic value of pregnancy in dairy cattle1. Journal of Dairy Science 89:3876-3885.  

7  Galvão KN, Federico P, De Vries A, Schuenemann GM. 2013. Economic comparison of reproductive programs for dairy herds using estrus detection, timed artificial insemination, or a combination. Journal of Dairy Science 96:2681-2693.  

8  Garritano S, Gemignani F, Voegele C, Nguyen-Dumont T, Le Calvez-Kelm F, De Silva D, Lesueur F, Landi S, Tavtigian SV. 2009. Determining the effectiveness of high resolution melting analysis for SNP genotyping and mutation scanning at the TP53 locus. Bio Med Central Genet 10:5.  

9  Kysel ’ ová J, Rychtáová J, Sztankóová Z, Czerneková V. 2012. Simultaneous identification of CSN3 and LGB genotypes in cattle by high-resolution melting curve analysis. Livestock Science 145:275-279.  

10  Lim D, Strucken EM, Choi BH, Chai HH, Cho YM, Jang GW, Kim T-H, Gondro C, Lee SH. 2016. Genomic footprints in selected and unselected beef cattle breeds in Korea. PLOS ONE 11:e0151324.  

11  Ortega MS, Kurian JJ, McKenna R, Hansen PJ. 2017. Characteristics of candidate genes associated with embryonic development in the cow: Evidence for a role for WBP1 in development to the blastocyst stage. PLOS ONE 12:e0178041.  

12  Penagaricano F, Khatib H. 2012. Association of milk protein genes with fertilization rate and early embryonic development in Holstein dairy cattle. The Journal of dairy research 79:47-52.  

13  Penagaricano F, Weigel KA, Khatib H. 2012. Genome-wide association study identifies candidate markers for bull fertility in Holstein dairy cattle. Animal genetics 43 Suppl 1:65-71.  

14  Penagaricano F, Weigel KA, Rosa GJ, Khatib H. 2013. Inferring quantitative trait pathways associated with bull fertility from a genome-wide association study. Frontiers in genetics 3:307.