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

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.

On the association between Chiari malformation type 1, bone mineral density and bone related genes

Background

Chiari malformation type 1 (C1M) is a neurological disease characterized by herniation of the cerebellar tonsils below the foramen magnum. Cranial bone constriction is suspected to be its main cause. To date, genes related to bone development (e.g. DKK1 or COL1A2) have been associated with C1M, while some bone diseases (e.g. Paget) have been found to cosegregate with C1M. Nevertheless, the association between bone mineral density (BMD) and C1M has not been investigated, yet. Here, we systematically investigate the association between C1M and BMD, and between bone related genes and C1M.

Methods

We have recruited a small cohort of C1M patients (12 unrelated patients) in whom we have performed targeted sequencing of an in-house bone-related gene panel and BMD determination through non-invasive DXA.

Results

In the search for association between the bone related genes and C1M we have found variants in more than one C1M patient in WNT16, CRTAP, MYO7A and NOTCH2. These genes have been either associated with craniofacial development in different ways, or previously associated with C1M (MYO7A). Regarding the potential link between BMD and C1M, we have found three osteoporotic patients and one patient who had high BMD, very close to the HBM phenotype values, although most patients had normal BMD.

Conclusions

Variants in bone related genes have been repeatedly found in some C1M cases. The relationship of bone genes with C1M deserves further study, to get a clearer estimate of their contribution to its etiology. No direct correlation between BMD and C1M was observed.

•

We used an in-house bone gene panel to investigate a small cohort of C1M patients.

•

Variants in WNT16, CRTAP, MYO7A and NOTCH2 were found in more than one C1M patient.

•

No clear relationship was found between C1M and BMD in this small C1M cohort.

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.

Wnt Pathway Extracellular Components and Their Essential Roles in Bone Homeostasis

The Wnt pathway is involved in several processes essential for bone development and homeostasis. For proper functioning, the Wnt pathway is tightly regulated by numerous extracellular elements that act by both activating and inhibiting the pathway at different moments. This review aims to describe, summarize and update the findings regarding the extracellular modulators of the Wnt pathway, including co-receptors, ligands and inhibitors, in relation to bone homeostasis, with an emphasis on the animal models generated, the diseases associated with each gene and the bone processes in which each member is involved. The precise knowledge of all these elements will help us to identify possible targets that can be used as a therapeutic target for the treatment of bone diseases such as osteoporosis.

Genetics and Genomics of SOST: Functional Analysis of Variants and Genomic Regulation in Osteoblasts

SOST encodes the sclerostin protein, which acts as a key extracellular inhibitor of the canonical Wnt pathway in bone, playing a crucial role in skeletal development and bone homeostasis. The objective of this work was to assess the functionality of two variants previously identified (the rare variant rs570754792 and the missense variant p.Val10Ile) and to investigate the physical interactors of the SOST proximal promoter region in bone cells. Through a promoter luciferase reporter assay we show that the minor allele of rs570754792, a variant located in the extended TATA box motif, displays a significant decrease in promoter activity. Likewise, through western blot studies of extracellular and intracellular sclerostin, we observe a reduced expression of the p.Val10Ile mutant protein. Finally, using a circular chromosome conformation capture assay (4C-seq) in 3 bone cell types (MSC, hFOB, Saos-2), we have detected physical interactions between the SOST proximal promoter and the ECR5 enhancer, several additional enhancers located between EVT4 and MEOX1 and a distant region containing exon 18 of DHX8. In conclusion, SOST presents functional regulatory and missense variants that affect its expression and displays physical contacts with far reaching genomic sequences, which may play a role in its regulation within bone cells.

Functional Assessment of Coding and Regulatory Variants From the DKK1 Locus

