Association study between growth hormone receptor (GHR ) gene polymorphisms and obesity in Korean population

Article information

J Exerc Rehabil Vol. 12, No. 6, 632-636, December, 2016
Publication date (electronic) : 2016 December 31
doi : https://doi.org/10.12965//jer.1632844.422
College of Nursing, Sungshin Women’s University, Seoul, Korea
*Corresponding author: Seung-Ae Yang, http://orcid.org/0000-0002-6235-2752, College of Nursing, Sungshin Women’s University, 2 Bomun-ro 34da-gil, Seongbuk-gu, Seoul 02844, Korea, Tel: +82-2-920-7728, Fax: +82-2-968-0560, E-mail: ewha63@sungshin.ac.kr
Received 2016 October 31; Accepted 2016 December 13.

Abstract

A main target of growth hormone (GH) is adipose tissue in human body. The GH secretion in obesity patients is impaired. It is needless to say that growth hormone receptor (GHR) is necessary in GH hormone signaling. The purpose of the present study is to examine the association between single nucleotide polymorphisms (SNPs) and the development of obesity. A total of 211 overweight/obese subjects with a body mass index (BMI) ≥23 kg/m2 and 157 nonoverweight/obese controls with a BMI of 18.5–23.0 kg/m2 were involved in this study. Seven SNPs including the rs6451620 (intron), rs4130114 (intron), rs4410646 (intron), rs6898743 (intron), rs4394131 (intron), rs6182 (Cys440Phe), and rs6184 (Pro579Thr) and rs2229765 SNPs of GHR gene were genotyped. Genotyping was performed using custom DNA chip. SNPStats was used to calculate the odds ratio, 95% confidence interval, and P-value. The link-age disequilibrium block and haplotypes among seven SNPs were determined using Haploview version 4.2. Dominant, recessive, and log-additive genetic models were conducted for genetic analyzing. Among tested SNPs in GHR gene, rs4410646 and rs6898743 showed significant association with obesity (rs4410646, P=0.02 in dominant model and P=0.036 in log-additive model; rs6898743, P=0.039 in dominant model and P=0.044 in log-additive model). In summary, these results suggest that GHR gene polymorphisms might play a role in the development of obesity in the Korean population.

INTRODUCTION

Obesity is a severe social problem in Korea (Kang, 2013). Not only social and economic factors (Chung et al., 2016) but also other various risk factors are important factors affecting the pathophysiologic and developing process of obesity. And life style and psychological factors are also related to hormone changes in obesity (Azuma et al., 2015; Kim et al., 2016; Nordkap et al., 2016).

Growth hormone (GH) is affected by psychological stress and many physiological stresses (He et al., 2015; Karaca et al., 2016). And GH action is related to body energy expenditure and body weight (Glad et al., 2015; Haruta et al., 2015; Liu et al., 2016). It is a major hormone in insulin resistance and body weight control. Therefore, GH is widely studied in many field of studies (Huang et al., 2016; Imran et al., 2016; Karaca et al., 2016).

However, hormone signaling of GH needs growth hormone receptor (GHR) (Carter-Su et al., 2016) because GH is not a steroid hormone (Carter-Su et al., 2016), and its main function is related to signaling of GHR into a cell. GHR deficiency is related with obesity, diabetes and cancer risk is likely appears. Therefore deficient GHR may affect GH function to be less modulation signals. It may show a kind of decreased pituitary hormone situation, which cause higher risk of metabolic syndrome (Khang et al., 2016), because the studies suggest that GH function is related with GHR function. GHR deficiency is related with adiposity. One of the famous example is GHR-null (GHR[−/−]) mouse (Sackmann-Sala et al., 2014), which showed insulin sensitivity and increased adiposity.

Insulin sensitivity is another problem in obesity, and GH is a major hormone in insulin sensitivity (Galescu et al., 2016; Poggiogalle et al., 2016; Suda et al., 2016). GH deficiency may lead to obesity, fatty liver, and sarcopenia. Those are closely related keywords to obesity, and these studies suggest that polymorphisms of GHR gene may associate with obesity. However, genetic effect such as common variation of functional gene may also affect its function (Barrie et al., 2017; Rahim et al., 2016). Therefore, in this study, we investigated the relationship between genetic polymorphisms of GHR gene and obese vs. nonobese in a Korean population, to compare polymorphism effect of GHR on obesity.

