Performance of FRAX in predicting fracture in the US postmenopausal women with varied race and genetic profiles
===============================================================================================================

* Qing Wu
* Xiangxue Xiao
* Yingke Xu

## Abstract

**Background** Whether the Fracture Risk Assessment Tool (FRAX) performed differently in estimating the 10-year fracture probability in women of different genetic profiling and race remained unclear.

**Methods** The genomic data in the Women’s Health Initiative study was analyzed (n=23,981). the genetic risk score (GRS) was calculated from 14 fracture-associated single nucleotide polymorphisms (SNPs) for each participant. FRAX without bone mineral density (BMD) was used to estimate fracture probability.

**Results** FRAX significantly overestimated the risk of major osteoporotic fracture (MOF) in the WHI study. The most enormous overestimation was observed in women with low GRS (predicted/observed ratio [POR]: 1.61, 95% CI: 1.45-1.79), in Asian women (POR: 3.5, 95% CI 2.48-4.81), and in African American women (POR: 2.59, 95% CI: 2.33-2.87). Compared to the low GRS group, the 10-year probability of MOF adjusted for the FRAX score was 21% and 30% higher in median GRS group and high GRS group, respectively. Asian, African American, and Hispanic women respectively had a 78%, 76%, and 56% lower hazard than Caucasian women after the FRAX score was adjusted for. The results were similar when for hip fractures.

**Conclusions** Our study suggested the FRAX performance varies significantly by both genetic profiling and race in postmenopausal women.

Keywords
*   genetic risk score (GRS)
*   bone mineral density (BMD)
*   single nucleotide polymorphism (SNP)
*   Fracture risk assessment tool (FRAX)

## 1. Introduction

Osteoporotic fracture continues to be a critical public health problem worldwide [1, 2]. One main reason is that the incidence of osteoporotic fracture increases exponentially throughout one’s life [3]. Approximately 40% of postmenopausal women will suffer at least one fracture in their lifetime [4-6]. Additionally, bone fractures often lead to devastating consequences, including functional decline, prolonged disability, and death [7]. With longevity increasing globally, the potentially high cumulative rate of osteoporosis and fractures, and the associated excess disability and mortality will lead to an inevitable increase in social and economic burdens worldwide [8, 9].

Furthermore, osteoporosis is a silent disease because bone loss occurs without any signs or symptoms [3]. Patients often do not aware that they have osteoporosis until a fracture occurs; thus, fracture prediction becomes critically important. Bone mineral density (BMD) can be used to stratify patients for fracture risk; however, it has low sensitivity [10]. Studies have found that most fractures occur in individuals who have a BMD above the threshold for osteoporosis [11, 12]. Because many factors other than BMD, such as age, gender, weight, height, smoking, alcohol consumption, medication, also contribute fracture risk, several algorithms have been developed to integrate these risk factors to assess fracture risk [16-20] more accurately. The Fracture Risk Assessment Tool (FRAX), which is the most widely used tool for fracture risk assessment, was developed by the Collaborating Centre for Metabolic Bone Diseases (Sheffield, UK). The FRAX is a computer-based program that computes the 10-year probability of major osteoporotic fracture (MOF, a composite of hip, humerus, forearm, and clinical vertebral fractures) and the hip fracture. The FRAX can be used with or without femoral neck BMD measurement [13]. Although FRAX improves fracture prediction over the BMD T-score method alone [14], the FRAX performance of predicting fracture risk varies in different study populations [15-17]. Hence, there is room for further improvement in fracture prediction.

FRAX was derived from nine cohorts and has been validated in 11 independent cohorts from around the world [14]. The US FRAX was calibrated from the data of the Rochester Epidemiology Project [18] cohort, which was composed predominantly of Caucasians [19]. For the US minorities, the FRAX estimates were adjusted basing on race-specific hip fracture incidence rates and race-specific mortality [20]. This adjustment was not empirically based. Racial/ethnic differences that influence fracture risk may not be adequately taken into account by US FRAX [21]. Additional studies are needed to examine the performance of FRAX in U. S. minorities.

Additionally, genetic profiling is an essential predisposition to bone deterioration and fragility fractures [22]. Genetic factors are also determinants of bone structure [23]. Although FRAX does not factor in genetic elements, mounting evidence shows that fracture susceptibility is genetically determined [24]. Virtually 50% of the variance in liability to fragility fracture is attributable to genetic elements [25]. With the advancement of genomic technologies in the past two decades, major genome-wide association studies (GWASs) and genome-wide meta-analyses have successfully identified numerous genetic loci associated with fracture [18, 26, 27]. To date, the largest genome-wide meta-analysis on fracture, which involved 32,961 participants, revealed 14 single nucleotide polymorphisms (SNPs) associated with fracture [18]. However, it remains unclear how these SNPs cause bone fragility and associated fracture. As the allelic frequency of these discovered SNPs featured high variability in the population, and each SNP is associated with small effect size, the contribution of any single SNP to fracture susceptibility is expected to be minimal [28]. The cumulative effects of many associated genetic variants possibly cause osteoporotic fracture [29, 30]. Thus polygenic scores summarized from risk alleles at each locus have commonly been employed to quantify the overall genetic effect contributing to fracture risk [31].

The performance of FRAX with different genetic profiling has not been reported in the literature. The performance of FRAX in minorities of the US was rarely studied. Thus our study aimed 1) to evaluate whether FRAX performs differently in estimating the 10-year absolute probability of MOF and hip fracture in postmenopausal women with different polygenic risk scores, and 2) assess FRAX performance in the prediction of MOF hip fracture in minority women. We also examined if the interaction of race and polygene score impacts the performance of FRAX in fracture prediction.

