Common Genetic Variants Are Associated with Accelerated Bone Mineral Density Loss after Hematopoietic Cell Transplantation

Abstract

Background

Bone mineral density (BMD) loss commonly occurs after hematopoietic cell transplantation (HCT). Hypothesizing that genetic variants may influence post-HCT BMD loss, we conducted a prospective study to examine the associations of single nucleotide polymorphisms (SNP) in bone metabolism pathways and acute BMD loss after HCT.

Methods and Findings

We genotyped 122 SNPs in 45 genes in bone metabolism pathways among 121 autologous and allogeneic HCT patients. BMD changes from pre-HCT to day +100 post-HCT were analyzed in relation to these SNPs in linear regression models. After controlling for clinical risk factors, we identified 16 SNPs associated with spinal or femoral BMD loss following HCT, three of which have been previously implicated in genome-wide association studies of bone phenotypes, including rs2075555 in COL1A1, rs9594738 in RANKL, and rs4870044 in ESR1. When multiple SNPs were considered simultaneously, they explained 5–35% of the variance in post-HCT BMD loss. There was a significant trend between the number of risk alleles and the magnitude of BMD loss, with patients carrying the most risk alleles having the greatest loss.

Conclusion

Our data provide the first evidence that common genetic variants play an important role in BMD loss among HCT patients similar to age-related BMD loss in the general population. This infers that the mechanism for post-HCT bone loss is a normal aging process that is accelerated during HCT. A limitation of our study comes from its small patient population; hence future larger studies are warranted to validate our findings.

About the Authors:

Song Yao

Affiliation: Department of Cancer Prevention and Control, Roswell Park Cancer Institute, Buffalo, New York, United States of America

Lara E. Sucheston

Affiliation: Department of Cancer Prevention and Control, Roswell Park Cancer Institute, Buffalo, New York, United States of America

Shannon L. Smiley

Affiliation: Department of Medicine, Roswell Park Cancer Institute, Buffalo, New York, United States of America

Warren Davis

Affiliation: Department of Cancer Prevention and Control, Roswell Park Cancer Institute, Buffalo, New York, United States of America

Jeffrey M. Conroy

Affiliation: Department of Cancer Genetics, Roswell Park Cancer Institute, Buffalo, New York, United States of America

Norma J. Nowak

Affiliations Department of Cancer Genetics, Roswell Park Cancer Institute, Buffalo, New York, United States of America, Department of Biochemistry, University at Buffalo, Buffalo, New York, United States of America

Christine B. Ambrosone

Affiliation: Department of Cancer Prevention and Control, Roswell Park Cancer Institute, Buffalo, New York, United States of America

Philip L. McCarthy Jr.

Contributed equally to this work with: Philip L. McCarthy Jr., Theresa Hahn

Affiliation: Department of Medicine, Roswell Park Cancer Institute, Buffalo, New York, United States of America

Theresa Hahn

Contributed equally to this work with: Philip L. McCarthy Jr., Theresa Hahn

* E-mail: [email protected]

Affiliation: Department of Medicine, Roswell Park Cancer Institute, Buffalo, New York, United States of America

Introduction

After the second decade of life, bone mineral density (BMD) gradually declines as a result of normal aging [1]. This process is accelerated in patients undergoing hematopoietic cell transplantation (HCT) [2], [3]. Rapid bone loss weakens bone strength and may permanently impair bone remodeling capability [4], [5], putting HCT survivors at high risk of premature fracture. We previously demonstrated that BMD loss within 4 months after HCT is equivalent to 7 to 17 years of normal bone aging [6]. Moreover, this loss occurs following autologous and allogeneic HCT with a similar incidence and severity.

In the general population, the inherited genetic background accounts for a large proportion of phenotypic variation in BMD loss and fracture risk [7], [8]. In numerous candidate gene association studies, a variety of genetic variants, mostly in the form of single nucleotide polymorphisms (SNP), have been associated with regulation of BMD and risk of fracture [9], with less than a handful of them confirmed in a recent large-scale meta-analysis performed by the Genetic Factors for Osteoporosis (GEFOS) Consortium [10]. A genome-wide association (GWA) approach has also been used to search for loci associated with various bone phenotypes, which identified a number of novel SNPs [11], [12], [13], [14], [15]. These genetic markers may have clinical significance in predicting risk of BMD loss and/or osteoporosis. Moreover, some may also be predictive for response to therapies targeting related molecules or pathways, such as a human monoclonal antibody against receptor activator of nuclear factor κB ligand (RANKL) (denosumab), which was recently approved by the FDA for postmenopausal women at risk for fracture [16].

Among HCT recipients, rapid BMD loss may be due to an accelerated bone aging process after exposure to high doses of corticosteroids, chemotherapy and/or total body irradiation. Indeed, we have previously demonstrated that bone loss occurs within three months of HCT in both allogeneic and autologous HCT patients [6]. Therefore, predisposing genetic variants for low BMD and osteoporosis in the general population may also predispose carriers to risk of severe post-HCT BMD loss. Recently, McClune and colleagues proposed comprehensive guidelines for the screening, prevention and treatment of osteoporosis after HCT. They recommend dual-energy X-ray absorptiometry (DXA) scans for all adult patients 1 year post-HCT or earlier for patients at high risk of bone loss [17]. Bisphosphonates are effective in reducing or reversing HCT-related bone loss in several small randomized studies; however, they rarely resolve bone loss to pre-HCT levels and have a significant side-effect, osteonecrosis of the jaw. We hypothesized that genetic variants are associated with BMD loss in HCT patients, and sought to identify a panel of genetic biomarkers that can be incorporated into the above guidelines to personalize prevention and treatment strategy.

Methods

Ethics statement: Stem cell transplant patients at the Blood and Marrow Transplantation Program included in this study are part of an institution-level Data Bank and Biorepository (DBBR) at Roswell Park Cancer Institute (RPCI). All participants have provided written consent for their data and biospecimens to be used for research purposes. This study has been reviewed and approved by RPCI Institutional Review Board.

Patients, BMD measurement, and DNA specimens

A detailed description of the patient population has been published previously [6]. Briefly, between 1/2006 and 1/2009, 206 consecutive adult (≥18 years) patients underwent their first HCT in the Blood and Marrow Transplantation Program at Roswell Park Cancer Institute (RPCI). As part of routine clinical care, the average BMD at lumbar spine (L2–L4) and dual femurs were measured by DXA scans at baseline (pre-HCT) and day +100 post-HCT. The median time for 197 patients (96%) who had a baseline DXA scan was 20 days before HCT, and 98 days for 146 patients (74%) who had a second DXA scan after HCT. Reasons for not obtaining a follow-up DXA scan were unstable medical status (n = 14), or disease relapse/early death (n = 37). Before HCT, patients were invited to participate in the DataBank and BioRepository (DBBR), a comprehensive data and biospecimen core facility at RPCI [18]. Upon their consent to the DBBR, blood samples were drawn and genomic DNA was extracted and stored for multidisplinary research. For this study, 121 patients including 67 autologous and 54 allogeneic HCT patients, consented to the DBBR and had adequate DNA available for genotyping.

Candidate gene and SNP selection and genotyping

Candidate genes from bone metabolism pathways were selected based on the following criteria: (1) genes that have been well studied and established in association with osteoporosis or other bone phenotypes in the general population; (2) key genes in bone metabolism pathways, including osteoblast or osteoclast differentiation, proliferation and activation; calcitronic hormone metabolism and activation; and bone matrix formation and regulation. As listed in Table S1, a total of 45 candidate genes were clustered into 4 groups based on underlying biological function, including (1) RANKL-RANK-OPG central signaling axis and regulating cytokines, (2) bone matrix proteins and regulating factors, (3) vitamin D receptor and metabolizing enzymes, and (4) steroid hormones and receptors.

SNPs for each candidate gene were selected based on the following criteria: (1) minor allele frequency (MAF) ≥0.10 in the population of European ancestry; (2) SNPs that have been well studied and established in association with osteoporosis or other bone traits in the general population, including those identified in recent GWA studies and the large prospective study Genetic Factors for Osteoporosis (GEFOS) Consortium [10]; (3) coding SNPs and SNPs in the 3′ or 5′ untranslated region (UTR); and (4) tagSNPs selected by the Tagger program [19] from the HapMap Project data [20], if no functional SNPs or well studied SNPs were available. In all, 158 SNPs were initially selected from the 45 candidate genes. To remove redundant SNPs due to high linkage disequilibrium (LD), pair-wise r2 was computed for SNPs in each gene and in genes clustered on the same chromosome by PLINK [21]; SNP_tools [22] was used for necessary data formatting. For SNP pairs in high LD (r2≥0.90), one SNP was randomly selected from each pair for further consideration. As a result, a total of 122 SNPs were included in the final analysis (see Table S1).

Genotyping of DNA extracted from peripheral blood was performed in the Genomics Shared Resource at RPCI utilizing the MassARRAY® technology and iPLEX Gold assay (Sequenom). Nine percent duplicate samples and in-house Coriell trio samples were included within plates and across plates for genotyping quality assurance. The average qualified call rate was 98.1% for each SNP and 98.3% for each DNA sample. There was no discordance among duplicate pairs, and no SNPs violated Mendelian inheritance or Hardy-Weinberg equilibrium.

Statistical analysis

All analyses were performed stratified by autologous and allogeneic HCTs since these two treatments, and the characteristics of the patients given these two treatments, are significantly different. Change in BMD was standardized to 100 days (BMD change divided by the number of days between the two DXA scans then times 100) and used as a continuous variable; thus, a negative value indicates BMD loss and a positive value indicates BMD gain. We used linear regression models to analyze the association of clinical factors and genetic variants with four separate outcomes: BMD changes at the spine and at the femur in autologous and allogeneic HCT patients. Assuming an additive genetic model, 54 patients in the small allogeneic HCT group and two-sided p-value of 0.05, we had 80% power to detect a minimum BMD change of 0.023–0.040 g/cm2 for a SNP with a minor allele frequency between 0.10–0.50. In the autologous HCT group which had a larger sample size of 67, we had 80% power to detect effect sizes <0.023 g/cm2.

Single SNP analysis.

SNPs associated with HCT-induced BMD change were identified as follows. First, each SNP was tested individually in a simple linear regression model without controlling for covariates; SNPs with p<0.10 were then added to regression models previously built with clinical factors for each of the four outcomes, including the cumulative dose of corticosteroids (in prednisone equivalents) between the two DXA scans for BMD loss in allogeneic HCT patients [6]. We measured the proportion of variance explained by all the covariates retained in a model, R2, and determined, ΔR2, an estimate of the proportion of variance that can be explained by the SNP(s) added to the base model, with adjustment for the number of covariates in the model. A SNP with a ΔR2 of 5% or higher was deemed significant. Effect size (i.e., absolute BMD change) per copy of the variant allele and 95% confidence intervals (CI) were calculated for each significant SNP, with and without adjustment for clinical risk factors. The rate of annualized BMD change between pre- and post-HCT DXA scans was compared and plotted by the genotypes of each significant SNP. To explore potential functionality of the significant SNPs, the SNP and CNV Annotation Database (SCAN) [23] was surveyed to explore whether it is a potential expression quantitative trait loci (eQTL) that has been associated with altered gene expression.

Multiple-SNP analysis.

We examined the combined effects of multiple variants using backward elimination. SNPs with a p≤0.05 were retained in the final models. Adjusted R2 was calculated for the complete models with both clinical and genetic factors and compared to the base models with only clinical factors, to determine the overall contribution of genetic factors. For SNPs retained in the model, the number of risk alleles was summed for each patient, and the annualized rate of BMD change was compared across patients by number of risk alleles. Although we measured change in BMD between a set interval of approximately 100 days, we presented annualized BMD loss rates for comparisons to the general population. Prior reports of BMD loss post-HCT demonstrated that reversal of BMD loss takes several years after HCT [3]. Multiple-SNP analysis was not performed for spinal BMD loss after allogeneic HCT because only one SNP was retained in the final comprehensive model. All analyses were two-tailed and were performed with SAS 9.2 (SAS Institute, Cary, NC).

Results

Individual SNP analysis

Table 1 summarizes the demographic and clinical characteristics of the study population separately by auto- and allo-HCT groups, which were statistically similar to the overall patient population described in our prior report [6]. Unadjusted p-values for associations between individual SNPs and BMD change between pre- and post-HCT are shown in Figure S1. After controlling for clinical risk factors, at least one SNP remained significant for each of the four outcomes, with a total of 16 significant SNPs for all four outcomes. SNP rs2075555 in the COL1A1 gene remained significant after Bonferroni correction for multiple comparisons (p = 0.037).

[Figure omitted. See PDF.]

Table 1. Characteristics of HCT patients with BMD measurement and genomic DNA samples.

https://doi.org/10.1371/journal.pone.0025940.t001

