Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment

Genetic Variation in RIN3 in the Belgian Population Supports Its Involvement in the Pathogenesis of Paget's Disease of Bone and Modifies the Age of Onset

Paget's disease of bone (PDB) is a common, late-onset bone disorder characterized by focal increase of bone turnover. Mutations in the SQSTM1 gene are found in up to 40% of patients and recent GWAS have led to novel associations with several loci. RIN3, the candidate gene located at the associated 14q32 locus, has recently been studied in a British cohort to elucidate its contribution to the pathogenesis. In this study, we performed a genetic screening of RIN3 in an unrelated cohort to validate these findings and to further explore genetic variation in this gene in the context of PDB. In our screening, we examined the 5' untranslated region (UTR), the exonic regions and the intron-exon boundaries of the gene in a control cohort and a patient cohort. Our findings show clustering of variation similar to the British cohort and support a protective role for common genetic variation (rs117068593, p.R279C) in the proline-rich region and a functionally relevant role for rare genetic variation in the domains that mediate binding and activation of its interaction partner, Rab5. Additive regression models, fitted for the common variants, validated the association of the rs117068593 variant with the disease (OR^+/+ 0.315; OR^+/- 0.562). In addition, our analyses revealed a potentially modifying effect of this variant on the age of onset of the disease. In conclusion, our findings support the involvement of genetic variation in RIN3 in PDB and suggest a role for RIN3 as a potential modifier of the age of onset of the disease.

High-density genotyping reveals signatures of selection related to acclimation and economically important traits in 15 local sheep breeds from Russia

Background

Domestication and centuries of selective breeding have changed genomes of sheep breeds to respond to environmental challenges and human needs. The genomes of local breeds, therefore, are valuable sources of genomic variants to be used to understand mechanisms of response to adaptation and artificial selection. As a step toward this we performed a high-density genotyping and comprehensive scans for signatures of selection in the genomes from 15 local sheep breeds reared across Russia.

Results

Results demonstrated that the genomes of Russian sheep breeds contain multiple regions under putative selection. More than 50% of these regions matched with intervals identified in previous scans for selective sweeps in sheep genomes. These regions contain well-known candidate genes related to morphology, adaptation, and domestication (e.g., KITLG, KIT, MITF, and MC1R), wool quality and quantity (e.g., DSG@, DSC@, and KRT@), growth and feed intake (e.g., HOXA@, HOXC@, LCORL, NCAPG, LAP3, and CCSER1), reproduction (e.g., CMTM6, HTRA1, GNAQ, UBQLN1, and IFT88), and milk-related traits (e.g., ABCG2, SPP1, ACSS1, and ACSS2). In addition, multiple genes that are putatively related to environmental adaptations were top-ranked in selected intervals (e.g., EGFR, HSPH1, NMUR1, EDNRB, PRL, TSHR, and ADAMTS5). Moreover, we observed that multiple key genes involved in human hereditary sensory and autonomic neuropathies, and genetic disorders accompanied with an inability to feel pain and environmental temperatures, were top-ranked in multiple or individual sheep breeds from Russia pointing to a possible mechanism of adaptation to harsh climatic conditions.

Conclusions

Our work represents the first comprehensive scan for signatures of selection in genomes of local sheep breeds from the Russian Federation of both European and Asian origins. We confirmed that the genomes of Russian sheep contain previously identified signatures of selection, demonstrating the robustness of our integrative approach. Multiple novel signatures of selection were found near genes which could be related to adaptation to the harsh environments of Russia. Our study forms a basis for future work on using Russian sheep genomes to spot specific genetic variants or haplotypes to be used in efforts on developing next-generation highly productive breeds, better suited to diverse Eurasian environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5537-0) contains supplementary material, which is available to authorized users.

Genome-scale Capture C promoter interactions implicate effector genes at GWAS loci for bone mineral density

Osteoporosis is a devastating disease with an essential genetic component. GWAS have discovered genetic signals robustly associated with bone mineral density (BMD), but not the precise localization of effector genes. Here, we carry out physical and direct variant to gene mapping in human mesenchymal progenitor cell-derived osteoblasts employing a massively parallel, high resolution Capture C based method in order to simultaneously characterize the genome-wide interactions of all human promoters. By intersecting our Capture C and ATAC-seq data, we observe consistent contacts between candidate causal variants and putative target gene promoters in open chromatin for ~ 17% of the 273 BMD loci investigated. Knockdown of two novel implicated genes, ING3 at ‘CPED1-WNT16’ and EPDR1 at ‘STARD3NL’, inhibits osteoblastogenesis, while promoting adipogenesis. This approach therefore aids target discovery in osteoporosis, here on the example of two relevant genes involved in the fate determination of mesenchymal progenitors, and can be applied to other common genetic diseases.

GWAS have identified numerous genetic loci for bone mineral density (BMD) and fracture risk. Here, the authors map these variants to putative target genes using ATAC-seq and Capture C of human osteoblasts and confirm ING3 and EPDR1 as BMD genes in in vitro osteoblast differentiation experiments.

LncRNA ZBTB40-IT1 modulated by osteoporosis GWAS risk SNPs suppresses osteogenesis

Previous genome-wide linkage and association studies have identified an osteoporosis-associated locus at 1p36 that harbors SNPs rs34920465 and rs6426749. The 1p36 locus also comprises the WNT4 gene with known role in bone metabolism and functionally unknown ZBTB40/lncRNA ZBTB40-IT1 genes. How these might interact to contribute to osteoporosis susceptibility is not known. In this study, we show that lncRNA ZBTB40-IT1 is able to suppress osteogenesis and promote osteoclastogenesis by regulating the expression of WNT4, RUNX2, OSX, ALP, COL1A1, OPG and RANKL in U-2OS and hFOB1.19 cell lines, whereas ZBTB40 plays an opposite role in bone metabolism. Treatment with parathyroid hormone significantly downregulates the expression of ZBTB40-IT1 in U-2OS cell lines. ZBTB40 can suppress ZBTB40-IT1 expression but has no response to parathyroid hormone treatment. Dual-luciferase assay and biotin pull-down assay demonstrate that osteoporosis GWAS lead SNPs rs34920465-G and rs6426749-C alleles can respectively bind transcription factors JUN::FOS and CREB1, and upregulate ZBTB40 and ZBTB40-IT1 expression. Our study discovers the critical role of ZBTB40 and lncRNA ZBTB40-IT1 in bone metabolism, and provides a mechanistic basis for osteoporosis GWAS lead SNPs rs34920465 and rs6426749.

Recent Advances in the Genetics of Fractures in Osteoporosis

Genetic susceptibility, together with old age, female sex, and low bone mineral density (BMD) are amongst the strongest determinants of fracture risk. Tmost recent large-scale genome-wide association study (GWAS) meta-analysis has yielded fifteen loci. This review focuses on the advances in the research of genetic determinants of fracture risk. We first discuss the genetic architecture of fracture risk, touching upon different methods and overall findings. We then discuss in a second paragraph the most recent advances in the field and focus on the genetics of fracture risk and also of other endophenotypes closely related to fracture risk such as bone mineral density (BMD). Application of state-of-the-art methodology such as Mendelian randzation in fracture GWAS are reviewed. The final part of this review touches upon potential future directions in genetic research of osteoporotic fractures.

Glycosyltransferase genes that cause monogenic congenital disorders of glycosylation are distinct from glycosyltransferase genes associated with complex diseases

Glycosylation of proteins, lipids and proteoglycans in human cells involves at least 167 identified glycosyltransferases (GTfs), and these orchestrate the biosynthesis of diverse types of glycoconjugates and glycan structures. Mutations in this part of the genome-the GTf-genome-cause more than 58 rare, monogenic congenital disorders of glycosylation (CDGs). They are also statistically associated with a large number of complex phenotypes, diseases or predispositions to complex diseases based on Genome-Wide Association Studies (GWAS). CDGs are extremely rare and often with severe medical consequences. In contrast, GWAS are likely to identify more common genetic variations and generally involve less severe and distinct traits. We recently confirmed that structural defects in GTf genes are extremely rare, which seemed at odds with the large number of GWAS pointing to GTf-genes. To resolve this issue, we surveyed the GTf-genome for reported CDGs and GWAS candidates; we found little overlap between the two groups of genes. Moreover, GTf-genes implicated by CDG or GWAS appear to constitute different classes with respect to their: (i) predicted roles in glycosylation pathways; (ii) potential for partial redundancy by closely homologous genes; and (iii) transcriptional regulation as evaluated by RNAseq data. Our analysis suggest that more complex traits are caused by dysregulation rather than structural deficiency of GTfs, which suggests that some glycosylation reactions may be predicted to be under tight regulation for fine-tuning of important biological functions.

Characterization of expression and alternative splicing of the gene cadherin-like and PC esterase domain containing 1 (Cped1)

Cadherin-like and PC-esterase domain containing 1 (CPED1) is an uncharacterized gene with no known function. Human genome wide association studies (GWAS) for bone mineral density (BMD) have repeatedly identified a significant locus on Chromosome 7 that contains the gene CPED1, but it remains unclear if this gene could be causative. While an open reading frame for this gene has been predicted, there has been no systematic exploration of expression or alternate splicing for CPED1 in humans or mice.Using mouse models, we demonstrate that Cped1 is alternately spliced whereby transcripts are generated with exon 3 or exons 16 and 17 removed. In calvarial-derived pre-osteoblasts, Cped1 utilizes the predicted promoter upstream of exon 1 as well as alternate promoters upstream of exon 3 and exon 12.Lastly, we have determined that some transcripts terminate at the end of exon 10 and therefore do not contain the cadherin like and the PC esterase domains.Together, these data suggest that multiple protein products may be produced by this gene, with some products either lacking or containing both the predicted functional domains. Our data provide a framework upon which future functional studies will be built to understand the role of this gene in bone biology.

Genome-wide survey of parent-of-origin effects on DNA methylation identifies candidate imprinted loci in humans

Genomic imprinting is an epigenetic mechanism leading to parent-of-origin silencing of alleles. So far, the precise number of imprinted regions in humans is uncertain. In this study, we leveraged genome-wide DNA methylation in whole blood measured longitudinally at three time points (birth, childhood and adolescence) and genome-wide association studies (GWAS) data in 740 mother-child duos from the Avon Longitudinal Study of parents and children to identify candidate imprinted loci. We reasoned that cis-meQTLs at genomic regions that were imprinted would show strong evidence of parent-of-origin associations with DNA methylation, enabling the detection of imprinted regions. Using this approach, we identified genome-wide significant cis-meQTLs that exhibited parent-of-origin effects (POEs) at 82 loci, 34 novel and 48 regions previously implicated in imprinting (3.7-10<P < 10-300). Using an independent dataset from the Brisbane Systems Genetic Study, we replicated 76 out of the 82 identified loci. POEs were remarkably consistent across time points and were so strong at some loci that methylation levels enabled good discrimination of parental transmissions at these and surrounding genomic regions. The implication is that parental allelic transmissions could be modelled at many imprinted (and linked) loci in GWAS of unrelated individuals given a combination of genetic and methylation data. Novel regions showing parent of origin effects on methylation will require replication using a different technology and further functional experiments to confirm that such effects arise through a genomic imprinting mechanism.

Transcriptional profiling of murine osteoblast differentiation based on RNA-seq expression analyses

Osteoblastic differentiation is a multistep process characterized by osteogenic induction of mesenchymal stem cells, which then differentiate into proliferative pre-osteoblasts that produce copious amounts of extracellular matrix, followed by stiffening of the extracellular matrix, and matrix mineralization by hydroxylapatite deposition. Although these processes have been well characterized biologically, a detailed transcriptional analysis of murine primary calvaria osteoblast differentiation based on RNA sequencing (RNA-seq) analyses has not previously been reported. Here, we used RNA-seq to obtain expression values of 29,148 genes at four time points as murine primary calvaria osteoblasts differentiate in vitro until onset of mineralization was clearly detectable by microscopic inspection. Expression of marker genes confirmed osteogenic differentiation. We explored differential expression of 1386 protein-coding genes using unsupervised clustering and GO analyses. 100 differentially expressed lncRNAs were investigated by co-expression with protein-coding genes that are localized within the same topologically associated domain. Additionally, we monitored expression of 237 genes that are silent or active at distinct time points and compared differential exon usage. Our data represent an in-depth profiling of murine primary calvaria osteoblast differentiation by RNA-seq and contribute to our understanding of genetic regulation of this key process in osteoblast biology.

Glucocorticoids suppress Wnt16 expression in osteoblasts in vitro and in vivo

Glucocorticoid-induced osteoporosis is a frequent complication of systemic glucocorticoid (GC) therapy and mainly characterized by suppressed osteoblast activity. Wnt16 derived from osteogenic cells is a key determinant of bone mass. Here, we assessed whether GC suppress bone formation via inhibiting Wnt16 expression. GC treatment with dexamethasone (DEX) decreased Wnt16 mRNA levels in murine bone marrow stromal cells (mBMSCs) time- and dose-dependently. Similarly, Wnt16 expression was also suppressed after DEX treatment in calvarial organ cultures. Consistently, mice receiving GC-containing slow-release prednisolone pellets showed lower skeletal Wnt16 mRNA levels and bone mineral density than placebo-treated mice. The suppression of Wnt16 by GCs was GC-receptor-dependent as co-treatment of mBMSCs with DEX and the GR antagonist RU-486 abrogated the GC-mediated suppression of Wnt16. Likewise, DEX failed to suppress Wnt16 expression in GR knockout-mBMSCs. In addition, Wnt16 mRNA levels were unaltered in bone tissue of GC-treated GR dimerization-defective GR^dim mice, suggesting that GCs suppress Wnt16 via direct DNA-binding mechanisms. Consistently, DEX treatment reduced Wnt16 promoter activity in MC3T3-E1 cells. Finally, recombinant Wnt16 restored DEX-induced suppression of bone formation in mouse calvaria. Thus, this study identifies Wnt16 as a novel target of GC action in GC-induced suppression of bone formation.

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (β = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (β = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.

An Osteoporosis Risk SNP at 1p36.12 Acts as an Allele-Specific Enhancer to Modulate LINC00339 Expression via Long-Range Loop Formation

Genome-wide association studies (GWASs) have reproducibly associated variants within intergenic regions of 1p36.12 locus with osteoporosis, but the functional roles underlying these noncoding variants are unknown. Through an integrative functional genomic and epigenomic analyses, we prioritized rs6426749 as a potential causal SNP for osteoporosis at 1p36.12. Dual-luciferase assay and CRISPR/Cas9 experiments demonstrate that rs6426749 acts as a distal allele-specific enhancer regulating expression of a lncRNA (LINC00339) (∼360 kb) via long-range chromatin loop formation and that this loop is mediated by CTCF occupied near rs6426749 and LINC00339 promoter region. Specifically, rs6426749-G allele can bind transcription factor TFAP2A, which efficiently elevates the enhancer activity and increases LINC00339 expression. Downregulation of LINC00339 significantly increases the expression of CDC42 in osteoblast cells, which is a pivotal regulator involved in bone metabolism. Our study provides mechanistic insight into how a noncoding SNP affects osteoporosis by long-range interaction, a finding that could indicate promising therapeutic targets for osteoporosis.

Multidimensional Bone Density Phenotyping Reveals New Insights Into Genetic Regulation of the Pediatric Skeleton

Osteoporosis is a complex disease with developmental origins. It is therefore important to understand the genetic contribution to pediatric areal bone mineral density (aBMD). Individual skeletal site phenotyping has been primarily used to identify pediatric aBMD loci. However, this approach is limited because there is a degree of aBMD discordance across skeletal sites. We therefore applied a novel multidimensional phenotyping approach to further understand the genetic regulation of pediatric aBMD. Our sample comprised a prospective, longitudinal cohort of 1293 children of European ancestry (52% female; up to seven annual measurements). Principal components analysis was applied to dual-energy X-ray absorptiometry-derived aBMD Z-scores for total hip, femoral neck, spine, and distal radius to generate multidimensional aBMD phenotypes (ie, principal component scores). We tested the association between a genetic score (percentage of bone lowering alleles at 63 loci) and each principal component. We also performed a genomewide association study (GWAS) using the multiethnic baseline data (n = 1885) to identify novel loci associated with these principal components. The first component (PC1) reflected a concordant phenotypic model of the skeleton (eg, higher loading score = higher BMD across all sites). In contrast, PC2 was discordant for distal radius versus spine and hip aBMD, and PC3 was discordant for spine versus distal radius and hip aBMD. The genetic score was associated with PC1 (beta = -0.05, p = 3.9 × 10^-10 ), but was not associated with discordant PC2 or PC3. Our GWAS discovered variation near CPED1 that associated with PC2 (rs67991850, p = 2.5 × 10^-11 ) and near RAB11FIP5 (rs58649746, p = 4.8 × 10^-9 ) that associated with PC3. In conclusion, an established bone fragility genetic summary score was associated with a concordant skeletal phenotype, but not discordant skeletal phenotypes. Novel associations were observed for the discordant multidimensional skeletal phenotypes that provide new biological insights into the developing skeleton. © 2017 American Society for Bone and Mineral Research.

Screening Gene Knockout Mice for Variation in Bone Mass: Analysis by μCT and Histomorphometry

PURPOSE OF REVIEW

The international mouse phenotyping consortium (IMPC) is producing defined gene knockout mouse lines. Here, a phenotyping program is presented that is based on micro-computed tomography (μCT) assessment of distal femur and vertebra. Lines with significant variation undergo a computer-based bone histomorphometric analysis.

RECENT FINDINGS

Of the 220 lines examined to date, approximately 15% have a significant variation (high or low) by μCT, most of which are not identified by the IMPC screen. Significant dimorphism between the sexes and bone compartments adds to the complexity of the skeletal findings. The μCT information that is posted at www.bonebase.org can group KOMP lines with similar morphological features. The histological data is presented in a graphic form that associates the cellular features with a specific anatomic group. The web portal presents a bone-centric view appropriate for the skeletal biologist/clinician to organize and understand the large number of genes that can influence skeletal health. Cataloging the relative severity of each variant is the first step towards compiling the dataset necessary to appreciate the full polygenic basis of degenerative bone disease.

Differential Expression Profiling of Long Noncoding RNA and mRNA during Osteoblast Differentiation in Mouse

Long noncoding RNAs (lncRNAs) are emerging as an important controller affecting metabolic tissue development, signaling, and function. However, little is known about the function and profile of lncRNAs in osteoblastic differentiation in mice. Here, we analyzed the RNA-sequencing (RNA-Seq) datasets obtained for 18 days in two-day intervals from neonatal mouse calvarial pre-osteoblast-like cells. Over the course of osteoblast differentiation, 4058 mRNAs and 3948 lncRNAs were differentially expressed, and they were grouped into 12 clusters according to the expression pattern by fuzzy c-means clustering. Using weighted gene coexpression network analysis, we identified 9 modules related to the early differentiation stage (days 2–8) and 7 modules related to the late differentiation stage (days 10–18). Gene ontology and KEGG pathway enrichment analysis revealed that the mRNA and lncRNA upregulated in the late differentiation stage are highly associated with osteogenesis. We also identified 72 mRNA and 89 lncRNAs as potential markers including several novel markers for osteoblast differentiation and activation. Our findings provide a valuable resource for mouse lncRNA study and improves our understanding of the biology of osteoblastic differentiation in mice.

Genetically Determined Later Puberty Impacts Lowered Bone Mineral Density in Childhood and Adulthood

Later puberty associates with lower areal bone mineral density (aBMD), and both are risk factors for osteoporosis. However, the association between puberty timing-associated genetic variants and aBMD during development, and the causal relationship between puberty timing and aBMD, remain uncharacterized. We constructed sex-specific polygenic risk scores (GRS) consisting of 333 genetic variants associated with later puberty in European-descent children in the Bone Mineral Density in Childhood Study (BMDCS), consisting of a longitudinal cohort with up to seven assessments (n = 933) and a cross-sectional cohort (n = 486). These GRS were tested for associations with age- and sex-specific aBMD Z-scores at the lumbar spine (LS), femoral neck (FN), total hip, and distal radius, accounting for clinical covariates using sex-stratified linear mixed models. The causal relationship between puberty timing and aBMD was tested in the BMDCS and in publicly available adult data (GEFOS consortium) using two-sample Mendelian randomization (MR). The puberty-delaying GRS was associated with later puberty and lower LS-aBMD in the BMDCS in both sexes (combined beta ± SE = -0.078 ± 0.024; p = 0.0010). In the MR framework, the puberty-delaying genetic instrument also supported a causal association with lower LS-aBMD and FN-aBMD in adults of both sexes. Our results suggest that pubertal timing is causal for diminished aBMD in a skeletal site- and sex-specific manner that tracks throughout life, potentially impacting later risk for osteoporosis, which should be tested in future studies. © 2017 American Society for Bone and Mineral Research.

Evaluation of shared genetic aetiology between osteoarthritis and bone mineral density identifies SMAD3 as a novel osteoarthritis risk locus

Osteoarthritis (OA) is a common complex disease with high public health burden and no curative therapy. High bone mineral density (BMD) is associated with an increased risk of developing OA, suggesting a shared underlying biology. Here, we performed the first systematic overlap analysis of OA and BMD on a genome wide scale. We used summary statistics from the GEFOS consortium for lumbar spine (n = 31,800) and femoral neck (n = 32,961) BMD, and from the arcOGEN consortium for three OA phenotypes (hip, n_cases=3,498; knee, n_cases=3,266; hip and/or knee, n_cases=7,410; n_controls=11,009). Performing LD score regression we found a significant genetic correlation between the combined OA phenotype (hip and/or knee) and lumbar spine BMD (r_g=0.18, P = 2.23 × 10⁻²), which may be driven by the presence of spinal osteophytes. We identified 143 variants with evidence for cross-phenotype association which we took forward for replication in independent large-scale OA datasets, and subsequent meta-analysis with arcOGEN for a total sample size of up to 23,425 cases and 236,814 controls. We found robustly replicating evidence for association with OA at rs12901071 (OR 1.08 95% CI 1.05–1.11, P_meta=3.12 × 10⁻¹⁰), an intronic variant in the SMAD3 gene, which is known to play a role in bone remodeling and cartilage maintenance. We were able to confirm expression of SMAD3 in intact and degraded cartilage of the knee and hip. Our findings provide the first systematic evaluation of pleiotropy between OA and BMD, highlight genes with biological relevance to both traits, and establish a robust new OA genetic risk locus at SMAD3.

Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects

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.

The RSPO3 gene as genetic markers for bone mass assessed by quantitative ultrasound in a population of young adults

Ultrasound bone mass measurement has been postulated as a valuable bone-health assessment tool for primary care. The aim of this study was to analyse the possible relationship between the SPTBN1, RSPO3, CCDC170, DKK1, GPATCH1, and TMEM135 genes, with calcaneal quantitative ultrasound (QUS) in a population of young adults. These genes were first associated with broadband ultrasound attenuation (BUA) in the GEFOS/GENOMOS study. A cross-sectional study was conducted on 575 individuals (mean age 20.41 ± 2.69). Bone mass at the right calcaneus was estimated by QUS. Six single-nucleotide polymorphisms (SNPs) in SPTBN1 (rs11898505), RSPO3 (rs7741021), CCDC170 (rs4869739), DKK1 (rs7902708), TMEM135 (rs597319), and GPATCH1 (rs10416265) were selected as genetic markers based on their previous association with calcaneal QUS. After adjusting for multiple confounding factors, the only significant association with QUS in our population was found for the rs7741021 SNP in the RSPO3 gene (P = 0.006) using the dominant model of inheritance. This suggests the possible implication of the RSPO3 gene in bone mass acquisition during early adulthood.

Increased Ca2+ signaling through CaV1.2 promotes bone formation and prevents estrogen deficiency-induced bone loss

While the prevalence of osteoporosis is growing rapidly with population aging, therapeutic options remain limited. Here, we identify potentially novel roles for CaV1.2 L-type voltage-gated Ca2+ channels in osteogenesis and exploit a transgenic gain-of-function mutant CaV1.2 to stem bone loss in ovariectomized female mice. We show that endogenous CaV1.2 is expressed in developing bone within proliferating chondrocytes and osteoblasts. Using primary BM stromal cell (BMSC) cultures, we found that Ca2+ influx through CaV1.2 activates osteogenic transcriptional programs and promotes mineralization. We used Prx1-, Col2a1-, or Col1a1-Cre drivers to express an inactivation-deficient CaV1.2 mutant in chondrogenic and/or osteogenic precursors in vivo and found that the resulting increased Ca2+ influx markedly thickened bone not only by promoting osteogenesis, but also by inhibiting osteoclast activity through increased osteoprotegerin secretion from osteoblasts. Activating the CaV1.2 mutant in osteoblasts at the time of ovariectomy stemmed bone loss. Together, these data highlight roles for CaV1.2 in bone and demonstrate the potential dual anabolic and anticatabolic therapeutic actions of tissue-specific CaV1.2 activation in osteoblasts.

A molecular roadmap of the AGM region reveals BMPER as a novel regulator of HSC maturation

Through transcriptional profiling of the mouse AGM region, McGarvey et al. identify potential niche regulators of HSC development. They show a new function of BMPER in regulating HSC maturation, likely via its modulation of BMP signalling.

In the developing embryo, hematopoietic stem cells (HSCs) emerge from the aorta-gonad-mesonephros (AGM) region, but the molecular regulation of this process is poorly understood. Recently, the progression from E9.5 to E10.5 and polarity along the dorso-ventral axis have been identified as clear demarcations of the supportive HSC niche. To identify novel secreted regulators of HSC maturation, we performed RNA sequencing over these spatiotemporal transitions in the AGM region and supportive OP9 cell line. Screening several proteins through an ex vivo reaggregate culture system, we identify BMPER as a novel positive regulator of HSC development. We demonstrate that BMPER is associated with BMP signaling inhibition, but is transcriptionally induced by BMP4, suggesting that BMPER contributes to the precise control of BMP activity within the AGM region, enabling the maturation of HSCs within a BMP-negative environment. These findings and the availability of our transcriptional data through an accessible interface should provide insight into the maintenance and potential derivation of HSCs in culture.

Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis

Osteoporosis is a common disease diagnosed primarily by measurement of bone mineral density (BMD). We undertook a genome-wide association study (GWAS) in 142,487 individuals from the UK Biobank to identify loci associated with BMD as estimated by quantitative ultrasound of the heel. We identified 307 conditionally independent single-nucleotide polymorphisms (SNPs) that attained genome-wide significance at 203 loci, explaining approximately 12% of the phenotypic variance. These included 153 previously unreported loci, and several rare variants with large effect sizes. To investigate the underlying mechanisms, we undertook (1) bioinformatic, functional genomic annotation and human osteoblast expression studies; (2) gene-function prediction; (3) skeletal phenotyping of 120 knockout mice with deletions of genes adjacent to lead independent SNPs; and (4) analysis of gene expression in mouse osteoblasts, osteocytes and osteoclasts. The results implicate GPC6 as a novel determinant of BMD, and also identify abnormal skeletal phenotypes in knockout mice associated with a further 100 prioritized genes.

Genetic Associations with Gestational Duration and Spontaneous Preterm Birth

BACKGROUND

Despite evidence that genetic factors contribute to the duration of gestation and the risk of preterm birth, robust associations with genetic variants have not been identified. We used large data sets that included the gestational duration to determine possible genetic associations.

METHODS

We performed a genomewide association study in a discovery set of samples obtained from 43,568 women of European ancestry using gestational duration as a continuous trait and term or preterm (<37 weeks) birth as a dichotomous outcome. We used samples from three Nordic data sets (involving a total of 8643 women) to test for replication of genomic loci that had significant genomewide association (P<5.0×10-8) or an association with suggestive significance (P<1.0×10-6) in the discovery set.

RESULTS

In the discovery and replication data sets, four loci (EBF1, EEFSEC, AGTR2, and WNT4) were significantly associated with gestational duration. Functional analysis showed that an implicated variant in WNT4 alters the binding of the estrogen receptor. The association between variants in ADCY5 and RAP2C and gestational duration had suggestive significance in the discovery set and significant evidence of association in the replication sets; these variants also showed genomewide significance in a joint analysis. Common variants in EBF1, EEFSEC, and AGTR2 showed association with preterm birth with genomewide significance. An analysis of mother-infant dyads suggested that these variants act at the level of the maternal genome.

CONCLUSIONS

In this genomewide association study, we found that variants at the EBF1, EEFSEC, AGTR2, WNT4, ADCY5, and RAP2C loci were associated with gestational duration and variants at the EBF1, EEFSEC, and AGTR2 loci with preterm birth. Previously established roles of these genes in uterine development, maternal nutrition, and vascular control support their mechanistic involvement. (Funded by the March of Dimes and others.).

MECHANISMS IN ENDOCRINOLOGY: Genetics of human bone formation

Throughout life, bone is continuously remodelled to be able to fulfil its multiple functions. The importance of strictly regulating the bone remodelling process, which is defined by the sequential actions of osteoclasts and osteoblasts, is shown by a variety of disorders with abnormalities in bone mass and strength. The best known and most common example of such a disorder is osteoporosis, which is marked by a decreased bone mass and strength that consequently results in an increased fracture risk. As osteoporosis is a serious health problem, a large number of studies focus on elucidating the aetiology of the disease as well as on the identification of novel therapeutic targets for the treatment of osteoporotic patients. These studies have demonstrated that a large amount of variation in bone mass and strength is often influenced by genetic variation in genes encoding important regulators of bone homeostasis. Throughout the years, studies into the genetic causes of osteoporosis as well as several rare monogenic disorders with abnormal high or low bone mass and strength have largely increased the knowledge on regulatory pathways important for bone resorption and formation. This review gives an overview of genes and pathways that are important for the regulation of bone formation and that are identified through their involvement in monogenic and complex disorders with abnormal bone mass. Furthermore, novel bone-forming strategies for the treatment of osteoporosis that resulted from these discoveries, such as antibodies against sclerostin, are discussed as well.

Omics analysis of human bone to identify genes and molecular networks regulating skeletal remodeling in health and disease

The skeleton is a metabolically active organ throughout life where specific bone cell activity and paracrine/endocrine factors regulate its morphogenesis and remodeling. In recent years, an increasing number of reports have used multi-omics technologies to characterize subsets of bone biological molecular networks. The skeleton is affected by primary and secondary disease, lifestyle and many drugs. Therefore, to obtain relevant and reliable data from well characterized patient and control cohorts are vital. Here we provide a brief overview of omics studies performed on human bone, of which our own studies performed on trans-iliacal bone biopsies from postmenopausal women with osteoporosis (OP) and healthy controls are among the first and largest. Most other studies have been performed on smaller groups of patients, undergoing hip replacement for osteoarthritis (OA) or fracture, and without healthy controls. The major findings emerging from the combined studies are: 1. Unstressed and stressed bone show profoundly different gene expression reflecting differences in bone turnover and remodeling and 2. Omics analyses comparing healthy/OP and control/OA cohorts reveal characteristic changes in transcriptomics, epigenomics (DNA methylation), proteomics and metabolomics. These studies, together with genome-wide association studies, in vitro observations and transgenic animal models have identified a number of genes and gene products that act via Wnt and other signaling systems and are highly associated to bone density and fracture. Future challenge is to understand the functional interactions between bone-related molecular networks and their significance in OP and OA pathogenesis, and also how the genomic architecture is affected in health and disease.

Bivariate genome-wide association meta-analysis of pediatric musculoskeletal traits reveals pleiotropic effects at the SREBF1/TOM1L2 locus

Bone mineral density is known to be a heritable, polygenic trait whereas genetic variants contributing to lean mass variation remain largely unknown. We estimated the shared SNP heritability and performed a bivariate GWAS meta-analysis of total-body lean mass (TB-LM) and total-body less head bone mineral density (TBLH-BMD) regions in 10,414 children. The estimated SNP heritability is 43% (95% CI: 34–52%) for TBLH-BMD, and 39% (95% CI: 30–48%) for TB-LM, with a shared genetic component of 43% (95% CI: 29–56%). We identify variants with pleiotropic effects in eight loci, including seven established bone mineral density loci: WNT4, GALNT3, MEPE, CPED1/WNT16, TNFSF11, RIN3, and PPP6R3/LRP5. Variants in the TOM1L2/SREBF1 locus exert opposing effects TB-LM and TBLH-BMD, and have a stronger association with the former trait. We show that SREBF1 is expressed in murine and human osteoblasts, as well as in human muscle tissue. This is the first bivariate GWAS meta-analysis to demonstrate genetic factors with pleiotropic effects on bone mineral density and lean mass.

Bone mineral density and lean skeletal mass are heritable traits. Here, Medina-Gomez and colleagues perform bivariate GWAS analyses of total body lean mass and bone mass density in children, and show genetic loci with pleiotropic effects on both traits.

Osteocyte-specific WNT1 regulates osteoblast function during bone homeostasis

Mutations in WNT1 cause osteogenesis imperfecta (OI) and early-onset osteoporosis, identifying it as a key Wnt ligand in human bone homeostasis. However, how and where WNT1 acts in bone are unclear. To address this mechanism, we generated late-osteoblast-specific and osteocyte-specific WNT1 loss- and gain-of-function mouse models. Deletion of Wnt1 in osteocytes resulted in low bone mass with spontaneous fractures similar to that observed in OI patients. Conversely, Wnt1 overexpression from osteocytes stimulated bone formation by increasing osteoblast number and activity, which was due in part to activation of mTORC1 signaling. While antiresorptive therapy is the mainstay of OI treatment, it has limited efficacy in WNT1-related OI. In this study, anti-sclerostin antibody (Scl-Ab) treatment effectively improved bone mass and dramatically decreased fracture rate in swaying mice, a model of global Wnt1 loss. Collectively, our data suggest that WNT1-related OI and osteoporosis are caused in part by decreased mTORC1-dependent osteoblast function resulting from loss of WNT1 signaling in osteocytes. As such, this work identifies an anabolic function of osteocytes as a source of Wnt in bone development and homoeostasis, complementing their known function as targets of Wnt signaling in regulating osteoclastogenesis. Finally, this study suggests that Scl-Ab is an effective genotype-specific treatment option for WNT1-related OI and osteoporosis.

A Genomewide Association Study Identifies Two Sex-Specific Loci, at SPTB and IZUMO3, Influencing Pediatric Bone Mineral Density at Multiple Skeletal Sites

Failure to achieve optimal bone mineral accretion during childhood and adolescence results in subsequent suboptimal peak bone mass, contributing to osteoporosis risk later in life. To identify novel genetic factors that influence pediatric bone mass at discrete skeletal sites, we performed a sex-stratified genomewide association study of areal bone mineral density (BMD) measured by dual-energy X-ray absorptiometry at the 1/3 distal radius, spine, total hip, and femoral neck in a cohort of 933 healthy European American children. We took forward signals with p < 5 × 10-5 and minor allele frequency (MAF) >5% into an independent cohort of 486 European American children in search of replication. In doing so, we identified five loci that achieved genome wide significance in the combined cohorts (nearest genes: CPED1, IZUMO3, RBFOX1, SPBT, and TBPL2), of which the last four were novel and two were sex-specific (SPTB in females and IZUMO3 in males), with all of them yielding associations that were particularly strong at a specific skeletal site. Annotation of potential regulatory function, expression quantitative trait loci (eQTL) effects and pathway analyses identified several potential target genes at these associated loci. This study highlights the importance of sex-stratified analyses at discrete skeletal sites during the critical period of bone accrual, and identifies novel loci for further functional follow-up to pinpoint key genes and better understand the regulation of bone development in children. © 2017 American Society for Bone and Mineral Research.

Identification of the gene-regulatory landscape in skeletal development and potential links to skeletal regeneration

A class of gene-regulatory elements called enhancers are the main mediators controlling quantitative, temporal and spatial gene expressions. In the course of evolution, the enhancer landscape of higher organisms such as mammals has become quite complex, exerting biological functions precisely and coordinately. In mammalian skeletal development, the master transcription factors Sox9, Runx2 and Sp7/Osterix function primarily through enhancers on the genome to achieve specification and differentiation of skeletal cells. Recently developed genome-scale analyses have shed light on multiple layers of gene regulations, uncovering not only the primary mode of actions of these transcription factors on skeletal enhancers, but also the relation of the epigenetic landscape to three-dimensional chromatin architecture. Here, we review findings on the emerging framework of gene-regulatory networks involved in skeletal development. We further discuss the power of genome-scale analyses to provide new insights into genetic diseases and regenerative medicine in skeletal tissues.

•

Skeletal development is coordinated by master transcription factors.

•

ChIP-seq analyses for the skeletal regulators identified their modes of actions.

•

Analyses of epigenetic landscape features distinct cell types in skeletal tissues.

•

Integrated analyses of the gene regulatory networks link to skeletal regeneration.

Bone Mass and Strength are Significantly Improved in Mice Overexpressing Human WNT16 in Osteocytes

Recently, we demonstrated that osteoblast-specific overexpression of human WNT16 increased both cortical and trabecular bone mass and structure in mice. To further identify the cell-specific role of Wnt16 in bone homeostasis, we created transgenic (TG) mice overexpressing human WNT16 in osteocytes using Dmp1 promoter (Dmp1-hWNT16 TG) on C57BL/6 (B6) background. We analyzed bone phenotypes and serum bone biomarkers, performed gene expression analysis and measured dynamic bone histomorphometry in Dmp1-hWNT16 TG and wild-type (WT) mice. Compared to WT mice, Dmp1-hWNT16 TG mice exhibited significantly higher whole-body, spine and femoral aBMD, BMC and trabecular (BV/TV, Tb.N, and Tb.Th) and cortical (bone area and thickness) parameters in both male and female at 12 weeks of age. Femur stiffness and ultimate force were also significantly improved in the Dmp1-hWNT16 TG female mice, compared to sex-matched WT littermates. In addition, female Dmp1-hWNT16 TG mice displayed significantly higher MS/BS, MAR and BFR/BS compared to the WT mice. Gene expression analysis demonstrated significantly higher mRNA level of Alp in both male and female Dmp1-hWNT16 TG mice and significantly higher levels of Osteocalcin, Opg and Rankl in the male Dmp1-hWNT16 TG mice in bone tissue compared to sex-matched WT mice. These results indicate that WNT16 plays a critical role for acquisition of both cortical and trabecular bone mass and strength. Strategies designed to use WNT16 as a target for therapeutic interventions will be valuable to treat osteoporosis and other low bone mass conditions.

A genetic factor associated with low final bone mineral density in children after a long-term glucocorticoids treatment

Treatment with glucocorticoids is associated with lower bone mineral density (BMD). We performed a genome-wide association study to analyze interactive effects between genotypes and cumulative dose of prednisone (PD) over 4.3 years of follow-up period on the final BMD Z-scores in 461 white children from the Childhood Asthma Management Program. No variants met the conventional criteria for genome-wide significance, and thus we looked for evidence of replication. The top 100-ranked single-nucleotide polymorphisms (SNPs) were then carried forward replication in 59 children with acute lymphoblastic leukemia (ALL) exposed to large fixed doses of PD as part of their chemotherapeutic regimen. Among them, rs6461639 (interaction P=1.88 × 10-5 in the CAMP population) showed a significant association with the final BMD Z-scores in the ALL population (P=0.016). The association of the ALL population was only present after correction for the anti-metabolite treatment arm (high vs low dose). We have identified a novel SNP, rs6461639, showing a significant effect on the final BMD Z-scores in two independent pediatric populations after long-term high-dose PD treatment.

RSPO3-LGR4 Regulates Osteogenic Differentiation Of Human Adipose-Derived Stem Cells Via ERK/FGF Signalling

The four R-spondins (RSPOs) and their three related receptors, LGR4, 5 and 6, have emerged as a major ligand-receptor system with critical roles in development and stem cell survival. However, the exact roles of the RSPO-LGR system in osteogenesis remain largely unknown. In the present study, we showed that RSPO3-shRNA increased the osteogenic potential of human adipose-derived stem cells (hASCs) significantly. Mechanistically, we demonstrated that RSPO3 is a negative regulator of ERK/FGF signalling. We confirmed that inhibition of the ERK1/2 signalling pathway blocked osteogenic differentiation in hASCs, and the increased osteogenic capacity observed after RSPO3 knockdown in hASCs was reversed by inhibition of ERK signalling. Further, silencing of LGR4 inhibited the activity of ERK signalling and osteogenic differentiation of hASCs. Most importantly, we found that loss of LGR4 abrogated RSPO3-regulated osteogenesis and RSPO3-induced ERK1/2 signalling inhibition. Collectively, our data show that ERK signalling works downstream of LGR4 and RSPO3 regulates osteoblastic differentiation of hASCs possibly via the LGR4-ERK signalling.

Integrating GWAS and Co-expression Network Data Identifies Bone Mineral Density Genes SPTBN1 and MARK3 and an Osteoblast Functional Module

Bone mineral density (BMD) is a highly heritable predictor of osteoporotic fracture. Genome-wide association studies (GWAS) for BMD have identified dozens of associations; yet, the genes responsible for most associations remain elusive. Here, we used a bone co-expression network to predict causal genes at BMD GWAS loci based on the premise that genes underlying a disease are often functionally related and functionally related genes are often co-expressed. By mapping genes implicated by BMD GWAS onto a bone co-expression network, we predicted and inferred the function of causal genes for 30 of 64 GWAS loci. We experimentally confirmed that two of the genes predicted to be causal, SPTBN1 and MARK3, are potentially responsible for the effects of GWAS loci on chromosomes 2p16.2 and 14q32.32, respectively. This approach provides a roadmap for the dissection of additional BMD GWAS associations. Furthermore, it should be applicable to GWAS data for a wide range of diseases.

Using Genetic Variation to Explore the Causal Effect of Maternal Pregnancy Adiposity on Future Offspring Adiposity: A Mendelian Randomisation Study

BACKGROUND

It has been suggested that greater maternal adiposity during pregnancy affects lifelong risk of offspring fatness via intrauterine mechanisms. Our aim was to use Mendelian randomisation (MR) to investigate the causal effect of intrauterine exposure to greater maternal body mass index (BMI) on offspring BMI and fat mass from childhood to early adulthood.

METHODS AND FINDINGS

We used maternal genetic variants as instrumental variables (IVs) to test the causal effect of maternal BMI in pregnancy on offspring fatness (BMI and dual-energy X-ray absorptiometry [DXA] determined fat mass index [FMI]) in a MR approach. This was investigated, with repeat measurements, from ages 7 to 18 in the Avon Longitudinal Study of Parents and Children (ALSPAC; n = 2,521 to 3,720 for different ages). We then sought to replicate findings with results for BMI at age 6 in Generation R (n = 2,337 for replication sample; n = 6,057 for total pooled sample). In confounder-adjusted multivariable regression in ALSPAC, a 1 standard deviation (SD, equivalent of 3.7 kg/m2) increase in maternal BMI was associated with a 0.25 SD (95% CI 0.21-0.29) increase in offspring BMI at age 7, with similar results at later ages and when FMI was used as the outcome. A weighted genetic risk score was generated from 32 genetic variants robustly associated with BMI (minimum F-statistic = 45 in ALSPAC). The MR results using this genetic risk score as an IV in ALSPAC were close to the null at all ages (e.g., 0.04 SD (95% CI -0.21-0.30) at age 7 and 0.03 SD (95% CI -0.26-0.32) at age 18 per SD increase in maternal BMI), which was similar when a 97 variant generic risk score was used in ALSPAC. When findings from age 7 in ALSPAC were meta-analysed with those from age 6 in Generation R, the pooled confounder-adjusted multivariable regression association was 0.22 SD (95% CI 0.19-0.25) per SD increase in maternal BMI and the pooled MR effect (pooling the 97 variant score results from ALSPAC with the 32 variant score results from Generation R) was 0.05 SD (95%CI -0.11-0.21) per SD increase in maternal BMI (p-value for difference between the two results = 0.05). A number of sensitivity analyses exploring violation of the MR results supported our main findings. However, power was limited for some of the sensitivity tests and further studies with relevant data on maternal, offspring, and paternal genotype are required to obtain more precise (and unbiased) causal estimates.

CONCLUSIONS

Our findings provide little evidence to support a strong causal intrauterine effect of incrementally greater maternal BMI resulting in greater offspring adiposity.

Quantitative imaging methods in osteoporosis

Osteoporosis is characterized by a decreased bone mass and quality resulting in an increased fracture risk. Quantitative imaging methods are critical in the diagnosis and follow-up of treatment effects in osteoporosis. Prior radiographic vertebral fractures and bone mineral density (BMD) as a quantitative parameter derived from dual-energy X-ray absorptiometry (DXA) are among the strongest known predictors of future osteoporotic fractures. Therefore, current clinical decision making relies heavily on accurate assessment of these imaging features. Further, novel quantitative techniques are being developed to appraise additional characteristics of osteoporosis including three-dimensional bone architecture with quantitative computed tomography (QCT). Dedicated high-resolution (HR) CT equipment is available to enhance image quality. At the other end of the spectrum, by utilizing post-processing techniques such as the trabecular bone score (TBS) information on three-dimensional architecture can be derived from DXA images. Further developments in magnetic resonance imaging (MRI) seem promising to not only capture bone micro-architecture but also characterize processes at the molecular level. This review provides an overview of various quantitative imaging techniques based on different radiological modalities utilized in clinical osteoporosis care and research.

Novel genetic loci underlying human intracranial volume identified through genome-wide association

Intracranial volume reflects the maximally attained brain size during development, and remains stable with loss of tissue in late life. It is highly heritable, but the underlying genes remain largely undetermined. In a genome-wide association study of 32,438 adults, we discovered five previously unknown loci for intracranial volume and confirmed two known signals. Four of the loci were also associated with adult human stature, but these remained associated with intracranial volume after adjusting for height. We found a high genetic correlation with child head circumference (ρgenetic = 0.748), which indicates a similar genetic background and allowed us to identify four additional loci through meta-analysis (Ncombined = 37,345). Variants for intracranial volume were also related to childhood and adult cognitive function, and Parkinson's disease, and were enriched near genes involved in growth pathways, including PI3K-AKT signaling. These findings identify the biological underpinnings of intracranial volume and their link to physiological and pathological traits.

Relationship between A163G osteoprotegerin gene polymorphism and other osteoporosis parameters in Roma and non-Roma postmenopausal women in eastern Slovakia

BACKGROUND

The study was focused on evaluating the possible correlation between biochemical, anthropometric, and genetic indicators of osteoporosis in postmenopausal women. The frequency of genotypes and differences in measured parameters were evaluated within two ethnically different groups of women in Slovakia.

METHODS

The study included 310 postmenopausal women divided into non-Roma and Roma groups. Based on results of densitometry, they were divided into control groups and women with osteoporosis and osteopenia. In all women, a genetic analysis of polymorphism of osteoprotegerin gene promotor region (A163G) was provided along with measurement of indicators of bone tissue metabolism.

RESULTS

There is a particularly low incidence of osteoporosis in Roma women. We found a correlation between bone mineral density (BMD), body mass index, and waist and hip circumference in women with osteoporosis and in Roma women with osteopenia. The frequency of the AG genotype was higher in non-Roma women with osteoporosis, but reached only 10.7% in Roma women with osteopenia. While the presence of the G allele in the non-Roma population was accompanied by higher BMD and markers of osteoformation, it was accompanied by significantly higher concentrations of parathyroid hormone in the Roma population.

CONCLUSION

The presence of the AG genotype has a different effect on bone metabolism in two ethnically diverse populations of women in Slovakia. In the general population, the presence of the G allele exhibited protective effects consistent with other studies, but in Roma population this appears to be the allele A. However, this requires a further study for confirmation and more detailed characterization of the differences between populations that have this work indicated.

An Emerging Regulatory Landscape for Skeletal Development

Skeletal development creates the physical framework that shapes our body and its actions. In the past two decades, genetic studies have provided important insights into the molecular processes at play, including the roles of signaling pathways and transcriptional effectors that coordinate an orderly, progressive emergence and expansion of distinct cartilage and bone cell fates in an invariant temporal and spatial pattern for any given skeletal element within that specific vertebrate species. Genome-scale studies have provided additional layers of understanding, moving from individual genes to the gene regulatory landscape, integrating regulatory information through cis-regulatory modules into cell type-specific gene regulatory programs. This review discusses our current understanding of the transcriptional control of mammalian skeletal development, focusing on recent genome-scale studies.

How rare bone diseases have informed our knowledge of complex diseases

Rare bone diseases, generally defined as monogenic traits with either autosomal recessive or dominant patterns of inheritance, have provided a rich database of genes and associated pathways over the past 2-3 decades. The molecular genetic dissection of these bone diseases has yielded some major surprises in terms of the causal genes and/or involved pathways. The discovery of genes/pathways involved in diseases such as osteopetrosis, osteosclerosis, osteogenesis imperfecta and many other rare bone diseases have all accelerated our understanding of complex traits. Importantly these discoveries have provided either direct validation for a specific gene embedded in a group of genes within an interval identified through a complex trait genome-wide association study (GWAS) or based upon the pathway associated with a monogenic trait gene, provided a means to prioritize a large number of genes for functional validation studies. In some instances GWAS studies have yielded candidate genes that fall within linkage intervals associated with monogenic traits and resulted in the identification of causal mutations in those rare diseases. Driving all of this discovery is a complement of technologies such as genome sequencing, bioinformatics and advanced statistical analysis methods that have accelerated genetic dissection and greatly reduced the cost. Thus, rare bone disorders in partnership with GWAS have brought us to the brink of a new era of personalized genomic medicine in which the prevention and management of complex diseases will be driven by the molecular understanding of each individuals contributing genetic risks for disease.

Genetics of pediatric bone strength

Osteoporosis is one of the most common chronic forms of disability in postmenopausal women and represents a major health burden around the world. Bone fragility is affected by bone mineral density (BMD), and, one of the most important factors in preventing osteoporosis is optimizing peak bone mass, which is achieved during growth in childhood and adolescence. BMD is a complex trait resulting from environmental and genetic factors. Genome-wide association studies have discovered robust genetic signals influencing BMD in adults, and similar studies have also been conducted to investigate the genetics of BMD in the pediatric setting. These latter studies have revealed that many adult osteoporosis-related loci also regulate BMD during growth. These investigations have the potential to profoundly impact public health and will allow for the eventual development of effective interventions for the prevention of osteoporosis.

Network Analysis Implicates Alpha-Synuclein (Snca) in the Regulation of Ovariectomy-Induced Bone Loss

The postmenopausal period in women is associated with decreased circulating estrogen levels, which accelerate bone loss and increase the risk of fracture. Here, we gained novel insight into the molecular mechanisms mediating bone loss in ovariectomized (OVX) mice, a model of human menopause, using co-expression network analysis. Specifically, we generated a co-expression network consisting of 53 gene modules using expression profiles from intact and OVX mice from a panel of inbred strains. The expression of four modules was altered by OVX, including module 23 whose expression was decreased by OVX across all strains. Module 23 was enriched for genes involved in the response to oxidative stress, a process known to be involved in OVX-induced bone loss. Additionally, module 23 homologs were co-expressed in human bone marrow. Alpha synuclein (Snca) was one of the most highly connected “hub” genes in module 23. We characterized mice deficient in Snca and observed a 40% reduction in OVX-induced bone loss. Furthermore, protection was associated with the altered expression of specific network modules, including module 23. In summary, the results of this study suggest that Snca regulates bone network homeostasis and ovariectomy-induced bone loss.

The case for genome-wide association studies of bone acquisition in paediatric and adolescent populations

Peak bone mass, the maximum amount of bone accrued at the end of the growth period, is an important predictor of future risk of osteoporosis and fracture. Hence, the contribution of genetic factors influencing bone accrual is of considerable interest to the osteoporosis research community. In this article, we review evidence that genetic factors play an important role in bone growth, describe the genetic loci implicated so far and briefly discuss lessons learned from the application of genome-wide association studies. Moreover, we attempt to make the case for genetic investigations of bone mineral density in paediatric and young adult populations, describing their potential to increase our knowledge of the process of bone metabolism throughout the life course, and in turn, identify novel targets for the pharmacological treatment of osteoporosis.

Using Mendelian randomization to investigate a possible causal relationship between adiposity and increased bone mineral density at different skeletal sites in children

BACKGROUND

Lean mass is positively associated with bone mineral density (BMD). However, the relationship between adiposity and BMD is more controversial. In particular, it is unclear if the observational association between the two reflects a causal effect of fat mass on BMD. Previous Mendelian randomization (MR) studies using variants in the FTO and MC4R genes as genetic instruments for adiposity have suggested that fat mass does indeed causally influence BMD. However, it is possible that these genetic variants pleiotropically influence lean mass and affect BMD through pathways independent of adiposity, invalidating one of the core assumptions of MR and complicating interpretation of the analysis.

METHODS

To investigate whether adiposity causally affects BMD, we investigated the relationship between fat mass and BMD at the skull (SK), upper limbs (UL) and lower limbs (LL), spine (SP) and pelvis (PE), using 32 body mass index (BMI)-associated SNPs, including a variant near ADCY3 that was strongly associated with fat but not lean mass in our sample. Dual-energy X-ray absorptiometry (DXA) scans and genetic data were available for 5221 subjects (mean age 9.9 years) from the Avon Longitudinal Study of Parents and Children. We performed a series of MR analyses involving single BMI-associated SNPs and allelic scores of these SNPs. We used new extensions of the MR method including MR Egger regression and multivariable MR, which are more robust to possible confounding effects due to horizontal pleiotropy and, in the case of multivariable MR, specifically account for the effect of lean mass in the analysis. Bidirectional Mendelian randomization analysis was also performed to examine whether BMD causally affected BMI and adiposity.

RESULTS

Observationally, fat mass was strongly positively related to BMD at all sites, but more weakly at the skull. Instrumental variables (IV) analyses using an allelic score of BMI SNPs suggested that fat mass was causally related to LL-BMD, UL-BMD, SP-BMD and PE-BMD but not SK-BMD. Multivariable MR, Egger regression and IV analyses involving the ADCY3 variant suggested a positive causal effect of adiposity on all sites except the skull, and that an effect was present even after taking lean mass into account. Finally, IV analyses using BMD allelic scores showed no evidence of reverse causality between BMD and fat mass.

CONCLUSIONS

Our results suggest that adiposity is causally related to increased BMD at all sites except the skull, perhaps reflecting positive effects of loading on bone formation at weighted but not unweighted sites. In contrast, we found no evidence for BMD causally affecting BMI or measures of adiposity. Our results illustrate how MR can be used profitably to investigate clinical questions relevant to osteoporosis.

Cdx-2 polymorphism in Vitamin D Receptor gene was associated with serum 25-hydroxyvitamin D levels, bone mineral density and fracture in middle-aged and elderly Chinese women

The aim of the current study was to examine the relationship between Cdx-2 polymorphism in the promoter region of the VDR gene and serum 25-hydroxyvitamin D (25(OH)D) levels, bone mineral density (BMD) and fracture in Chinese population. This was a cross-sectional study, which included 738 individuals (428 women and 310 men) aged 45 years or older. In women, the association of Cdx-2 polymorphism with serum 25(OH)D levels was significant adjusting for age, BMI, estimated glomerular filtration rate, menopausal status and season of blood collection (P = 0.002). Cdx-2 polymorphism was associated with lumbar spine BMD adjusted for age, BMI, menopausal status and serum 25(OH)D in women (P = 0.005). But it was not associated with femoral neck BMD or total hip BMD in women. In women, Cdx-2 polymorphism was also associated with fracture adjusted for age, BMI, menopausal status, serum 25(OH)D and total hip BMD (P = 0.03). Carriers of AA and AG genotypes was associated with a higher odds of fracture compared with the carriers of GG genotype (OR = 2.14, 95% CI 1.04-4.42 and OR = 1.90, 95% CI 1.03-3.51). In men, Cdx-2 polymorphism was not associated with serum 25(OH)D levels, BMD or fracture. Our results indicate that the association of Cdx-2 polymorphism in the VDR gene with serum 25(OH)D levels, BMD and fracture may have sex differences. Cdx-2 polymorphism in the VDR gene may affect the serum 25(OH)D concentrations and the risk of osteoporosis and fracture in middle-aged and elderly Chinese women.

Integrative Analysis of Genomics and Transcriptome Data to Identify Potential Functional Genes of BMDs in Females

Osteoporosis is known to be highly heritable. However, to date, the findings from more than 20 genome-wide association studies (GWASs) have explained less than 6% of genetic risks. Studies suggest that the missing heritability data may be because of joint effects among genes. To identify novel heritability for osteoporosis, we performed a system-level study on bone mineral density (BMD) by weighted gene coexpression network analysis (WGCNA), using the largest GWAS data set for BMD in the field, Genetic Factors for Osteoporosis Consortium (GEFOS-2), and a transcriptomic gene expression data set generated from transiliac bone biopsies in women. A weighted gene coexpression network was generated for 1574 genes with GWAS nominal evidence of association (p ≤ 0.05) based on dissimilarity measurement on the expression data. Twelve distinct gene modules were identified, and four modules showed nominally significant associations with BMD (p ≤ 0.05), but only one module, the yellow module, demonstrated a good correlation between module membership (MM) and gene significance (GS), suggesting that the yellow module serves an important biological role in bone regulation. Interestingly, through characterization of module content and topology, the yellow module was found to be significantly enriched with contractile fiber part (GO:044449), which is widely recognized as having a close relationship between muscle and bone. Furthermore, detailed submodule analyses of important candidate genes (HOMER1, SPTBN1) by all edges within the yellow module implied significant enrichment of functional connections between bone and cytoskeletal protein binding. Our study yielded novel information from system genetics analyses of GWAS data jointly with transcriptomic data. The findings highlighted a module and several genes in the model as playing important roles in the regulation of bone mass in females, which may yield novel insights into the genetic basis of osteoporosis. © 2016 American Society for Bone and Mineral Research.

Osteoporosis and Bone Mass Disorders: From Gene Pathways to Treatments

Genomic technologies have evolved rapidly contributing to the understanding of diseases. Genome-wide association studies (GWAS) and whole-exome sequencing have aided the identification of the genetic determinants of monogenic and complex conditions including osteoporosis and bone mass disorders. Overlap exists between the genes implicated in monogenic and complex forms of bone mass disorders, largely explained by the clustering of genes encoding factors in signaling pathways crucial for mesenchymal cell differentiation, skeletal development, and bone remodeling and metabolism. Numerous of the remaining discovered genes merit functional follow-up studies to elucidate their role in bone biology. The insight provided by genetic studies is serving the identification of biomarkers predictive of disease, redefining disease, response to treatment, and discovery of novel drug targets for skeletal disorders.

Bone Mass and Strength in School-Age Children Exhibit Sexual Dimorphism Related to Differences in Lean Mass: The Generation R Study

Bone strength, a key determinant of fracture risk, has been shown to display clear sexual dimorphism after puberty. We sought to determine whether sex differences in bone mass and hip bone geometry as an index of strength exist in school-age prepubertal children and the degree to which the differences are independent of body size and lean mass. We studied 3514 children whose whole-body and hip scans were measured using the same densitometer (GE-Lunar iDXA) at a mean age of 6.2 years. Hip dual-energy X-ray absorptiometry (DXA) scans underwent hip structural analyses (HSA) with derivation of bone strength indices. Sex differences in these parameters were assessed by regression models adjusted for age, height, ethnicity, weight, and lean mass fraction (LMF). Whole-body bone mineral density (BMD) and bone mineral content (BMC) levels were 1.3% and 4.3% higher in girls after adjustment by LMF. Independent of LMF, boys had 1.5% shorter femurs, 1.9% and 2.2% narrower shaft and femoral neck with 1.6% to 3.4% thicker cortices than girls. Consequent with this geometry configuration, girls observed 6.6% higher stresses in the medial femoral neck than boys. When considering LMF, the sexual differences on the derived bone strength indices were attenuated, suggesting that differences in muscle loads may reflect an innate disadvantage in bone strength in girls, as consequence of their lower muscular acquisition. In summary, we show that bone sexual dimorphism is already present at 6 years of age, with boys having stronger bones than girls, the relation of which is influenced by body composition and likely attributable to differential adaptation to mechanical loading. Our results support the view that early life interventions (ie, increased physical activity) targeted during the pre- and peripubertal stages may be of high importance, particularly in girls, because before puberty onset, muscle mass is strongly associated with bone density and geometry in children. © 2015 American Society for Bone and Mineral Research.

SNP- and haplotype-based genome-wide association studies for growth, carcass, and meat quality traits in a Duroc multigenerational population

Background

The aim of the present study was to compare the power of single nucleotide polymorphism (SNP)-based genome-wide association study (GWAS) and haplotype-based GWAS for quantitative trait loci (QTL) detection, and to detect novel candidate genes affecting economically important traits in a purebred Duroc population comprising seven-generation pedigree. First, we performed a simulation analysis using real genotype data of this population to compare the power (based on the null hypothesis) of the two methods. We then performed GWAS using both methods and real phenotype data comprising 52 traits, which included growth, carcass, and meat quality traits.

Results

In total, 836 animals were genotyped using the Illumina PorcineSNP60 BeadChip and 14 customized SNPs from regions of known candidate genes related to the traits of interest. The power of SNP-based GWAS was greater than that of haplotype-based GWAS in a simulation analysis. In real data analysis, a larger number of significant regions was obtained by SNP-based GWAS than by haplotype-based GWAS. For SNP-based GWAS, 23 genome-wide significant SNP regions were detected for 17 traits, and 120 genome-wide suggestive SNP regions were detected for 27 traits. For haplotype-based GWAS, 6 genome-wide significant SNP regions were detected for four traits, and 11 genome-wide suggestive SNP regions were detected for eight traits. All genome-wide significant SNP regions detected by haplotype-based GWAS were located in regions also detected by SNP-based GWAS. Four regions detected by SNP-based GWAS were significantly associated with multiple traits: on Sus scrofa chromosome (SSC) 1 at 304 Mb; and on SSC7 at 35–39 Mb, 41–42 Mb, and 103 Mb. The vertnin gene (VRTN) in particular, was located on SSC7 at 103 Mb and was significantly associated with vertebrae number and carcass lengths. Mapped QTL regions contain some candidate genes involved in skeletal formation (FUBP3; far upstream element binding protein 3) and fat deposition (METTL3; methyltransferase like 3).

Conclusion

Our results show that a multigenerational pig population is useful for detecting QTL, which are typically segregated in a purebred population. In addition, a novel significant region could be detected by SNP-based GWAS as opposed to haplotype-based GWAS.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0368-3) contains supplementary material, which is available to authorized users.