## 2. Experimental Section

### 2.1. Data Source

The WHI study is a nationwide, long-term health study that has focused on heart disease, breast and colorectal cancer, and fragility fractures in postmenopausal women. Between 1993 and 1998, the WHI has enrolled 161,808 women aged 50 to 79 years into one or more randomized Clinical Trials or to an Observational Study (OS), which were conducted at 40 clinical centers nationwide. Participants were provided by mail or telephone with questionnaires annually in the observational study, or semiannually in the clinical trials. The Institutional Review Board at each participating institution approved study protocols and consent forms [32].

### 2.2. Participants

The female participants included in the present study were combined from four WHI sub-studies: WHI Genomics and Randomized Trials Network (GARNET); National Heart Lung and Blood Institute (NHLBI); Population Architecture using Genomics and Epidemiology (PAGE); and Women’s Health Initiative Memory Study (WHIMS). These data were acquired through the database of Genotype and Phenotype (dbGap) ([https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000200.v12.p3](https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000200.v12.p3)) with the approval of the institutional review board at the University of Nevada, Las Vegas. The included participants were genotyped using either the Illumina (Illumina, San Diego, California) or Affymetrix 6.0 Array Set Platform (Affymetrix, Santa Clara, California). Participants who reported taking any medication known to influence osteoporosis, including bisphosphonates, calcitonin, parathyroid hormone, selective estrogen receptor modulators, luteinizing hormone-releasing hormone agents, and somatostatin agents, and participants who had incomplete information regarding risk factors included in FRAX, were excluded from the study sample. In total, 1,513 subjects were excluded from the present study, and there were 23,981 eligible participants from multiple racial backgrounds, genotype data, and adjudicated fracture outcomes available.

### 2.3. Outcomes: Incident Fractures

The primary outcomes are major MOF and hip fractures. The WHI participants were followed for 9 years on average from the baseline examination. The follow-up period was calculated from the enrollment (OS) or randomization (CT) to the time of the first fracture or death. People who did not experience a fracture or death were followed until the end of the initial WHI study. Self-reported fracture outcomes were identified by questionnaires semiannually for CT participants and annually for OS participants. Radiology reports were used to adjudicate all fractures in the CT, and hip fractures in the OS. Hip fractures were centrally or locally adjudicated using the same criteria. The agreement between central and local adjudication for hip fracture was 96%. Other types of fractures were locally adjudicated at clinical centers where BMD was not measured [32].

### 2.4. Genotyping

Genotype data produced from blood samples were acquired through dbGap. Genotype imputation was completed at the Sanger Imputation Server. The Haplotype Reference Consortium (HRC) reference panel and Positional Burrows-Wheeler Transform (PBWT) imputation algorithm were used for genotype imputation. All 14 fracture-associated SNPs, as reported by Estrada et al. [18], were successfully imputed.

### 2.5. Genetic risk scores

Genetic risk for fracture was quantified using a standardized metric described by Estrada et al. [18]. Briefly, this metric allows the composite assessment of genetic risk in complex traits by summarizing the genetic predisposition. Based on 14 fracture-associated SNPs discovered in the largest genome-wide meta-analysis [18], weighted GRS was calculated for each participant in this study as GRS = sum (x\_i*b_i); where x_i are individual’s genotype (0, 1, 2) for SNP i, and b_i are the effect size of this SNP. Linkage disequilibrium (LD) pruning was performed in advance to eliminate possible LD that existed between SNPs. None of the 14 SNPs were removed after pruning. To illustrate the different cumulative incidence of fracture in participants with different genetic profiles, we divided the participants into three GRS groups based on their weighted GRS, using distributions of 25%, 50%, and 25%.

### 2.6. Fracture probability

Since BMD measurement was unavailable for over 90% of the participants in this study sample, the existing FRAX score calculated by the FRAX without BMD (US FRAX version 3.0) for the 10-year probability of MOF and hip fracture in the data was used. The performance of FRAX with BMD according to race and GRS will be addressed in a different study. All predicted fracture probability in this study was estimated by FRAX without BMD, unless noted differently. The observed 10-year cumulative fracture incidence was assessed by race and GRS groups. The cumulative incidence function (CIF) was applied to derive the observed 10-year fracture probability for MOF, and hip fracture accounting for competing mortality risk [33].

### 2.7. Statistical Analysis

Demographic and baseline clinical characteristics are presented as mean ± SD for continuous variables or frequencies (%) for categorical variables. Differences between the two groups were examined by using Student’s t-test for continuous variables and by using chi-square tests for categorical variables, respectively. The ratio between FRAX predicted fracture probability and observed fracture probability (POR), with the corresponding 95% CI, was calculated for each group. Multivariable Cox proportional hazard models were employed to assess the effect GRS and race had on survival time to the first fracture or death with adjusting for baseline FRAX probability.

A sensitivity analysis was conducted on a small sample (N=14,722) in which participants who received interventions in any of the three clinical trials, namely Calcium and Vitamin D Trial (CAD), Hormone Therapy Trial (HT), and Dietary Modification Trial (DM) were excluded. Statistical analyses were performed using the SAS 9.4 (SAS Institute, Inc., Cary, NC, USA).

## 3. Results

This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation as well as the experimental conclusions that can be drawn.

### 3.1. Baseline characteristics

