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

Objective

Age at adiposity rebound (AR) is associated with obesity and Type 2 Diabetes in adults. The aim of the present study was to investigate the role of age at AR in adult fat mass, fat distribution and pubertal timing for a Swedish cohort.

Patients and Methods

This is a retrospective cohort study. Detailed growth charts were retrieved for the men participating in the population-based GOOD (Gothenburg Osteoporosis and Obesity Determinants) study (n = 573). Body composition was analysed using dual X-ray absorptiometry and computed tomography at 18–20 years of age. Age and BMI at AR were calculated using pediatric growth charts and AR was defined as the lowest BMI between 3 and 9 years of age.

Results

Subjects were divided into early (age at AR below 5.4 years of age), middle (age at AR 5.4 to 6.8 years of age) and late (age at AR after 6.8 years of age) age at AR tertiles. Subjects in the early age at AR tertile had higher young adult BMI (+8%), whole body fat mass (+34%) and amount of subcutaneous adipose tissue (+61%) than the subjects in the middle and late tertiles (p<0.01). The early age at AR tertile had an increased risk of obesity (Odds Ratio 4.1 [95% CI 1.2–13.9]) compared with the middle and late tertiles. In addition, the early age at AR tertile had Peak Height Velocity (PHV) 7 months earlier than the late tertile.

Conclusions

Early age at AR was associated with young adult obesity as a consequence of a high amount of subcutaneous adipose tissue in men. In addition we made the novel observation that early age at AR was associated with an early puberty in men.

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.

Author Summary

Genetic investigations on bone mineral density (BMD) variation in children allow the identification of factors determining peak bone mass and their influence on developing osteoporosis later in life. We ran a genome-wide association study (GWAS) for total body BMD based on 2,660 children of different ethnic backgrounds, followed by replication in an additional 12,066 individuals comprising children, young adults, and elderly populations. Our GWAS meta-analysis identified two independent signals in the 7q31.31 locus, arising from SNPs in the vicinity of WNT16, FAM3C, and C7orf58. These variants were also associated with skull BMD, a skeletal trait with much less environmental influence for which one of the signals displayed age-specific effects. Integration of functional studies in a Wnt16 knockout mouse model and gene expression profiles in human bone tissue provided additional evidence that WNT16 and C7orf58 underlie the described associations. All together our findings demonstrate the relevance of these factors for bone biology, the attainment of peak bone mass, and their likely impact on bone fragility later in life.

Context

Mutations in the low-density lipoprotein receptor-related protein 5 (LRP5) gene cause rare syndromes characterized by altered bone mineral density (BMD). More common LRP5 variants may affect osteoporosis risk in the general population.

Objective

To generate large-scale evidence on whether 2 common variants of LRP5 (Val667Met, Ala1330Val) and 1 variant of LRP6 (Ile1062Val) are associated with BMD and fracture risk.

Design and Setting

Prospective, multicenter, collaborative study of individual-level data on 37 534 individuals from 18 participating teams in Europe and North America. Data were collected between September 2004 and January 2007; analysis of the collected data was performed between February and May 2007. Bone mineral density was assessed by dual-energy x-ray absorptiometry. Fractures were identified via questionnaire, medical records, or radiographic documentation; incident fracture data were available for some cohorts, ascertained via routine surveillance methods, including radiographic examination for vertebral fractures.

Main Outcome Measures

Bone mineral density of the lumbar spine and femoral neck; prevalence of all fractures and vertebral fractures.

Results

The Met667 allele of LRP5 was associated with reduced lumbar spine BMD (n =25 052 [number of participants with available data]; 20-mg/cm2 lower BMD per Met667 allele copy; P=3.3 × 10−8), as was the Val1330 allele (n = 24 812; 14-mg/cm2 lower BMD per Val1330 copy; P=2.6 × 10−9). Similar effects were observed for femoral neck BMD, with a decrease of 11 mg/cm2 (P =3.8 × 10−5) and 8 mg/cm2 (P=5.0×10−6) for the Met667 and Val1330 alleles, respectively (n=25 193). Findings were consistent across studies for both LRP5 alleles. Both alleles were associated with vertebral fractures (odds ratio [OR], 1.26; 95% confidence interval [CI], 1.08–1.47 for Met667 [2001 fractures among 20 488 individuals] and OR, 1.12; 95% CI, 1.01–1.24 for Val1330 [1988 fractures among 20 096 individuals]). Risk of all fractures was also increased with Met667 (OR, 1.14; 95% CI, 1.05–1.24 per allele [7876 fractures among 31 435 individuals)]) and Val1330 (OR, 1.06; 95% CI, 1.01–1.12 per allele [7802 fractures among 31 199 individuals]). Effects were similar when adjustments were made for age, weight, height, menopausal status, and use of hormone therapy. Fracture risks were partly attenuated by adjustment for BMD. Haplotype analysis indicated that Met667 and Val1330 variants both independently affected BMD. The LRP6 Ile1062Val polymorphism was not associated with any osteoporosis phenotype. All aforementioned associations except that between Val1330 and all fractures and vertebral fractures remained significant after multiple-comparison adjustments.

