OPG and RANK polymorphisms are both associated with cortical bone mineral density: findings from a metaanalysis of the Avon longitudinal study of parents and children and gothenburg osteoporosis and obesity determinants cohorts

Endocrine Health and Health Care Disparities in the Pediatric and Sexual and Gender Minority Populations: An Endocrine Society Scientific Statement

Endocrine care of pediatric and adult patients continues to be plagued by health and health care disparities that are perpetuated by the basic structures of our health systems and research modalities, as well as policies that impact access to care and social determinants of health. This scientific statement expands the Society's 2012 statement by focusing on endocrine disease disparities in the pediatric population and sexual and gender minority populations. These include pediatric and adult lesbian, gay, bisexual, transgender, queer, intersex, and asexual (LGBTQIA) persons. The writing group focused on highly prevalent conditions-growth disorders, puberty, metabolic bone disease, type 1 (T1D) and type 2 (T2D) diabetes mellitus, prediabetes, and obesity. Several important findings emerged. Compared with females and non-White children, non-Hispanic White males are more likely to come to medical attention for short stature. Racially and ethnically diverse populations and males are underrepresented in studies of pubertal development and attainment of peak bone mass, with current norms based on European populations. Like adults, racial and ethnic minority youth suffer a higher burden of disease from obesity, T1D and T2D, and have less access to diabetes treatment technologies and bariatric surgery. LGBTQIA youth and adults also face discrimination and multiple barriers to endocrine care due to pathologizing sexual orientation and gender identity, lack of culturally competent care providers, and policies. Multilevel interventions to address these disparities are required. Inclusion of racial, ethnic, and LGBTQIA populations in longitudinal life course studies is needed to assess growth, puberty, and attainment of peak bone mass. Growth and development charts may need to be adapted to non-European populations. In addition, extension of these studies will be required to understand the clinical and physiologic consequences of interventions to address abnormal development in these populations. Health policies should be recrafted to remove barriers in care for children with obesity and/or diabetes and for LGBTQIA children and adults to facilitate comprehensive access to care, therapeutics, and technological advances. Public health interventions encompassing collection of accurate demographic and social needs data, including the intersection of social determinants of health with health outcomes, and enactment of population health level interventions will be essential tools.

Bone mineral density in high-level endurance runners: Part B-genotype-dependent characteristics

Purpose

Inter-individual variability in bone mineral density (BMD) exists within and between endurance runners and non-athletes, probably in part due to differing genetic profiles. Certainty is lacking, however, regarding which genetic variants may contribute to BMD in endurance runners and if specific genotypes are sensitive to environmental factors, such as mechanical loading via training.

Method

Ten single-nucleotide polymorphisms (SNPs) were identified from previous genome-wide and/or candidate gene association studies that have a functional effect on bone physiology. The aims of this study were to investigate (1) associations between genotype at those 10 SNPs and bone phenotypes in high-level endurance runners, and (2) interactions between genotype and athlete status on bone phenotypes.

Results

Female runners with P2RX7 rs3751143 AA genotype had 4% higher total-body BMD and 5% higher leg BMD than AC + CC genotypes. Male runners with WNT16 rs3801387 AA genotype had 14% lower lumbar spine BMD than AA genotype non-athletes, whilst AG + GG genotype runners also had 5% higher leg BMD than AG + GG genotype non-athletes.

Conclusion

We report novel associations between P2RX7 rs3751143 genotype and BMD in female runners, whilst differences in BMD between male runners and non-athletes with the same WNT16 rs3801387 genotype existed, highlighting a potential genetic interaction with factors common in endurance runners, such as high levels of mechanical loading. These findings contribute to our knowledge of the genetic associations with BMD and improve our understanding of why some runners have lower BMD than others.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00421-021-04789-z.

Level and change in bone microarchitectural parameters and their relationship with previous fracture and established bone mineral density loci

BACKGROUND

Osteoporosis is characterised by a reduction of bone mineral density (BMD) and predisposition to fracture. Bone microarchitecture, measured by high resolution peripheral quantitative computed tomography (HR-pQCT), has been related to fragility fractures and BMD and has been the subject of large-scale genome-wide analysis. We investigated whether fracture was related to baseline values and longitudinal changes in bone microarchitecture and whether bone microarchitecture was associated with established BMD loci.

METHODS

115 males and 99 females (aged 72-81 at baseline) from the Hertfordshire Cohort Study (HCS) were analysed. Fracture history was determined in 2011-2012 by self-report and vertebral fracture assessment. Participants underwent HR-pQCT scans of the distal radius and tibia in 2011-2012 and 2017. Previous fracture in relation to baseline values and changes in tibial HR-pQCT parameters was examined using sex-adjusted logistic regression with and without adjustment for age, sociodemographic, lifestyle and clinical characteristics; baseline values and changes in parameters associated with previous fracture were then examined in relation to four established BMD loci after adjustment for sex and age.

