Osteoporosis is characterized by low bone mineral density (BMD). The advancement of highthroughput technologies and integrative approaches provided an opportunity for deciphering the mechanisms underlying osteoporosis. Here, we generated genomic, transcriptomic, methylomic, and metabolomic datasets from 119 subjects with high (n = 61) and low (n = 58) BMDs. By adopting sparse multiple discriminative canonical correlation analysis, we identified an optimal multi-omics biomarker panel with 74 differentially expressed genes (DEGs), 75 differentially methylated CpG sites (DMCs), and 23 differential metabolic products (DMPs). By linking genetic data, we identified 199 targeted BMD-associated expression/methylation/metabolite quantitative trait loci (eQTLs/meQTLs/ metaQTLs). The reconstructed networks/pathways showed extensive biomarker interactions, and a substantial proportion of these biomarkers were enriched in RANK/RANKL, MAPK/TGF-b, and WNT/b-catenin pathways and G-protein-coupled receptor, GTP-binding/GTPase, telomere/mitochondrial activities that are essential for bone metabolism. Five biomarkers (FADS2, ADRA2A, FMN1, RABL2A, SPRY1) revealed causal effects on BMD variation. Our study provided an innovative framework and insights into the pathogenesis of osteoporosis.

A major challenge in translating findings from genome-wide association studies (GWAS) to biological mechanisms is pinpointing functional variants because only a very small percentage of variants associated with a given trait actually impact the trait. We used an extensive epigenetics, transcriptomics, and genetics analysis of the TBX15/WARS2 neighbourhood to prioritize this region’s best-candidate causal variants for the genetic risk of osteoporosis (estimated bone density, eBMD) and obesity (waist-hip ratio or waist circumference adjusted for body mass index). TBX15 encodes a transcription factor that is important in bone development and adipose biology. Manual curation of 692 GWAS-derived variants gave eight strong candidates for causal SNPs that modulate TBX15 transcription in subcutaneous adipose tissue (SAT) or osteoblasts, which highly and specifically express this gene. None of these SNPs were prioritized by Bayesian fine-mapping. The eight regulatory causal SNPs were in enhancer or promoter chromatin seen preferentially in SAT or osteoblasts at TBX15 intron-1 or upstream. They overlap strongly predicted, allele-specific transcription factor binding sites. Our analysis suggests that these SNPs act independently of two missense SNPs in TBX15 . Remarkably, five of the regulatory SNPs were associated with eBMD and obesity and had the same trait-increasing allele for both. We found that WARS2 obesity-related SNPs can be ascribed to high linkage disequilibrium with TBX15 intron-1 SNPs. Our findings from GWAS index, proxy, and imputed SNPs suggest that a few SNPs, including three in a 0.7-kb cluster, act as causal regulatory variants to fine-tune TBX15 expression and, thereby, affect both obesity and osteoporosis risk.

DNA methylation is an important epigenetic modification that contributes to the lineage commitment and specific functions of different cell types. In this study, we compared ENCODE-generated genome-wide DNA methylation profiles of human osteoblast with 21 other types of human cells in order to identify osteoblast-specific methylation events. For most of the cell strains, data from two isogenic replicates were included, resulting in a total of 51 DNA methylation datasets. We identified 852 significant osteoblast-specific differentially methylated CpGs (DMCs) and 295 significant differentially methylated regions (DMRs). Significant DMCs/DMRs were not enriched in CpG islands (CGIs) and promoters, but more strongly enriched in CGI shores/shelves and in gene body and intergenic regions. The genes associated with significant DMRs were highly enriched in biological processes related to transcriptional regulation and critical for regulating bone metabolism and skeletal development under physiologic and pathologic conditions. By integrating the DMR data with the extensive gene expression and chromatin epigenomics data, we observed complex, context-dependent relationships between DNA methylation, chromatin states, and gene expression, suggesting diverse DNA methylation-mediated regulatory mechanisms. Our results also highlighted a number of novel osteoblast-relevant genes. For example, the integrated evidences from DMR analysis, histone modification and RNA-seq data strongly support that there is a novel isoform of neurexin-2 (NRXN2) gene specifically expressed in osteoblast. NRXN2 was known to function as a cell adhesion molecule in the vertebrate nervous system, but its functional role in bone is completely unknown and thus worth further investigation. In summary, we reported a comprehensive analysis of osteoblast-specific DNA methylation profiles and revealed novel insights into the epigenetic basis of osteoblast differentiation and activity.

Background DNA methylation is an important heritable epigenetic mark that plays a crucial role in transcriptional regulation and the pathogenesis of various human disorders. The commonly used DNA methylation measurement approaches, e.g., Illumina Infinium HumanMethylation-27 and -450 BeadChip arrays (27 K and 450 K arrays) and reduced representation bisulfite sequencing (RRBS), only cover a small proportion of the total CpG sites in the human genome, which considerably limited the scope of the DNA methylation analysis in those studies. Results We proposed a new computational strategy to impute the methylation value at the unmeasured CpG sites using the mixture of regression model (MRM) of radial basis functions, integrating information of neighboring CpGs and the similarities in local methylation patterns across subjects and across multiple genomic regions. Our method achieved a better imputation accuracy over a set of competing methods on both simulated and empirical data, particularly when the missing rate is high. By applying MRM to an RRBS dataset from subjects with low versus high bone mineral density (BMD), we recovered methylation values of ~ 300 K CpGs in the promoter regions of chromosome 17 and identified some novel differentially methylated CpGs that are significantly associated with BMD. Conclusions Our method is well applicable to the numerous methylation studies. By expanding the coverage of the methylation dataset to unmeasured sites, it can significantly enhance the discovery of novel differential methylation signals and thus reveal the mechanisms underlying various human disorders/traits.

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 finemapping 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.