The study included a total of 23,981 women for analysis. During an average of 12 years of follow-up, 5,555 (23.23%) women died, and 1,637 (6.9%) women sustained at least one MOF during the follow-up. **Table 1**. compares the baseline characteristics of women with MOF and women without an MOF during the follow-up. Weighted GRS was significantly higher in women who sustained an MOF than in those who did not (p<.0001). Women who sustained an MOF were also older, had lower body mass index (BMI), more alcohol consumption, a higher prevalence of prior fractures, and more hip fractures in their family history. FRAX scores of both MOF and hip were significantly higher in women with a fracture incidence (p<.0001). The mean of GRS between races are significantly different (Supplemental Figure 1).

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/01/06/2020.01.03.20016444/T1)

Table 1. 
Baseline Characteristics of 23,981 women according to whether they sustained a MOF during follow-up

![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/06/2020.01.03.20016444/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/01/06/2020.01.03.20016444/F1)

Figure 1. 
Crude (unadjusted) 10-year cumulative incidence of major osteoporotic (A) and hip fracture (B) stratified by the GRS group, including competing mortality risk. The difference in the cumulative incidence rates among different GRS groups was tested by using Gray’s test, p-value<0.01 indicating a significant difference between the groups.

### 3.2. Performance of FRAX in predicting MOF and hip fracture

The crude 10-year cumulative incidence of MOF and hip fracture by the GRS group is shown in **Figure 1**. Significant between-group differences were observed for both MOF (p<.001) and hip fracture (p<.001). The incidences of MOF and hip fracture were greater in the high GRS group. The crude 10-year cumulative incidence of MOF and hip fracture by race is shown in **Figure 2**. Significant between-group differences were observed for both MOF (p<.001) and hip fracture (p<.001). The incidence of MOF and hip fracture were higher in Caucasian women.

![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/06/2020.01.03.20016444/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/01/06/2020.01.03.20016444/F2)

Figure 2. 
Crude (unadjusted) 10-year cumulative incidence of major osteoporotic (A) and hip fracture (B) stratified by race, including competing mortality risk. The difference in the cumulative incidence rates among different racial groups was tested by using Gray’s test, p-value<0.01 indicating a significant difference between the groups.

The predicted versus observed 10-year probability of MOF by GRS groups, with accounting for competing mortality, are shown in **Figure 3A**. The 10-year MOF probability derived from FRAX significantly overestimated risk across all GRS groups. The greatest overestimation by FRAX was observed in women who had low GRS, in which the 10-year predicted probability of MOF was 6.02% versus observed 3.74%, with a corresponding predicted/observed ratio (POR) of 1.61 (95% CI, 1.45-1.79), followed by the high GRS group with a POR of 1.38 (95% CI, 1.27-1.50); and in the median GRS group, the POR was 1.40 (95% CI, 1.32-1.49). For hip fracture outcome, where the 10-year predicted probability calculated by FRAX overestimated fracture risk in all GRS groups, however, with similar POR across the three GRS groups (**Figure 3B**).

![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/06/2020.01.03.20016444/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/01/06/2020.01.03.20016444/F3)

Figure 3. 
Observed versus predicted 10-year major osteoporotic fracture (A) and hip fracture (B) probability stratified by the GRS group. The dotted line indicates a relative ratio of 1 (reference line), ratio >1 indicates that FRAX overestimates fracture probability.

The predicted versus observed 10-year probability of MOF by racial groups, with competing mortality risk accounted for are shown in **Figure 4A**. The 10-year probability of MOF calculated by FRAX significantly overestimated fracture risk in most racial groups, and the greatest overestimation was observed in Asian women. In Asian women, the predicted 10-year probability of MOF was 7.26% versus observed 2.03%, and the POR was 3.5 (95% CI 2.48-4.81). In African American women, the predicted 10-year probability of MOF was 3.79% as opposed to observed 1.46%, with the POR being 2.59 (95% CI 2.33-2.87). The 10-year probability of hip fracture estimated without BMD overestimated risk in all racial groups except American Indians. The 10-year predicted probability of hip fracture was in this group was 1.75% as opposed to observed 1.91%, with the POR being 0.91 (95% CI 0.46-1.62) (**Figure 4B**).

![Figure 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/06/2020.01.03.20016444/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/01/06/2020.01.03.20016444/F4)

Figure 4. 
Observed versus predicted 10-year major osteoporotic fracture (A) and hip fracture (B) probability stratified by race. The dotted line indicates a relative ratio of 1 (reference line), ratio >1 indicates that FRAX overestimates fracture probability.

### 3.3. Race/ethnicity and the fracture outcome

In the multivariate Cox proportional hazard model, after adjusting for baseline FRAX probability in the model, weighted GRS calculated from 14 fracture-related SNPs was significantly associated with subsequent MOF. Compared to the low GRS group, the 10-year probability of MOF was 21% higher for women with medium genetic risk (HR=1.21, 95% CI 1.05-1.39) and 30% higher for women with high genetic risk (HR=1.30, 95% CI 1.12-1.50). Similar findings with hip fracture outcomes were observed. Compared to the low GRS group, the 10-year probability of hip fracture was 27% higher for women in the medium GRS group (HR=1.27, 95% CI 1.04-1.55) and 46% higher for women in the high GRS group (HR=1.46, 95% CI 1.17-1.80) (**Table 2**).

View this table:
[Table 2.](http://medrxiv.org/content/early/2020/01/06/2020.01.03.20016444/T2)

Table 2. 
Hazard Ratios (HR) with 95% confidence interval (CI) for outcomes of incidence fracture according to the GRS group, adjusted for FRAX score: Results of multivariate Cox proportional hazard model.

After controlling for baseline fracture probability estimated by FRAX, race remained a significant predictor of subsequent MOF and hip fracture. Compared to Caucasian women, Asian women had a 78% lower hazard of MOF (HR=0.22, 95% CI 0.12-0.41) and hip fracture (HR=0.22, 95% CI 0.09-0.52). Similarly, the FRAX-adjusted hazard ratio of MOF and hip fracture for African-American women and Hispanic women were also significantly lower (**Table 3**).

View this table:
[Table 3.](http://medrxiv.org/content/early/2020/01/06/2020.01.03.20016444/T3)

Table 3. 
Hazard Ratios (HR) with 95% confidence interval (CI) for outcomes of incidence fracture according to race, adjusted for FRAX score: Results of multivariate Cox proportional hazard model.

### 3.4. GRS and the fracture outcome

The potential impact of GRS on the estimated probabilities of MOF and hip fractures across different racial groups was also assessed. When adjusted for FRAX probability and race, high GRS was associated with an increased probability of MOF (GRS high vs. low: HR=1.08, 95% CI 0.92-1.25) and of hip fracture (HR=1.17, 95% CI 0.93-1.46) (**Table 2**); however, the increase was not statistically significant in both outcomes. When adjusted for the baseline FRAX probability and GRS simultaneously, the impact of race on the estimated probabilities MOF and hip fracture was slightly attenuated but remained statistically significant. Compared to Caucasian women, American Indians, Asians, African American, and Hispanic women had a 61%, 78%, 76%, and 57% lower hazard of MOF, respectively. Similar findings were observed with hip fracture outcomes (**Table 3**).

### 3.5. Sensitivity analysis

We conducted a sensitivity analysis in which we excluded subjects who received intervention in either of the three clinical trials. When adjusted for the FRAX probability, we observed an increased HR of MOF with GRS (GRS high as opposed to low: HR=1.39, 95% CI 1.10 – 1.74); whereas the impact of GRS on the estimated probability of hip fracture attenuated slightly (HR=1.41, 95% CI 1.03-1.95), but remained significant. However, when adjusted for both race and FRAX probability, the association between GRS and hip fracture was not significant (Supplementary Table 1). When controlling for GRS and FRAC probability, the effects of race on the estimated probabilities of MOF and hip fracture remained significant. Compare to Caucasian women, the race and FRAX-adjusted hazard of MOF was 90%, 78%, and 66% lower in Asian, African-American, and Hispanic women, respectively (**Supplementary Table 2**).

## 4. Discussion

The present study found that FRAX overestimated the risk of fracture in women aged 50-79 years, and the degree of overestimation by FRAX in the low GRS group is greater than high genetic risk groups in both outcomes of MOF and hip fracture. In the multivariate analysis, genetic profiling was further demonstrated to be a significant predictor of MOF and hip fracture, independent of FRAX probability.

Genetic factors that influence osteoporotic fracture risk have been long recognized. Genetics are determinants of bone structure and thus a predisposition to fragility. Hereditary factors contribute almost half of the variance in fracture susceptibility [25]. However, genetic factors are not counted in the FRAX or any other existing clinical fracture risk assessment models. Since FRAX is the most commonly used fracture prediction model, determining if the performance of FRAX varies with different genetic profiling has become crucially important. The largest and most updated GWAS meta-analysis has identified 14 SNPs that are significantly associated with fracture risk at a genome-wide significant level [25]. Although these individual SNPs have modest effect size on fracture risk, the GRS, as summarized from these individual risk SNPs, enables us to examine if FRAX performance varies with different genetic risk factors. The varied prediction performance of FRAX by GRS observed in our study suggests that the accuracy of FRAX can be improved by incorporating genetic profiling. Several studies suggested that genetic profiling may help improve the accuracy of various fracture prediction models. For example, GRS of 39 SNPs increased the precision of nonvertebral fracture prediction in postmenopausal Korean women [34]. Additionally, a profiling of 63 SNPs improved the accuracy of non-trauma fracture prediction [29]. One of our recent studies on the US older men also found that GRS is one of the most important variables in MOF prediction models developed by the gradient boosting approach (manuscript under review).

The present study also provides compelling evidence that FRAX overestimates the risk of MOF and hip fracture in women 50-79 years old, across all racial groups, but especially in minorities. In Asian, African-American, and Hispanic women, the observed probability of fracture, in terms of both MOF and hip fracture, was significantly lower than the risk estimated by FRAX, indicating that the FRAX did not adequately capture racial and ethnic differences of fracture risk. Additionally, our multivariate analysis demonstrated that race is a significant predictor of MOF and hip fracture independent of the cumulative fracture risk estimated by FRAX, suggesting that FRAX does not have adequate adjustment for racial difference. Racial and ethnic difference that influences fracture risk not being adjusted for adequately in the FRAX has long been a concern [13]. As we know, the US FRAX was calibrated from the REP data, composed predominantly of Caucasians. For non-Caucasian minorities, the FRAX estimates were adjusted based on race-specific hip fracture incidence rates and race-specific mortality [46]. This adjustment for minorities in FRAX is not empirically based, thus making the prediction accuracy of FRAX increasingly uncertain, especially for MOF, a composite of hip, humerus, forearm, and clinical vertebral fractures. The current FRAX adjustment model, based on race-specific hip fracture incidence rate and race-specific mortality, remains likely to be inadequate for MOF risk estimation in minorities. In this study, we observed that the overestimated risk for MOF by FRAX was much higher than that for hip fracture, at least validated that the US FRAX has not adjusted race adequately for MOF. A prior study conducted on the same WHI sample assessing the accuracy of FRAX without BMD in predicting fracture also demonstrated that the FRAX has significant lower sensitivity in identifying incidence fractures in African-American and Hispanic women [35]. Another study on 2266 postmenopausal women who participated in the Hong Kong Osteoporosis study revealed that the predictive accuracy of FRAX with BMD was not substantially different from the model with BMD alone [36]. Considering the generally lower incidence of fracture in Asians than in Caucasians and the prominent effect of BMD in fracture prediction in the Asian population, the absence of BMD in the present study may explain the significant overestimation of fracture in this racial group. Besides, inconsistent findings regarding the performance of FRAX without BMD was reported in several other studies. Leslie et al. observed that the fracture probability estimated without BMD overestimate risk among the general population [37], which are consistent with findings from the present study. Other studies have reported underestimation of fracture risk by FRAX [15, 38, 39], but their methodologies were found to be problematic lately because they either comparing incidence with probabilities or failed to take the competing mortality risk into account [40].

When both FRAX probability and race were adjusted simultaneously in the multivariate model, the effect of GRS was reduced, which be due to the following reasons. First, genetic profiling regarding osteoporosis or osteoporotic fracture varies in different racial groups; the effect of race and GRS on fracture could be overlapping (See Supplemental Figure 1). Second, the genetic effect on fracture probability may not be fully captured by the limited number of discovered risk SNPs.

With more fracture-related genetic components being discovered in the future, a larger effect of GRS on fracture risk prediction should be foreseen.

Limitations to this study are acknowledged. First, the WHI data we used only included women 50-79 years, so our findings may not apply to men or to women who are not in the study age range. Second, rare genetic variants with high effect size were not included in the present study, because risk SNPs used in this analysis were identified from a GWAS meta-analysis study, which likely discovered common variants, not rare variants [18]. The limited number of fracture-associated SNPs may not capture all genetic risk, which partially explained the reduced effect of GRS in the model when both FRAX probability and race were included. Third, our study only focuses on FRAX without BMD because the BMD measurement was unavailable for most of the study subjects. The performance of FRAX with BMD will be examined in a future study. Finally, the sample size of Asian and American Indian subjects were very small in this study; the results may, therefore, be underpowered.

## 5. Conclusions

To the best of our knowledge, this is the first study to assess FRAX performance in the prediction of MOF and hip fractures in groups with different genetic profiling and of various races. Our findings suggested genetic profiling of an individual should be considered in fracture prediction, as genetic factors have been demonstrated to be a significant risk factor for osteoporotic fracture, independent of FRAX. Our results also demonstrated that FRAX performed differently in different races, and thus the effect of race in osteoporotic fracture prediction has not heretofore adequately been taken accounted for by existing FRAX models. Fully integrating genetic profiling and racial factors into the existing fracture assessment model is very likely to improve the accuracy of fracture prediction. Thus, developing racial/ethnic-specific, individualized fracture risk assessment models will provide more accurate fracture risk assessment. Further studies, especially these including men, larger sample of minorities, and more comprehensive fracture-associated genetic variants, are warranted.

## Data Availability

The data/analyses presented in the current publication are based on the use of study data downloaded from the dbGaP web site, under phs000200

[https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study\_id=phs000200.v12.p3](https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000200.v12.p3) 

## Supplementary Materials

The following are available online at [www.mdpi.com/xxx/s1](http://www.mdpi.com/xxx/s1), Figure S1: Box plot for the distribution of GRS in different race groups (N=23,981), Table S1: Hazard Ratios (HR) with 95% confidence interval (CI) for outcomes of Incidence fracture according to GRS group, adjusted for FRAX score: results from sensitivity analysis after excluding subjects who received intervention in either of the three clinical trials. (N=14,722), Table S2: Hazard Ratios (HR) with 95% confidence interval (CI) for outcomes of incidence fracture according to race, adjusted for FRAX score: results from sensitivity analysis after excluding subjects who received intervention in either of the three clinical trials. (N=14,722)

## Author Contributions

conceptualization, Q.W.; methodology, Q.W., X.X. and Y.X.; software, Y.X.; validation, Q.W., X.X. and Y.X.; formal analysis, Y.X.; investigation, Q.W., X.X. and Y.X.; resources, Q.W.; data curation, X.X. and Y.X.; writing—original draft preparation, Q.W., X.X.; writing—review and editing, Q.W. X.X. and Y.X; visualization, X.X. and Y.X; supervision, Q.W.; project administration, Q.W.; funding acquisition, Q.W.

## Funding

This research was funded by a grant from the National Institute of General Medical Sciences (GR08954), a grant from the National Institute on Minority Health and Health Disparities of the National Institutes of Health (R15MD010475).

## Conflicts of Interest

The authors declare no conflict of interest.

## Acknowledgments

The data/analyses presented in the current publication are based on the use of study data downloaded from the dbGaP web site, under phs000200 ([https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000200.v12.p3](https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000200.v12.p3)). The research and analysis described in the current publication were supported by the Genome Acquisition to Analytics (GAA) Research Core of the Personalized Medicine Center of Biomedical Research Excellence in the Nevada Institute of Personalized Medicine, and the National Supercomputing Institute at the University of Nevada Las Vegas provided facilities for bioinformatical analysis in this study.

*   Received January 3, 2020.
*   Revision received January 3, 2020.
*   Accepted January 6, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Johnell, O. & Kanis, J. A. (2004) An estimate of the worldwide prevalence, mortality and disability associated with hip fracture, Osteoporosis international : a journal established as result of cooperation between the European Foundation for Osteoporosis and the National Osteoporosis Foundation of the USA. 15, 897–902.
    
    

2.  2.Johnell, O. & Kanis, J. A. (2006) An estimate of the worldwide prevalence and disability associated with osteoporotic fractures, Osteoporosis international : a journal established as result of cooperation between the European Foundation for Osteoporosis and the National Osteoporosis Foundation of the USA. 17, 1726–33.
    
    

3.  3.Sözen, T., Özışık, L. & Başaran, N. Ç. (2017) An overview and management of osteoporosis, Eur J Rheumatol. 4, 46–56.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 

4.  4.Reginster, J. Y. & Burlet, N. (2006) Osteoporosis: a still increasing prevalence, Bone. 38, S4–9.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16455317&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 

5.  5.Watts, N. B., Bilezikian, J. P., Camacho, P. M., Greenspan, S. L., Harris, S. T., Hodgson, S. F., Kleerekoper, M., Luckey, M. M., McClung, M. R., Pollack, R. P. & Petak, S. M. (2010) American Association of Clinical Endocrinologists Medical Guidelines for Clinical Practice for the diagnosis and treatment of postmenopausal osteoporosis, Endocrine practice : official journal of the American College of Endocrinology and the American Association of Clinical Endocrinologists. 16 Suppl 3, 1–37.
    
    

6.  6.Wright, N. C., Looker, A. C., Saag, K. G., Curtis, J. R., Delzell, E. S., Randall, S. & Dawson-Hughes, B. (2014) The recent prevalence of osteoporosis and low bone mass in the United States based on bone mineral density at the femoral neck or lumbar spine, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 29, 2520–2526.
    
    

7.  7.Dyer, S. M., Crotty, M., Fairhall, N., Magaziner, J., Beaupre, L. A., Cameron, I. D., Sherrington, C. & Fragility Fracture Network Rehabilitation Research Special Interest, G. (2016) A critical review of the long-term disability outcomes following hip fracture, BMC Geriatr. 16, 158–158.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 

8.  8.Sànchez-Riera, L., Carnahan, E., Vos, T., Veerman, L., Norman, R., Lim, S. S., Hoy, D., Smith, E., Wilson, N., Nolla, J. M., Chen, J. S., Macara, M., Kamalaraj, N., Li, Y., Kok, C., Santos-Hernández, C. & March, L. (2014) The global burden attributable to low bone mineral density, Annals of the Rheumatic Diseases. 73, 1635.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjk6IjczLzkvMTYzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAxLzA2LzIwMjAuMDEuMDMuMjAwMTY0NDQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

9.  9.Harvey, N., Dennison, E. & Cooper, C. (2010) Osteoporosis: impact on health and economics, Nature Reviews Rheumatology. 6, 99–105.
    
    

10. 10.Gnudi, S. & Malavolta, N. (2003) Comparison Between T-Score-Based Diagnosis of Osteoporosis and Specific Skeletal Site Measurements: Prognostic Value for Predicting Fracture Risk, Journal of Clinical Densitometry. 6, 267–273.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14514997&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 

11. 11.Schuit, S. C., van der Klift, M., Weel, A. E., de Laet, C. E., Burger, H., Seeman, E., Hofman, A., Uitterlinden, A. G., van Leeuwen, J. P. & Pols, H. A. (2004) Fracture incidence and association with bone mineral density in elderly men and women: the Rotterdam Study, Bone. 34, 195–202.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bone.2003.10.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14751578&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000188861100021&link_type=ISI) 

12. 12.Siris, E. S., Chen, Y. T., Abbott, T. A., Barrett-Connor, E., Miller, P. D., Wehren, L. E. & Berger, M. L. (2004) Bone mineral density thresholds for pharmacological intervention to prevent fractures, Arch Intern Med. 164, 1108–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archinte.164.10.1108&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15159268&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221623600008&link_type=ISI) 

13. 13.Kanis, J. A. (2008) Assessment of osteoporosis at the primary health-care level. Technical Report, [http://www.shefacuk/FRAX](http://www.shefacuk/FRAX).
    
    

14. 14.Kanis, J. A., Oden, A., Johnell, O., Johansson, H., De Laet, C., Brown, J., Burckhardt, P., Cooper, C., Christiansen, C., Cummings, S., Eisman, J. A., Fujiwara, S., Gluer, C., Goltzman, D., Hans, D., Krieg, M. A., La Croix, A., McCloskey, E., Mellstrom, D., Melton, L. J., 3rd, Pols, H., Reeve, J., Sanders, K., Schott, A. M., Silman, A., Torgerson, D., van Staa, T., Watts, N. B. & Yoshimura, N. (2007) The use of clinical risk factors enhances the performance of BMD in the prediction of hip and osteoporotic fractures in men and women, Osteoporosis international : a journal established as result of cooperation between the European Foundation for Osteoporosis and the National Osteoporosis Foundation of the USA. 18, 1033–46.
    
    

15. 15.Sornay-Rendu, E., Munoz, F., Delmas, P. D. & Chapurlat, R. D. (2010) The FRAX tool in French women: How well does it describe the real incidence of fracture in the OFELY cohort?, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 25, 2101–7.
    
    

16. 16.Crandall, C. J., Schousboe, J. T., Morin, S. N., Lix, L. M. & Leslie, W. (2019) Performance of FRAX and FRAX-Based Treatment Thresholds in Women Aged 40 Years and Older: The Manitoba BMD Registry, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 34, 1419–1427.
    
    

17. 17.Briot, K., Paternotte, S., Kolta, S., Eastell, R., Felsenberg, D., Reid, D. M., Glüer, C.-C. & Roux, C. (2014) FRAX®: Prediction of Major Osteoporotic Fractures in Women from the General Population: The OPUS Study, PLOS ONE. 8, e83436.
    
    

18. 18.Estrada, K. Styrkarsdottir,  U. Evangelou,  E. Hsu,  Y. H. Duncan,  E. L. Ntzani,  E. E. Oei,  L. Albagha,  O. M. Amin,  N. Kemp,  J. P. Koller,  D. L. Li,  G. Liu,  C. T. Minster,  R. L. Moayyeri,  A. Vandenput,  L. Willner,  D. Xiao,  S. M. Yerges-Armstrong,  L. M. Zheng,  H. F. Alonso,  N. Eriksson,  J. Kammerer,  C. M. Kaptoge,  S. K. Leo,  P. J. Thorleifsson,  G. Wilson,  S. G. Wilson,  J. F. Aalto,  V. Alen,  M. Aragaki,  A. K. Aspelund,  T. Center,  J. R. Dailiana,  Z. Duggan,  D. J. Garcia,  M. Garcia-Giralt,  N. Giroux,  S. Hallmans,  G. Hocking,  L. J. Husted,  L. B. Jameson,  K. A. Khusainova,  R. Kim,  G. S. Kooperberg,  C. Koromila,  T. Kruk,  M. Laaksonen,  M. Lacroix,  A. Z. Lee,  S. H. Leung,  P. C. Lewis,  J. R. Masi,  L. Mencej-Bedrac,  S. Nguyen,  T. V. Nogues,  X. Patel,  M. S. Prezelj,  J. Rose,  L. M. Scollen,  S. Siggeirsdottir,  K. Smith,  A. V. Svensson,  O. Trompet,  S. Trummer,  O. van Schoor,  N. M. Woo,  J. Zhu,  K. Balcells,  S. Brandi,  M. L. Buckley,  B. M. Cheng,  S. Christiansen,  C. Cooper,  C. Dedoussis,  G. Ford,  I. Frost,  M. Goltzman,  D. Gonzalez-Macias,  J. Kahonen,  M. Karlsson,  M. Khusnutdinova,  E. Koh,  J. M. Kollia,  P. Langdahl,  B. L. Leslie,  W. D. Lips,  P. Ljunggren,  O. Lorenc,  R. S. Marc,  J. Mellstrom,  D. Obermayer-Pietsch,  B. Olmos,  J. M. Pettersson-Kymmer,  U. Reid,  D. M. Riancho,  J. A. Ridker,  P. M. Rousseau,  F. Slagboom, P. E. Tang, N. L., et al. (2012) Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture, Nature genetics. 44, 491–501.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2249&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22504420&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 

19. 19.Melton, L. J., 3rd (1996) History of the Rochester Epidemiology Project, Mayo Clinic proceedings. 71, 266–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4065/71.3.266&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8594285&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996TZ98900008&link_type=ISI) 

20. 20.Dawson-Hughes, B., Tosteson, A. N., Melton, L. J., 3rd, Baim, S., Favus, M. J., Khosla, S. & Lindsay, R. L. (2008) Implications of absolute fracture risk assessment for osteoporosis practice guidelines in the USA, Osteoporosis Int. 19, 449–58.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00198-008-0559-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18292975&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253995100006&link_type=ISI) 

