|
|
||||||||
CLINICAL STUDIES |
1 Clinical and Molecular Osteoporosis Research Unit, Department of Clinical Sciences Malmö, Lund University, 205 02 Malmö, Sweden2 Department of Orthopaedics, Malmö University Hospital, 205 02 Malmö, Sweden3 Department of Orthopaedics, Karolinska University Hospital, S-171 76 Stockholm, Sweden4 Medical Genetics Unit, Department of Clinical Sciences Malmö, Lund University, 205 20 Malmö, Sweden
(Correspondence should be addressed to K Åkesson; Email: kristina.akesson{at}med.lu.se)
| Abstract |
|---|
|
|
|---|
Design: Population-based association study.
Methods: Four BMP2 SNPs studied. Total body fat and lean mass measured by DEXA in two cohorts: PEAK-25 women aged 25 (±0.00) (n=993) and osteoporosis prospective risk assessment (OPRA) women aged 75 (±0.00) years (n=1001).
Results: We found no association between BMP2 SNPs and fat or lean mass, however, we observed consistent although non-significant trends. Polymorphisms, rs235767 and Ser37Ala, exerted opposing effects on most parameters of soft tissue and bone mass in both cohorts. This relationship appeared to be age specific with large differences between alleles observed (fat mass; Ser37Ala: 14.3% (PEAK-25), –3.5% (OPRA)). These initial results appear to suggest that alleles exerting a beneficial effect in young women may subsequently contribute to phenotypes associated with osteoporosis risk in elderly women.
Conclusions: While further analyses in other comparative populations are necessary, in this study of almost 2000 women we observed interesting, although non-significant trends, regarding the effects of variation in the BMP2 gene on parameters of body mass. Although the exact nature of the relationship remains uncertain, we suggest that the mechanisms are influenced by age and environmental factors.
| Introduction |
|---|
|
|
|---|
Both the fat and lean mass components of soft tissue influence bone mass (3). At its simplest level, the association between body mass and BMD may be mediated through the effects of skeletal loading and via the action of oestrogen, which is metabolised in peripheral adipose tissue (4). More recently, the complexity of the relationship between fat and bone has been recognised. Osteoblasts and adipocytes share a common mesenchymal stem cell lineage with a degree of plasticity between them (5), depending on environmental and physiological signals (6, 7, 8). Evidence suggests both shared and individual genes contribute to the overlapping phenotypes of BMD and body weight (9). One such pleiotropic candidate is the bone morphogenetic protein-2 (BMP2), which is implicated in both osteoblastogenesis and adipogenesis from common progenitor cells. The pathway favoured depends on complex interactions including BMP2 concentration and several other regulators (10). Evidence from in vitro studies also suggests that muscle satellite cells are able to differentiate into osteoblasts following BMP2 exposure (11).
BMP2 gene variants have variously been linked to BMD and osteoporotic fracture, and in our cohorts with ultrasound parameters (12). Since adipocytes may be influenced by BMP2, we hypothesised that variation in BMP2 might contribute to the distribution of fat mass and hence indirectly contribute to bone mass. To investigate the potential relationship between BMP2 gene variants and body composition, we used two population-based Swedish cohorts, the PEAK-25 study of young women and the OPRA study of elderly women.
| Materials and methods |
|---|
|
|
|---|
The study subjects were Caucasian participants in one of two population-based studies. The Malmö Osteoporosis Prospective Risk Assessment (OPRA) study consists of women all aged 75 years, collected between 1995 and 1999. Women were invited to follow-up visits 1, 3 and 5 years after inclusion and we report here on measurements made at the 5-year follow-up visit.
The PEAK-25 cohort consists of young women living in Malmö, Sweden, between 1999 and 2003. The women, all of whom were aged 25 years at inclusion, were randomly selected from the city files. Full details of both cohorts have been reported previously (12, 13). All participants gave informed consent and the studies were approved by the Lund University Ethics Committee.
Measurements
Total body scans from dual-energy X-ray absorptiometry (PEAK-25: Lunar Prodigy, OPRA: Lunar DPX-L (Lunar Corporation, Madison, WI, USA)) were used to analyse body composition; fat mass (g), lean mass (g) and BMD (g/cm2). Lumbar spine (LS) and femoral neck (FN) BMD were also measured. The participants included in the reported analysis comprised n=1001 (OPRA) and n=993 (PEAK-25), for whom a DXA measurement and BMP2 genotype result was available.
BMP2 SNPs and genotyping
In addition to the potentially functional Ser37Ala SNP (rs2273073) originally implicated in low BMD, we selected an additional three SNPs to encompass the gene and its surroundings (http://www.ensembl.org). In the absence of available data on tagging SNPs at that time, SNPs were selected on the following criteria: relative position to the gene, minor allele frequency in Caucasian populations and distance between the polymorphisms. These were (in chromosomal order) rs235710 (10 kb upstream), rs2273073 (Ser37Ala) in exon 2, rs235767 in intron 2 and rs235754 (8 kb downstream) (12).
The polymorphisms were all genotyped blind by the 5' nuclease Taqman assay (Applied Biosystems, Foster City, California, USA) with 1% routinely repeated to ensure agreement of results. A further 3% of all PEAK-25 samples and 4.5% of all OPRA samples underwent additional checks for the Ser37Ala polymorphism, with concurrence of results. Full details of the methodology have been described previously (12).
Statistical analysis
Statistical analyses were performed using SPSS v14 (SPSS Inc., Chicago). Hardy–Weinberg equilibrium (HWE) was calculated by
2-test. Linkage disequilibrium between SNPs (Lewontin's D' and r2) was calculated using the EH, 2BY2 (14) and Haploview programs. Haplotypes were derived from Ser37Ala and rs235767 using PHASE v2.02 (15). Two common haplotypes accounted for >98% of alleles at the BMP2 locus: haplotype 1 Ser/A (61.5%) and haplotype 2 Ser/C (36.5%). Haplotype analyses were performed, according to whether individuals possessed zero, one or two copies of the haplotype.
Non-normally distributed variables were log transformed prior to analysis (<0.95, Shapiro–Wilk). Genotype-specific differences in variables were investigated using GLM-ANOVA, correcting for confounders including height, oestrogen use, years since menarche and smoking, but not age, since all individuals within each cohort were of the same age (PEAK-25: 25.5 years, range 25.03–25.99 and OPRA: 75.2 years, range 75.0–75.8). Multiple regression analysis determined independent predictors of fat and lean mass, using a model containing genotype or haplotype, height, smoking, age at menarche and contraceptive pill use as variables. The P values are presented uncorrected. Significance was not retained when the adjusted level of significance equivalent to P=0.05, based on a modified Bonferroni correction, was applied.
| Results |
|---|
|
|
|---|
2=5.9; P=0.05). Neither cohort displayed any genotype-related differences in the general clinical characteristics.
|
8%) and less fat (25%) tissue mass, which is reflected in the lower ratio of fat to lean mass. Smoking incidence was higher in the young women (26% vs 13.8%) and the age of onset of menarche was younger (12.7 vs 14.0 years). There were no genotype-related differences in height, weight or BMI for either cohort.
|
In both the PEAK-25 and OPRA cohorts, there was no significant association with fat mass for either genotypes or haplotypes, neither could we detect any differences in the fat–lean mass ratio. Categorising individuals according to whether they fell into the lowest or highest quintile of BMI did not appreciably alter the results.
In the young women, regression analysis identified rs235767 as a predictor of lean mass (P=0.008) although this decreased to P=0.06 when height was additionally included in the model. Although individuals carrying one or two copies of the variant C allele had incrementally lower values, the association was not significant. Individuals carrying two copies of the common haplotype 1 (Ser-A) had the highest lean mass values when compared with one copy (intermediate) and zero copies (lowest), P=0.04. As before, when height was taken into consideration, the relationship was abolished. The data did not alter appreciably when BMI category was taken into consideration. Haplotype 2 (Ser-C) was a non-significant predictor of lean mass (P=0.081) and values decreased with haplotype copy number.
In the elderly women, BMP2 genotype was not an independent predictor of lean mass at baseline although a borderline association was identified between rs235767 and lean mass (P=0.052) at follow-up. Individuals who were homozygous for the variant allele had higher lean mass at both time points but this was not incremental and was not significant either before or after correction. There were no significant genotype differences for the ratios of fat to lean mass or weight change between 20 and 75 years. Analysing only those in the highest and lowest quintiles of BMI did not appreciably alter the results. Exclusion of 69 women with reported diabetes (based on the known effects on adipose tissue from diabetes) gave essentially the same results. Regression analysis identified haplotype 1 as an independent predictor of lean mass at follow-up (P=0.037) but not baseline. Similarly, haplotype 2 had a non-significant trend (P=0.08). Although non-significant, individuals carrying one or two copies of haplotype 1 had the lowest values but not in a dose-dependent manner, while individuals carrying one or two copies of haplotype 2 had the highest values and a gene dose effect was observed.
SNP and age effects on body composition
SNP rs235767 showed consistent dose dependent, although non-significant trends in the direction of effect across the parameters measured. For BMI, weight, lean mass, fat mass and BMD at all sites, CC homozygotes had the lowest values for all measures (Fig. 1A PEAK-25). SNP Ser37Ala also showed consistent non-significant trends of effect, however, in this instance the Ala/Ala homozygotes, tended to display the highest mean values for the measures of body composition and BMD at the spine (Fig. 1A). (For the other SNPs, the relationship between genotype and phenotype was more variable (data not shown)).
|
As in the young women, rs235767 showed consistently similar, non-significant trends in the direction of effect across the parameters measured at baseline. It is interesting to note, however, that the direction of the effect is the opposite of what is seen in the younger women, i.e. the homozygous variant individuals tended to have the highest mean values. With the exception of lean mass, the relationship was dose dependent (Fig. 1B). As far as can be compared, given the absence of homozygotes for the rare Ala allele in the OPRA cohort, the Ser37Ala trends, as observed in the young women, were the opposite of those for rs235767 (Fig. 1B OPRA). Individuals with two copies of haplotype 1 tended to have the lowest values and individuals with two copies of haplotype 2 tended to have the highest values (data not shown).
Although it is not possible to make direct comparisons between the young and older women, the results suggest that in young adulthood, SNP Ser37Ala seems to be associated with higher fat and lean mass levels (for heterozygotes,
7% greater fat and 3% greater lean mass compared with Ser/Ser homozygotes) and higher BMD levels. In old age, this relationship is reversed and the variant allele is associated with reduced fitness in terms of these measurements (Fig. 2A and B).
|
5 years later. It is apparent that in very elderly women, rs235767 CC homozygotes and Ser37Ala heterozygotes loose more fat mass than lean mass during this period compared with the other genotypes.
|
| Discussion |
|---|
|
|
|---|
Several genes that have been studied in relation to BMD, e.g. VDR, ER and LRP5, have pleiotropic effects and are implicated in a range of diseases including obesity and insulin resistance (17). The BMP2 gene has received some attention with regard to BMD (18, 19), however, despite a role in adipocyte differentiation both in vitro and in vivo (10) it has not yet been studied with regard to body composition. In the present study, we have investigated the relationship between variation in the BMP2 gene and the components of body composition: fat and lean mass in 25- and 75-year-old women.
Our results indicate that while BMP2 genotype is not significantly associated with absolute measures of fat mass, lean mass (or indeed bone mass), there is nevertheless a possible, albeit complex relationship between them which appears to be altered by age and possibly by other genetic and environmental interactions.
Two polymorphisms, rs235767 and the potentially functional Ser37Ala, exerted opposing effects on almost all the parameters of soft tissue and bone mass measured in both cohorts of women. We also observed that alleles that could be considered beneficial to bone mass in the young women subsequently contributed to phenotypes associated with osteoporosis risk in the elderly women. Opposing effects of polymorphisms, such as we have observed, have previously been reported for PPAR-
(20) and MTHFR genotypes (21); however, we are unaware of other studies that have reported, as we do, age-related reversal of the relationship between genotype and phenotype. We cannot fully explain our observations in the context of the present study. The results could have arisen simply by chance – replication in comparative populations is necessary to refute or confirm our findings. In the absence of detailed dietary information in these cohorts, we also cannot exclude the possibility that small differences in BMI could have been contributed to by differences in the diet.
It could, however, reflect age-related environmental exposures, given that there are 55–60 years between the childhood and young adolescence of the two cohorts, or gene–microenvironment interactions. An alternative explanation could be antagonistic pleiotropy (whereby genes, advantageous during times of reproductive activity have harmful effects later in life), which is thought to be uncommon in humans, but has been observed in animal models of reproduction, immunity and body composition (22). The literature provides evidence to suggest that our findings are not the result of statistical anomaly. BMPs control stem cell commitment to various lineages and in vitro studies indicate that adipocyte differentiation depends on BMP2 dosage, precursor cell type and cell stage (23, 24) and the presence of other regulators such as Schnurri-2 and PPAR-
(10, 25, 26). The balance between osteogenic and adipogenic effects mediated by BMP2 appear to be reciprocal and influenced by ageing, such that ageing activates adipogenic and suppresses osteogenic characters in stem cells (27, 28, 29) an observation that may contribute indirectly to the explanation of our findings.
Twin studies indicate that the similarity between monozygotic twins diminish with age, suggesting that gene penetrance is susceptible to external influences (30, 31). By studying two cohorts from opposite ends of the age spectrum, we have identified relationships that might have otherwise been overlooked as a methodological artefact between individually published studies and therein lies one of the strengths of our study. The narrow age range and similar ethnic background of our cohorts also eliminate at least some of the confounding variation. Our data also highlight the importance of considering the effect of gene variation, not as static, but in the context of constantly changing interactions with the microenvironment. A weakness in our study is that at the time when genotyping was completed, no functional SNPs apart from Ser37Ala had been identified. Despite having chosen our SNPs to broadly encompass the gene region and inclusion of the potentially functional Ser37Ala polymorphism, it is not possible to capture all variation within the gene. However, 76% of the genetic variation in the locus is represented by five possible haplotypes identified in our population and inclusion of tag SNPs from the conserved haplotype block covering the gene, would have captured only an additional 15%. Until replication studies have been performed in comparative age-distributed populations; however, we should be cautious in the conclusions we draw from these potentially interesting observations.
To conclude, in this study consisting of two distinct populations of almost 2000 women, we suggest that polymorphisms of the BMP2 gene may be linked with measures of body composition, fat and lean mass. Interestingly, the effects are not uniform, indicating that age and environmental factors may play important roles through mechanisms that are not possible to fully explore in the current setting.
| Acknowledgements |
|---|
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |