Identification of novel RANK polymorphisms and their putative association with low BMD among postmenopausal women

Genetic polymorphisms and their influence on therapeutic response to alendronate-a pilot study

Introduction: Osteoporosis has a strong genetic contribution, and several genes have been shown to influence bone mineral density. Variants in the human genome are considered important causes of differences in drug responses observed in clinical practice. In terms of bone mineral density, about 26–53% of patients do not respond to amino-bisphosphonate therapies, of which alendronate is the most widely used. Material and method: The current study is prospective, observational, analytical, longitudinal and cohort type. It included 25 postmenopausal women treated with alendronate for 1 year. Bone mineral density at lumbar spine and proximal femur was measured and bone turnover markers (C-terminal telopeptide of type I collagen and procollagen 1N-terminal propeptide) were evaluated at 0 and 12 months of treatment. Six single nucleotide polymorphisms in osteoporosis-candidate genes were genotyped (FDPS rs2297480, LRP5 rs3736228, SOST rs1234612, VKORC1 rs9934438, GGPS1 rs10925503 and RANKL rs2277439). Treatment response was evaluated by percentage changes in bone mineral density and bone turnover markers. Results: The heterozygous CT of FDPS rs2297480 showed lower increases in BMD values in the lumbar spine region and the homozygous CC of the GGPS1 rs10925503 showed lower increases in terms of BMD at the total hip region. No association was found for LRP5 rs3736228, SOST rs1234612, VKORC1 rs9934438 and RANKL rs2277439. Conclusions: Romanian postmenopausal women with osteoporosis carrying the CT genotype of FDPS rs2297480 or the CC genotype of GGPS1 rs10925503 could have an unsatisfactory response to alendronate treatment. Key words: osteoporosis; genetic polymorphism; alendronate; bone mineral density; bone turnover markers,

RANKL/RANK/OPG Polymorphisms and Heel Quantitative Ultrasound in Young Adults

BACKGROUND

The receptor activator of the nuclear factor-kappa B ligand (RANKL), the receptor activator of nuclear factor-kappa B (RANK), and the osteoprotegerin (OPG) signaling pathway play an important role in the regulation of bone remodeling and osteoclast differentiation. Quantitative ultrasound (QUS) is a relatively recent and noninvasive method providing structural information on microstructure, bone elasticity, and connectivity. However, in contrast to bone mineral density measurements, the possible association of the RANKL/RANK/OPG pathway with heel QUS has not been analyzed.

OBJECTIVES

The aim of this study was to assess, for the first time, the contribution of the RANKL/RANK/OPG pathway genes in the genetic background of heel QUS parameters.

METHODS

Ten single-nucleotide polymorphisms (SNPs) of RANKL (rs9594759, rs12585014, rs7988338, rs2148073), RANK (rs1805034, rs12458117, rs3018362), and OPG (rs4355801, rs3102735, rs2073618) were selected as genetic markers and genotyped using Open Array technology in 575 self-reported Caucasian individuals aged 18-25. Bone mass in the right calcaneus was estimated with QUS to obtain the broadband ultrasound attenuation (BUA) measurement (dB/MHz). Linear regression analyses were performed to test the possible association between the SNPs and BUA.

RESULTS

Linear regression analysis of all the tested SNPs revealed no significant association with the BUA parameter after adjusting for age, gender, weight, height, physical activity, and calcium intake. The lowest p-value was observed for the rs9594759 RANKL polymorphism and heel QUS (p = .06; b* = -.075, 95% CI [-0.960, 0.028]).

CONCLUSION

Our results suggest that the polymorphism of the RANKL, RANK, and OPG genes does not make a significant genetic contribution to heel ultrasound measurements in a population of young Caucasian adults. Further studies replicating the results in independent populations are needed to support these initial findings.

Polymorphisms in genes in the RANKL/RANK/OPG pathway are associated with bone mineral density at different skeletal sites in post-menopausal women

Association between 22 single nucleotide polymorphisms (SNPs) in the TNFSF11, TNFRSF11A, and TNFRSF11B genes in the RANKL/RANK/OPG pathway with bone mineral density (BMD) in 881 post-menopausal women. Our results suggest that TNFSF11 and TNFRSF11A, but not TNFRSF11B, genetic polymorphisms influence BMD mainly in the femoral neck in post-menopausal Chinese women.

INTRODUCTION

The aim of this study was to assess the relationship of polymorphisms in the TNFSF11, TNFRSF11A, and TNFRSF11B genes in the RANKL/RANK/OPG pathway with bone mineral density (BMD) in a cohort of Chinese post-menopausal women.

METHODS

A cross-sectional study was conducted in 881 post-menopausal women aged 50-89 years. All participants underwent lumbar spinal (LS) and femoral neck (FN) BMD evaluation by dual-energy X-ray absorptiometry. Twenty-two TNFSF11, TNFRSF11A, and TNFRSF11B SNPs were genotyped. We tested whether a single SNP or a haplotype was associated with BMD variations.

RESULTS

Two SNPs in the TNFSF11 gene (rs2277439 and rs2324851) and one in the TNFRSF11A gene (rs7239261) were found to be significantly associated with FN BMD (p = 0.014, 0.013, and 0.047, respectively). Haplotype TGACGT of TNFSF11 rs9525641-rs2277439-rs2324851-rs2875459-rs2200287-rs9533166 was a genetic risk factor toward a lower FN BMD (beta = -0.1473; p = 0.01126). In contrary, haplotype TAGCGT of TNFSF11 rs9525641-rs2277439-rs2324851-rs2875459-rs2200287-rs9533166 was genetic protective factor for LS BMD (beta = 0.3923; p = 0.04917).

CONCLUSIONS

Our findings suggest that TNFSF11 and TNFRSF11A, but not TNFRSF11B, genetic polymorphisms influence BMD mainly in the femoral neck in post-menopausal Chinese women. This contributes to the understanding of the role of genetic variation in this pathway in determining bone health.

Multiple gene polymorphisms can improve prediction of nonvertebral fracture in postmenopausal women

Clinical risk factors (CRFs), with or without bone mineral density (BMD), are used to determine the risk of osteoporotic fracture (OF), which has a heritable component. In this study we investigated whether genetic profiling can additionally improve the ability to predict OF. Using 1229 unrelated Korean postmenopausal women, 39 single-nucleotide polymorphisms (SNPs) in 30 human genomic loci were tested for association with osteoporosis-related traits, such as BMD, osteoporosis, vertebral fracture (VF), nonvertebral fracture (NVF), and any fracture. To estimate the effects of genetic profiling, the genetic risk score (GRS) was calculated using five prediction models: (Model I) GRSs only; (Model II) BMD only; (Model III) CRFs only; (Model IV) CRFs and BMD; and (Model V) CRFs, BMD, and GRS. A total of 21 SNPs within 19 genes associated with one or more osteoporosis-related traits and were included for GRS calculation. GRS associated with BMD before and after adjustment for CRFs (p ranging from <0.001 to 0.018). GRS associated with NVF before and after adjustment for CRFs and BMD (p ranging from 0.017 to 0.045), and with any fracture after adjustment for CRFs and femur neck BMD (p = 0.049). In terms of predicting NVF, the area under the receiver operating characteristic curve (AUC) for Model I was 0.55, which was lower than the AUCs of Models II (0.60), III (0.64), and IV (0.65). Adding GRS to Model IV (in Model V) increased the AUC to 0.67, and improved the accuracy of NVF classification by 11.5% (p = 0.014). In terms of predicting any fracture, the AUC of Model V (0.68) was similar to that of Model IV (0.68), and Model V did not significantly improve the accuracy of any fracture classification (p = 0.39). Thus, genetic profiling may enhance the accuracy of NVF predictions and help to delineate the intervention threshold.

