Pharmacogenomics of osteoporosis: a pathway approach

Pharmacogenetics of Osteoporosis: A Pathway Analysis of the Genetic Influence on the Effects of Antiresorptive Drugs

Osteoporosis is a skeletal disorder defined by a decreased bone mineral density (BMD) and an increased susceptibility to fractures. Bisphosphonates and selective oestrogen receptor modulators (SERM) are among the most widely used drugs. They inhibit bone resorption by targeting the mevalonate and oestrogen pathways, respectively. The aim of this study was to determine if common variants of genes in those pathways influence drug responses. We studied 192 women treated with oral aminobisphosphonates and 51 with SERMs. Genotypes at 154 SNPs of the mevalonate pathway and 806 in the oestrogen pathway were analyzed. Several SNPs located in genes FDPS and FNTA were associated with the bisphosphonate-induced changes in hip bone mineral density (BMD), whereas polymorphisms of the PDSS1, CYP19A1, CYP1A1, and CYP1A2 genes were associated with SERM-induced changes in spine BMD. After multivariate analyses, genotypes combining genes FDPS and FNTA showed a stronger association with bisphosphonate response (r = 0.34; p = 0.00009), whereas the combination of CYP19A1 and PDSS1 genotypes was associated with the response to SERMs (r = 0.62, p = 0.0003). These results suggest that genotyping genes in these pathways may help predict the response to antiresorptive drugs and hence make personalized therapeutic choices.

Parathyroid hormone signaling in mature osteoblasts/osteocytes protects mice from age-related bone loss

Aging is accompanied by osteopenia, characterized by reduced bone formation and increased bone resorption. Osteocytes, the terminally differentiated osteoblasts, are regulators of bone homeostasis, and parathyroid hormone (PTH) receptor (PPR) signaling in mature osteoblasts/osteocytes is essential for PTH-driven anabolic and catabolic skeletal responses. However, the role of PPR signaling in those cells during aging has not been investigated. The aim of this study was to analyze the role of PTH signaling in mature osteoblasts/osteocytes during aging. Mice lacking PPR in osteocyte (Dmp1-PPRKO) display an age-dependent osteopenia characterized by a significant decrease in osteoblast activity and increase in osteoclast number and activity. At the molecular level, the absence of PPR signaling in mature osteoblasts/osteocytes is associated with an increase in serum sclerostin and a significant increase in osteocytes expressing 4-hydroxy-2-nonenals, a marker of oxidative stress. In Dmp1-PPRKO mice there was an age-dependent increase in p16Ink4a/Cdkn2a expression, whereas it was unchanged in controls. In vitro studies demonstrated that PTH protects osteocytes from oxidative stress-induced cell death. In summary, we reported that PPR signaling in osteocytes is important for protecting the skeleton from age-induced bone loss by restraining osteoclast's activity and protecting osteocytes from oxidative stresses.

Postmenopausal Osteoporosis reference genes for qPCR expression assays

Osteoporosis (OP) is a multifactorial disease influenced by genetic factors in more than half of the cases. In spite of the efforts to clarify the relationship among genetic factors and susceptibility to develop OP, many genetic associations need to be further functionally validated. Besides, some limitations as the choice of stably expressed reference genes (RG) should be overcome to ensure the quality and reproducibility of gene expression assays. To our knowledge, a validation study for RG in OP is still missing. We compared the expression levels, using polymerase chain reaction quantitative real time (qPCR) of 10 RG (G6PD, B2M, GUSB, HSP90, EF1A, RPLP0, GAPDH, ACTB, 18 S and HPRT1) to assess their suitability in OP analysis by using GeNorm, Normfinder, BestKeeper and RefFinder programs. A minimal number of two RG was recommended by GeNorm to obtain a reliable normalization. RPLP0 and B2M were identified as the most stable genes in OP studies while ACTB, 18 S and HPRT1 were inadequate for normalization in our data set. Moreover, we showed the dramatic effects of suboptimal RG choice on the quantification of a target gene, highlighting the importance in the identification of the most appropriate reference gene to specific diseases. We suggest the use of RPLP0 and B2M as the most stable reference genes while we do not recommend the use of the least stable reference genes HPRT1, 18 S and ACTB in OP expression assays using PBMC as biological source. Additionally, we emphasize the importance of individualized and careful choice in software and reference genes selection.

Association between polymorphisms in GGPS1 and RANKL genes and postmenopausal osteoporosis in Romanian women

Objectives: This study aimed to assess the relationship between bone mineral density, fragility fractures, fracture risk and polymorphisms of two osteoporosis-candidate genes (GGPS1 and RANKL) in Romanian women with postmenopausal osteoporosis. Methods: An analytical, prospective, transversal, observational, case-control study on 364 postmenopausal women, of which 228 were previously diagnosed with osteoporosis, was carried out between June 2016 and August 2017 in Cluj Napoca, Romania. Clinical data and blood samples were collected from all study participants. Polymorphisms in GGPS1 and RANKL genes were genotyped using TaqMan SNP Genotyping assays, run on a QuantStudio 3 real-time PCR machine. Results: The CT genotype in GGPS1 rs10925503 was associated with significant lower bone mineral density values at lumbar spine and femoral neck sites and a higher fracture risk compared to controls. No significant association was found between genotypes of RANKL rs2277439 with bone mineral density or fracture risk compared to the healthy controls. Conclusions: Our study showed a strong association between low bone mineral density and genotype CT of GGPS1 rs10925503 polymorphisms. No association was found for RANKL rs2277439 polymorphism.

Evaluation of selected traditional Chinese medical extracts for bone mineral density maintenance: A mechanistic study

Objective

To investigate the development of a minimal traditional Chinese medicine (TCM) formula using selected TCM ingredients and evaluating their biological activity with bone-specific in vitro tests. Finally, determining if the minimal formula can maintain bone mineral density (BMD) in a low bone mass (LBM)/osteoporosis (OP) model system.

Methods and results

Sixteen different TCM plant extracts were tested for estrogenic, osteogenic and osteoclastic activities. Despite robust activation of the full-length estrogen receptors α and β by Psoralea corylifolia and Epimedium brevicornu, these extracts do not activate the isolated estrogen ligand binding domains (LBD) of either ERα or ERβ; estrogen (17-β estradiol) fully activates the LBD of ERα and ERβ. E. brevicornu and Drynaria fortunei extracts activated cyclic AMP response elements (CRE) individually and when combined these ingredients stimulated the production of osteoblastic markers Runx2 and Bmp4 in MC3T3-E1 cells. E. brevicornu, Salvia miltiorrhiza, and Astragalus onobrychis extracts inhibited the Il-1β mediated activation of NF-κβ and an E. brevicornu/D. fortunei combination inhibited the development of osteoclasts from precursor cells. Further, a minimal formula containing the E. brevicornu/D. fortunei combination with or without a third ingredient (S. miltiorrhiza, Angelica sinensis, or Lycium barbarum) maintained bone mineral density (BMD) similar to an estradiol-treated control group in the ovariectomized rat; a model LBM/OP system.

Conclusion

A minimal formula consisting of TCM plant extracts that activate CRE and inhibit of NF-κβ activation, but do not behave like estrogen, maintain BMD in a LBM/OP model system.

Polymorphisms in key bone modulator cytokines genes influence bisphosphonates therapy in postmenopausal women

Osteoporosis is a multifactorial and debilitating disease resulting from decreased bone mineral density (BMD) and loss of tissue microarchitecture. Ineffective therapies may lead to bone fractures and subsequent death. Single nucleotide polymorphisms (SNPs) in key immune regulator genes have been associated with therapeutic response to bisphosphonates, which are the first therapeutic line of choice for osteoporosis. However, cytokine pathways and their relation with therapeutic adhesion remain to be fully elucidated. Aimed at better understanding these processes, we investigated the response to bisphosphonate therapy in postmenopausal women and four SNPs in key proinflammatory cytokines genes: IL23R +2284 (C>A) (rs10889677), IL17A +672 (G>A) (rs7747909), IL12B +1188 (T>G) (rs3212227) and INF-γ -1616 (G>A) (rs2069705). A total of 69 patients treated with bisphosphonate were followed for a period of 1 up to 4 years, genotyped and compared according to their changes in bone mineral density (BMD) and level of biochemical markers during their treatment. The INF-γ -1616 G/G associated with increased BMD values in femoral neck (GG/AA, p = 0.016) and decreased BMD values in total hip (GG/GA, p = 0.019; GG/AA, p = 0.011). In relation to biochemical markers, INF-γ -1616 SNP associated with increased alkaline phosphatase (GG/AA; p < 0.0001) and parathyroid hormone levels (AA/GA; p = 0.017). Vitamin D values changes were related to IL17A +672 (GG/GA, p = 0.034) and to IL12B +1188 (TT/TG, p = 0.046) SNPs. Besides, significant differences in changes of calcium levels correlated with IL23R +2284 (CC/CA, p = 0.016) genotypes. Altogether, we suggest that these polymorphisms may play an important role for therapeutic decisions in osteoporosis treatment.