21. 21.(2007) WHO scientific group on the assessment of osteoporosis at the primary health-care Level in Summary Meeting Report, Brussels, Belgium, 5-7 May 2004 (the World Health Organization Collaborating Centre for Metabolic Bone Diseases, ed)
    
    

22. 22.Stewart, T. L. & Ralston, S. H. (2000) Role of genetic factors in the pathogenesis of osteoporosis, The Journal of endocrinology. 166, 235–45.
    
    [Abstract](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoiam9lIjtzOjU6InJlc2lkIjtzOjk6IjE2Ni8yLzIzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAxLzA2LzIwMjAuMDEuMDMuMjAwMTY0NDQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

23. 23.Mikkola, T. M., Sipilä, S., Rantanen, T., Sievänen, H., Suominen, H., Kaprio, J., Koskenvuo, M., Kauppinen, M. & Heinonen, A. (2008) Genetic and Environmental Influence on Structural Strength of Weight-Bearing and Non–Weight-Bearing Bone: A Twin Study, Journal of Bone and Mineral Research. 23, 492–498.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1359/jbmr.071205&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18072876&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 

24. 24.Koromani, F., Trajanoska, K., Rivadeneira, F. & Oei, L. (2019) Recent Advances in the Genetics of Fractures in Osteoporosis, Front Endocrinol (Lausanne). 10, 337–337.
    
    

25. 25.Michaëlsson, K., Melhus, H., Ferm, H., Ahlbom, A. & Pedersen, N. L. (2005) Genetic Liability to Fractures in the Elderly, Archives of Internal Medicine. 165, 1825–1830.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archinte.165.16.1825&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16157825&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231834900003&link_type=ISI) 

