Osteoporosis is a highly prevalent disorder characterized by low bone mineral density and an increased risk of fracture, termed osteoporotic fracture. Notably, bone mineral density, osteoporosis and osteoporotic fracture are highly heritable; however, determining the genetic architecture, and especially the underlying genomic and molecular mechanisms, of osteoporosis in vivo in humans is still challenging. In addition to susceptibility loci identified in genome-wide association studies, advances in various omics technologies, including genomics, transcriptomics, epigenomics, proteomics and metabolomics, have all been applied to dissect the pathogenesis of osteoporosis. However, each technology individually cannot capture the entire view of the disease pathology and thus fails to comprehensively identify the underlying pathological molecular mechanisms, especially the regulatory and signalling mechanisms. A change to the status quo calls for integrative multi-omics and inter-omics analyses with approaches in 'systems genetics and

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.

Microsatellites are a major type of molecular markers in genetics studies. Their mutational dynamics are not clear. We investigated the patterns and characteristics of 97 mutation events unambiguously identified, from 53 multigenerational pedigrees with 630 subjects, at 362 autosomal dinucleotide microsatellite loci. A size-dependent mutation bias (in which long alleles are biased toward contraction, whereas short alleles are biased toward expansion) is observed. There is a statistically significant negative relationship between the magnitude (repeat numbers changed during mutation) and direction (contraction or expansion) of mutations and standardized allele size. Contrasting with earlier findings in humans, most mutation events (63%) in our study are multistep events that involve changes of more than one repeat unit. There was no correlation between mutation rate and recombination rate. Our data indicate that mutational dynamics at microsatellite loci are more complicated than the generalized stepwise mutation models.

Several genome-wide scans have been performed to detect loci that regulate BMD, but these have yielded inconsistent results, with limited replication of linkage peaks in different studies. In an effort to improve statistical power for detection of these loci, we performed a meta-analysis of genome-wide scans in which spine or hip BMD were studied. Evidence was gained to suggest that several chromosomal loci regulate BMD in a site-specific and sex-specific manner.Introduction: BMD is a heritable trait and an important predictor of osteoporotic fracture risk. Several genome-wide scans have been performed in an attempt to detect loci that regulate BMD, but there has been limited replication of linkage peaks between studies. In an attempt to resolve these inconsistencies, we conducted a collaborative meta-analysis of genome-wide linkage scans in which femoral neck BMD (FN-BMD) or lumbar spine BMD (LS-BMD) had been studied. Materials and Methods: Data were accumulated from nine genome-wide scans involving 11,842 subjects. Data were analyzed separately for LS-BMD and FN-BMD and by sex. For each study, genomic bins of 30 cM were defined and ranked according to the maximum LOD score they contained. While various densitometers were used in different studies, the ranking approach that we used means that the results are not confounded by the fact that different measurement devices were used. Significance for high average rank and heterogeneity was obtained through Monte Carlo testing. Results: For LS-BMD, the quantitative trait locus (QTL) with greatest significance was on chromosome 1p13.3-q23.3 (p ‫ס‬ 0.004), but this exhibited high heterogeneity and the effect was specific for women. Other significant LS-BMD QTLs were on chromosomes 12q24. 31-qter, 3p25.3-p22.1, 11p12-q13.3, and 1q32-q42.3, including one on 18p11-q12.3 that had not been detected by individual studies. For FN-BMD, the strongest QTL was on chromosome 9q31.1-q33.3 (p ‫ס‬ 0.002). Other significant QTLs were identified on chromosomes 17p12-q21.33, 14q13.1-q24.1, 9q21.32-q31.1, and 5q14.3-q23.2. There was no correlation in average ranks of bins between men and women and the loci that regulated BMD in men and women and at different sites were largely distinct.

Bone mineral density (BMD) is a major risk factor for osteoporosis. Circulating monocytes may serve as early progenitors of osteoclasts and produce a wide variety of factors important to bone metabolism. However, little is known about the roles of circulating monocytes in relation to the pathophysiology of osteoporosis. Using the Affymetrix HG-U133A GeneChip® array, we performed a comparative gene expression study of circulating monocytes in subjects with high and low BMD. We identified in total 66 differentially expressed genes including some novel as well as some already known to be relevant to bone metabolism. Three genes potentially contributing to bone metabolism, CCR3 (chemokine receptor 3), HDC (histidine decarboxylase, i.e. the histamine synthesis enzyme), and GCR (glucocorticoid receptor), were confirmed by quantitative real-time reverse transcriptase-PCR as up-regulated in subjects with lower BMD. In addition, significant negative correlation was observed between expression levels of the genes and BMD Z-scores. These three genes and/or their products mediate monocyte chemotaxis, histamine production, and/or sensitivity to glucocorticoids. Our results suggest a novel pathophysiological mechanism for osteoporosis that is characterized by increased recruitment of circulating monocyte into bone, enhanced monocyte differentiation into osteoclasts, as well as osteoclast stimulation via monocyte functional changes. This is the first in vivo microarray study of osteoporosis in humans. The results may contribute to identification of new genes and their functions for osteoporosis and suggest genetic markers to discern individuals at higher risk to osteoporosis with an aim for preventive intervention and treatment.