Race-specific FRAX models are evidence-based and support equitable care: a response to the ASBMR Task Force report on Clinical Algorithms for Fracture Risk.

Task Force on 'Clinical Algorithms for Fracture Risk' commissioned by the American Society for Bone and Mineral Research (ASBMR) Professional Practice Committee has recommended that FRAX® models in the US do not include adjustment for race and ethnicity. This position paper finds that an agnostic model would unfairly discriminate against the Black, Asian and Hispanic communities and recommends the retention of ethnic and race-specific FRAX models for the US, preferably with updated data on fracture and death hazards. In contrast, the use of intervention thresholds based on a fixed bone mineral density unfairly discriminates against the Black, Asian and Hispanic communities in the US. This position of the Working Group on Epidemiology and Quality of Life of the International Osteoporosis Foundation (IOF) is endorsed both by the IOF and the European Society for Clinical and Economic Aspects of Osteoporosis, Osteoarthritis and Musculoskeletal Diseases (ESCEO).

A meta-analysis of previous falls and subsequent fracture risk in cohort studies.

The relationship between self-reported falls and fracture risk was estimated in an international meta-analysis of individual-level data from 46 prospective cohorts. Previous falls were associated with an increased fracture risk in women and men and should be considered as an additional risk factor in the FRAX® algorithm. INTRODUCTION: Previous falls are a well-documented risk factor for subsequent fracture but have not yet been incorporated into the FRAX algorithm. The aim of this study was to evaluate, in an international meta-analysis, the association between previous falls and subsequent fracture risk and its relation to sex, age, duration of follow-up, and bone mineral density (BMD). METHODS: The resource comprised 906,359 women and men (66.9% female) from 46 prospective cohorts. Previous falls were uniformly defined as any fall occurring during the previous year in 43 cohorts; the remaining three cohorts had a different question construct. The association between previous falls and fracture risk (any clinical fracture, osteoporotic fracture, major osteoporotic fracture, and hip fracture) was examined using an extension of the Poisson regression model in each cohort and each sex, followed by random-effects meta-analyses of the weighted beta coefficients. RESULTS: Falls in the past year were reported in 21.4% of individuals. During a follow-up of 9,102,207 person-years, 87,352 fractures occurred of which 19,509 were hip fractures. A previous fall was associated with a significantly increased risk of any clinical fracture both in women (hazard ratio (HR) 1.42, 95% confidence interval (CI) 1.33-1.51) and men (HR 1.53, 95% CI 1.41-1.67). The HRs were of similar magnitude for osteoporotic, major osteoporotic fracture, and hip fracture. Sex significantly modified the association between previous fall and fracture risk, with predictive values being higher in men than in women (e.g., for major osteoporotic fracture, HR 1.53 (95% CI 1.27-1.84) in men vs. HR 1.32 (95% CI 1.20-1.45) in women, P for interaction = 0.013). The HRs associated with previous falls decreased with age in women and with duration of follow-up in men and women for most fracture outcomes. There was no evidence of an interaction between falls and BMD for fracture risk. Subsequent risk for a major osteoporotic fracture increased with each additional previous fall in women and men. CONCLUSIONS: A previous self-reported fall confers an increased risk of fracture that is largely independent of BMD. Previous falls should be considered as an additional risk factor in future iterations of FRAX to improve fracture risk prediction.

Adjusting conventional FRAX estimates of fracture probability according to the number of prior falls in the preceding year.

A greater propensity to falling is associated with higher fracture risk. This study provides adjustments to FRAX-based fracture probabilities accounting for the number of prior falls. INTRODUCTION: Prior falls increase subsequent fracture risk but are not currently directly included in the FRAX tool. The aim of this study was to quantify the effect of the number of prior falls on the 10-year probability of fracture determined with FRAX®. METHODS: We studied 21,116 women and men age 40 years or older (mean age 65.7 ± 10.1 years) with fracture probability assessment (FRAX®), self-reported falls for the previous year, and subsequent fracture outcomes in a registry-based cohort. The risks of death, hip fracture, and non-hip major osteoporotic fracture (MOF-NH) were determined by Cox proportional hazards regression for fall number category versus the whole population (i.e., an average number of falls). Ten-year probabilities of hip fracture and major osteoporotic fracture (MOF) were determined according to the number of falls from the hazards of death and fracture incorporated into the FRAX model for the UK. The probability ratios (number of falls vs. average number of falls) provided adjustments to conventional FRAX estimates of fracture probability according to the number of falls. RESULTS: Compared with the average number of falls, the hazard ratios for hip fracture, MOF-NH and death were lower than unity in the absence of a fall history. Hazard ratios increased progressively with an increasing number of reported falls. The probability ratio rose progressively as the number of reported falls increased. Probability ratios decreased with age, an effect that was more marked the greater the number of prior falls. CONCLUSION: The probability ratios provide adjustments to conventional FRAX estimates of fracture probability according to the number of prior falls.

