IL21R and PTH may underlie variation of femoral neck bone mineral density as revealed by a genome-wide association study

Twelve years of GWAS discoveries for osteoporosis and related traits: advances, challenges and applications

Osteoporosis is a common skeletal disease, affecting ~200 million people around the world. As a complex disease, osteoporosis is influenced by many factors, including diet (e.g. calcium and protein intake), physical activity, endocrine status, coexisting diseases and genetic factors. In this review, we first summarize the discovery from genome-wide association studies (GWASs) in the bone field in the last 12 years. To date, GWASs and meta-analyses have discovered hundreds of loci that are associated with bone mineral density (BMD), osteoporosis, and osteoporotic fractures. However, the GWAS approach has sometimes been criticized because of the small effect size of the discovered variants and the mystery of missing heritability, these two questions could be partially explained by the newly raised conceptual models, such as omnigenic model and natural selection. Finally, we introduce the clinical use of GWAS findings in the bone field, such as the identification of causal clinical risk factors, the development of drug targets and disease prediction. Despite the fruitful GWAS discoveries in the bone field, most of these GWAS participants were of European descent, and more genetic studies should be carried out in other ethnic populations to benefit disease prediction in the corresponding population.

Comprehensive Genetic Analysis of 128 Candidate Genes in a Cohort With Idiopathic, Severe, or Familial Osteoporosis

Context

The genetic bases of osteoporosis (OP), a disorder with high heritability, are poorly understood at an individual level. Cases of idiopathic or familial OP have long puzzled clinicians as to whether an actionable genetic cause could be identified.

Objective

We performed a genetic analysis of 28 cases of idiopathic, severe, or familial osteoporosis using targeted massively parallel sequencing.

Design

Targeted sequencing of 128 candidate genes was performed using Illumina NextSeq. Variants of interest were confirmed by Sanger sequencing or SNP array.

Patients and Setting

Thirty-seven patients in an academic tertiary hospital participated (54% male; median age, 44 years; 86% with fractures), corresponding to 28 sporadic or familial cases.

Main Outcome Measure

The identification of rare stop-gain, indel, splice site, copy-number, or nonsynonymous variants altering protein function.

Results

Altogether, we identified 28 variants of interest, but only 3 were classified as pathogenic or likely pathogenic variants: COL1A2 p.(Arg708Gln), WNT1 p.(Gly169Asp), and IDUA p.(His82Gln). An association of variants in different genes was found in 21% of cases, including a young woman with severe OP bearing WNT1, PLS3, and NOTCH2 variants. Among genes of uncertain significance analyzed, a potential additional line of evidence has arisen for GWAS candidates GPR68 and NBR1, warranting further studies.

Conclusions

While we hope that continuing efforts to identify genetic predisposition to OP will lead to improved and personalized care in the future, the likelihood of identifying actionable pathogenic variants in intriguing cases of idiopathic or familial osteoporosis is seemingly low.

Bone metabolism genes variation and response to bisphosphonate treatment in women with postmenopausal osteoporosis

Introduction

Long-term treatment is used in patients with osteoporosis, and bisphosphonates (BPs) are the most commonly prescribed medications. However, in some patients this therapy is not effective, cause different side effects and complications. Unfortunately, at least one year is needed to identify and confirm an ineffectiveness of BPs therapy on bone mineral density (BMD). Among other factors, a response to BPs therapy may also be explained by genetic factors. The aim of this study was to analyze the influence of SOST, PTH, FGF2, FDPS, GGPS1, and LRP5 gene variants on the response to treatment with aminobisphosphonates.

Materials and methods

Women with postmenopausal osteoporosis were included to this study if they used aminobisphosphonates for at least 12 months. Exclusion criteria were: persistence on BPs therapy less than 80%, bone metabolic diseases, diseases deemed to affect bone metabolism, malignant tumours, using of any medications influencing BMD. The study protocol was approved by the local ethics committee. The BMD at the lumbar spine and femoral neck were measured using dual x-ray absorptiometry (GE Lunar) before and at least 12 months after treatment with BPs. According to BMD change, patients were divided in two groups–responders and non-responders to BPs terapy. Polymorphic variants in SOST, PTH, FGF2, FDPS, GGPS1, and LRP5 genes were determined using PCR analysis with TaqMan probes (Thermo Scientific).

Results

In total, 201 women with BPs therapy were included in the study. No statistically significant differences were observed in age, age at menopause, weight, height, BMI and baseline BMD levels between responders (122 subjects) and non-responders (79 subjects).

As single markers, the SOST rs1234612 T/T (OR = 2.3; P = 0.02), PTH rs7125774 T/T (OR = 2.8, P = 0.0009), FDPS rs2297480 G/G (OR = 29.3, P = 2.2×10⁻⁷), and GGPS1 rs10925503 C/C+C/T (OR = 2.9; P = 0.003) gene variants were over-represented in non-responders group. No significant association between FGF2 rs6854081 and LRP5 rs3736228 gene variants and response to BPs treatment was observed. The carriers of T-T-G-C allelic combination (constructed from rs1234612, rs7125774, rs2297480, and rs10925503) were predisposed to negative response to BPs treatment (OR = 4.9, 95% CI 1.7–14.6, P = 0.005). The C-C-T-C combination was significantly over-represented in responders (OR = 0.1, 95% CI 0.1–0.5, P = 0.006).

Conclusions

Our findings highlight the importance of identified single gene variants and their allelic combinations for pharmacogenetics of BPs therapy of osteoporosis. Complex screening of these genetic markers could be used as a new strategy for personalized antiresorptive therapy.

P4 medicine and osteoporosis: a systematic review

BACKGROUND

Osteoporosis is the most frequent bone metabolic disease. In order to improve early detection, prediction, prevention, diagnosis, and treatment of the disease, a new model of P4 medicine (personalized, predictive, preventive, and participatory medicine) could be applied. The aim of this work was to systematically review the publications of four different types of "omics" studies related to osteoporosis, in order to discover novel predictive, preventive, diagnostic, and therapeutic targets for better management of the geriatric population.

METHODS

To systematically search the PubMed database, we created specific groups of criteria for four different types of "omics" information on osteoporosis: genomic, transcriptomic, proteomic, and metabolomic. We then analyzed the intersections between them in order to find correlations and common pathways or molecules with important roles in osteoporosis, and with a potential application in disease prediction, prevention, diagnosis, or treatment.

RESULTS

Altogether, 180 publications of "omics" studies in the field of osteoporosis were found and reviewed at first selection. After introducing the inclusion and exclusion criteria (the secondary selection), 46 papers were included in the systematic review.

CONCLUSIONS

The intersection of reviewed papers identified five genes (ESR1, IBSP, CTNNB1, SOX4, and IDUA) and processes like the Wnt pathway, JAK/STAT signaling, and ERK/MAPK, which should be further validated for their predictive, diagnostic, or other clinical value in osteoporosis. Such molecular insights will enable us to fit osteoporosis into the P4 strategy and could increase the effectiveness of disease prediction and prevention, with a decrease in morbidity in the geriatric population.

A genome-wide association study meta-analysis of clinical fracture in 10,012 African American women

Background

Osteoporosis is a major public health problem associated with excess disability and mortality. It is estimated that 50–70% of the variation in osteoporotic fracture risk is attributable to genetic factors. The purpose of this hypothesis-generating study was to identify possible genetic determinants of fracture among African American (AA) women in a GWAS meta-analysis.

Methods

Data on clinical fractures (all fractures except fingers, toes, face, skull or sternum) were analyzed among AA female participants in the Women's Health Initiative (WHI) (N = 8155), Cardiovascular Health Study (CHS) (N = 504), BioVU (N = 704), Health ABC (N = 651), and the Johnston County Osteoarthritis Project (JoCoOA) (N = 291). Affymetrix (WHI) and Illumina (Health ABC, JoCoOA, BioVU, CHS) GWAS panels were used for genotyping, and a 1:1 ratio of YRI:CEU HapMap haplotypes was used as an imputation reference panel. We used Cox proportional hazard models or logistic regression to evaluate the association of ~ 2.5 million SNPs with fracture risk, adjusting for ancestry, age, and geographic region where applicable. We conducted a fixed-effects, inverse variance-weighted meta-analysis. Genome-wide significance was set at P < 5 × 10^− 8.

Results

One SNP, rs12775980 in an intron of SVIL on chromosome 10p11.2, reached genome-wide significance (P = 4.0 × 10^− 8). Although this SNP has a low minor allele frequency (0.03), there was no evidence for heterogeneity of effects across the studies (I² = 0). This locus was not reported in any previous osteoporosis-related GWA studies. We also interrogated previously reported GWA-significant loci associated with fracture or bone mineral density in our data. One locus (SMOC1) generalized, but overall there was not substantial evidence of generalization. Possible reasons for the lack of generalization are discussed.

Conclusion

This GWAS meta-analysis of fractures in African American women identified a potentially novel locus in the supervillin gene, which encodes a platelet-associated factor and was previously associated with platelet thrombus formation in African Americans. If validated in other populations of African descent, these findings suggest potential new mechanisms involved in fracture that may be particularly important among African Americans.

•

This was a hypothesis-generating GWAS for fracture in African Americans.

•

One potentially novel locus (SVIL) was identified at GWA-significant levels.

•

SVIL has been associated with platelet thrombus formation in African-Americans.

The genetics of bone mass and susceptibility to bone diseases

Osteoporosis is characterized by low bone mass and an increased risk of fracture. Genetic factors, environmental factors and gene-environment interactions all contribute to a person's lifetime risk of developing an osteoporotic fracture. This Review summarizes key advances in understanding of the genetics of bone traits and their role in osteoporosis. Candidate-gene approaches dominated this field 20 years ago, but clinical and preclinical genetic studies published in the past 5 years generally utilize more-sophisticated and better-powered genome-wide association studies (GWAS). High-throughput DNA sequencing, large genomic databases and improved methods of data analysis have greatly accelerated the gene-discovery process. Linkage analyses of single-gene traits that segregate in families with extreme phenotypes have led to the elucidation of critical pathways controlling bone mass. For example, components of the Wnt-β-catenin signalling pathway have been validated (in both GWAS and functional studies) as contributing to various bone phenotypes. These notable advances in gene discovery suggest that the next decade will witness cataloguing of the hundreds of genes that influence bone mass and osteoporosis, which in turn will provide a roadmap for the development of new drugs that target diseases of low bone mass, including osteoporosis.

Genome-wide association reveals QTL for growth, bone and in vivo carcass traits as assessed by computed tomography in Scottish Blackface lambs

BACKGROUND

Improving meat quality including taste and tenderness is critical to the protection and development of markets for sheep meat. Phenotypic selection for such measures of meat quality is constrained by the fact that these parameters can only be measured post-slaughter. Carcass composition has an impact on meat quality and can be measured on live animals using advanced imaging technologies such as X-ray computed tomography (CT). Since carcass composition traits are heritable, they are potentially amenable to improvement through marker-assisted and genomic selection. We conducted a genome-wide association study (GWAS) on about 600 Scottish Blackface lambs for which detailed carcass composition phenotypes, including bone, fat and muscle components, had been captured using CT and which were genotyped for ~40,000 single nucleotide polymorphisms (SNPs) using the Illumina OvineSNP50 chip.

RESULTS

We confirmed that the carcass composition traits were heritable with moderate to high (0.19-0.78) heritabilities. The GWAS analyses revealed multiple SNPs and quantitative trait loci (QTL) that were associated with effects on carcass composition traits and were significant at the genome-wide level. In particular, we identified a region on ovine chromosome 6 (OAR6) associated with bone weight and bone area that harboured SNPs with p values of 5.55 × 10(-8) and 2.63 × 10(-9), respectively. The same region had effects on fat area, fat density, fat weight and muscle density. We identified plausible positional candidate genes for these OAR6 QTL. We also detected a SNP that reached the genome-wide significance threshold with a p value of 7.28 × 10(-7) and was associated with muscle density on OAR1. Using a regional heritability mapping approach, we also detected regions on OAR3 and 24 that reached genome-wide significance for bone density.

CONCLUSIONS

We identified QTL on OAR1, 3, 24 and particularly on OAR6 that are associated with effects on muscle, fat and bone traits. Based on available evidence that indicates that these traits are genetically correlated with meat quality traits, these associated SNPs have potential applications in selective breeding for improved meat quality. Further research is required to determine whether the effects associated with the OAR6 QTL are caused by a single gene or several closely-linked genes.

