Identification of Single Nucleotide Polymorphisms in AMBP gene and their associations with economic traits in Korean cattle (Hanwoo)

Review
  Jung-Min Han1Hong Sik Kong1*

Abstract

The alpha-1-microglobulin/Bikunin precursor (AMBP) gene encodes alpha-1-microglobulin (A1M), an immunosuppressive lipocalin, and Bikunin is located on bovine chromosome 8 and single nucleotide polymorphisms (SNPs) in AMBP gene, a gene in the QTL region, are associated with meat quality traits in swine. This study aims to investigate the AMBP gene in association to economic traits in Korean cattle (Hanwoo). Through DNA sequencing, we detected a total of 14 (C119A, T148C, T179C, T241A, G4500T, T4605C, G9457A, T9638C, T9648G, C9662T, G9675A, G9936A, T9939C, G9998A) in different regions. We investigated the effect of these mutations by determining their association and variation with economic traits in Hanwoo. Our statistical association analysis revealed that T148C, T179C, T241A, G9457A, C9662T, G9936A, T9939C and G9998A were significantly associated with carcass weight, longissimus dorsi muscle area, backfat thickness and marbling score. Therefore, the AMBP gene may thus serve as a genetic marker for economic traits and as a basis for molecular breeding in Hanwoo.

Keyword



Introduction

Hanwoo (Bos taurus) is a Korean breed of cattle that has been in used as draft animal for 5,000 years (Kim et al., 2011). Hanwoo may have developed independently from an ancestor of the European and African species (Yoon et al., 2005). In Korea, consumer demand is driving efforts to increase meat production and produce higher quality meat (Chung, 2005; Seong et al., 2011). All beef carcasses in Korea are graded by the Korea Institute for Animal Products Quality Evaluation (KAPE; www.ekape.or.kr) both in qualitative and quantative terms. Under KAPE standards (Kim at al., 2011; Choy et al., 2012), the major quantitative components that determine carcass grades are backfat thickness (BF), longissimus dorsi muscle area (LMA), carcass weight (CW) and marbling score (MS).

