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−14; LOC285735, rs271170, p = 2.7×10−12; OPG, rs7839059, p = 1.2×10−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.

Author Summary

Osteoporosis is a common highly heritable skeletal disease characterized by reduced bone mineral density (BMD) and deteriorated bone microstructure, resulting in an increased risk of fracture. Most previous genetic epidemiology studies have focused on the genetics of the complex trait BMD, not being able to separate genetic determinants of the trabecular and cortical bone compartments and bone microstructure. The trabecular and cortical BMDs can be analysed separately by computed tomography. Therefore, we performed separate genome-wide association studies for trabecular and cortical BMDs, demonstrating that the genetic determinants of cortical and trabecular BMDs differ. Genetic variants in the RANKL, LOC285735, OPG, and ESR1 loci were associated with cortical BMD, while a genetic variant in the FMN2/GREM2 locus was associated with trabecular BMD. Two of these are novel bone-related loci. Follow-up analyses of bone microstructure demonstrated that a genetic variant in the RANKL locus is associated with cortical porosity and that the FMN2/GREM2 locus is associated with trabecular number and thickness. We propose that a genetic variant in the RANKL locus influences cortical BMD via effects on cortical porosity, and that a genetic variant in the FMN2/GREM2 locus influences trabecular BMD and fracture risk via effects on both trabecular number and thickness.

A mendelian randomization study based on data from multiple cohorts conducted by Karani Santhanakrishnan Vimaleswaran and colleagues re-examines the causal nature of the relationship between vitamin D levels and obesity.

Background

Obesity is associated with vitamin D deficiency, and both are areas of active public health concern. We explored the causality and direction of the relationship between body mass index (BMI) and 25-hydroxyvitamin D [25(OH)D] using genetic markers as instrumental variables (IVs) in bi-directional Mendelian randomization (MR) analysis.

Methods and Findings

We used information from 21 adult cohorts (up to 42,024 participants) with 12 BMI-related SNPs (combined in an allelic score) to produce an instrument for BMI and four SNPs associated with 25(OH)D (combined in two allelic scores, separately for genes encoding its synthesis or metabolism) as an instrument for vitamin D. Regression estimates for the IVs (allele scores) were generated within-study and pooled by meta-analysis to generate summary effects.

Associations between vitamin D scores and BMI were confirmed in the Genetic Investigation of Anthropometric Traits (GIANT) consortium (n = 123,864). Each 1 kg/m2 higher BMI was associated with 1.15% lower 25(OH)D (p = 6.52×10−27). The BMI allele score was associated both with BMI (p = 6.30×10−62) and 25(OH)D (−0.06% [95% CI −0.10 to −0.02], p = 0.004) in the cohorts that underwent meta-analysis. The two vitamin D allele scores were strongly associated with 25(OH)D (p≤8.07×10−57 for both scores) but not with BMI (synthesis score, p = 0.88; metabolism score, p = 0.08) in the meta-analysis. A 10% higher genetically instrumented BMI was associated with 4.2% lower 25(OH)D concentrations (IV ratio: −4.2 [95% CI −7.1 to −1.3], p = 0.005). No association was seen for genetically instrumented 25(OH)D with BMI, a finding that was confirmed using data from the GIANT consortium (p≥0.57 for both vitamin D scores).

Conclusions

On the basis of a bi-directional genetic approach that limits confounding, our study suggests that a higher BMI leads to lower 25(OH)D, while any effects of lower 25(OH)D increasing BMI are likely to be small. Population level interventions to reduce BMI are expected to decrease the prevalence of vitamin D deficiency.

Please see later in the article for the Editors' Summary

Editors' Summary

Background