The Regulatory Roles of MicroRNAs in Bone Remodeling and Perspectives as Biomarkers in Osteoporosis

MicroRNAs are involved in many cellular and molecular activities and played important roles in many biological and pathological processes, such as tissue formation, cancer development, diabetes, neurodegenerative diseases, and cardiovascular diseases. Recently, it has been reported that microRNAs can modulate the differentiation and activities of osteoblasts and osteoclasts, the key cells that are involved in bone remodeling process. Meanwhile, the results from our and other research groups showed that the expression profiles of microRNAs in the serum and bone tissues are significantly different in postmenopausal women with or without fractures compared to the control. Therefore, it can be postulated that microRNAs might play important roles in bone remodeling and that they are very likely to be involved in the pathological process of postmenopausal osteoporosis. In this review, we will present the updated research on the regulatory roles of microRNAs in osteoblasts and osteoclasts and the expression profiles of microRNAs in osteoporosis and osteoporotic fracture patients. The perspective of serum microRNAs as novel biomarkers in bone loss disorders such as osteoporosis has also been discussed.

The analysis of DKK1 polymorphisms in relation to skeletal phenotypes and bone response to alendronate treatment in Chinese postmenopausal women

Aim: To investigate the correlation between DKK1 polymorphisms with bone phenotypes and response to alendronate treatment. Materials & methods: Five tag single nucleotide polymorphisms of DKK1 were analyzed in 639 Chinese postmenopausal women with osteoporosis or osteopenia. Bone mineral density (BMD), β-CTX and ALP were measured before and after alendronate treatment. Results: Genotypes at rs1896367, rs1528877 and rs2241529 correlated to baseline BMD (p < 0.05). rs1528877 and rs2241529 polymorphisms correlated to baseline β-CTX levels (p < 0.05). rs2241529 polymorphisms of DKK1 had a small influence on the skeletal response to alendronate treatment (p < 0.05). Conclusion: DKK1 polymorphisms may correlate to baseline BMD and serum β-CTX levels, but present a weak effect on the response to alendronate.

Polymorphisms Associated With Low Bone Mass and High Risk of Atraumatic Fracture

Osteoporosis is a serious disease characterized by high morbidity and mortality due to atraumatic fractures. In the pathogenesis of osteoporosis, except environment and internal factors, such as hormonal imbalance and genetic background, are also in play. In this study candidate genes for osteoporosis were classified according to metabolic or hormonal pathways, which regulate bone mineral density and bone quality (estrogen, RANKL/RANK/OPG axis, mevalonate, the canonical circuit and genes regulating the vitamin D system). COL1A1 and/or COL1A2 genes, which encode formation of the procollagen 1 molecule, were also studied. Mutations in these genes are well-known causes of the inborn disease ‘osteogenesis imperfecta’. In addition to this, polymorphisms in COL1A1 and/or COL1A2 have been found to be associated with parameters of bone quality in adult subjects. The authors discuss the perspectives for the practical utilization of pharmacogenetics (identification of single candidate genes using PCR) and pharmacogenomics (using genome wide association studies (GWAS) to choose optimal treatment for osteoporosis). Potential predictors of antiresorptive therapy efficacy include the following well established genes: ER, FDPS, Cyp19A1, VDR, Col1A1, and Col1A2, as well as the gene for the canonical (Wnt) pathway. Unfortunately, the positive outcomes seen in most association studies have not been confirmed by other researchers. The controversial results could be explained by the use of different methodological approaches in individual studies (different sample size, homogeneity of investigated groups, ethnic differences, or linkage disequilibrium between genes). The key pitfall of association studies is the low variability (7-10 %) of bone phenotypes associated with the investigated genes. Nevertheless, the identification of new genes and the verification of their association with bone density and/or quality (using both PCR and GWAS), remain a great challenge in the optimal prevention and treatment of osteoporosis.

LRP5 polymorphisms and response to alendronate treatment in Chinese postmenopausal women with osteoporosis

AIM

To investigate the association between LRP5 gene polymorphisms and response to alendronate in Chinese osteoporotic women.

MATERIALS & METHODS

Six hundred and thirty nine Chinese postmenopausal women with osteopenia or osteoporosis were included and received alendronate treatment. The A1330V polymorphism of LRP5 was investigated. Bone mineral density (BMD) and bone turnover markers (ALP and β-isomerized carboxy-telopeptide of type I collagen [β-CTX]) were measured before and after treatment. The correlation of LRP5 polymorphisms with changes in BMD and bone turnover biomarkers were analyzed after treatment.

