The heritability of bone mineral density, ultrasound of the calcaneus and hip axis length: a study of postmenopausal twins

The chromatin remodeling factor Arid1a cooperates with Jun/Fos to promote osteoclastogenesis by epigenetically upregulating Siglec15 expression

Osteoporosis is characterized by an imbalance between osteoclast-mediated bone resorption and osteoblast-related bone formation, particularly increased osteoclastogenesis. However, the mechanisms by which epigenetic factors regulate osteoclast precursor differentiation during osteoclastogenesis remain poorly understood. Here, we show that the specific knockout of the chromatin remodeling factor Arid1a in bone marrow-derived macrophages (BMDMs) results in increased bone mass. The loss of Arid1a in BMDM inhibits cell-cell fusion and maturation of osteoclast precursors, thereby suppressing osteoclast differentiation. Mechanistically, Arid1a increases the chromatin access in the gene promoter region of sialic acid-binding Ig-like lectin 15 (Siglec15) by transcription factor Jun/Fos, which results in the upregulation of Siglec15 and promotion of osteoclast differentiation. However, the loss of Arid1a reprograms the chromatin structure to restrict Siglec15 expression in osteoclast precursors, thereby inhibiting BMDM differentiation into mature osteoclasts. Deleting Arid1a after ovariectomy (a model for postmenopausal bone loss) alleviated bone loss and maintained bone mass. In summary, epigenetic reprogramming mediated by Arid1a loss suppresses osteoclast differentiation and may serve as a promising therapeutic strategy for treating bone loss diseases.

Genetic and environmental determinants of bone quality: a cross-sectional analysis of the Hungarian Twin Registry

There is abundant evidence that bone mineral content is highly heritable, while the heritability of bone quality (i.e. trabecular bone score [TBS] and quantitative ultrasound index [QUI]) is rarely investigated. We aimed to disentangle the role of genetic, shared and unique environmental factors on TBS and QUI among Hungarian twins. Our study includes 82 twin (48 monozygotic, 33 same-sex dizygotic) pairs from the Hungarian Twin Registry. TBS was determined by DXA, QUI by calcaneal bone ultrasound. To estimate the genetic and environmental effects, we utilized ACE-variance decomposition. For the unadjusted model of TBS, an AE model provided the best fit with > 80% additive genetic heritability. Adjustment for age, sex, BMI and smoking status improved model fit with 48.0% of total variance explained by independent variables. Furthermore, there was a strong dominant genetic effect (73.7%). In contrast, unadjusted and adjusted models for QUI showed an AE structure. Adjustments improved model fit and 25.7% of the total variance was explained by independent variables. Altogether 70-90% of the variance in QUI was related to additive genetic influences. We found a strong genetic heritability of bone quality in unadjusted models. Half of the variance of TBS was explained by age, sex and BMI. Furthermore, the adjusted model suggested that the genetic component of TBS could be dominant or an epistasis could be present. In contrast, independent variables explained only a quarter of the variance of QUI and the additive heritability explained more than half of all the variance.

Identification of Transcripts with Shared Roles in the Pathogenesis of Postmenopausal Osteoporosis and Cardiovascular Disease

Epidemiological evidence suggests existing comorbidity between postmenopausal osteoporosis (OP) and cardiovascular disease (CVD), but identification of possible shared genes is lacking. The skeletal global transcriptomes were analyzed in trans-iliac bone biopsies (n = 84) from clinically well-characterized postmenopausal women (50 to 86 years) without clinical CVD using microchips and RNA sequencing. One thousand transcripts highly correlated with areal bone mineral density (aBMD) were further analyzed using bioinformatics, and common genes overlapping with CVD and associated biological mechanisms, pathways and functions were identified. Fifty genes (45 mRNAs, 5 miRNAs) were discovered with established roles in oxidative stress, inflammatory response, endothelial function, fibrosis, dyslipidemia and osteoblastogenesis/calcification. These pleiotropic genes with possible CVD comorbidity functions were also present in transcriptomes of microvascular endothelial cells and cardiomyocytes and were differentially expressed between healthy and osteoporotic women with fragility fractures. The results were supported by a genetic pleiotropy-informed conditional False Discovery Rate approach identifying any overlap in single nucleotide polymorphisms (SNPs) within several genes encoding aBMD- and CVD-associated transcripts. The study provides transcriptional and genomic evidence for genes of importance for both BMD regulation and CVD risk in a large collection of postmenopausal bone biopsies. Most of the transcripts identified in the CVD risk categories have no previously recognized roles in OP pathogenesis and provide novel avenues for exploring the mechanistic basis for the biological association between CVD and OP.

Identification of an novel genetic variant associated with osteoporosis: insights from the Taiwan Biobank Study

Purpose

The purpose of this study was to identify new independent significant SNPs associated with osteoporosis using data from the Taiwan Biobank (TWBB).

Material and Methods

The dataset was divided into discovery (60%) and replication (40%) subsets. Following data quality control, genome-wide association study (GWAS) analysis was performed, adjusting for sex, age, and the top 5 principal components, employing the Scalable and Accurate Implementation of the Generalized mixed model approach. This was followed by a meta-analysis of TWBB1 and TWBB2. The Functional Mapping and Annotation (FUMA) platform was used to identify osteoporosis-associated loci. Manhattan and quantile-quantile plots were generated using the FUMA platform to visualize the results. Independent significant SNPs were selected based on genome-wide significance (P < 5 × 10-8) and independence from each other (r2 < 0.6) within a 1 Mb window. Positional, eQTL(expression quantitative trait locus), and Chromatin interaction mapping were used to map SNPs to genes.

Results

A total of 29 084 individuals (3154 osteoporosis cases and 25 930 controls) were used for GWAS analysis (TWBB1 data), and 18 918 individuals (1917 cases and 17 001 controls) were utilized for replication studies (TWBB2 data). We identified a new independent significant SNP for osteoporosis in TWBB1, with the lead SNP rs76140829 (minor allele frequency = 0.055, P-value = 1.15 × 10-08). Replication of the association was performed in TWBB2, yielding a P-value of 6.56 × 10-3. The meta-analysis of TWBB1 and TWBB2 data demonstrated a highly significant association for SNP rs76140829 (P-value = 7.52 × 10-10). In the positional mapping of rs76140829, 6 genes (HABP2, RP11-481H12.1, RNU7-165P, RP11-139 K1.2, RP11-57H14.3, and RP11-214 N15.5) were identified through chromatin interaction mapping in mesenchymal stem cells.

Conclusions

Our GWAS analysis using the Taiwan Biobank dataset unveils rs76140829 in the VTI1A gene as a key risk variant associated with osteoporosis. This finding expands our understanding of the genetic basis of osteoporosis and highlights the potential regulatory role of this SNP in mesenchymal stem cells.

A longitudinal genome-wide association study of bone mineral density mean and variability in the UK Biobank

Bone mineral density (BMD) is an essential predictor of osteoporosis and fracture. We conducted a genome-wide trajectory analysis of BMD and analyzed the BMD change.

PURPOSE

This study aimed to identify the genetic architecture and potential biomarkers of BMD.

METHODS

Our analysis included 141,261 white participants from the UK Biobank with heel BMD phenotype data. We used a genome-wide trajectory analysis tool, TrajGWAS, to conduct a genome-wide association study (GWAS) of BMD. Then, we validated our findings in previously reported BMD genetic associations and performed replication analysis in the Asian participants. Finally, gene-set enrichment analysis (GSEA) of the identified candidate genes was conducted using the FUMA platform.

RESULTS

A total of 52 genes associated with BMD trajectory mean were identified, of which the top three significant genes were WNT16 (P = 1.31 × 10-126), FAM3C (P = 4.18 × 10-108), and CPED1 (P = 8.48 × 10-106). In addition, 114 genes associated with BMD within-subject variability were also identified, such as AC092079.1 (P = 2.72 × 10-13) and RGS7 (P = 4.72 × 10-10). The associations for these candidate genes were confirmed in the previous GWASs and replicated successfully in the Asian participants. GSEA results of BMD change identified multiple GO terms related to skeletal development, such as SKELETAL SYSTEM DEVELOPMENT (Padjusted = 2.45 × 10-3) and REGULATION OF OSSIFICATION (Padjusted = 2.45 × 10-3). KEGG enrichment analysis showed that these genes were mainly enriched in WNT SIGNALING PATHWAY.

CONCLUSIONS

Our findings indicated that the CPED1-WNT16-FAM3C locus plays a significant role in BMD mean trajectories and identified several novel candidate genes contributing to BMD within-subject variability, facilitating the understanding of the genetic architecture of BMD.

Investigating the association between serum ADAM/ADAMTS levels and bone mineral density by mendelian randomization study

PURPOSE