The DKK1 gene encodes an extracellular inhibitor of the Wnt pathway with an important role in bone tissue development, bone homeostasis, and different critical aspects of bone biology. Several BMD genome‐wide association studies (GWASs) have consistently found association with SNPs in the DKK1 genomic region. For these reasons, it is important to assess the functionality of coding and regulatory variants in the gene. Here, we have studied the functionality of putative regulatory variants, previously found associated with BMD in different studies by others and ourselves, and also six missense variants present in the general population. Using a Wnt‐pathway‐specific luciferase reporter assay, we have determined that the variants p.Ala41Thr, p.Tyr74Phe, p.Arg120Leu, and p.Ser157Ile display a reduced DKK1 inhibitory capacity as compared with WT. This result agrees with the high‐bone‐mass (HBM) phenotype of two women from our cohort who carried mutations p.Tyr74Phe or p.Arg120Leu. On the other hand, by means of a circularized chromosome conformation capture‐ (4C‐) sequencing experiment, we have detected that the region containing 24 BMD‐GWA variants, located 350‐kb downstream of DKK1, interacts both with DKK1 and the LNCAROD (LncRNA‐activating regulator of DKK1, AKA LINC0148) in osteoblastic cells. In conclusion, we have shown that some rare coding variants are partial loss‐of‐function mutations that may lead to a HBM phenotype, whereas the common SNPs associated with BMD in GWASs belong to a putative long‐range regulatory region, through a yet unknown mechanism involving LNCAROD. © 2020 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Alterations in DNA methylation profiles in cancellous bone of postmenopausal women with osteoporosis

Osteoporosis is characterized by systemic microarchitecture impairment and bone loss, which ultimately lead to fragility fractures. This disease is most common in older people, especially in postmenopausal women. Cancellous bone is affected by osteoporosis earlier than cortical bone, and DNA methylation microarray analysis of the hip cancellous bone of patients with osteoarthritis revealed differential methylation. In view of the important role of cancellous bone in bone development, we examined genome‐wide DNA methylation profiles in the cancellous bone from patients with postmenopausal osteoporosis versus healthy postmenopausal women using Illumina 850K methylation microarray analysis. Under a threshold of P < 0.05, we obtained a total of 8973 differentially methylated genes, such as SOX6, ACE, SYK and TGFB3. Under a threshold of P < 0.05 and |△β| > 0.2, a total of 17 and 34 key differentially methylated genes were further identified at the promoter region and cytosine‐ phosphate‐ guanine (CpG) islands (such as PRKCZ, GNA11 and COL4A1), respectively. PLEKHA2, PLEKHB1, PNPLA7, SCD, MGST3 and TSNAX were the most common differentially methylated genes at both the promoter region and CpG islands. Five important signaling pathways, including the calcium signaling pathway, the cyclic guanosine phospho‐protein kinase G (cGMP‐PKG) signaling pathway, endocytosis, the Rap1 signaling pathway and the AMPK signaling pathway were identified. Our study may be suitable as a basis for exploring the mechanisms underlying osteoporosis in postmenopausal women.

Twist1 Inactivation in Dmp1-Expressing Cells Increases Bone Mass but Does Not Affect the Anabolic Response to Sclerostin Neutralization

Wnt signaling plays a major role in bone metabolism. Advances in our understanding of secreted regulators of Wnt have yielded several therapeutic targets to stimulate osteoanabolism—the most promising of which is the Wnt inhibitor sclerostin. Sclerostin antibody recently gained approval for clinical use to treat osteoporosis, but the biology surrounding sclerostin antagonism is still incompletely understood. Numerous factors regulate the efficacy of sclerostin inhibition on bone formation, a process known as self-regulation. In previous communications we reported that the basic helix-loop-helix transcription factor Twist1—a gene know to regulate skeletal development—is highly upregulated among the osteocyte cell population in mice treated with sclerostin antibody. In this communication, we tested the hypothesis that preventing Twist1 upregulation by deletion of Twist1 from late-stage osteoblasts and osteocytes would increase the efficacy of sclerostin antibody treatment, since Twist1 is known to restrain osteoblast activity in many models. Twist1-floxed loss-of-function mice were crossed to the Dmp1-Cre driver to delete Twist1 in Dmp1-expressing cells. Conditional Twist1 deletion was associated with a mild but significant increase in bone mass, as assessed by dual energy x-ray absorptiometry (DXA) and microCT (µCT) for many endpoints in both male and female mice. Biomechanical properties of the femur were not affected by conditional mutation of Twist1. Sclerostin antibody improved all bone properties significantly, regardless of Twist1 status, sex, or endpoint examined. No interactions were detected when Twist1 status and antibody treatment were examined together, suggesting that Twist1 upregulation in the osteocyte population is not an endogenous mechanism that restrains the osteoanabolic effect of sclerostin antibody treatment. In summary, Twist1 inhibition in the late-stage osteoblast/osteocyte increases bone mass but does not affect the anabolic response to sclerostin neutralization.

Cortical and Trabecular Bone Analysis of Patients With High Bone Mass From the Barcelona Osteoporosis Cohort Using 3-Dimensional Dual-Energy X-ray Absorptiometry: A Case-Control Study

High bone mass (HBM), a rare phenotype, can be detected by dual-energy X-ray absorptiometry (DXA) scanning. Measurements with peripheral quantitative computed tomography at the tibia have found increased trabecular bone mineral density and changes in cortical bone density and structure, all of which lead to increased bone strength. However, no studies on cortical and trabecular bone have been performed at the femur. The recently developed 3-dimensional (3D)-DXA software algorithm quantifies the trabecular and cortical volumetric bone mineral density (vBMD) and the anatomical distribution of cortical thickness using routine hip DXA scans. We analyzed the femurs of 15 women with HBM and 15 controls from the Barcelona Osteoporosis (BARCOS) cohort using the 3D-DXA technique. The mean vBMD of proximal femur was 29.7% higher in HBM cases than in controls for the integral bone, 41.3% higher for the trabecular bone, and 7.3% higher for the cortical bone (p < 0.001). No differences in bone size were detected between cases and controls. Patients with HBM had a thicker cortex and higher trabecular and cortical vBMDs, as measured by 3D-DXA at the femur and compared to controls; bone size was similar in both groups. To the best of our knowledge, this is the first description of trabecular and cortical characteristics of the hip in patients with HBM.

Distinct DNA methylation profiles in bone and blood of osteoporotic and healthy postmenopausal women

DNA methylation affects expression of associated genes and may contribute to the missing genetic effects from genome-wide association studies of osteoporosis. To improve insight into the mechanisms of postmenopausal osteoporosis, we combined transcript profiling with DNA methylation analyses in bone. RNA and DNA were isolated from 84 bone biopsies of postmenopausal donors varying markedly in bone mineral density (BMD). In all, 2529 CpGs in the top 100 genes most significantly associated with BMD were analyzed. The methylation levels at 63 CpGs differed significantly between healthy and osteoporotic women at 10% false discovery rate (FDR). Five of these CpGs at 5% FDR could explain 14% of BMD variation. To test whether blood DNA methylation reflect the situation in bone (as shown for other tissues), an independent cohort was selected and BMD association was demonstrated in blood for 13 of the 63 CpGs. Four transcripts representing inhibitors of bone metabolism-MEPE, SOST, WIF1, and DKK1-showed correlation to a high number of methylated CpGs, at 5% FDR. Our results link DNA methylation to the genetic influence modifying the skeleton, and the data suggest a complex interaction between CpG methylation and gene regulation. This is the first study in the hitherto largest number of postmenopausal women to demonstrate a strong association among bone CpG methylation, transcript levels, and BMD/fracture. This new insight may have implications for evaluation of osteoporosis stage and susceptibility.

SNPs in bone-related miRNAs are associated with the osteoporotic phenotype

Biogenesis and function of microRNAs can be influenced by genetic variants in the pri-miRNA sequences leading to phenotypic variability. This study aims to identify single nucleotide polymorphisms (SNPs) affecting the expression levels of bone-related mature microRNAs and thus, triggering an osteoporotic phenotype. An association analysis of SNPs located in pri-miRNA sequences with bone mineral density (BMD) was performed in the OSTEOMED2 cohort (n = 2183). Functional studies were performed for assessing the role of BMD-associated miRNAs in bone cells. Two SNPs, rs6430498 in the miR-3679 and rs12512664 in the miR-4274, were significantly associated with femoral neck BMD. Further, we measured these BMD-associated microRNAs in trabecular bone from osteoporotic hip fractures comparing to non-osteoporotic bone by qPCR. Both microRNAs were found overexpressed in fractured bone. Increased matrix mineralization was observed after miR-3679-3p inhibition in human osteoblastic cells. Finally, genotypes of rs6430498 and rs12512664 were correlated with expression levels of miR-3679 and miR-4274, respectively, in osteoblasts. In both cases, the allele that generated higher microRNA expression levels was associated with lower BMD values. In conclusion, two osteoblast-expressed microRNAs, miR-3679 and miR-4274, were associated with BMD; their overexpression could contribute to the osteoporotic phenotype. These findings open new areas for the study of bone disorders.

Mutations in Known Monogenic High Bone Mass Loci Only Explain a Small Proportion of High Bone Mass Cases

High bone mass (HBM) can be an incidental clinical finding; however, monogenic HBM disorders (eg, LRP5 or SOST mutations) are rare. We aimed to determine to what extent HBM is explained by mutations in known HBM genes. A total of 258 unrelated HBM cases were identified from a review of 335,115 DXA scans from 13 UK centers. Cases were assessed clinically and underwent sequencing of known anabolic HBM loci: LRP5 (exons 2, 3, 4), LRP4 (exons 25, 26), SOST (exons 1, 2, and the van Buchem's disease [VBD] 52-kb intronic deletion 3'). Family members were assessed for HBM segregation with identified variants. Three-dimensional protein models were constructed for identified variants. Two novel missense LRP5 HBM mutations ([c.518C>T; p.Thr173Met], [c.796C>T; p.Arg266Cys]) were identified, plus three previously reported missense LRP5 mutations ([c.593A>G; p.Asn198Ser], [c.724G>A; p.Ala242Thr], [c.266A>G; p.Gln89Arg]), associated with HBM in 11 adults from seven families. Individuals with LRP5 HBM (∼prevalence 5/100,000) displayed a variable phenotype of skeletal dysplasia with increased trabecular BMD and cortical thickness on HRpQCT, and gynoid fat mass accumulation on DXA, compared with both non-LRP5 HBM and controls. One mostly asymptomatic woman carried a novel heterozygous nonsense SOST mutation (c.530C>A; p.Ser177X) predicted to prematurely truncate sclerostin. Protein modeling suggests the severity of the LRP5-HBM phenotype corresponds to the degree of protein disruption and the consequent effect on SOST-LRP5 binding. We predict p.Asn198Ser and p.Ala242Thr directly disrupt SOST binding; both correspond to severe HBM phenotypes (BMD Z-scores +3.1 to +12.2, inability to float). Less disruptive structural alterations predicted from p.Arg266Cys, p.Thr173Met, and p.Gln89Arg were associated with less severe phenotypes (Z-scores +2.4 to +6.2, ability to float). In conclusion, although mutations in known HBM loci may be asymptomatic, they only account for a very small proportion (∼3%) of HBM individuals, suggesting the great majority are explained by either unknown monogenic causes or polygenic inheritance.

Unexplained high BMD in DXA-scanned patients is generalized throughout the skeleton and characterized by thicker cortical and trabecular bone

Unexplained high bone mineral density (BMD) is a rare condition and the mechanisms responsible are yet to be described in detail. The aim of the study was to identify patients with unexplained high BMD from a local DXA database and compare their radiological phenotype with an age- and a gender-matched group of population-based controls. We defined high BMD as a DXA Z-score ≥ + 2.5 at the total hip and lumbar spine. We characterized the findings as "unexplained" if no osteodegenerative changes, bone metabolic disease, or arthritis at the hip or lumbar spine was observed. All participants were investigated with high-resolution peripheral quantitative computed tomography (HR-pQCT), QCT, DXA, fasting blood samples, a 24-h urine sample, and questionnaires. The DXA database contained data on 25,118 patients. Initially, 138 (0.55%) potential participants with high BMD were identified, and during the study ten additional cases were identified from new DXA scans. Sixty-seven patients accepted to participate in the study, and among these we identified 15 women and one man with unexplained high BMD. These 15 women had higher BMD throughout the skeleton relative to controls, similar area/volume at the hip and the distal extremities, a higher number of trabeculae, which was thicker than in the controls, and a higher finite element estimated bone strength. The 15 women were heavier and had a higher fat mass then controls. We conclude that patients with unexplained high BMD have a generalized high BMD phenotype throughout their skeleton, which is characterized with a denser microarchitecture.