Email updates

Keep up to date with the latest news and content from BMC Genetics and BioMed Central.

Open Access Research article

Growth hormone-releasing hormone (GHRH) polymorphisms associated with carcass traits of meat in Korean cattle

Hyun Sub Cheong1, Du-Hak Yoon2, Lyoung Hyo Kim1, Byung Lae Park1, Yoo Hyun Choi1, Eui Ryong Chung3, Yong Min Cho2, Eng Woo Park2, Il-Cheong Cheong2, Sung-Jong Oh2, Sung-Gon Yi4, Taesung Park4 and Hyoung Doo Shin1*

Author Affiliations

1 Department of Genetic Epidemiology, SNP Genetics, Inc., Seoul, 153-803, Korea

2 National Livestock Research Institute, RDA, 441-706, Korea

3 Department of Biotechnology, Sangi University, Wonju, Kangwon Do, 220-702, Korea

4 Department of Statistics, Seoul National University, Seoul, 151-747, Korea

For all author emails, please log on.

BMC Genetics 2006, 7:35  doi:10.1186/1471-2156-7-35


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2156/7/35


Received:1 November 2005
Accepted:3 June 2006
Published:3 June 2006

© 2006 Cheong et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Cold carcass weight (CW) and longissimus muscle area (EMA) are the major quantitative traits in beef cattle. In this study, we found several polymorphisms of growth hormone-releasing hormone (GHRH) gene and examined the association of polymorphisms with carcass traits (CW and EMA) in Korean native cattle (Hanwoo).

Results

By direct DNA sequencing in 24 unrelated Korean cattle, we identified 12 single nucleotide polymorphisms within the 9 kb full gene region, including the 1.5 kb promoter region. Among them, six polymorphic sites were selected for genotyping in our beef cattle (n = 428) and five marker haplotypes (frequency > 0.1) were identified. Statistical analysis revealed that -4241A>T showed significant associations with CW and EMA.

Conclusion

Our findings suggest that polymorphisms in GHRH might be one of the important genetic factors that influence carcass yield in beef cattle. Sequence variation/haplotype information identified in this study would provide valuable information for the production of a commercial line of beef cattle.

Background

The successful application of marker-assisted selection in the commercial animal population will depend on the identification of genes, including identification of genes underlying quantitative traits, exploration of genetic polymorphisms that are involved in different phenotypes of quantitative traits, and understanding how these genes/polymorphisms interact with the environment or with other genes affecting economic traits.

The growth hormone (GH) is essential for post-natal growth and general metabolism, and also plays an important role in lactation. Current knowledge indicates that GH exerts a key influence in nutrient use [1], mammary development [2], and growth [3]. There have been several reports of association between quantitative traits in cattle, such as growth performance and carcass merit, and polymorphisms in the GH gene [4-6].

The regulation of GH synthesis and secretion is multifactorial, but the predominant regulators of GH are the hypothalamic hormones, GH-releasing hormone (GHRH), GH secretagogue (GHS), and somatostatin (SS) [7]. In spite of the functional importance of GHRH in the regulation of GH, only one PCR-restriction fragment-length polymorphism (RFLP) [8] has been reported in cattle.

In this study, we examined GHRH as one of candidate genes in meat production. We performed extensive screening of GHRH by direct sequencing to detect polymorphisms and examined genetic association with the carcass traits. Here, we present 12 polymorphisms identified in GHRH and the results of an association study with meat quantity in Korean native cattle (Hanwoo).

Results and discussion

By direct DNA sequencing, 12 polymorphisms were identified in GHRH: one in 5'UTR and 11 in introns. The locations and allele frequencies of polymorphism are shown in Table 1 and Figure 1. By pair-wise linkage analysis with DNA from the 24 unrelated Korean cattle, which were used for direct sequencing, we have found that two sets of polymorphisms were in absolute LDs (|D'| = 1 and r2 = 1). Several sets of polymorphisms in complete LDs (|D'| = 1 and r2 ≠ 1) were also identified (Figure 1 and Table 2).

thumbnailFigure 1. Map of SNPs in GHRH on chromosome 13. The exons are marked by black block, and 5' and 3' UTRs indicated by open blocks. First base of translational site is denoted as nucleotide +1. Asterisks (*) indicate polymorphisms genotyped in a larger Korean cattle (n = 428). The minor allele frequency based on 24 sequencing samples only, which is different with minor allele frequency of absolutely linked SNP genotyped in larger population.

Table 1. Genotype and allele frequencies of 12 polymorphisms detected in GHRH

Table 2. Linkage disequilibrium coefficient (|D'| and r2) among GHRH SNPs

Among identified polymorphisms, six SNPs (-4241A>T, -3195T>A, -618T>A, +114C>A, +2042A>G and +2279C>T) were selected for larger-scale genotyping based on LDs (only one polymorphism if there are absolute LDs (r2 = 1)) and frequencies (> 0.1). Five major haplotypes (freq. > 0.1) were constructed (Table 3). Minor allele frequencies of SNP are shown in Table 1.

Table 3. Haplotypes and frequencies of GHRH among Korean native cattle.

Associations of GHRH polymorphisms with carcass traits were analyzed using the mixed effect model with sire and age as covariates. Sire was treated as a random effect and age a fixed effect. The obtained P values were corrected for multiple testing by the effective number of independent marker loci (5.35) in GHRH. In addition, the permutation tests were performed for adjusting P values and for controlling the false discovery rate (FDR). Among five common haplotypes (freq. > 0.1) identified, ht3 was not used for further analysis because it was almost completely (> 93%) tagged by a single SNP, +2042A>G (Table 3). Haplotypes (ht6 – ht9) with frequencies less than 0.1 were not analyzed either.

By statistical analyses,-4241A>T showed significant associations with cold carcass weight (CW) and longissimus muscle area (EMA) and the genetic effects of -4241A>T were in gene dose dependent manner. CW and EMA were highest in "T" allele homozygotes (CW = 321.7 ± 36.4 and EMA = 78.2 ± 9.3), intermediate in "A/T" heterozygotes (CW = 315.3 ± 33.1 and EMA = 75.7 ± 8.7), and lowest in "A" allele homozygotes (CW = 306.2 ± 32.5 and EMA = 74.1 ± 8.2) (PCOR = 0.025 and 0.046 for CW and EMA, respectively). The simple corrected PCOR and the permutation based PWY provided consistent results, though there were some slight differences. The QSAM represents the FDR value showing the consistent results. Similar mixed model was fit for the haplotype association analysis. However, none of the haplotypes were found to be significantly associated with CW and EMA (Table 4).

Table 4. Association analyses of the GHRH polymorphisms with carcass traits (CW and EMA) among Korean native cattle

Linear growth in vertebrate organisms is highly dependent on the GH [9]. GH is a lipolytic hormone, activating lipase and thereby mobilizing fat from adipose tissue [10]. While the primary action of GH is to stimulate skeletal and visceral growth, it has important metabolic actions as well. As a consequence of these actions, deficiency of GH can result in lowered growth rate [11], delayed bone maturation [12], decreased body mass [13], and hypoglycemia [14].

In a recent study, a GHRH polymorphism (AluI) was associated with yield traits (the average daily gain and expected progeny difference for fat thickness) in landrace pigs [15]. In this study, we also found that polymorphisms in GHRH, one of the predominant regulators of GH releasing [16], were associated with carcass yield traits (CW and EMA) in Korean native cattle.

Although the mechanisms involved in the association of alternative genotypes in the UTR and intronic SNPs with CW and EMA are not currently understood, the crucial role of the non-coding portion of genomes is now widely acknowledged. Polymorphisms within introns can affect gene function by affecting both the splice donor or acceptor site or regions nearby and regulatory motifs within introns. And UTRs are involved in many post-transciptional regulatory pathways that control mRNA localisation, stability, translation efficiency and initiation of protein synthesis. The post-transcriptional events play an important, yet incompletely understood, role in regulatory gene expression and cellular behaviour; many of the identified cis-acting elements for translational regulation occur within the UTR [17].

Conclusion

We have identified 12 polymorphisms in GHRH, and six polymorphic sites were selected for genotyping in Korean native cattle. Statistical analysis revealed that GHRH -4241A>T showed significant association with carcass traits. Replication of our finding in an independent dataset and/or functional validation of polymorphisms should be performed in the future.

Methods

Animals and phenotypic data

The Korean native cattle genomic DNA samples were obtained from 428 steers produced from 76 sires used in progeny testing program of National Livestock Research Institute (NLRI) of Korea. The dams were inseminated randomly with young sires. All steers were fed for 731.39 ± 16.53 days period under tightly controlled feeding program in Daekwanryeong and Namwon branch of NLRI. Live weights were determined before slaughter. Mean of live weights was 539.93 ± 51.96 kg. Yield grades for carcasses were determined by CW and EMA. After a 24-h chill, CW was measured, and then the left side of each carcass was cut between the last rib and the first lumbar vertebrae to determine EMA. EMA was determined at the surface using a grid [18]. Means of carcass traits were 311.47 ± 33.39 kg (CW) and 75.16 ± 8.62 cm2 (EMA).

Sequencing analysis of GHRH

We have sequenced the 9 kb full genome, including the promoter region (1.5 kb), to discover variants in 24 unrelated Korean native cattle using the ABI PRISM 3730 DNA analyzer (Applied Biosystems, Foster City, CA). Nineteen primer sets for the amplification and sequencing analysis were designed based on GenBank sequences (Ref. Genome seq.; AF242855 released on July 30, 2000). Primer information is available on website [19]. Sequence variants were verified by chromatograms.

Genotyping by single-base extension (SBE) and electrophoresis

For genotyping of polymorphic sites, amplifying and extension primers were designed for single-base extension (SBE) [20]. Primer extension reactions were performed with the SNaPshot ddNTP Primer Extension Kit (Applied Biosystems). To clean up the primer extension reaction, one unit of SAP (shrimp alkaline phosphatase) was added to the reaction mixture, and the mixture was incubated at 37°C for 1 hour, followed by 15 min at 72°C for enzyme inactivation. The DNA samples, containing extension products, and Genescan 120 Liz size standard solution was added to Hi-Di formamide (Applied Biosystems) according to the recommendations of the manufacturer. The mixture was incubated at 95°C for 5 min, followed by 5 min on ice, and then electrophoresis was performed using the ABI Prism 3100 Genetic Analyzer. The results were analyzed using the program of ABI Prism GeneScan and Genotyper (Applied Biosystems). Probe information is available on website [19].

Statistics

The X2 tests were used to determine whether individual variants were in equilibrium at each locus in the population (Hardy-Weinberg equilibrium). We examined a widely used measure of linkage disequilibrium between all pairs of biallelic loci, Lewontin's D' (|D'|) [21], and r2. Haplotypes and their frequencies were inferred using the algorithm developed by Stephens et al [22]. Phase probabilities of each site were calculated for each individual using this software. Association analyses with carcass traits (CW and EMA) were performed using a mixed effect model treating "sire" as a random effect. Age at slaughter was also included in the model. Other covariates were not available for this analysis. We fit a full model that includes all six SNPs in the model. We think the full model is more appropriate for controlling the closely linked SNPs more effectively. The effective number of independent marker loci in GHRH was calculated to correct for multiple testing. The effective number in GHRH was calculated using the software SNPSpD [23], which is based on the spectral decomposition (SpD) of matrices of pair-wise LD between SNPs. The resulting number of independent marker loci was applied to correct for multiple testing [24]. In addition, we performed the permutation test by controlling the P values by the Westfall and Young's method [25] and by controlling the FDR [26]. For the haplotype analyses, we fit the model including four haplotypes with the same covariates and performed the permutation test in a similar manner.

Authors' contributions

HS Cheong and HD Shin carried out the statistical analyses and drafted the manuscript. LY Kim, BL Park and YH Choi carried out the sequencing, SNP discovery and genotyping. DH Yoon, ER Chung and EW Park participated in the animal and phenotypic data collection. YM Cho participated in the design of the study. IC Cheong and SJ Oh participated in its design and coordination and helped to draft the manuscript. SG Yi and T Park carried out further statistical analyses. All authors read and approved the final manuscript.

Acknowledgements

This work was supported by a grant from BioGreen 21 Program, Rural Development Administration, Republic of Korea. The work of T Park was partially supported by the Natioanal Research Laboratory Program of Korea Science and Engineering Foundation.

References

  1. Bauman DE: Bovine somatotropin: review of an emerging animal technology.

    J Dairy Sci 1992, 75(12):3432-3451. PubMed Abstract | Publisher Full Text OpenURL

  2. Sejrsen K, Foldager J, Sorensen MT, Akers RM, Bauman DE: Effect of exogenous bovine somatotropin on pubertal mammary development in heifers.

    J Dairy Sci 1986, 69(6):1528-1535. PubMed Abstract OpenURL

  3. Bass JJ, Oldham JM, Hodgkinson SC, Fowke PJ, Sauerwein H, Molan P, Breier BH, Gluckman PD: Influence of nutrition and bovine growth hormone (GH) on hepatic GH binding, insulin-like growth factor-I and growth of lambs.

    J Endocrinol 1991, 128(2):181-186. PubMed Abstract OpenURL

  4. Falaki M, Gengler N, Sneyers M, Prandi A, Massart S, Formigoni A, Burny A, Portetelle D, Renaville R: Relationships of polymorphisms for growth hormone and growth hormone receptor genes with milk production traits for Italian Holstein-Friesian bulls.

    J Dairy Sci 1996, 79(8):1446-1453. PubMed Abstract | Publisher Full Text OpenURL

  5. Grochowska R, Sorensen P, Zwierzchowski L, Snochowski M, Lovendahl P: Genetic variation in stimulated GH release and in IGF-I of young dairy cattle and their associations with the leucine/valine polymorphism in the GH gene.

    J Anim Sci 2001, 79(2):470-476. PubMed Abstract | Publisher Full Text OpenURL

  6. Lagziel A, Lipkin E, Ezra E, Soller M, Weller JI: An MspI polymorphism at the bovine growth hormone (bGH) gene is linked to a locus affecting milk protein percentage.

    Anim Genet 1999, 30(4):296-299. PubMed Abstract | Publisher Full Text OpenURL

  7. Lin-Su K, Wajnrajch MP: Growth Hormone Releasing Hormone (GHRH) and the GHRH Receptor.

    Rev Endocr Metab Disord 2002, 3(4):313-323. PubMed Abstract | Publisher Full Text OpenURL

  8. Moody DE, Pomp D, Barendse W: Rapid communication: restriction fragment length polymorphism in amplification products of the bovine growth hormone-releasing hormone gene.

    J Anim Sci 1995, 73(12):3789. PubMed Abstract | Publisher Full Text OpenURL

  9. McDowell RS, Elias KA, Stanley MS, Burdick DJ, Burnier JP, Chan KS, Fairbrother WJ, Hammonds RG, Ingle GS, Jacobsen NE: Growth hormone secretagogues: characterization, efficacy, and minimal bioactive conformation.

    Proc Natl Acad Sci U S A 1995, 92(24):11165-11169. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Heffernan M, Summers RJ, Thorburn A, Ogru E, Gianello R, Jiang WJ, Ng FM: The effects of human GH and its lipolytic fragment (AOD9604) on lipid metabolism following chronic treatment in obese mice and beta(3)-AR knock-out mice.

    Endocrinology 2001, 142(12):5182-5189. PubMed Abstract | Publisher Full Text OpenURL

  11. de Graaff LC, Mulder PG, Hokken-Koelega AC: Body proportions before and during growth hormone therapy in children with chronic renal failure.

    Pediatr Nephrol 2003, 18(7):679-684. PubMed Abstract | Publisher Full Text OpenURL

  12. Johannsson G, Bengtsson BA: Growth hormone and the acquisition of bone mass.

    Horm Res 1997, 48 Suppl 5:72-77. PubMed Abstract OpenURL

  13. Bonert VS, Elashoff JD, Barnett P, Melmed S: Body mass index determines evoked growth hormone (GH) responsiveness in normal healthy male subjects: diagnostic caveat for adult GH deficiency.

    J Clin Endocrinol Metab 2004, 89(7):3397-3401. PubMed Abstract | Publisher Full Text OpenURL

  14. Perras B, Schultes B, Schwaiger R, Metz C, Wesseler W, Born J, Fehm HL: Growth hormone-releasing hormone facilitates hypoglycemia-induced release of cortisol.

    Regul Pept 2002, 110(1):85-91. PubMed Abstract | Publisher Full Text OpenURL

  15. Franco MM, Antunes RC, Silva HD, Goulart LR: Association of PIT1, GH and GHRH polymorphisms with performance and carcass traits in Landrace pigs.

    J Appl Genet 2005, 46(2):195-200. PubMed Abstract | Publisher Full Text OpenURL

  16. Vance ML: Growth-hormone-releasing hormone.

    Clin Chem 1990, 36(3):415-420. PubMed Abstract OpenURL

  17. Conklin D, Jonassen I, Aasland R, Taylor WR: Association of nucleotide patterns with gene function classes: application to human 3' untranslated sequences.

    Bioinformatics 2002, 18(1):182-189. PubMed Abstract | Publisher Full Text OpenURL

  18. APGS: Report of business for animal products grading. Animal products grading system, National Livestock Co-operatives Federation, Korea.

    1995.

  19. Primer and probe information. [http://www.snp-genetics.com/reference/Supplementary information to GHRH.doc] webcite

  20. Vreeland WN, Meagher RJ, Barron AE: Multiplexed, high-throughput genotyping by single-base extension and end-labeled free-solution electrophoresis.

    Anal Chem 2002, 74(17):4328-4333. PubMed Abstract | Publisher Full Text OpenURL

  21. Hedrick PW: Gametic disequilibrium measures: proceed with caution.

    Genetics 1987, 117(2):331-341. PubMed Abstract | Publisher Full Text OpenURL

  22. Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data.

    Am J Hum Genet 2001, 68(4):978-989. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Nyholt DR: A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other.

    Am J Hum Genet 2004, 74(4):765-769. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Nyholt DR: A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other.

    Am J Hum Genet 2004, 74():765-769. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Westfall PH, Young SS, John Wiley & Sons: Resampling-based multiple testing: examples and methods for p-value adjustment.

    1993.

  26. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing.

    JRSS series B 1995, 57:289-300. OpenURL