RESULTS

Previous fracture was related to: higher trabecular area (fully-adjusted odds ratio [95% CI] per SD greater baseline value: 2.18 [1.27,3.73], p = 0.005); lower total volumetric BMD (0.53 [0.34,0.84], p = 0.007), cortical area (0.53 [0.30,0.95], p = 0.032), cortical BMD (0.56 [0.36,0.88], p = 0.011) and cortical thickness (0.45 [0.27,0.77], p = 0.004); and greater declines in trabecular BMD (p = 0.001). Associations were robust in sex- and fully-adjusted analysis. Relationships between BMD loci and these HR-pQCT parameters were weak: rs3801387 (WNT16) was related to decline in trabecular BMD (p = 0.011) but no other associations were significant (p > 0.05).

CONCLUSION

Baseline values of HR-pQCT parameters and greater decline in trabecular BMD were associated with fracture. Change in trabecular BMD was associated with WNT16 which has been demonstrated to influence bone health in murine models and human genome-wide association studies (GWAS).

Opportunities and Challenges in Functional Genomics Research in Osteoporosis: Report From a Workshop Held by the Causes Working Group of the Osteoporosis and Bone Research Academy of the Royal Osteoporosis Society on October 5th 2020

The discovery that sclerostin is the defective protein underlying the rare heritable bone mass disorder, sclerosteosis, ultimately led to development of anti-sclerostin antibodies as a new treatment for osteoporosis. In the era of large scale GWAS, many additional genetic signals associated with bone mass and related traits have since been reported. However, how best to interrogate these signals in order to identify the underlying gene responsible for these genetic associations, a prerequisite for identifying drug targets for further treatments, remains a challenge. The resources available for supporting functional genomics research continues to expand, exemplified by "multi-omics" database resources, with improved availability of datasets derived from bone tissues. These databases provide information about potential molecular mediators such as mRNA expression, protein expression, and DNA methylation levels, which can be interrogated to map genetic signals to specific genes based on identification of causal pathways between the genetic signal and the phenotype being studied. Functional evaluation of potential causative genes has been facilitated by characterization of the "osteocyte signature", by broad phenotyping of knockout mice with deletions of over 7,000 genes, in which more detailed skeletal phenotyping is currently being undertaken, and by development of zebrafish as a highly efficient additional in vivo model for functional studies of the skeleton. Looking to the future, this expanding repertoire of tools offers the hope of accurately defining the major genetic signals which contribute to osteoporosis. This may in turn lead to the identification of additional therapeutic targets, and ultimately new treatments for osteoporosis.

Factors influencing peak bone mass gain

Bone mass is a key determinant of osteoporosis and fragility fractures. Epidemiologic studies have shown that a 10% increase in peak bone mass (PBM) at the population level reduces the risk of fracture later in life by 50%. Low PBM is possibly due to the bone loss caused by various conditions or processes that occur during adolescence and young adulthood. Race, gender, and family history (genetics) are responsible for the majority of PBM, but other factors, such as physical activity, calcium and vitamin D intake, weight, smoking and alcohol consumption, socioeconomic status, age at menarche, and other secondary causes (diseases and medications), play important roles in PBM gain during childhood and adolescence. Hence, the optimization of lifestyle factors that affect PBM and bone strength is an important strategy to maximize PBM among adolescents and young people, and thus to reduce the low bone mass or osteoporosis risk in later life. This review aims to summarize the available evidence for the common but important factors that influence bone mass gain during growth and development and discuss the advances of developing high PBM.

Associations of LRP5 Gene With Bone Mineral Density, Bone Turnover Markers, and Fractures in the Elderly With Osteoporosis

Objective: The study aimed to explore the associations of rs4988300 and rs634008 in the low-density lipoprotein receptor-related protein 5 (LRP5) gene with bone mineral density (BMD), bone turnover markers (BTM), and fractures in elderly patients with osteoporosis (OP). Methods: Our study included 328 unrelated OP patients with or without fractures. Genomic DNA was extracted for genotyping. BTM levels were assessed by electrochemiluminescence (ECL). Dual-energy X-ray absorptiometry (DXA) was employed to measure BMD in the lumbar spine (LS) and proximal femur. Basic features between the OP and fracture groups were analyzed using the t-test. The Chi-square test was performed to analyze the differences in allele and genotype frequencies. The associations of single-nucleotide polymorphisms (SNPs) with BMD and BTM in the subgroups were investigated by the analysis of covariance (ANCOVA) adjusted for confounding factors. Results: In both females and males, individuals with fractures exhibited higher BTM levels and lower BMD values than those with OP (P < 0.05). The allele and genotype frequencies of rs4988300 in the subgroups were significantly different (P < 0.05). In both females and males suffering from OP, participants with rs4988300 GG or rs634008 TT presented lower procollagen I N-terminal propeptide (PINP) levels (P < 0.05). Women with OP carrying rs4988300 GG exhibited lower BMD values at FN and TH (P < 0.05). In both females and males with fractures, individuals carrying rs4988300 GG genotype or rs634008 TT genotype exhibited lower PINP levels and BMD values at FN and TH than those with other genotypes (P < 0.05). Conclusions: Rs4988300 and rs634008 polymorphisms in the LRP5 gene are associated with bone phenotypes in the elderly with OP or fractures.