Table 2 summarizes the unadjusted and adjusted BMD change associated with each of the 16 SNPs significant in individual SNP tests. Annualized rate of BMD change by genotype are shown in Figure 1. There was a significant dose-response trend between rate of BMD loss and number of risk alleles at the spine and the femur in autologous HCT and at the femur in allogeneic HCT. The functions of these genes, SNPs and published literature are summarized in Table S2. Among these, three SNPs, including rs2075555 in COL1A1 gene [11], [24], rs9594738 in RANKL gene [13] and rs4870044 in ESR1 gene [13], were significantly associated with bone phenotypes in GWA studies of the general population. Moreover, rs2075555 and rs9594738 were associated with altered gene expression in the SCAN (p<0.0001). Another 9 significant SNPs in our analysis were previously linked to bone phenotypes in candidate gene association studies, including rs4588 in GC, rs6256 in PTH, rs419598 in IL1RN, rs1801131 in MTHFR, rs1800896 in IL10, and rs1061624 in TNFRSF1B, rs1042358 in ALOX12, rs759330 in BGLAP, and rs2235579 in CLCN7. Four of these SNPs (rs759330, rs2235579, rs1800896 and rs1061624) were also associated with altered gene expression in the SCAN. The other 4 SNPs, including rs3787557 and rs2296241 in CYP24A1, rs25645 in CSF3, and rs1321080 in RUNX2, were selected as tagSNPs or SNPs in functional genomic regions, with no known associations with bone phenotypes. Two of these SNPs (rs2296241 and rs25645) were associated with altered gene expression in the SCAN.

[Figure omitted. See PDF.]

Figure 1. Annualized rate of bone mineral density change after hematopoietic cell transplantation by genotype for SNPs significant in the multivariable models including clinical risk factors.

Mean and standard error of the annualized rate of BMD change is plotted by genotype for significant SNPs in the multivariable models including clinical risk factors. If the frequency of the homozygous variant genotype was ≤5%, it was combined with the heterozygous genotype (includes rs3787557 [CYP24A1], rs6265 [PTH], rs1801131 [MTHFR] and rs2075555 [COL1A1]). Please refer to Table S2 for a description of functions of these genes and SNPs. Abbreviation: BMD, bone mineral density; HCT, hematopoietic cell transplantation; SNP, single nucleotide polymorphism.

https://doi.org/10.1371/journal.pone.0025940.g001

[Figure omitted. See PDF.]

Table 2. Change in bone mineral density per copy of the variant allele for individual SNPs which were significant in the multivariable models including clinical risk factors1.

https://doi.org/10.1371/journal.pone.0025940.t002

Multiple-SNP analysis

The multiple-SNP models are shown in Table 3. For the spinal BMD change after autologous HCT, three SNPs (rs759330 in BGLAP, rs4588 in GC, and rs3787557 in CYP24A1) remained in the final model and explained 16% of the variance in BMD loss, in addition to the 10% determined by the clinical risk factors. Similarly, SNPs retained in the final models explained an additional 35%, 5% and 29% of the variance in femoral BMD loss after autologous HCT, spinal BMD loss after allogeneic HCT, and femoral BMD loss after allogeneic HCT, respectively.

[Figure omitted. See PDF.]

Table 3. Multiple SNP models of bone mineral density change after hematopoietic cell transplantation.

https://doi.org/10.1371/journal.pone.0025940.t003

When the annualized rate of BMD loss was compared by the number of risk alleles, there was a significant association between higher number of risk alleles and higher BMD loss (Figure 2). In addition, the multiple-SNP models identified patients at the highest risk of BMD loss. Patients carrying the highest number of risk alleles had the most severe BMD loss.

[Figure omitted. See PDF.]

Figure 2. Annualized rate of bone mineral density change after hematopoietic cell transplantation by the number of risk alleles of the SNPs significant in the multiple SNP models.

Mean and standard error of annualized rate of BMD change is plotted by the number of risk alleles in significant SNPs from the multiple SNP models as shown in Table 3. The p-values for a trend between the number of risk alleles and BMD change were ≤0.003 for all three outcomes. The proportion of HCT patients carrying the specified number of risk alleles is shown. The multiple SNP analysis was not performed for spinal BMD change after allogeneic HCT since only one SNP was significant in the multiple SNP model. Please refer to Table S2 for a description of these genes and SNPs. Abbreviation: BMD, bone mineral density; HCT, hematopoietic cell transplantation; SNP.

https://doi.org/10.1371/journal.pone.0025940.g002

Discussion

