Promoter 2 -1025 T/C polymorphism in the RUNX2 gene is associated with femoral neck bmd in Spanish postmenopausal women

Identifying rare variants in genes related to bone phenotypes in a cohort of postmenopausal women

Rare genetic variants in genes previously described to be involved in bone monogenic disorders were identified in postmenopausal women split into two groups according to extreme bone mineral density (BMD) values and lumbar spine Z-scores. A pathogenic variant in COL1A2 gene found in a woman with low BMD highlights the overlap between osteogenesis imperfecta and osteoporosis, which may share their genetic etiology. Other variants were not clearly associated with the extreme BMD, suggesting that there is little contribution of rare variants to postmenopausal osteoporosis.

PURPOSE

We aimed to evaluate whether extreme values of bone mineral density (BMD) in a population-based cohort of postmenopausal women (BARCOS) could be determined by rare genetic variants in genes related to monogenic bone disorders.

METHODS

A panel of 127 genes related to different skeletal phenotypes was designed. Massive sequencing by targeted capture of these genes was performed in 104 DNA samples from those women of the BARCOS cohort that exhibited the highest (HZ group) and lowest (LZ group) LS Z-scores, ranging from + 0.70 to + 3.80 and from - 2.35 to - 4.26, respectively. 5'UTR, 3'UTR, splice region, missense, nonsense, and short indel variants with MAF < 0.01 were annotated with CADD version 1.6 and considered in the analysis.

RESULTS

After filtering those variants with CADD > 25 and present only in one of the groups (either LZ or HZ), six variants were detected, most of which (5/6) were in the LZ group in TCIRG1, COL1A2, SEC24D, LRP4, and ANO5 genes, while only one, in the LMNA gene, was in the HZ group. According to the ClinVar database, the COL1A2 variant, causative of a recessive form of osteogenesis imperfecta, is described as pathogenic, while the other variants are considered of uncertain significance (VUS).

CONCLUSION

The variant identified in COL1A2 in a woman from the LZ group highlights the genetic overlap between monogenic diseases such as osteogenesis imperfecta and complex diseases like osteoporosis. However, the other variants were not clearly associated with the extreme BMD, suggesting that there is little contribution of rare variants to postmenopausal osteoporosis.

Polymorphisms in the Runx2 and osteocalcin genes affect BMD in postmenopausal women: a systematic review and meta-analysis

PURPOSE

Runx2 and osteocalcin have pivotal roles in bone homeostasis. Polymorphism of these two genes could alter the function of osteoblasts and consequently bone mineral density (BMD). Attempts to understand the relationship between these polymorphisms and BMD in postmenopausal women across a variety of populations have yielded inconsistent results. This meta-analysis seeks to define the relationship between these polymorphisms with BMD in postmenopausal women.

METHODS

Eligible studies were identified from three electronic databases. Data were extracted from the eligible studies (4 studies on Runx2 and 6 studies on osteocalcin), and associations of Runx2 T > C and osteocalcin HindIII polymorphisms with BMD in postmenopausal women were assessed using standard difference in means (SDM) and 95% confidence intervals (CI) as statistical measures.

RESULTS

A significant difference in the lumbar spine (LS) BMD in postmenopausal women was observed between the TT and CC homozygotes for the Runx2 T > C (SDM = -0.445, p-value = 0.034). The mutant genotypes (CC) showed significantly lower LS BMD in comparison to wild type genotypes under recessive model of genetic analysis (TC + TT vs. CC: SDM = -0.451, p-value = 0.032). For osteocalcin, HindIII polymorphism, the mutant genotypes (HH) was associated with significantly higher BMD for both LS and femoral neck (FN) than the wild type (hh) homozygotes (SDM = 0.152, p-value = 0.008 and SDM = 0.139, p-value = 0.016 for LS and FN, respectively). There was no association between total hip (TH) BMD and the osteocalcin HindIII polymorphism.

CONCLUSIONS

Runx2 T > C and osteocalcin HindIII polymorphisms influence the level of BMD in postmenopausal women and may be used as predictive markers of osteoporosis.

Genetic analysis in a familial case with high bone mineral density suggests additive effects at two loci

Osteoporosis is the most common bone disease, characterized by a low bone mineral density (BMD) and increased risk of fracture. At the other end of the BMD spectrum, some individuals present strong, fracture‐resistant, bones. Both osteoporosis and high BMD are heritable and their genetic architecture encompasses polygenic inheritance of common variants and some cases of monogenic highly penetrant variants in causal genes. We have investigated the genetics of high BMD in a family segregating this trait in an apparently Mendelian dominant pattern. We searched for rare causal variants by whole‐exome sequencing in three affected and three nonaffected family members. Using this approach, we have identified 38 rare coding variants present in the proband and absent in the three individuals with normal BMD. Although we have found four variants shared by the three affected members of the family, we have not been able to relate any of these to the high‐BMD phenotype. In contrast, we have identified missense variants in two genes, VAV3 and ADGRE5, each shared by two of out of three affected members, whose loss of function fits with the phenotype of the family. In particular, the proband, a woman displaying the highest BMD (sum Z‐score = 7), carries both variants, whereas the other two affected members carry one each. VAV3 encodes a guanine‐nucleotide‐exchange factor with an important role in osteoclast activation and function. Although no previous cases of VAV3 mutations have been reported in humans, Vav3 knockout (KO) mice display dense bones, similarly to the high‐BMD phenotype present in our family. The ADGRE5 gene encodes an adhesion G protein‐coupled receptor expressed in osteoclasts whose KO mouse displays increased trabecular bone volume. Combined, these mouse and human data highlight VAV3 and ADGRE5 as novel putative high‐BMD genes with additive effects, and potential therapeutic targets for osteoporosis. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

The Polymorphism at PLCB4 Promoter (rs6086746) Changes the Binding Affinity of RUNX2 and Affects Osteoporosis Susceptibility: An Analysis of Bioinformatics-Based Case-Control Study and Functional Validation

PURPOSE

Genome-wide association studies have identified numerous genetic variants that are associated with osteoporosis risk; however, most of them are present in the non-coding regions of the genome and the functional mechanisms are unknown. In this study, we aimed to investigate the potential variation in runt domain transcription factor 2 (RUNX2), which is an osteoblast-specific transcription factor that normally stimulates bone formation and osteoblast differentiation, regarding variants within RUNX2 binding sites and risk of osteoporosis in postmenopausal osteoporosis (PMOP).

METHODS

We performed bioinformatics-based prediction by combining whole genome sequencing and chromatin immunoprecipitation sequencing to screen functional SNPs in the RUNX2 binding site using data from the database of Taiwan Biobank; Case-control studies with 651 postmenopausal women comprising 107 osteoporosis patients, 290 osteopenia patients, and 254 controls at Tri-Service General Hospital between 2015 and 2019 were included. The subjects were examined for bone mass density and classified into normal and those with osteopenia or osteoporosis by T-scoring with dual-energy X-ray absorptiometry. Furthermore, mRNA expression and luciferase reporter assay were used to provide additional evidence regarding the associations identified in the association analyses. Chi-square tests and logistic regression were mainly used for statistical assessment.

RESULTS

Through candidate gene approaches, 3 SNPs in the RUNX2 binding site were selected. A novel SNP rs6086746 in the PLCB4 promoter was identified to be associated with osteoporosis in Chinese populations. Patients with AA allele had higher risk of osteoporosis than those with GG+AG (adjusted OR = 6.89; 95% confidence intervals: 2.23-21.31, p = 0.001). Moreover, the AA genotype exhibited lower bone mass density (p < 0.05). Regarding mRNA expression, there were large differences in the correlation between PLCB4 and different RUNX2 alleles (Cohen's q = 0.91). Functionally, the rs6086746 A allele reduces the RUNX2 binding affinity, thus enhancing the suppression of PLCB4 expression (p < 0.05).

CONCLUSIONS

Our results provide further evidence to support the important role of the SNP rs6086746 in the etiology of osteopenia/osteoporosis, thereby enhancing the current understanding of the susceptibility to osteoporosis. We further studied the mechanism underlying osteoporosis regulation by PLCB4.

Functional characterization of the C7ORF76 genomic region, a prominent GWAS signal for osteoporosis in 7q21.3

Genome-wide association studies (GWAS) have repeatedly identified genetic variants associated with bone mineral density (BMD) and osteoporotic fracture in non-coding regions of C7ORF76, a poorly studied gene of unknown function. The aim of the present study was to elucidate the causality and molecular mechanisms underlying the association. We re-sequenced the genomic region in two extreme BMD groups from the BARCOS cohort of postmenopausal women to search for functionally relevant variants. Eight selected variants were tested for association in the complete cohort and 2 of them (rs4342521 and rs10085588) were found significantly associated with lumbar spine BMD and nominally associated with osteoporotic fracture. cis-eQTL analyses of these 2 SNPs, together with SNP rs4727338 (GWAS lead SNP in Estrada et al., Nat Genet. 44:491-501, 2012), performed in human primary osteoblasts, disclosed a statistically significant influence on the expression of the proximal neighbouring gene SLC25A13 and a tendency on the distal SHFM1. We then studied the functionality of a putative upstream regulatory element (UPE), containing rs10085588. Luciferase reporter assays showed transactivation capability with a strong allele-dependent effect. Finally, 4C-seq experiments in osteoblastic cell lines showed that the UPE interacted with different tissue-specific enhancers and a lncRNA (LOC100506136) in the region. In summary, this study is the first one to analyse in depth the functionality of C7ORF76 genomic region. We provide functional regulatory evidence for the rs10085588, which may be a causal SNP within the 7q21.3 GWAS signal for osteoporosis.

Common and rare variants of WNT16, DKK1 and SOST and their relationship with bone mineral density

Numerous GWAS and candidate gene studies have highlighted the role of the Wnt pathway in bone biology. Our objective has been to study in detail the allelic architecture of three Wnt pathway genes: WNT16, DKK1 and SOST, in the context of osteoporosis. We have resequenced the coding and some regulatory regions of these three genes in two groups with extreme bone mineral density (BMD) (n = ∼50, each) from the BARCOS cohort. No interesting novel variants were identified. Thirteen predicted functional variants have been genotyped in the full cohort (n = 1490), and for ten of them (with MAF > 0.01), the association with BMD has been studied. We have found six variants nominally associated with BMD, of which 2 WNT16 variants predicted to be eQTLs for FAM3C (rs55710688, in the Kozak sequence and rs142005327, within a putative enhancer) withstood multiple-testing correction. In addition, two rare variants in functional regions (rs190011371 in WNT16b 3′UTR and rs570754792 in the SOST TATA box) were found only present in three women each, all with BMD below the mean of the cohort. Our results reinforce the higher importance of regulatory versus coding variants in these Wnt pathway genes and open new ways for functional studies of the relevant variants.

β‑Ecdysterone promotes autophagy and inhibits apoptosis in osteoporotic rats

Osteoporosis is an aging process of skeletal tissues with characteristics of reductions in bone mass and microarchitectural deterioration of bone tissue. The present study aimed to investigate the effects of glucocorticoid-induced osteoporosis on osteoblasts and to examine the roles of β-ecdysterone (β-Ecd) involved. In the present study, an in vivo model of osteoporosis was established through the subcutaneous implantation of prednisolone (PRED) into Sprague-Dawley rats, with or without a subcutaneous injection of β-Ecd (5 or 10 mg/kg body weight). Expression of Beclin-1 and microtubule-associated protein 1A/1B-light chain 3I/II and apoptosis in lumbar vertebrae tissues was measured by immunofluorescence and TUNEL assays, respectively. Serum concentration of calcium and phosphorus, and the activity of tartrate-resistant acid phosphatase (TRAP) and alkaline phosphatase (ALP) were measured by biochemical assay. Reverse transcription-quantitative polymerase chain reaction and western blotting was used for detect the expression of related genes and proteins. PRED treatment inhibited bone formation by decreasing bone mineral density, and suppressing the expression of Runt-related transcription factor 2 and bone morphogenetic protein 2, while enhancing the activity of alkaline phosphatase, upregulating the expression of receptor activator of nuclear factor-κB ligand, and increasing the serum content of calcium, phosphorus and tartrate-resistant acid phosphatase in rats. Additionally, PRED was revealed to inhibit autophagy through the downregulation of Beclin-1, autophagy protein 5 and microtubule-associated protein 1A/1B-light chain 3I/II expression, whereas it induced the apoptosis, through the activation of caspase-3 and the suppression of apoptosis regulator BCL2 expression. Notably, the PRED-induced alterations in bone formation, autophagy and apoptosis were revealed to be attenuated by β-Ecd administration. In conclusion, the findings of the present study suggested that β-Ecd may be a promising candidate for the development of therapeutic strategies for the treatment of osteoporosis, through the induction of autophagy and the inhibition of apoptosis in vivo.

Increased Runx2 expression associated with enhanced Wnt signaling in PDLLA internal fixation for fracture treatment

Poly-D-L lactide (PDLLA) biodegradable implants to heal fractures are widely applied in orthopedic surgeries. However, whether the process of fracture healing is regulated differently when PDLLA is used compared with traditional metal materials remains unclear. Runt-related transcription factor 2 (Runx2) and canonical Wnt signaling are essential and may interact reciprocally in the regulation of osteogenesis during bone repair. In the present study, a rat femoral open osteotomy model was used to compare the curative efficacy of a PDLLA rod and Kirschner wire under intramedullary fixation for fracture treatment. The dynamic expression of Runx2 and key components of the canonical Wnt signaling in callus tissue during fracture healing was also investigated. The results of the current study indicate that at weeks 4 and 6 following fixation, the callus bone structural parameters of microCT were significantly improved by PDLLA rod compared to that of Kirschner wire. In addition, at weeks 4 and 6 after fixation, the protein and mRNA expression of Runx2 and the positive regulators of canonical Wnt signaling, such as Wnts and β-catenin, were significantly increased. However, the protein and mRNA expression levels of the negative regulators of canonical Wnt signaling, such as glycogen synthase kinase-3β, were significantly decreased in callus tissue when treated with PDLLA rod compared with Kirschner wire. Collectively, these data indicate that compared to the traditional metal material, using PDLLA internal fixation for fracture treatment may further improve bone formation, which is associated with the increased expression of Runx2 and the enhancement of canonical Wnt signaling.

Osterix and RUNX2 are Transcriptional Regulators of Sclerostin in Human Bone

Sclerostin, encoded by the SOST gene, works as an inhibitor of the Wnt pathway and therefore is an important regulator of bone homeostasis. Due to its potent action as an inhibitor of bone formation, blocking sclerostin activity is the purpose of recently developed anti-osteoporotic treatments. Two bone-specific transcription factors, RUNX2 and OSX, have been shown to interact and co-ordinately regulate the expression of bone-specific genes. Although it has been recently shown that sclerostin is targeted by OSX in mice, there is currently no information of whether this is also the case in human cells. We have identified SP-protein family and AML1 consensus binding sequences at the human SOST promoter and have shown that OSX, together with RUNX2, binds to a specific region close to the transcription start site. Furthermore, we show that OSX and RUNX2 activate SOST expression in a co-ordinated manner in vitro and that SOST expression levels show a significant positive correlation with OSX/RUNX2 expression levels in human bone. We also confirmed previous results showing an association of several SOST/RUNX2 polymorphisms with bone mineral density.

Novel hedgehog agonists promote osteoblast differentiation in mesenchymal stem cells

Hedgehog (Hh) family members are involved in multiple cellular processes including proliferation, migration, differentiation, and cell fate determination. Recently, the novel Hh agonists Hh-Ag 1.3 and 1.7 were identified in a high-throughput screening of small molecule compounds that activate the expression of Gli1, a target of Hh signaling. This study demonstrates that Hh-Ag 1.3 and 1.7 strongly activate the expression of endogenous Gli1 and promote osteoblast differentiation in the mesenchymal stem cell line C3H10T1/2. Both compounds stimulated alkaline phosphatase activity in a dose-dependent manner, and induced osteoblast marker gene expression in C3H10T1/2 cells, which indicated that they had acquired an osteoblast identity. Of the markers, the expression of osterix/Sp7, a downstream target of runt-related transcription factor (Runx)2, was induced by Hh-Ag 1.7, which also rescued the osteoblast differentiation defect of RD-127, a mesenchymal cell line from Runx2-deficient mice. Hh-Ags also activated canonical Wnt signaling and synergized with low doses of BMP-2 to enhance osteoblastic potential. Thus, Hh-Ag 1.7 could be useful for bone healing in individuals with abnormalities in osteogenesis, such as osteoporosis patients and the elderly, and can contribute to the development of novel therapeutics for the treatment of bone fractures and defects.

The association between genetic variants of RUNX2, ADIPOQ and vertebral fracture in Korean postmenopausal women

Contrary to the traditional belief that obesity acts as a protective factor for bone, recent epidemiologic studies have shown that body fat might be a risk factor for osteoporosis and bone fracture. Accordingly, we evaluated the association between the phenotypes of osteoporosis or vertebral fracture and variants of obesity-related genes, peroxisome proliferator-activated receptor-gamma (PPARG), runt-related transcription factor 2 (RUNX2), leptin receptor (LEPR), and adiponectin (ADIPOQ). In total, 907 postmenopausal healthy women, aged 60-79 years, were included in this study. BMD and biomarkers of bone health and adiposity were measured. We genotyped for four single nucleotide polymorphisms (SNPs) from four genes (PPARG, RUNX2, LEPR, ADIPOQ). A general linear model for continuous dependent variables and a logistic regression model for categorical dependent variables were used to analyze the statistical differences among genotype groups. Compared with the TT subjects at rs7771980 in RUNX2, C-carrier (TC + CC) subjects had a lower vertebral fracture risk after adjusting for age, smoking, alcohol, total calorie intake, total energy expenditure, total calcium intake, total fat intake, weight, body fat. Odds ratio (OR) and 95% interval (CI) for the vertebral fracture risk was 0.55 (95% CI 0.32-0.94). After adjusting for multiple variables, the prevalence of vertebral fracture was highest in GG subjects at rs1501299 in ADIPOQ (p = 0.0473). A high calcium intake (>1000 mg/day) contributed to a high bone mineral density (BMD) in GT + TT subjects at rs1501299 in ADIPOQ (p for interaction = 0.0295). Even if the mechanisms between obesity-related genes and bone health are not fully established, the results of our study revealed the association of certain SNPs from obesity-related genes with BMD or vertebral fracture risk in postmenopausal Korean women.

Transcription factor Runx2 in the low bone mineral density of girls with adolescent idiopathic scoliosis

OBJECTIVE

The molecular mechanism of low bone mass in girls with adolescent idiopathic scoliosis (AIS) has not been ascertained. Runx2 is a critical transcription factor regulating osteoblast differentiation and maturation. The present study aimed to explore the possible relationship between Runx2 expression in osteoblasts and bone mineral density (BMD) in subjects with AIS.

METHODS

Twenty-two girls with AIS scheduled to corrective surgery with iliac crest as donor site of autograft for spinal fusion were recruited. The BMD of lumbar spine and femoral neck were assessed by dual-energy X-ray absorptiometry, then patients were divided into two groups with either normal or reduced BMD. Cancellous bone was harvested from their iliac crests for primary culture of osteoblasts. mRNA and protein expression of Runx2 were assayed by reverse transcription-polymerase chain reaction and western blotting, respectively. Results were compared between the two groups and correlated with BMD.

RESULTS

AIS patients with normal BMD showed comparable maturity and body mass index but significant lower Cobb angle of main curve than those of patients with reduced BMD. The mean BMD of lumbar spine and femoral neck were 0.993 g/m(2) and 0.911 g/m(2) in patients with normal BMD, and were 0.757 g/m(2) and 0.733 g/m(2) in those with reduced BMD, respectively. The differences were significant between two groups (P < 0.05). The relative mean mRNA and protein expression of Runx2 were 0.49 ± 0.12 and 0.062 ± 0.020 in AIS with normal BMD, 0.35 ± 0.12 and 0.042 ± 0.006 in AIS with reduced BMD, respectively. Significantly lower Runx2 mRNA and protein expression were found in patients with AIS patients with reduced BMD than in those with normal BMD (P < 0.05). After controlling for age, weight and body mass index, positive correlations were found between Runx2 expression of both mRNA and protein and BMD of lumbar spine and femoral neck.

CONCLUSION

The abnormal expression of Runx2 in patients with AIS and reduced BMD indicates abnormal regulation of differentiation of their osteoblasts. Runx2 may play an important role in the pathogenesis of reduced BMD in patients with AIS.

Genetic analysis of high bone mass cases from the BARCOS cohort of Spanish postmenopausal women

The aims of the study were to establish the prevalence of high bone mass (HBM) in a cohort of Spanish postmenopausal women (BARCOS) and to assess the contribution of LRP5 and DKK1 mutations and of common bone mineral density (BMD) variants to a HBM phenotype. Furthermore, we describe the expression of several osteoblast-specific and Wnt-pathway genes in primary osteoblasts from two HBM cases. A 0.6% of individuals (10/1600) displayed Z-scores in the HBM range (sum Z-score >4). While no mutation in the relevant exons of LRP5 was detected, a rare missense change in DKK1 was found (p.Y74F), which cosegregated with the phenotype in a small pedigree. Fifty-five BMD SNPs from Estrada et al. [NatGenet 44:491-501,2012] were genotyped in the HBM cases to obtain risk scores for each individual. In this small group of samples, Z-scores were found inversely related to risk scores, suggestive of a polygenic etiology. There was a single exception, which may be explained by a rare penetrant genetic variant, counterbalancing the additive effect of the risk alleles. The expression analysis in primary osteoblasts from two HBM cases and five controls suggested that IL6R, DLX3, TWIST1 and PPARG are negatively related to Z-score. One HBM case presented with high levels of RUNX2, while the other displayed very low SOX6. In conclusion, we provide evidence of lack of LRP5 mutations and of a putative HBM-causing mutation in DKK1. Additionally, we present SNP genotyping and expression results that suggest additive effects of several genes for HBM.

Identification of putative target genes of the transcription factor RUNX2

Comparisons of the genomes of Neandertals and Denisovans with present-day human genomes have suggested that the gene RUNX2, which encodes a transcription factor, may have been positively selected during early human evolution. Here, we overexpress RUNX2 in ten human cell lines and identify genes that are directly or indirectly affected by RUNX2 expression. We find a number of genes not previously known to be affected by RUNX2 expression, in particular BIRC3, genes encoded on the mitochondrial genome, and several genes involved in bone and tooth formation. These genes are likely to provide inroads into pathways affected by RUNX2 and potentially by the evolutionary changes that affected RUNX2 in modern humans.

Analyses of RANK and RANKL in the post-GWAS context: functional evidence of vitamin D stimulation through a RANKL distal region

Over the past decade, many genome-wide association studies (GWAs) and meta-analyses have identified genes and regions involved in osteoporotic phenotypes. Nevertheless, the large majority of these results were not tested at any functional level. GWA-associated single-nucleotide polymorphisms (SNPs) near candidate genes such as RANK and RANKL suggest that these SNPs and/or other variants nearby may be involved in bone phenotype determination. This study focuses on SNPs along these two genes, which encode proteins with a well-established role in the bone remodeling equilibrium. Thirty-three SNPs, chosen for their location in evolutionary conserved regions or replicated from previous studies, were genotyped in the BARCOS cohort of 1061 postmenopausal women and tested for association with osteoporotic phenotypes. SNP rs9594738, which lies 184 kb upstream of the RANKL gene, was the only SNP found to be associated with a bone phenotype (dominant model: beta coefficient = -0.034, p = 1.5 × 10(-4) , for lumbar spine bone mineral density). Functional experiments exploring a distal region (DR) of 831 bp that harbors this SNP in a centered position (nt 470) demonstrated its capacity to inhibit the RANKL promoter in reporter gene assays. Remarkably, this DR inhibition was significantly reduced in the presence of vitamin D. In conclusion, the GWA-associated SNP rs9594738 lies in a region involved in transcription regulation through which vitamin D could be regulating RANKL expression and bone mineral density.

Mesenchymal stem cells from osteoporotic patients feature impaired signal transduction but sustained osteoinduction in response to BMP-2 stimulation

Osteoporotic fractures show reduced callus formation and delayed bone healing. Cellular sources of fracture healing are mesenchymal stem cells (MSC) that differentiate into osteoblasts by stimulation with osteoinductive cytokines, such as BMP-2. We hypothesized that impaired signal transduction and reduced osteogenic differentiation capacity in response to BMP-2 may underlie the delayed fracture healing. Therefore, MSC were isolated from femoral heads of healthy and osteoporotic patients. Grouping was carried out by bone mineral densitometry in an age-matched manner. MSC were stimulated with BMP-2. Signal transduction was assessed by western blotting of pSMAD1/5/8 and pERK1/2 as well as by quantitative RT-PCR of Runx-2, Dlx5, and Osteocalcin. Osteogenic differentiation was assessed by quantifying Alizarin Red staining. Osteoporotic MSC featured an accurate phosphorylation pattern of SMAD1/5/8 but a significantly reduced activation of ERK1/2 by BMP-2 stimulation. Furthermore, osteoporotic MSC showed significantly reduced basal expression levels of Runx-2 and Dlx5. However, Runx-2, Dlx5, and Osteocalcin expression showed adequate up-regulation due to BMP-2 stimulation. The global osteogenic differentiation in standard osteogenic differentiation media was reduced in osteoporotic MSC. Nevertheless, osteoporotic MSC were shown to feature an adequate induction of osteogenic differentiation due to BMP-2 stimulation. Taken together, we here demonstrate osteoporosis associated alterations in BMP-2 signaling but sustained specific osteogenic differentiation capacity in response to BMP-2. Therefore, BMP-2 may represent a promising therapeutic agent for the treatment of fractures in osteoporotic patients.

Polyalanine repeat polymorphism in RUNX2 is associated with site-specific fracture in post-menopausal females

Runt related transcription factor 2 (RUNX2) is a key regulator of osteoblast differentiation. Several variations within the RUNX2 gene have been found to be associated with significant changes in BMD, which is a major risk factor for fracture. In this study we report that an 18 bp deletion within the polyalanine tract (17A>11A) of RUNX2 is significantly associated with fracture. Carriers of the 11A allele were found to be nearly twice as likely to have sustained fracture. Within the fracture category, there was a significant tendency of 11A carriers to present with fractures of distal radius and bones of intramembranous origin compared to bones of endochondral origin (p = 0.0001). In a population of random subjects, the 11A allele was associated with decreased levels of serum collagen cross links (CTx, p = 0.01), suggesting decreased bone turnover. The transactivation function of the 11A allele showed a minor quantitative decrease. Interestingly, we found no effect of the 11A allele on BMD at multiple skeletal sites. These findings suggest that the 11A allele is a biologically relevant polymorphism that influences serum CTx and confers enhanced fracture risk in a site-selective manner related to intramembranous bone ossification.

Glutamine repeat variants in human RUNX2 associated with decreased femoral neck BMD, broadband ultrasound attenuation and target gene transactivation

RUNX2 is an essential transcription factor required for skeletal development and cartilage formation. Haploinsufficiency of RUNX2 leads to cleidocranial displaysia (CCD) a skeletal disorder characterised by gross dysgenesis of bones particularly those derived from intramembranous bone formation. A notable feature of the RUNX2 protein is the polyglutamine and polyalanine (23Q/17A) domain coded by a repeat sequence. Since none of the known mutations causing CCD characterised to date map in the glutamine repeat region, we hypothesised that Q-repeat mutations may be related to a more subtle bone phenotype. We screened subjects derived from four normal populations for Q-repeat variants. A total of 22 subjects were identified who were heterozygous for a wild type allele and a Q-repeat variant allele: (15Q, 16Q, 18Q and 30Q). Although not every subject had data for all measures, Q-repeat variants had a significant deficit in BMD with an average decrease of 0.7SD measured over 12 BMD-related parameters (p = 0.005). Femoral neck BMD was measured in all subjects (-0.6SD, p = 0.0007). The transactivation function of RUNX2 was determined for 16Q and 30Q alleles using a reporter gene assay. 16Q and 30Q alleles displayed significantly lower transactivation function compared to wild type (23Q). Our analysis has identified novel Q-repeat mutations that occur at a collective frequency of about 0.4%. These mutations significantly alter BMD and display impaired transactivation function, introducing a new class of functionally relevant RUNX2 mutants.

Common polymorphisms rather than rare genetic variants of the Runx2 gene are associated with femoral neck BMD in Spanish women

RUNX2 is a transcription factor essential for osteoblast differentiation and skeletal morphogenesis. Its mutation creates cleidocranial dysplasia (CCD), a disorder characterized by skeletal abnormalities and bone mineral density (BMD) alterations. The purpose of the present study has been to clarify whether polymorphisms affecting this gene could be associated with changes in BMD in women. To that end, we performed an association study of BMD values from 776 women with two single nucleotide polymorphisms (SNPs) located at P2 promoter (-1025 T>C) and at exon 2 (+198 G>A), and with a deletion polymorphism (17Ala>11Ala), also located at exon 2. We found an association of -1025 T>C SNP with femoral neck BMD (FN-BMD), being the women of TC/CC genotype who have higher BMD than women of TT genotype (P = 0.006). This association was independent of age, weight, menopausal status, or hormone replacement therapy (HRT) use as shown by regression analysis. When women of highest versus lowest quartile of BMD were compared, this association became more evident (P = 0.002), extending also to +198 G>A SNP (GA/AA women with higher FN-BMD; P < 0.05). In addition, we describe herein three novel rare variants in the polyglutamine domain of RUNX2 protein: an in-frame insertion and two deletions in exon 2, resulting in the insertions of 7 and deletions of 7 and 5 glutamines, respectively. These variants do not produce CCD, increased frequency of bone fracture, or BMD alterations. In conclusion, common polymorphisms in Runx2 are associated with FN-BMD. Nevertheless, rare variants that modify the polyglutamine domain of RUNX2 neither have any effect on BMD nor produce the CCD phenotype. These results underscore the significance of polymorphisms in the 5'-region of Runx2 in the determination of FN-BMD.

Analysis of three functional polymorphisms in relation to osteoporosis phenotypes: replication in a Spanish cohort

Osteoporosis is a complex disease involving many putative genetic factors. Association analysis of functional SNPs in candidate genes is an important tool for their identification. However, this approach is affected by limited power, population stratification, and other drawbacks that lead to discordant results. Replication in independent cohorts is essential. We performed association analyses of three functional polymorphisms previously associated with bone phenotypes--namely, Ala222Val in MTHFR, Ile1062Val in LRP6, and -13910C>T in LCT--in a cohort of 944 postmenopausal Spanish women, all of them with lumbar spine (LS) bone mineral density (BMD) data and most with femoral neck (FN) BMD and fracture data. We found significant differences between genotypes only for the MTHFR polymorphism and vertebral factures, with an OR of 2.27 (95% CI 1.17-4.38) for the TT vs. CC/CT genotypes, P = 0.018. We present genotype and allele frequency data for LCT -13910C>T for a Spanish population, where the T allele (conferring lactase persistence) has a frequency of 38.6%. Genotype frequencies were consistent with observed clines in Europe and with the prevalence of lactase nonpersistence. The LCT -13910C>T polymorphism was significantly associated with height and weight, such that T allele carriers were 0.88 cm taller (95% CI 0.08-1.59 cm, P = 0.032, adjusted by age) than CC individuals and TT homozygotes were 1.91 kg heavier than CC/CT individuals (95% CI 0.11-3.71 kg, P = 0.038, adjusted by age). In conclusion, no significant association was observed between the studied polymorphisms and LS BMD or FN BMD in postmenopausal Spanish women, and only MTHFR Ala222Val was associated with vertebral fractures.

Polymorphisms and haplotypes across the osteoprotegerin gene associated with bone mineral density and osteoporotic fractures

Osteoprotegerin plays a key role in bone remodelling. We studied the association between 24 polymorphisms and haplotypes on the OPG gene and bone mineral density and fractures. After multiple-testing correction, one SNP and two block-haplotypes were significantly associated with FN BMD. Two other block-haplotypes were associated with fracture.

INTRODUCTION AND HYPOTHESIS

Osteoprotegerin (OPG) plays a key role in bone remodelling. Here we studied the association between polymorphisms and haplotypes on the OPG gene and bone mineral density (BMD) and fractures.

METHODS

Twenty-four single nucleotide polymorphisms (SNPs) were selected to cover six haplotypic blocks and were genotyped in 964 postmenopausal Spanish women. Haplotypes were established with HaploStats. Association was analysed by GLM (for BMD) and logistic regression (for fractures) both at single SNP and haplotype levels.

RESULTS

Upon adjustment for multiple testing (p < 0.0073), one of the SNPs (SNP #17, rs1032129) remained significantly associated with FN BMD (p = 0.001). Four block-haplotypes stood multiple-testing correction. Two remained associated with FN BMD and two with fracture. The association of block-4 haplotype "AC" (of SNPs #18 and #17) with FN BMD (p = 0.0002) was stronger than that of SNP#17 alone and was the best result overall. A global assessment of the results indicated that all the alleles and haplotypes with a protective effect, at p < 0.05, belonged to a frequent long-range haplotype.

CONCLUSIONS

In conclusion, these results provide a detailed picture of the involvement of common variants and haplotypes of the OPG gene in bone phenotypes.

Genetics of osteoporosis: accelerating pace in gene identification and validation

Osteoporosis is characterized by low bone mineral density and structural deterioration of bone tissue, leading to an increased risk of fractures. It is the most common metabolic bone disorder worldwide, affecting one in three women and one in eight men over the age of 50. In the past 15 years, a large number of genes have been reported as being associated with osteoporosis. However, only in the past 4 years we have witnessed an accelerated pace in identifying and validating osteoporosis susceptibility loci. This increase in pace is mostly due to large-scale association studies, meta-analyses, and genome-wide association studies of both single nucleotide polymorphisms and copy number variations. A comprehensive review of these developments revealed that, to date, at least 15 genes (VDR, ESR1, ESR2, LRP5, LRP4, SOST, GRP177, OPG, RANK, RANKL, COLIA1, SPP1, ITGA1, SP7, and SOX6) can be reasonably assigned as confirmed osteoporosis susceptibility genes, whereas, another >30 genes are promising candidate genes. Notably, confirmed and promising genes are clustered in three biological pathways, the estrogen endocrine pathway, the Wnt/beta-catenin signaling pathway, and the RANKL/RANK/OPG pathway. New biological pathways will certainly emerge when more osteoporosis genes are identified and validated. These genetic findings may provide new routes toward improved therapeutic and preventive interventions of this complex disease.

Dose-dependent effects of Runx2 on bone development

Runx2 controls the commitment of mesenchymal cells to the osteoblastic lineage. Distinct promoters, designated P1 and P2, give rise to functionally similar Runx2-II and Runx2-I isoforms. We postulate that this dual promoter gene structure permits temporal and spatial adjustments in the amount of Runx2 isoforms necessary for optimal bone development. To evaluate the gene dose-dependent effect of Runx2 isoforms on bone development, we intercrossed selective Runx2-II(+/-) with nonselective Runx2-II(+/-)/Runx2-I(+/-) mice to create compound mutant mice: Runx2-II(+/-), Runx2-II(+/-)/Runx2-I(+/-), Runx2-II(-/-), Runx2-II(-/-)/Runx2-I(+/-), Runx2-II(-/-)/Runx2-I(-/-). Analysis of the different Runx2-deficient genotypes showed gene dose-dependent differences in the level of expression of the Runx2 isoforms. In addition, we found that Runx2-I is predominately expressed in the perichondrium and proliferating chondrocytes, whereas Runx2-II is expressed in hypertrophic chondrocytes and metaphyseal osteoblasts. Newborn mice showed impaired development of a mineralized skeleton, bone length, and widening of the hypertrophic zone that were proportionate to the reduction in total Runx2 protein expression. Osteoblast differentiation ex vivo was also proportionate to total amount of Runx2 expression that correlated with reduced Runx2 binding to the osteocalcin promoter by quantitative chromatin immunoprecipitation analysis. Functional analysis of P1 and P2 promoters showed differential regulation of the two promoters in osteoblastic cell lines. These findings support the possibility that the total amount of Runx2 derived from two isoforms and the P1 and P2 promoters, by regulating the time, place, and amount of Runx2 in response to changing environmental cues, impacts on bone development.

Association of a RUNX2 promoter polymorphism with bone mineral density in postmenopausal Korean women

Osteoporosis is characterized by impaired osteoblastogenesis. Bone mineral density (BMD) is a major determinant of bone strength. RUNX2 is an osteoblast-specific transcription factor involved in osteoblast differentiation and ossification. To determine whether RUNX2 is associated with BMD in an ethnically distinct population, we investigated SNPs within the two RUNX2 promoters (P1 and P2) using the Illuminar GoldenGate system in 729 postmenopausal Korean women. Subjects bearing the minor homozygote genotype (CC) at the RUNX2 -1025 T > C SNP (rs7771980) located in P2 showed a significant association with reduced lumbar spine BMD (p = 0.02) and BMDs at proximal femur sites (trochanter, p = 0.05; total femur, p = 0.04) compared with subjects carrying the major homozygote genotype (TT) or the heterozygote genotype (TC), respectively. These results present an interesting genotype association complementary to the previously reported association of BMD with the RUNX2 -1025 T > C P2 SNP in Spanish and Australian cohorts. Therefore, we suggest that the RUNX2 P2 polymorphism (-1025 T > C) may be a useful genetic marker for bone metabolism and may play an important role in BMD in postmenopausal Korean women.

Quantitative trait loci, genes, and polymorphisms that regulate bone mineral density in mouse

This is an in silico analysis of data available from genome-wide scans. Through analysis of QTL, genes and polymorphisms that regulate BMD, we identified 82 BMD QTL, 191 BMD-associated (BMDA) genes, and 83 genes containing known BMD-associated polymorphisms (BMDAP). The catalogue of all BMDA/BMDAP genes and relevant literatures are provided. In total, there are substantially more BMDA/BMDAP genes in regions of the genome where QTL have been identified than in non-QTL regions. Among 191 BMDA genes and 83 BMDAP genes, 133 and 58 are localized in QTL regions, respectively. The difference was still noticeable for the chromosome distribution of these genes between QTL and non-QTL regions. These results have allowed us to generate an integrative profile of QTL, genes, polymorphisms that determine BMD. These data could facilitate more rapid and comprehensive identification of causal genes underlying the determination of BMD in mouse and provide new insights into how BMD is regulated in humans.

A haplotype-based analysis of the LRP5 gene in relation to osteoporosis phenotypes in Spanish postmenopausal women

LRP5 encodes the low-density lipoprotein receptor-related protein 5, a transmembrane protein involved in Wnt signaling. LRP5 is an important regulator of osteoblast growth and differentiation, affecting bone mass in vertebrates. Whether common variations in LRP5 are associated with normal BMD variation or osteoporotic phenotypes is of great relevance. We used a haplotype-based approach to search for common disease-associated variants in LRP5 in a cohort of 964 Spanish postmenopausal women. Twenty-four SNPs were selected, covering the LRP5 region, including the missense changes p.V667M and p.A1330V. The SNPs were genotyped and evaluated for association with BMD at the lumbar spine (LS) or femoral neck (FN) and with osteoporotic fracture, at single SNP and haplotype levels, by regression methods. Association with LS BMD was found for SNP 1, rs312009, located in the 5'-flanking region (p = 0.011, recessive model). SNP 6, rs2508836, in intron 1, was also associated with BMD, both at LS (p = 0.025, additive model) and FN (p = 0.031, recessive model). Two polymorphisms were associated with fracture: SNP 11, rs729635, in intron 1, and SNP 15, rs643892, in intron 5 (p = 0.007 additive model and p = 0.019 recessive model, respectively). Haplotype analyses did not provide additional information, except for haplotype "GC" of the block located at the 3'end of the gene. This haplotype spans intron 22 and the 3' untranslated region and was associated with FN BMD (p = 0.029, one copy of the haplotype versus none). In silico analyses showed that SNP 1 (rs312009) lies in a putative RUNX2 binding site. Electro-mobility shift assays confirmed RUNX2 binding to this site.

Molecular genetic studies of gene identification for osteoporosis

This review comprehensively summarizes the most important and representative molecular genetics studies of gene identification for osteoporosis published up to the end of September 2007. It is intended to constitute a sequential update of our previously published reviews covering the available data up to the end of 2004. Evidence from candidate gene-association studies, genome-wide linkage and association studies, as well as functional genomic studies (including gene-expression microarray and proteomics) on osteogenesis and osteoporosis, are reviewed separately. Studies of transgenic and knockout mice models relevant to osteoporosis are summarized. The major results of all studies are tabulated for comparison and ease of reference. Comments are made on the most notable findings and representative studies for their potential influence and implications on our present understanding of genetics of osteoporosis. The format adopted by this review should be ideal for accommodating future new advances and studies.