M. Carola Zillikens

Twenty bone-mineral-density loci identified by large-scale meta-analysis of genome-wide association studies

Bone mineral density (BMD) is a heritable complex trait used in the clinical diagnosis of osteoporosis and the assessment of fracture risk. We performed meta-analysis of five genome-wide association studies of femoral neck and lumbar spine BMD in 19,195 subjects of Northern European descent. We identified 20 BMD loci that reached genome-wide significance (GWS; P < 5 × 10−8), of which 13 map to regions not previously associated with this trait: 1p31.3 (GPR177), 2p21 (SPTBN1), 3p22 (CTNNB1), 4q21.1 (MEPE), 5q14 (MEF2C), 7p14 (STARD3NL), 7q21.3 (FLJ42280), 11p11.2 (LRP4, ARHGAP1, F2), 11p14.1 (DCDC5), 11p15 (SOX6), 16q24 (FOXL1), 17q21 (HDAC5) and 17q12 (CRHR1). The meta-analysis also confirmed at GWS level seven known BMD loci on 1p36 (ZBTB40), 6q25 (ESR1), 8q24 (TNFRSF11B), 11q13.4 (LRP5), 12q13 (SP7), 13q14 (TNFSF11) and 18q21 (TNFRSF11A). The many SNPs associated with BMD map to genes in signaling pathways with relevance to bone metabolism and highlight the complex genetic architecture that underlies osteoporosis and variation in BMD.

NRXN3 Is a Novel Locus for Waist Circumference: A Genome-Wide Association Study from the CHARGE Consortium

Central abdominal fat is a strong risk factor for diabetes and cardiovascular disease. To identify common variants influencing central abdominal fat, we conducted a two-stage genome-wide association analysis for waist circumference (WC). In total, three loci reached genome-wide significance. In stage 1, 31,373 individuals of Caucasian descent from eight cohort studies confirmed the role of FTO and MC4R and identified one novel locus associated with WC in the neurexin 3 gene [NRXN3 (rs10146997, p = 6.4×10−7)]. The association with NRXN3 was confirmed in stage 2 by combining stage 1 results with those from 38,641 participants in the GIANT consortium (p = 0.009 in GIANT only, p = 5.3×10−8 for combined analysis, n = 70,014). Mean WC increase per copy of the G allele was 0.0498 z-score units (0.65 cm). This SNP was also associated with body mass index (BMI) [p = 7.4×10−6, 0.024 z-score units (0.10 kg/m2) per copy of the G allele] and the risk of obesity (odds ratio 1.13, 95% CI 1.07–1.19; p = 3.2×10−5 per copy of the G allele). The NRXN3 gene has been previously implicated in addiction and reward behavior, lending further evidence that common forms of obesity may be a central nervous system-mediated disorder. Our findings establish that common variants in NRXN3 are associated with WC, BMI, and obesity.

A meta-analysis of the association of fracture risk and body mass index in women.

Several recent studies suggest that obesity may be a risk factor for fracture. The aim of this study was to investigate the association between body mass index (BMI) and future fracture risk at different skeletal sites. In prospective cohorts from more than 25 countries, baseline data on BMI were available in 398,610 women with an average age of 63 (range, 20–105) years and follow up of 2.2 million person‐years during which 30,280 osteoporotic fractures (6457 hip fractures) occurred. Femoral neck BMD was measured in 108,267 of these women. Obesity (BMI ≥ 30 kg/m2) was present in 22%. A majority of osteoporotic fractures (81%) and hip fractures (87%) arose in non‐obese women. Compared to a BMI of 25 kg/m2, the hazard ratio (HR) for osteoporotic fracture at a BMI of 35 kg/m2 was 0.87 (95% confidence interval [CI], 0.85–0.90). When adjusted for bone mineral density (BMD), however, the same comparison showed that the HR for osteoporotic fracture was increased (HR, 1.16; 95% CI, 1.09–1.23). Low BMI is a risk factor for hip and all osteoporotic fracture, but is a protective factor for lower leg fracture, whereas high BMI is a risk factor for upper arm (humerus and elbow) fracture. When adjusted for BMD, low BMI remained a risk factor for hip fracture but was protective for osteoporotic fracture, tibia and fibula fracture, distal forearm fracture, and upper arm fracture. When adjusted for BMD, high BMI remained a risk factor for upper arm fracture but was also a risk factor for all osteoporotic fractures. The association between BMI and fracture risk is complex, differs across skeletal sites, and is modified by the interaction between BMI and BMD. At a population level, high BMI remains a protective factor for most sites of fragility fracture. The contribution of increasing population rates of obesity to apparent decreases in fracture rates should be explored.