High risk of lumbar spine osteoporosis with the RANK rs3018362 polymorphism

Osteoporosis is characterized by reduced bone mineral density (BMD) and quality, increasing the risk of fractures. A large number of genes involved in bone metabolism have been implicated in the genesis of osteoporosis; these include RANK and RANKL. Polymorphisms of these genes have been implicated in osteoporosis. The aim of this study was to determine the association of the RANK rs3018362 and RANKL rs12585014 polymorphisms with risk of osteoporosis. Four hundred Mexican women aged 40 years old or above were genotyped by real-time PCR and several demographic and risk factors were explored. The GA and AA genotypes of the rs3018362 polymorphism were associated with a high risk of osteoporosis in the dominant model (p=.0062; OR = 2.16, 95% CI: 1.24-3.78). In summary, the rs3018362 polymorphism in the RANK gene seems to be associated with osteoporosis of the lumbar spine while the RANKL rs12585014 is not, although more studies are needed to confirm these results.

pQCT bone geometry and strength: population epidemiology and concordance in Australian children aged 11-12 years and their parents

OBJECTIVES

To describe the epidemiology and concordance of bone health in a population-based sample of Australian parent-child dyads at child age 11-12 years.

DESIGN

Population-based cross-sectional study (the Child Health CheckPoint) nested between waves 6 and 7 of the Longitudinal Study of Australian Children (LSAC).

SETTING

Assessment centres in seven cities around Australia, February 2015-March 2016.

PARTICIPANTS

of all participating CheckPoint families (n=1874), bone data were available for 1222 dyads (1271 children, 50% girls; 1250 parents, 86% mothers).

OUTCOME MEASURES

Peripheral quantitative CT (pQCT) of the non-dominant leg scanned at the 4% (distal) and 66% (mid-calf) tibial sites. Stratec XCT 2000 software generated estimates of bone density, geometry and polar stress-strain index.Parent-child concordance were assessed using Pearson's correlation coefficients and multivariable linear regression models. Percentiles were determined using survey weights. Survey weights and methods accounted for LSAC's complex sampling, stratification and clustering within postcodes.

RESULTS

Concordances were greater for the geometric pQCT parameters (periosteal circumference 0.38, 95% CI 0.33 to 0.43; endosteal circumference 0.42, 95% CI 0.37 to 0.47; total cross-sectional area 0.37, 95% CI 0.32 to 0.42) than density (cortical density 0.25, 95% CI 0.19 to 0.30). Mother-child and father-child values were similar. Relationships attenuated only slightly on adjustment for age, sex and body mass index. Percentiles and concordance are presented for the whole sample and by sex.

CONCLUSIONS

There is strong parent-child concordance in bone geometry and, to a lesser extent, density even before the period of peak adolescent bone deposition. This geometrical concordance suggests that future intergenerational bone studies could consider using pQCT rather than the more commonly used dual X-ray absorptiometry (DXA).

The RANKL rs12585014 polymorphism is associated with age at menarche in postmenopausal women with hip fracture

The RANK/RANKL/OPG signaling is important in the regulation of bone turnover. The aim of the present work was to analyze the rs3018362 and rs12585014 polymorphisms in the RANK and RANKL genes, as well as risk factors in postmenopausal women. Women with hip fracture, with femoral neck osteoporosis and controls (n = 646) were recruited. From these, 303 women who fulfill the inclusion criteria were genotyped using real-time PCR with TaqMan probes. There were no associations of the rs3018362 and rs12585014 with osteoporosis or fracture. When women were divided by age at menarche, the rs12585014 GG genotype was strongly associated with age at menarche >13 years [p = .00774, OR = 6.429 (1.907-21.103)] in women with hip fracture. Significant differences in risk factors such as body mass index, age at menopause, use of estrogens, the presence of hypertension, and diabetes mellitus were found. Carrying the GG genotype of rs12585014 entails a higher risk of having menarche later (>13 years), which could involves a greater risk of fractures. The rs3018362 and rs12585014 do not seem to be associated with hip osteoporosis or hip fracture in Mexican women.

SNPs in the vicinity of P2X7R, RANK/RANKL/OPG and Wnt signalling pathways and their association with bone phenotypes in academy footballers

CONTEXT

Genotype plays an important role in influencing bone phenotypes, such as bone mineral density, but the role of genotype in determining responses of bone to exercise has yet to be elucidated.

OBJECTIVE

To determine whether 10 SNPs associated with genes in the vicinity of P2X7R, RANK/RANKL/OPG and Wnt Signalling Pathways are associated with bone phenotypes in elite academy footballers (Soccer players) and to determine whether these genotypes are associated with training induced changes in bone. Design, participants, and methods: 99 elite academy footballers volunteered to participate. Peripheral computed tomography of the tibia (4%, 14%, 38% and 66% sites) was performed immediately before and 12 weeks after an increase in football training volume. Genotypes were determined using proprietary fluorescence-based competitive allele-specific PCR assays.

RESULTS

No significant genotype by time interactions were shown for any of the SNPs analysed (P > .05). A main effect of genotype was shown. SOST SNP rs1877632 (trabecular density), P2X7R SNPs rs1718119 (cortical thickness and CSA), rs3751143 (SSI, CSA, cortical CSA and periosteal circumference) RANK/RANKL/OPG SNPs rs9594738 (periosteal circumference), rs1021188 (cortical thickness and CSA) and rs9594759 (cortical density) were associated with bone phenotypes (P < .05).

CONCLUSIONS

No association was shown between P2X7R, RANK/RANKL/OPG and Wnt Signalling SNPs and a change in bone phenotypes following 12 weeks of increased training volume in elite academy footballers. However, SNPs were associated with bone phenotypes pre training. These data highlight the complexity of the interaction between SNPs in the vicinity of the RANK/RANKL/OPG, P2X7R and Wnt metabolic regulatory pathways and bone phenotypes in elite academy footballers.

Increased detection of genetic loci associated with risk predictors of osteoporotic fracture using a pleiotropic cFDR method

Although GWAS have been successful in identifying some osteoporosis associated loci, the findings explain only a small fraction of the total genetic variance. In this study we use a recently developed novel pleiotropic conditional false discovery rate (cFDR) method to identify novel genetic loci associated with two risk traits for osteoporotic fracture (the clinical outcome and end result of osteoporosis), Height (HT) and Femoral Neck (FNK) BMD. The cFDR method allows us to improve the detection of associated variants by incorporating any potentially shared genetic mechanisms between the two associated traits. We analyzed the summary statistics from two GWAS meta-analyses for single nucleotide polymorphisms (SNPs) that are associated with HT and FNK BMD. Using the cFDR method, we show enrichment in the identification of SNPs associated with each trait conditioned on their strength of association with the second trait. The findings revealed 18 SNPs that are associated with both HT and FNK BMD, 4 of which had not previously been reported to play a role in bone health. The novel SNPs located at KIF1B and the intergenic region between FERD3L and TWISTNB are noteworthy as these genes may be associated with processes that are functionally important in bone metabolism. By leveraging GWAS results from related phenotypes we identified several novel loci that may contribute to the proportion of variability explained for each trait, although we cannot speculate about these potential contributions to heritability based on this analysis alone.

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.

Osteoporosis and Bone Mass Disorders: From Gene Pathways to Treatments

Genomic technologies have evolved rapidly contributing to the understanding of diseases. Genome-wide association studies (GWAS) and whole-exome sequencing have aided the identification of the genetic determinants of monogenic and complex conditions including osteoporosis and bone mass disorders. Overlap exists between the genes implicated in monogenic and complex forms of bone mass disorders, largely explained by the clustering of genes encoding factors in signaling pathways crucial for mesenchymal cell differentiation, skeletal development, and bone remodeling and metabolism. Numerous of the remaining discovered genes merit functional follow-up studies to elucidate their role in bone biology. The insight provided by genetic studies is serving the identification of biomarkers predictive of disease, redefining disease, response to treatment, and discovery of novel drug targets for skeletal disorders.

Genetics of Bone Mass in Childhood and Adolescence: Effects of Sex and Maturation Interactions

We aimed to determine if adult bone mineral density (BMD) susceptibility loci were associated with pediatric bone mass and density, and if sex and pubertal stage influenced any association. We analyzed prospective areal BMD (aBMD) and bone mineral content (BMC) data from the Bone Mineral Density in Childhood Study (n = 603, European ancestry, 54% female). Linear mixed models were used to assess if 77 single-nucleotide polymorphisms (SNPs) near known adult BMD susceptibility loci interacted with sex and pubertal stage to influence the aBMD/BMC; adjusting for age, BMI, physical activity, and dietary calcium. The strongest main association was observed between an SNP near C7orf58 and distal radius aBMD. However, this association had a significant sex • SNP interaction, revealing a significant association only in females (b = -0.32, p = 1.8 × 10(-6)). Furthermore, the C12orf23 locus had significant interactions with both sex and pubertal stage, revealing associations in females during Tanner stage I for total hip aBMD (b = 0.24, p = 0.001) and femoral neck aBMD (b = 0.27, p = 3.0 × 10(-5)). In contrast, the sex • SNP interactions for loci near LRP5 and WNT16 uncovered associations that were only in males for total body less head BMC (b = 0.22, p = 4.4 × 10(-4)) and distal radius aBMD (b = 0.27, p = 0.001), respectively. Furthermore, the LRP5 locus interacted with both sex and pubertal stage, demonstrating associations that were exclusively in males during Tanner V for total hip aBMD (b = 0.29, p = 0.003). In total, significant sex • SNP interactions were found at 15 loci; pubertal stage • SNP interactions at 23 loci and 19 loci interacted with both sex and pubertal stage. In conclusion, variants originally associated with adult BMD influence bone mass in children of European ancestry, highlighting the fact that many of these loci operate early in life. However, the direction and magnitude of associations for a large number of SNPs only became evident when accounting for sex and maturation.

Novel candidate genes putatively involved in stress fracture predisposition detected by whole-exome sequencing

While genetic factors in all likelihood contribute to stress fracture (SF) pathogenesis, a few studies focusing on candidate genes have previously been reported. The objective of this study is to gain better understanding on the genetic basis of SF in a gene-naive manner. Exome sequence capture followed by massive parallel sequencing of two pooled DNA samples from Israeli combat soldiers was employed: cases with high grade SF and ethnically matched healthy controls. The resulting sequence variants were individually verified using the Sequenom™ platform and the contribution of the genetic alterations was validated in a second cohort of cases and controls. In the discovery set that included DNA pool of cases (n = 34) and controls (n = 60), a total of 1174 variants with >600 reads/variant/DNA pool were identified, and 146 (in 127 genes) of these exhibited statistically significant (P < 0·05) different rates between SF cases and controls after multiple comparisons correction. Subsequent validation of these 146 sequence variants individually in a total of 136 SF cases and 127 controls using the Sequenom™ platform validated 20/146 variants. Of these, three missense mutations (rs7426114, rs4073918, rs3752135 in the NEB, SLC6A18 and SIGLEC12 genes, respectively) and three synonymous mutations (rs2071856, rs2515941, rs716745 in the ELFN2, GRK4, LRRC55 genes) displayed significant different rates in SF cases compared with controls. Exome sequencing seemingly unravelled novel candidate genes as involved in SF pathogenesis and predisposition.

RANK/RANKL/OPG pathway: genetic associations with stress fracture period prevalence in elite athletes

CONTEXT

The RANK/RANKL/OPG signalling pathway is important in the regulation of bone turnover, with single nucleotide polymorphisms (SNPs) in genes within this pathway associated with bone phenotypic adaptations.

OBJECTIVE

To determine whether four SNPs associated with genes in the RANK/RANKL/OPG signalling pathway were associated with stress fracture injury in elite athletes.

DESIGN, PARTICIPANTS, AND METHODS

Radiologically confirmed stress fracture history was reported in 518 elite athletes, forming the Stress Fracture Elite Athlete (SFEA) cohort. Data were analysed for the whole group and were sub-stratified into male and cases of multiple stress fracture groups. Genotypes were determined using proprietary fluorescence-based competitive allele-specific PCR assays.

RESULTS

SNPs rs3018362 (RANK) and rs1021188 (RANKL) were associated with stress fracture injury (P<0.05). 8.1% of the stress fracture group and 2.8% of the non-stress fracture group were homozygote for the rare allele of rs1021188. Allele frequency, heterozygotes and homozygotes for the rare allele of rs3018362 were associated with stress fracture period prevalence (P<0.05). Analysis of the male only group showed 8.2% of rs1021188 rare allele homozygotes had suffered a stress fracture whilst 2.5% of the non-stress fracture group were homozygous. In cases of multiple stress fractures, homozygotes for the rare allele of rs1021188 and individuals possessing at least one copy of the rare allele of rs4355801 (OPG) were shown to be associated with stress fracture injury (P<0.05).

CONCLUSIONS

The data support an association between SNPs in the RANK/RANKL/OPG signalling pathway and the development of stress fracture injury. The association of rs3018362 (RANK) and rs1021188 (RANKL) with stress fracture injury susceptibility supports their role in the maintenance of bone health and offers potential targets for therapeutic interventions.

Association analysis between g.18873C>T and g.27522G>A genetic polymorphisms of OPG and bone mineral density in Chinese postmenopausal women

Several studies report that the OPG is an important candidate gene in the pathogenesis of osteoporosis. This study aimed to detect the potential association of OPG gene polymorphisms with osteoporosis in postmenopausal women. We recruited 928 subjects containing 463 with primary postmenopausal osteoporosis and 465 healthy volunteers as controls. The BMD of neck hip, lumbar spine (L_2–4), and total hip were assessed by dual-energy X-ray absorptiometry (DEXA). Through the created restriction site-polymerase chain reaction (CRS-PCR), PCR-restriction fragment length polymorphism (PCR-RFLP), and DNA sequencing methods, the g.18873C>T and g.27522G>A have been investigated. As for g.18873C>T, our data indicated that subjects with CC genotype have significantly higher BMD value than those of CT and TT genotypes (all P values < 0.05). As for g.27522G>A, the BMD values of subjects with GG genotype were significantly higher than those of GA and AA genotypes (all P values < 0.05). Our findings suggest that the OPG g.18873C>T and g.27522G>A genetic polymorphisms are associated with the decreased risk for osteoporosis in Chinese postmenopausal women.

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.

The influence of the genetic and non-genetic factors on bone mineral density and osteoporotic fractures in Chinese women

To investigate the effects of genetic and non-genetic factors on bone mineral densities (BMDs) and osteoporotic fractures. This was a cross-sectional study to investigate the relationships between 18 SNPs and non-genetic factors with BMDs and osteoporotic fractures in 1012 Chinese Han women. Five SNPs in genes GPR177, CTNNB1, MEF2C, SOX6, and TNFRSF11B were associated with L1-4 or total hip BMDs. rs11898505 in SPTBN1 gene was associated with osteoporotic fractures. Subjects carrying the largest number of risk alleles (highest 10 %) not only had lower BMD values as compared to those carrying the least number of risk alleles (lowest 10 %), they also had a higher risk of fracture [P = 0.002, OR = 2.252, 95 %CI (1.136, 4.463)]. Results from multivariate stepwise regression analysis revealed that age [P < 0.001, OR = 1.038, 95 % CI (1.018, 1.058)], number of falls in a year [P < 0.001, OR = 2.347, 95 % CI (1.459, 3.774)], the G risk allele in rs11898505 [P = 0.023, OR = 1.559, 95 % CI (1.062, 2.290)], and the L1-4 BMD [P = 0.017, OR = 0.286, 95 % CI (0.102, 0.798)] were associated with the occurrence of osteoporotic fractures. Genetic (rs11898505) and non-genetic factors (age, number of falls in a year and L1-4 BMD) could work in concert to contribute to the risk of osteoporotic fractures.

The relationship between adipose tissue and bone metabolism

OBJECTIVES

The authors have set out to evaluate the literature relevant to the dynamic regulation of adipogenesis and osteogenesis.

DESIGN AND METHODS

A detailed search of the past and recent literature was conducted on Pubmed using a combination of keywords including: adipogenesis, bone marrow, hematopoiesis, mesenchymal stromal/stem cell, and osteogenesis.

RESULTS

Throughout one's lifespan, the bone marrow microenvironment provides a unique niche for mesenchymal stromal/stem cells (BMSCs) and hematopoietic stem cells (HSCs). The marrow changes as a function of biological age and pathophysiology. Historically, clinical biochemistry has observed these changes from an HSC and hematological perspective. Nevertheless, these changes also reflect the balance between BMSC adipogenic and osteogenic processes which can display an inverse or reciprocal relationship. Multiple hormonal factors and nuclear hormone receptor ligands and drugs are responsible for BMSC lineage selection. Data from a number of laboratories now implicates endocrine feedback loops between extramedullary adipose depots and the central nervous system.

CONCLUSIONS

This concise review provides a perspective on the mechanisms regulating BMSC differentiation in the context of biological aging, obesity, and osteoporosis.

"Single nucleotide polymorphisms of the OPG/RANKL system genes in primary hyperparathyroidism and their relationship with bone mineral density"

Background

Primary hyperparathyroidism (PHPT) affects mainly cortical bone. It is thought that parathyroid hormone (PTH) indirectly regulates the activity of osteoclasts by means of the osteoprotegerin/ligand of the receptor activator of nuclear factor-κβ (OPG/RANKL) system. Several studies have confirmed that OPG (osteoprotegerin) and RANKL (ligand of the receptor activator of nuclear factor-κβ) loci are determinants of bone mineral density (BMD) in the general population. The aim of this study is to analyze the relationship between fractures and BMD and the rs3102735 (163 A/G), rs3134070 (245 T/G) and rs2073618 (1181 G/C) SNPs of the OPG and the rs2277438 SNP of the RANKL, in patients with sporadic PHPT.

Methods