Our analysis identified a number of SNPs in bone metabolism pathways predictive for acute BMD loss following HCT. Because our candidate genes and SNPs have been previously implicated in age-related bone loss, osteoporosis and/or bone phenotypes in the general population, our findings support the hypothesis that HCT-related BMD loss is an acceleration of the normal bone aging process. Moreover, SNPs showed greater effects in HCT patients in our study than those identified in the general population for bone phenotypes. Four SNPs accounted for 35% of variance in femoral BMD change after autologous HCT, whereas the 12 most significant loci in a GWA study accounted for only about 3% of variance in hip and spine BMD in a healthy Icelandic population [13]. The stronger effect may be driven by potent exposures during HCT.

Among the 16 SNPs significant in our analysis, only rs9594738 in RANKL gene was related to BMD loss at both the spine and the femur in allogeneic patients. This marginal overlap could be attributed to differences in exposures related to allogeneic versus autologous transplant procedures, as well as to differences in susceptibility to HCT-induced BMD loss at these two anatomic sites. The trabecular interconnection of the lumbar spine is much higher than that of the femur [25]. Thus, the difference in bone microarchitecture may explain the lack of overlap of SNPs associated with BMD and risk of fracture at the two sites in previous studies in the general populations [10] and the observations that femoral bone was more susceptible to HCT-induced bone loss than spinal bone [6], [26].

A limitation of our study comes from the relatively small sample size and the large number of tests performed: 122 SNPs in relation to 4 quantitative outcomes in 121 patients stratified by autologous and allogeneic HCT. The statistical power was thus limited. Instead of using a Bonferroni correction, which would be over conservative, we propose a hierarchy of evidence in our findings based on curated knowledge of the SNPs in the literature.

The first level of evidence includes three SNPs implicated in two GWA studies of bone phenotypes. We found that rs2075555 in the intronic region of COL1A1 is associated with femoral BMD loss after allogeneic HCT. The gene encodes for type I collagen, a major matrix protein for bone formation and this SNP has been previously associated with femoral neck width in women and shaft width in men [11]. Two other SNPs, rs9594738 near RANKL and rs4870044 near ESR1 from a GWA study of BMD in an Icelandic population [13], and a meta-analysis of GWA studies of BMD14, were also significant in our study. RANKL plays a central role in coupling the activities of osteoblasts and osteoclasts; and ESR1 encodes for estrogen receptor alpha, which is key in estrogen-regulated activities.

The second level of evidence includes 9 SNPs, which have been previously studied in relation to various bone phenotypes using a candidate gene approach: rs4588 in GC, rs6256 in PTH, rs419598 in IL1RN, rs1801131 in MTHFR, rs1800896 in IL10, and rs1061624 in TNFRSF1B, and rs1042358 in ALOX12, rs759330 in BGLAP, and rs2235579 in CLCN7 (Table S2). The first 5 SNPs are among the most widely studied candidates for osteoporosis and/or other bone phenotypes.

The third level of evidence includes 4 SNPs which have not been previously studied in relation to bone phenotypes: rs3787557 and rs2296241 in CYP24A1, rs25645 in CSF3, and rs1321080 in RUNX2. These SNPs are either tagSNPs or located in potentially functional regions, and these genes were selected for their importance in bone metabolism (Table S2).

As discussed above, we relied on a priori knowledge to rank our confidence in our results: those previously implicated in GWA studies were deemed most reliable, followed by those implicated in candidate gene studies, and lastly those without support from the literature. Although this knowledge-based hierarchy of evidence is not statistically stringent, and we cannot exclude the possibility of false positive (or negative) findings in our results, it aligns with our exploratory purpose and provides support that genetic variants contribute to post-HCT bone loss. Moreover, our findings may provide clues for future studies. Thus, in addition to rs2075555 in COL1A1, which remained significant in our analysis under the most stringent threshold of significance, we also reported other SNPs that may be interesting. Our study was designed to evaluate only one post-HCT BMD measurement (about +100 days); thus we do not know how the genetic variants described in our study affect long-term BMD post-HCT. We chose day +100 for post-HCT BMD measurement for three reasons. First, BMD loss occurs at the highest rate within the first 3–6 months after HCT [3], thus a true gene-exposure influence would be strongest in the time period closest to the exposure. Second, patients who were osteoporotic at their day +100 DXA scan were treated with a bisphosphonate soon afterwards which may distort the association between genotypes and BMD loss. Third, as an exploratory study, BMD loss at day +100 was an outcome immediately available, while we continue to follow-up patients with the plan to assess genotypes and long-term BMD loss.