Primary hyperparathyroidism and fracture probability.

The incidence of hip and major osteoporotic fracture was increased in patients with primary hyperparathyroidism even in patients not referred for parathyroidectomy. The risk of death was also increased which attenuated an effect on fracture probabilities. The findings argue for widening the indications for parathyroidectomy in mild primary hyperparathyroidism. INTRODUCTION: Primary hyperparathyroidism (PHPT) is associated with an increase in the risk of fracture. In FRAX, the increase in risk is assumed to be mediated by low bone mineral density (BMD). However, the risk of death is also increased and its effect on fracture probability is not known. OBJECTIVE: The aim of this study was to determine whether PHPT affects hip fracture and major osteoporotic fracture risk independently of bone mineral density (BMD) and whether this and any increase in mortality affects the assessment of fracture probability. METHODS: A register-based survey of patients with PHPT and matched controls in Denmark were identified from hospital registers. The incidence of death, hip fracture, and major osteoporotic fracture were determined for computing fracture probabilities excluding time after parathyroidectomy. The gradient of risk for fracture for differences in BMD was determined in a subset of patients and in BMD controls. The severity of disease was based on serum calcium and parathyroid hormone levels. RESULTS: We identified 6884 patients with biochemically confirmed PHPT and 68,665 matched population controls. On follow-up, excluding time after parathyroidectomy in those undergoing surgery, patients with PHPT had a higher risk of death (+52%), hip fracture (+48%), and major osteoporotic fracture (+36%) than population controls. At any given age, average 10-year probabilities of fracture were higher in patients with PHPT than population controls. The gradient of fracture risk with differences in BMD was similar in cases and controls. Results were similar when confined to patients not undergoing parathyroidectomy. Fracture probability decreased with the severity of disease due to an increase in mortality rather than fracture risk. CONCLUSION: The risk of hip and other major osteoporotic fracture is increased in PHPT irrespective of the disease severity. Fracture probability was attenuated due to the competing effect of mortality. The increased fracture risk in patients treated conservatively argues for widening the indications for parathyroidectomy in mild PHPT.

Adjusting conventional FRAX estimates of fracture probability according to the number of prior fractures.

The risk of a recurrent fragility fracture is high following a first fracture and higher still with more than one prior fracture. This study provides adjustments to FRAX-based fracture probabilities accounting for the number of prior fractures. INTRODUCTION: Prior fractures increase subsequent fracture risk. The aim of this study was to quantify the effect of the number of prior fractures on the 10-year probability of fracture determined with FRAX®. METHODS: The study used data from the Reykjavik Study fracture register that documented prospectively all fractures at all skeletal sites in a large sample of the population of Iceland. Ten-year probabilities of hip fracture and major osteoporotic fracture (MOF) were determined according to the number of prior osteoporotic fractures over a 20-year interval from the hazards of death and fracture. Fracture probabilities were also computed for a prior osteoporotic fracture irrespective of the number of previous fractures. The probability ratios provided adjustments to conventional FRAX estimates of fracture probability according to the number of prior fractures. RESULTS: Probability ratios to adjust 10-year FRAX probabilities of a hip fracture and MOF increased with the number of prior fractures but decreased with age in both men and women. Probability ratios were similar in men and women and for hip fracture and MOF. Mean probability ratios according to the number of prior fractures for all scenarios were 0.95, 1.08, 1.21 and 1.35, for 1,2, 3 and 4 or more prior fractures, respectively. Thus, a simple rule of thumb is to downward adjust FRAX-based fracture probabilities by 5% in the presence of a single prior fracture and to uplift probabilities by 10, 20 and 30% with a history of 2, 3 and 4 or more prior fractures, respectively. CONCLUSION: The probability ratios provide adjustments to conventional FRAX estimates of fracture probability according to the number of prior fractures.