Genome-wide Survey of Runs of Homozygosity Identifies Recessive Loci for Bone Mineral Density in Caucasian and Chinese Populations

Runs of homozygosity (ROHs), in which both parental alleles are identical, have been proposed to have recessive effects on multiple human complex diseases. Osteoporosis is a common complex disease characterized by low bone mineral density (BMD), which is highly heritable. And recessive loci that contribute to BMD variations have been identified. In this study, we performed a genome-wide ROHs association study using our SNP array data from three GWAS samples including 4,900 subjects from Caucasian and Chinese populations. Significant results were further subjected to replication in 3,747 additional subjects. ROHs associated with BMD were also tested for associations with osteoporotic fractures in a GWAS fracture sample. Combining results from all the samples, we identified 697 autosomal regions with ROHs. Among these, we detected genome-wide significant associations between BMD and 6 ROHs, including ROH1q31.3, 1p31.1, 3q26.1, 11q12.1, 21q22.1 and 15q22.3 (combined P = 6.29 × 10(-5)-3.17 × 10(-8)). Especially, ROH1p31.1 was found to be associated with increased risk of osteoporotic hip fractures (odds ratio [OR] = 3.71, P = 0.032). To investigate the functional relevance of the identified ROHs, we performed cis-expression quantitative trait locus (eQTL) analysis in lymphoblast cell lines. Three ROHs, including ROH1p31.1, 11q12.1, and 15q22.3, were found to be significantly associated with mRNA expression levels of their nearby genes (PeQTL < 0.05). In summary, our findings reveal that ROHs could play as recessive-acting determinants contributing to the pathogenesis of osteoporosis. Further molecular and functional studies are needed to explore and clarify the potential mechanism.

Parathyroid hormone and its receptor gene polymorphisms: implications in osteoporosis and in fracture healing

Parathyroid glands secrete parathyroid hormone (PTH) which plays multiple roles in calcium homeostasis and in bone remodeling. Secretion of PTH is regulated by extracellular calcium levels and other humoral factors including 1α,25(OH)2D3. PTH regulates gene expression and induces biological effects directly and indirectly. The human gene encoding PTH is located on chromosome 11. In this review, we study the diverse PTH along with its receptor gene polymorphisms and their association with osteoporosis and fracture healing. Genetic factors are associated with osteoporosis by influencing bone mineral density (BMD), bone turnover, calcium homeostasis, and susceptibility to osteoporotic fractures. Polymorphisms in genes encoding PTH may contribute to genetic regulation of BMD and thus susceptibility to fracture risk. PTH stimulates the proliferation of osteoprogenitor cells, production of alkaline phosphatise, and bone matrix proteins that contribute to hard callus formation and increases strength at the site of fractured bone. During remodeling, PTH promotes osteoclastogenesis restoring the original shape, structure, and mechanical strength of the bone. Some PTH polymorphisms have shown an association with fracture risk. Further research is needed to elucidate the relative importance of PTH genetics and the mechanisms of genetic contributions to gene-gene interactions in the pathogenesis of osteoporosis and in fracture healing.

Genetic variants in adult bone mineral density and fracture risk genes are associated with the rate of bone mineral density acquisition in adolescence

Previous studies have identified 63 single-nucleotide polymorphisms (SNPs) associated with bone mineral density (BMD) in adults. These SNPs are thought to reflect variants that influence bone maintenance and/or loss in adults. It is unclear whether they affect the rate of bone acquisition during adolescence. Bone measurements and genetic data were available on 6397 individuals from the Avon Longitudinal Study of Parents and Children at up to five follow-up clinics. Linear mixed effects models with smoothing splines were used for longitudinal modelling of BMD and its components bone mineral content (BMC) and bone area (BA), from 9 to 17 years. Genotype data from the 63 adult BMD associated SNPs were investigated individually and as a genetic risk score in the longitudinal model. Each additional BMD lowering allele of the genetic risk score was associated with lower BMD at age 13 [per allele effect size, 0.002 g/cm(2) (SE = 0.0001, P = 1.24 × 10(-38))] and decreased BMD acquisition from 9 to 17 years (P = 9.17 × 10(-7)). This association was driven by changes in BMC rather than BA. The genetic risk score explained ∼2% of the variation in BMD at 9 and 17 years, a third of that explained in adults (6%). Genetic variants that putatively affect bone maintenance and/or loss in adults appear to have a small influence on the rate of bone acquisition through adolescence.

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.

Gene-gene interaction between RBMS3 and ZNF516 influences bone mineral density

Osteoporosis is characterized by low bone mineral density (BMD), a highly heritable trait that is determined, in part, by the actions and interactions of multiple genes. Although an increasing number of genes have been identified to have independent effects on BMD, few studies have been performed to identify genes that interact with one another to affect BMD. In this study, we performed gene-gene interaction analyses in selected candidate genes in individuals with extremely high versus low hip BMD (20% tails of the distributions), in two independent U.S. Caucasian samples. The first sample contained 916 unrelated subjects with extreme hip BMD Z-scores selected from a population composed of 2286 subjects. The second sample consisted of 400 unrelated subjects with extreme hip BMD Z-scores selected from a population composed of 1000 subjects. Combining results from these two samples, we found one interacting gene pair (RBMS3 versus ZNF516) which, even after Bonferroni correction for multiple testing, showed consistently significant effects on hip BMD. RMBS3 harbored two single-nucleotide polymorphisms (SNPs), rs6549904 and rs7640046, both of which had significant interactions with an SNP, rs4891159, located on ZNF516 (p = 7.04 × 10(-11) and 1.03 × 10(-10) ). We further validated these results in two additional samples of Caucasian and African descent. The gene pair, RBMS3 versus ZNF516, was successfully replicated in the Caucasian sample (p = 8.07 × 10(-3) and 2.91 × 10(-3) ). For the African sample, a significant interaction was also detected (p = 0.031 and 0.043), but the direction of the effect was opposite to that observed in the three Caucasian samples. By providing evidence for genetic interactions underlying BMD, this study further delineates the genetic architecture of osteoporosis.

Polymorphisms in the inflammatory genes CIITA, CLEC16A and IFNG influence BMD, bone loss and fracture in elderly women

Osteoclast activity and the fine balance between bone formation and resorption is affected by inflammatory factors such as cytokines and T lymphocyte activity, mediated by major histocompatibility complex (MHC) molecules, in turn regulated by the MHC class II transactivator (MHC2TA). We investigated the effect of functional polymorphisms in the MHC2TA gene (CIITA), and two additional genes; C-type lectin domain 16A (CLEC16A), in linkage disequilibrium with CIITA and Interferon-γ (IFNG), an inducer of CIITA; on bone density, bone resorption markers, bone loss and fracture risk in 75 year-old women followed for up to 10 years (OPRA n = 1003) and in young adult women (PEAK-25 n = 999). CIITA was associated with BMD at age 75 (lumbar spine p = 0.011; femoral neck (FN) p = 0.049) and age 80 (total body p = 0.015; total hip p = 0.042; FN p = 0.028). Carriers of the CIITA rs3087456(G) allele had 1.8–3.4% higher BMD and displayed increased rate of bone loss between age 75 and 80 (FN p = 0.013; total hip p = 0.030; total body p = 3.8E⁻⁵). Despite increasing bone loss, the rs3087456(G) allele was protective against incident fracture overall (p = 0.002), osteoporotic fracture and hip fracture. Carriers of CLEC16A and IFNG variant alleles had lower BMD (p<0.05) and ultrasound parameters and a lower risk of incident fracture (CLEC16A, p = 0.011). In 25-year old women, none of the genes were associated with BMD. In conclusion, variation in inflammatory genes CIITA, CLEC-16A and INFG appear to contribute to bone phenotypes in elderly women and suggest a role for low-grade inflammation and MHC class II expression for osteoporosis pathogenesis.

Clinical review: Genome-wide association studies of skeletal phenotypes: what we have learned and where we are headed

CONTEXT

The primary goals of genome-wide association studies (GWAS) are to discover new molecular and biological pathways involved in the regulation of bone metabolism that can be leveraged for drug development. In addition, the identified genetic determinants may be used to enhance current risk factor profiles.