It is of interest to investigate the association of donor genotype with BMD loss in allogeneic HCT recipients, considering their hematopoietic cells, and thus osteoclast precursors, are replaced by donors' hematopoietic cells once complete donor chimerism has been achieved. Unfortunately, donor DNA is not available from the allogeneic HCT patients in this study. However, given that the acute BMD loss we observed occurred quickly following transplantation, prior to complete donor chimerism, and that the acute BMD loss in our study was similar between autologous and allogeneic HCT recipients (ie., in the absence and presence of donor genomes), any effect of donor genotype may only be evidenced in chronic BMD loss with longer follow-up. This question warrants further investigation.

In summary, given the continuous bone loss due to normal aging, accelerated post-HCT bone loss may never recover without aggressive therapy, and HCT patients may be at much higher risk of fracture compared to their age-matched peers. Considering the well established relationship between low BMD and fracture risk in the general population, we support the routine post HCT surveillance guidelines developed by McClune et al [17]. Further studies are needed to determine the impact and timing of different prevention and/or intervention strategies.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Log-transformed p-values for associations between bone mineral density change and individual SNPs without adjustment for clinical risk factors. Individual SNP (gene) remained significant in the multivariable models including clinical risk factors are labeled in circles. A number of other SNPs significant in univariate analysis but dropped out of multivariate analysis were not labeled. Please refer to Table S2 for description of functions of these genes and SNPs. Abbreviation: HCT, hematopoietic cell transplantation; SNP, single nucleotide polymorphism.

https://doi.org/10.1371/journal.pone.0025940.s001

(TIF)

Table S1.

Selected 122 single nucleotide polymorphisms from 46 candidate genes in bone metabolism pathways. 1SNPs that have been previously studied in the large prospective study Genetic Markers for Osteoporosis (GENOMOS) Consortium were labeled as “GENOMOS SNPs”. SNPs that have been identified in genome wide association (GWA) studies were labeled as “GWAS SNP”. TagSNPs selected to represent genetic variations in genes that have not been previous studied and have no common SNPs in functional regions were labeled as “tagSNP”. Functional regions include 3′ and 5′ untranslated region (UTR) and codon. For coding SNPs, amino acid changes are shown.

https://doi.org/10.1371/journal.pone.0025940.s002

(DOC)

Table S2.

Functions and summary of the published literature for genes and SNPs significant in the individual SNP models including clinical risk factors

https://doi.org/10.1371/journal.pone.0025940.s003

(DOC)

Acknowledgments

The authors are grateful to the RPCI patients who agreed to participate in this study, the BMT nursing and laboratory staff, the nurse practitioners, physician assistants and the entire clinical team who cares for these patients, and the BMT data managers for assistance with data collection.

Author Contributions

Conceived and designed the experiments: SY LES SLS CBA PLM TH. Performed the experiments: SY WD JMC NJN. Analyzed the data: SY LES SLS CBA PLM TH. Contributed reagents/materials/analysis tools: WD JMC NJN CBA. Wrote the paper: SY LES SLS WD JMC NJN CBA PLM TH.

Citation: Yao S, Sucheston LE, Smiley SL, Davis W, Conroy JM, Nowak NJ, et al. (2011) Common Genetic Variants Are Associated with Accelerated Bone Mineral Density Loss after Hematopoietic Cell Transplantation. PLoS ONE6(10): e25940. https://doi.org/10.1371/journal.pone.0025940

References

1. Benjamin RM (2010) Bone health: preventing osteoporosis. J Am Diet Assoc 110: 498.RM Benjamin2010Bone health: preventing osteoporosis.J Am Diet Assoc110498

2. Kelly PJ, Atkinson K, Ward RL, Sambrook PN, Biggs JC, et al. (1990) Reduced bone mineral density in men and women with allogeneic bone marrow transplantation. Transplantation 50: 881–883.PJ KellyK. AtkinsonRL WardPN SambrookJC Biggs1990Reduced bone mineral density in men and women with allogeneic bone marrow transplantation.Transplantation50881883

3. Schulte CM, Beelen DW (2004) Bone loss following hematopoietic stem cell transplantation: a long-term follow-up. Blood 103: 3635–3643.CM SchulteDW Beelen2004Bone loss following hematopoietic stem cell transplantation: a long-term follow-up.Blood10336353643

4. Tauchmanova L, Serio B, Del Puente A, Risitano AM, Esposito A, et al. (2002) Long-lasting bone damage detected by dual-energy x-ray absorptiometry, phalangeal osteosonogrammetry, and in vitro growth of marrow stromal cells after allogeneic stem cell transplantation. J Clin Endocrinol Metab 87: 5058–5065.L. TauchmanovaB. SerioA. Del PuenteAM RisitanoA. Esposito2002Long-lasting bone damage detected by dual-energy x-ray absorptiometry, phalangeal osteosonogrammetry, and in vitro growth of marrow stromal cells after allogeneic stem cell transplantation.J Clin Endocrinol Metab8750585065

