- Research article
- Open Access
Growth hormone-releasing hormone (GHRH) polymorphisms associated with carcass traits of meat in Korean cattle
BMC Genetics volume 7, Article number: 35 (2006)
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).
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.
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.
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 , mammary development , and growth . 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) . In spite of the functional importance of GHRH in the regulation of GH, only one PCR-restriction fragment-length polymorphism (RFLP)  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).
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.
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 PCORand the permutation based PWYprovided consistent results, though there were some slight differences. The QSAMrepresents 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).
Linear growth in vertebrate organisms is highly dependent on the GH . GH is a lipolytic hormone, activating lipase and thereby mobilizing fat from adipose tissue . 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 , delayed bone maturation , decreased body mass , and hypoglycemia .
In a recent study, a GHRH polymorphism (Alu I) was associated with yield traits (the average daily gain and expected progeny difference for fat thickness) in landrace pigs . In this study, we also found that polymorphisms in GHRH, one of the predominant regulators of GH releasing , 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 .
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.
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 . 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 . 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) . 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 .
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'|) , and r2. Haplotypes and their frequencies were inferred using the algorithm developed by Stephens et al . 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 , 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 . In addition, we performed the permutation test by controlling the P values by the Westfall and Young's method  and by controlling the FDR . For the haplotype analyses, we fit the model including four haplotypes with the same covariates and performed the permutation test in a similar manner.
Bauman DE: Bovine somatotropin: review of an emerging animal technology. J Dairy Sci. 1992, 75 (12): 3432-3451.
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.
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.
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.
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.
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. 10.1046/j.1365-2052.1999.00472.x.
Lin-Su K, Wajnrajch MP: Growth Hormone Releasing Hormone (GHRH) and the GHRH Receptor. Rev Endocr Metab Disord. 2002, 3 (4): 313-323. 10.1023/A:1020949507265.
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-
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. 10.1073/pnas.92.24.11165.
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. 10.1210/en.142.12.5182.
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.
Johannsson G, Bengtsson BA: Growth hormone and the acquisition of bone mass. Horm Res. 1997, 48 Suppl 5: 72-77.
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. 10.1210/jc.2003-032213.
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. 10.1016/S0167-0115(02)00163-5.
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.
Vance ML: Growth-hormone-releasing hormone. Clin Chem. 1990, 36 (3): 415-420.
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. 10.1093/bioinformatics/18.1.182.
APGS: Report of business for animal products grading. Animal products grading system, National Livestock Co-operatives Federation, Korea. 1995
Primer and probe information. [http://www.snp-genetics.com/reference/Supplementary information to GHRH.doc]
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. 10.1021/ac0258094.
Hedrick PW: Gametic disequilibrium measures: proceed with caution. Genetics. 1987, 117 (2): 331-341.
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. 10.1086/319501.
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. 10.1086/383251.
Nyholt DR: A simple correction for multiple testing for single-nucleotidepolymorphisms in linkage disequilibrium with each other. Am J Hum Genet. 2004, 74 (): 765-769. 10.1086/383251.
Westfall PH, Young SS, John Wiley & Sons: Resampling-based multiple testing: examples andmethods for p-value adjustment. 1993,
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical andpowerful approach to multiple testing. JRSS series B. 1995, 57: 289-300.
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.
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.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Cheong, H.S., Yoon, DH., Kim, L.H. et al. Growth hormone-releasing hormone (GHRH) polymorphisms associated with carcass traits of meat in Korean cattle. BMC Genet 7, 35 (2006). https://doi.org/10.1186/1471-2156-7-35
- Growth Hormone
- Beef Cattle
- Carcass Weight
- Carcass Trait
- Allele Homozygote