Earlier reports indicate that single nucleotide polymorphisms (SNPs) in alpha-1-microglobulin/Bikunin precursor (AMBP) gene, a gene in the QTL region (http://www.animalgenome.org/cgi-bin/QTLdb/BT/ link_genom2snpb?chromos=21&upcm=73.81&locm=76.81),are associated with meat quality traits in swine(Cinaretal.,2012). AMBP gene is located on the bovine chromosome 8 and consists of ten exons. It is a single-copy gene and is expressed exclusively in the liver (Grewal et al., 2005; Tyagi et al., 2002), encoding two plasma glycoproteins: alpha-1-microglobulin (A1M) and Bikunin (Grewal et al., 2005; Tyagi et al., 2002; Allhorn et al., 2002). Although A1M and Bikunin originate from AMBP gene cleavage by a furin-like protease, the structural and functional relationship between the two mature molecules remains unknown (Tyagi et al., 2002). The first 6 exons code for A1M, exons 8 and 9 code for each Kunitz domain and exons 7 and 10 encode the N- and C-terminal ends of Bikunin, respectively (Fries and Blom, 2000).

AMBP gene is reported to assumed to be cis-regulated and was part of the glycolysis pathway (Heidt et al., 2013). A1M is to play a role in the regulation of development (Logdberg et al., 2000), cell growth (Logdberg and Akerstrom, 1981; Logdberg et al., 1986), metabolism (Allhorn et al., 2002) and the immune response (Wester et al., 2000; Santin and Cannas, 1999). Bikunin is a urinary trypsin inhibitor belonging to the Kunitz-type superfamily of protease inhibitors (Fries and Blom, 2000). It is thought to participate in the control of such events as endothelial cell growth or extracellular matrix stabilization (Tyagi et al., 2002).

The aim to the study was to detection of novel polymorphisms within AMBP gene by sequencing and to determine the genetic association between the SNPs and economic traits, such meat quality and carcass traits

Materials and methods

Sample collection and Genomic DNA Extraction

A total of 105 Korean cattle (Hanwoo) steers obtained from Korea Institute for Animal Products Quality Evaluation (KAPE) were used as samples. Genomic DNA from tissue samples was extracted using a Toyobo MagExtractor Kit, as described in the manufacturer's instructions (Toyobo, Japan). The DNA was quantified by NanoDrop ND-1000 UV - Vis Spectrophotometer (NanoDrop Technologies, USA). For each specimen, the weight at slaughter was recorded as the CW. BF and LMA were calculated at the 12th and 13th rib interface. MS was estimated on a cross section at the 12th and 13thrib interface and scored on a scale from 1 to 9 associated with the most marbling.

Primer design and Polymerase Chain Reaction

Primers used in the study were based on the nucleotide sequence of the bovine AMBP gene available in the GenBank database (Accession No. AC_000165). Detailed information on oligonucleotide primers, amplicon size and corresponding annealing temperature is depicted in Table 1.

Table 1. Primer sequences, location and product size for SNP analysis of AMBP gene.

http://dam.zipot.com:8080/sites/jabg/images/N0270020304/Table_JABG_02_03_04_T1.jpg

PCR amplification was performed using genomic DNA (20~200 ng/㎕), 10 pmole of each primer, 2 ㎕ of 10X buffer (10 mM Tris-HCl, 50 mM KCl, 1.5 mM MgCl2,pH8.3),40㎛ of dNTPs and 0.5 units Taq polymerase (Genetbio, Korea). Thermal cycling began with an initial denaturation step of 94℃ for 5 min, followed by 35 cycles of 94℃ for 1 min, optimized annealing temperature 63.5℃ for 30 s, 72℃ for 30 s and a final extension at 72℃ for 5 min using a Peltier Thermal Cycler 200 (MJ Research, USA).

Sequencing analysis

The amplified PCR product of each sample was excised from a 2.0% agarose gel and purified with AccuPrep®Gel PurificationKit(Bioneer™, Korea). The purified products were sequenced using BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, USA). Sequencing reactions were purified by Ethanol and run on ABI 3130 XL Genetic Analyzer (Applied Biosystems, USA) to determine nucleotide sequences. Detection of SNPs were performed with Seq Man II software (DNA Star, USA).

Statistical analysis

Allele and genotype frequencies were calculated by simple counting. A linkage disequilibrium (LD) block of polymorphisms was tested by Haploview version 3.32 (Barrett et al., 2005). The haplotypes and their frequencies were calculated with the expectation-maximization algorithm (Stephens et al., 2001). Statistical analysis of the association between the genotypes (or haplotypes) of A1M and economic traits were evaluated with the least square method (GLM implemented in the SAS software package: SAS institute, USA). The general linear model used in this study was

Yij=μ+Gi+βA+Si+eij,

Where Yijis the observation of the performance traits, μ is the overall mean for each trait, Giis the ith level of the fixed factor genotype, β is the regression coefficient for age, A is the age, Si is the it sex effect and eij is the random error for the ijth individual.

Results

