| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Medical Research Council Epidemiology Resource Center (E.M.D., H.E.S., C.C.) and Human Genetics Division (S.R., A.V., I.N.M.D.), University of Southampton, Southampton General Hospital, Southampton, United Kingdom SO16 6YD
Address all correspondence and requests for reprints to: Dr. Elaine Dennison, Medical Research Council Resource Center, Southampton General Hospital, Southampton, United Kingdom SO16 6YD. E-mail: emd{at}mrc.soton.ac.uk.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
| Subjects and Methods |
|---|
|
|
|---|
Bone mineral density (BMD) was measured in each subject at baseline and follow-up 4 yr later by dual energy x-ray absorptiometry at the lumbar spine and proximal femur (15). From these measurements we were able to derive annual bone loss rates both as absolute values and as a percentage of the original bone density. Short-term in vivo measurement precision for the Hologic QDR 1000 (Vertec Scientific, Ltd., Calleva Park, Reading, UK), expressed as the coefficient of variation, was 1.1% for lumbar spine BMD and 1.8% for femoral neck BMD; these figures were obtained by scanning six volunteers who were not part of the study, undergoing five scans on the same day, getting on and off the table between examinations. The long-term (4-yr) precision for the QDR 1000 using a spine phantom was less than 1%. Thoracolumbar radiographs taken at baseline were assessed for osteoarthritis according to the Kellgren and Lawrence score grading system (16).
Genomic DNA was extracted from 5 ml potassium-EDTA venous blood, quantified by picoGreen assay, and concentrations were equalized. Long-term stock DNA aliquots were made, and working 96-well plates of DNA dilutions to 7 ng/µl were prepared. Degenerate oligo primer amplifications were made from dilution plates to conserve stock DNA, and 96- or 384-well PCRs were performed from degenerate oligo primer-DNA, representing 0.07 ng of the original genomic DNA. A 384-well, allele-specific PCR (using flanking primers and internal allele-specific primers) in conjunction with a 384-well, microplate array, diagonal gel electrophoresis; fluorescent image screening; and PhoretiX gel image analysis were used for SNP analyses (17). Two single nucleotide variants in the GH gene (coded GHV1 and GHV4) were analyzed. GHV1 was in the promoter region of the GH gene, whereas GHV4 lay in the fourth exon downstream of this. GHV1 represents an adenine to guanine translocation at position 5157 in GenBank sequence J03071; GHV4 is located in intron 4 of the GH gene at position 6331 in GenBank sequence J03071 and represents a thymidine to adenine transversion.
To characterize the circulating GH profile of a subgroup of our cohort in more detail, 35 men and 36 women were admitted to hospital for overnight rest, and an indwelling iv cannula was inserted. At 0730 h, a venous blood sample was drawn for measurement of serum GH concentration, and samples were drawn every 20 min thereafter for 24 h. Standard meals were taken at 0800, 1230, and 1800 h, and a normal daily routine was encouraged within the confines of the hospital. At 0600 h, an additional blood sample was drawn for the measurement of IGF-I, IGF-binding protein-1 and -3, and GH-binding protein. An overnight urine collection was also made. Women taking hormone replacement therapy and those with diabetes mellitus were excluded.
GH was measured using a chemiluminescence technique (Nichols Institute, San Juan Capistrano, CA). The within-assay coefficients of variation (CVs) were 5.5%, 6.8%, and 10.5% at serum concentrations of 0.4, 10.0, and 18.6 mU/liter, respectively, and the between-assay coefficients of variation were 12.1%, 12.3%, and 9.0% at serum concentrations of 3.3, 6.3, and 18.0 mU/liter, respectively. The sensitivity of the assay was 0.05 mU/liter. Standards in the assay were calibrated against the International Reference Preparation (80/505).
GH profiles were analyzed to derive the following values. The circulatory profile of GH was characterized by estimating the concentration at or below which 5%, 50%, or 95% of the time is spent from the total of 72 values measured in each subject (11).
Serum IGF-I was measured using an in-house polyclonal RIA with acid-alcohol extraction. The sensitivity of the assay was 0.07 U/ml. The within-assay CVs were 11.3%, 6.5%, and 4.7% at serum concentrations of 0.23, 1.23, and 3.53 U/ml, respectively; the between-assay CVs were 10.5%, 12.1%, and 5.1% at concentrations of 0.38, 0.99, and 3.54 U/ml, respectively. IGF-binding protein-3 was measured using an RIA kit (Diagnostic Systems Laboratories, Inc., Webster, TX). The sensitivity of the assay was 0.9 ng/ml. The within-assay CVs were 8.1% and 5.4% at serum concentrations of 2200 and 7800 ng/ml, respectively, whereas the between-assay CVs were 9.8% and 4.8% at serum concentrations of 2200 and 8500 ng/ml, respectively. Estradiol was measured by a competitive RIA using 125I tracer (Diagnostic Products Corp., Los Angeles, CA). The between-assay CVs were 10%, 7.9%, and 7.7% at 230, 1003, and 1456 pmol/liter, respectively.
Serum osteocalcin was measured in a fasting blood sample using an RIA with antiserum raised to human osteocalcin and human osteocalcin as the standard and tracer (Nichols Institute). The intraassay CV was 3%; the interassay CV was 5% at 4.8 mg/liter, and the detection limit was 0.3 mg/liter. A 2-h morning urine sample was also obtained, after overnight fasting. Type I collagen cross-linked N-telopeptide (NTX/Cr) was measured using an ELISA (Ostex, Inc., Seattle, WA). The intraassay CV was 3.8%, and the interassay coefficient of variation was 6.2%. Creatinine was quantified in urine samples by the Jaffe procedure, and cross-link excretion was expressed as a ratio to urinary creatinine concentration.
The normality of variables was assessed, and variables were transformed as required. All statistical analyses were carried out using STATA (release 7.0, Stata Corporation, College Station, TX). ANOVA and linear regression were used to examine the associations between genotype and both BMD and bone loss rate at the lumbar spine and proximal femur. Multiple linear regression was performed to adjust for the effects of adult lifestyle variables that influence bone density and bone loss. Some genotypes were rare; hence, we performed some analyses dichotomizing for the presence/absence of an allele. GH haplotypes were inferred using the PHASE package (18). Bone phenotypes were related to phased GH gene haplotypes using multiple linear regression based on haplotype count variables (coded using STATA). The magnitude (D') and the significance of pairwise linkage disequilibrium (LD) between GH1 polymorphisms were computed by the 2LD program (36). D' is defined as D/Dmax, where D = f(AB) pu, with f(AB) being the frequency of the AB haplotype, and p and u being the allele frequencies at the two loci. Dmax is minimum [p(1 u), (1 p) u], when D > 0, or minimum [pu, (1 p)(1 u)], when D < 0. D' ranges from 1 to 1, 1, and 1, indicating complete LD, and 0, indicating complete lack of LD. The significance of LD was determined from
2, with
2 being N D2/p(1 p) u(1 u), where N is the sample size.
All subjects gave informed written consent, and permission to conduct the study was granted by the East and North Hertfordshire ethical committees.
| Results |
|---|
|
|
|---|
|
Table 2
shows the baseline bone density and annualized bone loss rate among men and women in the study according to allelic variation in the GHV1 and GHV4 SNPs. Observations were similar among men and women, and data are presented for both sexes combined after adjusting for age, gender, calcium intake, body weight, social class, smoking status, alcohol consumption, physical activity, and osteoarthritis (Kellgren and Lawrence score) (16). There were no statistically significant associations between either genetic marker with baseline bone density at the lumbar spine or proximal femur. There were, however, significant associations of each marker with the rate of bone loss at the lumbar spine and proximal femur. The 2 allele at the GHV1 locus was associated with significantly greater loss rates at the spine (P = 0.03) and hip (P = 0.04), whereas the 2 allele at the GHV4 locus was associated with a significantly greater loss rate at the spine (P = 0.008).
|
Haplotype analysis revealed a total of seven commonly occurring haplotypes for the two SNPs of interest (GHV1 and GHV4): 1 2 1; 2 2 1; 2 2 2; 1 2 2; 2 1 2; 2 1 1; and 1 1 1. Haplotype 1 (1 2 1) was set as the reference. Significant associations were observed between haplotypes containing the 2 allele at GHV1 and/or GHV4, and reduced bone density (as well as accelerated bone loss) at the lumbar spine, replicating our previous findings. Results at the proximal femur were again weaker and were not statistically significant. Hence, our findings suggest that the genotypes GHV1 and GHV4, or polymorphic genes in LD, are acting as allelically associated markers within the same haplotype block and are associated with accelerated bone loss rate in the seventh decade.
Thirty-five men and 36 women also underwent detailed characterization of their circulating GH profiles, as previously described. We have previously observed differing profiles of GH concentration according to gender in this cohort (19) and hence analyzed relationships with genotype in each gender separately (Table 3
). Among men, the skeletally unfavorable 2 allele at the GHV1 locus was associated with lower trough (P = 0.03), median (P = 0.06), and peak (P = 0.07) 24-h GH concentrations, as displayed graphically in Fig. 1
. In women, although analyses were also limited by the small sample size, there did not appear to be any association between GH genotype and circulating GH concentration. Similarly, there were no relationships observed between GHV4 genotype and circulating GH concentration. The sexual dimorphism observed with regard to the relationships between GHV1 genotype and GH concentration was mirrored in relationships between bone loss rate and GH concentration; relationships were significant in men, but not in women, at the total femoral region with trough (r = 0.39; P = 0.02) and median (r = 0.43; P = 0.01) concentrations. Similar relationships were observed at the femoral neck. We found no evidence of gender-GHV1 allelic interaction in the prediction of BMD or bone loss rate at any site.
|
|
Finally, we investigated the possible interaction between genotype and early environment as predictors of bone mass and bone loss (Table 4
). These analyses were performed in view of the associations previously demonstrated between another genotype associated with bone mass (the VDR gene) and early life (12). The table shows the annualized change in lumbar spine bone density according to weight at 1 yr and GHV1 genotype. Individuals in the lowest third of the distribution of weight at 1 yr showed a strong relationship between GHV1 genotype and bone loss rate (difference in mean bone loss rate, 0.94% between 11 and 22 genotypes), whereas these relationships were not apparent for subjects in the intermediate or higher thirds of the infant weight distribution. Test for interaction was significant (P = 0.04, after adjusting for potential confounding variables).
|
| Discussion |
|---|
|
|
|---|
There were a number of limitations to this study. The individuals studied were selected because we have accurate records of their early life; these subjects were all born in Hertfordshire, UK, and still live there. They have, however, previously been shown to have anthropometric and lifestyle characteristics similar to those of the general population (20). Not all subjects who participated in the baseline study (224 men and 186 women) were seen at follow-up 4 yr later, but we found no difference in baseline body build, lifestyle, and BMD between those who did and did not participate in the whole study (15). The lesser relationships observed at the lumbar spine may reflect arthritic changes; we attempted to adjust for this by including an osteoarthritis score according to the methods of Kellgren and Lawrence (16) in our regression models, although another approach may have been to examine relationships with the L1 vertebra only, because this vertebra is often uninvolved in osteoarthritic change. The study performed was small in comparison with many genetic studies, and the allelic frequency of some genotypes limited power for some analyses. Hence, although this study had 77% power to detect a significant difference in lumbar spine BMD change between the two homozygotes, 11 and 22, of the GHV4 genotype, this figure fell to 36% power when looking for a significant difference in femoral neck BMD change between the two homozygotes, 11 and 22, of the GHV4 genotype. Hence, our findings are all the more striking and require validation in a larger study powered to examine genetic effects on bone mass.
GH stimulates linear growth in childhood and bone remodeling throughout life. It stimulates chondrocytes in the growth plate to secrete IGF-I, which, in turn, signals the chondrocytes to differentiate, leading to cartilage formation and linear growth. In addition, GH probably has effects independent of IGF-I. Although the role of GH in the risk of osteoporosis is unclear, GH deficiency is associated with low adult BMD (19, 20, 21, 22, 23), although this is often difficult to correct with GH treatment (24, 23). Case-control studies also show lower circulating IGF-I levels in patients with osteoporosis than in normal controls (28). The pulsatile nature of GH secretion is controlled by two hypothalamic hormones with opposing effects: GHRH and somatostatin or GH release-inhibiting hormone (29). The biological function of this pulsatile pattern is unknown, but the amplitude of GH peaks correlate with height during childhood (30, 31), suggesting that the GH concentration may be an important determinant of skeletal growth.
The GH1 gene is located on chromosome 17q23; although little is known about the role of polymorphisms in the gene, case reports of mutations in the GH receptor have been described (32, 33, 34). Other work in animals has focused on the effects of GH oversecretion on bone modeling and remodeling in transgenic mice (35); this reveals that continual systemic stimulation by GH results in an initial increase in bone mass, but that bone quality is unfavorably affected. The proximal promoter region of the GH1 gene is highly polymorphic, containing at least 15 SNPs manifesting in 40 different haplotypes. Recently reported functional studies using an in vitro assay of reporter gene expression have suggested that 12 of these haplotypes are associated with significantly reduced expression, whereas 10 are associated with increased levels (14). Although measurements of skeletal status were not available in that study, a relationship was reported between an individuals GH1 promoter haplotype combination and adult height. The results of both studies suggest that this locus warrants closer scrutiny in explaining the complex interindividual differences in skeletal size, bone density, and bone loss.
Separate studies have related GHV1 to GH transcription in vitro (14) and GHV4 to growth retardation (13). We have examined both markers in this study; their strong LD (>60%; P < 1020) explains the similarity of their associations with bone loss. It also permits comparison between our study of osteoporosis phenotypes and GH secretion and previous reports of GH gene function in vitro and GH1 sequence polymorphism in growth retardation. Additional analyses in this cohort demonstrate strong LD not only across the GH gene, but also extending to the adjacent placental GH2 and lactogen gene cluster, and moderate LD as far as the angiotensin-converting enzyme gene locus situated an estimated 1 megabase away (34). The most obvious hypothetical chain of causality is that GH gene diversity influences GH expression and circulating levels (quantitative on allele), with consequent unfavorable effects on skeletal health, most notably in homozygotes for low secretion. However, strong LD across not only GH1, but the entire GH1-CSH1-GH2-CSH2 gene cluster, means that identification and proof of a causal site(s) (as opposed to proxy markers and risk haplotype(s)) may be difficult.
We have previously reported data suggesting that the intrauterine environment, using birth weight as a marker, may modulate the relationship between the VDR gene and adult bone size and volumetric density (12). This study provides evidence of a significant interaction between the early environment and a second genotypic locus, such that individuals who carry the unfavorable genotype, who also grew poorly in fetal life and infancy, are susceptible to accelerated bone loss in the seventh decade. Additional studies to identify the mechanisms for this association are required.
In conclusion, our study implicates, for the first time, common diversity in the GH1 region in predisposing to osteoporosis via effects on the level of GH expression. The interaction with infant weight suggests that early environmental influences may modulate the effect of GH1 genotype on bone loss some 7 decades later. Larger replication studies, with expansion of the range of phenotypes examined and comprehensive studies of common diversity in the GH1 region, need to be established
| Acknowledgments |
|---|
| Footnotes |
|---|
The Southampton Genetic Epidemiology Research Group consists of: Xiao-he Chen, Tom R. Gaunt, Himanshu Patel, Shu Ye, Faiza Tabassum, Nigel Arden, Avan Aihie Sayer, David I. W. Phillips, and David J. P. Barker.
Abbreviations: BMD, Bone mineral density; CV, coefficients of variation; LD, linkage disequilibrium; NTX/Cr, type I collagen cross-linked N-telopeptide; SNP, single nucleotide polymorphism; VDR, vitamin D receptor.
Received January 29, 2004.
Accepted July 8, 2004.
| References |
|---|
|
|
|---|
gene polymorphisms and bone mineral density: haplotype analysis in women from the United Kingdom. J Bone Miner Res 16:128134[CrossRef][Medline]
This article has been cited by other articles:
![]() |
E. M. DENNISON, H. E. SYDDALL, K. A. JAMESON, A. A. SAYER, T. R. GAUNT, S. RODRIGUEZ, I. N.M. DAY, C. COOPER, M. A. LIPS, and the Hertfordshire Cohort Study Group A Study of Relationships Between Single Nucleotide Polymorphisms from the Growth Hormone-Insulin-like Growth Factor Axis and Bone Mass: the Hertfordshire Cohort Study J Rheumatol, July 1, 2009; 36(7): 1520 - 1526. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Giustina, G. Mazziotti, and E. Canalis Growth Hormone, Insulin-Like Growth Factors, and the Skeleton Endocr. Rev., August 1, 2008; 29(5): 535 - 559. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Giordano, M. Godi, S. Mellone, A. Petri, D. Vivenza, L. Tiradani, Y. Carlomagno, D. Ferrante, T. Arrigo, G. Corneli, et al. A Functional Common Polymorphism in the Vitamin D-Responsive Element of the GH1 Promoter Contributes to Isolated Growth Hormone Deficiency J. Clin. Endocrinol. Metab., March 1, 2008; 93(3): 1005 - 1012. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. H. D. Bassett, P. J. O'Shea, S. Sriskantharajah, B. Rabier, A. Boyde, P. G. T. Howell, R. E. Weiss, J.-P. Roux, L. Malaval, P. Clement-Lacroix, et al. Thyroid Hormone Excess Rather Than Thyrotropin Deficiency Induces Osteoporosis in Hyperthyroidism Mol. Endocrinol., May 1, 2007; 21(5): 1095 - 1107. [Abstract] [Full Text] [PDF] |
||||
![]() |
B Walpole, T D Noakes, M Collins, and R J Trent Growth hormone 1 (GH1) gene and performance and post-race rectal temperature during the South African Ironman triathlon * Commentary Br. J. Sports Med., February 1, 2006; 40(2): 145 - 150. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |