|
|
||||||||
1School of Medicine and Pharmacology, University of Western Australia; 2Department of Endocrinology and Diabetes, Sir Charles Gairdner Hospital; 3Western Australian Institute of Medical Research; and 4School of Exercise, Biomedical and Health Sciences, Edith Cowan University, Perth, Western Australia
Submitted 16 November 2004 ; accepted in final form 16 December 2004
| ABSTRACT |
|---|
|
|
|---|
estradiol; fracture; bone density
It is not known if there is an association between CYP19 polymorphisms, estradiol concentration, and measures of bone mass and vertebral fracture in older women. On the basis of previous studies in younger postmenopausal women and the demonstrated importance of endogenous estradiol concentration in determining BMD and fracture risk in older women (5), we hypothesize that the CYP19 TTTA microsatellite is associated with estrogen concentration and consequently bone density, strength, and vertebral fracture in women many years past the menopause. Therefore, we examined the association of this TTTA microsatellite of the CYP19 gene in an unselected, elderly female population with circulating estradiol concentration, measures of bone mass and strength, and vertebral deformities.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Study design. The study design consisted of a preliminary evaluation of the association of microsatellite repeats on hip BMD measured by dual-energy X-ray absorptiometry (DXA). Once the association of a 171-bp repeat had been identified at the total hip site, an a priori hypothesis testing approach was applied to the biochemical data, spine BMD, quantitative ultrasound (QUS), and peripheral quantitative computed tomography (pQCT) measurements and vertebral fracture.
Quantitative QUS, BMD, and pQCT. QUS of the calcaneum of the left foot was measured twice in 1,129 of these subjects using a Lunar Achilles Ultrasound machine (Lunar, Madison, WI). The manufacturer's quality assurance methods were employed. The average measurement of the speed of sound (SOS), broadband ultrasound (BUA), and stiffness was reported. The coefficient of variation (CV) for SOS and BUA using the manufacturer's standards were 0.43 and 1.59%, respectively. BMD measurements of the hip on 1,065 subjects were carried out using an Hologic Acclaim 4500A detector fan beam densitometer (Hologic, Waltham, MA), 12 mo after the commencement of the study. The CV at the total hip site was 1.0% (9). Lumbar spine BMD measurements were available for 211 subjects. The CV at the spine was 1.1% (9). pQCT bone structure and density were measured on 1,048 subjects at 5 yr at the radius at a site 4% of the length of the tibial shaft to the ankle, using a Stratec XCT 2000 pQCT (Stratec Medizintechnik, Pforzheim, Germany). The voxel size was set at 150 µm in the x- and y-direction and 1,000 µm in the z-direction, which increased the scan time to 5 min. Trabecular and cortical BMD was ascertained using the peel mode 1 algorithm of the manufacturer's analysis software package, version 5.50. For the tibia, the CV error for total BMD was 1.5%, for trabecular BMD was 2.9%, and for cortical BMD was 2.5%. The cross-sectional area of cortical bone was measured using a periosteal bone density threshold of 710 mg/cm3 and an endosteal threshold of 169 mg/cm3 determined as the lowest density that consistently allowed delineation of an endosteal surface in these patients. A previously validated biomechanical parameter, the stress strain index (SSI), was calculated as the product of the section modulus and cortical density normalized to the maximal physiological cortical density of human bones (1,200 mg/cm3) for the polar moment and the bending moments in the x- and y-direction, where the y-direction is the widest part of the radius and the x-direction is perpendicular to this (1).
Morphometric X-ray absorptiometry analysis. Dual-energy high-definition lateral morphometric X-ray absorptiometry (MXA) scans of vertebra T4-L4 (Hologic QDR 4500 A) were available for 976 of the subjects who had a CYP19 genotype. The positions of six reference markers for each vertebra were placed at the corners and in the middle of the upper and lower surface of each vertebra as described by the operator's manual. The mean CV in MXA measurements of the anterior and posterior heights of vertebrae L4 to T12 was 4.8% determined by the one operator who undertook the scans and point placement on duplicate measurements from 30 subjects. Vertebral heights reference data were calculated according to a modification of the procedure described by McCloskey et al. (17).
Demographic, anthropometric, and lifestyle factors. All subjects completed a lifestyle questionnaire that included age at menopause and smoking history. A quantitative food frequency questionnaire (11) was used to ascertain the average daily alcohol and calcium intake. Weight and height were measured. Body mass index (BMI) was calculated as height/(weight)2.
Biochemical assessment. A randomly selected subgroup of 242 subjects had a blood and urine sample collected after an overnight fast at baseline. The urine samples were analyzed for creatinine, calcium, and phosphorus using routine methods (BM/Hitachi 747 Analyser, Boehringer Mannheim, Mannheim, Germany). Urine deoxyprydinoline was measured by HPLC (20) and corrected for creatinine excretion. The blood samples were analyzed for alkaline phosphatase, creatinine, calcium, and phosphorus using routine methods (BM/Hitachi 747 Analyser). Serum osteocalcin was determined using RIA techniques as previously described (14, 19). Serum estradiol was measured by RIA (Orion Diagnostica, Espoo, Finland) on 1,146 subjects. This assay had a sensitivity of 5 pmol/l and an interassay CV of 20% at 18 pmol/l and 6.6% at 100 pmol/l and an intra-assay CV of 5% at 100 pmol/l. Sex hormone-binding globulin (SHBG) was measured on these subjects using an immunochemiluminometric assay (Imulite, Los Angeles, CA). The intra- and interassay CV were 7.1 and 6.8%, respectively, at 24 nmol/l. Serum intact parathyroid hormone was measured by immunochemiluminometric assay (27) with intra- and interassay CV of 3.6 and 6.2%, respectively.
Genotype analysis.
Genomic DNA was extracted and purified from EDTA whole blood samples. The tetranucleotide [TTTA]n repeat begins at the 682 bp of the human CYP19 gene and has been reported to be repeated up to 24 times (18). This region was amplified by PCR using previously described TET-labeled primers (18). Amplification conditions were as follows: 94°C for 4 min and then 94°C (30 s); 55°C (30 s); 72°C (30 s) for 5 cycles; 94°C (5 s); 55°C (30 s); 72°C (45 s) for 20 cycles; 94°C (5 s); 55°C (45 s); and 72°C (1 min 20 s) for 15 cycles followed by a 7-min extension. The PCR product was electrophoresed in preheated (55°C) 6% polyacrylamide gel for 50 min at 70 W. Samples were visualized on the Hitachi FMBIO (Tokyo, Japan) using a 585-nm filter. The smallest allele has previously been reported to not be part of the [TTTA]n polymorphism and is due to a 3 base insertion/deletion
50 bp upstream of the [TTTA]n repeat (8). This insertion/deletion polymorphism has been reported to be in complete linkage disequilibrium with a [TTTA]n of 7 in a caucasian population (2).
Statistical analysis.
Hardy Weinberg equilibrium for the CYP19 TTTA allele was tested using the Genepop program (http://wbiomed.curtin.edu.au/genepop/genepop_op1.html). A preliminary analysis of the association of CYP19 TTTA alleles with low total hip BMD (T
1) was performed using the CLUMP program, which uses a Monte Carlo simulation technique to compare the departure of observed values from expected values conditional on the marginal totals. This approach overcomes the problems of multiple comparisons (22). All subsequent statistical analyses were performed using SPSS Windows version 11 (SPSS, Chicago, IL). The subjects were divided into two groups depending on whether they had at least one allele with a TTTA allele of 7 repeats that included the 3-bp insertion 50 bp from the TTTA microsatellite (A2), which was a cutoff chosen based on the preliminary analysis of the hip BMD data and applied to all other analyses. The data were also analyzed on the basis of a dichotomized allele length of [TTTA]n of 10 or longer, or less than 10, as previously described (6). Significant differences in QUS, height, weight, and biochemical parameters between the groups were determined using unpaired Student's t-test. Significant differences in age, number of years since menopause, alcohol consumption, cigarette smoking, and calcium intake between the groups were examined using the Mann-Whitney U-test, as these variables were not normally distributed. The free estradiol index (FEI) was calculated as the molar ratio of serum estradiol concentration divided by serum sex binding globulin concentration. The log10 values for the deoxypyridonoline creatinine ratio and the FEI were calculated to normalize the distribution of these variables. QUS parameters were compared between groups using Student's t-test. BMD and pQCT were measured at 1 and 5 yr, respectively, after commencement of calcium treatment, therefore the estimated means of BMD and pQCT parameters were calculated after adjustment for calcium supplementation or placebo and the groups were compared by Wald t-test using the general linear modeling (GLM) procedure. In some analysis, BMD, QUS, and pQCT values were adjusted for the FEI and means were compared using GLM. Those subjects with a T score for hip BMD of
1 SD below the premenopausal mean (
0.855) were classified as being osteopenic and those subjects with a T score for hip BMD of
2.5 SD below the premenopausal mean (
0.675) were classified as being osteoporotic according to World Health Organization (WHO) criteria (12) applied to a reference population of young normal women collected at our laboratory. The effect of a CYP19 TTTA A2 on the distribution into the normal compared with the osteopenic and osteoporotic groups was calculated using logistic regression analysis, adjusting for the FEI and the results were reported as the odds ratio (OR) and 95% confidence intervals (CI).
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
Relationship between CYP19 [TTTA]n alleles and osteopenia and osteoporosis. In this study population, 34.5% had normal BMD, 52.7% were osteopenic, and 12.2% were osteoporotic at the total hip site according to WHO criteria. Those subjects with an A2 allele were less likely to be osteopenic (OR 0.72, 95% CI 0.520.99). The decreased risk of being osteoporotic with an A2 allele did not achieve statistical significance (OR 0.64, 95% CI 0.381.05). When the osteopenic and osteoporotic (i.e., low BMD) subjects were considered together, an A2 allele was associated with a decreased risk of low BMD (OR 0.70, 95% CI 0.520.95). After adjustment for the FEI, those subjects with an A2 allele were no longer less likely to have low BMD (OR 0.92, 95% CI 0.651.31).
The analysis based on low and high TTTAn repeats demonstrated that those subjects with a low repeat number allele were less likely to be in the osteoporotic group (OR 0.60, 95% CI 0.380.95). The association of low repeat numbers with osteopenia did not reach statistical significance (OR 0.78, 95% CI 0.581.04) compared with those subjects who had low repeat number alleles. When the osteopenic and osteoporotic groups were considered together, a low number of repeats was associated with a decreased risk of low BMD (OR 0.74, 95% CI 0.560.98).
Relationship between CYP19 [TTTA]n alleles and vertebral fracture.
Vertebral fracture, determined by MXA, was assessed on 976 subjects. Of these, 16.1% had one or more anterior wedge fractures (AW), 17.8% had one or more crush fractures, and 5.6% had one or more central collapse fractures. A combination of AW and crush fractures occurred in 3.7% of subjects. The various types of fracture were subsequently analyzed together, with 33.5% of subjects having one or more vertebral fractures. A strong interactive effect between WHO-defined bone density groups and vertebral fracture with an A2 allele was observed (Fig. 6). In the osteoporotic group, the presence of an A2 allele was associated with a lower risk of vertebral fracture (OR 0.27, CI 0.090.79). This was not observed in osteopenic or normal subjects (Fig. 6). There was no difference in the prevalence of vertebral fractures in the osteoporotic subjects who had an A2 allele compared with those subjects who had an A2 allele and had a hip BMD T score greater than 2.5 (
2 = 1.65, P = 0.2; Fig. 6). The high prevalence of vertebral fractures in the osteoporotic subjects who did not have an A2 allele was highly significantly different from those subjects who did not have an A2 allele and who had a hip BMD T score greater than 2.5 (
2 = 11.26, P = 0.001; Fig. 6). After adjustment for the FEI, the association in the osteoporotic group between the A2 allele and vertebral fracture remained significant (OR 0.21, CI 0.070.68). When a high repeat-containing allele group was compared with the low repeat allele group, no significant difference in vertebral fracture was observed, either considered as a whole group (OR 0.90, CI 0.691.18), or in the osteoporotic group alone (OR 0.53, CI 0.231.24), for the short compared with the long repeat allele groups.
|
| DISCUSSION |
|---|
|
|
|---|
A previous study examining the association of the CYP19 TTTA microsatellite with bone mass and fracture risk reported that a short TTTA repeat was associated with a higher risk of osteoporosis at the spine and vertebral fracture (16). This is opposite to the pattern observed in our study, where the shorter TTTA repeat was associated with a higher BMD. The ethnically homogeneous Italian population studied by Masi et al. (16) is considerably different from the population reported in the present study of postmenopausal women of Northern European origin. Our population was much older (75 compared with 57 yr of age) and had a much lower rate of WHO-defined osteoporosis (12.6 compared with 52.8%). As the prior clinical diagnosis of osteoporosis was an exclusion from this study, this is likely to account for the relatively low prevalence of DXA-defined osteoporosis in this population. Another population-based study that used the Danish Osteoporosis Prevention Study (DOPS) failed to find an association between CYP19 TTTA microsatellite and BMD in young postmenopausal women, but it did find an interaction with this polymorphism to the BMD response to hormone replacement therapy (29). The reason why our population-based study differed from the DOPS population-based study in the association of BMD with the TTTA microsatellite polymorphism may be due to the difference in ages of the two study cohorts or may be due to the difference in the analysis methods. The variation in results between this and other studies is similar to the discrepant results observed with the association of breast cancer and the CYP19 TTTA microsatellite polymorphism, where some studies show an increased breast cancer risk associated with a [TTTA]12 allele (7, 15), association with other allele lengths (2), or decreased breast cancer risk associated with the [TTTA]12 allele (23). The present study was comprehensive enough to allow a priori testing for effects on bone density at the spine site and QUS at the heel, together with effects on biochemical parameters and vertebral fracture, after the initial determination of the significant association of a TTTA A2 allele on hip BMD. The DOPS study analysis relied on a strategy of comparing "short" and "long" allele lengths (29), which may have obscured differences associated with specific TTTA repeat lengths. In our study, we found the comparison of short and long alleles to be a less powerful analytic method than determining the precise allele that was associated with a difference in bone phenotypes.
The differences in DXA and pQCT BMD, calcaneal QUS, and SSI parameters observed in this study are associated with a lower mean deoxypyridinoline creatinine ratio and a higher FEI. This indicates that a CYP19 A2 allele is associated with decreased bone resorption, perhaps as a result of the higher endogenous estrogen concentration associated with this allele. It has been previously reported that the low endogenous concentrations of estrogen present in late menopausal women are associated with decreased bone resorption and that suppression of this endogenous estrogen production results in increased bone resorption, but not bone formation (10), consistent with our observations.
In this study we demonstrated, in those subjects who had osteoporosis, that the presence of an A2 allele was associated with a greatly reduced risk of vertebral fracture. The CYP19 TTTA microsatellite has been previously reported to be associated with vertebral fracture (16). The increased prevalence of vertebral fracture with an A2 allele was independent of the FEI. Aromatase has direct target organ-specific effects that can result in increased levels of estrogen in target tissues, without substantially affecting circulating estrogen concentrations (24). Such effects would therefore not be assessed by measurement of the FEI so it is possible that the TTTA microsatellite polymorphism is associated with other factors affecting aromatase activity not accounted for by the FEI measurement. We observed that the FEI was associated with between 35 and 39% of the difference observed in the QUS parameters, between 42 and 67% of the difference for hip BMD and accounted for all of the difference in lumbar spine BMD, between the A2 allele present and absent groups. It is important to note that circulating estradiol concentration was measured in a sensitive and reproducible assay. Nevertheless, the circulating concentrations in the two CYP19 TTTA microsatellite groups were not different. An association with circulating estrogen was only evident after adjustment for SHBG, which was slightly but not significantly lower in the A2 group. We are aware that these data do not completely support the concept that an A2 allele has its effects on bone as a result of effects on circulating bioactive estrogen.
The CYP19 TTTA repeat is present in an intronic region that is not associated with gene regulation or with posttranslational expression. Therefore, it is likely that the association we observed is the result of the CYP19 TTTA microsatellite polymorphism being in linkage disequilibrium with a potentially causative polymorphism. This could also account for the reversed association between TTTA allele lengths and bone parameters observed between this and other studies. It is not clear why the association we observed with higher bone mass was with the TTTAn = 7 allele and not also with the TTTAn = 7(3) allele. This suggests that the 3 base deletion/insertion polymorphism that results in the TTTAn = 7(3) allele is associated with a putative causative polymorphic site independently of the TTTAn = 7 allele. A previous association of the insertion/deletion polymorphism with premenopausal breast cancer risk has been published, indicating that this is a possibility (28). Further studies are therefore required to ascertain potentially functional polymorphisms in the CYP19 gene and to validate their biological function in regulating estrogen synthesis and bone structure. An example of such a polymorphsim is the C-T substitution at base pair 1558 in exon 10 (26), which has been reported to be in linkage disequilibrium with the TTTA microsatellite polymorphism (15).
In conclusion, a significant association of the CYP19 TTTA microsatellite with BMD, bone strength, and vertebral fracture was observed in elderly women. This may be associated with endogenous estrogen concentration and subsequent effects on bone turnover. Given that associations between polymorphisms in the CYP19 gene and bone-related phenotypes have now been observed in a number of studies, further investigation to elucidate the functional polymorphisms in the CYP19 gene and their mechanisms of action should be pursued.
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. L. Cote, W. Yoo, A. S. Wenzlaff, G. M. Prysak, S. K. Santer, G. B. Claeys, A. L. Van Dyke, S. J. Land, and A. G. Schwartz Tobacco and estrogen metabolic polymorphisms and risk of non-small cell lung cancer in women Carcinogenesis, April 1, 2009; 30(4): 626 - 635. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. W. Bell, B. W. Brannigan, K. Matsuo, D. M. Finkelstein, R. Sordella, J. Settleman, T. Mitsudomi, and D. A. Haber Increased Prevalence of EGFR-Mutant Lung Cancer in Women and in East Asian Populations: Analysis of Estrogen-Related Polymorphisms Clin. Cancer Res., July 1, 2008; 14(13): 4079 - 4084. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Czajka-Oraniec, W. Zgliczynski, A. Kurylowicz, M. Mikula, and J. Ostrowski Association between gynecomastia and aromatase (CYP19) polymorphisms Eur. J. Endocrinol., May 1, 2008; 158(5): 721 - 727. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. H. Olson, E. V. Bandera, and I. Orlow Variants in Estrogen Biosynthesis Genes, Sex Steroid Hormone Levels, and Endometrial Cancer: A HuGE Review Am. J. Epidemiol., February 1, 2007; 165(3): 235 - 245. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ferrari Single Gene Mutations and Variations Affecting Bone Turnover and Strength: a Selective 2006 Update IBMS BoneKEy, December 1, 2006; 3(12): 11 - 29. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |