Search

At present, there have only been a few DNA sequencing-based studies to explore the genetic determinants of bone mineral density (BMD). We carried out the largest whole genome sequencing analysis to date for femoral neck and spine BMD (n = 4981), with one of the highest average sequencing depths implemented thus far at 22×, in a multiethnic sample (58% Caucasian and 42% African American) from the Louisiana Osteoporosis Study (LOS). The LOS samples were combined with summary statistics from the GEFOS consortium and several independent samples of various ethnicities to perform GWAS meta-analysis (n = 44 506). We identified 31 and 30 genomic risk loci for femoral neck and spine BMD, respectively. The findings substantiate many previously reported susceptibility loci (e.g. WNT16 and ESR1) and reveal several others that are either novel or have not been widely replicated in GWAS for BMD, including two for femoral neck (IGF2 and ZNF423) and one for spine (SIPA1). Although we were not able to uncover ethnicity specific differences in the genetic determinants of BMD, we did identify several loci which demonstrated sex-specific associations, including two for women (PDE4D and PIGN) and three for men (TRAF3IP2, NFIB and LYSMD4). Gene-based rare variant association testing detected MAML2, a regulator of the Notch signaling pathway, which has not previously been suggested, for association with spine BMD. The findings provide novel insights into the pathophysiological mechanisms of osteoporosis., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Transcriptome-wide association studies (TWAS) systematically investigate the association of genetically predicted gene expression with disease risk, providing an effective approach to identify novel susceptibility genes. Osteoporosis is the most common metabolic bone disease, associated with reduced bone mineral density (BMD) and increased risk of osteoporotic fractures, whereas genetic factors explain approximately 70% of the variance in phenotypes associated with bone. BMD is commonly assessed using dual-energy X-ray absorptiometry (DXA) to obtain measurements (g/cm 2 ) of areal BMD. However, quantitative computed tomography (QCT) measured 3D volumetric BMD (vBMD) (g/cm 3 ) has important advantages compared with DXA since it can evaluate cortical and trabecular microstructural features of bone quality, which can be used to directly predict fracture risk. Here, we performed the first TWAS for volumetric BMD (vBMD) by integrating genome-wide association studies (GWAS) data from two independent cohorts, namely the Framingham Heart Study (FHS, n = 3298) and the Osteoporotic Fractures in Men (MrOS, n = 4641), with tissue-specific gene expression data from the Genotype-Tissue Expression (GTEx) project. We first used stratified linkage disequilibrium (LD) score regression approach to identify 12 vBMD-relevant tissues, for which vBMD heritability is enriched in tissue-specific genes of the given tissue. Focusing on these tissues, we subsequently leveraged GTEx expression reference panels to predict tissue-specific gene expression levels based on the genotype data from FHS and MrOS. The associations between predicted gene expression levels and vBMD variation were then tested by MultiXcan, an innovative TWAS method that integrates information available across multiple tissues. We identified 70 significant genes associated with vBMD, including some previously identified osteoporosis-related genes such as LYRM2 and NME8, as well as some novel loci such as DNAAF2 and SPAG16. Our findings provide novel insights into the pathophysiological mechanisms of osteoporosis and highlight several novel vBMD-associated genes that warrant further investigation., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Objective: To examine the contribution of skeletal mass index (SMI) as a mediator in the relationship between menarcheal age and hip/spine bone mineral density (BMD) in premenopausal women by race/ethnicity., Methods: The data of 4,329 participants (age ≥ 18; mean age=35.7 ± 9.5) of Whites (n = 2,543), African Americans (n = 1,236), and Asians (n = 550) enrolled from October 2011 to January 2019 from the Louisiana Osteoporosis Study were analyzed. After adjustment for physiological and behavioral factors, multivariable linear regression analyses were conducted to evaluate each component of the proposed mediation models, and mediation was verified by the bootstrapping resampling approach., Results: Premenopausal women with early menarcheal age tended to have higher SMI and BMD than women with normal menarcheal age among all races/ethnicities included. Women with late menarcheal age were, however, more likely to have a lower SMI than their counterparts with normal menarcheal age (r = -0.212, 95% CI = [-0.321 to -0.103] for White women; r = -0.181, 95% CI = [-0.410 to -0.008] for African-American women; r = -0.174, 95% CI = [-0.343 to -0.006] for Asian women). Similar results were found for both spine and hip BMD. SMI fully mediated the difference in BMD due to different menarcheal ages among Whites, African Americans, and Asian women with early menarcheal age; however, no mediating effects were observed for Asian women with late menarcheal age., Conclusions: SMI, as a full mediator, affected the relationship between menarcheal age and BMD among premenopausal women, and the mediating effects varied by race/ethnicity. To prevent or slow down the loss of hip/spine BMD and the development of osteoporosis, measures aiming at minimizing the risk for muscle mass loss should be recommended, especially for White and African-American women with late menarcheal age., Competing Interests: Financial disclosure/conflicts of interest: None reported., (Copyright © 2021 by The North American Menopause Society.)

Context: Although metabolic profiles appear to play an important role in menopausal bone loss, the functional mechanisms by which metabolites influence bone mineral density (BMD) during menopause are largely unknown., Objective: We aimed to systematically identify metabolites associated with BMD variation and their potential functional mechanisms in peri- and postmenopausal women., Design and Methods: We performed serum metabolomic profiling and whole-genome sequencing for 517 perimenopausal (16%) and early postmenopausal (84%) women aged 41 to 64 years in this cross-sectional study. Partial least squares regression and general linear regression analysis were applied to identify BMD-associated metabolites, and weighted gene co-expression network analysis was performed to construct co-functional metabolite modules. Furthermore, we performed Mendelian randomization analysis to identify causal relationships between BMD-associated metabolites and BMD variation. Finally, we explored the effects of a novel prominent BMD-associated metabolite on bone metabolism through both in vivo/in vitro experiments., Results: Twenty metabolites and a co-functional metabolite module (consisting of fatty acids) were significantly associated with BMD variation. We found dodecanoic acid (DA), within the identified module causally decreased total hip BMD. Subsequently, the in vivo experiments might support that dietary supplementation with DA could promote bone loss, as well as increase the osteoblast and osteoclast numbers in normal/ovariectomized mice. Dodecanoic acid treatment differentially promoted osteoblast and osteoclast differentiation, especially for osteoclast differentiation at higher concentrations in vitro (eg,10, 100 μM)., Conclusions: This study sheds light on metabolomic profiles associated with postmenopausal osteoporosis risk, highlighting the potential importance of fatty acids, as exemplified by DA, in regulating BMD., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Endocrine Society. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Bone mineral density (BMD) and lean body mass (LBM) not only have a considerable heritability each, but also are genetically correlated. However, common genetic determinants shared by both traits are largely unknown. In the present study, we performed a bivariate genome-wide association study (GWAS) meta-analysis of hip BMD and trunk lean mass (TLM) in 11,335 subjects from 6 samples, and performed replication in estimated heel BMD and TLM in 215,234 UK Biobank (UKB) participants. We identified 2 loci that nearly attained the genome-wide significance (GWS, p < 5.0 × 10 -8 ) level in the discovery GWAS meta-analysis and that were successfully replicated in the UKB sample: 11p15.2 (lead SNP rs12800228, discovery p = 2.88 × 10 -7 , replication p = 1.95 × 10 -4 ) and 18q21.32 (rs489693, discovery p = 1.67 × 10 -7 , replication p = 1.17 × 10 -3 ). The above 2 pleiotropic loci may play a pleiotropic role for hip BMD and TLM development. So our findings provide useful insights that further enhance our understanding of genetic interplay between BMD and LBM.

Aiming to investigate whether genetic risk factors (GRFs) for fracture and bone mineral density (BMD) identified from people of European descent can help improve the prediction of osteoporotic fracture (OF) risk and BMD in Chinese populations, we built assessment models for femoral neck (FN)-fracture prediction and BMD value prediction using 700 elderly Chinese Han subjects and 1,620 unrelated Chinese Han subjects, respectively. 17 fracture-associated genes and 82 FN-BMD associated genes identified in people of European descent were used to build a logistic regression model with clinical risk factors (CRFs) for FN-fracture prediction in Chinese. Meanwhile 107 BMD-associated genes from people of European descent were used to build a multiple linear regression model with CRFs for BMD prediction in Chinese. A Lasso algorithm was employed for informative SNP selection to construct the genetic risk score (GRS) with ten-fold cross-validation. The results showed that, adding fracture GRF and FN-BMD GRF to the model with CRFs, the area under the receiver operating characteristic curve (AUC) decrease from 0.653 to 0.587 and 0.588, respectively, for FN fracture prediction. 62.3% and 61.8% of the risk variation were explained by the Model with CRFs and fracture GRF and by the Model with CRFs and FN-BMD GRF, respectively, as compared to 65.5% in the Model with CRFs only. The net reclassification improvement (NRI) index in the reclassification analysis is 0.56% (P = 0.57) and 1.13% (P = 0.29), respectively. There is no significant difference either between the performance of the model with CRFs and that of the model with both CRFs and GRF for BMD prediction. We concluded that, in the current study, GRF of fracture identified in people of European descent does not contributes to improve the fracture prediction in Chinese; and GRF of BMD from people of European descent cannot help improve the accuracy of the fracture prediction in Chinese perhaps partially because GRF of BMD from people of European descent may not contribute to BMD prediction in Chinese. This study highlights the limited utility of the current genetics studies largely focused on people of European descent for disease or risk factor prediction in other ethnic groups, and calls for more and larger scale studies focused on other ethnic groups.

Osteoporosis is a common metabolic bone disease. The impact of global blood plasma proteins on the risk of osteoporosis remains elusive now. We performed a large-scale polygenic risk score (PRS) analysis to evaluate the potential effects of blood plasma proteins on the development of osteoporosis in 2286 Caucasians, including 558 males and 1728 females. Bone mineral density (BMD) and bone areas at ulna & radius, hip, and spine were measured using Hologic 4500W DXA. BMD/bone areas values were adjusted for age, sex, height, and weight as covariates. Genome-wide SNP genotyping of 2286 Caucasian subjects was performed using Affymetrix Human SNP Array 6.0. The 267 blood plasma proteins-associated SNP loci and their genetic effects were obtained from recently published genome-wide association study (GWAS) using a highly multiplexed aptamer-based affinity proteomics platform. The polygenetic risk score (PRS) of study subjects for each blood plasma protein was calculated from the genotypes data of the 2286 Caucasian subjects by PLINK software. Pearson correlation analysis of individual PRS values and BMD/bone area value was performed using R. Additionally, gender-specific analysis also was performed by Pearson correlation analysis. 267 blood plasma proteins were analyzed in this study. For BMD, we observed association signals between 41 proteins and BMD, mainly including whole body total BMD versus Factor H (p value = 9.00 × 10 -3 ), whole body total BMD versus BGH3 (p value = 1.40 × 10 -2 ), spine total BMD versus IGF-I (p value = 2.15 × 10 -2 ), and spine total BMD versus SAP (p value = 3.90 × 10 -2 ). As for bone areas, association evidence was observed between 45 blood plasma proteins and bone areas, such as ferritin versus spine area (p value = 1.90 × 10 -2 ), C4 versus hip area (p value = 1.25 × 10 -2 ), and hemoglobin versus right ulna and radius area (p value = 2.70 × 10 -2 ). Our study results suggest the modest impact of blood plasma proteins on the variations of BMD/bone areas, and identify several candidate blood plasma proteins for osteoporosis.

Osteoporosis is a skeletal disorder characterized by low bone mineral density (BMD) and deterioration of bone microarchitecture. To identify novel genetic loci underlying osteoporosis, an effective strategy is to focus on scanning of variants with high potential functional impacts. Enhancers play a crucial role in regulating cell-type-specific transcription. Therefore, single-nucleotide polymorphisms (SNPs) located in enhancers (enhancer-SNPs) may represent strong candidate functional variants. Here, we performed a targeted analysis for potential functional enhancer-SNPs that may affect gene expression and biological processes in bone-related cells, specifically, osteoblasts, and peripheral blood monocytes (PBMs), using five independent cohorts (n = 5905) and the genetics factors for osteoporosis summary statistics, followed by comprehensive integrative genomic analyses of chromatin states, transcription, and metabolites. We identified 15 novel enhancer-SNPs associated with femoral neck and lumbar spine BMD, including 5 SNPs mapped to novel genes (e.g., rs10840343 and rs10770081 in IGF2 gene) and 10 novel SNPs mapped to known BMD-associated genes (e.g., rs2941742 in ESR1 gene, and rs10249092 and rs4342522 in SHFM1 gene). Interestingly, enhancer-SNPs rs10249092 and rs4342522 in SHFM1 were tightly linked, but annotated to different enhancers in PBMs and osteoblasts, respectively, suggesting that even tightly linked SNPs may regulate the same target gene and contribute to the phenotype variation in cell-type-specific manners. Importantly, ten enhancer-SNPs may also regulate BMD variation by affecting the serum metabolite levels. Our findings revealed novel susceptibility loci that may regulate BMD variation and provided intriguing insights into the genetic mechanisms of osteoporosis.

Osteoporosis is a prevalent bone metabolic disease, mainly caused by excessive bone resorption (by osteoclasts) over bone formation (by osteoblasts). Identifying the key transcription factors and understanding the regulatory network influencing osteoclastogenesis will be helpful to explore the potential biological mechanism for osteoporosis. In our study, peripheral blood monocyte (PBM) was used as a cell model for bone mineral density (BMD) research. PBMs serve as progenitors of osteoclasts and produce important cytokines for osteoclastogenesis. In our study, via exon arrays, gene expression profiles of PBMs were analyzed between high versus low hip BMD groups. Transcription factors for differentially expressed genes were then predicted based on the enrichment analysis. We found that 591 genes were differentially expressed between the two BMD groups (nominally significant, raw p value < 0.05). For high BMD subjects, 482 genes were up-regulated and 109 genes were down-regulated. We then found 29 potential transcription factors for up-regulated genes and nine transcription factors for down-regulated genes. Among these transcription factors, HMGA1 and NFKB2 were differentially expressed between high versus low BMD groups. In addition, their regulation types with their target genes were consistent with the information from public databases. Our findings of key transcription factors and their target genes for osteoporosis were further validated by GWAS analysis. Overall, we predicted important transcription factors for osteoporosis. We were also able to infer the regulatory mechanism that exists between transcription factors and target genes in bone metabolism.

Objective: The relationships between early growth parameters and bone mineral density (BMD) remain elusive now. In this study, we performed a large scale polygenic risk score (PRS) analysis to evaluate the potential impact of early growth parameters on the variations of BMD., Methods: We used 2286 Caucasian subjects as cohort 1 and 3404 Framingham Heart Study (FHS) subjects as cohort 2 in this study. BMD at ulna & radius, hip and spine were measured using dual energy X-ray absorptiometry. BMD values were adjusted for age, sex, height and weight as covariates. Genome-wide single-nucleotide polymorphism (SNP) genotyping of the 2286 Caucasian subjects was performed using Affymetrix Human SNP Array 6.0. The GWAS datasets of early growth parameters were driven from the Early Growth Genetics Consortium, including birth weight (BW), birth head circumference (BHC), childhood body mass index (CBMI), pubertal height growth related indexes and tanner stage. Polygenic Risk Score (PRSice) and linkage disequilibrium (LD) score regression analysis were conducted to assess the genetic correlation between early growth parameters and BMD., Results: We detected significant genetic correlations in cohort 1, such as total spine BMD vs. CBMI (p value = 1.51 × 10 -4 , rg = 0.4525), right ulna and radius BMD vs. CBMI (p value = 1.51 × 10 -4 , rg = 0.4399) and total body BMD vs. tanner stage (p value = 7.00 × 10 -4 , rg = -0.0721). For cohort 2, significant correlations were observed for total spine BMD vs. height change standard deviation score (SDS) between 8 years and adult (denoted as PGF + PGM) (p value = 3.97 × 10 -4 , rg = -0.1425), femoral neck BMD vs. the timing of peak height velocity by looking at the height change SDS between age 14 years and adult (denoted as PTF + PTM) (p value = 7.04 × 10 -4 , rg = -0.2185), and total spine BMD vs. PTF + PTM (p value = 6.86 × 10 -4 , rg = -0.2180)., Conclusion: Our study results suggest that some early growth parameters could affect the variations of BMD., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Several studies indicated bone mineral density (BMD) and alcohol intake might share common genetic factors. The study aimed to explore potential SNPs/genes related to both phenotypes in US Caucasians at the genome-wide level. A bivariate genome-wide association study (GWAS) was performed in 2069 unrelated participants. Regular drinking was graded as 1, 2, 3, 4, 5, or 6, representing drinking alcohol never, less than once, once or twice, three to six times, seven to ten times, or more than ten times per week respectively. Hip, spine, and whole body BMDs were measured. The bivariate GWAS was conducted on the basis of a bivariate linear regression model. Sex-stratified association analyses were performed in the male and female subgroups. In males, the most significant association signal was detected in SNP rs685395 in DYNC2H1 with bivariate spine BMD and alcohol drinking (P = 1.94 × 10 -8 ). SNP rs685395 and five other SNPs, rs657752, rs614902, rs682851, rs626330, and rs689295, located in the same haplotype block in DYNC2H1 were the top ten most significant SNPs in the bivariate GWAS in males. Additionally, two SNPs in GRIK4 in males and three SNPs in OPRM1 in females were suggestively associated with BMDs (of the hip, spine, and whole body) and alcohol drinking. Nine SNPs in IL1RN were only suggestively associated with female whole body BMD and alcohol drinking. Our study indicated that DYNC2H1 may contribute to the genetic mechanisms of both spine BMD and alcohol drinking in male Caucasians. Moreover, our study suggested potential pleiotropic roles of OPRM1 and IL1RN in females and GRIK4 in males underlying variation of both BMD and alcohol drinking.

Aiming to identify genomic variants associated with osteoporosis, we performed a genome-wide association meta-analysis of bone mineral density (BMD) at Ward's triangle of the hip in 7175 subjects from 6 samples. We performed in silico replications with femoral neck, trochanter, and inter-trochanter BMDs in 6912 subjects from the Framingham heart study (FHS), and with forearm, femoral neck and lumbar spine BMDs in 32965 subjects from the GEFOS summary results. Combining the evidence from all samples, we identified 2 novel loci for areal BMD: 1q43 (rs1414660, discovery p=1.20×10(-8), FHS p=0.05 for trochanter BMD; rs9287237, discovery p=3.55×10(-7), FHS p=9.20×10(-3) for trochanter BMD, GEFOS p=0.02 for forearm BMD, nearest gene FMN2) and 2q32.2 (rs56346965, discovery p=7.48×10(-7), FHS p=0.10 for inter-trochanter BMD, GEFOS p=0.02 for spine BMD, nearest gene NAB1). The two lead SNPs rs1414660 and rs56346965 are eQTL sites for the genes GREM2 and NAB1 respectively. Functional annotation of GREM2 and NAB1 illustrated their involvement in BMP signaling pathway and in bone development. We also replicated three previously reported loci: 5q14.3 (rs10037512, discovery p=3.09×10(-6), FHS p=8.50×10(-3), GEFOS p=1.23×10(-24) for femoral neck BMD, nearest gene MEF2C), 6q25.1 (rs3020340, discovery p=1.64×10(-6), GEFOS p=1.69×10(-3) for SPN-BMD, nearest gene ESR1) and 7q21.3 (rs13310130, discovery p=8.79×10(-7), GEFOS p=2.61×10(-7) for spine BMD, nearest gene SHFM1). Our findings provide additional insights that further enhance our understanding of bone development, osteoporosis, and fracture pathogenesis., (Copyright © 2016 Elsevier Inc. All rights reserved.)

To identify genetic variants that influence bone mineral density (BMD) in East Asians, we performed a quantitative trait analysis of lumbar spine, total hip and femoral neck BMD in a Korean population-based cohort (N=2729) and follow-up replication analysis in a Chinese Han population and two Caucasian populations (N=1547, 2250 and 987, respectively). From the meta-analysis of the stage 1 discovery analysis and stage 2 replication analysis, we identified four BMD loci that reached near-genome-wide significance level (P<5×10(-7)). One locus on 1q23 (UHMK1, rs16863247, P=4.1×10(-7) for femoral neck BMD and P=3.2×10(-6) for total hip BMD) was a novel BMD signal. Interestingly, rs16863247 was very rare in Caucasians (minor allele frequency<0.01), indicating that this association could be specific to East Asians. In gender specific analysis, rs1160574 on 1q32 (KCNH1) was associated with femoral neck BMD (P=2.1×10(-7)) in female subjects. rs9371538 in the known BMD region on 6q25 ESR1 was associated with lumbar spine BMD (P=5.6×10(-9)). rs7776725 in the known BMD region on 7q31 WTN16 was associated with total hip BMD (P=8.6×10(-9)). In osteoblasts, endogenous UHMK1 expression was increased during differentiation and UHMK1 knockdown decreased its differentiation, while UHMK1 overexpression increased its differentiation. In osteoclasts, endogenous UHMK1 expression was decreased during differentiation and UHMK1 knockdown increased its differentiation, while UHMK1 overexpression decreased its differentiation. In conclusion, our genome-wide association study identified the UHMK1 gene as a novel BMD locus specific to East Asians. Functional studies suggest a role of UHMK1 on regulation of osteoblasts and osteoclasts., (Copyright © 2016 Elsevier Inc. All rights reserved.)

To identify susceptibility genes for osteoporosis, we conducted an integrative analysis that combined epigenomic elements and previous genome-wide association studies (GWASs) data, followed by validation at population and functional levels, which could identify common regulatory elements and predict new susceptibility genes that are biologically meaningful to osteoporosis. By this approach, we found a set of distinct epigenomic elements significantly enriched or depleted in the promoters of osteoporosis-associated genes, including 4 transcription factor binding sites, 27 histone marks, and 21 chromatin states segmentation types. Using these epigenomic marks, we performed reverse prediction analysis to prioritize the discovery of new candidate genes. Functional enrichment analysis of all the prioritized genes revealed several key osteoporosis related pathways, including Wnt signaling. Genes with high priority were further subjected to validation using available GWASs datasets. Three genes were significantly associated with spine bone mineral density, including BDNF, PDE4D, and SATB2, which all closely related to bone metabolism. The most significant gene BDNF was also associated with osteoporotic fractures. RNA interference revealed that BDNF knockdown can suppress osteoblast differentiation. Our results demonstrated that epigenomic data could be used to indicate common epigenomic marks to discover additional loci with biological functions for osteoporosis.

Unlabelled: Osteoporosis is mainly characterized by low bone mineral density (BMD), and can be attributed to excessive bone resorption by osteoclasts. Migration of circulating monocytes from blood to bone is important for subsequent osteoclast differentiation and bone resorption. Identification of those genes and pathways related to osteoclastogenesis and BMD will contribute to a better understanding of the pathophysiological mechanisms of osteoporosis. In this study, we applied the LC-nano-ESI-MS(E) (Liquid Chromatograph-nano-Electrospray Ionization-Mass Spectrometry) for quantitative proteomic profiling in 33 female Caucasians with discordant BMD levels, with 16 high vs. 17 low BMD subjects. Protein quantitation was accomplished by label-free measurement of total ion currents collected from MS(E) data. Comparison of protein expression in high vs. low BMD subjects showed that ITGA2B (p=0.0063) and GSN (p=0.019) were up-regulated in the high BMD group. Additionally, our protein-RNA integrative analysis showed that RHOA (p=0.00062) differentially expressed between high vs. low BMD groups. Network analysis based on multiple tools revealed two pathways: "regulation of actin cytoskeleton" (p=1.13E-5, FDR=3.34E-4) and "leukocyte transendothelial migration" (p=2.76E-4, FDR=4.71E-3) that are functionally relevant to osteoporosis. Consistently, ITGA2B, GSN and RHOA played crucial roles in these two pathways respectively. All together, our study strongly supported the contribution of the genes ITGA2B, GSN and RHOA and the two pathways to osteoporosis risk., Biological Significance: Mass spectrometry based quantitative proteomics study integrated with network analysis identified novel genes and pathways related to osteoporosis. The results were further verified in multiple level studies including protein-RNA integrative analysis and genome wide association studies., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Protein phosphorylation regulates a wide variety of cellular processes. Thus, we hypothesize that single-nucleotide polymorphisms (SNPs) that may modulate protein phosphorylation could affect osteoporosis risk. Based on a previous conventional genome-wide association (GWA) study, we conducted a three-stage meta-analysis targeting phosphorylation-related SNPs (phosSNPs) for femoral neck (FN)-bone mineral density (BMD), total hip (HIP)-BMD, and lumbar spine (LS)-BMD phenotypes. In stage 1, 9593 phosSNPs were meta-analyzed in 11,140 individuals of various ancestries. Genome-wide significance (GWS) and suggestive significance were defined by α = 5.21 × 10(-6) (0.05/9593) and 1.00 × 10(-4), respectively. In stage 2, nine stage 1-discovered phosSNPs (based on α = 1.00 × 10(-4)) were in silico meta-analyzed in Dutch, Korean, and Australian cohorts. In stage 3, four phosSNPs that replicated in stage 2 (based on α = 5.56 × 10(-3), 0.05/9) were de novo genotyped in two independent cohorts. IDUA rs3755955 and rs6831280, and WNT16 rs2707466 were associated with BMD phenotypes in each respective stage, and in three stages combined, achieving GWS for both FN-BMD (p = 8.36 × 10(-10), p = 5.26 × 10(-10), and p = 3.01 × 10(-10), respectively) and HIP-BMD (p = 3.26 × 10(-6), p = 1.97 × 10(-6), and p = 1.63 × 10(-12), respectively). Although in vitro studies demonstrated no differences in expressions of wild-type and mutant forms of IDUA and WNT16B proteins, in silico analyses predicts that WNT16 rs2707466 directly abolishes a phosphorylation site, which could cause a deleterious effect on WNT16 protein, and that IDUA phosSNPs rs3755955 and rs6831280 could exert indirect effects on nearby phosphorylation sites. Further studies will be required to determine the detailed and specific molecular effects of these BMD-associated non-synonymous variants., (© 2015 American Society for Bone and Mineral Research.)

Runs of homozygosity (ROHs), in which both parental alleles are identical, have been proposed to have recessive effects on multiple human complex diseases. Osteoporosis is a common complex disease characterized by low bone mineral density (BMD), which is highly heritable. And recessive loci that contribute to BMD variations have been identified. In this study, we performed a genome-wide ROHs association study using our SNP array data from three GWAS samples including 4,900 subjects from Caucasian and Chinese populations. Significant results were further subjected to replication in 3,747 additional subjects. ROHs associated with BMD were also tested for associations with osteoporotic fractures in a GWAS fracture sample. Combining results from all the samples, we identified 697 autosomal regions with ROHs. Among these, we detected genome-wide significant associations between BMD and 6 ROHs, including ROH1q31.3, 1p31.1, 3q26.1, 11q12.1, 21q22.1 and 15q22.3 (combined P = 6.29 × 10(-5)-3.17 × 10(-8)). Especially, ROH1p31.1 was found to be associated with increased risk of osteoporotic hip fractures (odds ratio [OR] = 3.71, P = 0.032). To investigate the functional relevance of the identified ROHs, we performed cis-expression quantitative trait locus (eQTL) analysis in lymphoblast cell lines. Three ROHs, including ROH1p31.1, 11q12.1, and 15q22.3, were found to be significantly associated with mRNA expression levels of their nearby genes (PeQTL < 0.05). In summary, our findings reveal that ROHs could play as recessive-acting determinants contributing to the pathogenesis of osteoporosis. Further molecular and functional studies are needed to explore and clarify the potential mechanism., (© 2015 American Society for Bone and Mineral Research.)

Integration of multiple profiling data and construction of functional gene networks may provide additional insights into the molecular mechanisms of complex diseases. Osteoporosis is a worldwide public health problem, but the complex gene-gene interactions, post-transcriptional modifications and regulation of functional networks are still unclear. To gain a comprehensive understanding of osteoporosis etiology, transcriptome gene expression microarray, epigenomic miRNA microarray and methylome sequencing were performed simultaneously in 5 high hip BMD (Bone Mineral Density) subjects and 5 low hip BMD subjects. SPIA (Signaling Pathway Impact Analysis) and PCST (Prize Collecting Steiner Tree) algorithm were used to perform pathway-enrichment analysis and construct the interaction networks. Through integrating the transcriptomic and epigenomic data, firstly we identified 3 genes (FAM50A, ZNF473 and TMEM55B) and one miRNA (hsa-mir-4291) which showed the consistent association evidence from both gene expression and methylation data; secondly in network analysis we identified an interaction network module with 12 genes and 11 miRNAs including AKT1, STAT3, STAT5A, FLT3, hsa-mir-141 and hsa-mir-34a which have been associated with BMD in previous studies. This module revealed the crosstalk among miRNAs, mRNAs and DNA methylation and showed four potential regulatory patterns of gene expression to influence the BMD status. In conclusion, the integration of multiple layers of omics can yield in-depth results than analysis of individual omics data respectively. Integrative analysis from transcriptomics and epigenomic data improves our ability to identify causal genetic factors, and more importantly uncover functional regulation pattern of multi-omics for osteoporosis etiology.

MicroRNAs (miRNAs) are critical post-transcriptional regulators. Based on a previous genome-wide association (GWA) scan, we conducted a polymorphism in microRNA target sites (poly-miRTS)-centric multistage meta-analysis for lumbar spine (LS)-, total hip (HIP)- and femoral neck (FN)-bone mineral density (BMD). In stage I, 41 102 poly-miRTSs were meta-analyzed in seven cohorts with a genome-wide significance (GWS) α = 0.05/41 102 = 1.22 × 10(-6). By applying α = 5 × 10(-5) (suggestive significance), 11 poly-miRTSs were selected, with FGFRL1 rs4647940 and PRR5 rs3213550 as top signals for FN-BMD (P = 7.67 × 10(-6) and 1.58 × 10(-5)) in gender-combined sample. In stage II in silico replication (two cohorts), FGFRL1 rs4647940 was the only signal marginally replicated for FN-BMD (P = 5.08 × 10(-3)) at α = 0.10/11 = 9.09 × 10(-3). PRR5 rs3213550 was also selected based on biological significance. In stage III de novo genotyping replication (two cohorts), FGFRL1 rs4647940 was the only signal significantly replicated for FN-BMD (P = 7.55 × 10(-6)) at α = 0.05/2 = 0.025 in gender-combined sample. Aggregating three stages, FGFRL1 rs4647940 was the single stage I-discovered and stages II- and III-replicated signal attaining GWS for FN-BMD (P = 8.87 × 10(-12)). Dual-luciferase reporter assays demonstrated that FGFRL1 3' untranslated region harboring rs4647940 appears to be hsa-miR-140-5p's target site. In a zebrafish microinjection experiment, dre-miR-140-5p is shown to exert a dramatic impact on craniofacial skeleton formation. Taken together, we provided functional evidence for a novel FGFRL1 poly-miRTS rs4647940 in a previously known 4p16.3 locus, and experimental and clinical genetics studies have shown both FGFRL1 and hsa-miR-140-5p are important for bone formation., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Background: Osteoporosis is caused by excessive bone resorption (by osteoclasts) over bone formation (by osteoblasts). Monocytes are important to osteoporosis by serving as progenitors of osteoclasts and produce cytokines for osteoclastogenesis., Aim: To identify osteoporosis-related genes, we performed microarray analyses of monocytes using Affymetrix 1.0 ST arrays in 42 (including 16 pre- and 26 postmenopausal) high hip BMD (bone mineral density) vs. 31 (including 15 pre- and 16 postmenopausal) low hip BMD Caucasian female subjects. Here, high vs. low BMD is defined as belonging to top vs. bottom 30% of BMD values in population., Method: Differential gene expression analysis in high vs. low BMD subjects was conducted in the total cohort as well as pre- and post-menopausal subjects. Focusing on the top differentially expressed genes identified in the total, the pre- and the postmenopausal subjects (with a p <5E-03), we performed replication of the findings in 3 independent datasets of microarray analyses of monocytes (total N = 125)., Results: We identified (in the 73 subjects) and successfully replicated in all the 3 independent datasets 2 genes, DAXX and PLK3. Interestingly, both genes are apoptosis induction genes and both down-regulated in the low BMD subjects. Moreover, using the top 200 genes identified in the meta-analysis across all of the 4 microarray datasets, GO term enrichment analysis identified a number of terms related to induction of apoptosis, for which the majority of component genes are also down-regulated in the low BMD subjects. Overall, our result may suggest that there might be a decreased apoptosis activity of monocytes in the low BMD subjects., Conclusion: Our study for the first time suggested a decreased apoptosis rate (hence an increased survival) of monocytes, an important osteoclastogenic cell, as a novel mechanism for osteoporosis.