A Disintegrin and Metalloproteinase (ADAM) and A Disintegrin and Metalloproteinase with Thrombospondin Motif (ADAMTS) have been reported potentially involved in bone metabolism and related to bone mineral density. This Mendelian Randomization (MR) analysis was performed to determine whether there are causal associations of serum ADAM/ADAMTS with BMD in rid of confounders.

METHODS

The genome-wide summary statistics of four site-specific BMD measurements were obtained from studies in individuals of European ancestry, including forearm (n = 8,143), femoral neck (n = 32,735), lumbar spine (n = 28,498) and heel (n = 426,824). The genetic instrumental variables for circulating levels of ADAM12, ADAM19, ADAM23, ADAMTS5 and ADAMTS6 were retrieved from the latest genome-wide association study of European ancestry (n = 5336 ~ 5367). The estimated causal effect was given by the Wald ratio for each variant, the inverse-variance weighted model was used as the primary approach to combine estimates from multiple instruments, and sensitivity analyses were conducted to assess the robustness of MR results. The Bonferroni-corrected significance was set at P < 0.0025 to account for multiple testing, and a lenient threshold P < 0.05 was considered to suggest a causal relationship.

RESULTS

The causal effects of genetically predicted serum ADAM/ADAMTS levels on BMD measurements at forearm, femoral neck and lumbar spine were not statistically supported by MR analyses. Although causal effect of ADAMTS5 on heel BMD given by the primary MR analysis (β = -0.006, -0.010 to 0.002, P = 0.004) failed to reach Bonferroni-corrected significance, additional MR approaches and sensitivity analyses indicated a robust causal relationship.

CONCLUSION

Our study provided suggestive evidence for the causal effect of higher serum levels of ADAMTS5 on decreased heel BMD, while there was no supportive evidence for the associations of ADAM12, ADAM19, ADAM23, and ADAMTS6 with BMD at forearm, femoral neck and lumbar spine in Europeans.

UBAP2 plays a role in bone homeostasis through the regulation of osteoblastogenesis and osteoclastogenesis

Osteoporosis is a condition characterized by decreased bone mineral density (BMD) and reduced bone strength, leading to an increased risk of fractures. Here, to identify novel risk variants for susceptibility to osteoporosis-related traits, an exome-wide association study is performed with 6,485 exonic single nucleotide polymorphisms (SNPs) in 2,666 women of two Korean study cohorts. The rs2781 SNP in UBAP2 gene is suggestively associated with osteoporosis and BMD with p-values of 6.1 × 10-7 (odds ratio = 1.72) and 1.1 × 10-7 in the case-control and quantitative analyzes, respectively. Knockdown of Ubap2 in mouse cells decreases osteoblastogenesis and increases osteoclastogenesis, and knockdown of ubap2 in zebrafish reveals abnormal bone formation. Ubap2 expression is associated with E-cadherin (Cdh1) and Fra1 (Fosl1) expression in the osteclastogenesis-induced monocytes. UBAP2 mRNA levels are significantly reduced in bone marrow, but increased in peripheral blood, from women with osteoporosis compared to controls. UBAP2 protein level is correlated with the blood plasma level of the representative osteoporosis biomarker osteocalcin. These results suggest that UBAP2 has a critical role in bone homeostasis through the regulation of bone remodeling.

Exome-wide screening identifies novel rare risk variants for bone mineral density

Bone mineral density (BMD) is an independent risk factor of osteoporosis-related fractures. We performed gene-based burden tests to assess the association between rare variants and BMD, and identified several BMD candidate genes.

PURPOSE

BMD is highly heritable and a major predictor of osteoporotic fractures, but its genetic basis remains unclear. We aimed to identify rare risk variants contributing to BMD.

METHODS

Utilizing the newly released UK Biobank 200,643 exome dataset, we conducted a gene-based exome-wide association study in males and females, respectively. First, 100,639 males and 117,338 females with BMD values were included in the polygenic risk scores (PRS) analysis. Among individuals with lower 30% PRS, cases were individuals with top 10% BMD, and individuals with bottom 10% BMD were the controls. Considering the effects of vitamin D (VD), individuals with the highest 30% VD concentration were selected for VD-BMD analysis. After quality control, 741 males and 697 females were included in the BMD analysis, and 717 males and 708 females were included in the VD-BMD analysis. The variants were annotated by ANNOVAR software, then BMD and VD-BMD qualified variants were imported into the SKAT R-package to perform gene-based burden tests, respectively.

RESULTS

The gene-based burden test of the exonic variants identified genome-wide candidate associations in ANKRD18A (P = 1.60 × 10^-5, P_{Bonferroni adjust} = 2.11 × 10^-3), C22orf31 (P = 3.49 × 10^-4, P_{Bonferroni adjust} = 3.17 × 10^-2), and SPATC1L (P = 1.09 × 10^-5, P_{Bonferroni adjust} = 8.80 × 10^-3). For VD-BMD analysis, three genes were associated with BMD, such as NIPAL1 (P = 1.06 × 10^-3, P_{Bonferroni adjust} = 3.91 × 10^-2).

CONCLUSIONS

Our study suggested that rare variants contribute to BMD, providing new sights for broadening the genetic structure of BMD.

Lack of causal association between heart failure and osteoporosis: a Mendelian randomization study

Objectives

Heart failure (HF) has been implicated in osteoporosis. However, causality remains unestablished. Here, we sought to assess causal associations of genetic liability to HF with osteoporosis using Mendelian randomization (MR) analyses.

Methods

Independent single nucleotide polymorphisms associated with HF at genome-wide significance were derived from a large genome-wide association study (GWAS) (including up to 977,323 individuals). We obtained summary statistics for forearm (FA) bone mineral density (BMD) (n = 8,143), femoral neck (FN) BMD (n = 32,735), lumbar spine (LS) BMD (n = 28,498), heel (HE) BMD (n = 426,824), and fracture (n = 1,214,434) from other GWAS meta-analyses. Inverse variance weighted (IVW) and several supplementary methods were performed to calculate the MR estimates.

Results

Genetically determined HF has no causal effect on FA-BMD (odds ratio (OR) 1.17; 95% confidence interval (CI) 0.82, 1.66; P = 0.389), FN-BMD (OR 1.01; 95% CI 0.85, 1.19; P = 0.936), LS-BMD (OR 0.96; 95% CI 0.80, 1.17; P = 0.705), HE-BMD (OR 1.01; 95% CI 0.90, 1.13; P = 0.884), and fracture risk (OR 1.00; 95% CI 0.92, 1.10; P = 0.927). Complementary analyses returned broadly consistent results.

Conclusion

This MR study provides genetic evidence that HF may not lead to an increased risk of reduced BMDs or fracture.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01385-8.

Does Aging Activate T-cells to Reduce Bone Mass and Quality?

PURPOSE OF REVIEW

Aging leads to decline in bone mass and quality starting at age 30 in humans. All mammals undergo a basal age-dependent decline in bone mass. Osteoporosis is characterized by low bone mass and changes in bone microarchitecture that increases the risk of fracture. About a third of men over the age of 50 years are osteoporotic because they have higher than basal bone loss. In women, there is an additional acute decrement in bone mass, atop the basal rate, associated with loss of ovarian function (menopause) causing osteoporosis in about half of the women. Both genetics and environmental factors such as smoking, chronic infections, diet, microbiome, and metabolic disease can modulate basal age-dependent bone loss and eventual osteoporosis. Here, we review recent studies on the etiology of age-dependent decline in bone mass and propose a mechanism that integrates both genetic and environmental factors.

RECENT FINDINGS

Recent findings support that aging and menopause dysregulate the immune system leading to sterile low-grade inflammation. Both animal models and human studies demonstrate that certain kinds of inflammation, in both men and women, mediate bone loss. Senolytics, meant to block a wide array of age-induced effects by preventing cellular senescence, have been shown to improve bone mass in aged mice. Based on a synthesis of the recent data, we propose that aging activates long-lived tissue resident memory T-cells to become senescent and proinflammatory, leading to bone loss. Targeting this population may represent a promising osteoporosis therapy. Emerging data indicates that there are several mechanisms that lead to sterile low-grade chronic inflammation, inflammaging, that cause age- and estrogen-loss dependent osteoporosis in men and women.

The Impact of Environmental and Genetic Factors on Bone Quality in Monozygotic and Dizygotic Twins

The purpose of the research was to assess the genetic and environmental influences on bone properties. One hundred thirty-two pairs of twins (99/33 monozygotic/dizygotic) underwent anthropometric measurements and phalangeal quantitative ultrasound (DBM Sonic 1200, Igea, Italy) measuring the amplitude speed of sound (AD-SoS, m/s). The mean age was 16.78 ± 12.35 years for monozygotic twins and 14.30 ± 8 years for dizygotic. Interpair and intrapair correlations between twins were calculated. In the groups of monozygotic and dizygotic twins, Ad-SoS correlated significantly with age (r = 0.56−0.73, p < 0.05), weight (r = 0.73−0.78, p < 0.05), and height (r = 0.80−0.81, p < 0.05). The strongest intrapair correlation (r = 0.99−0.998) was noted in monozygotic females for Ad-SoS, weight, and height. There was a statistically significant correlation between the intrapair difference of Ad-SoS and age but only in the groups of monozygotic and dizygotic females (r = 0.281, r2 = 0.079, and p = 0.028; r = 0.544, r2 = 0.296, and p = 0.01, respectively). After age adjustment, it was estimated that 28.62% of Ad-SoS in women and 13.2% of Ad-SoS in men was explained by genetic influence, leading to the conclusion that Ad-SoS changed with age, weight, and height. The strongest correlation between pairs of twins was observed in monozygotic twins. The differences in bone values between female twins arose with age, which indicated the role of environmental factors.

Genetic Risk for Osteoporosis and the Benefit of Adherence to Healthy Lifestyles

Objectives: We aimed to explore how healthy lifestyles and genetic factors influence the risk of Osteoporosis (OP). Methods: In this prospective cohort study, we first performed a genome-wide association study (GWAS) of estimated bone mineral density (eBMD) and constructed the genetic risk score (GRS) based on the effect of single nucleotide polymorphism (SNP) on eBMD. We then assessed the effect of three-level GRS and adherence to healthy lifestyles on the risk of OP and fracture, respectively. Finally, we assessed the joint effects of GRS and lifestyle on the OP and fracture risk. Results: People with higher GRS have a lower risk of OP and fracture. Negative associations were detected between healthy lifestyle factors and the risk of OP and fracture. Compare with the group with high GRS and favorable lifestyles, the group with low GRS and unfavorable lifestyles had a high Hazard Ratio (HR). Conclusion: The findings suggest that adherence to healthy lifestyles can reduce the risk of OP and fracture in people with different genetic risks.

Lobular Difference in Heritability of Brain Atrophy among Elderly Japanese: A Twin Study

Background and Objectives: Brain atrophy is related to cognitive decline. However, the heritability of brain atrophy has not been fully investigated in the Eastern Asian population. Materials and Methods: Brain imaging of 74 Japanese twins registered in the Osaka University Twin Registry was conducted with voxel-based morphometry SPM12 and was processed by individual voxel-based morphometry adjusting covariates (iVAC) toolbox. The atrophy of the measured lobes was obtained by comparing the focal volume to the average of healthy subjects. Classical twin analysis was used to measure the heritability of its z-scores. Results: The heritability of brain atrophy ranged from 0.23 to 0.97, depending upon the lobes. When adjusted to age, high heritability was reported in the frontal, frontal-temporal, and parietal lobes, but the heritability in other lobes was lower than 0.70. Conclusions: This study revealed a relatively lower heritability in brain atrophy compared to other ethnicities. This result suggests a significant environmental impact on the susceptibility of brain atrophy the Japanese. Therefore, environmental factors may have more influence on the Japanese than in other populations.

Identification of Known and Novel Long Noncoding RNAs Potentially Responsible for the Effects of Bone Mineral Density (BMD) Genomewide Association Study (GWAS) Loci

Osteoporosis, characterized by low bone mineral density (BMD), is the most common complex disease affecting bone and constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 associations influencing BMD. It has been shown that perturbations to long noncoding RNAs (lncRNAs) influence BMD and the activities of bone cells; however, the extent to which lncRNAs are involved in the genetic regulation of BMD is unknown. Here, we combined the analysis of allelic imbalance (AI) in human acetabular bone fragments with a transcriptome-wide association study (TWAS) and expression quantitative trait loci (eQTL) colocalization analysis using data from the Genotype-Tissue Expression (GTEx) project to identify lncRNAs potentially responsible for GWAS associations. We identified 27 lncRNAs in bone that are located in proximity to a BMD GWAS association and harbor single-nucleotide polymorphisms (SNPs) demonstrating AI. Using GTEx data we identified an additional 31 lncRNAs whose expression was associated (false discovery rate [FDR] correction < 0.05) with BMD through TWAS and had a colocalizing eQTL (regional colocalization probability [RCP] > 0.1). The 58 lncRNAs are located in 43 BMD associations. To further support a causal role for the identified lncRNAs, we show that 23 of the 58 lncRNAs are differentially expressed as a function of osteoblast differentiation. Our approach identifies lncRNAs that are potentially responsible for BMD GWAS associations and suggest that lncRNAs play a role in the genetics of osteoporosis. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Increased production of 8-oxo-7,8-dihydroguanine in human urine, a novel biomarker of osteoporosis

Osteoporosis is a worldwide disease that seriously affects the quality of life and survival rate of the elderly. The detection of bone biomarkers will provide supplementary information of bone mineral density, contributing to the accurate diagnosis of osteoporosis and better health care for prevention. This study aimed to investigate the efficacy of oxidative stress markers-8-oxo-7,8-dihydroguanine (8-oxoGsn) and 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodGsn) in the assessment of osteoporosis. We conducted a cross-sectional study among menopausal women with a mean (standard deviation) age of 62.967 (7.798) years old (n = 151). Participants were recruited for the bone mineral density (BMD) assessment, blood and urinary samples. Urinary 8-oxo-7,8-dihydro-2'-deoxyguanosine and 8-oxo-7,8-dihydro-guanine concentrations were measured by ultra performance liquid chromatography and tandem mass spectrometry (UPLC-MS/MS). The urinary 8-oxoGsn/Cre value differed significantly between normal and osteoporotic participants (p < 0.001), while the 8-oxodGsn/Cre value did not (p = 0.720). Even after adjusting for the age and body mass index, the BMD was still associated with urinary 8-oxoGsn/Cre value. ROC analysis showed that 8-oxoGsn has a strong diagnostic value for osteoporosis (AUC =0.744). The results show for the first time that 8-oxoGsn may be a biomarker for the future diagnosis of osteoporosis in women.

The genetic overlap between osteoporosis and craniosynostosis

Osteoporosis is the most prevalent bone condition in the ageing population. This systemic disease is characterized by microarchitectural deterioration of bone, leading to increased fracture risk. In the past 15 years, genome-wide association studies (GWAS), have pinpointed hundreds of loci associated with bone mineral density (BMD), helping elucidate the underlying molecular mechanisms and genetic architecture of fracture risk. However, the challenge remains in pinpointing causative genes driving GWAS signals as a pivotal step to drawing the translational therapeutic roadmap. Recently, a skull BMD-GWAS uncovered an intriguing intersection with craniosynostosis, a congenital anomaly due to premature suture fusion in the skull. Here, we recapitulate the genetic contribution to both osteoporosis and craniosynostosis, describing the biological underpinnings of this overlap and using zebrafish models to leverage the functional investigation of genes associated with skull development and systemic skeletal homeostasis.

Association between alcohol and bone mineral density in a Mendelian randomization study: the Dong-gu study

INTRODUCTION

Many previous studies have reported a positive relationship between alcohol and bone mineral density (BMD). However, the causality between alcohol and BMD has not been fully evaluated.

MATERIALS AND METHODS

This study enrolled 8892 participants from the Dong-gu study. Mendelian randomization (MR) using two-stage least-squared regression was used to evaluate the association between the genetically predicted amount of alcohol consumption per day and BMD. The aldehyde dehydrogenase 2 (ALDH2) rs671 polymorphism was used as instrumental variables for alcohol consumption. Age, smoking history, and BMI were adjusted in the multivariate model.

RESULTS

Self-reported alcohol consumption was positively related to total hip and lumbar spine BMD in both sexes. In multivariate Mendelian randomization analysis, the genetically predicted amount of alcohol consumption was positively associated with both total hip and lumbar spine BMD in men. Total hip BMD and lumbar spine BMD increased by 0.004 g/cm² (95% confidence interval [CI] 0.002-0.007) and 0.007 g/cm² (95% CI 0.004-0.011) with doubling of alcohol consumption. However, in women, genetically predicted alcohol consumption was not significantly associated with BMD.

CONCLUSION

In our MR study, genetically predicted alcohol consumption was positively associated with BMD in men. This result suggests that the association between alcohol consumption and BMD is causal.

Human genotype-to-phenotype predictions: Boosting accuracy with nonlinear models

Genotype-to-phenotype prediction is a central problem of human genetics. In recent years, it has become possible to construct complex predictive models for phenotypes, thanks to the availability of large genome data sets as well as efficient and scalable machine learning tools. In this paper, we make a threefold contribution to this problem. First, we ask if state-of-the-art nonlinear predictive models, such as boosted decision trees, can be more efficient for phenotype prediction than conventional linear models. We find that this is indeed the case if model features include a sufficiently rich set of covariates, but probably not otherwise. Second, we ask if the conventional selection of single nucleotide polymorphisms (SNPs) by genome wide association studies (GWAS) can be replaced by a more efficient procedure, taking into account information in previously selected SNPs. We propose such a procedure, based on a sequential feature importance estimation with decision trees, and show that this approach indeed produced informative SNP sets that are much more compact than when selected with GWAS. Finally, we show that the highest prediction accuracy can ultimately be achieved by ensembling individual linear and nonlinear models. To the best of our knowledge, for some of the phenotypes that we consider (asthma, hypothyroidism), our results are a new state-of-the-art.

Identification of New Genes and Loci Associated With Bone Mineral Density Based on Mendelian Randomization

Bone mineral density (BMD) is a complex and highly hereditary trait that can lead to osteoporotic fractures. It is estimated that BMD is mainly affected by genetic factors (about 85%). BMD has been reported to be associated with both common and rare variants, and numerous loci related to BMD have been identified by genome-wide association studies (GWAS). We systematically integrated expression quantitative trait loci (eQTL) data with GWAS summary statistical data. We mainly focused on the loci, which can affect gene expression, so Summary data-based Mendelian randomization (SMR) analysis was implemented to investigate new genes and loci associated with BMD. We identified 12,477 single-nucleotide polymorphisms (SNPs) regulating 564 genes, which are associated with BMD. The genetic mechanism we detected could make a contribution in the density of BMD in individuals and play an important role in understanding the pathophysiology of cataclasis.

Using "-omics" Data to Inform Genome-wide Association Studies (GWASs) in the Osteoporosis Field

PURPOSE OF REVIEW

Osteoporosis constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 loci influencing bone mineral density (BMD); however, few of the causal genes have been identified. Here, we review approaches that use "-omics" data and genetic- and systems genetics-based analytical strategies to facilitate causal gene discovery.

RECENT FINDINGS

The bone field is beginning to adopt approaches that are commonplace in other disease disciplines. The slower progress has been due in part to the lack of large-scale "omics" data on bone and bone cells. This is however changing, and approaches such as eQTL colocalization, transcriptome-wide association studies (TWASs), network, and integrative approaches are beginning to provide significant insight into the genes responsible for BMD GWAS associations. The use of "-omics" data to inform BMD GWASs has increased in recent years, leading to the identification of novel regulators of BMD in humans. The ultimate goal will be to use this information to develop more effective therapies to treat and ultimately prevent osteoporosis.

Level and change in bone microarchitectural parameters and their relationship with previous fracture and established bone mineral density loci

BACKGROUND

Osteoporosis is characterised by a reduction of bone mineral density (BMD) and predisposition to fracture. Bone microarchitecture, measured by high resolution peripheral quantitative computed tomography (HR-pQCT), has been related to fragility fractures and BMD and has been the subject of large-scale genome-wide analysis. We investigated whether fracture was related to baseline values and longitudinal changes in bone microarchitecture and whether bone microarchitecture was associated with established BMD loci.

METHODS

115 males and 99 females (aged 72-81 at baseline) from the Hertfordshire Cohort Study (HCS) were analysed. Fracture history was determined in 2011-2012 by self-report and vertebral fracture assessment. Participants underwent HR-pQCT scans of the distal radius and tibia in 2011-2012 and 2017. Previous fracture in relation to baseline values and changes in tibial HR-pQCT parameters was examined using sex-adjusted logistic regression with and without adjustment for age, sociodemographic, lifestyle and clinical characteristics; baseline values and changes in parameters associated with previous fracture were then examined in relation to four established BMD loci after adjustment for sex and age.

RESULTS

Previous fracture was related to: higher trabecular area (fully-adjusted odds ratio [95% CI] per SD greater baseline value: 2.18 [1.27,3.73], p = 0.005); lower total volumetric BMD (0.53 [0.34,0.84], p = 0.007), cortical area (0.53 [0.30,0.95], p = 0.032), cortical BMD (0.56 [0.36,0.88], p = 0.011) and cortical thickness (0.45 [0.27,0.77], p = 0.004); and greater declines in trabecular BMD (p = 0.001). Associations were robust in sex- and fully-adjusted analysis. Relationships between BMD loci and these HR-pQCT parameters were weak: rs3801387 (WNT16) was related to decline in trabecular BMD (p = 0.011) but no other associations were significant (p > 0.05).

CONCLUSION

Baseline values of HR-pQCT parameters and greater decline in trabecular BMD were associated with fracture. Change in trabecular BMD was associated with WNT16 which has been demonstrated to influence bone health in murine models and human genome-wide association studies (GWAS).

Strong Genetic Effects on Bone Mineral Density in Multiple Locations with Two Different Techniques: Results from a Cross-Sectional Twin Study

Background and Objectives: Previous studies have demonstrated that risk of hip fracture is at least partly heritable. The aim of this study was to determine the magnitude of the genetic component of bone mineral density (BMD), using both X-ray and ultrasound assessment at multiple sites. Materials and Methods: 216 adult, healthy Hungarian twins (124 monozygotic, MZ, 92 dizygotic, DZ; mean age 54.2 ± 14.3 years), recruited from the Hungarian Twin Registry with no history of oncologic disease underwent cross-sectional BMD studies. We measured BMD, T- and Z-scores with dual energy X-ray absorptiometry (DEXA) at multiple sites (lumbar spine, femoral neck, total hip and radius). Quantitative bone ultrasound (QUS) was also performed, resulting in a calculated value of estimated bone mineral density (eBMD) in the heel bone. Heritability was calculated using the univariate ACE model. Results: Bone density had a strong genetic component at all sites with estimates of heritability ranging from 0.613 to 0.838 in the total sample. Lumbar BMD and calcaneus eBMD had major genetic components with estimates of 0.828 and 0.838 respectively, and least heritable (0.653) at the total hip. BMD of the radius had also a strong genetic component with an estimate of 0.806. No common environmental effect was found. The remaining variance was influenced by unique environment (0.162 to 0.387). In females only, slightly higher additive genetic estimates were found, especially in the case of the femoral neck and total hip. Conclusion: Bone mineral density is strongly heritable, especially in females at all locations using both DEXA and QUS, which may explain the importance of family history as a risk factor for bone fractures. Unshared environmental effects account for the rest of the variance with slight differences in magnitude across various bone regions, supporting the role of lifestyle in preventing osteoporotic fractures with various efficacy in different bone regions.

A Meta-Analysis of the Transferability of Bone Mineral Density Genetic Loci Associations From European to African Ancestry Populations

Genetic studies of bone mineral density (BMD) largely have been conducted in European populations. We therefore conducted a meta-analysis of six independent African ancestry cohorts to determine whether previously reported BMD loci identified in European populations were transferable to African ancestry populations. We included nearly 5000 individuals with both genetic data and assessments of BMD. Genotype imputation was conducted using the 1000G reference panel. We assessed single-nucleotide polymorphism (SNP) associations with femoral neck and lumbar spine BMD in each cohort separately, then combined results in fixed effects (or random effects if study heterogeneity was high, I2 index >60) inverse variance weighted meta-analyses. In secondary analyses, we conducted locus-based analyses of rare variants using SKAT-O. Mean age ranged from 12 to 68 years. One cohort included only men and another cohort included only women; the proportion of women in the other four cohorts ranged from 52% to 63%. Of 56 BMD loci tested, one locus, 6q25 (C6orf97, p = 8.87 × 10-4 ), was associated with lumbar spine BMD and two loci, 7q21 (SLC25A13, p = 2.84 × 10-4 ) and 7q31 (WNT16, p = 2.96 × 10-5 ), were associated with femoral neck BMD. Effects were in the same direction as previously reported in European ancestry studies and met a Bonferroni-adjusted p value threshold, the criteria for transferability to African ancestry populations. We also found associations that met locus-specific Bonferroni-adjusted p value thresholds in 11q13 (LRP5, p < 2.23 × 10-4 ), 11q14 (DCDC5, p < 5.35 × 10-5 ), and 17p13 (SMG6, p < 6.78 × 10-5 ) that were not tagged by European ancestry index SNPs. Rare single-nucleotide variants in AKAP11 (p = 2.32 × 10-2 ), MBL2 (p = 4.09 × 10-2 ), MEPE (p = 3.15 × 10-2 ), SLC25A13 (p = 3.03 × 10-2 ), STARD3NL (p = 3.35 × 10-2 ), and TNFRSF11A (p = 3.18 × 10-3 ) were also associated with BMD. The majority of known BMD loci were not transferable. Larger genetic studies of BMD in African ancestry populations will be needed to overcome limitations in statistical power and to identify both other loci that are transferable across populations and novel population-specific variants. © 2020 American Society for Bone and Mineral Research (ASBMR).