26. 26.Kim, S. K. (2018) Identification of 613 new loci associated with heel bone mineral density and a polygenic risk score for bone mineral density, osteoporosis and fracture, PLOS ONE. 13, e0200785.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0200785&link_type=DOI) 

27. 27.Morris, J. A., Kemp, J. P., Youlten, S. E., Laurent, L., Logan, J. G., Chai, R. C., Vulpescu, N. A., Forgetta, V., Kleinman, A., Mohanty, S. T., Sergio, C. M., Quinn, J., Nguyen-Yamamoto, L., Luco, A.-L., Vijay, J., Simon, M.-M., Pramatarova, A., Medina-Gomez, C., Trajanoska, K., Ghirardello, E. J., Butterfield, N. C., Curry, K. F., Leitch, V. D., Sparkes, P. C., Adoum, A.-T., Mannan, N. S., Komla-Ebri, D. S. K., Pollard, A. S., Dewhurst, H. F., Hassall, T. A. D., Beltejar, M.-J. G., Agee, M., Alipanahi, B., Auton, A., Bell, R. K., Bryc, K., Elson, S. L., Fontanillas, P., Furlotte, N. A., McCreight, J. C., Huber, K. E., Litterman, N. K., McIntyre, M. H., Mountain, J. L., Noblin, E. S., Northover, C. A. M., Pitts, S. J., Sathirapongsasuti, J. F., Sazonova, O. V., Shelton, J. F., Shringarpure, S., Tian, C., Tung, J. Y., Vacic, V., Wilson, C. H., Adams, D. J., Vaillancourt, S. M., Kaptoge, S., Baldock, P., Cooper, C., Reeve, J., Ntzani, E. E., Evangelou, E., Ohlsson, C., Karasik, D., Rivadeneira, F., Kiel, D. P., Tobias, J. H., Gregson, C. L., Harvey, N. C., Grundberg, E., Goltzman, D., Adams, D. J., Lelliott, C. J., Hinds, D. A., Ackert-Bicknell, C. L., Hsu, Y.-H., Maurano, M. T., Croucher, P. I., Williams, G. R., Bassett, J. H. D., Evans, D. M., Richards, J. B. & and Me Research, T. (2019) An atlas of genetic influences on osteoporosis in humans and mice, Nature genetics. 51, 258–266.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0302-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30598549&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 

28. 28.Bush, W. S. & Moore, J. H. (2012) Chapter 11: Genome-wide association studies, PLoS Comput Biol. 8, e1002822–e1002822.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1002822&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23300413&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=23300413&link_type=ISI) 

29. 29.Ho-Le, T. P., Center, J. R., Eisman, J. A., Nguyen, H. T. & Nguyen, T. V. (2017) Prediction of Bone Mineral Density and Fragility Fracture by Genetic Profiling, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 32, 285–293.
    
    

30. 30.Eriksson, J., Evans, D. S., Nielson, C. M., Shen, J., Srikanth, P., Hochberg, M., McWeeney, S., Cawthon, P. M., Wilmot, B., Zmuda, J., Tranah, G., Mirel, D. B., Challa, S., Mooney, M., Crenshaw, A., Karlsson, M., Mellstrom, D., Vandenput, L., Orwoll, E. & Ohlsson, C. (2015) Limited clinical utility of a genetic risk score for the prediction of fracture risk in elderly subjects, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 30, 184–94.
    
    

31. 31.Gibson, G. (2019) On the utilization of polygenic risk scores for therapeutic targeting, PLOS Genetics. 15, e1008060.
    
    

32. 32.(1998) Design of the Women’s Health Initiative clinical trial and observational study. The Women’s Health Initiative Study Group, Control Clin Trials. 19, 61–109.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0197-2456(97)00078-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9492970&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000071850900006&link_type=ISI) 