We enrolled 298 Caucasian patients with PHPT and 328 healthy volunteers in a cross-sectional study. We analyzed anthropometric data, history of fractures or renal lithiasis, biochemical determinants including markers for bone remodelling, BMD measurements in the lumbar spine, total hip, femoral neck and distal radius, and genotyping for the SNPs to be studied.

Results

Regarding the age of diagnosis, BMI, menopause status, frequency of fractures or renal lithiasis, we found no differences between genotypes in any of the SNPs studied in the PHPT group. Significant lower BMD in the distal radius with similar PTH levels was found in the minor allele homozygotes (GG) compared to heterozygotes and major allele homozygotes in both OPG rs3102735 (163 A/G) and OPG rs3134070 (245 T/G) SNPs in those with PHPT compared to control subjects. We found no differences between genotypes of the OPG rs2073618 (1181 G/C) SNP with regard to BMD in the PHPT subjects. In the evaluation of rs2277438 SNP of the RANKL in PHPT patients, we found a non significant trend towards lower BMD in the 1/3 distal radius and at total hip in the minor allele homocygotes (GG) genotype group versus heterocygotes and major allele homocygotes (AA).

Conclusions

Our study provides the first evaluation of the relationship between SNPs of the OPG/RANK system and sporadic PHPT. Subjects with PHPT and minor homocygote genotype (GG) for the OPG rs3102735 (163 A/G) and OPG rs3134070 (245 T/G) SNPs have lower BMD in the distal radius, and this association does not appear to be mediated by differences in PTH serum levels.

Influence of polymorphisms in the RANKL/RANK/OPG signaling pathway on volumetric bone mineral density and bone geometry at the forearm in men

We sought to determine the influence of single-nucleotide polymorphisms (SNPs) in RANKL, RANK, and OPG on volumetric bone mineral density (vBMD) and bone geometry at the radius in men. Pairwise tag SNPs (r (2) ≥ 0.8) for RANKL (n = 8), RANK (n = 44), and OPG (n = 22) and five SNPs near RANKL and OPG strongly associated with areal BMD in genomewide association studies were previously genotyped in men aged 40-79 years in the European Male Ageing Study (EMAS). Here, these SNPs were analyzed in a subsample of men (n = 589) who had peripheral quantitative computed tomography (pQCT) performed at the distal (4%) and mid-shaft (50%) radius. Estimated parameters were total and trabecular vBMD (mg/mm(3)) and cross-sectional area (mm(2)) at the 4% site and cortical vBMD (mg/mm(3)); total, cortical, and medullary area (mm(2)); cortical thickness (mm); and stress strain index (SSI) (mm(3)) at the 50% site. We identified 12 OPG SNPs associated with vBMD and/or geometric parameters, including rs10505348 associated with total vBMD (β [95% CI] = 9.35 [2.12-16.58], P = 0.011), cortical vBMD (β [95% CI] = 5.62 [2.10-9.14], P = 0.002), cortical thickness (β [95% CI] = 0.08 [0.03-0.13], P = 0.002), and medullary area (β [95% CI] = -2.90 [-4.94 to -0.86], P = 0.005) and rs2073618 associated with cortical vBMD (β [95% CI] = -4.30 [-7.78 to -0.82], P = 0.015) and cortical thickness (β [95% CI] = -0.08 [-0.13 to -0.03], P = 0.001). Three RANK SNPs were associated with vBMD, including rs12956925 associated with trabecular vBMD (β [95% CI] = -7.58 [-14.01 to -1.15], P = 0.021). There were five RANK SNPs associated with geometric parameters, including rs8083511 associated with distal radius cross-sectional area (β [95% CI] = 8.90 [0.92-16.88], P = 0.029). No significant association was observed between RANKL SNPs and pQCT parameters. Our findings suggest that genetic variation in OPG and RANK influences radius vBMD and geometry in men.

TNFRSF11B gene haplotype and its association with bone mineral density variations in postmenopausal Mexican-Mestizo women

OBJECTIVE

Osteoporosis is a complex health disease characterized by low bone mineral density (BMD), which is determined by an interaction of genetics with metabolic and environmental factors. The tumor necrosis factor receptor superfamily, member 11b (TNFRSF11B) gene, has been investigated in relation to BMD. Three polymorphisms in/nearby TNFRSF11B have been associated with BMD variations in some populations. The aim of this study was to investigate the possible association among three SNPs of TNFRSF11B and their haplotypes with the presence of BMD variations in postmenopausal Mexican Mestizo women.

SUBJECTS AND METHODS