Improved prediction of fracture risk leveraging a genome-wide polygenic risk score

BACKGROUND

Accurately quantifying the risk of osteoporotic fracture is important for directing appropriate clinical interventions. While skeletal measures such as heel quantitative speed of sound (SOS) and dual-energy X-ray absorptiometry bone mineral density are able to predict the risk of osteoporotic fracture, the utility of such measurements is subject to the availability of equipment and human resources. Using data from 341,449 individuals of white British ancestry, we previously developed a genome-wide polygenic risk score (PRS), called gSOS, that captured 25.0% of the total variance in SOS. Here, we test whether gSOS can improve fracture risk prediction.

METHODS

We examined the predictive power of gSOS in five genome-wide genotyped cohorts, including 90,172 individuals of European ancestry and 25,034 individuals of Asian ancestry. We calculated gSOS for each individual and tested for the association between gSOS and incident major osteoporotic fracture and hip fracture. We tested whether adding gSOS to the risk prediction models had added value over models using other commonly used clinical risk factors.

RESULTS

A standard deviation decrease in gSOS was associated with an increased odds of incident major osteoporotic fracture in populations of European ancestry, with odds ratios ranging from 1.35 to 1.46 in four cohorts. It was also associated with a 1.26-fold (95% confidence interval (CI) 1.13-1.41) increased odds of incident major osteoporotic fracture in the Asian population. We demonstrated that gSOS was more predictive of incident major osteoporotic fracture (area under the receiver operating characteristic curve (AUROC) = 0.734; 95% CI 0.727-0.740) and incident hip fracture (AUROC = 0.798; 95% CI 0.791-0.805) than most traditional clinical risk factors, including prior fracture, use of corticosteroids, rheumatoid arthritis, and smoking. We also showed that adding gSOS to the Fracture Risk Assessment Tool (FRAX) could refine the risk prediction with a positive net reclassification index ranging from 0.024 to 0.072.

CONCLUSIONS

We generated and validated a PRS for SOS which was associated with the risk of fracture. This score was more strongly associated with the risk of fracture than many clinical risk factors and provided an improvement in risk prediction. gSOS should be explored as a tool to improve risk stratification to identify individuals at high risk of fracture.

Evaluating and Strengthening the Evidence for Nutritional Bone Research: Ready to Break New Ground?

A healthy diet is essential to attain genetically determined peak bone mass and maintain optimal skeletal health across the adult lifespan. Despite the importance of nutrition for bone health, many of the nutritional requirements of the skeleton across the lifespan remain underexplored, poorly understood, or controversial. With increasingly aging populations, combined with rapidly changing diets and lifestyles globally, one anticipates large increases in the prevalence of osteoporosis and incidence of osteoporotic fractures. Robust, transparent, and reproducible nutrition research is a cornerstone for developing reliable public health recommendations to prevent osteoporosis and osteoporotic fractures. However, nutrition research is often criticized or ignored by healthcare professionals due to the overemphasis of weak science, conflicting, confusing or implausible findings, industry interests, common misconceptions, and strong opinions. Conversely, spurious research findings are often overemphasized or misconstrued by the media or prominent figures especially via social media, potentially leading to confusion and a lack of trust by the general public. Recently, reforms of the broader discipline of nutrition science have been suggested and promoted, leading to new tools and recommendations to attempt to address these issues. In this perspective, we provide a brief overview of what has been achieved in the field on nutrition and bone health, focusing on osteoporosis and osteoporotic fractures. We discuss what we view as some of the challenges, including inherent difficulties in assessing diet and its change, disentangling complex interactions between dietary components and between diet and other factors, selection of bone-related outcomes for nutrition studies, obtaining evidence with more unbiased designs, and perhaps most importantly, ensuring the trust of the public and healthcare professionals. This perspective also provides specific recommendations and highlights new developments and future opportunities for scientists studying nutrition and bone health. © 2021 American Society for Bone and Mineral Research (ASBMR).

LGR4 Gene Polymorphisms Are Associated With Bone and Obesity Phenotypes in Chinese Female Nuclear Families

OBJECTIVE

The current study was conducted to determine whether peak bone mineral density (BMD) and obesity phenotypes are associated with certain LGR4 gene polymorphisms found in Chinese nuclear families with female children.

METHODS

A total of 22 single nucleotide polymorphisms (SNPs) located in and around the LGR4 gene were identified in 1,300 subjects who were members of 390 Chinese nuclear families with female children. Then, BMD readings of the femoral neck, total hip, and lumbar spine as well as measurements of the total lean mass (TLM), total fat mass (TFM), and trunk fat mass were obtained via dual-energy X-ray absorptiometry. The quantitative transmission disequilibrium test was used to analyze the associations between specific SNPs and LGR4 haplotypes and peak BMD as well as between LGR4 haplotypes and TLM, percent lean mass, TFM, percent fat mass, trunk fat mass, and body mass index (BMI).

RESULTS

Here, rs7936621 was significantly associated with the BMD values for the total hip and lumbar spine, while rs10835171 and rs6484295 were associated with the trunk fat mass and BMI, respectively. Regarding the haplotypes, we found significant associations between GAA in block 2 and trunk fat mass and BMI, between AGCGT in block 3 and total hip BMD, between TGCTCC in block 5 and femoral neck BMD, and between TACTTC in block 5 and both lumbar spine and femoral neck BMD (all P-values < 0.05).

CONCLUSION

Genetic variations of the LGR4 gene are related to peak BMD, BMI, and trunk fat mass.

MEPE loss-of-function variant associates with decreased bone mineral density and increased fracture risk

A major challenge in genetic association studies is that most associated variants fall in the non-coding part of the human genome. We searched for variants associated with bone mineral density (BMD) after enriching the discovery cohort for loss-of-function (LoF) mutations by sequencing a subset of the Nord-Trøndelag Health Study, followed by imputation in the remaining sample (N = 19,705), and identified ten known BMD loci. However, one previously unreported variant, LoF mutation in MEPE, p.(Lys70IlefsTer26, minor allele frequency [MAF] = 0.8%), was associated with decreased ultradistal forearm BMD (P-value = 2.1 × 10⁻¹⁸), and increased osteoporosis (P-value = 4.2 × 10⁻⁵) and fracture risk (P-value = 1.6 × 10⁻⁵). The MEPE LoF association with BMD and fractures was further evaluated in 279,435 UK (MAF = 0.05%, heel bone estimated BMD P-value = 1.2 × 10⁻¹⁶, any fracture P-value = 0.05) and 375,984 Icelandic samples (MAF = 0.03%, arm BMD P-value = 0.12, forearm fracture P-value = 0.005). Screening for the MEPE LoF mutations before adulthood could potentially prevent osteoporosis and fractures due to the lifelong effect on BMD observed in the study. A key implication for precision medicine is that high-impact functional variants missing from the publicly available cosmopolitan panels could be clinically more relevant than polygenic risk scores.

Bone mineral density (BMD) is associated with fracture risk and many genetic loci with small effect sizes have been discovered by genome-wide association studies (GWAS). Here, the authors discover a large-effect rare loss-of-function genetic variant for BMD in the MEPE gene in the Norwegian HUNT study which replicates in the UK Biobank.

Comprehensive Genetic Analysis of 128 Candidate Genes in a Cohort With Idiopathic, Severe, or Familial Osteoporosis

Context

The genetic bases of osteoporosis (OP), a disorder with high heritability, are poorly understood at an individual level. Cases of idiopathic or familial OP have long puzzled clinicians as to whether an actionable genetic cause could be identified.

Objective

We performed a genetic analysis of 28 cases of idiopathic, severe, or familial osteoporosis using targeted massively parallel sequencing.

Design

Targeted sequencing of 128 candidate genes was performed using Illumina NextSeq. Variants of interest were confirmed by Sanger sequencing or SNP array.

Patients and Setting

Thirty-seven patients in an academic tertiary hospital participated (54% male; median age, 44 years; 86% with fractures), corresponding to 28 sporadic or familial cases.

Main Outcome Measure

The identification of rare stop-gain, indel, splice site, copy-number, or nonsynonymous variants altering protein function.

Results

Altogether, we identified 28 variants of interest, but only 3 were classified as pathogenic or likely pathogenic variants: COL1A2 p.(Arg708Gln), WNT1 p.(Gly169Asp), and IDUA p.(His82Gln). An association of variants in different genes was found in 21% of cases, including a young woman with severe OP bearing WNT1, PLS3, and NOTCH2 variants. Among genes of uncertain significance analyzed, a potential additional line of evidence has arisen for GWAS candidates GPR68 and NBR1, warranting further studies.

Conclusions

While we hope that continuing efforts to identify genetic predisposition to OP will lead to improved and personalized care in the future, the likelihood of identifying actionable pathogenic variants in intriguing cases of idiopathic or familial osteoporosis is seemingly low.

The Musculoskeletal Knowledge Portal: Making Omics Data Useful to the Broader Scientific Community

The development of high-throughput genotyping technologies and large biobank collections, complemented with rapid methodological advances in statistical genetics, has enabled hypothesis-free genome-wide association studies (GWAS), which have identified hundreds of genetic variants across many loci associated with musculoskeletal conditions. Similarly, basic scientists have valuable molecular cellular and animal data based on musculoskeletal disease that would be enhanced by being able to determine the human translation of their findings. By integrating these large-scale human genomic musculoskeletal datasets with complementary evidence from model organisms, new and existing genetic loci can be statistically fine-mapped to plausibly causal variants, candidate genes, and biological pathways. Genes and pathways identified using this approach can be further prioritized as drug targets, including side-effect profiling and the potential for new indications. To bring together these big data, and to realize the vision of creating a knowledge portal, the International Federation of Musculoskeletal Research Societies (IFMRS) established a working group to collaborate with scientists from the Broad Institute to create the Musculoskeletal Knowledge Portal (MSK-KP)(http://mskkp.org/). The MSK consolidates omics datasets from humans, cellular experiments, and model organisms into a central repository that can be accessed by researchers. The vision of the MSK-KP is to enable better understanding of the biological mechanisms underlying musculoskeletal disease and apply this knowledge to identify and develop new disease interventions. © 2020 American Society for Bone and Mineral Research (ASBMR).

A single genetic locus associated with pediatric fractures: A genome-wide association study on 3,230 patients

The understanding of the biological and environmental risk factors of fractures in pediatrics is limited. Previous studies have reported that fractures involve heritable traits, but the genetic factors contributing to the risk of fractures remain elusive. Furthermore, genetic influences specific to immature bone have not been thoroughly studied. Therefore, the aim of the present study was to identify genetic variations that are associated with fractures in early childhood. The present study used a prospective Northern Finland Birth Cohort (year 1986; n=9,432). The study population was comprised of 3,230 cohort members with available genotype data. A total of 48 members of the cohort (1.5%) had in-hospital treated bone fractures during their first 6 years of life. Furthermore, individuals without fracture (n=3,182) were used as controls. A genome-wide association study (GWAS) was performed using a frequentist association test. In the GWAS analysis, a linear regression model was fitted to test for additive effects of single-nucleotide polymorphisms (SNPs; genotype dosage) adjusting for sex and performing population stratification using genotypic principal components. Using the GWAS analysis, the present study identified one locus with a significant association with fractures during childhood on chromosome 10 (rs112635931) and six loci with a suggested implication. The lead SNP rs112635931 was located near proline- and serine-rich 2 (PROSER2) antisense RNA 1 (PROSER2-AS1) and PROSER2, thus suggesting that these may be novel candidate genes associated with the risk of pediatric fractures.

A genome-wide scan for pleiotropy between bone mineral density and nonbone phenotypes

Osteoporosis is the most common metabolic bone disorder globally and is characterized by skeletal fragility and microarchitectural deterioration. Genetic pleiotropy occurs when a single genetic element is associated with more than one phenotype. We aimed to identify pleiotropic loci associated with bone mineral density (BMD) and nonbone phenotypes in genome-wide association studies. In the discovery stage, the NHGRI-EBI Catalog was searched for genome-wide significant associations (P value < 5 × 10-8), excluding bone-related phenotypes. SNiPA was used to identify proxies of the significantly associated single nucleotide polymorphisms (SNPs) (r 2 = 1). We then assessed putative genetic associations of this set of SNPs with femoral neck (FN) and lumbar spine (LS) BMD data from the GEFOS Consortium. Pleiotropic variants were claimed at a false discovery rate < 1.4 × 10-3 for FN-BMD and < 1.5 × 10-3 for LS-BMD. Replication of these genetic markers was performed among more than 400 000 UK Biobank participants of European ancestry with available genetic and heel bone ultrasound data. In the discovery stage, 72 BMD-related pleiotropic SNPs were identified, and 12 SNPs located in 11 loci on 8 chromosomes were replicated in the UK Biobank. These SNPs were associated, in addition to BMD, with 14 different phenotypes. Most pleiotropic associations were exhibited by rs479844 (AP5B1, OVOL1 genes), which was associated with dermatological and allergic diseases, and rs4072037 (MUC1 gene), which was associated with magnesium levels and gastroenterological cancer. In conclusion, 12 BMD-related genome-wide significant SNPs showed pleiotropy with nonbone phenotypes. Pleiotropic associations can deepen the genetic understanding of bone-related diseases by identifying shared biological mechanisms with other diseases or traits.

Targeted next generation sequencing of nine osteoporosis-related genes in the Wnt signaling pathway among Chinese postmenopausal women

PURPOSE

This study aimed to explore the association between low-frequency and rare variants of Wnt signaling genes and postmenopausal osteoporosis (OP) by the next generation sequencing (NGS) technology.

METHODS

We performed targeted NGS of nine Wnt signaling genes in 400 Chinese postmenopausal women, including 226 cases with decreased bone mineral density (BMD) and 174 controls with normal values. Proxy External Controls Association Test (ProxECAT) and logistic regression analysis were performed by data from internal cases (n = 226) and Genome Aggregation Database (gnomAD) East Asian controls (n = 9435).

RESULTS

The genomic region of interest (ROI) of 94 functional low-frequency and rare variants was associated with OP risk (P < 0.05). The LGR6 gene was associated significantly with OP risk and BMD measurements (BMD, T-score and Z-score) (adjusted-P < 0.05) after adjusting for confounders. The allele A of rs199693693 (K82N) in LRP6 and G of novel variant 1: 202287949 (R840G) in LGR6 were associated with higher BMD, T-score, and Z-score (all adjusted-P < 0.05). ProxECAT showed that LGR4 was significantly different between the internal cases and the external controls (all adjusted-P < 0.05). Logistic regression analysis revealed that the allele G of rs765778410 (T645A) (OR = 26.16, 95% CI: 4.36-156.95, adjusted-P value = 0.026) in LGR6 and A of rs61370283 (L987M) (OR = 15.39, 95% CI: 2.98-79.55, adjusted-P value = 0.037) in LRP5 were associated with increased risk of postmenopausal OP.

CONCLUSION

The LGR4 and LGR6 genes and four potential functional rare variants associate with postmenopausal OP risk. These results highlight the significance of rare functional variants in postmenopausal OP genetics and provide new insights into the potential mutations in this field.

Identification of PIEZO1 polymorphisms for human bone mineral density

Bone mineral density (BMD) is a key indicator for diagnosis and treatment for osteoporosis; the reduction of BMD could increase the risk of osteoporotic fracture. It was very recently found that Piezo1 mediated mechanically evoked responses in bone and further participated in bone formation in mice. Here, we performed cross phenotype meta-analysis for human BMD at lumbar spine (LS), femoral neck (FN), distal radius/forearm (FA) and heel and screened out 14 top SNPs for PIEZO1, these SNPs were overlapped with putative enhancers, DNase-I hypersensitive sites and active promoter flanking regions. We found that the signal of the best SNP rs62048221 was mainly from heel ultrasound estimated BMD (-0.02 SD per T allele, P = 8.50E-09), where calcaneus supported most of the mechanical force of body when standing, walking and doing physical exercises. Each copy of the effect allele T of SNP rs62048221 was associated with a decrease of 0.0035 g/cm² BMD (P = 4.6E-27, SE = 0.0003) in UK Biobank data within 477,760 samples. SNP rs62048221 was located at the enhancer region (HEDD enhancer ID 2331049) of gene PIEZO1, site-directed ChIP assays in human mesenchymal stem cells (hMSCs) showed significant enrichment of H3K4me1 and H3K27ac in this region, luciferase assays showed that rs62048221 could significantly affect the activity of the enhancer where it resides. Our results first suggested that SNP rs62048221 might mediate the PIEZO1 expression level via modulating the activity of cis-regulatory elements and then further affect the BMD.

BMD-Related Genetic Risk Scores Predict Site-Specific Fractures as Well as Trabecular and Cortical Bone Microstructure

CONTEXT

It is important to identify patients at highest risk of fractures.

OBJECTIVE

To compare the separate and combined performances of bone-related genetic risk scores (GRSs) for prediction of forearm, hip and vertebral fractures separately, as well as of trabecular and cortical bone microstructure parameters separately.

DESIGN, SETTING, AND PARTICIPANTS

Using 1103 single nucleotide polymorphisms (SNPs) independently associated with estimated bone mineral density of the heel (eBMD), we developed a weighted GRS for eBMD and determined its contribution to fracture prediction beyond 2 previously developed GRSs for femur neck BMD (49 SNPs) and lumbar spine BMD (48 SNPs). Associations between these GRSs and forearm (ncases = 1020; ncontrols = 2838), hip (ncases = 1123; ncontrols = 2630) and vertebral (ncases = 288; ncontrols = 1187) fractures were evaluated in 3 Swedish cohorts. Associations between the GRSs and trabecular and cortical bone microstructure parameters (n = 426) were evaluated in the MrOS Sweden cohort.

RESULTS

We found that eBMDGRS was the only significant independent predictor of forearm and vertebral fractures while both FN-BMDGRS and eBMDGRS were significant independent predictors of hip fractures. The eBMDGRS was the major GRS contributing to prediction of trabecular bone microstructure parameters while both FN-BMDGRS and eBMDGRS contributed information for prediction of cortical bone microstructure parameters.

CONCLUSIONS

The eBMDGRS independently predicts forearm and vertebral fractures while both FN-BMDGRS and eBMDGRS contribute independent information for prediction of hip fractures. We propose that eBMDGRS captures unique information about trabecular bone microstructure useful for prediction of forearm and vertebral fractures. These findings may facilitate personalized medicine to predict site-specific fractures as well as cortical and trabecular bone microstructure separately.

Identification of novel variants and candidate genes associated with porcine bone mineral density using genome-wide association study

Pig leg weakness not only causes huge economic losses for producers but also affects animal welfare. However, genes with large effects on pig leg weakness have not been identified and suitable methods to study porcine leg weakness are urgently needed. Bone mineral density (BMD) is an important indicator for determining leg soundness in pigs. Increasing pig BMD is likely to improve pig leg soundness. In this study, porcine BMD was measured using an ultrasound bone densitometer in a population with 212 Danish Landrace pigs and 537 Danish Yorkshires. After genotyping all the individuals using GeneSeek Porcine 50K SNP chip, genetic parameter estimation was performed to evaluate the heritability of BMD. Genome-wide association study and haplotype analysis were also performed to identify the variants and candidate genes associated with porcine BMD. The results showed that the heritability of BMD was 0.21 in Landrace and 0.31 in Yorkshire. Five single-nucleotide polymorphisms on chromosome 6 identified were associated with porcine BMD at suggestive significance level. Two candidate quantitative trait loci (74.47 to 75.33 Mb; 80.20 to 83.83 Mb) and three potential candidate genes (ZBTB40, CNR2, and Lin28a) of porcine BMD were detected in this study.

Characterisation of genetic regulatory effects for osteoporosis risk variants in human osteoclasts

Background

Osteoporosis is a complex disease with a strong genetic contribution. A recently published genome-wide association study (GWAS) for estimated bone mineral density (eBMD) identified 1103 independent genome-wide significant association signals. Most of these variants are non-coding, suggesting that regulatory effects may drive many of the associations. To identify genes with a role in osteoporosis, we integrate the eBMD GWAS association results with those from our previous osteoclast expression quantitative trait locus (eQTL) dataset.

Results

We identify sixty-nine significant cis-eQTL effects for eBMD GWAS variants after correction for multiple testing. We detect co-localisation of eBMD GWAS and osteoclast eQTL association signals for 21 of the 69 loci, implicating a number of genes including CCR5, ZBTB38, CPE, GNA12, RIPK3, IQGAP1 and FLCN. Summary-data-based Mendelian Randomisation analysis of the eBMD GWAS and osteoclast eQTL datasets identifies significant associations for 53 genes, with TULP4 presenting as a strong candidate for pleiotropic effects on eBMD and gene expression in osteoclasts. By performing analysis using the GARFIELD software, we demonstrate significant enrichment of osteoporosis risk variants among high-confidence osteoclast eQTL across multiple GWAS P value thresholds. Mice lacking one of the genes of interest, the apoptosis/necroptosis gene RIPK3, show disturbed bone micro-architecture and increased osteoclast number, highlighting a new biological pathway relevant to osteoporosis.

Conclusion

We utilise a unique osteoclast eQTL dataset to identify a number of potential effector genes for osteoporosis risk variants, which will help focus functional studies in this area.

Mendelian Randomization Identifies CpG Methylation Sites With Mediation Effects for Genetic Influences on BMD in Peripheral Blood Monocytes

Osteoporosis is mainly characterized by low bone mineral density (BMD) and is an increasingly serious public health concern. DNA methylation is a major epigenetic mechanism that may contribute to the variation in BMD and may mediate the effects of genetic and environmental factors of osteoporosis. In this study, we performed an epigenome-wide DNA methylation analysis in peripheral blood monocytes of 118 Caucasian women with extreme BMD values. Further, we developed and implemented a novel analytical framework that integrates Mendelian randomization with genetic fine mapping and colocalization to evaluate the causal relationships between DNA methylation and BMD phenotype. We identified 2,188 differentially methylated CpGs (DMCs) between the low and high BMD groups and distinguished 30 DMCs that may mediate the genetic effects on BMD. The causal relationship was further confirmed by eliminating the possibility of horizontal pleiotropy, linkage effect and reverse causality. The fine-mapping analysis determined 25 causal variants that are most likely to affect the methylation levels at these mediator DMCs. The majority of the causal methylation quantitative loci and DMCs reside within cell type-specific histone mark peaks, enhancers, promoters, promoter flanking regions and CTCF binding sites, supporting the regulatory potentials of these loci. The established causal pathways from genetic variant to BMD phenotype mediated by DNA methylation provide a gene list to aid in designing future functional studies and lead to a better understanding of the genetic and epigenetic mechanisms underlying the variation of BMD.

Use of Mendelian Randomization to Examine Causal Inference in Osteoporosis

Epidemiological studies have identified many risk factors for osteoporosis, however it is unclear whether these observational associations reflect true causal effects, or the effects of latent confounding or reverse causality. Mendelian randomization (MR) enables causal relationships to be evaluated, by examining the relationship between genetic susceptibility to the risk factor in question, and the disease outcome of interest. This has been facilitated by the development of two-sample MR analysis, where the exposure and outcome are measured in different studies, and by exploiting summary result statistics from large well-powered genome-wide association studies that are available for thousands of traits. Though MR has several inherent limitations, the field is rapidly evolving and at least 14 methodological extensions have been developed to overcome these. The present paper aims to discuss some of the limitations in the MR analytical framework, and how this method has been applied to the osteoporosis field, helping to reinforce conclusions about causality, and discovering potential new regulatory pathways, exemplified by our recent MR study of sclerostin.

The genetic architecture of osteoporosis and fracture risk

Osteoporosis and fracture risk are common complex diseases, caused by an interaction of numerous disease susceptibility genes and environmental factors. With the advances in genomic technologies, large-scale genome-wide association studies (GWAS) have been performed which have broadened our understanding of the genetic architecture and biological mechanisms of complex disease. Currently, more than ~90 loci have been found associated with DXA derived bone mineral density (BMD), over ~500 loci with heel estimated BMD and several others with other less widely available bone parameters such as bone geometry, shape, and microarchitecture. Notably, several of the pathways identified by the GWAS efforts correspond to pathways that are currently targeted for the treatment of osteoporosis. Overall, tremendous progress in the field of the genetics of osteoporosis has been achieved with the discovery of WNT16, EN1, DAAM2, and GPC6 among others. Assessment of the function and biological mechanisms of the remaining genes may further untangle the complex genetic landscape of osteoporosis and fracture risk. With this review we aimed to provide a general overview of the existing GWAS studies on osteoporosis traits and fracture risk.

COL1A1, CCDC170, and ESR1 single nucleotide polymorphisms associated with distal radius fracture in postmenopausal Mexican women

Objective: The aim of this study was to analyze the genetic association of five ESR1 single nucleotide polymorphisms (SNPs) (rs3020331, rs851982, rs1999805, rs2234693, rs3020404), four COL1A1 SNPs (rs1800012, rs2075555, rs2412298, rs1107946), and two SNPs on the CCDC170 gene (rs9479055, rs4870044) with distal radius fracture (DRF) in a group of postmenopausal Mexican women.Methods: A case-control study was conducted. Cases (n = 182) were women above the age of 38 years with low-energy DRF, and controls (n = 201) were women without. Analysis was done through real-time polymerase chain reaction. Frequencies and Hardy-Weinberg equilibrium were calculated. A multivariate analysis including bone mass index, age, menarche, and menopause as covariables was carried out. Finally, haplotype and linkage disequilibrium (LD) analyses were performed.Results: COL1A1 rs1107946 was strongly associated with DRF. Both CCDC170 SNPs showed strong association with DRF. For the ESR1 gene, four SNPs (rs2234693, 3020404, rs3020331, and rs851982) showed very strong association with DRF. Additionally, the region between the latter two showed strong LD.Conclusions: A strong association of DRF with variants in these genes was found, including haplotypes and a region with strong LD on ESR1. The results suggest that these SNPs could be useful to detect the population at risk of presenting DRF among Mexican perimenopausal women.

An atlas of genetic influences on osteoporosis in humans and mice

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.

Integrative genomic analysis for the functional roles of ITPKC in bone mineral density

Osteoporosis is defined by low bone mineral density (BMD), which is mainly due to the imbalances in osteoclast and osteoblast activity. Previous studies indicated that early activation of osteoclasts relies on calcium entry through store-operated calcium (SOC) entry, and several genes, including STIM1, ORAI1, and ITPKC, are known as key regulators of SOC entry. However, the relationships between STIM1, ORAI1, ITPKC, and human BMD are still unclear. In order to investigate the plausible associations between these genes and BMD, we conducted a meta-analysis of genes expression and BMD using the publicly available GEO database. We further recruited 1044 subjects and tested associations between polymorphisms in these genes and BMD. Clinical information (including age, sex, and BMI) was collected and used for the analysis. Our results indicated that ITPKC gene expression was significantly associated with BMD. Furthermore, we found that one ITPKC SNP (rs2607420) was significantly associated with lumbar spine BMD. Through bioinformatics analysis, rs2607420 was found to be very likely to participate in the regulation of ITPKC expression. Our findings suggest that ITPKC is a susceptibility gene for BMD, and rs2607420 may play an important role in the regulation of this gene.

Genome-wide identification of m6A-associated SNPs as potential functional variants for bone mineral density

This study investigated the effect of the N⁶-methyladenosine (m⁶A)-associated SNPs on bone mineral density (BMD) and found plenty of m⁶A-SNPs that were associated with BMD. This study increases our understanding on the regulation patterns of SNP and may provide new clues for further detection of functional mechanism underlying the associations between SNPs and osteoporosis.

INTRODUCTION

m⁶A plays critical roles in many fundamental biological processes and a variety of diseases. The m⁶A-associated SNPs may be potential functional variants for BMD. The aim of this study was to investigate the effect of the genome-wide m⁶A-SNPs on BMD.

METHODS

We examined the association of m⁶A-SNPs with femoral neck (FN) and lumbar spine (LS) BMD in 32,961 individuals and quantitative heel ultrasounds (eBMD) in 142,487 individuals. Furthermore, we performed expression quantitative trait locus (eQTL) analyses for the m⁶A-SNPs using whole genome data of about 10.5 million SNPs and 21,323 mRNAs from 43 Chinese individuals, as well as public available data. Differential expression analyses were also performed to support the identified genes.

RESULTS

We found 138, 125, and 993 m⁶A-SNPs which were associated with FN-BMD, LS-BMD, and eBMD (P < 0.05), respectively. The associations of rs11614913 (P = 8.92 × 10^-10) in MIR196A2 and rs1110720 (P = 2.05 × 10^-10) in ESPL1 with LS-BMD reached the genome-wide significance level. In addition, a total of 24 m⁶A-SNPs were significantly associated with eBMD (P < 5.0 × 10^-8). Further eQTL analyses showed that 47 of these BMD-associated m⁶A-SNPs were associated with expressions of the 46 corresponding local genes. Moreover, the expressions of 26 of these genes were associated with BMD.

CONCLUSION

The present study represents the first effort of investigating the associations and the mechanisms underlying the link between m⁶A-SNPs and BMD. The results suggested that m⁶A-SNP may play important roles in the pathology of osteoporosis.

Nonlinear association between bone mineral density and all-cause mortality: the Dong-gu study

There was a U-shaped association between hip BMD and all-cause mortality, with the lowest mortality in the 90th percentile in males. However, there was an inverse linear relationship in females. In contrast, the association between lumbar spine BMD and mortality was less evident in males, with no association in females.

INTRODUCTION

Bone mineral density (BMD) is reported inversely associated with mortality. Although some previous studies provided evidence for nonlinear associations, these were not adequately assessed in most previous works.

METHODS

We evaluated the nonlinear relationship between BMD and mortality in Asians. Our study involved 8629 participants in the Dong-gu study from 2007 to 2010. Cox proportional hazard regression was used to calculate hazard ratios (HRs) according to BMD categories after adjusting for potential confounders. During a follow-up of 6.7 ± 1.4 years, 712 participants died.

RESULTS

There was a U-shaped association between hip BMD and all-cause mortality, with the lowest mortality in the 90th percentile in males. However, there was an inverse linear relationship in females. In males, compared with the 75th to 95th percentile group, the < 2.5th percentile group had a 3.89 (95% CI 2.41-6.28)-fold higher risk and the 2.5th to 5th percentile group had a 2.51 (95% CI 1.25-5.04)-fold higher risk. The HR was 2.51 (95% CI 1.25, 5.04) in the > 97.5th percentile group. In females, compared with that in the 75th to 95th percentile group, the HR was 2.33 (95% CI 1.24, 4.39) in the < 2.5th percentile group. In contrast, the association between lumbar spine BMD and mortality was less evident in males, with no association in females.

CONCLUSION

In conclusion, this study shows that the association between BMD and mortality varies by gender and that high and low BMD are predictors of all-cause mortality in males.

Genetic profiling of decreased bone mineral density in an independent sample of Caucasian women

Genetic risk of low bone mineral density in women remains unclear. This study found that a large percentage of Caucasian women have a high genetic risk of osteoporosis, and genetic risk scores are significantly associated with BMD variation in a bone healthy sample of Caucasian women.

INTRODUCTION

We aimed to examine the distribution of risk alleles in an independent sample and to determine if such genetic components are associated with bone mineral density (BMD) variation in the sample.

METHODS

Existing genotype data of 1205 women in the cross-sectional Genomic Wide Scans for Female Osteoporosis Gene Study (GWSFO) were analyzed. Multi-loci genetic risk scores (GRSs) based on 62 BMD-associated single nucleotide polymorphisms (SNPs) were calculated. Regression analysis was employed to assess the association between GRSs and BMD. To examine the effect of SNPs clustered within key pathways associated with the development of osteoporosis, subtype weighted GRS specific to WNT signaling (6 SNPs), RANK-RANKL-OPG (3 SNPs), and mesenchymal stem differentiation (3 SNPs) were generated for analysis.

RESULTS

The unweighted GRS ranged from 48 to 80. One third of the women carried 66% risk alleles. After adjusting for age, height, and body weight, each unit increase of weighted GRS was associated with a decrease in BMD of 0.097 at femur (p < 0.0001) and 0.110 (p < 0.0001) at lumbar spine. The weighted GRS accounted for only 3.17-4.52% of BMD variance. The WNT signaling pathway GRS (6 SNPs) and the RANK-RANKL-OPG signaling pathway GRS (3 SNPs) both were significantly associated with decreased BMD at femur neck (p = 0.0004 and p = 0.0063, respectively) and lumbar spine (p < 0.0001 and p = 0.0001, respectively), while the mesenchymal stem cell differentiation pathway (3 SNPs) GRSs were associated only with the lumbar spine BMD (p = 0.045).

CONCLUSIONS

A substantially large percentage of healthy Caucasian women have a high genetic risk of osteoporosis. Weighted GRS was significantly associated with decreased BMD. The contribution of subtype GRS to the BMD variation differs by specific biological pathway and skeletal regions.

Identification of 613 new loci associated with heel bone mineral density and a polygenic risk score for bone mineral density, osteoporosis and fracture

Low bone mineral density (BMD) leads to osteoporosis, and is a risk factor for bone fractures, including stress fractures. Using data from UK Biobank, a genome-wide association study identified 1,362 independent SNPs that clustered into 899 loci of which 613 are new. These data were used to train a genetic algorithm using 22,886 SNPs as predictors and showing a correlation with heel bone mineral density of 0.415. Combining this genetic algorithm with height, weight, age and sex resulted in a correlation with heel bone mineral density of 0.496. Individuals with low scores (2.2% of total) showed a change in BMD of -1.16 T-score units, an increase in risk for osteoporosis of 17.4 fold and an increase in risk for fracture of 1.87 fold. Genetic predictors could assist in the identification of individuals at risk for osteoporosis or fractures.

Genetic architecture: the shape of the genetic contribution to human traits and disease

Genetic architecture describes the characteristics of genetic variation that are responsible for heritable phenotypic variability. It depends on the number of genetic variants affecting a trait, their frequencies in the population, the magnitude of their effects and their interactions with each other and the environment. Defining the genetic architecture of a complex trait or disease is central to the scientific and clinical goals of human genetics, which are to understand disease aetiology and aid in disease screening, diagnosis, prognosis and therapy. Recent technological advances have enabled genome-wide association studies and emerging next-generation sequencing studies to begin to decipher the nature of the heritable contribution to traits and disease. Here, we describe the types of genetic architecture that have been observed, how architecture can be measured and why an improved understanding of genetic architecture is central to future advances in the field.

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.