Obesity—having an unhealthy amount of body fat—is increasing worldwide. In the US, for example, a third of the adult population is now obese. Obesity is defined as having a body mass index (BMI, an indicator of body fat calculated by dividing a person's weight in kilograms by their height in meters squared) of more than 30.0 kg/m2. Although there is a genetic contribution to obesity, people generally become obese by consuming food and drink that contains more energy than they need for their daily activities. Thus, obesity can be prevented by having a healthy diet and exercising regularly. Compared to people with a healthy weight, obese individuals have an increased risk of developing diabetes, heart disease and stroke, and tend to die younger. They also have a higher risk of vitamin D deficiency, another increasingly common public health concern. Vitamin D, which is essential for healthy bones as well as other functions, is made in the skin after exposure to sunlight but can also be obtained through the diet and through supplements.

Why Was This Study Done?

Observational studies cannot prove that obesity causes vitamin D deficiency because obese individuals may share other characteristics that reduce their circulating 25-hydroxy vitamin D [25(OH)D] levels (referred to as confounding). Moreover, observational studies cannot indicate whether the larger vitamin D storage capacity of obese individuals (vitamin D is stored in fatty tissues) lowers their 25(OH)D levels or whether 25(OH)D levels influence fat accumulation (reverse causation). If obesity causes vitamin D deficiency, monitoring and treating vitamin D deficiency might alleviate some of the adverse health effects of obesity. Conversely, if low vitamin D levels cause obesity, encouraging people to take vitamin D supplements might help to control the obesity epidemic. Here, the researchers use bi-directional “Mendelian randomization” to examine the direction and causality of the relationship between BMI and 25(OH)D. In Mendelian randomization, causality is inferred from associations between genetic variants that mimic the influence of a modifiable environmental exposure and the outcome of interest. Because gene variants do not change over time and are inherited randomly, they are not prone to confounding and are free from reverse causation. Thus, if a lower vitamin D status leads to obesity, genetic variants associated with lower 25(OH)D concentrations should be associated with higher BMI, and if obesity leads to a lower vitamin D status, then genetic variants associated with higher BMI should be associated with lower 25(OH)D concentrations.

What Did the Researchers Do and Find?

The researchers created a “BMI allele score” based on 12 BMI-related gene variants and two “25(OH)D allele scores,” which are based on gene variants that affect either 25(OH)D synthesis or breakdown. Using information on up to 42,024 participants from 21 studies, the researchers showed that the BMI allele score was associated with both BMI and with 25(OH)D levels among the study participants. Based on this information, they calculated that each 10% increase in BMI will lead to a 4.2% decrease in 25(OH)D concentrations. By contrast, although both 25(OH)D allele scores were strongly associated with 25(OH)D levels, neither score was associated with BMI. This lack of an association between 25(OH)D allele scores and obesity was confirmed using data from more than 100,000 individuals involved in 46 studies that has been collected by the GIANT (Genetic Investigation of Anthropometric Traits) consortium.

What Do These Findings Mean?

These findings suggest that a higher BMI leads to a lower vitamin D status whereas any effects of low vitamin D status on BMI are likely to be small. That is, these findings provide evidence for obesity as a causal factor in the development of vitamin D deficiency but not for vitamin D deficiency as a causal factor in the development of obesity. These findings suggest that population-level interventions to reduce obesity should lead to a reduction in the prevalence of vitamin D deficiency and highlight the importance of monitoring and treating vitamin D deficiency as a means of alleviating the adverse influences of obesity on health.

Additional Information

Please access these Web sites via the online version of this summary at http://dx.doi.org/10.1371/journal.pmed.1001383.

The US Centers for Disease Control and Prevention provides information on all aspects of overweight and obesity (in English and Spanish); a data brief provides information about the vitamin D status of the US population

The World Health Organization provides information on obesity (in several languages)

The UK National Health Service Choices website provides detailed information about obesity and a link to a personal story about losing weight; it also provides information about vitamin D

The International Obesity Taskforce provides information about the global obesity epidemic

The US Department of Agriculture's ChooseMyPlate.gov website provides a personal healthy eating plan; the Weight-control Information Network is an information service provided for the general public and health professionals by the US National Institute of Diabetes and Digestive and Kidney Diseases (in English and Spanish)

The US Office of Dietary Supplements provides information about vitamin D (in English and Spanish)

MedlinePlus has links to further information about obesity and about vitamin D (in English and Spanish)

Wikipedia has a page on Mendelian randomization (note: Wikipedia is a free online encyclopedia that anyone can edit; available in several languages)

Overview and details of the collaborative large-scale genetic association study (D-CarDia) provide information about vitamin D and the risk of cardiovascular disease, diabetes and related traits

Clinical relevance of a genetic predisposition to elevated blood pressure was quantified during the transition from childhood to adulthood in a population-based Finnish cohort (N=2,357). Blood pressure was measured at baseline in 1980 (age 3–18 years) and in follow-ups in 1983, 1986, 2001 and 2007. Thirteen single nucleotide polymorphisms associated with blood pressure were genotyped and three genetic risk scores associated with systolic and diastolic blood pressure and their combination were derived for all participants. Effects of the genetic risk score were 0.47 mmHg for systolic and 0.53 mmHg for diastolic blood pressure (both p<0.01). The combination genetic risk score was associated with diastolic blood pressure from age 9 onwards (β=0.68 mmHg, p=0.015). Replications in 1194 participants of the Bogalusa Heart Study showed essentially similar results. The participants in the highest quintile of the combination genetic risk score had a 1.82-fold risk of hypertension in adulthood (p<0.0001) compared with the lowest quintile, independent of a family history of premature hypertension. These findings show that genetic variants are associated with preclinical blood pressure traits in childhood, individuals with several susceptibility alleles have on average a 0.5 mmHg higher blood pressure and this trajectory continues from childhood to adulthood.

Depressive mood is often preceded by sleep problems, suggesting that they increase the risk of depression. Sleep problems can also reflect prodromal symptom of depression, thus temporal precedence alone is insufficient to confirm causality. The authors applied recently introduced statistical causal-discovery algorithms that can estimate causality from cross-sectional samples in order to infer the direction of causality between the two sets of symptoms from a novel perspective. Two common-population samples were used; one from the Young Finns study (690 men and 997 women, average age 37.7 years, range 30–45), and another from the Wisconsin Longitudinal study (3101 men and 3539 women, average age 53.1 years, range 52–55). These included three depression questionnaires (two in Young Finns data) and two sleep problem questionnaires. Three different causality estimates were constructed for each data set, tested in a benchmark data with a (practically) known causality, and tested for assumption violations using simulated data. Causality algorithms performed well in the benchmark data and simulations, and a prediction was drawn for future empirical studies to confirm: for minor depression/dysphoria, sleep problems cause significantly more dysphoria than dysphoria causes sleep problems. The situation may change as depression becomes more severe, or more severe levels of symptoms are evaluated; also, artefacts due to severe depression being less well presented in the population data than minor depression may intervene the estimation for depression scales that emphasize severe symptoms. The findings are consistent with other emerging epidemiological and biological evidence.

Background

Concern has been recently raised about possible adverse cardio-metabolic effects of high selenium status, such as increased risks of diabetes and hyperlipidemia. However, most of the evidence comes from selenium-replete populations such as the US.

Objectives

To examine cross-sectional and longitudinal associations of serum selenium with cardiovascular risk factors in Finland where selenium levels were among the lowest in the world until the early 1980s before the implementation of a nationwide selenium fertilization program.

Methods

Serum selenium was measured in 1,235 young Finns aged 3-18 years at baseline in 1980 (pre-fertilization), and in a subgroup (N=262) at the 6-year follow-up (1986, post-fertilization). During the 27-year follow-up, serum lipids, blood pressure, BMI, and smoking were assessed five times (1980, 1983, 1986, 2001, and 2007).

Results

Mean (±SD) serum selenium concentrations were 74.3±14.0 ng/mL in 1980 and 106.6±12.5 ng/mL in 1986 (average increase 32.3 ng/mL; 95% CI: 30.3 to 34.3, p<0.0001). In univariate and multivariable cross-sectional models in 1980 and 1986, increased serum selenium levels were consistently associated with increased total, HDL- and LDL-cholesterol. However, the average longitudinal changes in lipids were −0.20 mmol/L (95% CI: −0.30 to −0.10, p<0.0001) for total cholesterol, 0.06 mmol/L (95% CI: 0.03 to 0.10, p<0.0001) for HDL-cholesterol, and −0.23 mmol/L (95% CI: −0.31 to −0.14, p<0.0001) for LDL-cholesterol. Selenium measured in 1986 was not associated with lipids assessed in 2001 and 2007.

Conclusions

Cross-sectional findings from the Young Finns study corroborate positive associations of selenium status with serum lipids. However, longitudinal evidence does not support the causality of this link.

Background

Genome-wide association (GWA) studies have identified several susceptibility loci for metabolic syndrome (MetS) component traits, but have had variable success in identifying susceptibility loci to the syndrome as an entity. We conducted a GWA study on MetS and its component traits in four Finnish cohorts consisting of 2637 MetS cases and 7927 controls, both free of diabetes, and followed the top loci in an independent sample with transcriptome and NMR-based metabonomics data. Furthermore, we tested for loci associated with multiple MetS component traits using factor analysis and built a genetic risk score for MetS.

Methods and Results

A previously known lipid locus, APOA1/C3/A4/A5 gene cluster region (SNP rs964184), was associated with MetS in all four study samples (P=7.23×10−9 in meta-analysis). The association was further supported by serum metabolite analysis, where rs964184 associated with various VLDL, TG, and HDL metabolites (P=0.024-1.88×10−5). Twenty-two previously identified susceptibility loci for individual MetS component traits were replicated in our GWA and factor analysis. Most of these associated with lipid phenotypes and none with two or more uncorrelated MetS components. A genetic risk score, calculated as the number of alleles in loci associated with individual MetS traits, was strongly associated with MetS status.

Conclusions

Our findings suggest that genes from lipid metabolism pathways have the key role in the genetic background of MetS. We found little evidence for pleiotropy linking dyslipidemia and obesity to the other MetS component traits such as hypertension and glucose intolerance.

Rationale

Asthma has substantial morbidity and mortality and a strong genetic component, but identification of genetic risk factors is limited by availability of suitable studies.

Objectives

To test if population-based cohorts with self-reported physician-diagnosed asthma and genome-wide association (GWA) data could be used to validate known associations with asthma and identify novel associations.

Methods

The APCAT (Analysis in Population-based Cohorts of Asthma Traits) consortium consists of 1,716 individuals with asthma and 16,888 healthy controls from six European-descent population-based cohorts. We examined associations in APCAT of thirteen variants previously reported as genome-wide significant (P<5x10−8) and three variants reported as suggestive (P<5×10−7). We also searched for novel associations in APCAT (Stage 1) and followed-up the most promising variants in 4,035 asthmatics and 11,251 healthy controls (Stage 2). Finally, we conducted the first genome-wide screen for interactions with smoking or hay fever.

Main Results

We observed association in the same direction for all thirteen previously reported variants and nominally replicated ten of them. One variant that was previously suggestive, rs11071559 in RORA, now reaches genome-wide significance when combined with our data (P = 2.4×10−9). We also identified two genome-wide significant associations: rs13408661 near IL1RL1/IL18R1 (PStage1+Stage2 = 1.1x10−9), which is correlated with a variant recently shown to be associated with asthma (rs3771180), and rs9268516 in the HLA region (PStage1+Stage2 = 1.1x10−8), which appears to be independent of previously reported associations in this locus. Finally, we found no strong evidence for gene-environment interactions with smoking or hay fever status.

Conclusions

Population-based cohorts with simple asthma phenotypes represent a valuable and largely untapped resource for genetic studies of asthma.

Background

Pregnancy-related cardiovascular conditions are associated with both poorer pregnancy outcomes and cardiovascular disease later in life. Little is known about the relationship between preconception cardiovascular risk factor levels and pregnancy complications.

Methods

Data from the Cardiovascular Risk in Young Finns Study were linked with birth registry data for 1142 primiparous women. Age-standardized levels of total cholesterol, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, triglycerides, blood pressure, insulin, and glucose from the study visit prior to last menstrual period were calculated. These factors were examined as predictors of gestational age, preterm birth (<37 weeks), birthweight, low birthweight (<2500 g), small-for-gestational-age (weight <10th percentile for gestational age), hypertensive disorders of pregnancy, and gestational diabetes, using linear and Poisson regression with adjustment for age, body mass index, smoking, and socioeconomic status.

Results

Higher triglycerides were associated with a higher risk of hypertensive disorders (adjusted risk ratio [aRR]= 1.42 [95% confidence interval (CI) = 0.90–2.23]), pre-eclampsia (1.70 [1.08–2.65]), and gestational diabetes (1.68 [1.25–2.25]). After removing women with pregnancy complications (n=30), the estimated aRR for the association between systolic blood pressure and preterm birth was 1.23 (95% CI= 0.99–1.54); for HDL-c and low birthweight, 0.97 (0.73–1.28); for diastolic blood pressure and small-for-gestational-age, 0.98 (0.81–1.20); and for systolic blood pressure and small-for-gestational-age, 1.18 (0.97–1.45).

Conclusions

High lipid levels before pregnancy predict an increased risk of pre-eclampsia and gestational diabetes. Reported associations between these pregnancy complications and later cardiovascular disease of the mother are probably explained, at least in part, by maternal conditions that precede pregnancy. Interventions to improve cardiovascular health before pregnancy may reduce risk of pregnancy complications.

Association testing of multiple correlated phenotypes offers better power than univariate analysis of single traits. We analyzed 6,600 individuals from two population-based cohorts with both genome-wide SNP data and serum metabolomic profiles. From the observed correlation structure of 130 metabolites measured by nuclear magnetic resonance, we identified 11 metabolic networks and performed a multivariate genome-wide association analysis. We identified 34 genomic loci at genome-wide significance, of which 7 are novel. In comparison to univariate tests, multivariate association analysis identified nearly twice as many significant associations in total. Multi-tissue gene expression studies identified variants in our top loci, SERPINA1 and AQP9, as eQTLs and showed that SERPINA1 and AQP9 expression in human blood was associated with metabolites from their corresponding metabolic networks. Finally, liver expression of AQP9 was associated with atherosclerotic lesion area in mice, and in human arterial tissue both SERPINA1 and AQP9 were shown to be upregulated (6.3-fold and 4.6-fold, respectively) in atherosclerotic plaques. Our study illustrates the power of multi-phenotype GWAS and highlights candidate genes for atherosclerosis.

Author Summary

In this study, we aim to identify novel genetic variants for metabolism, characterize their effects on nearby genes, and show that the nearby genes are associated with metabolism and atherosclerosis. To discover new genetic variants, we use an alternative approach to traditional genome-wide association studies: we leverage the information in phenotype covariance to increase our statistical power. We identify variants at seven novel loci and then show that our top signals drive expression of nearby genes AQP9 and SERPINA1 in multiple tissues. We demonstrate that AQP9 and SERPINA1 gene expression, in turn, is associated with metabolite levels. Finally, we show that the genes are associated with atherosclerosis using mouse atherosclerotic lesion size (AQP9) as well as tissue from healthy human arteries and atherosclerotic plaques (AQP9 and SERPINA1). This study illustrates that multivariate analysis of correlated metabolites can boost power for gene discovery substantially. Further functional work will need to be performed to elucidate the biological role of SERPINA1 and AQP9 in atherosclerosis.

Stature is a classical and highly heritable complex trait, with 80%–90% of variation explained by genetic factors. In recent years, genome-wide association studies (GWAS) have successfully identified many common additive variants influencing human height; however, little attention has been given to the potential role of recessive genetic effects. Here, we investigated genome-wide recessive effects by an analysis of inbreeding depression on adult height in over 35,000 people from 21 different population samples. We found a highly significant inverse association between height and genome-wide homozygosity, equivalent to a height reduction of up to 3 cm in the offspring of first cousins compared with the offspring of unrelated individuals, an effect which remained after controlling for the effects of socio-economic status, an important confounder (χ2 = 83.89, df = 1; p = 5.2×10−20). There was, however, a high degree of heterogeneity among populations: whereas the direction of the effect was consistent across most population samples, the effect size differed significantly among populations. It is likely that this reflects true biological heterogeneity: whether or not an effect can be observed will depend on both the variance in homozygosity in the population and the chance inheritance of individual recessive genotypes. These results predict that multiple, rare, recessive variants influence human height. Although this exploratory work focuses on height alone, the methodology developed is generally applicable to heritable quantitative traits (QT), paving the way for an investigation into inbreeding effects, and therefore genetic architecture, on a range of QT of biomedical importance.

Author Summary

Studies investigating the extent to which genetics influences human characteristics such as height have concentrated mainly on common variants of genes, where having one or two copies of a given variant influences the trait or risk of disease. This study explores whether a different type of genetic variant might also be important. We investigate the role of recessive genetic variants, where two identical copies of a variant are required to have an effect. By measuring genome-wide homozygosity—the phenomenon of inheriting two identical copies at a given point of the genome—in 35,000 individuals from 21 European populations, and by comparing this to individual height, we found that the more homozygous the genome, the shorter the individual. The offspring of first cousins (who have increased homozygosity) were predicted to be up to 3 cm shorter on average than the offspring of unrelated parents. Height is influenced by the combined effect of many recessive variants dispersed across the genome. This may also be true for other human characteristics and diseases, opening up a new way to understand how genetic variation influences our health.

Background

The object of this study was to identify temperament patterns in the Finnish population, and to determine the relationship between these profiles and life habits, socioeconomic status, and health.

Methods/Principal Findings

A cluster analysis of the Temperament and Character Inventory subscales was performed on 3,761 individuals from the Northern Finland Birth Cohort 1966 and replicated on 2,097 individuals from the Cardiovascular Risk in Young Finns study. Clusters were formed using the k-means method and their relationship with 115 variables from the areas of life habits, socioeconomic status and health was examined.

Results

Four clusters were identified for both genders. Individuals from Cluster I are characterized by high persistence, low extravagance and disorderliness. They have healthy life habits, and lowest scores in most of the measures for psychiatric disorders. Cluster II individuals are characterized by low harm avoidance and high novelty seeking. They report the best physical capacity and highest level of income, but also high rate of divorce, smoking, and alcohol consumption. Individuals from Cluster III are not characterized by any extreme characteristic. Individuals from Cluster IV are characterized by high levels of harm avoidance, low levels of exploratory excitability and attachment, and score the lowest in most measures of health and well-being.

Conclusions

This study shows that the temperament subscales do not distribute randomly but have an endogenous structure, and that these patterns have strong associations to health, life events, and well-being.

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

Author Summary

Bone traits are highly dependent on genetic factors. To date, numerous genetic loci for bone mineral density (BMD) and only one locus for osteoporotic fracture have been previously identified to be genome-wide significant. Cortical bone has been reported to be an important determinant of bone strength; so far, no genome-wide association studies (GWAS) have been performed for cortical bone thickness (CBT) of the tibial and radial diaphysis or BMD at forearm, a skeletal site rich in cortical bone. Therefore, we performed two separated meta-analyses of GWAS for cortical thickness of the tibia in 3 independent cohorts of 5,878 men and women, and for forearm BMD in 5 cohorts of 5,672 individuals. We identified the 7q31 locus, which contains WNT16, to be associated with CBT and BMD. Four SNPs from this locus were then tested in 2,023 osteoporotic fracture cases and 3,740 controls. One of these SNPs was genome-wide significant, and two were genome-wide suggestive, for forearm fracture. Generating a mouse with targeted disruption of Wnt16, we also demonstrated that mice lacking this protein had substantially thinner bone cortices and reduced bone strength than their wild-type littermates. These findings highlight WNT16 as a clinically relevant member of the Wnt signaling pathway and increase our understanding of the etiology of osteoporosis-related phenotypes and fracture.