An Integration of Genome-Wide Association Study and Gene Expression Profiling to Prioritize the Discovery of Novel Susceptibility Loci for Osteoporosis-Related Traits

Osteoporosis is a complex disorder and commonly leads to fractures in elderly persons. Genome-wide association studies (GWAS) have become an unbiased approach to identify variations in the genome that potentially affect health. However, the genetic variants identified so far only explain a small proportion of the heritability for complex traits. Due to the modest genetic effect size and inadequate power, true association signals may not be revealed based on a stringent genome-wide significance threshold. Here, we take advantage of SNP and transcript arrays and integrate GWAS and expression signature profiling relevant to the skeletal system in cellular and animal models to prioritize the discovery of novel candidate genes for osteoporosis-related traits, including bone mineral density (BMD) at the lumbar spine (LS) and femoral neck (FN), as well as geometric indices of the hip (femoral neck-shaft angle, NSA; femoral neck length, NL; and narrow-neck width, NW). A two-stage meta-analysis of GWAS from 7,633 Caucasian women and 3,657 men, revealed three novel loci associated with osteoporosis-related traits, including chromosome 1p13.2 (RAP1A, p = 3.6×10−8), 2q11.2 (TBC1D8), and 18q11.2 (OSBPL1A), and confirmed a previously reported region near TNFRSF11B/OPG gene. We also prioritized 16 suggestive genome-wide significant candidate genes based on their potential involvement in skeletal metabolism. Among them, 3 candidate genes were associated with BMD in women. Notably, 2 out of these 3 genes (GPR177, p = 2.6×10−13; SOX6, p = 6.4×10−10) associated with BMD in women have been successfully replicated in a large-scale meta-analysis of BMD, but none of the non-prioritized candidates (associated with BMD) did. Our results support the concept of our prioritization strategy. In the absence of direct biological support for identified genes, we highlighted the efficiency of subsequent functional characterization using publicly available expression profiling relevant to the skeletal system in cellular or whole animal models to prioritize candidate genes for further functional validation.

Mild hyponatremia as a risk factor for fractures: the Rotterdam Study.

Recent studies suggest that mild hyponatremia is associated with fractures, but prospective studies are lacking. We studied whether hyponatremia is associated with fractures, falls, and/or bone mineral density (BMD). A total of 5208 elderly subjects with serum sodium assessed at baseline were included from the prospective population‐based Rotterdam Study. The following data were analyzed: BMD, vertebral fractures (mean follow‐up 6.4 years), nonvertebral fractures (7.4 years), recent falls, comorbidity, medication, and mortality. Hyponatremia was detected in 399 subjects (7.7%, 133.4 ± 2.0 mmol/L). Subjects with hyponatremia were older (73.5 ± 10.3 years versus 70.0 ± 9.0 years, p < .001), had more recent falls (23.8% versus16.4%, p < .01), higher type 2 diabetes mellitus prevalence (22.2% versus 10.3%, p < .001), and more often used diuretics (31.1% versus 15.0%, p < .001). Hyponatremia was not associated with lower BMD but was associated with increased risk of incident nonvertebral fractures [hazard ratio (HR) =1.39, 95% confidence interval (CI) 1.11–1.73, p = .004] after adjustment for age, sex, and body mass index. Further adjustments for disability index, use of diuretics, use of psycholeptics, recent falls, and diabetes did not modify results. In the fully adjusted model, subjects with hyponatremia also had increased risk of vertebral fractures at baseline [odds ratio (OR) = 1.78, 95% CI 1.04–3.06, p = .037] but not at follow‐up. Finally, all‐cause mortality was higher in subjects with hyponatremia (HR = 1.21, 95% CI 1.03–1.43, p = .022). It is concluded that mild hyponatremia in the elderly is associated with an increased risk of vertebral fractures and incident nonvertebral fractures but not with BMD. Increased fracture risk in hyponatremia also was independent of recent falls, pointing toward a possible effect on bone quality.

Femoral neck BMD is a strong predictor of hip fracture susceptibility in elderly men and women because it detects cortical bone instability: The Rotterdam Study

We studied HSA measurements in relation to hip fracture risk in 4806 individuals (2740 women). Hip fractures (n = 147) occurred at the same absolute levels of bone instability in both sexes. Cortical instability (propensity of thinner cortices in wide diameters to buckle) explains why hip fracture risk at different BMD levels is the same across sexes.

High Bone Mineral Density and Fracture Risk in Type 2 Diabetes as Skeletal Complications of Inadequate Glucose Control: The Rotterdam Study