We investigated genetic variations within the bovine AMBP gene using DNA sequencing technologies. By comparing the sequencing results with the DNA sequence of Bos taurus AMBP gene, we identified 14 SNPs. 5 SNPs (C119A, T148C, T179C, G4500T, G9457A) were detected in the exon regions and 9 SNPs (T241A, T4605C, T9638C, T9648G, C9662T, G9675A, G9936A, T9939C, G9998A) were detected in the intron regions (Fig. 1, Table 2). The T148C was missense mutation. T148C changes amino acid from valine to alanine (Table 2). The locations, allele frequency, and genotype frequency for Hanwoo are shown in Table 3. The frequency of T148C and T179C polymorphisms in allele C was higher than that for allele T,and 0.488, 0.476 and 0.036 displayed the CC, CT and TT genotypes, respectively. The frequency of T148C and T179C are the same in genotype. And the genotype frequencies were as follows : 0.194 CC, 0.548 CA and 0.258 AA for the C119A; 0.310 AA, 0.595 AT and 0.095 TT for the T241A; 0.270 CC, 0.568 CT and 0.162 TT for the C4500T; 0.417 TT, 0.555 TC and 0.028 CC for the T4605C; 0.091 GG, 0.352 GA and 0.557AA for the G9457A; 0.110 TT, 0.305 TC and 0.585 CC for the T9638C; 0.582 TT, 0.316 TG and 0.101 GG for the T9648G; 0.812 CC, 0.163 CT and 0.025 TT for the C9662T; 0.570 GG, 0.329 GA and 0.101 AA for the G9675A; 0.110 GG, 0.341 GA and 0.549 AA for the G9936A; 0.549 TT, 0.352 TC and 0.099 CC for the T9939C; 0.553 GG, 0.351 GA and 0.096 AA for the G9998A (Table 3).

http://dam.zipot.com:8080/sites/jabg/images/N0270020304/Figure_JABG_02_03_04_F1.jpg

Figure 1. AMBP Gene structure and detected SNPs in the Hanwoo cattle.

Table 2. Location and amino acid change of 14 polymorphism detected in AMBP gene

http://dam.zipot.com:8080/sites/jabg/images/N0270020304/Table_JABG_02_03_04_T2.jpg

Table 3. Genotype and allele frequencies of SNPs in Hanwoo on AMBP gene.

http://dam.zipot.com:8080/sites/jabg/images/N0270020304/Table_JABG_02_03_04_T3.jpg

Pairwise LD analysis with the 3 polymorphisms parsed into 1LD block in A1M region and 6 polymorphisms parsed into 2LD blocks in Bikunin region. Three major haplotypes were constructed in A1M region and two major haplotypes were constructed in Bikunin region (frequency < 0.1) (Fig. 2).

http://dam.zipot.com:8080/sites/jabg/images/N0270020304/Figure_JABG_02_03_04_F2.jpg