Association of genetic polymorphisms of RANK, RANKL and OPG with bone mineral density in Chinese peri- and postmenopausal women

OBJECTIVES

To explore the influence of 14 single nucleotide polymorphisms (SNPs) in receptor activator of nuclear factor-kappa B (RANK), RANK ligand (RANKL) and osteoprotegerin (OPG) on bone mineral density (BMD) in a Chinese female population.

DESIGN AND METHODS

A cross-sectional study was conducted in 108 perimenopausal and 127 postmenopausal women aged 43-65 years. All participants underwent lumbar spinal and nondominant femoral BMD evaluation by dual energy X-ray absorptiometry. Fourteen RANK, RANKL and OPG genotypes were determined by chip-based MALDI-TOF mass spectrometry. The differences between the BMDs of the RANK genotypes were analyzed.

RESULTS

Five SNPs (rs6993813, rs4355801, rs1032129 and rs2073618 in OPG and rs3018362 in RANK) were significantly associated with BMD or with BMD adjusted for body weight or years since menopause, mostly at the femoral neck but also partly at the total hip (p<0.05). The risk allele frequencies observed in our sample were different from those found in Europeans but the effects of these risk alleles on BMD values had the same direction in our cohort as in Europeans, except for rs3018362 with G as the risk allele, which was contrary to other studies. None of the SNPs in RANKL were associated with BMD at any anatomical site.

CONCLUSIONS

Our findings suggest that OPG and RANK but not RANKL genetic polymorphisms influence BMD mainly in the femoral neck in peri- and postmenopausal Chinese women. This contributes to the understanding of the role of genetic variation in this pathway in determining bone health.

Genetic variation in the TNFRSF11A gene encoding RANK is associated with susceptibility to Paget's disease of bone

RANK (receptor activator of nuclear factor-κB), encoded by TNFRSF11A, is a key protein in osteoclastogenesis. TNFRSF11A mutations cause Paget's disease of bone (PDB)-like diseases (ie, familial expansile osteolysis, expansile skeletal hyperphosphatasia, and early-onset PDB) and an osteoclast-poor form of osteopetrosis. However, no TNFRSF11A mutations have been found in classic PDB, neither in familial nor in isolated cases. To investigate the possible relationship between TNFRSF11A polymorphisms and sporadic PDB, we conducted an association study including 32 single-nucleotide polymorphisms (SNPs) in 196 Belgian sporadic PDB patients and 212 control individuals. Thirteen SNPs and 3 multimarker tests (MMTs) turned out to have a p value of between .036 and 3.17 × 10(-4) , with the major effect coming from females. Moreover, 6 SNPs and 1 MMT withstood the Bonferroni correction (p < .002). Replication studies were performed for 2 nonsynonymous SNPs (rs35211496 and rs1805034) in a Dutch and a British cohort. Interestingly, both SNPs resulted in p values ranging from .013 to 8.38 × 10(-5) in both populations. Meta-analysis over three populations resulted in p = .002 for rs35211496 and p = 1.27 × 10(-8) for rs1805034, again mainly coming from the female subgroups. In an attempt to identify the underlying causative SNP, we performed functional studies for the coding SNPs as well as resequencing efforts of a 31-kb region harboring a risk haplotype within the Belgian females. However, neither approach resulted in significant evidence for the causality of any of the tested genetic variants. Therefore, further studies are needed to identify the real cause of the increased risk to develop PDB shown to be present within TNFRSF11A.

Genetics of osteoporosis