EVIDENCE ACQUISITION

There have been more than 40 published GWAS on skeletal phenotypes, predominantly focused on dual-energy x-ray absorptiometry-derived bone mineral density (BMD) of the hip and spine.

EVIDENCE SYNTHESIS

Sixty-six BMD loci have been replicated across all the published GWAS, confirming the highly polygenic nature of BMD variation. Only seven of the 66 previously reported genes (LRP5, SOST, ESR1, TNFRSF11B, TNFRSF11A, TNFSF11, PTH) from candidate gene association studies have been confirmed by GWAS. Among 59 novel BMD GWAS loci that have not been reported by previous candidate gene association studies, some have been shown to be involved in key biological pathways involving the skeleton, particularly Wnt signaling (AXIN1, LRP5, CTNNB1, DKK1, FOXC2, HOXC6, LRP4, MEF2C, PTHLH, RSPO3, SFRP4, TGFBR3, WLS, WNT3, WNT4, WNT5B, WNT16), bone development: ossification (CLCN7, CSF1, MEF2C, MEPE, PKDCC, PTHLH, RUNX2, SOX6, SOX9, SPP1, SP7), mesenchymal-stem-cell differentiation (FAM3C, MEF2C, RUNX2, SOX4, SOX9, SP7), osteoclast differentiation (JAG1, RUNX2), and TGF-signaling (FOXL1, SPTBN1, TGFBR3). There are still 30 BMD GWAS loci without prior molecular or biological evidence of their involvement in skeletal phenotypes. Other skeletal phenotypes that either have been or are being studied include hip geometry, bone ultrasound, quantitative computed tomography, high-resolution peripheral quantitative computed tomography, biochemical markers, and fractures such as vertebral, nonvertebral, hip, and forearm.

CONCLUSIONS

Although several challenges lie ahead as GWAS moves into the next generation, there are prospects of new discoveries in skeletal biology. This review integrates findings from previous GWAS and provides a roadmap for future directions building on current GWAS successes.

Generation and characterization of human anti-human IL-21 neutralizing monoclonal antibodies

Interleukin-21 (IL-21) is a type I four-helical bundle cytokine that exerts a variety of significant effects on many hematopoietic cells, including T and B lymphocytes and natural killer cells. IL-21 is produced predominantly by CD4+ T cells and natural killer T cells and, when aberrantly overexpressed, appears to play important roles in a wide variety of autoimmune disorders. To generate potential therapeutic reagents capable of inhibiting IL-21 for clinical use, we immunized human immunoglobulin transgenic mice with IL-21 and then identified and cloned a panel of human anti-human IL-21 binding monoclonal antibodies. IL-21 neutralizing and IL-21-binding, non-neutralizing antibodies were assigned to distinct epitope "bins" based on surface plasmon resonance competition studies. The most potent neutralizing antibodies had extremely high (sub pM) affinity for IL-21 and were able to block IL-21 activity in various biological assays using either an IL-21R-transfected pre-B-cell line or primary human B cells, and their neutralizing activity was, in some cases, superior to that of a soluble form of the high affinity heterodimeric IL-21 receptor. Characterization of this panel of IL-21 antibodies provided the basis for the selection of a therapeutic candidate antibody capable of inhibiting IL-21 activity for the treatment of autoimmune and inflammatory diseases.

Are bone mineral density loci associated with hip osteoporotic fractures? A validation study on previously reported genome-wide association loci in a Chinese population

Osteoporosis is a heritable disease characterized mainly by low bone mineral density (BMD) and/or osteoporotic fractures (OF). Most genome-wide association studies on osteoporosis have focused on BMD, whereas little effort has been expended to identify genetic variants directly linked to OF. To determine whether BMD-loci are also associated with OF risk, we performed a validation study to examine 23 BMD-loci reported by recent genome-wide association studies for association with hip OF risk. Our sample consisted of 700 elderly Chinese Han subjects, 350 with hip OF and 350 healthy matched controls. We identified four BMD-loci that were significantly associated with hip OF in this Chinese population, including 7q21 (FLJ42280, P = 1.17 × 10(-4) for rs4729260; P = 0.008 for rs7781370), 6p21 (MHC, P = 0.004 for rs3130340), 13q14 (TNFSF11, P = 0.012 for rs9533090; P = 0.018 for rs9594759; P = 0.020 for rs9594738; P = 0.044 for rs9594751), and 18q21 (TNFRSF11A, P = 0.015 for rs884205). The SNP rs4729260 at 7q21 remained significantly associated, even after conservative Bonferroni's correction. Our results further highlight the importance of these loci in the pathogenesis of osteoporosis, and demonstrate that it is feasible and useful to use OF as the direct phenotype to conduct genetic studies, to enhance our understanding of the genetic architecture of osteoporosis.