33. 33.Zhang, M.-J., Zhang, X. & Scheike, T. H. (2008) Modeling cumulative incidence function for competing risks data, Expert Rev Clin Pharmacol. 1, 391–400.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19829754&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F06%2F2020.01.03.20016444.atom) 

34. 34.Lee, S. H., Lee, S. W., Ahn, S. H., Kim, T., Lim, K. H., Kim, B. J., Cho, E. H., Kim, S. W., Kim, T. H., Kim, G. S., Kim, S. Y., Koh, J. M. & Kang, C. (2013) Multiple gene polymorphisms can improve prediction of nonvertebral fracture in postmenopausal women, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 28, 2156–64.
    
    

35. 35.Crandall, C. J., Larson, J., LaCroix, A., Cauley, J. A., LeBoff, M. S., Li, W., LeBlanc, E. S., Edwards, B. J., Manson, J. E. & Ensrud, K. (2019) Predicting Fracture Risk in Younger Postmenopausal Women: Comparison of the Garvan and FRAX Risk Calculators in the Women’s Health Initiative Study, Journal of general internal medicine. 34, 235–242.
    
    

36. 36.Cheung, E., Cheung, C. L., Kung, A. W. & Tan, K. C. (2014) Possible FRAX-based intervention thresholds for a cohort of Chinese postmenopausal women, Osteoporosis international : a journal established as result of cooperation between the European Foundation for Osteoporosis and the National Osteoporosis Foundation of the USA. 25, 1017–23.
    
    

37. 37.Leslie, W. D., Morin, S. N., Lix, L. M., Niraula, S., McCloskey, E. V., Johansson, H., Harvey, N. C. & Kanis, J. A. (2019) Performance of FRAX in Women with Breast Cancer Initiating Aromatase Inhibitor Therapy: A Registry-Based Cohort Study, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 34, 1428–1435.
    
    

38. 38.Bolland, M. J., Siu, A. T., Mason, B. H., Horne, A. M., Ames, R. W., Grey, A. B., Gamble, G. D. & Reid, R. (2011) Evaluation of the FRAX and Garvan fracture risk calculators in older women, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 26, 420–7.
    
    

39. 39.Pluskiewicz, W., Adamczyk, P., Franek, E., Leszczynski, P., Sewerynek, E., Wichrowska, H., Napiorkowska, L., Kostyk, T., Stuss, M., Stepien-Klos, W., Golba, K. S. & Drozdzowska, B. (2010) Ten-year probability of osteoporotic fracture in 2012 Polish women assessed by FRAX and nomogram by Nguyen et al.—Conformity between methods and their clinical utility, Bone. 46, 1661–1667.
    
    

40. 40.Kanis, J. A., Oden, A., Johansson, H. & McCloskey, E. (2012) Pitfalls in the external validation of FRAX, Osteoporosis international : a journal established as result of cooperation between the European Foundation for Osteoporosis and the National Osteoporosis Foundation of the USA. 23, 423–31.