Development of a novel method to measure bone marrow fat fraction in older women using high-resolution peripheral quantitative computed tomography.

Bone marrow adipose tissue (BMAT) has been implicated in a number of conditions associated with bone deterioration and osteoporosis. Several studies have found an inverse relationship between BMAT and bone mineral density (BMD), and higher levels of BMAT in those with prevalent fracture. Magnetic resonance imaging (MRI) is the gold standard for measuring BMAT, but its use is limited by high costs and low availability. We hypothesized that BMAT could also be accurately quantified using high-resolution peripheral quantitative computed tomography (HR-pQCT). METHODS: In the present study, a novel method to quantify the tibia bone marrow fat fraction, defined by MRI, using HR-pQCT was developed. In total, 38 postmenopausal women (mean [standard deviation] age 75.9 [3.1] years) were included and measured at the same site at the distal (n = 38) and ultradistal (n = 18) tibia using both MRI and HR-pQCT. To adjust for partial volume effects, the HR-pQCT images underwent 0 to 10 layers of voxel peeling to remove voxels adjacent to the bone. Linear regression equations were then tested for different degrees of voxel peeling, using the MRI-derived fat fractions as the dependent variable and the HR-pQCT-derived radiodensity as the independent variables. RESULTS: The most optimal HR-pQCT derived model, which applied a minimum of 4 layers of peeled voxel and with more than 1% remaining marrow volume, was able to explain 76% of the variation in the ultradistal tibia bone marrow fat fraction, measured with MRI (p < 0.001). CONCLUSION: The novel HR-pQCT method, developed to estimate BMAT, was able to explain a substantial part of the variation in the bone marrow fat fraction and can be used in future studies investigating the role of BMAT in osteoporosis and fracture prediction.

Menopausal hormone therapy reduces the risk of fracture regardless of falls risk or baseline FRAX probability-results from the Women's Health Initiative hormone therapy trials.

In a combined analysis of 25,389 postmenopausal women aged 50-79 years, enrolled in the two Women's Health Initiative hormone therapy trials, menopausal hormone therapy vs. placebo reduced the risk of fracture regardless of baseline FRAX fracture probability and falls history. INTRODUCTION: The aim of this study was to determine if the anti-fracture efficacy of menopausal hormone therapy (MHT) differed by baseline falls history or fracture risk probability as estimated by FRAX, in a combined analysis of the two Women's Health Initiative (WHI) hormone therapy trials. METHODS: A total of 25,389 postmenopausal women aged 50-79 years were randomized to receive MHT (n = 12,739) or matching placebo (n = 12,650). At baseline, questionnaires were used to collect information on falls history, within the last 12 months, and clinical risk factors. FRAX 10-year probability of major osteoporotic fracture (MOF) was calculated without BMD. Incident clinical fractures were verified using medical records. An extension of Poisson regression was used to investigate the relationship between treatment and fractures in (1) the whole cohort; (2) those with prior falls; and (3) those without prior falls. The effect of baseline FRAX probability on efficacy was investigated in the whole cohort. RESULTS: Over 4.3 ± 2.1 years (mean ± SD), MHT (vs. placebo) significantly reduced the risk of any clinical fracture (hazard ratio [HR] 0.72 [95% CI, 0.65-0.78]), MOF (HR 0.60 [95% CI, 0.53-0.69]), and hip fracture (0.66 [95% CI, 0.45-0.96]). Treatment was effective in reducing the risk of any clinical fracture, MOF, and hip fracture in women regardless of baseline FRAX MOF probability, with no evidence of an interaction between MHT and FRAX (p > 0.30). Similarly, there was no interaction (p > 0.30) between MHT and prior falls. CONCLUSION: In the combined WHI trials, compared to placebo, MHT reduces fracture risk regardless of FRAX probability and falls history in postmenopausal women.

Low-Frequency Synonymous Coding Variation in CYP2R1 Has Large Effects on Vitamin D Levels and Risk of Multiple Sclerosis.