5. Lee WY, Cho SW, Oh ES, Oh KW, Lee JM, et al. (2002) The effect of bone marrow transplantation on the osteoblastic differentiation of human bone marrow stromal cells. J Clin Endocrinol Metab 87: 329–335.WY LeeSW ChoES OhKW OhJM Lee2002The effect of bone marrow transplantation on the osteoblastic differentiation of human bone marrow stromal cells.J Clin Endocrinol Metab87329335

6. Yao S, Smiley SL, West K, Lamonica D, Battiwalla M, et al. (2010) Accelerated bone mineral density loss occurs with similar incidence and severity, but with different risk factors, after autologous versus allogeneic hematopoietic cell transplantation. Biol Blood Marrow Transplant 16: 1130–1137.S. YaoSL SmileyK. WestD. LamonicaM. Battiwalla2010Accelerated bone mineral density loss occurs with similar incidence and severity, but with different risk factors, after autologous versus allogeneic hematopoietic cell transplantation.Biol Blood Marrow Transplant1611301137

7. Peacock M, Turner CH, Econs MJ, Foroud T (2002) Genetics of osteoporosis. Endocr Rev 23: 303–326.M. PeacockCH TurnerMJ EconsT. Foroud2002Genetics of osteoporosis.Endocr Rev23303326

8. Eisman JA (1999) Genetics of osteoporosis. Endocr Rev 20: 788–804.JA Eisman1999Genetics of osteoporosis.Endocr Rev20788804

9. Liu YZ, Liu YJ, Recker RR, Deng HW (2003) Molecular studies of identification of genes for osteoporosis: the 2002 update. J Endocrinol 177: 147–196.YZ LiuYJ LiuRR ReckerHW Deng2003Molecular studies of identification of genes for osteoporosis: the 2002 update.J Endocrinol177147196

10. Richards JB, Kavvoura FK, Rivadeneira F, Styrkarsdottir U, Estrada K, et al. (2009) Collaborative meta-analysis: associations of 150 candidate genes with osteoporosis and osteoporotic fracture. Ann Intern Med 151: 528–537.JB RichardsFK KavvouraF. RivadeneiraU. StyrkarsdottirK. Estrada2009Collaborative meta-analysis: associations of 150 candidate genes with osteoporosis and osteoporotic fracture.Ann Intern Med151528537

11. Kiel DP, Demissie S, Dupuis J, Lunetta KL, Murabito JM, et al. (2007) Genome-wide association with bone mass and geometry in the Framingham Heart Study. BMC Med Genet 8: Suppl 1S14.DP KielS. DemissieJ. DupuisKL LunettaJM Murabito2007Genome-wide association with bone mass and geometry in the Framingham Heart Study.BMC Med Genet8Suppl 1S14

12. Richards JB, Rivadeneira F, Inouye M, Pastinen TM, Soranzo N, et al. (2008) Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study. Lancet 371: 1505–1512.JB RichardsF. RivadeneiraM. InouyeTM PastinenN. Soranzo2008Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study.Lancet37115051512

13. Styrkarsdottir U, Halldorsson BV, Gretarsdottir S, Gudbjartsson DF, Walters GB, et al. (2008) Multiple genetic loci for bone mineral density and fractures. N Engl J Med 358: 2355–2365.U. StyrkarsdottirBV HalldorssonS. GretarsdottirDF GudbjartssonGB Walters2008Multiple genetic loci for bone mineral density and fractures.N Engl J Med35823552365

14. Rivadeneira F, Styrkarsdottir U, Estrada K, Halldorsson BV, Hsu YH, et al. (2009) Twenty bone-mineral-density loci identified by large-scale meta-analysis of genome-wide association studies. Nat Genet 41: 1199–1206.F. RivadeneiraU. StyrkarsdottirK. EstradaBV HalldorssonYH Hsu2009Twenty bone-mineral-density loci identified by large-scale meta-analysis of genome-wide association studies.Nat Genet4111991206

15. Styrkarsdottir U, Halldorsson BV, Gretarsdottir S, Gudbjartsson DF, Walters GB, et al. (2009) New sequence variants associated with bone mineral density. Nat Genet 41: 15–17.U. StyrkarsdottirBV HalldorssonS. GretarsdottirDF GudbjartssonGB Walters2009New sequence variants associated with bone mineral density.Nat Genet411517