One thousand unrelated postmenopausal women of Mexican-Mestizo ethnic origin, who attended the outpatient clinic for routine, general medical evaluation, were invited and 750 women accepted to participate in the study. A structured questionnaire for risk factors was applied and BMD was measured in total hip and lumbar spine by dual-energy X-ray absorptiometry. DNA was obtained from blood leukocytes. Three single-nucleotide polymorphisms in TNFRSF11B gene were studied: rs4355801, rs2073618, and rs6993813. Real-time PCR allelic discrimination was used for genotyping. Deviations from Hardy-Weinberg equilibrium were tested. Pairwise linkage disequilibrium between single nucleotide polymorphisms was calculated by direct correlation r(2), and haplotype analysis was conducted.

RESULTS

Of the subjects, 31% had osteoporosis, 45.1% had osteopenia, and 23.9% had normal BMD. Genotype and allele distributions showed no significant differences; however, A-G-T haplotype was associated with variations in femoral neck BMD (P=0.022).

CONCLUSIONS

In our study population, analysis of the haplotypes of TNFRSF11B is a better genetic marker for variations in BMD.

The ESR1 (6q25) locus is associated with calcaneal ultrasound parameters and radial volumetric bone mineral density in European men

Purpose

Genome-wide association studies (GWAS) have identified 6q25, which incorporates the oestrogen receptor α gene (ESR1), as a quantitative trait locus for areal bone mineral density (BMD_a) of the hip and lumbar spine. The aim of this study was to determine the influence of this locus on other bone health outcomes; calcaneal ultrasound (QUS) parameters, radial peripheral quantitative computed tomography (pQCT) parameters and markers of bone turnover in a population sample of European men.

Methods

Eight single nucleotide polymorphisms (SNP) in the 6q25 locus were genotyped in men aged 40–79 years from 7 European countries, participating in the European Male Ageing Study (EMAS). The associations between SNPs and measured bone parameters were tested under an additive genetic model adjusting for centre using linear regression.

Results

2468 men, mean (SD) aged 59.9 (11.1) years had QUS measurements performed and bone turnover marker levels measured. A subset of 628 men had DXA and pQCT measurements. Multiple independent SNPs showed significant associations with BMD using all three measurement techniques. Most notably, rs1999805 was associated with a 0.10 SD (95%CI 0.05, 0.16; p = 0.0001) lower estimated BMD at the calcaneus, a 0.14 SD (95%CI 0.05, 0.24; p = 0.004) lower total hip BMD_a, a 0.12 SD (95%CI 0.02, 0.23; p = 0.026) lower lumbar spine BMD_a and a 0.18 SD (95%CI 0.06, 0.29; p = 0.003) lower trabecular BMD at the distal radius for each copy of the minor allele. There was no association with serum levels of bone turnover markers and a single SNP which was associated with cortical density was also associated with cortical BMC and thickness.

Conclusions

Our data replicate previous associations found between SNPs in the 6q25 locus and BMD_a at the hip and extend these data to include associations with calcaneal ultrasound parameters and radial volumetric BMD.

Genome-wide association meta-analysis of cortical bone mineral density unravels allelic heterogeneity at the RANKL locus and potential pleiotropic effects on bone

Previous genome-wide association (GWA) studies have identified SNPs associated with areal bone mineral density (aBMD). However, this measure is influenced by several different skeletal parameters, such as periosteal expansion, cortical bone mineral density (BMD(C)) cortical thickness, trabecular number, and trabecular thickness, which may be under distinct biological and genetic control. We have carried out a GWA and replication study of BMD(C), as measured by peripheral quantitative computed tomography (pQCT), a more homogenous and valid measure of actual volumetric bone density. After initial GWA meta-analysis of two cohorts (ALSPAC n = 999, aged ∼15 years and GOOD n = 935, aged ∼19 years), we attempted to replicate the BMD(C) associations that had p<1×10(-5) in an independent sample of ALSPAC children (n = 2803) and in a cohort of elderly men (MrOS Sweden, n = 1052). The rs1021188 SNP (near RANKL) was associated with BMD(C) in all cohorts (overall p = 2×10(-14), n = 5739). Each minor allele was associated with a decrease in BMD(C) of ∼0.14SD. There was also evidence for an interaction between this variant and sex (p = 0.01), with a stronger effect in males than females (at age 15, males -6.77mg/cm(3) per C allele, p = 2×10(-6); females -2.79 mg/cm(3) per C allele, p = 0.004). Furthermore, in a preliminary analysis, the rs1021188 minor C allele was associated with higher circulating levels of sRANKL (p<0.005). We show this variant to be independent from the previously aBMD associated SNP (rs9594738) and possibly from a third variant in the same RANKL region, which demonstrates important allelic heterogeneity at this locus. Associations with skeletal parameters reflecting bone dimensions were either not found or were much less pronounced. This finding implicates RANKL as a locus containing variation associated with volumetric bone density and provides further insight into the mechanism by which the RANK/RANKL/OPG pathway may be involved in skeletal development.