Osteoporosis is a common disease with a strong genetic component characterized by reduced bone mass, defects in the microarchitecture of bone tissue, and an increased risk of fragility fractures. Twin and family studies have shown high heritability of bone mineral density (BMD) and other determinants of fracture risk such as ultrasound properties of bone, skeletal geometry, and bone turnover. Osteoporotic fractures also have a heritable component, but this reduces with age as environmental factors such as risk of falling come into play. Susceptibility to osteoporosis is governed by many different genetic variants and their interaction with environmental factors such as diet and exercise. Notable successes in identification of genes that regulate BMD have come from the study of rare Mendelian bone diseases characterized by major abnormalities of bone mass where variants of large effect size are operative. Genome-wide association studies have also identified common genetic variants of small effect size that contribute to regulation of BMD and fracture risk in the general population. In many cases, the loci and genes identified by these studies had not previously been suspected to play a role in bone metabolism. Although there has been extensive progress in identifying the genes and loci that contribute to the regulation of BMD and fracture over the past 15 yr, most of the genetic variants that regulate these phenotypes remain to be discovered.

Genetic variation in the RANKL/RANK/OPG signaling pathway is associated with bone turnover and bone mineral density in men

The aim of this study was to determine if single-nucleotide polymorphisms (SNPs) in RANKL, RANK, and OPG influence bone turnover and bone mineral density (BMD) in men. Pairwise tag SNPs (r(2) > or = 0.8) were selected for RANKL, RANK, and OPG and their 10-kb flanking regions. Selected tag SNPs plus five SNPs near RANKL and OPG, associated with BMD in published genome-wide association studies (GWAS), were genotyped in 2653 men aged 40 to 79 years of age recruited for participation in a population-based study of male aging, the European Male Ageing Study (EMAS). N-terminal propeptide of type I procollagen (PINP) and C-terminal cross-linked telopeptide of type I collagen (CTX-I) serum levels were measured in all men. BMD at the calcaneus was estimated by quantitative ultrasound (QUS) in all men. Lumbar spine and total-hip areal BMD (BMD(a)) was measured by dual-energy X-ray absorptiometry (DXA) in a subsample of 620 men. Multiple OPG, RANK, and RANKL SNPs were associated with bone turnover markers. We also identified a number of SNPs associated with BMD, including rs2073618 in OPG and rs9594759 near RANKL. The minor allele of rs2073618 (C) was associated with higher levels of both PINP (beta = 1.83, p = .004) and CTX-I (beta = 17.59, p = 4.74 x 10(-4)), and lower lumbar spine BMD(a) (beta = -0.02, p = .026). The minor allele of rs9594759 (C) was associated with lower PINP (beta = -1.84, p = .003) and CTX-I (beta = -27.02, p = 6.06 x 10(-8)) and higher ultrasound BMD at the calcaneus (beta = 0.01, p = .037). Our findings suggest that genetic variation in the RANKL/RANK/OPG signaling pathway influences bone turnover and BMD in European men.

Gene-gene interactions in RANK/RANKL/OPG system influence bone mineral density in postmenopausal women

Receptor activator of nuclear factor kappaB (RANK) is one of the proteins in regulation of osteoclastogenesis via RANK/RANKL/OPG. Gene that codes for RANK protein (TNFRSF11A) was associated with osteoporotic fractures in a recent genome-wide association study. As variations in the RANK gene could alter its expression and activity, the aim of our study was to evaluate the influence of four RANK gene polymorphisms on bone mineral density (BMD) and biochemical markers. We evaluated 467 postmenopausal women and 117 elderly men. All subjects were genotyped for the presence of RANK polymorphisms -670G>C, +34694C>T, +34901G>A and +35966insdelC. BMD and biochemical markers were measured. Significant associations of +35966insdelC with BMD at lumbar spine (BMD-ls), total hip (BMD-th) and femoral neck (BMD-fn) were found in postmenopausal women (p=0.020, 0.024 and 0.034), but not in men. Significant gene-gene interaction was proved for two RANK polymorphisms in combination with OPG and RANKL polymorphisms studied previously in postmenopausal women. Firstly, RANK/RANKL (+34901G>A/-290C>T) combination was associated with BMD-fn, BMD-th and BMD-ls (p=0.034, 0.016 and 0.050), and secondly, RANK/OPG combination (+35966insdelC/K3N) showed influence on BMD-fn and BMD-ls (p=0.043 and 0.039). Our results suggest that gene-gene interactions between RANK and OPG, and RANK and RANKL influence BMD in postmenopausal women.

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.