DDR2 gene, playing an essential role in regulating osteoblast differentiation and chondrocyte maturation, may influence bone mineral density (BMD) and osteoporosis, but the genetic variations actually leading to the association remain to be elucidated. Therefore, the aim of this study was to investigate whether the genetic variants in DDR2 are associated with BMD and fracture risk. This study was performed in three samples from two ethnicities, including 1,300 Chinese Han subjects, 700 Chinese Han subjects (350 with osteoporotic hip fractures and 350 healthy controls) and 2,286 US white subjects. Twenty-eight SNPs in DDR2 were genotyped and tested for associations with hip BMD and fractures. We identified 3 SNPs in DDR2 significantly associated with hip BMD in the Chinese population after multiple testing adjustments, which were rs7521233 (P = 1.06×10-4, β: -0.018 for allele C), rs7553831 (P = 1.30×10-4, β: -0.018 for allele T), and rs6697469 (P = 1.59×10-3, β: -0.015 for allele C), separately. These three SNPs were in high linkage disequilibrium. Haplotype analyses detected two significantly associated haplotypes, including one haplotype in block 2 (P = 9.54×10-4, β: -0.016) where these three SNPs located. SNP rs6697469 was also associated with hip fractures (P = 0.043, OR: 1.42) in the Chinese population. The effect on fracture risk was consistent with its association with lower BMD. However, in the white population, we didn't observe significant associations with hip BMD. eQTL analyses revealed that SNPs associated with BMD also affected DDR2 mRNA expression levels in Chinese. Our findings, together with the prior biological evidence, suggest that DDR2 could be a new candidate for osteoporosis in Chinese population. Our results also reveal an ethnic difference, which highlights the need for further genetic studies in each ethnic group.

Low bone mineral density (BMD) is a risk factor for osteoporosis. Osteoporosis is more prevalent in females than in males. So far, the pathophysiological mechanisms underlying osteoporosis are unclear. Peripheral blood monocytes (PBMs) are precursors of bone-resorbing osteoclasts. This study aims to identify PBM-expressed proteins (genes) influencing hip BMD in humans. We utilized three independent study cohorts (N=34, 29, 40), including premenopausal Caucasians with discordant hip BMD. We studied PBM proteome-wide protein expression profiles in cohort 1 and identified 57 differentially expressed proteins (DEPs) between low vs. high BMD subjects. One protein gelsolin (GSN), after validation by Western blotting, was subject to follow-up. We compared GSN mRNA level in PBM between low vs. high BMD subjects in cohorts 2 and 3. We genotyped SNPs across GSN in 2286 unrelated Caucasians (cohort 4) and 1627 Chinese (cohort 5) and tested their association with hip BMD in females and males, respectively. We discovered and validated that GSN protein expression level in PBM was down-regulated 3.0-fold in low vs. high BMD subjects (P<0.05). Down-regulation of GSN in PBM in low BMD subjects was also observed at mRNA level in both cohort 2 and cohort 3. We identified that SNP rs767770 was significantly associated with hip BMD in female Caucasians (P=0.0003) only. Integrating analyses of the datasets at DNA, RNA, and protein levels from female Caucasians substantiated that GSN is highly significant for hip BMD (P=0.0001). We conclude that GSN is a significant gene influencing hip BMD in female Caucasians., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10(-8)) level: 14q24.2 (rs227425, P-value 3.98 × 10(-13), SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10(-9), CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n = 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis.

Bone mineral density (BMD) is a major index for diagnosing osteoporosis. PhosSNPs are nonsynonymous SNPs that affect protein phosphorylation. The relevance and significance of phosSNPs to BMD and osteoporosis is unknown. This study aimed to identify and characterize phosSNPs significant for BMD in humans. We conducted a pilot genomewide phosSNP association study for BMD in three independent population samples, involving ∼5000 unrelated individuals. We identified and replicated three phosSNPs associated with both spine BMD and hip BMD in Caucasians. Association with hip BMD for one of these phosSNPs, ie, rs6265 (major/minor allele: G/A) in BDNF gene, was also suggested in Chinese. Consistently in both ethnicities, individuals carrying the AA genotype have significantly lower hip BMD than carriers of the GA and GG genotypes. Through in vitro molecular and cellular studies, we found that compared to osteoblastic cells transfected with wild-type BDNF-Val66 (encoded with allele G at rs6265), transfection of variant BDNF-Met66 (encoded with allele A at rs6265) significantly decreased BDNF protein phosphorylation (at amino acid residue T62), expression of osteoblastic genes (OPN, BMP2, and ALP), and osteoblastic activity. The findings are consistent with and explain our prior observations in general human populations. We conclude that phosSNP rs6265, by regulating BDNF protein phosphorylation and osteoblast differentiation, influences hip BMD in humans. This study represents our first endeavor to dissect the functions of phosSNPs in bone, which might stimulate extended large-scale studies on bone or similar studies on other human complex traits and diseases., (© 2013 American Society for Bone and Mineral Research.)

Osteoporosis is characterized by low bone mineral density (BMD), a highly heritable trait that is determined, in part, by the actions and interactions of multiple genes. Although an increasing number of genes have been identified to have independent effects on BMD, few studies have been performed to identify genes that interact with one another to affect BMD. In this study, we performed gene-gene interaction analyses in selected candidate genes in individuals with extremely high versus low hip BMD (20% tails of the distributions), in two independent U.S. Caucasian samples. The first sample contained 916 unrelated subjects with extreme hip BMD Z-scores selected from a population composed of 2286 subjects. The second sample consisted of 400 unrelated subjects with extreme hip BMD Z-scores selected from a population composed of 1000 subjects. Combining results from these two samples, we found one interacting gene pair (RBMS3 versus ZNF516) which, even after Bonferroni correction for multiple testing, showed consistently significant effects on hip BMD. RMBS3 harbored two single-nucleotide polymorphisms (SNPs), rs6549904 and rs7640046, both of which had significant interactions with an SNP, rs4891159, located on ZNF516 (p = 7.04 × 10(-11) and 1.03 × 10(-10) ). We further validated these results in two additional samples of Caucasian and African descent. The gene pair, RBMS3 versus ZNF516, was successfully replicated in the Caucasian sample (p = 8.07 × 10(-3) and 2.91 × 10(-3) ). For the African sample, a significant interaction was also detected (p = 0.031 and 0.043), but the direction of the effect was opposite to that observed in the three Caucasian samples. By providing evidence for genetic interactions underlying BMD, this study further delineates the genetic architecture of osteoporosis., (Copyright © 2013 American Society for Bone and Mineral Research.)

Bone mineral density (BMD) measured at the femoral neck (FN) is the most important risk phenotype for osteoporosis and has been used as a reference standard for describing osteoporosis. The specific genes influencing FN BMD remain largely unknown. To identify such genes, we first performed a genome-wide association (GWA) analysis for FN BMD in a discovery sample consisting of 983 unrelated white subjects. We then tested the top significant single-nucleotide polymorphisms (SNPs; 175 SNPs with p < 5 x 10(-4)) for replication in a family-based sample of 2557 white subjects. Combing results from these two samples, we found that two genes, parathyroid hormone (PTH) and interleukin 21 receptor (IL21R), achieved consistent association results in both the discovery and replication samples. The PTH gene SNPs, rs9630182, rs2036417, and rs7125774, achieved p values of 1.10 x 10(-4), 3.24 x 10(-4), and 3.06 x 10(-4), respectively, in the discovery sample; p values of 6.50 x 10(-4), 5.08 x 10(-3), and 5.68 x 10(-3), respectively, in the replication sample; and combined p values of 3.98 x 10(-7), 9.52 x 10(-6), and 1.05 x 10(-5), respectively, in the total sample. The IL21R gene SNPs, rs8057551, rs8061992, and rs7199138, achieved p values of 1.51 x 10(-4), 1.53 x 10(-4), and 3.88 x 10(-4), respectively, in the discovery sample; p values of 2.36 x 10(-3), 6.74 x 10(-3), and 6.41 x 10(-3), respectively, in the replication sample; and combined p values of 2.31 x 10(-6), 8.62 x 10(-6), and 1.41 x 10(-5), respectively, in the total sample. The effect size of each SNP was approximately 0.11 SD estimated in the discovery sample. PTH and IL21R both have potential biologic functions important to bone metabolism. Overall, our findings provide some new clues to the understanding of the genetic architecture of osteoporosis., ((c) 2010 American Society for Bone and Mineral Research.)

Both bone mineral density (BMD) and lean body mass (LBM) are important physiological measures with strong genetic determination. Besides, BMD and LBM might have common genetic factors. Aiming to identify pleiotropic genomic loci underlying BMD and LBM, we performed bivariate genome-wide association study meta-analyses of femoral neck bone mineral density and LBM at arms and legs, and replicated in the large-scale UK Biobank cohort sample. Combining the results from discovery meta-analysis and replication sample, we identified three genomic loci at the genome-wide significance level (p < 5.0 × 10–8): 2p23.2 (lead SNP rs4477866, discovery p = 3.47 × 10–8, replication p = 1.03 × 10–4), 16q12.2 (rs1421085, discovery p = 2.04 × 10–9, replication p = 6.47 × 10–14) and 18q21.32 (rs11152213, discovery p = 3.47 × 10–8, replication p = 6.69 × 10–6). Our findings not only provide useful insights into lean mass and bone mass development, but also enhance our understanding of the potential genetic correlation between BMD and LBM. [ABSTRACT FROM AUTHOR]

publisher‐imprint‐name Springer volume‐issue‐count 1 issue‐article‐count 20 issue‐toc‐levels 0 issue‐pricelist‐year 2020 issue‐copyright‐holder The Author(s) issue‐copyright‐year 2020 article‐contains‐esm Yes article‐numbering‐style Unnumbered article‐registration‐date‐year 2020 article‐registration‐date‐month 2 article‐registration‐date‐day 17 article‐toc‐levels 0 toc‐levels 0 volume‐type Regular journal‐product ArchiveJournal numbering‐style Unnumbered article‐grants‐type OpenChoice metadata‐grant OpenAccess abstract‐grant OpenAccess bodypdf‐grant OpenAccess bodyhtml‐grant OpenAccess bibliography‐grant OpenAccess esm‐grant OpenAccess online‐first false pdf‐file‐reference BodyRef/PDF/40169_2020_Article_272.pdf pdf‐type Typeset target‐type OnlinePDF issue‐type Regular article‐type OriginalPaper journal‐subject‐primary Medicine & Public Health journal‐subject‐secondary Medicine/Public Health, general journal‐subject‐collection Medicine open‐access true --> Backgrounds: To explore the genetic correlation between schizophrenia (SCZ) and osteoporosis (OP). Design, setting, participants, measurements: We conducted a trans‐ethnic two‐stage genetic correlation analysis of OP and SCZ, totally invoking 2286 Caucasia subjects in discovery stage and 4124 Chinese subjects in replication stage. The bone mineral density (BMD) and bone area values of ulna & radius, hip and spine were measured using Hologic 4500W dual energy X‐ray absorptiometry machine. SCZ was diagnosed according to DSM‐IV criteria. For the genome‐wide association study (GWAS) of Caucasian OP, Chinese OP and Chinese SCZ, SNP genotyping was performed using Affymetrix SNP 6.0 array. For the GWAS of Caucasian SCZ, SNP genotyping was conducted using the Affymetrix 5.0 array, Affymetrix 6.0 array and Illumina 550 K array. Polygenetic risk scoring (PRS) analysis was conducted by PRSice software. Also, Linkage disequilibrium score regression (LD Score regression) analysis was performed to evaluate the genetic correlation between OP and SCZ. Multi‐trait analysis of GWAS (MTAG) was performed to detect novel candidate genes for osteoporosis and SCZ. Results: In the Caucasia discovery samples, significant genetic correlations were observed for ulna & radius BMD vs. SCZ (P value = 0.010), ulna & radius area vs. SCZ (P value = 0.031). In the Chinese replication samples, we observed significant correlation for ulna & radius area vs. SCZ (P value = 0.019). In addition, LD Score regression also identified significant genetic correlations between SCZ and bone phenotypes in Caucasian and Chinese sample respectively. MTAG analysis identified several novel candidate genes, such as CTNNA2 (MTAG P value = 2.24 × 10−6) for SCZ and FADS2 (MTAG P value = 2.66 × 10−7) for osteoporosis. Conclusions: Our study results support the overlapped genetic basis for osteoporosis and SCZ, and provide novel clues for elucidating the biological mechanism of increased osteoporosis risk in SCZ patients. [ABSTRACT FROM AUTHOR]

Context: Though genome-wide association studies (GWASs) have identified hundreds of genetic variants associated with osteoporosis related traits, such as bone mineral density (BMD) and fracture, it remains a challenge to interpret their biological functions and underlying biological mechanisms.Objective: Integrate diverse expression quantitative trait loci and splicing quantitative trait loci data with several powerful GWAS datasets to identify novel candidate genes associated with osteoporosis.Design, Setting, and Participants: Here, we conducted a transcriptome-wide association study (TWAS) for total body BMD (TB-BMD) (n = 66 628 for discovery and 7697 for validation) and fracture (53 184 fracture cases and 373 611 controls for discovery and 37 857 cases and 227 116 controls for validation), respectively. We also conducted multi-SNP-based summarized mendelian randomization analysis to further validate our findings.Results: In total, we detected 88 genes significantly associated with TB-BMD or fracture through expression or ribonucleic acid splicing. Summarized mendelian randomization analysis revealed that 78 of the significant genes may have potential causal effects on TB-BMD or fracture in at least 1 specific tissue. Among them, 64 genes have been reported in previous GWASs or TWASs for osteoporosis, such as ING3, CPED1, and WNT16, as well as 14 novel genes, such as DBF4B, GRN, TMUB2, and UNC93B1.Conclusions: Overall, our findings provide novel insights into the pathogenesis mechanisms of osteoporosis and highlight the power of a TWAS to identify and prioritize potential causal genes. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Introduction: Osteoporosis is a prevalent bone metabolic disease characterized by bone fragility. As a key pathophysiological mechanism, the disease is caused by excessive bone resorption (by osteoclasts) over bone formation (by osteoblasts). Peripheral blood monocytes (PBMs) is a major systemic cell model for bone metabolism by serving as progenitors of osteoclasts and producing cytokines important for osteoclastogenesis. Protein-coding genes for osteoporosis have been widely studied by mRNA analyses of PBMs in high versus low hip bone mineral density (BMD) subjects. However, long noncoding RNAs (lncRNAs), which account for a large proportion of human transcriptome, have seldom been studied.Methods: In this study, microarray analyses of monocytes were performed using Affymetrix exon 1.0 ST arrays in 73 Caucasian females (age: 47-56). LncRNA profile was generated by re-annotating exon array for lncRNAs detection, which yielded 12,007 lncRNAs mapped to the human genome.Results: 575 lncRNAs were differentially expressed between the two groups. In the high BMD subjects, 309 lncRNAs were upregulated and 266 lncRNAs were downregulated (nominally significant, raw p-value < 0.05). To investigate the relationship between mRNAs and lncRNAs, we used two approaches to predict the target genes of lncRNAs and found that 26 candidate lncRNAs might regulate mRNA expression. The majority of these lncRNAs were further validated to be potentially correlated with BMD by GWAS analysis.Conclusion: Overall, our findings for the first time reported the lncRNAs profiles for osteoporosis and suggested the potential regulatory mechanism of lncRNAs on protein-coding genes in bone metabolism. [ABSTRACT FROM AUTHOR]

Summary Mitochondrial DNA (mtDNA) variants are involved in the pathogenesis of human complex diseases, especially for age-related disorders, including osteoporosis. However, the role of mtDNA variants in osteoporosis is largely unknown. In this study, we performed a mitochondria-wide association study for osteoporosis in a large sample of 2286 unrelated Caucasian subjects. A total of 445 mtSNPs were genotyped and 72 mtSNPs survived the quality control. We first examined association between mtSNPs and bone mineral density (BMD), and identified that an mtSNP, mt4823 within the ND2 gene, was strongly associated with hip BMD ( P= 2.05 × 10−4), even after Bonferroni correction. The C allele of mt4823 was associated with reduced hip BMD and the effect size (β) was ∼0.044. Another SNP mt15885 within the MT-CYB gene was associated both with spine ( P= 1.66 × 10−3) and hip BMD ( P= 0.023). The T allele of mt15885 had a protective effect on spine (β= 0.064) and hip BMD (β= 0.038). Next, we classified subjects into the nine common European haplogroups and conducted association analyses. Subjects classified as haplogroup X had significantly lower hip BMD values than others ( P= 0.040). Our results highlighted the importance of mtDNA variants in osteoporosis. [ABSTRACT FROM AUTHOR]

Bivariate genome-wide association analysis identified three pleiotropic loci underlying osteoporosis and obesity Osteoporosis and obesity are two inheritable diseases indexed by low bone mineral density (BMD) and high body mass index (BMI), and are correlated. 2 Abbreviations: BMD, bone mineral density; BMI, body mass index; GWA, genome-wide association. [Extracted from the article]

Background: Individuals’ peak bone mineral density (BMD) achieved and maintained at ages 20–40 years is the most powerful predictor of low bone mass and osteoporotic fractures later in life. The aim of this study was to identify metabolomic factors associated with peak BMD variation in US Caucasian women. Methods: A total of 136 women aged 20–40 years, including 65 subjects with low and 71 with high hip BMD, were enrolled. The serum metabolites were assessed using a liquid chromatography-mass spectrometry (LC-MS) method. The partial least-squares discriminant analysis (PLS-DA) method and logistic regression models were used, respectively, to examine the associations of metabolomic profiles and individual metabolites with BMD. Results: The low and high BMD groups could be differentiated by the detected serum metabolites using PLS-DA (Ppermutation = 0.008). A total of 14 metabolites, including seven amino acids and amino acid derivatives, five lipids (including three bile acids), and two organic acids, were significantly associated with the risk for low BMD. Most of these metabolites are novel in that they have never been linked with BMD in humans earlier. The prediction model including the newly identified metabolites significantly improved the classification of the groups with low and high BMD. The area under the receiver operating characteristic curve without and with metabolites were 0.88 (95% CI: 0.83–0.94) and 0.97 (95% CI: 0.94–0.99), respectively (P for the difference = 0.0004). Conclusion: Metabolomic profiling may improve the risk prediction of osteoporosis among Caucasian women. Our findings also suggest the potential importance of the metabolism of amino acids and bile acids in bone health. [ABSTRACT FROM AUTHOR]

In the present study, aiming to identify loci associated with osteoporosis, we conducted a joint association study of 2 independent genome-wide association meta-analyses of femoral neck and lumbar spine bone mineral densities (BMDs): 1) an in-house study of 6 samples involving 7484 subjects, and 2) the GEFOS-seq study of 7 samples involving 32,965 subjects. The in-house samples were imputed by the 1000 genomes project phase 3 reference panel. SNP-based association test was applied to 7,998,108 autosomal SNPs in each meta-analysis, and for each SNP the 2 association signals were then combined for joint analysis and for mutual replication. Combining the evidence from both studies, we identified 2 novel loci associated with BMDs at the genome-wide significance level (α = 5.0 × 10 −8 ): 20p12.1 ( rs73100693 p = 2.65 × 10 −8 , closest gene MACROD2 ) and 20q13.33 ( rs2380128 p = 3.44 × 10 −8 , OSBPL2 ). We also replicated 7 loci that were reported by two recent studies on heel and total body BMD. Our findings provide useful insights that enhance our understanding of bone development, osteoporosis and fracture pathogenesis. [ABSTRACT FROM AUTHOR]