Vitamin D insufficiency is common, correctable, and influenced by genetic factors, and it has been associated with risk of several diseases. We sought to identify low-frequency genetic variants that strongly increase the risk of vitamin D insufficiency and tested their effect on risk of multiple sclerosis, a disease influenced by low vitamin D concentrations. We used whole-genome sequencing data from 2,619 individuals through the UK10K program and deep-imputation data from 39,655 individuals genotyped genome-wide. Meta-analysis of the summary statistics from 19 cohorts identified in CYP2R1 the low-frequency (minor allele frequency = 2.5%) synonymous coding variant g.14900931G>A (p.Asp120Asp) (rs117913124[A]), which conferred a large effect on 25-hydroxyvitamin D (25OHD) levels (-0.43 SD of standardized natural log-transformed 25OHD per A allele; p value = 1.5 × 10-88). The effect on 25OHD was four times larger and independent of the effect of a previously described common variant near CYP2R1. By analyzing 8,711 individuals, we showed that heterozygote carriers of this low-frequency variant have an increased risk of vitamin D insufficiency (odds ratio [OR] = 2.2, 95% confidence interval [CI] = 1.78-2.78, p = 1.26 × 10-12). Individuals carrying one copy of this variant also had increased odds of multiple sclerosis (OR = 1.4, 95% CI = 1.19-1.64, p = 2.63 × 10-5) in a sample of 5,927 case and 5,599 control subjects. In conclusion, we describe a low-frequency CYP2R1 coding variant that exerts the largest effect upon 25OHD levels identified to date in the general European population and implicates vitamin D in the etiology of multiple sclerosis.

A decade of FRAX: how has it changed the management of osteoporosis?

The fracture risk assessment tool, FRAX®, was released in 2008 and provides country-specific algorithms for estimating individualized 10-year probability of hip and major osteoporotic fracture (hip, clinical spine, distal forearm, and proximal humerus). Since its release, 71 models have been made available for 66 countries covering more than 80% of the world population. The website receives approximately 3 million visits annually. Following independent validation, FRAX has been incorporated into more than 80 guidelines worldwide. The application of FRAX in assessment guidelines has been heterogeneous with the adoption of several different approaches in setting intervention thresholds. Whereas most guidelines adopt a case-finding strategy, the case for FRAX-based community screening in the elderly is increasing. The relationship between FRAX and efficacy of intervention has been explored and is expected to influence treatment guidelines in the future.

The gut microbiota is a major regulator of androgen metabolism in intestinal contents.

Androgens exert important effects both in androgen-responsive tissues and in the intestinal tract. To determine the impact of the gut microbiota (GM) on intestinal androgen metabolism, we measured unconjugated (free) and glucuronidated androgen levels in intestinal contents from the small intestine, with a low bacterial density, and from cecum and colon, with a high bacterial density. Using a specific, sensitive gas chromatography-tandem mass spectrometry method, we detected high levels of glucuronidated testosterone (T) and dihydrotestosterone (DHT) in small intestinal content of mice of both sexes, whereas in the distal intestine we observed remarkably high levels of free DHT, exceeding serum levels by >20-fold. Similarly, in young adult men high levels of unconjugated DHT, >70-fold higher than in serum, were detected in feces. In contrast to mice with a normal GM composition, germ-free mice had high levels of glucuronidated T and DHT, but very low free DHT levels, in the distal intestine. These findings demonstrate that the GM is involved in intestinal metabolism and deglucuronidation of DHT and T, resulting in extremely high free levels of the most potent androgen, DHT, in the colonic content of young and healthy mice and men.

Pubertal timing and adult fracture risk in men: A population-based cohort study.

BACKGROUND: Puberty is a critical period for bone mass accrual, and late puberty in boys is associated with reduced bone mass in adult men. The role of variations in pubertal timing within the normal range for adult fracture risk in men is, however, unknown. We, therefore, assessed the association between age at peak height velocity (PHV), an objective measure of pubertal timing, and fracture risk in adult men. METHODS AND FINDINGS: In the BMI Epidemiology Study Gothenburg, 31,971 Swedish men born between January 1, 1945, and December 31, 1961, with detailed growth data (height and weight) available from centrally archived school healthcare records and the conscription register were followed until December 31, 2016. Age at PHV was calculated according to a modified infancy-childhood-puberty model, and fracture information was retrieved from the Swedish National Patient Register. The mean ± SD age at PHV was 14.1 ± 1.1 years. In total, 5,872 men (18.4%) sustained at least 1 fracture after 20 years of age and 5,731 men (17.9%) sustained a non-vertebral fracture after 20 years of age during a mean ± SD follow-up of 37.3 ± 11.7 years. Cox proportional hazards models adjusted for birth year and country of origin revealed that age at PHV was associated with the risk of any fracture and non-vertebral fracture. Participants with age at PHV in the highest tertile (after 14.5 years of age) were at greater risk of any fracture (hazard ratio [HR] 1.15, 95% confidence interval [CI] 1.08-1.22, P < 0.001) and non-vertebral fracture (HR 1.16, 95% CI 1.09-1.24, P < 0.001) compared with those with age at PHV in the lowest tertile (at 13.6 years of age or younger). Additional adjustments for birthweight, childhood BMI, adult educational level, and young adult height did not attenuate the associations between age at PHV and adult fracture risk. Limitations of this study include the inability to adjust for important risk factors for fracture, inadequate power to assess the relation between pubertal timing and specific fracture types, and the limited generalizability to other populations. CONCLUSIONS: In this study, we observed that late pubertal timing was associated with increased adult fracture risk in men. These findings suggest that information on pubertal timing might aid in the identification of those men at greatest risk of fracture.