MATERIALS AND METHODS

Study subjects

Three hundred sixty-eight subjects who participated in a general health check-up program were enrolled in this study. The biochemical characteristics including fasting plasma glucose, fasted glycated hemoglobin, total cholesterol, and low-density lipoprotein cholesterol were measured. The values of biochemical test are shown in Table 1. Body mass index (BMI) is calculated based on the weight (kg) divided by the square of height (m) (Yang, 2015). All subjects were classified according to the classification of Korean Society for the Study of Obesity (underweight, BMI<18 kg/m2; normal, BMI 18 to <23 kg/m2; moderately obese, BMI 23 to <25 kg/m2; obesity I, BMI 25 to <30 kg/m2; obesity II, BMI≥30 kg/m2). Two hundred eleven subjects with a BMI≥23 kg/m2 were included in overweight/obese group (mean age±standard deviation, 44.7±6.4 years) and 157 subjects with a BMI of 18.5–23.0 kg/m2 in nonoverweight/obese control group (mean age±standard deviation, 43.6±6.2 years).

Clinical data of subjects included in the study

SNP selection and genotyping

Genomic DNAs were extracted from peripheral bloods of all subjects by Roche DNA extraction kit (Roche Diagnostics Corp., Indianapolis, IN, USA) and quantified using nanodrop machine. We selected seven SNPs including the rs6451620 (intron), rs4130114 (intron), rs4410646 (intron), rs6898743 (intron), rs4394131 (intron), rs6182 (Cys440Phe), and rs6184 (Pro579Thr) and rs2229765 in GHR gene. Each SNP was genotyped using custom DNA chip. After that, each genotype of SNPs was analyzed using targeted genotyping analysis software GCOS software (Affymetrix).

Statistical analysis