16. Perrone M (June 01, 2010) FDA clears Amgen's bone-strengthening drug Prolia. M. PerroneJune 01, 2010FDA clears Amgen's bone-strengthening drug Prolia.Available: http://www.biosciencetechnology.com/News/FeedsAP/2010/06/fda-clears-amgens-bone-strengthening-drug-prolia/ Accessed 2011 September 14. Available: http://www.biosciencetechnology.com/News/FeedsAP/2010/06/fda-clears-amgens-bone-strengthening-drug-prolia/ Accessed 2011 September 14.

17. McClune BL, Polgreen LE, Burmeister LA, Blaes AH, Mulrooney DA, et al. (2011) Screening, prevention and management of osteoporosis and bone loss in adult and pediatric hematopoietic cell transplant recipients. Bone Marrow Transplant 46: 1–9.BL McCluneLE PolgreenLA BurmeisterAH BlaesDA Mulrooney2011Screening, prevention and management of osteoporosis and bone loss in adult and pediatric hematopoietic cell transplant recipients.Bone Marrow Transplant4619

18. Ambrosone CB, Nesline MK, Davis W (2006) Establishing a cancer center data bank and biorepository for multidisciplinary research. Cancer Epidemiol Biomarkers Prev 15: 1575–1577.CB AmbrosoneMK NeslineW. Davis2006Establishing a cancer center data bank and biorepository for multidisciplinary research.Cancer Epidemiol Biomarkers Prev1515751577

19. de Bakker PI, Yelensky R, Pe'er I, Gabriel SB, Daly MJ, et al. (2005) Efficiency and power in genetic association studies. Nat Genet 37: 1217–1223.PI de BakkerR. YelenskyI. Pe'erSB GabrielMJ Daly2005Efficiency and power in genetic association studies.Nat Genet3712171223

20. The International HapMap Consortium (2003) The International HapMap Project. Nature 426: 789–796.The International HapMap Consortium2003The International HapMap Project.Nature426789796

21. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.S. PurcellB. NealeK. Todd-BrownL. ThomasMA Ferreira2007PLINK: a tool set for whole-genome association and population-based linkage analyses.Am J Hum Genet81559575

22. Chen B, Wilkening S, Drechsel M, Hemminki K (2009) SNP_tools: A compact tool package for analysis and conversion of genotype data for MS-Excel. BMC Res Notes 2: 214.B. ChenS. WilkeningM. DrechselK. Hemminki2009SNP_tools: A compact tool package for analysis and conversion of genotype data for MS-Excel.BMC Res Notes2214

23. Gamazon ER, Zhang W, Konkashbaev A, Duan S, Kistner EO, et al. (2010) SCAN: SNP and copy number annotation. Bioinformatics 26: 259–262.ER GamazonW. ZhangA. KonkashbaevS. DuanEO Kistner2010SCAN: SNP and copy number annotation.Bioinformatics26259262

24. Murabito JM, Rosenberg CL, Finger D, Kreger BE, Levy D, et al. (2007) A genome-wide association study of breast and prostate cancer in the NHLBI's Framingham Heart Study. BMC Med Genet 8: Suppl 1S6.JM MurabitoCL RosenbergD. FingerBE KregerD. Levy2007A genome-wide association study of breast and prostate cancer in the NHLBI's Framingham Heart Study.BMC Med Genet8Suppl 1S6

25. Amling M, Herden S, Posl M, Hahn M, Ritzel H, et al. (1996) Heterogeneity of the skeleton: comparison of the trabecular microarchitecture of the spine, the iliac crest, the femur, and the calcaneus. J Bone Miner Res 11: 36–45.M. AmlingS. HerdenM. PoslM. HahnH. Ritzel1996Heterogeneity of the skeleton: comparison of the trabecular microarchitecture of the spine, the iliac crest, the femur, and the calcaneus.J Bone Miner Res113645

26. Yao S, McCarthy PL, Dunford LM, Roy DM, Brown K, et al. (2008) High prevalence of early-onset osteopenia/osteoporosis after allogeneic stem cell transplantation and improvement after bisphosphonate therapy. Bone Marrow Transplant 41: 393–398.S. YaoPL McCarthyLM DunfordDM RoyK. Brown2008High prevalence of early-onset osteopenia/osteoporosis after allogeneic stem cell transplantation and improvement after bisphosphonate therapy.Bone Marrow Transplant41393398