Increased adipose tissue aromatase activity improves insulin sensitivity and reduces adipose tissue inflammation in male mice.

Females are, in general, more insulin sensitive than males. To investigate whether this is a direct effect of sex-steroids (SS) in white adipose tissue (WAT), we developed a male mouse model overexpressing the aromatase enzyme, converting testosterone (T) to estradiol (E2), specifically in WAT (Ap2-arom mice). Adipose tissue E2 levels were increased while circulating SS levels were unaffected in male Ap2-arom mice. Importantly, male Ap2-arom mice were more insulin sensitive compared with WT mice and exhibited increased serum adiponectin levels and upregulated expression of Glut4 and Irs1 in WAT. The expression of markers of macrophages and immune cell infiltration was markedly decreased in WAT of male Ap2-arom mice. The adipogenesis was enhanced in male Ap2-arom mice, supported by elevated Pparg expression in WAT and enhanced differentiation of preadipocyte into mature adipocytes. In summary, increased adipose tissue aromatase activity reduces adipose tissue inflammation and improves insulin sensitivity in male mice. We propose that estrogen increases insulin sensitivity via a local effect in WAT on adiponectin expression, adipose tissue inflammation, and adipogenesis.

Meta-analysis of genome-wide association studies identifies two loci associated with circulating osteoprotegerin levels.

Osteoprotegerin (OPG) is involved in bone homeostasis and tumor cell survival. Circulating OPG levels are also important biomarkers of various clinical traits, such as cancers and atherosclerosis. OPG levels were measured in serum or in plasma. In a meta-analysis of genome-wide association studies in up to 10 336 individuals from European and Asian origin, we discovered that variants >100 kb upstream of the TNFRSF11B gene encoding OPG and another new locus on chromosome 17q11.2 were significantly associated with OPG variation. We also identified a suggestive locus on chromosome 14q21.2 associated with the trait. Moreover, we estimated that over half of the heritability of OPG levels could be explained by all variants examined in our study. Our findings provide further insight into the genetic regulation of circulating OPG levels.

Hydroxysteroid (17ß)-dehydrogenase 1-deficient female mice present with normal puberty onset but are severely subfertile due to a defect in luteinization and progesterone production.

Hydroxysteroid (17ß)-dehydrogenase type 1 (HSD17B1) catalyzes the conversion of low active 17-ketosteroids, androstenedione (A-dione) and estrone (E1) to highly active 17-hydroxysteroids, testosterone (T) and E2, respectively. In this study, the importance of HSD17B1 in ovarian estrogen production was determined using Hsd17b1 knockout (HSD17B1KO) mice. In these mice, the ovarian HSD17B enzyme activity was markedly reduced, indicating a central role of HSD17B1 in ovarian physiology. The lack of Hsd17b activity resulted in increased ovarian E1:E2 and A-dione:T ratios, but we also observed reduced progesterone concentration in HSD17B1KO ovaries. Accordingly with the altered steroid production, altered expression of Star, Cyp11a1, Lhcgr, Hsd17b7, and especially Cyp17a1 was observed. The ovaries of HSD17B1KO mice presented with all stages of folliculogenesis, while the corpus luteum structure was less defined and number reduced. Surprisingly, bundles of large granular cells of unknown origin appeared in the stroma of the KO ovaries. The HSD17B1KO mice presented with severe subfertility and failed to initiate pseudopregnancy. However, the HSD17B1KO females presented with normal estrous cycle defined by vaginal smears and normal puberty appearance. This study indicates that HSD17B1 is a key enzyme in ovarian steroidogenesis and has a novel function in initiation and stabilization of pregnancy.

