Osteoporotic fracture is among the most common and costly of diseases. While reasonably heritable, its genetic determinants have remained elusive. Forearm fractures are the most common clinically recognized osteoporotic fractures with a relatively high heritability. To establish an atlas of the genetic determinants of forearm fractures, we performed genome-wide association analyses including 100,026 forearm fracture cases. We identified 43 loci, including 26 new fracture loci. Although most fracture loci associated with bone mineral density, we also identified loci that primarily regulate bone quality parameters. Functional studies of one such locus, at TAC4, revealed that Tac4-/- mice have reduced mechanical bone strength. The strongest forearm fracture signal, at WNT16, displayed remarkable bone-site-specificity with no association with hip fractures. Tall stature and low body mass index were identified as new causal risk factors for fractures. The insights from this atlas may improve fracture prediction and enable therapeutic development to prevent fractures.

Appendicular lean mass (ALM) associates with mobility and bone mineral density (BMD). While associations between gut microbiota composition and ALM have been reported, previous studies rely on relatively small sample sizes. Here, we determine the associations between prevalent gut microbes and ALM in large discovery and replication cohorts with information on relevant confounders within the population-based Norwegian HUNT cohort (n = 5196, including women and men). We show that the presence of three bacterial species - Coprococcus comes, Dorea longicatena, and Eubacterium ventriosum - are reproducibly associated with higher ALM. When combined into an anabolic species count, participants with all three anabolic species have 0.80 kg higher ALM than those without any. In an exploratory analysis, the anabolic species count is positively associated with femoral neck and total hip BMD. We conclude that the anabolic species count may be used as a marker of ALM and BMD. The therapeutic potential of these anabolic species to prevent sarcopenia and osteoporosis needs to be determined.

Recent preclinical studies reported that the BCL-2 inhibitor venetoclax can impair bone growth. A strategy to prevent such a side effect of this promising anticancer drug is highly desired. Earlier in vitro and in vivo studies suggested that the mitochondrial peptide humanin has the potential to prevent drug-induced growth impairment.

Excess vitamin A has been associated with decreased cortical bone thickness and increased fracture risk. While most studies in rodents have employed high dosages of vitamin A for short periods of time, we investigated the bone phenotype in mice after longer exposure to more clinically relevant doses. For 1, 4 and 10 weeks, mice were fed a control diet (4.5 µg retinyl acetate/g chow), a diet modeled from the human upper tolerable limit (UTL; 20 µg retinyl acetate/g chow) and a diet three times UTL (supplemented; 60 µg retinyl acetate/g chow). Time-dependent decreases in periosteal circumference and bone mineral content were noted with the supplemented dose. These reductions in cortical bone resulted in a significant time-dependent decrease of predicted strength and a non-significant trend toward reduced bone strength as analyzed by three-point bending. Trabecular bone in tibiae and vertebrae remained unaffected when vitamin A was increased in the diet. Dynamic histomorphometry demonstrated that bone formation was substantially decreased after 1 week of treatment at the periosteal site with the supplemental dose. Increasing amount of vitamin A decreased endocortical circumference, resulting in decreased marrow area, a response associated with enhanced endocortical bone formation. In the presence of bisphosphonate, vitamin A had no effect on cortical bone, suggesting that osteoclasts are important, even if effects on bone resorption were not detected by osteoclast counting, genes in cortical bone or analysis of serum TRAP5b and CTX. In conclusion, our results indicate that even clinically relevant doses of vitamin A have a negative impact on the amount of cortical bone.

We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10-9 , adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2  > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

Osteoporosis is a common aging-related disease diagnosed primarily using bone mineral density (BMD). We assessed genetic determinants of BMD as estimated by heel quantitative ultrasound in 426,824 individuals, identifying 518 genome-wide significant loci (301 novel), explaining 20% of its variance. We identified 13 bone fracture loci, all associated with estimated BMD (eBMD), in ~1.2 million individuals. We then identified target genes enriched for genes known to influence bone density and strength (maximum odds ratio (OR) = 58, P = 1 × 10-75) from cell-specific features, including chromatin conformation and accessible chromatin sites. We next performed rapid-throughput skeletal phenotyping of 126 knockout mice with disruptions in predicted target genes and found an increased abnormal skeletal phenotype frequency compared to 526 unselected lines (P < 0.0001). In-depth analysis of one gene, DAAM2, showed a disproportionate decrease in bone strength relative to mineralization. This genetic atlas provides evidence linking associated SNPs to causal genes, offers new insight into osteoporosis pathophysiology, and highlights opportunities for drug development.

Androgens have important cardiometabolic actions in males, but their metabolic role in females is unclear. To determine the physiologic androgen receptor (AR)-dependent actions of androgens on atherogenesis in female mice, we generated female AR-knockout (ARKO) mice on an atherosclerosis-prone apolipoprotein E (apoE)-deficient background. After 8 weeks on a high-fat diet, but not on a normal chow diet, atherosclerosis in aorta was increased in ARKO females (+59% vs. control apoE-deficient mice with intact AR gene). They also displayed increased body weight (+18%), body fat percentage (+62%), and hepatic triglyceride levels, reduced insulin sensitivity, and a marked atherogenic dyslipidemia (serum cholesterol, +52%). Differences in atherosclerosis, body weight, and lipid levels between ARKO and control mice were abolished in mice that were ovariectomized before puberty, consistent with a protective action of ovarian androgens mediated via the AR. Furthermore, the AR agonist dihydrotestosterone reduced atherosclerosis (-41%; thoracic aorta), subcutaneous fat mass (-44%), and cholesterol levels (-35%) in ovariectomized mice, reduced hepatocyte lipid accumulation in hepatoma cells in vitro, and regulated mRNA expression of hepatic genes pivotal for lipid homeostasis. In conclusion, we demonstrate that the AR protects against diet-induced atherosclerosis in female mice and propose that this is mediated by modulation of body composition and lipid metabolism.

In a previous population-based study on 3369 European men with self-reported prostate cancer (PCa), it was shown that androgen receptor (AR) haplotype designated H2 was associated with high levels of serum PSA (prostate-specific antigen) concentration, and, at the same time, with low risk for PCa. The aim of this study was to replicate this finding in other cohorts, with registry-based cancer diagnosis.

WISP2 is a novel adipokine, most highly expressed in the adipose tissue and primarily in undifferentiated mesenchymal cells. As a secreted protein, it is an autocrine/paracrine activator of canonical WNT signaling and, as an intracellular protein, it helps to maintain precursor cells undifferentiated. To examine effects of increased WISP2 in vivo, we generated an aP2-WISP2 transgenic (Tg) mouse. These mice had increased serum levels of WISP2, increased lean body mass and whole body energy expenditure, hyperplastic brown/white adipose tissues and larger hyperplastic hearts. Obese Tg mice remained insulin sensitive, had increased glucose uptake by adipose cells and skeletal muscle in vivo and ex vivo, increased GLUT4, increased ChREBP and markers of adipose tissue lipogenesis. Serum levels of the novel fatty acid esters of hydroxy fatty acids (FAHFAs) were increased and transplantation of Tg adipose tissue improved glucose tolerance in recipient mice supporting a role of secreted FAHFAs. The growth-promoting effect of WISP2 was shown by increased BrdU incorporation in vivo and Tg serum increased mesenchymal precursor cell proliferation in vitro. In contrast to conventional canonical WNT ligands, WISP2 expression was inhibited by BMP4 thereby allowing normal induction of adipogenesis. WISP2 is a novel secreted regulator of mesenchymal tissue cellularity.

Leptin is an adipocyte-secreted hormone, the circulating levels of which correlate closely with overall adiposity. Although rare mutations in the leptin (LEP) gene are well known to cause leptin deficiency and severe obesity, no common loci regulating circulating leptin levels have been uncovered. Therefore, we performed a genome-wide association study (GWAS) of circulating leptin levels from 32,161 individuals and followed up loci reaching P<10(-6) in 19,979 additional individuals. We identify five loci robustly associated (P<5 × 10(-8)) with leptin levels in/near LEP, SLC32A1, GCKR, CCNL1 and FTO. Although the association of the FTO obesity locus with leptin levels is abolished by adjustment for BMI, associations of the four other loci are independent of adiposity. The GCKR locus was found associated with multiple metabolic traits in previous GWAS and the CCNL1 locus with birth weight. Knockdown experiments in mouse adipose tissue explants show convincing evidence for adipogenin, a regulator of adipocyte differentiation, as the novel causal gene in the SLC32A1 locus influencing leptin levels. Our findings provide novel insights into the regulation of leptin production by adipose tissue and open new avenues for examining the influence of variation in leptin levels on adiposity and metabolic health.

Lipoprotein synthesis is controlled by estrogens, but the exact mechanisms underpinning this regulation and the role of the hepatic estrogen receptor α (ERα) in cholesterol physiology are unclear. Utilizing a mouse model involving selective ablation of ERα in the liver, we demonstrate that hepatic ERα couples lipid metabolism to the reproductive cycle. We show that this receptor regulates the synthesis of cholesterol transport proteins, enzymes for lipoprotein remodeling, and receptors for cholesterol uptake. Additionally, ERα is indispensable during proestrus for the generation of high-density lipoproteins efficient in eliciting cholesterol efflux from macrophages. We propose that a specific interaction with liver X receptor α (LXRα) mediates the broad effects of ERα on the hepatic lipid metabolism.

Cancer of unknown primary site (CUP) is a fatal cancer diagnosed through metastases at various organs. Little is known about germline genetics of CUP which appears worth of a search in view of reported familial associations in CUP. In the present study, samples from CUP patients were identified from 2 Swedish biobanks and a German clinical trial, totaling 578 CUP patients and 7628 regionally matched controls. Diagnostic data specified the organ where metastases were diagnosed. We carried out a genome-wide association study on CUP cases and controls. In the whole sample set, 6 loci reached an allelic p-value in the range of 10-7 and were supported by data from the three centers. Three associations were located next to non-coding RNA genes. rs2660852 flanked 5'UTR of LTA4H (leukotriene A4 hydrolase), rs477145 was intronic to TIAM1 (T-cell lymphoma invasion and metastases) and rs2835931 was intronic to KCNJ6 (potassium channel, inwardly rectifying subfamily J, member 6). In analysis of subgroups of CUP patients (smokers, non-smokers and CUP with liver metastases) genome-wide significant associations were noted. For patients with liver metastases associations on chromosome 6 and 11, the latter including a cluster of genes DHCR7 and NADSYN1, encoding key enzymes in cholesterol and NAD synthesis, and KRTAP5-7, encoding a keratin associated protein. This first GWAS on CUP provide preliminary evidence that germline genes relating to inflammation (LTA4H), metastatic promotion (TIAM1) in association with lipid metabolic disturbance (chromosome 11 cluster) may contribute to the risk of CUP.

Generalized osteoporosis in postmenopausal rheumatoid arthritis (RA) is caused both by estrogen deficiency and by the inflammatory disease. The relative importance of each of these factors is unknown. The aim of this study was to establish a murine model of osteoporosis in postmenopausal RA, and to evaluate the relative importance and mechanisms of menopause and arthritis-related osteoporosis. To mimic postmenopausal RA, DBA/1 mice were ovariectomized, followed by the induction of type II collagen-induced arthritis. After the mice had been killed, paws were collected for histology, one femur for bone mineral density (BMD) and sera for analyses of markers of bone resorption (RatLaps; type I collagen cross-links, bone formation (osteocalcin) and cartilage destruction (cartilage oligomeric matrix protein), and for the evaluation of antigen-specific and innate immune responsiveness. Ovariectomized mice displayed more severe arthritis than sham-operated controls. At termination of the experiment, arthritic control mice and non-arthritic ovariectomized mice displayed trabecular bone losses of 26% and 22%, respectively. Ovariectomized mice with arthritis had as much as 58% decrease in trabecular BMD. Interestingly, cortical BMD was decreased by arthritis but was not affected by hormonal status. In addition, markers of bone resorption and cartilage destruction were increased in arthritic mice, whereas markers of bone formation were increased in ovariectomized mice. This study demonstrates that the loss of endogenous estrogen and inflammation contribute additively and equally to osteoporosis in experimental postmenopausal polyarthritis. Markers of bone remodeling and bone marrow lymphocyte phenotypes indicate different mechanisms for the development of osteoporosis caused by ovariectomy and arthritis in this model.

The gut microbiota modulates host metabolism and development of immune status. Here we show that the gut microbiota is also a major regulator of bone mass in mice. Germ-free (GF) mice exhibit increased bone mass associated with reduced number of osteoclasts per bone surface compared with conventionally raised (CONV-R) mice. Colonization of GF mice with a normal gut microbiota normalizes bone mass. Furthermore, GF mice have decreased frequency of CD4(+) T cells and CD11b(+) /GR 1 osteoclast precursor cells in bone marrow, which could be normalized by colonization. GF mice exhibited reduced expression of inflammatory cytokines in bone and bone marrow compared with CONV-R mice. In summary, the gut microbiota regulates bone mass in mice, and we provide evidence for a mechanism involving altered immune status in bone and thereby affected osteoclast-mediated bone resorption. Further studies are required to evaluate the gut microbiota as a novel therapeutic target for osteoporosis.

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.

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⁻¹⁴; LOC285735, rs271170, p = 2.7×10⁻¹²; OPG, rs7839059, p = 1.2×10⁻¹⁰; and ESR1/C6orf97, rs6909279, p = 1.1×10⁻⁹). 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⁻⁹). 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.

Approaches exploiting trait distribution extremes may be used to identify loci associated with common traits, but it is unknown whether these loci are generalizable to the broader population. In a genome-wide search for loci associated with the upper versus the lower 5th percentiles of body mass index, height and waist-to-hip ratio, as well as clinical classes of obesity, including up to 263,407 individuals of European ancestry, we identified 4 new loci (IGFBP4, H6PD, RSRC1 and PPP2R2A) influencing height detected in the distribution tails and 7 new loci (HNF4G, RPTOR, GNAT2, MRPS33P4, ADCY9, HS6ST3 and ZZZ3) for clinical classes of obesity. Further, we find a large overlap in genetic structure and the distribution of variants between traits based on extremes and the general population and little etiological heterogeneity between obesity subgroups.

We hypothesized that bone resorption acts to increase bone strength through stimulation of periosteal expansion. Hence, we examined whether bone resorption, as reflected by serum β-C-telopeptides of type I collagen (CTX), is positively associated with periosteal circumference (PC), in contrast to inverse associations with parameters related to bone remodeling such as cortical bone mineral density (BMDC ). CTX and mid-tibial peripheral quantitative computed tomography (pQCT) scans were available in 1130 adolescents (mean age 15.5 years) from the Avon Longitudinal Study of Parents and Children (ALSPAC). Analyses were adjusted for age, gender, time of sampling, tanner stage, lean mass, fat mass, and height. CTX was positively related to PC (β=0.19 [0.13, 0.24]) (coefficient=SD change per SD increase in CTX, 95% confidence interval)] but inversely associated with BMDC (β=-0.46 [-0.52,-0.40]) and cortical thickness [β=-0.11 (-0.18, -0.03)]. CTX was positively related to bone strength as reflected by the strength-strain index (SSI) (β=0.09 [0.03, 0.14]). To examine the causal nature of this relationship, we then analyzed whether single-nucleotide polymorphisms (SNPs) within key osteoclast regulatory genes, known to reduce areal/cortical BMD, conversely increase PC. Fifteen such genetic variants within or proximal to genes encoding receptor activator of NF-κB (RANK), RANK ligand (RANKL), and osteoprotegerin (OPG) were identified by literature search. Six of the 15 alleles that were inversely related to BMD were positively related to CTX (p<0.05 cut-off) (n=2379). Subsequently, we performed a meta-analysis of associations between these SNPs and PC in ALSPAC (n=3382), Gothenburg Osteoporosis and Obesity Determinants (GOOD) (n=938), and the Young Finns Study (YFS) (n=1558). Five of the 15 alleles that were inversely related to BMD were positively related to PC (p<0.05 cut-off). We conclude that despite having lower BMD, individuals with a genetic predisposition to higher bone resorption have greater bone size, suggesting that higher bone resorption is permissive for greater periosteal expansion.

Boys present with higher proportions of immature/naïve CD5+ B cells than girls up to 3 years of age. Boys also have higher fractions of regulatory T cells (Tregs) in early infancy, but the mechanisms for these sex-related differences are unknown. In the prospective FARMFLORA follow-up study of 23 boys and 25 girls, we investigated if these immunological differences remained at 8 years of age. We also examined if testosterone or dihydrotestosterone (DHT) levels at birth and at 8 years of age were associated with immune maturation. Immunological variables and androgen levels were examined and measured in blood samples obtained at birth, 3-5 days and at 8 years of age. Boys had higher proportions of CD5+ and immature/transitional CD24hiCD38hi B cells, whereas girls had higher fractions of B cells with a memory phenotype at 8 years of age. School-aged boys also presented with higher frequencies of Tregs, and a greater capacity to produce T-cell-associated cytokines. Among boys, higher cord blood DHT levels were associated with higher proportions of CD5+ B cells in early infancy and at 8 years of life. These results suggest that DHT actions in utero might be involved in the mechanism for delayed peripheral B-cell maturation in boys.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.