Association of IL-15 polymorphisms with bone mineral density in postmenopausal Korean women

Interleukin-15 (IL-15) has been suggested to participate in bone metabolism by stimulating osteoclast differentiation and mediating inflammatory bone loss. This study investigated the effect of IL-15 gene polymorphisms on the bone mineral density (BMD) and bone fracture rates of postmenopausal women. Sequencing of the IL-15 gene in 24 Koreans revealed 16 single-nucleotide polymorphisms (SNPs), of which five were selected for further study. Postmenopausal Korean women (n = 844) were genotyped for these SNPs, and their BMDs and risk of fractures were assessed. It was found that the +20A > G, +13467C > A, +13653A > T, and +13815A > T IL-15 gene polymorphisms were significantly associated with the BMD of the lumbar spine and femoral neck and that their effects were gene-dose dependent. BMD was reduced when the minor allele of +13467A and +13653T or the common allele of +20A and +13815A was present. Haplotype (ht) analyses revealed that ht1 (GCAT) and ht2 (AATA) were associated with BMD of the lumbar spine and femoral neck. However, there was no association between the risk of fracture and IL-15 SNPs or hts. These results suggest that the +20A > G, +13467C > A, +13653A > T, and +13815A > T SNPs in the IL-15 gene affect BMD and, thus, could be genetic markers of osteoporosis.

Association analyses of RANKL/RANK/OPG gene polymorphisms with femoral neck compression strength index variation in Caucasians

Femoral neck compression strength index (fCSI), a novel phenotypic parameter that integrates bone density, bone size, and body size, has significant potential to improve hip fracture risk assessment. The genetic factors underlying variations in fCSI, however, remain largely unknown. Given the important roles of the receptor activator of the nuclear factor-kappaB ligand/receptor activator of the nuclear factor-kappaB/osteoprotegerin (RANKL/RANK/OPG) pathway in the regulation of bone remodeling, we tested the associations between RANKL/RANK/OPG polymorphisms and variations in fCSI as well as its components (femoral neck bone mineral density [fBMD], femoral neck width [FNW], and weight). This was accomplished with a sample comprising 1873 subjects from 405 Caucasian nuclear families. Of the 37 total SNPs studied in these three genes, 3 SNPs, namely, rs12585014, rs7988338, and rs2148073, of RANKL were significantly associated with fCSI (P = 0.0007, 0.0007, and 0.0005, respectively) after conservative Bonferroni correction. Moreover, the three SNPs were approximately in complete linkage disequilibrium. Haplotype-based association tests corroborated the single-SNP results since haplotype 1 of block 1 of the RANKL gene achieved an even more significant association with fCSI (P = 0.0003) than any of the individual SNPs. However, we did not detect any significant associations of these genes with fBMD, FNW, or weight. In summary, our findings suggest that the RANKL gene may play an important role in variation in fCSI, independent of fBMD and non-fBMD components.

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.

Genetics of osteoporosis

Osteoporosis is a frequent skeletal disorder, particularly among postmenopausal women. It affects approximately 30% of women and 12% of men above 50 years of age. It is characterized by reduced bone mass and alterations in bone microarchitecture that result in impaired bone strength and a propensity to fracture. Decreased bone mass is the consequence of an imbalance in the bone remodeling process, resulting from complex interactions between acquired and genetic factors. The former include physical activity, nutrition and other lifestyle habits, as well as the skeletal effects of some diseases and drug therapies. Genetic factors have been extensively studied during the past 15 years. We will review some important studies that exemplify the advances and the difficulties in this research field.

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.