Genetic determinants of trabecular and cortical volumetric bone mineral densities and bone microstructure.

Most previous genetic epidemiology studies within the field of osteoporosis have focused on the genetics of the complex trait areal bone mineral density (aBMD), not being able to differentiate genetic determinants of cortical volumetric BMD (vBMD), trabecular vBMD, and bone microstructural traits. The objective of this study was to separately identify genetic determinants of these bone traits as analysed by peripheral quantitative computed tomography (pQCT). Separate GWA meta-analyses for cortical and trabecular vBMDs were performed. The cortical vBMD GWA meta-analysis (n = 5,878) followed by replication (n = 1,052) identified genetic variants in four separate loci reaching genome-wide significance (RANKL, rs1021188, p = 3.6×10⁻¹4; LOC285735, rs271170, p = 2.7×10⁻¹²; OPG, rs7839059, p = 1.2×10⁻¹°; and ESR1/C6orf97, rs6909279, p = 1.1×10⁻9). The trabecular vBMD GWA meta-analysis (n = 2,500) followed by replication (n = 1,022) identified one locus reaching genome-wide significance (FMN2/GREM2, rs9287237, p = 1.9×10⁻9). High-resolution pQCT analyses, giving information about bone microstructure, were available in a subset of the GOOD cohort (n = 729). rs1021188 was significantly associated with cortical porosity while rs9287237 was significantly associated with trabecular bone fraction. The genetic variant in the FMN2/GREM2 locus was associated with fracture risk in the MrOS Sweden cohort (HR per extra T allele 0.75, 95% confidence interval 0.60-0.93) and GREM2 expression in human osteoblasts. In conclusion, five genetic loci associated with trabecular or cortical vBMD were identified. Two of these (FMN2/GREM2 and LOC285735) are novel bone-related loci, while the other three have previously been reported to be associated with aBMD. The genetic variants associated with cortical and trabecular bone parameters differed, underscoring the complexity of the genetics of bone parameters. We propose that a genetic variant in the RANKL locus influences cortical vBMD, at least partly, via effects on cortical porosity, and that a genetic variant in the FMN2/GREM2 locus influences GREM2 expression in osteoblasts and thereby trabecular number and thickness as well as fracture risk.

Meta-analysis of genome-wide scans for total body BMD in children and adults reveals allelic heterogeneity and age-specific effects at the WNT16 locus.

To identify genetic loci influencing bone accrual, we performed a genome-wide association scan for total-body bone mineral density (TB-BMD) variation in 2,660 children of different ethnicities. We discovered variants in 7q31.31 associated with BMD measurements, with the lowest P = 4.1 × 10(-11) observed for rs917727 with minor allele frequency of 0.37. We sought replication for all SNPs located ± 500 kb from rs917727 in 11,052 additional individuals from five independent studies including children and adults, together with de novo genotyping of rs3801387 (in perfect linkage disequilibrium (LD) with rs917727) in 1,014 mothers of children from the discovery cohort. The top signal mapping in the surroundings of WNT16 was replicated across studies with a meta-analysis P = 2.6 × 10(-31) and an effect size explaining between 0.6%-1.8% of TB-BMD variance. Conditional analyses on this signal revealed a secondary signal for total body BMD (P = 1.42 × 10(-10)) for rs4609139 and mapping to C7orf58. We also examined the genomic region for association with skull BMD to test if the associations were independent of skeletal loading. We identified two signals influencing skull BMD variation, including rs917727 (P = 1.9 × 10(-16)) and rs7801723 (P = 8.9 × 10(-28)), also mapping to C7orf58 (r(2) = 0.50 with rs4609139). Wnt16 knockout (KO) mice with reduced total body BMD and gene expression profiles in human bone biopsies support a role of C7orf58 and WNT16 on the BMD phenotypes observed at the human population level. In summary, we detected two independent signals influencing total body and skull BMD variation in children and adults, thus demonstrating the presence of allelic heterogeneity at the WNT16 locus. One of the skull BMD signals mapping to C7orf58 is mostly driven by children, suggesting temporal determination on peak bone mass acquisition. Our life-course approach postulates that these genetic effects influencing peak bone mass accrual may impact the risk of osteoporosis later in life.