RESULTS

After 12 months of treatment, participants with CC and CT genotypes had a larger increase in lumbar spine BMD and a larger decrease in serum β-CTX and ALP levels than those with TT genotype (all p < 0.001). No significant genotype-treatment interaction was found in hip BMD.

CONCLUSION

The A1330V polymorphism of LRP5 is possibly correlated with response to alendronate treatment in Chinese women with osteoporosis, and the TT genotype could possibly predict a weak response to alendronate.

Pharmacogenomics of osteoporotic fractures

Osteoporosis is a prevalent disease that typically reduces bone strength and predisposes to fractures. It is a multifactorial disorder resulting from the interaction of genetic and acquired factors. Candidate gene studies and, more recently, genome-wide studies have identified a number of polymorphisms significantly associated with bone mass and fractures. Anti-resorptive drugs, which inhibit the differentiation and activity of osteoclasts, are frequently used to treat patients with osteoporosis.Several candidate gene studies have explored the association of genetic factors with drug response, including some common polymorphisms of the gene encoding FDPS (Farnesyl diphosphate synthase), an enzyme that is the main target of aminobisphosphonates. Although scarce data are available, interesting opportunities are open for a better understanding of the pharmacogenetics of osteoporosis and osteoporotic fractures. They include the reanalysis of data already available from epidemiological studies and clinical trials, as well as obtaining pharmacogenetic data in new studies. However, based upon the experience with previous genome-wide association studies, large collaborative efforts would be likely needed to obtain meaningful results.

Parathyroid hormone (PTH)/PTH-related peptide type 1 receptor (PPR) signaling in osteocytes regulates anabolic and catabolic skeletal responses to PTH

Parathyroid hormone (PTH) is the only Food and Drug Administration-approved anabolic agent to treat osteoporosis; however, the cellular targets of PTH action in bone remain controversial. PTH modulates bone turnover by binding to the PTH/PTH-related peptide (PTHrP) type 1 receptor (PPR), a G-protein-coupled receptor highly expressed in bone and kidneys. Osteocytes, the most abundant cells in adult bone, also express PPR. However, the physiological relevance of PPR signaling in osteocytes remains to be elucidated. Toward this goal, we generated mice with PPR deletion in osteocytes (Ocy-PPRKO). Skeletal analysis of these mice revealed a significant increase in bone mineral density and trabecular and cortical bone parameters. Osteoblast activities were reduced in these animals, as demonstrated by decreased collagen type I α1 mRNA and receptor activator of NF-κB ligand (RANKL) expression. Importantly, when subjected to an anabolic or catabolic PTH regimen, Ocy-PPRKO animals demonstrated blunted skeletal responses. PTH failed to suppress SOST/Sclerostin or induce RANKL expression in Ocy-PPRKO animals compared with controls. In vitro, osteoclastogenesis was significantly impaired in Ocy-PPRKO upon PTH administration, indicating that osteocytes control osteoclast formation through a PPR-mediated mechanism. Taken together, these data indicate that PPR signaling in osteocytes is required for bone remodeling, and receptor signaling in osteocytes is needed for anabolic and catabolic skeletal responses.

Genome-wide approaches for identifying genetic risk factors for osteoporosis

Osteoporosis, the most common type of bone disease worldwide, is clinically characterized by low bone mineral density (BMD) and increased susceptibility to fracture. Multiple genetic and environmental factors and gene-environment interactions have been implicated in its pathogenesis. Osteoporosis has strong genetic determination, with the heritability of BMD estimated to be as high as 60%. More than 80 genes or genetic variants have been implicated in risk of osteoporosis by hypothesis-free genome-wide studies. However, these genes or genetic variants can only explain a small portion of BMD variation, suggesting that many other genes or genetic variants underlying osteoporosis risk await discovery. Here, we review recent progress in genome-wide studies of osteoporosis and discuss their implications for medicine and the major challenges in the field.

The future of pharmacogenetics for osteoporosis

The possibility to predict the outcome of medical treatments, both in terms of efficacy and development of adverse effects, is the main goal of modern personalized medicine. The principal aim of pharmacogenetics is to design specific predictive genetic tests, to be performed prior to any drug treatment, and to tailor the therapy for each patient based on the results of these tests. Few pharmacogenetic tests are today validated and commonly applied in clinical practice, and none in the area of osteoporosis and bone disorders. Surely, the complex regulation of bone metabolism and the involvement of numerous different molecular pathways makes it difficult to individuate responsible genes and polymorphisms involved in the modulation of anti-osteoporotic drug response and, subsequently, in designing specific predictive analyses.