VKORC1 common variation and bone mineral density in the Third National Health and Nutrition Examination Survey

Osteoporosis, defined by low bone mineral density (BMD), is common among postmenopausal women. The distribution of BMD varies across populations and is shaped by both environmental and genetic factors. Because the candidate gene vitamin K epoxide reductase complex subunit 1 (VKORC1) generates vitamin K quinone, a cofactor for the gamma-carboxylation of bone-related proteins such as osteocalcin, we hypothesized that VKORC1 genetic variants may be associated with BMD and osteoporosis in the general population. To test this hypothesis, we genotyped six VKORC1 SNPs in 7,159 individuals from the Third National Health and Nutrition Examination Survey (NHANES III). NHANES III is a nationally representative sample linked to health and lifestyle variables including BMD, which was measured using dual energy x-ray absorptiometry (DEXA) on four regions of the proximal femur. In adjusted models stratified by race/ethnicity and sex, SNPs rs9923231 and rs9934438 were associated with increased BMD (p = 0.039 and 0.024, respectively) while rs8050894 was associated with decreased BMD (p = 0.016) among non-Hispanic black males (n = 619). VKORC1 rs2884737 was associated with decreased BMD among Mexican-American males (n = 795; p = 0.004). We then tested for associations between VKORC1 SNPs and osteoporosis, but the results did not mirror the associations observed between VKORC1 and BMD, possibly due to small numbers of cases. This is the first report of VKORC1 common genetic variation associated with BMD, and one of the few reports available that investigate the genetics of BMD and osteoporosis in diverse populations.

Genetic epidemiology of age-related osteoporosis and its clinical applications

Osteoporosis is an important and complex disorder that is highly prevalent worldwide. This disease poses a major challenge to modern medicine and its treatment is associated with high costs. Numerous studies have endeavored to decipher the pathogenesis of this disease. The clinical assessment of patients often incorporates information about a family history of osteoporotic fractures. Indeed, the observation of an increased risk of fracture in an individual with a positive parental history of hip fracture provides strong evidence for the heritability of osteoporosis. The onset and progression of osteoporosis are generally controlled by multiple genetic and environmental factors, as well as interactions between them, with rare cases determined by a single gene. In an attempt to identify the genetic markers of complex diseases such as osteoporosis, there has been a move away from traditional linkage mapping studies and candidate gene association studies to higher-density genome-wide association studies. The advent of high-throughput technology enables genotyping of millions of DNA markers in the human genome, and consequently the identification and characterization of causal variants and loci that underlie osteoporosis. This Review presents an overview of the major findings since 2007 and clinical applications of these genome-wide linkage and association studies.

Molecular genetic studies of gene identification for osteoporosis: the 2009 update

Osteoporosis is a complex human disease that results in increased susceptibility to fragility fractures. It can be phenotypically characterized using several traits, including bone mineral density, bone size, bone strength, and bone turnover markers. The identification of gene variants that contribute to osteoporosis phenotypes, or responses to therapy, can eventually help individualize the prognosis, treatment, and prevention of fractures and their adverse outcomes. Our previously published reviews have comprehensively summarized the progress of molecular genetic studies of gene identification for osteoporosis and have covered the data available to the end of September 2007. This review represents our continuing efforts to summarize the important and representative findings published between October 2007 and November 2009. The topics covered include genetic association and linkage studies in humans, transgenic and knockout mouse models, as well as gene-expression microarray and proteomics studies. Major results are tabulated for comparison and ease of reference. Comments are made on the notable findings and representative studies for their potential influence and implications on our present understanding of the genetics of osteoporosis.

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.