A genome-wide association meta-analysis of circulating sex hormone-binding globulin reveals multiple Loci implicated in sex steroid hormone regulation.

Sex hormone-binding globulin (SHBG) is a glycoprotein responsible for the transport and biologic availability of sex steroid hormones, primarily testosterone and estradiol. SHBG has been associated with chronic diseases including type 2 diabetes (T2D) and with hormone-sensitive cancers such as breast and prostate cancer. We performed a genome-wide association study (GWAS) meta-analysis of 21,791 individuals from 10 epidemiologic studies and validated these findings in 7,046 individuals in an additional six studies. We identified twelve genomic regions (SNPs) associated with circulating SHBG concentrations. Loci near the identified SNPs included SHBG (rs12150660, 17p13.1, p = 1.8 × 10(-106)), PRMT6 (rs17496332, 1p13.3, p = 1.4 × 10(-11)), GCKR (rs780093, 2p23.3, p = 2.2 × 10(-16)), ZBTB10 (rs440837, 8q21.13, p = 3.4 × 10(-09)), JMJD1C (rs7910927, 10q21.3, p = 6.1 × 10(-35)), SLCO1B1 (rs4149056, 12p12.1, p = 1.9 × 10(-08)), NR2F2 (rs8023580, 15q26.2, p = 8.3 × 10(-12)), ZNF652 (rs2411984, 17q21.32, p = 3.5 × 10(-14)), TDGF3 (rs1573036, Xq22.3, p = 4.1 × 10(-14)), LHCGR (rs10454142, 2p16.3, p = 1.3 × 10(-07)), BAIAP2L1 (rs3779195, 7q21.3, p = 2.7 × 10(-08)), and UGT2B15 (rs293428, 4q13.2, p = 5.5 × 10(-06)). These genes encompass multiple biologic pathways, including hepatic function, lipid metabolism, carbohydrate metabolism and T2D, androgen and estrogen receptor function, epigenetic effects, and the biology of sex steroid hormone-responsive cancers including breast and prostate cancer. We found evidence of sex-differentiated genetic influences on SHBG. In a sex-specific GWAS, the loci 4q13.2-UGT2B15 was significant in men only (men p = 2.5 × 10(-08), women p = 0.66, heterogeneity p = 0.003). Additionally, three loci showed strong sex-differentiated effects: 17p13.1-SHBG and Xq22.3-TDGF3 were stronger in men, whereas 8q21.12-ZBTB10 was stronger in women. Conditional analyses identified additional signals at the SHBG gene that together almost double the proportion of variance explained at the locus. Using an independent study of 1,129 individuals, all SNPs identified in the overall or sex-differentiated or conditional analyses explained ~15.6% and ~8.4% of the genetic variation of SHBG concentrations in men and women, respectively. The evidence for sex-differentiated effects and allelic heterogeneity highlight the importance of considering these features when estimating complex trait variance.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk.

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)

Eight common genetic variants associated with serum DHEAS levels suggest a key role in ageing mechanisms.

Dehydroepiandrosterone sulphate (DHEAS) is the most abundant circulating steroid secreted by adrenal glands--yet its function is unknown. Its serum concentration declines significantly with increasing age, which has led to speculation that a relative DHEAS deficiency may contribute to the development of common age-related diseases or diminished longevity. We conducted a meta-analysis of genome-wide association data with 14,846 individuals and identified eight independent common SNPs associated with serum DHEAS concentrations. Genes at or near the identified loci include ZKSCAN5 (rs11761528; p = 3.15 × 10(-36)), SULT2A1 (rs2637125; p =  2.61 × 10(-19)), ARPC1A (rs740160; p =  1.56 × 10(-16)), TRIM4 (rs17277546; p =  4.50 × 10(-11)), BMF (rs7181230; p = 5.44 × 10(-11)), HHEX (rs2497306; p =  4.64 × 10(-9)), BCL2L11 (rs6738028; p = 1.72 × 10(-8)), and CYP2C9 (rs2185570; p = 2.29 × 10(-8)). These genes are associated with type 2 diabetes, lymphoma, actin filament assembly, drug and xenobiotic metabolism, and zinc finger proteins. Several SNPs were associated with changes in gene expression levels, and the related genes are connected to biological pathways linking DHEAS with ageing. This study provides much needed insight into the function of DHEAS.