Figure 2. Haplotypes and LD coefficients in AMBP gene. (a) LD blocks comprised 3 polymorphisms of A1M region are shown using Haploview and haplotype frequencies. (b) LD blocks among 8 polymorphisms of Bikunin region are shown using Haploview and haplotype frequencies. The color code on the Haploview plot follows the standard color scheme: shades of red (|D’|< 1, LOD≥2); bright purple (|D'|< 1, LOD< 2).

To investigate the effects of SNPs, we analyzed the association of AMBP genotypes to determine the effects on variations in economic traits. Both the T148C and the T179C SNPs were associated with a significant effect on BF, LMA and MS. Hanwoo with the CC genotype displayed a significantly thicker BF than cattle with genotype TC, Hanwoo with the TT had a larger LMA than those with CC, Hanwoo with TT had a higher MS than those with TC (P< 0.05). The T241A significantly affected MS. Hanwoo with TT and TA had a higher MS than those with AA (P<0.05). The G9457A significantly affected BF and CW. Hanwoo with GG genotype displayed a significantly thicker BF than cattle with GA, Hanwoo with AA had a heavier CW than those with GG (P<0.05). The C9662T significantly affected CW. Hanwoo with CT had heavier CW than those with CC (P<0.05). And G9936A, T9939C, G9998A significantly affected BF (P<0.05) (Table 4).

Table 4. Least square meam and standard error analysis for economic traits in AMBP gene

http://dam.zipot.com:8080/sites/jabg/images/N0270020304/Table_JABG_02_03_04_T4.jpg

CW=carcass weight; LMA=longissimus dorsi muscle area; BF=backfat thickness; MS=marbling score.

a,b,c Different superscripts within columns differ significantly (P <0.05).

Discussion

AMBP gene is located on bovine chromosome 8 and SNPs in AMBP gene, a gene in the QTL region, are associated with meat quality traits in swine (Cinar et al., 2012). Therefore, this study aims to investigate the AMBP gene in association to economic traits in Hanwoo.

We show 14 SNPs were detected in the AMBP gene. The T148C displayed a change in amino acid from valine to alanine (Table 2). Based on our statistical association analysis, T148C and T179C were the same in genotype. They also significantly affect BF, LMA and MS. The T241A significantly affect MS, the G9457A significantly affect BF and CW, the C9662T significantly affect CW and G9936A, T9939C, G9998A significantly affect BF (P<0.05) (Table 4). And 3 polymorphisms parsed into LD block in A1M region. Previous studies have reported that several SNPs in the QTL region that influence carcass traits (CW, LMA, MS) are mapped at the region on BTA8 where AMBP gene is located (McClure et al., 2010). In pig, QTLs for meat quality traits (meat color, water content, drip loss and pH and conductivity) have been mapped in the region on SSC1 where AMBP gene in located (Liu et al., 2009). The SNP of AMBP gene was associated with meat quality traits in pig (P<0.05) (Cinar et al., 2012). And a partial protein sequence alignment of amino acid sequences from cattle, pig, human and rat, the alanine altered by T148C is conserved in all species (Fig. 3). This conservation suggests that T148C may alter economic traits in Hanwoo. These results provide supporting evidence that the AMBP gene may be a marker for carcass traits in Hanwoo.

http://dam.zipot.com:8080/sites/jabg/images/N0270020304/Figure_JABG_02_03_04_F3.jpg

Figure 3. AMBP Gene structure and detected SNPs in the Hanwoo cattle. Bos : Bos taurus (Cattle), Sus : Sus scrofa (Pig), Homo : Homo sapiens (Human), Mus : Mus musculus (Rat). * : same sequence, red frame : indicating residue of valine associated with missense T148C mutation.

These results suggest that the SNPs on the AMBP gene influence economic traits in Hanwoo. And it also provide evidence that the AMBP is a candidate gene for economic traits.

References

1 Atapattu NSBM, Saranga K, Abeywickrama LM, Gunawardane WWDA and Munasinghe M. 2016. Contribution of Backyard Chicken Egg Production to Fulfill Household Nutrient Requirement: A Study in Southern Sri Lanka.Int. J. Livst. Res (7). 

2 Sri Lanka poultry firm ups export of parent stock chicks. 2017. Economynext. Available at :https:// economynext.com/Sri_Lanka_poultry_firm_ups_export_of_parent_stock_chiks-3-7758-.html [Accessed 25 Aug.2018].  

3 Bett RC, Bhuiyan AKFH, Khan MS, Silva GLLP, Thuy LT, Sarker SC, Abeykoon MND, Nguyen TTH, Sadef S, Kariuki E, Baltenweck I, Poole J, Mwai O and Ibrahim MNM. 2014. Indigenous chicken production in the South and South East Asia. Livestock Research for Rural Development. 26:229. 

4 Bogahawatta C. 1987. Seasonal, Retail, and Wholesale Price Analysis of Some Major Agricultural Commodities. Marketing & Food Policy Division, Agrarian Research and Training Institute, Sri Lanka Available at: http://dl.nsf.ac.lk/ohs/harti/22935.pdf 

5 Bulumulla PBAIK. 2013. Phenotypic and genetic characterization of Sri Lankan indigenous chicken (Gallus gallusdomesticus) and Ceylon Jungle Fowl (Gallus lafayetti), and evaluation of indigenous chicken production systems, (M.Phil Thesis). Postgraduate Institute of Agriculture, University of Peradeniya, Sri Lanka. 

6 Chandrasiri AND. 2002 Country Report on The State of Animal Genetic Resources (Sri Lanka) Submitted to the Food and Agriculture Organization of the United Nations. Veterinary Research Institute, Peradeniya, Sri Lanka. . Available at: http://www.fao.org/docrep/pdf/010/a1250e/annexes/CountryReports/SriLanka.pdf 

7 Chantong W and Kaneene B. 2011. Poultry raising systems and highly pathogenic avian influenza outbreaks in Thailand: The situation, associations, and impacts. Southeast Asian J Trop Med Public Health. 42(3): 596–608. 

8 CIC. 2017. Livestock Solutions-Poultry. Available at: http://www.cic.lk/our-segments/livestocksolutions/# livestock-2 [ Accessed 30th Aug. 2018]. 

9 Department of Animal Production and Health (DAPH), Sri Lanka.2017. Annual Report,2017. Peradeniya, Available at:http://daph.gov.lk/web/images/content_image/publications/annual_reports/ Annual_Report_2017.pdf [Accessed 30 Aug. 2018]. 

10 Department of animal production and health(DAPH) Sri Lanka. 2016. Poultry Industry-Key statistics 2007- 2016. Peradeniya. Available at: http://www.daph.gov.lk/web/images/content_image/Livestock_stat/key_stat/ keystats_2016_Poultry.pdf . [Accessed 25 Aug 2018] 

11 Department of Animal Production and Health (DAPH), Sri Lanka.2015.Livestock statistical bulletin 2015.Peradeniya. Available at: http://www.daph.gov.lk/web/images/content_image/publications/ other_publications/stat_bulletin_2015.PDF. 

12 Department of Census and Statistics, Sri Lanka. National Livestock statistics 2005-2017. Available at: http://www.statistics.gov.lk/agriculture/Livestock/LivestockStatistics.html [Accessed 25 Aug. 2018]. 

13 Gajaweer, CJ. 2015. Recent Trends in Poultry Industry in Sri Lanka (master’s thesis). Department of Agricultural Economics and Extension, Faculty of Agriculture, University of Ruhuna, Sri Lanka. 

14 Gamage D, Gamage DVS de S, Jeyaruban MG and Gunaratne, SP. 1997. On-farm Reciprocal Recurrent Selection (pedigree) of CPRS Strains. Terminal report CARP project 12/149/121. 

15 Gamage DVS de S and Gunaratne SP. 2009. Smallholder Poultry Production, Livestock for rural development and poverty reduction:Sri Lankan experience. Hector Kobbekaduwa Agrarian Research and Training Institute, Colombo-7, Sri Lanka. ISBN: 978-955-612-107-0 61-95. 

16 Gamage DVS ds S, Disanayaka KGJS and Gunaratne SP. 2005.Research Report: A Study on Performance of Backyard Poultry Management System in Western Province, Sri Lanka. 

17 Gunaratne SP, Chandrasiri ADN, Hemalatha M and Roberts JA. 1993.The Productivity and Nutrition of Village Chicken in Sri Lanka. In: Newcastle Disease in Village Chickens. Control with Thermostable Oral vaccines, Ed. P B Spradbrow. Proceedings No. 39, Australian Centre for International Agriculture Research, Canberra. 

18 Iddamalgoda, A., Sugiyama, M., Oguri, K., Arahata, K., Premasiri, H. P. 1998. A study on development of Sri Lankan poultry industry-importance of rural poultry sector. Research Bulletin-Faculty Of Agriculture Gifu University, 87-96. Available at http://repository.lib.gifu-u.ac.jp/bitstream/123456789/5544/1/ KJ00000707791.pdf 

19 Jeyaruban MG, Gamage DVS de S and Gunaratne SP. 1996.Comparative Study on the Performance of Newly Evolved Strains at Central Poultry Research Station with Some Exotic Breeds. Terminal report CARP project. 

20 Kumudinie DLN, Kumara MKCN, Wijesinghe R, Seneviratne DGSD, Jayasinghe M. 2017. Poultry Sector Forecast 2017. Deparment of Animal Production and Health, Gatambe, Peradeniya, Sri Lanka. 

21  Liyanage RP, Dematawewa CMB and Silva P. 2015. Comparative Study on morphological and morphometric features of village chicken in Sri Lanka. Trop. Agric. Res. 26(2):261-273. 

22 Niroshan TD, Demetawewa CMB, Gamage DVS de S and Srikaran S. 2004. Egg Production Characteristics of Different Ecotypes of Village Chicken in Sri Lanka. Sri Lanka Association for the Advancement of Science, Proceedings of the 60lh Annual Sessions, Colombo.  

23 Rajapaksha, RRB.2014. Egg quality traits and functional properties of eggs of some backyard chicken ecotypes in Sri Lanka (Master thesis). Department of Animal Science, Faculty of Agriculture, University of Peradeniya.  

24  Samaraweera AM, Silva P, Abeykone NDF, Ibrahim MNM, Okeyo AM and Han JL. 2014. Population Genetic Structure of Sri Lankan Backyard Chicken Flocks: Implication for Conservation and Genetic Improvement Programs. Proceedings of the 10th World Congress of Genetics Applied to Livestock Production,Vancouver, Canada. 

25  Sanjeewa MN, Liyanage RP, Vidanarachchi JK and Silva GLLP. 2011. Association between egg production and body morphology of some village chicken ecotypes in Sri Lanka. In : Proceedings of the Peradeniya University Research Sessions, University of Peradeniya, Sri Lanka 

26  Silva P, Abeykone NDF, Samaraweera AM, Han JL, Ibrahim MNM and Okeyo AM.2014. Genetic Diversity and Adaptability Exist among Backyard Poultry Populations in Sri Lanka. Proceedings of the 0th World Congress of Genetics Applied to Livestock Production, Vancouver, Canada. 

27   

28  Silva P, Liyanage RP, Senadheera S. and Dematawewa CMB. 2016. Monograph on indigenous chicken in Sri Lanka. UNEP-GEF-ILRI FAnGR Asia project, University of Peradeniya, Sri Lanka. 

29  Silva,P, Dematawewa M. and Chandrasiri N. 2010. Farm animal genetic resources.Silva, P. (Edited), Indigenous animal genetic resources in Sri Lanka –status, potential and opportunities, UNEP-GEF-ILRI FAnGR Asia Project publication.pp. 1-46. 

30  Sugiyama M, Iddamalgoda A, Oguri K and Kamiya N. 2003. Development of livestock sector in Asia: an analysis of present situation of livestock sector and its importance for future development. Gifu City Women’s College, Gifu, Japan, 52. Available at http://www.gifu cwc.ac.jp/tosyo/kiyo/zenbun52/ Development_sugiyama.pdf. [Accessed 28 Aug. 2018]. 

31 Te Velde DW. 2008. The global financial crisis and developing countries. ODI Background Note. Overseas Development Institute. London. Available at https://www.odi.org/sites/odi.org.uk/files/odi-assets/ publications-opinion-files/3339.pdf.  

32 Wickramarachchi AR., Herath HMLK, Jayasinghe-Mudalige UK, Edirisinghe JC, Udugama JMM, Lokuge LDMN and Wijesuriya W. 2017. An Analysis of Price Behavior of Major Poultry Products in Sri Lanka. J. Agric. Sci. – Sri Lanka.12(2):138-148. DOI:https://jas.sljol.info/article/10.4038/jas.v12i2.8231/  

33  Wickramarathne, SHG. 1999. Village chicken production in Sri Lanka, (M.Phil. Thesis).Postgraduate Institute of Agriculture, University of Peradeniya, Sri Lanka 

34 Wickramaratne SHG, Gunaratne SP and Chandrasiri AND.1993. Growth Performance of Village Chicken under Different Feeding Systems. Trop. Agric. Res. 5: 359-367.  

35  Wijewantha P, Kothalawala KACHA, Wickramasinghe UKM, Rannasinghe RKK, 2010. Livestock Statistical Bulletin, Department of Animal Production and Health., Available at: http://www.daph .gov.lk/ web/mages/ df/ stb201031.pdf 

36  Wong et al., 2004. International Chicken Polymorphism Map Consortium. A genetic variation map for chicken with 2.8 million single nucleotide polymorphisms.Nature,432:717-722