To evaluate the odds ratio (OR), 95% confidence interval (CI), and P-value, SNPStats (http://bioinfo.iconcologia.net/index.php) and IBM SPSS Statistics ver. 23.0 (IBM Co., Armonk, NY, USA) were used. Genetic models (dominant [major homozygous vs. heterozygous+minor homozygous], recessive [major homozygous+heterozygous vs. minor homozygous], and log-additive [major homozygous vs. heterozygous vs. minor homozygous] models) were applied (Kim et al., 2014; Zare et al., 2013). The haplotype analysis in linkage disequilibrium (LD) was performed using Haploview 4.2. The P-value less than 0.05 was considered as statistically significant.

RESULTS

To find genetic factors involved in the development of obesity, we performed an association study between polymorphisms of GHR gene and overweight/obese. Firstly, we calculated Hardy-Weinberg equilibrium (HWE) in tested polymorphisms of GHR gene. The genotype distributions of the tested seven SNPs of GHR gene were consistent with the HWE (rs6451620, P=1.00; rs4130114, P=1.00; rs4410646, P=0.16; rs6898743, P=1.00; rs4394131, P=1.00; rs6182, and rs6184, P=1.00; rs2229765, P=1.00, data not shown).

Tables 2 and 3 showed genotype and allele distributions of seven SNPs in GHR gene. Among tested seven SNPs of GHR gene, two SNPs (rs4410646 and rs6898743) revealed the significant association with obesity (rs4410646, OR, 1.68, 95% CI, 1.09–2.60, P=0.02 in dominant model [A/A genotype vs. A/C genotype+C/C genotype] and OR, 1.38, 95% CI, 1.02–1.87, P=0.036 in log-additive model [A/A genotype vs. A/C genotype vs. C/C genotype] rs6898743, OR, 0.63, 95% CI, 0.41–0.98, P=0.039 in dominant model [C/C genotype vs. C/G genotype+G/G genotype] and OR, 0.73, 95% CI, 0.54–0.99, P=0.044 in log-additive model [C/C genotype vs. C/G genotype vs. G/G gentoype]).

Genotype analysis of polymorphisms in growth hormone receptor (GHR ) gene between overweight/obesity and control subjects

Allele analysis of polymorphisms in growth hormone receptor (GHR ) gene between overweight/obesity and control subjects

However, the other five SNPs did not show any statistical significant association with susceptibility of obesity.

In haplotype analysis, we found two haplotypes (GC and TA haplotype) of one LD block in GHR gene, which consisted of rs6182 and rs6184. But two haplotypes of GHR gene showed no statistical significant difference between the control group and the overweight/obese group (GC haplotype, P=0.22 and TA haplotype, P=0.30) (Table 4).

Haplotype analysis of polymorphisms of growth hormone receptor (GHR ) gene between overweight/obesity and control subjects

DISCUSSION

In present study, seven SNPs in GHR gene were studied, however, only two intron SNPs of them were weakly associated with obesity (rs4410646 and rs6898743). The two SNPs were not similar in major/minor allele frequencies, however, both of them showed association to obesity in the dominant and log-additive models, and not in recessive model. Among the other SNPs, four SNPs (rs4130114, rs4394131, rs6182, and rs6184) displayed similar major/minor allele frequencies, and all of them showed no association with obesity.

Association between rs4410646 SNP of GHR gene and obesity was previously found in previous study (Mong et al., 2010). They reported that rs4410646 was significantly associated with the body composition (P=0.0044) and blood pressure factor scores (P=0.00017) (Mong et al., 2010). They also found rs7703713 SNP may affect IGF-1 activity factor score. However, the SNP was not included in our study. Interestingly, in their study, C genotype was associated with less obesity. It was on the contrary to our study. Our result was that C genotype was associated with overweight/obese group.

In addition, the association between rs6898743 SNP of GHR gene and cancer were studied, and rs6898743 was associated with esophageal adenocarcinoma (McElholm et al., 2010; Ong et al., 2014). But there was no obesity study was found.

In general, our result may show only weak association of GHR SNPs and obesity susceptibility. However, some presumable results may support our results. A variation in GHR gene was studied in various anthropometry (Glad et al., 2015), and relationship to central obesity was found. However, in their discussion, deficit in GHR does not always directly show GH deficiency, because they found responses to exogeneous GH in deficient GHR. They suggested the complex and continuous lifelong relationship among GH sensitivity, GH secretion, and other metabolic mechanisms in normal GH–IGF1 axis.

Another example of complex GH sensitivity may be the white adipose tissue of GHR-null mice (Sackmann-Sala et al., 2014). They presented the insulin, leptin, adiponectin, and protein changes in enhanced insulin sensitivity (Sackmann-Sala et al., 2014). In addition, Beattie et al. (2002) reported that decreased binding affinity may not reflect biological activities in the protein of living organism.

To our knowledge, this is the first study to report the relationship between SNPs in the GHR gene and obesity by BMI. Our findings support that GHR may play a role in pathophysiologic signaling of obesity or other related diseases. And further study in a large numbers with more anthropometry will be needed to verify our findings.

Notes

CONFLICT OF INTEREST

No potential conflict of interest relevant to this article was reported.

References

Azuma, K., Adachi, Y., Hayashi, H., & Kubo, K.Y.Chronic psychological stress as a risk factor of osteoporosis. J UOEH, (2015). 37, 245–253.
Barrie, E.S., Hartmann, K., Lee, S.H., Frater, J.T., Seweryn, M., Wang, D., & Sadee, W.The CHRNA5/CHRNA3/CHRNB4 nicotinic receptor regulome: genomic architecture, regulatory variants, and clinical associations. Hum Mutat, (2017). 38, 112–119.
Beattie, J., Phillips, K., Shand, J.H., Brocklehurst, S., Flint, D.J., & Allan, G.J.Sensitivity of hybrid ovine/rat GH receptors to oGH and rat GH in transfected FDC-P1 mouse myeloid cells in vitro. Mol Cell Biochem, (2002). 238, 137–143.
Carter-Su, C., Schwartz, J., & Argetsinger, L.S.Growth hormone signaling pathways. Growth Horm IGF Res, (2016). 28, 11–15.
Chung, W., Lee, S., Lim, S.J., & Kim, J.Modifying effects of education on the association between lifestyle behaviors and the risk of obesity: evidence from South Korea. BMC Public Health, (2016). 16, 1100.
Galescu, O.A., Crocker, M.K., Altschul, A.M., Marwitz, S.E., Brady, S.M., & Yanovski, J.A.A pilot study of the effects of niacin administration on free fatty acid and growth hormone concentrations in children with obesity. Pediatr Obes, (2016). Sep. 21. [Epub]. https://doi.org/10.1111/ijpo.12184.
Glad, C.A., Carlsson, L.M., Melander, O., Almgren, P., Sjöström, L., Nilsson, S., Larsson, I., Svensson, P.A., & Johannsson, G.The GH receptor exon 3 deleted/full-length polymorphism is associated with central adiposity in the general population. Eur J Endocrinol, (2015). 172, 123–128.
Haruta, I., Fuku, Y., Kinoshita, K., Yoneda, K., Morinaga, A., Amitani, M., Amitani, H., Asakawa, A., Sugawara, H., Takeda, Y., Bowers, C.Y., & Inui, A.One-year intranasal application of growth hormone releasing peptide-2 improves body weight and hypoglycemia in a severely emaciated anorexia nervosa patient. J Cachexia Sarcopenia Muscle, (2015). 6, 237–241.
He, F.Q., Fang, G., Wang, B., Guo, X.J., & Guo, C.L.Perinatal stress effects on later anxiety and hormone secretion in male mandarin voles. Behav Neurosci, (2015). 129, 789–800.
Huang, C.C., Chou, D., Yeh, C.M., & Hsu, K.S.Acute food deprivation enhances fear extinction but inhibits long-term depression in the lateral amygdala via ghrelin signaling. Neuropharmacology, (2016). 101, 36–45.
Imran, S.A., Tiemensma, J., Kaiser, S.M., Vallis, M., Doucette, S., Abidi, E., Yip, C.E., De Tugwell, B., Siddiqi, F., & Clarke, D.B.Morphometric changes correlate with poor psychological outcomes in patients with acromegaly. Eur J Endocrinol, (2016). 174, 41–50.
Kang, Y.S.Obesity associated hypertension: new insights into mechanism. Electrolyte Blood Press, (2013). 11, 46–52.
Karaca, Z., Tanrıverdi, F., Ünlühızarcı, K., & Kelestimur, F.GH and pituitary hormone alterations after traumatic brain injury. Prog Mol Biol Transl Sci, (2016). 138, 167–191.
Khang, A.R., Ku, E.J., Kim, Y.A., Roh, E., Bae, J.H., Oh, T.J., Kim, S.W., Shin, C.S., Kim, S.Y., & Kim, J.H.Sex differences in the prevalence of metabolic syndrome and its components in hypopituitary patients: comparison with an age- and sex-matched nationwide control group. Pituitary, (2016). 19, 573–581.
Kim, J.H., Chang, S.M., Bae, J.N., Cho, S.J., Lee, J.Y., Kim, B.S., & Cho, M.J.Mental-physical comorbidity in Korean adults: results from a Nationwide General Population Survey in Korea. Psychiatry Investig, (2016). 13, 496–503.
Kim, S.K., Park, H.J., Chung, J.H., Kim, J.W., Seok, H., Lew, B.L., & Sim, W.Y.Association between interleukin 18 polymorphisms and alopecia areata in Koreans. J Interferon Cytokine Res, (2014). 34, 349–353.
Liu, Z., Cordoba-Chacon, J., Kineman, R.D., Cronstein, B.N., Muzumdar, R., Gong, Z., Werner, H., & Yakar, S.Growth hormone control of hepatic lipid metabolism. Diabetes, (2016). 65, 3598–3609.
McElholm, A.R., McKnight, A.J., Patterson, C.C., Johnston, B.T., Hardie, L.J., & Murray, L.J.Finbar Group. A population-based study of IGF axis polymorphisms and the esophageal inflammation, metaplasia, adenocarcinoma sequence. Gastroenterology, (2010). 139, 204–212.e3.
Mong, J.L., Ng, M.C., Guldan, G.S., Tam, C.H., Lee, H.M., Ma, R.C., So, W.Y., Wong, G.W., Kong, A.P., Chan, J.C., & Waye, M.M.Associations of the growth hormone receptor (GHR) gene polymorphisms with adiposity and IGF-I activity in adolescents. Clin Endocrinol (Oxf), (2010). 73, 313–322.
Nordkap, L., Jensen, T.K., Hansen, ÅM, Lassen, T.H., Bang, A.K., Joensen, U.N., Blomberg Jensen, M., Skakkebæk, N.E., & Jørgensen, N.Psychological stress and testicular function: a cross-sectional study of 1,215 Danish men. Fertil Steril, (2016). 105, 174–187.e1–2.
Ong, J., Salomon, J., te Morsche, R.H., Roelofs, H.M., Witteman, B.J., Dura, P., Lacko, M., & Peters, W.H.Polymorphisms in the insulin-like growth factor axis are associated with gastrointestinal cancer. PLoS One, (2014). 9, e90916.
Poggiogalle, E., Lubrano, C., Gnessi, L., Mariani, S., Lenzi, A., & Donini, L.M.Fatty liver index associates with relative sarcopenia and GH/IGF- 1 status in obese subjects. PLoS One, (2016). 11, e0145811.
Rahim, M., El Khoury, L.Y., Raleigh, S.M., Ribbans, W.J., Posthumus, M., Collins, M., & September, A.V.Human genetic variation, sport and exercise medicine, and achilles tendinopathy: role for angiogenesis-associated genes. OMICS, (2016). 20, 520–527.
Sackmann-Sala, L., Berryman, D.E., Lubbers, E.R., Zhang, H., Vesel, C.B., Troike, K.M., Gosney, E.S., List, E.O., & Kopchick, J.J.Age-related and depot-specific changes in white adipose tissue of growth hormone receptor-null mice. J Gerontol A Biol Sci Med Sci, (2014). 69, 34–43.
Suda, K., Matsumoto, R., Fukuoka, H., Iguchi, G., Hirota, Y., Nishizawa, H., Bando, H., Yoshida, K., Odake, Y., Takahashi, M., Sakaguchi, K., Ogawa, W., & Takahashi, Y.The influence of type 2 diabetes on serum GH and IGF-I levels in hospitalized Japanese patients. Growth Horm IGF Res, (2016). 29, 4–10.
Yang, S.A.Association between exonic polymorphism (rs629849, Gly1619Arg) of IGF2R gene and obesity in Korean population. J Exerc Rehabil, (2015). 11, 282–286.
Zare, M., Shahtaheri, S.J., Mehdipur, P., Abedinejad, M., & Zare, S.The influence of CYP1A1 and GSTM1 polymorphism on the concentration of urinary 1-hydroxypyrene in cPAHs exposed Iranian anode plant workers. Mol Cell Toxicol, (2013). 9, 303–309.

Article information Continued

Table 1

Clinical data of subjects included in the study

Characteristic Control (n=157) Overweight/obesity (n=211)
Age (yr) 43.6±6.2 44.8±6.4
Body mass index (kg/m2) 21.4±1.2 25.6±2.0
Fasting plasma glucose (mg/dL) 89.7±11.4 93.8±14.8
HbA1c (%) 5.3±0.4 5.5±0.7
Total cholesterol (mg/dL) 186.0±27.7 196.9±34.1
LDL-C (mg/dL) 108.8±28.3 119.0±32.3

Values are presented as mean±standard deviation.

HbA1c, fasted glycated hemoglobin; LDL-C, low-density lipoprotein cholesterol.

Table 2

Genotype analysis of polymorphisms in growth hormone receptor (GHR ) gene between overweight/obesity and control subjects

SNPs Genotype Control Overweight/obese Models OR (95% CI) P-value
rs6451620 intron G/G 52 (33.5) 82 (38.9) Dominant 0.77 (0.49–1.20) 0.24
A/G 76 (49.0) 102 (48.3) Recessive 0.68 (0.38–1.23) 0.21
A/A 27 (17.4) 27 (12.8) Log-additive 0.79 (0.58–1.08) 0.14

rs4130114 intron T/T 128 (83.7) 179 (87.3) Dominant 0.77 (0.42–1.40) 0.39
G/T 24 (15.7) 26 (12.7) Recessive 0.00 (0.00–NA) NA
G/G 1 (0.6) 0 (0) Log-additive 0.75 (0.42–1.35) 0.33

rs4410646 intron A/A 71 (45.2) 68 (32.2) Dominant 1.68 (1.09–2.60) 0.02*
A/C 63 (40.1) 103 (48.8) Recessive 1.31 (0.74–2.32) 0.35
C/C 23 (14.7) 40 (19.0) Log-additive 1.38 (1.02–1.87) 0.036*

rs6898743 intron C/C 52 (33.1) 96 (45.7) Dominant 0.63 (0.41–0.98) 0.039*
C/G 77 (49.0) 86 (41.0) Recessive 0.71 (0.40–1.28) 0.26
G/G 28 (17.8) 28 (13.3) Log-additive 0.73 (0.54–0.99) 0.044*

rs4394131 intron T/T 127 (81.4) 170 (81.0) Dominant 1.06 (0.62–1.83) 0.82
A/T 28 (17.9) 38 (18.1) Recessive 1.51 (0.13–17.25) 0.73
A/A 1 (0.6) 2 (1.0) Log-additive 1.08 (0.65–1.78) 0.78

rs6182 Cys440Phe G/G 134 (85.3) 187 (89.0) Dominant 0.70 (0.37–1.31) 0.26
G/T 22 (14.0) 23 (10.9) Recessive 0.00 (0.00–NA) 0.13
T/T 1 (0.6) 0 (0) Log-additive 0.67 (0.37–1.23) 0.20

rs6184 Pro579Thr C/C 133 (85.3) 179 (88.6) Dominant 0.74 (0.39–1.39) 0.35
A/C 23 (14.7) 23 (11.4) Recessive NA NA
A/A 0 (0) 0 (0) Log-additive NA NA

Values are presented as number (%).

SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval; NA, not applicable.

*

P<0.05 significant association.

Table 3

Allele analysis of polymorphisms in growth hormone receptor (GHR ) gene between overweight/obesity and control subjects

SNPs Allele Control Overweight/obese OR (95% CI) P-value
rs6451620 intron G 180 (58.1) 266 (63.0) 1
A 130 (41.9) 156 (37.0) 0.812 (0.602–1.096) 0.17

rs4130114 intron T 280 (91.5) 384 (93.7) 1
G 26 (8.5) 26 (6.3) 0.729 (0.414–1.283) 0.27

rs4410646 intron A 205 (65.3) 239 (56.6) 1
C 109 (34.7) 183 (43.4) 1.440 (1.065–1.947) 0.018*

rs6898743 intron C 181 (57.6) 278 (66.2) 1
G 133 (42.4) 142 ( 33.8) 0.695 (0.514–0.940) 0.018*

rs4394131 intron T 282 (90.4) 378 (90.0) 1
A 30 (9.6) 42 (10.0) 1.044 (0.638–1.710) 0.86

rs6182 Cys440Phe G 290 (92.4) 397 (94.5) 1
T 24 (7.6) 23 (5.5) 0.700 (0.387–1.265) 0.24

rs6184 Pro579Thr C 289 (92.6) 381 (94.3) 1
A 23 (7.4) 23 (5.7) 0.759 (0.417–1.379) 0.37

Values are presented as number (%).

SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval; NA, not applicable.

*

P<0.05 significant association.

Table 4

Haplotype analysis of polymorphisms of growth hormone receptor (GHR ) gene between overweight/obesity and control subjects

Haplotype Frequency Control Overweight/obese Chi-square P-value
+ +
GC 0.933 289.0 25.0 397.9 24.1 1.482 0.22
TA 0.062 23.0 291.0 23.0 399.0 1.079 0.30

Haplotypes consist of rs6182 and rs6184.