OBJECTIVE Individuals with type 2 diabetes have increased fracture risk despite higher bone mineral density (BMD). Our aim was to examine the influence of glucose control on skeletal complications. RESEARCH DESIGN AND METHODS Data of 4,135 participants of the Rotterdam Study, a prospective population-based cohort, were available (mean follow-up 12.2 years). At baseline, 420 participants with type 2 diabetes were classified by glucose control (according to HbA1c calculated from fructosamine), resulting in three comparison groups: adequately controlled diabetes (ACD; n = 203; HbA1c <7.5%), inadequately controlled diabetes (ICD; n = 217; HbA1c ≥7.5%), and no diabetes (n = 3,715). Models adjusted for sex, age, height, and weight (and femoral neck BMD) were used to test for differences in bone parameters and fracture risk (hazard ratio [HR] [95% CI]). RESULTS The ICD group had 1.1–5.6% higher BMD, 4.6–5.6% thicker cortices, and −1.2 to −1.8% narrower femoral necks than ACD and ND, respectively. Participants with ICD had 47–62% higher fracture risk than individuals without diabetes (HR 1.47 [1.12–1.92]) and ACD (1.62 [1.09–2.40]), whereas those with ACD had a risk similar to those without diabetes (0.91 [0.67–1.23]). CONCLUSIONS Poor glycemic control in type 2 diabetes is associated with fracture risk, high BMD, and thicker femoral cortices in narrower bones. We postulate that fragility in apparently “strong” bones in ICD can result from microcrack accumulation and/or cortical porosity, reflecting impaired bone repair.

Genome-wide meta-analysis of observational studies shows common genetic variants associated with macronutrient intake

Background: Macronutrient intake varies substantially between individuals, and there is evidence that this variation is partly accounted for by genetic variants. Objective: The objective of the study was to identify common genetic variants that are associated with macronutrient intake. Design: We performed 2-stage genome-wide association (GWA) meta-analysis of macronutrient intake in populations of European descent. Macronutrients were assessed by using food-frequency questionnaires and analyzed as percentages of total energy consumption from total fat, protein, and carbohydrate. From the discovery GWA (n = 38,360), 35 independent loci associated with macronutrient intake at P < 5 × 10−6 were identified and taken forward to replication in 3 additional cohorts (n = 33,533) from the DietGen Consortium. For one locus, fat mass obesity-associated protein (FTO), cohorts with Illumina MetaboChip genotype data (n = 7724) provided additional replication data. Results: A variant in the chromosome 19 locus (rs838145) was associated with higher carbohydrate (β ± SE: 0.25 ± 0.04%; P = 1.68 × 10−8) and lower fat (β ± SE: −0.21 ± 0.04%; P = 1.57 × 10−9) consumption. A candidate gene in this region, fibroblast growth factor 21 (FGF21), encodes a fibroblast growth factor involved in glucose and lipid metabolism. The variants in this locus were associated with circulating FGF21 protein concentrations (P < 0.05) but not mRNA concentrations in blood or brain. The body mass index (BMI)–increasing allele of the FTO variant (rs1421085) was associated with higher protein intake (β ± SE: 0.10 ± 0.02%; P = 9.96 × 10−10), independent of BMI (after adjustment for BMI, β ± SE: 0.08 ± 0.02%; P = 3.15 × 10−7). Conclusion: Our results indicate that variants in genes involved in nutrient metabolism and obesity are associated with macronutrient consumption in humans. Trials related to this study were registered at clinicaltrials.gov as NCT00005131 (Atherosclerosis Risk in Communities), NCT00005133 (Cardiovascular Health Study), NCT00005136 (Family Heart Study), NCT00005121 (Framingham Heart Study), NCT00083369 (Genetic and Environmental Determinants of Triglycerides), NCT01331512 (InCHIANTI Study), and NCT00005487 (Multi-Ethnic Study of Atherosclerosis).

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 (r2 = 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 study of aging

Human longevity and healthy aging show moderate heritability (20%-50%). We conducted a meta-analysis of genome-wide association studies from 9 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium for 2 outcomes: (1) all-cause mortality, and (2) survival free of major disease or death. No single nucleotide polymorphism (SNP) was a genome-wide significant predictor of either outcome (p < 5 × 10(-8)). We found 14 independent SNPs that predicted risk of death, and 8 SNPs that predicted event-free survival (p < 10(-5)). These SNPs are in or near genes that are highly expressed in the brain (HECW2, HIP1, BIN2, GRIA1), genes involved in neural development and function (KCNQ4, LMO4, GRIA1, NETO1) and autophagy (ATG4C), and genes that are associated with risk of various diseases including cancer and Alzheimers disease. In addition to considerable overlap between the traits, pathway and network analysis corroborated these findings. These findings indicate that variation in genes involved in neurological processes may be an important factor in regulating aging free of major disease and achieving longevity.