Conclusions

Common LRP5 variants are consistently associated with BMD and fracture risk across different white populations. The magnitude of the effect is modest. LRP5 may be the first gene to reach a genome-wide significance level (a conservative level of significance [herein, unadjusted P<10−7] that accounts for the many possible comparisons in the human genome) for a phenotype related to osteoporosis.

Cortical bone dimensions are important determinants of bone strength. Gender differences in cortical bone size caused by greater periosteal expansion in males than in females during the pubertal growth spurt are well established both in humans and in experimental animal models. However, the mechanism by which gender influences cortical bone size is still a matter of investigation. The role of androgens and estrogen in pubertal bone growth has been examined in human disorders as well as animal models, such as gonadectomized or sex steroid receptor knockout mice. Based on the findings that growth hormone (GH) and insulin-like growth factor I (IGF-I) are major regulators of postnatal skeletal growth, we and others have predicted that sex hormones interact with the GH/IGF-I axis to regulate cortical bone size. However, studies conflict as to whether estrogen and androgens impact cortical bone size through the canonical pathway, through GH without IGF-I mediation, through IGF-I without GH stimulation, or independent of GH/IGF-I. We review recent data on the impact of sex steroids and components of the GH/IGF axis on sexual dimorphism in bone size. While the GH/IGF-I axis is a major player in regulating peak bone size, the relative contribution of GH/IGF-dependent mechanisms to sex differences in cortical bone size remains to be established.

Growth hormone and insulin-like growth factor-1 (IGF-1) are essential for the achievement of normal longitudinal bone growth and bone mass. Preclinical studies using various knockout models have shown that both endocrine (mostly liver-derived) IGF-1 and bone-derived IGF-1 contribute to normal longitudinal skeletal growth and cortical bone size. Since bone size is an important determinant of bone strength, and hence fracture risk, we evaluated the predictive role of serum IGF-1 for fracture risk. The population-based Osteoporotic Fractures in Men Sweden cohort recently showed for the first time that older men with low serum IGF-1 levels have an increased fracture risk, especially for the two most important fracture types, hip and vertebral fractures. This association between serum IGF-1 and incident fracture risk is partly mediated via bone mineral density. Future studies are required to identify the mechanisms by which endocrine and local IGF-1 regulate skeletal growth and bone size. In addition, possible mediators affecting the impact of IGF-1 on fractures in men remain to be elucidated.

Testosterone concentrations in men are associated with cardiovascular morbidity, osteoporosis, and mortality and are affected by age, smoking, and obesity. Because of serum testosterone's high heritability, we performed a meta-analysis of genome-wide association data in 8,938 men from seven cohorts and followed up the genome-wide significant findings in one in silico (n = 871) and two de novo replication cohorts (n = 4,620) to identify genetic loci significantly associated with serum testosterone concentration in men. All these loci were also associated with low serum testosterone concentration defined as <300 ng/dl. Two single-nucleotide polymorphisms at the sex hormone-binding globulin (SHBG) locus (17p13-p12) were identified as independently associated with serum testosterone concentration (rs12150660, p = 1.2×10−41 and rs6258, p = 2.3×10−22). Subjects with ≥3 risk alleles of these variants had 6.5-fold higher risk of having low serum testosterone than subjects with no risk allele. The rs5934505 polymorphism near FAM9B on the X chromosome was also associated with testosterone concentrations (p = 5.6×10−16). The rs6258 polymorphism in exon 4 of SHBG affected SHBG's affinity for binding testosterone and the measured free testosterone fraction (p<0.01). Genetic variants in the SHBG locus and on the X chromosome are associated with a substantial variation in testosterone concentrations and increased risk of low testosterone. rs6258 is the first reported SHBG polymorphism, which affects testosterone binding to SHBG and the free testosterone fraction and could therefore influence the calculation of free testosterone using law-of-mass-action equation.

Author Summary

Testosterone is the most important testicular androgen in men. Low serum testosterone concentrations are associated with cardiovascular morbidity, metabolic syndrome, type 2 diabetes mellitus, atherosclerosis, osteoporosis, sarcopenia, and increased mortality risk. Thus, there is growing evidence that serum testosterone is a valuable biomarker of men's overall health status. Studies in male twins indicate that there is a strong heritability of serum testosterone. Here we perform a large-scale genome-wide association study to examine the effects of common genetic variants on serum testosterone concentrations. By examining 14,429 men, we show that genetic variants in the sex hormone-binding globulin (SHBG) locus and on the X chromosome are associated with a substantial variation in serum testosterone concentrations and increased risk of low testosterone. The reported associations may now be used in order to better understand the functional background of recently identified disease associations related to low testosterone. Importantly, we identified the first known genetic variant, which affects SHBG's affinity for binding testosterone and the free testosterone fraction and could therefore influence the calculation of free testosterone. This finding suggests that individual-based SHBG-testosterone affinity constants are required depending on the genotype of this single-nucleotide polymorphism.

Background

Transgenic mice with low levels of global insulin-like growth factor-I (IGF-I) throughout their life span, including pre- and postnatal development, have increased longevity. This study investigated whether specific deficiency of liver-derived, endocrine IGF-I is of importance for life span.

Methods and Findings

Serum IGF-I was reduced by approximately 80% in mice with adult, liver-specific IGF-I inactivation (LI-IGF-I-/- mice), and body weight decreased due to reduced body fat. The mean life span of LI-IGF-I-/- mice (n = 84) increased 10% vs. control mice (n = 137) (Cox's test, p<0.01), mainly due to increased life span (16%) of female mice [LI-IGF-I-/- mice (n = 31): 26.7±1.1 vs. control (n = 67): 23.0±0.7 months, p<0.001]. Male LI-IGF-I-/- mice showed only a tendency for increased longevity (p = 0.10). Energy expenditure, measured as oxygen consumption during and after submaximal exercise, was increased in the LI-IGF-I-/- mice. Moreover, microarray and RT-PCR analyses showed consistent regulation of three genes (heat shock protein 1A and 1B and connective tissue growth factor) in several body organs in the LI-IGF-I-/- mice.

Conclusions

Adult inactivation of liver-derived, endocrine IGF-I resulted in moderately increased mean life span. Body weight and body fat decreased in LI-IGF-I-/- mice, possibly due to increased energy expenditure during exercise. Genes earlier reported to modulate stress response and collagen aging showed consistent regulation, providing mechanisms that could underlie the increased mean life span in the LI-IGF-I-/- mice.

Background

Vitamin D is crucial for maintaining musculoskeletal health. Recently, vitamin D insufficiency has been linked to a number of extraskeletal disorders, including diabetes, cancer, and cardiovascular disease. Determinants of circulating 25-hydroxyvitamin D (25-OH D) include sun exposure and dietary intake, but its high heritability suggests that genetic determinants may also play a role.

Methods

We performed a genome-wide association study of 25-OH D among ∼30,000 individuals of European descent from 15 cohorts. Five cohorts were designated as discovery cohorts (n=16,125), five as in silico replication cohorts (n=9,366), and five as de novo replication cohorts (n=8,378). Association results were combined using z-score-weighted meta-analysis. Vitamin D insufficiency was defined as 25-OH D <75 nmol/L or <50 nmol/L.

Findings

Variants at three loci reached genome-wide significance in the discovery cohorts, and were confirmed in the replication cohorts: 4p12 (overall P=1.9 × 10-109 for rs2282679, in GC); 11q12 (P=2.1 × 10-27 for rs12785878, near DHCR7); 11p15 (P=3.3 × 10-20 for rs10741657, near CYP2R1). Variants at an additional locus (20q13, CYP24A1) were genome-wide significant in the pooled sample (P=6.0 × 10-10 for rs6013897). A genotype score was constructed using the three confirmed variants. Those in the top quartile of genotype scores had 2- to 2.5-fold elevated odds of vitamin D insufficiency (P≤1 × 10-26).

Interpretation

Variants near genes involved in cholesterol synthesis (DHCR7), hydroxylation (CYP2R1, CYP24A1), and vitamin D transport (GC) influence vitamin D status. Genetic variation at these loci identifies individuals of European descent who have substantially elevated risk of vitamin D insufficiency.

Androgens are important regulators of bone mass but the relative importance of testosterone (T) versus dihydrotestosterone (DHT) for the activation of the androgen receptor (AR) in bone is unknown. 5α-reductase is responsible for the irreversible conversion of T to the more potent AR activator DHT. There are two well established isoenzymes of 5α-reductase (type 1 and type 2), encoded by separate genes (Srd5a1 and Srd5a2). 5α-reductase type 2 is predominantly expressed in male reproductive tissues whereas 5α-reductase type 1 is highly expressed in liver and moderately expressed in several other tissues including bone. The aim of the present study was to investigate the role of 5α-reductase type 1 for bone mass using Srd5a1−/− mice. Four-month-old male Srd5a1−/− mice had reduced trabecular bone mineral density (−36%, p<0.05) and cortical bone mineral content (−15%, p<0.05) but unchanged serum androgen levels compared with wild type (WT) mice. The cortical bone dimensions were reduced in the male Srd5a1−/− mice as a result of a reduced cortical periosteal circumference compared with WT mice. T treatment increased the cortical periosteal circumference (p<0.05) in orchidectomized WT mice but not in orchidectomized Srd5a1−/− mice. Male Srd5a1−/− mice demonstrated a reduced forelimb muscle grip strength compared with WT mice (p<0.05). Female Srd5a1−/− mice had slightly increased cortical bone mass associated with elevated circulating levels of androgens. In conclusion, 5α-reductase type 1 inactivated male mice have reduced bone mass and forelimb muscle grip strength and we propose that these effects are due to lack of 5α-reductase type 1 expression in bone and muscle. In contrast, the increased cortical bone mass in female Srd5a1−/− mice, is an indirect effect mediated by elevated circulating androgen levels.

Introduction

Postmenopausal patients with rheumatoid arthritis (RA) are often treated with corticosteroids. Loss of estrogen, the inflammatory disease and exposure to corticosteroids all contribute to the development of osteoporosis. Therefore, our aim was to investigate if addition of the selective estrogen receptor modulator raloxifene, or estradiol, could prevent loss of bone mineral density in ovariectomized and dexamethasone treated mice with collagen-induced arthritis (CIA).

Methods

Female DBA/1-mice were ovariectomized or sham-operated, and CIA was induced. Treatment with dexamethasone (Dex) (125 μg/d), estradiol (E2) (1 μg/d) or raloxifene (Ral) (120 μg/day) alone, or the combination of Dex + E2 or Dex + Ral, was started after disease onset, and continued until termination of the experiments. Arthritic paws were collected for histology and one of the femoral bones was used for measurement of bone mineral density.

Results

Dex-treatment alone protected against arthritis and joint destruction, but had no effect on osteoporosis in CIA. However, additional treatment with either Ral or E2 resulted in completely preserved bone mineral density.

Conclusions

Addition of raloxifene or estradiol to dexamethasone-treatment in experimental postmenopausal polyarthritis prevents generalized bone loss.

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.

Author Summary

Dehydroepiandrosterone sulphate (DHEAS), mainly secreted by the adrenal gland, is the most abundant circulating steroid in humans. It shows a significant physiological decline after the age of 25 and diminishes about 95% by the age of 85 years, which has led to speculation that a relative DHEAS deficiency may contribute to the development of common age-related diseases or diminished longevity. Twin- and family-based studies have shown that there is a substantial genetic effect with heritability estimate of 60%, but no specific genes regulating serum DHEAS concentration have been identified to date. Here we take advantage of recent technical and methodological advances to examine the effects of common genetic variants on serum DHEAS concentrations. By examining 14,846 Caucasian individuals, we show that eight common genetic variants are associated with serum DHEAS concentrations. Genes at or near these genetic variants include BCL2L11, ARPC1A, ZKSCAN5, TRIM4, HHEX, CYP2C9, BMF, and SULT2A1. These genes have various associations with steroid hormone metabolism—co-morbidities of ageing including type 2 diabetes, lymphoma, actin filament assembly, drug and xenobiotic metabolism, and zinc finger proteins—suggesting a wider functional role for DHEAS than previously thought.

RHO family proteins are important for the function of inflammatory cells. They are modified with a 20-carbon geranylgeranyl lipid in a process catalyzed by protein geranylgeranyltransferase type I (GGTase-I). Geranylgeranylation is viewed as essential for the membrane targeting and activity of RHO proteins. Consequently, inhibiting GGTase-I to interfere with RHO protein activity has been proposed as a strategy to treat inflammatory disorders. However, here we show that mice lacking GGTase-I in macrophages develop severe joint inflammation resembling erosive rheumatoid arthritis. The disease was initiated by the GGTase-I–deficient macrophages and was transplantable and reversible in bone marrow transplantation experiments. The cells accumulated high levels of active GTP-bound RAC1, CDC42, and RHOA, and RAC1 remained associated with the plasma membrane. Moreover, GGTase-I deficiency activated p38 and NF-κB and increased the production of proinflammatory cytokines. The results challenge the view that geranylgeranylation is essential for the activity and localization of RHO family proteins and suggest that reduced geranylgeranylation in macrophages can initiate erosive arthritis.

Background

Collagen-induced arthritis (CIA) is an often-used murine model for human rheumatoid arthritis (RA). Earlier studies have shown potent anti-arthritic effects with the female sex hormone estradiol and the selective estrogen receptor modulator (SERM) raloxifene in CIA in DBA/1-mice. B10.Q-ncf1*/*mice are B10.Q mice with a mutated Ncf1 gene. In B10.Q-ncf1*/*mice, CIA develops as a chronic relapsing disease, which more accurately mimics human RA. We investigated the role of endogenous and exogenous sex steroids and raloxifene in the course of this model of chronic arthritis. We also examined whether treatment would prevent the development of inflammation-triggered generalized osteoporosis.

Methods

Female B10.Q-ncf1*/*mice were sham-operated or ovariectomized, and CIA was induced. 22 days later, when 30% of the mice had developed arthritis, treatment with raloxifene, estradiol or vehicle was started, and the clinical disease was evaluated continuously. Treatment was continued until day 56 after immunization. At termination of the experiment (day 73), bone mineral density (BMD) was analyzed, paws were collected for histological examination, and sera were analyzed for markers of cartilage turnover and pro-inflammatory cytokines.

Results

Raloxifene and estradiol treatment, as well as endogenous estrogen, decreased the frequency of arthritis, prevented joint destruction and countered generalized osteoporosis. These effects were associated with lower serum levels of the pro-inflammatory cytokine IL-6.

Conclusions

This is the first study to show that raloxifene and estradiol can ameliorate established erosive arthritis and inflammation-triggered osteoporosis in this chronic arthritis model. We propose that treatment with raloxifene could be a beneficial addition to the treatment of postmenopausal RA.

Context

Regulation of fat mass appears to be associated with immune functions. Studies of knockout mice show that endogenous interleukin (IL)-6 can suppress mature-onset obesity.

Objective

To systematically investigate associations of single nucleotide polymorphisms (SNPs) near the IL-6 (IL6) and IL-6 receptor (IL6R) genes with body fat mass, in support for our hypothesis that variants of these genes can be associated with obesity.

Design and Study Subjects

The Gothenburg Osteoporosis and Obesity Determinants (GOOD) study is a population-based cross-sectional study of 18-20 years old men (n=1 049), from the Gothenburg area (Sweden). Major findings were confirmed in two additional cohorts consisting of elderly men from the Osteoporotic Fractures in Men (MrOS) Sweden (n=2 851) and MrOS US (n=5 611) multicenter population-based studies.

Main Outcome

The genotype distributions and their association with fat mass in different compartments, measured with dual-energy X-ray absorptiometry (DXA).

Results

Out of 18 evaluated tag single nucleotide polymorphisms (SNPs) near the IL6 and IL6R genes, a recently identified SNP rs10242595 G/A [minor allele frequency (MAF) = 29%] 3′ of the IL6 gene was negatively associated with the primary outcome total body fat mass (effect size -0.11 standard deviation (SD) units/A allele, P=0.02). This negative association with fat mass was also confirmed in the combined MrOS Sweden and MrOS US cohorts (effect size -0.05 SD units/A allele; P=0.002). When all three cohorts were combined (n= 8 927, Caucasian subjects), rs10242595*A showed a negative association with total body fat mass (effect size -0.05 SD units/A allele, P<0.0002). Furthermore, the rs10242595*A was associated with low body mass index [(BMI, effect size -0.03, P<0.001)] and smaller regional fat masses. None of the other SNPs investigated in the GOOD study were reproducibly associated with body fat.

Conclusions

The IL6 gene polymorphism rs10242595*A is associated with decreased fat mass in three combined cohorts of 8 927 Caucasian men.

Previous genome-wide association (GWA) studies have identified SNPs associated with areal bone mineral density (aBMD). However, this measure is influenced by several different skeletal parameters, such as periosteal expansion, cortical bone mineral density (BMDC) cortical thickness, trabecular number, and trabecular thickness, which may be under distinct biological and genetic control. We have carried out a GWA and replication study of BMDC, as measured by peripheral quantitative computed tomography (pQCT), a more homogenous and valid measure of actual volumetric bone density. After initial GWA meta-analysis of two cohorts (ALSPAC n = 999, aged ∼15 years and GOOD n = 935, aged ∼19 years), we attempted to replicate the BMDC associations that had p<1×10−5 in an independent sample of ALSPAC children (n = 2803) and in a cohort of elderly men (MrOS Sweden, n = 1052). The rs1021188 SNP (near RANKL) was associated with BMDC in all cohorts (overall p = 2×10−14, n = 5739). Each minor allele was associated with a decrease in BMDC of ∼0.14SD. There was also evidence for an interaction between this variant and sex (p = 0.01), with a stronger effect in males than females (at age 15, males −6.77mg/cm3 per C allele, p = 2×10−6; females −2.79 mg/cm3 per C allele, p = 0.004). Furthermore, in a preliminary analysis, the rs1021188 minor C allele was associated with higher circulating levels of sRANKL (p<0.005). We show this variant to be independent from the previously aBMD associated SNP (rs9594738) and possibly from a third variant in the same RANKL region, which demonstrates important allelic heterogeneity at this locus. Associations with skeletal parameters reflecting bone dimensions were either not found or were much less pronounced. This finding implicates RANKL as a locus containing variation associated with volumetric bone density and provides further insight into the mechanism by which the RANK/RANKL/OPG pathway may be involved in skeletal development.

Author Summary

Previous studies that have identified genetic polymorphisms involved in bone density have used a technique that cannot differentiate between cortical and trabecular bone. We have carried out the first genome-wide association study using a bone scanning method that can differentiate between the constituent parts of bone. We found a genetic variant (rs1021188) near the RANKL gene that was associated with the density of cortical bone in the three cohorts that we studied (ranging in age from 15 to 78 years old). We also found that this variant may have a more prominent effect on cortical bone density in males than females. In addition, the minor C allele of rs1021188 was associated with higher circulating levels of free RANKL. Although the RANKL gene has been previously identified as being important for bone structure (albeit with a different SNP showing association), we show for the first time that this may be primarily due to its influence on the density of cortical bone, rather than the size of the bone or other bone features.

Introduction

Our objective was to determine the frequency of and factors associated with prevalent vertebral compression fractures in female systemic lupus erythematosus (SLE) patients attending rheumatologists in western Sweden.

Methods

In this cross sectional study 150 women were included. They were examined with x-ray of thoracic and lumbar spine (Th4 to L4). A reduction of at least 20% of any vertebral height, assessed by Genant's semiquantitative method, was defined as a fracture. Bone mineral density (BMD) was measured by dual energy x-ray absorptiometry (DXA).

Results

Median patient age was 47 years (20 to 82) and disease duration 11 years (1 to 41). Only 6 (4%) women had a history of clinical compressions whereas 43 (29%) had at least one radiological fracture each. The patients with at least one fracture at any site were characterized by older age (P < 0.001), being postmenopausal (P < 0.01), higher Systemic Lupus International Collaborative Clinics Damage Index (P < 0.05), lower BMD total hip and femoral neck (P < 0.05), more peripheral fractures (P < 0.01), medication with bisphosphonates (P <0.05) and calcium and vitamin D3 (P < 0.05). There were no significant differences regarding current or cumulative glucocorticosteroid dose between the groups. In logistic regression analyses high age remained as a risk factor of at least one vertebral fracture at any site whereas low BMD in total hip was associated with vertebral fracture in the lumbar spine.

Conclusions

Radiological compression fractures are common but seldom diagnosed in SLE patients. High age and low BMD in total hip, but not in spine, was associated with vertebral fractures.

OBJECTIVE

The amount of visceral adipose tissue is a risk factor for the metabolic syndrome. It is unclear how BMI changes during childhood and adolescence predict adult fat distribution. We hypothesized that there are critical periods during development for the prediction of adult subcutaneous and visceral fat mass by BMI changes during childhood and adolescence.

RESEARCH DESIGN AND METHODS

Detailed growth charts were retrieved for the men participating in the population-based Gothenburg Osteoporosis and Obesity Determinants (GOOD) Study (n = 612). Body composition was analyzed using dual-energy X-ray absorptiometry and adipose tissue areas using abdominal computed tomography at 18 to 20 years of age.

RESULTS

The main finding in the present study was that subjects with increases in BMI Z score of more than 1 SD during adolescence had, independent of prepubertal BMI, both larger subcutaneous (+138%; P < 0.001) and visceral adipose tissue areas (+91%; P < 0.001) than subjects with unchanged BMI Z-score. In contrast, subjects with increases in BMI Z score of more than 1 SD during late childhood had a larger amount of adult subcutaneous adipose tissue (+83%; P < 0.001) than subjects with unchanged BMI Z score but an unaffected amount of visceral adipose tissue. BMI changes during adolescence predict both visceral and subcutaneous adipose tissue of the abdomen, whereas BMI changes during late childhood predict only the subcutaneous adipose tissue.

CONCLUSIONS

The amount of visceral adipose tissue in young adult men was associated with BMI changes specifically during adolescence, whereas the amount of subcutaneous adipose tissue was associated with BMI changes during both late childhood and adolescence.

Introduction

Dichloroacetate (DCA) has been in clinical use for the treatment of lactacidosis and inherited mitochondrial disorders. It has potent anti-tumor effects both in vivo and in vitro, facilitating apoptosis and inhibiting proliferation. The pro-apoptotic and anti-proliferative properties of DCA prompted us to investigate the effects of this compound in arthritis.

Methods

In the present study, we used DCA to treat murine collagen type II (CII)-induced arthritis (CIA), an experimental model of rheumatoid arthritis. DBA/1 mice were treated with DCA given in drinking water.

Results

Mice treated with DCA displayed much slower onset of CIA and significantly lower severity (P < 0.0001) and much lower frequency (36% in DCA group vs. 86% in control group) of arthritis. Also, cartilage and joint destruction was significantly decreased following DCA treatment (P = 0.005). Moreover, DCA prevented arthritis-induced cortical bone mineral loss. This clinical picture was also reflected by lower levels of anti-CII antibodies in DCA-treated versus control mice, indicating that DCA affected the humoral response. In contrast, DCA had no effect on T cell- or granulocyte-mediated responses. The beneficial effect of DCA was present in female DBA/1 mice only. This was due in part to the effect of estrogen, since ovariectomized mice did not benefit from DCA treatment to the same extent as sham-operated controls (day 30, 38.7% of ovarectomized mice had arthritis vs. only 3.4% in sham-operated group).

Conclusion

Our results indicate that DCA delays the onset and alleviates the progression of CIA in an estrogen-dependent manner.

Current genome-wide association studies (GWAS) are moving towards the use of large cohorts of primary cell lines to study a disease of interest and to assign biological relevance to the genetic signals identified. Here, we use a panel of human osteoblasts (HObs) to carry out a transcriptomic survey, similar to recent studies in lymphoblastoid cell lines (LCLs). The distinct nature of HObs and LCLs is reflected by the preferential grouping of cell type–specific genes within biologically and functionally relevant pathways unique to each tissue type. We performed cis-association analysis with SNP genotypes to identify genetic variations of transcript isoforms, and our analysis indicates that differential expression of transcript isoforms in HObs is also partly controlled by cis-regulatory genetic variants. These isoforms are regulated by genetic variants in both a tissue-specific and tissue-independent fashion, and these associations have been confirmed by RT–PCR validation. Our study suggests that multiple transcript isoforms are often present in both tissues and that genetic control may affect the relative expression of one isoform to another, rather than having an all-or-none effect. Examination of the top SNPs from a GWAS of bone mineral density show overlap with probeset associations observed in this study. The top hit corresponding to the FAM118A gene was tested for association studies in two additional clinical studies, revealing a novel transcript isoform variant. Our approach to examining transcriptome variation in multiple tissue types is useful for detecting the proportion of genetic variation common to different cell types and for the identification of cell-specific isoform variants that may be functionally relevant, an important follow-up step for GWAS.

Author Summary

The transcriptome of any given cell type is a complex program of controlled gene expression underlying its biological function. An additional layer of molecular complexity involving individual genetic variation can modulate the transcriptome within the same tissue type, conferring potential phenotypic differences between individuals at the cellular level. This study highlights common and unique aspects of the transcriptome between the well-characterized lymphoblastoid cell lines from the International HapMap Project and those of a cultured primary cell type, human osteoblasts. We observe that inter-individual genetic variation can regulate transcript isoform expression in tissue-specific and tissue-independent manners, indicating that genetic differences among individuals can alter the transcriptome in one or more tissues, ultimately leading to altered biological functions within the lymphoblasts and/or osteoblasts. Pursuant to this, genome wide association studies on bone mineral density (BMD) have identified a number of significant loci and polymorphisms highly linked to the BMD quantitative phenotype. A small proportion of these polymorphisms overlap with our highly significant SNPs regulating the osteoblast transcriptome, revealing a potential molecular basis for this phenotype at the transcriptional level. This study highlights the importance of examining the differing transcriptomes and cis-regulatory mechanisms governing the biological and functional roles of varied tissue types.

Northwestern Europeans are among the tallest of human populations. The increase in body height in these people appears to have reached a plateau, suggesting the ubiquitous presence of an optimal environment in which genetic factors may have exerted a particularly strong influence on human growth. Therefore, we performed a genome-wide association study (GWAS) of body height using 2.2 million markers in 10 074 individuals from three Dutch and one German population-based cohorts. Upon genotyping, the 12 most significantly height-associated single nucleotide polymorphisms (SNPs) from this GWAS in 6912 additional individuals of Dutch and Swedish origin, a genetic variant (rs6717918) on chromosome 2q37.1 was found to be associated with height at a genome-wide significance level (Pcombined = 3.4 × 10−9). Notably, a second SNP (rs6718438) located ∼450 bp away and in strong LD (r2 = 0.77) with rs6717918 was previously found to be suggestive of a height association in 29 820 individuals of mainly northwestern European ancestry, and the over-expression of a nearby natriuretic peptide precursor type C (NPPC) gene, has been associated with overgrowth and skeletal anomalies. We also found a SNP (rs10472828) located on 5p14 near the natriuretic peptide receptor 3 (NPR3) gene, encoding a receptor of the NPPC ligand, to be associated with body height (Pcombined = 2.1 × 10−7). Taken together, these results suggest that variation in the C-type natriuretic peptide signaling pathway, involving the NPPC and NPR3 genes, plays an important role in determining human body height.

Background

Hormone-sensitive lipase (HSL) is expressed predominantly in adipose tissue, where it plays an important role in catecholamine-stimulated hydrolysis of stored tri- and diglycerides, thus mobilizing fatty acids. HSL exhibits broad substrate specificity and besides acylglycerides it hydrolyzes cholesteryl esters, retinyl esters and lipoidal esters. Despite its role in fatty acid mobilization, HSL null mice have been shown to be resistant to diet-induced obesity.

Methodology/Principal Findings

Following a high-fat diet (HFD) regimen, energy expenditure, measured using indirect calorimetry, was increased in HSL null mice. White adipose tissue of HSL null mice was characterized by reduced mass and reduced protein expression of PPARγ, a key transcription factor in adipogenesis, and stearoyl-CoA desaturase 1, the expression of which is known to be positively correlated to the differentiation state of the adipocyte. The protein expression of uncoupling protein-1 (UCP-1), the highly specific marker of brown adipocytes, was increased 7-fold in white adipose tissue of HSL null mice compared to wildtype littermates. Transmission electron microscopy revealed an increase in the size of mitochondria of white adipocytes of HSL null mice. The mRNA expression of pRb and RIP140 was decreased in isolated white adipocytes, while the expression of UCP-1 and CPT1 was increased in HSL null mice compared to wildtype littermates. Basal oxygen consumption was increased almost 3-fold in white adipose tissue of HSL null mice and was accompanied by increased uncoupling activity.

Conclusions

These data suggest that HSL is involved in the determination of white versus brown adipocytes during adipocyte differentiation The exact mechanism(s) underlying this novel role of HSL remains to be elucidated, but it seems clear that HSL is required to sustain normal expression levels of pRb and RIP140, which both promote differentiation into the white, rather than the brown, adipocyte lineage.

Background

Estrogens and androgens have extensive effects on the immune system, for example they suppress both T and B lymphopoiesis in thymus and bone marrow. Submandibular glands are sexually dimorphic in rodents, resulting in larger granular convoluted tubules in males compared to females. The aim of the present experiments was to investigate the estrogenic and androgenic effects of 4-estren-3α,17β-diol (estren) on thymus, bone marrow and submandibular glands, and compare the effects to those of 17β-estradiol (E2) and 5α-dihydrotestosterone (DHT), respectively. Estrogen receptors (ERs) were blocked by treatment of mice with the ER-antagonist ICI 182,780; also, knock-out mice lacking one or both ERs were used.

Results

As expected, the presence of functional ERs was mandatory for all the effects of E2. Similar to DHT-treatment, estren-treatment resulted in decreased thymus weight, as well as decreased frequency of bone marrow B cells. Treatment with estren or DHT also resulted in a shift in submandibular glands towards an androgen phenotype. All the effects of estren and DHT were independent of ERs.

Conclusion

Our study is the first to show that estren has similar effects as the androgen DHT on lymphopoiesis in thymus and bone marrow, and on submandibular glands, and that these effects are independent of estrogen receptors. This supports the hypothesis of estren being able to signal through the androgen receptor.