Causal inference in health and disease: a review of the principles and applications of Mendelian randomization.

Mendelian randomization (MR) is a genetic epidemiological technique that uses genetic variation to infer causal relationships between modifiable exposures and outcome variables. Conventional observational epidemiological studies are subject to bias from a range of sources; MR analyses can offer an advantage in that they are less prone to bias as they use genetic variants inherited at conception as "instrumental variables" which are proxies of an exposure. However, as with all research tools, MR studies must be carefully designed to yield valuable insights into causal relationships between exposures and outcomes, and to avoid biased or misleading results that undermine the validity of the causal inferences drawn from the study. In this review we outline Mendel's laws of inheritance, the assumptions and principles that underlie MR, MR study designs and methods, and how MR analyses can be applied, and reported. Using the example of serum phosphate concentrations on liability to kidney stone disease we illustrate how MR estimates may be visualized and, finally, we contextualize MR in bone and mineral research including exemplifying how this technique could be employed to inform clinical studies and future guidelines concerning bone mineral density and fracture risk. This review provides a framework to enhance understanding of how MR may be used to triangulate evidence and progress research in bone and mineral metabolism as we strive to infer causal effects in health and disease.

Mendelian randomization is a powerful tool used by researchers to understand how factors like lifestyle choices or medical conditions affect our health. Some research studies can be hampered by inaccurate information or measurements which mean the results are unreliable. Mendelian randomization uses genetic information to evade this problem. Because a person's genetic information does not change the studies can provide more reliable conclusions about how different factors influence our health. This review provides an overview of Mendelian randomization principles, its applications, and how researchers in bone and mineral research have used it to uncover cause-and-effect relationships.

Enhanced fracture risk prediction: a novel multi-trait genetic approach integrating polygenic scores of fracture-related traits.

The novel metaPGS, integrating multiple fracture-related genetic traits, surpasses traditional polygenic scores in predicting fracture risk. Demonstrating a robust association with incident fractures, this metaPGS offers significant potential for enhancing clinical fracture risk assessment and tailoring prevention strategies. INTRODUCTION: Current polygenic scores (PGS) have limited predictive power for fracture risk. To improve genetic prediction, we developed and evaluated a novel metaPGS combining genetic information from multiple fracture-related traits. METHODS: We derived individual PGS from genome-wide association studies of 16 fracture-related traits and employed an elastic-net logistic regression model to examine the association between the 16 PGSs and fractures. An optimal metaPGS was constructed by combining 11 significant individual PGSs selected by the elastic regularized regression model. We evaluated the predictive power of the metaPGS alone and in combination with clinical risk factors recommended by guidelines. The discrimination ability of metaPGS was assessed using the concordance index. Reclassification was assessed using net reclassification improvement (NRI) and integrated discrimination improvement (IDI). RESULTS: The metaPGS had a significant association with incident fractures (HR 1.21, 95% CI 1.18-1.25 per standard deviation of metaPGS), which was stronger than previously developed bone mineral density (BMD)-related individual PGSs. Models with PGS_FNBMD, PGS_TBBMD, and metaPGS had slightly higher but statistically non-significant c-index than the base model (0.640, 0.644, 0.644 vs. 0.638). However, the reclassification analysis showed that compared to the base model, the model with metaPGS improves the reclassification of fracture. CONCLUSIONS: The metaPGS is a promising approach for stratifying fracture risk in the European population, improving fracture risk prediction by combining genetic information from multiple fracture-related traits.

Murine IRF8 Mutation Offers New Insight into Osteoclast and Root Resorption.

Interferon regulatory factor 8 (IRF8), a transcription factor expressed in immune cells, functions as a negative regulator of osteoclasts and helps maintain dental and skeletal homeostasis. Previously, we reported that a novel mutation in the IRF8 gene increases susceptibility to multiple idiopathic cervical root resorption (MICRR), a form of tooth root resorption mediated by increased osteoclast activity. The IRF8 G388S variant in the highly conserved C-terminal motif is predicted to alter the protein structure, likely impairing IRF8 function. To investigate the molecular basis of MICRR and IRF8 function in osteoclastogenesis, we generated Irf8 knock-in (KI) mice using CRISPR/Cas9 technique modeling the human IRF8G388S mutation. The heterozygous (Het) and homozygous (Homo) Irf8 KI mice showed no gross morphological defects, and the development of hematopoietic cells was unaffected and similar to wild-type (WT) mice. The Irf8 KI Het and Homo mice showed no difference in macrophage gene signatures important for antimicrobial defenses and inflammatory cytokine production. Consistent with the phenotype observed in MICRR patients, Irf8 KI Het and Homo mice demonstrated significantly increased osteoclast formation and resorption activity in vivo and in vitro when compared to WT mice. The oral ligature-inserted Het and Homo mice displayed significantly increased root resorption and osteoclast-mediated alveolar bone loss compared to WT mice. The increased osteoclastogenesis noted in KI mice is due to the inability of IRF8G388S mutation to inhibit NFATc1-dependent transcriptional activation and downstream osteoclast specific transcripts, as well as its impact on autophagy-related pathways of osteoclast differentiation. This translational study delineates the IRF8 domain important for osteoclast function and provides novel insights into the IRF8 mutation associated with MICRR. IRF8G388S mutation mainly affects osteoclastogenesis while sparing immune cell development and function. These insights extend beyond oral health and significantly advance our understanding of skeletal disorders mediated by increased osteoclast activity and IRF8's role in osteoclastogenesis.

The Effect of Asfotase Alfa on Plasma and Urine Pyrophosphate Levels and Pseudofractures in a Patient With Adult-Onset Hypophosphatasia.

Hypophosphatasia (HPP) is an inherited disease caused by variants of the ALPL gene encoding tissue-nonspecific alkaline phosphatase. Adult-onset HPP (adult HPP), known as a mild form of HPP, develops symptoms involving osteomalacia after the age of 18 years. Asfotase alfa (AA) is a modulated recombinant human alkaline phosphatase (ALP) that has been established as a first-line therapy for severe forms of HPP, such as perinatal and infantile forms. We described a 64-year-old female who presented with pseudofractures in bilateral femur diaphyses and impaired mobility. Low serum ALP activity and a high concentration of urine phosphoethanolamine indicated the diagnosis of HPP, which was confirmed by the identification of a homozygous variant in the ALPL gene (c.319G > A; p.Val107Ile). An in vitro transfection experiment to measure the ALP activity of this novel variant protein was performed, resulting in 40% of the residual enzymatic activity compared with the wild type. AA was initiated to facilitate the union of pseudofracture and to improve mobility. After 6 months, radiographic images revealed the disappearance of fracture lines, and improvement of ambulatory ability was confirmed by the 6-minute walk test (525 to 606 m). The EQ-5D-5L index was also improved (0.757 to 0.895). Within a follow-up period, the levels of urine pyrophosphate corrected by urine creatinine (uPPi/Cre) declined in parallel with the level of plasma PPi (plasma PPi: 6.34 to 1.04 µM, uPPi/Cre: 226.8 to 75.4 nmol/mg). The beneficial effect of AA on pseudofracture healing in adult HPP was presented, although the application of AA should be restricted to patients exhibiting relatively severe manifestations. In addition, a novel pathogenic variant of the ALPL gene was identified with the supportive result of functional analysis. Furthermore, when monitoring patients with HPP treated with AA, uPPi/Cre might be a convenient substitute for plasma PPi, which requires immediate filtration after blood sampling. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC. on behalf of American Society for Bone and Mineral Research.

Circulating lipocalin-2 and features of metabolic syndrome in community-dwelling older women: A cross-sectional study.

Lipocalin-2 (LCN2) is released by several cell types including osteoblasts and adipocytes and has been suggested as a marker of renal dysfunction, metabolic syndrome (MetS) and type 2 diabetes (T2D). Whether LCN2 is linked to these diseases in older women remains unknown. This study investigated whether LCN2 is related to features of MetS and T2D in older women. This cross-sectional study included 705 non-diabetic women (mean age 75.1 ± 2.6 years) for MetS analysis and 76 women (mean age 75.4 ± 2.8 years) with T2D. Total circulating LCN2 levels were analysed using a two-step chemiluminescent microparticle monoclonal immunoassay. MetS was determined by a modified National Cholesterol Education Program Adult Treatment Panel III classification. Multivariable-adjusted logistic regression analysis was used to assess odds ratios between LCN2 quartiles and MetS. Women in the highest LCN2 quartile had approximately 3 times greater risk for MetS compared to women in the lowest quartile (OR 3.05; 95%CI 1.86-5.02). Women with T2D or MetS scores of ≥ 3 had higher LCN2 levels compared to women with a MetS score of 0 (p < 0.05). Higher LCN2 correlated with higher body mass index, fat mass, triglycerides and glycated haemoglobin and lower high-density lipoprotein cholesterol and estimated glomerular filtration rate (p < 0.05). Higher circulating levels of LCN2 are associated with worsened cardio-metabolic risk factors and increased odds of MetS and T2D in older women. Whether it can be used as a biomarker for identifying those at risk for MetS and T2D should be explored further.

Ethnic disparities in fracture risk assessment using polygenic scores.

Whether the PGS developed using data from European ancestry is predictive of fracture risk for minorities remains unclear. This study demonstrated that PGSs based on common BMD-related genetic variants discovered in the European ancestry cohort are predictive of fracture risk in people of Asian but not African ancestry. PURPOSE: Large-scale genome-wide association studies (GWAS) on bone mineral density (BMD) have been conducted predominantly in European cohorts. Genetic models based on common variants associated with BMD have been evaluated using almost exclusively European data, which could potentially exacerbate health disparities due to different linkage disequilibrium among different ethnic groups. METHODS: UK Biobank (UKB) is a large-scale population-based observational study starting in 2006 that recruited 502,617 individuals aged between 40 and 69 years with genotypic and phenotypic data available. Based on the summary statistics of two GWAS studies of femoral neck BMD and total body BMD, we derived four PGSs and assessed the association between each PGS and prevalent/incident fractures within each ethnic group separately using Multivariate logistic regressions and Cox proportional hazard models. All models were adjusted for age, sex, and the first four principal components. RESULTS: We assessed four PGSs derived from European cohorts. Significant associations were observed between PGSs and fracture in European and Asian cohorts but not in the African cohort. Of all four PGSs, [Formula: see text] performed the best. A standard deviation decreases in [Formula: see text] were associated with an increased hazard ratio (HR) of 1.24 (1.22-1.27), 1.28 (0.83-1.99), and 1.34 (1.10-1.64) in European, African, and Asian ancestry, respectively. A low BMD-related PGS is associated with up to 2.35- and 4.31-fold increased fracture risk in European and Asian populations. CONCLUSIONS: These results showed that PGSs based on common BMD-related genetic variants discovered in the European ancestry cohort are predictive of fracture risk in people of Asian but not African ancestry.

The clinical utility of the BMD-related comprehensive genome-wide polygenic score in identifying individuals with a high risk of osteoporotic fractures.

The potential of bone mineral density (BMD)-related genome-wide polygenic score (PGS) in identifying individuals with a high risk of fractures remains unclear. This study suggests that an efficient PGS enables the identification of strata with up to a 1.5-fold difference in fracture incidence. Incorporating PGS into clinical diagnosis is anticipated to increase the population-level screening benefits. PURPOSE: This study sought to construct genome-wide polygenic scores for femoral neck and total body BMD and to estimate their potential in identifying individuals with a high risk of osteoporotic fractures. METHODS: Genome-wide polygenic scores were developed and validated for femoral neck and total body BMD. We externally tested the PGSs, both by themselves and in combination with available clinical risk factors, in 455,663 European ancestry individuals from the UK Biobank. The predictive accuracy of the developed genome-wide PGS was also compared with previously published restricted PGS employed in fracture risk assessment. RESULTS: For each unit decrease in PGSs, the genome-wide PGSs were associated with up to 1.17-fold increased fracture risk. Out of four studied PGSs, [Formula: see text] (HR: 1.03; 95%CI 1.01-1.05, p = 0.001) had the weakest and the [Formula: see text] (HR: 1.17; 95%CI 1.15-1.19, p < 0.0001) had the strongest association with an incident fracture. In the reclassification analysis, compared to the FRAX base model, the models with [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text] improved the reclassification of fracture by 1.2% (95% CI, 1.0 to 1.3%), 0.2% (95% CI, 0.1 to 0.3%), 1.4% (95% CI, 1.3 to 1.5%), and 2.2% (95% CI, 2.1 to 2.4%), respectively. CONCLUSIONS: Our findings suggested that an efficient PGS estimate enables the identification of strata with up to a 1.7-fold difference in fracture incidence. Incorporating PGS information into clinical diagnosis is anticipated to increase the benefits of screening programs at the population level.

Identification of a Functional Susceptibility Variant for Adolescent Idiopathic Scoliosis that Upregulates Early Growth Response 1 (EGR1)-Mediated UNCX Expression.

Adolescent idiopathic scoliosis (AIS) is a serious health problem affecting 3% of live births all over the world. Many loci associated with AIS have been identified by previous genome wide association studies, but their biological implication remains mostly unclear. In this study, we evaluated the AIS-associated variants in the 7p22.3 locus by combining in silico, in vitro, and in vivo analyses. rs78148157 was located in an enhancer of UNCX, a homeobox gene and its risk allele upregulated the UNCX expression. A transcription factor, early growth response 1 (EGR1), transactivated the rs78148157-located enhancer and showed a higher binding affinity for the risk allele of rs78148157. Furthermore, zebrafish larvae with UNCX messenger RNA (mRNA) injection developed body curvature and defective neurogenesis in a dose-dependent manner. rs78148157 confers the genetic susceptibility to AIS by enhancing the EGR1-regulated UNCX expression. © 2022 American Society for Bone and Mineral Research (ASBMR).

Comparisons Between Different Anti-osteoporosis Medications on Postfracture Mortality: A Population-Based Study.

CONTEXT: Osteoporosis is becoming a global epidemic in aging societies. Anti-osteoporotic medications can prevent fractures, and their pleiotropic effect on mortality is interesting but not well compared among each other. OBJECTIVE: To provide real-world evidence on the pleiotropic effect of different anti-osteoporotic medications on all-cause mortality, stratified by fracture site, sex, and age. METHODS: This longitudinal population-based postfracture cohort study, included mega-data from subjects ≥40 years of age with osteoporotic fracture who used anti-osteoporotic medications as recorded in Taiwan's National Health Insurance Research Database from 2009 to 2017 and followed until 2018. A multivariate Cox proportional hazards model with immortal time bias was used to assess the relationship between fracture sites and mortality stratified by anti-osteoporosis medication. RESULTS: A total of 46 729 subjects with an average age of 74.45 years (80.0% female) and a mean follow-up period of 4.73 years were enrolled. In the total fracture group, compared with raloxifene and bazedoxifene, we found that alendronate/risedronate (hazard ratio [HR] 0.83; 95% CI, 0.79-0.88), denosumab (HR 0.86; 95% CI, 0.81-0.91), and zoledronic acid (HR 0.78; 95% CI, 0.73-0.84) resulted in significantly lower mortality. Similar trends were observed in the hip, vertebral, or nonhip/nonvertebral fracture groups. Subjects receiving long-acting zoledronic acid showed the lowest mortality in the subanalysis according to sex or age over 65 years. CONCLUSION: This real-world mega-data study suggests that the usage of osteoporotic medication, especially a long-acting regimen, may lower postfracture mortality.

The Effects of Acute High-Intensity Interval Exercise and Hyperinsulinemic-Euglycemic Clamp on Osteoglycin Levels in Young and Middle-Aged Men.

Osteoglycin (OGN) is a leucine-rich proteoglycan that has been implicated in the regulation of glucose in animal models. However, its relationship with glucose control in humans is unclear. We examined the effect of high-intensity interval exercise (HIIE) and hyperinsulinemic-euglycemic clamp on circulating levels of OGN as well as whether circulating OGN levels are associated with markers of glycemic control and cardio-metabolic health. Serum was analyzed for OGN (ELISA) levels from 9 middle-aged obese men (58.1 ± 2.2 years, body mass index [BMI] = 33.1 ± 1.4 kg∙m-2, mean ± SEM) and 9 young men (27.8 ± 1.6 years, BMI = 24.4 ± 0.08 kg∙m-2) who previously completed a study involving a euglycemic-hyperinsulinemic clamp at rest and after HIIE (4x4 minutes cycling at approximately 95% peak heart rate (HRpeak), interspersed with 2 minutes of active recovery). Blood pressure, body composition (dual-energy X-ray absorptiometry), and insulin sensitivity (hyperinsulinemic-euglycemic clamp) were assessed. Serum OGN was higher in the young cohort compared with the middle-aged cohort (65.2 ± 10.1 ng/mL versus 36.5 ± 4. 5 ng/mL, p ≤ 0.05). Serum OGN was unaffected by acute HIIE but decreased after the insulin clamp compared with baseline (~-27%, p = 0.01), post-exercise (~-35%, p = 0.01), and pre-clamp (~-32%, p = 0.02) time points, irrespective of age. At baseline, lower circulating OGN levels were associated with increased age, BMI, and fat mass, whereas higher OGN levels were related to lower fasting glucose. Higher OGN levels were associated with a higher glucose infusion rate. Exercise had a limited effect on circulating OGN. The mechanisms by which OGN affects glucose regulation should be explored in the future. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Identification of Known and Novel Long Noncoding RNAs Potentially Responsible for the Effects of Bone Mineral Density (BMD) Genomewide Association Study (GWAS) Loci.

Osteoporosis, characterized by low bone mineral density (BMD), is the most common complex disease affecting bone and constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 associations influencing BMD. It has been shown that perturbations to long noncoding RNAs (lncRNAs) influence BMD and the activities of bone cells; however, the extent to which lncRNAs are involved in the genetic regulation of BMD is unknown. Here, we combined the analysis of allelic imbalance (AI) in human acetabular bone fragments with a transcriptome-wide association study (TWAS) and expression quantitative trait loci (eQTL) colocalization analysis using data from the Genotype-Tissue Expression (GTEx) project to identify lncRNAs potentially responsible for GWAS associations. We identified 27 lncRNAs in bone that are located in proximity to a BMD GWAS association and harbor single-nucleotide polymorphisms (SNPs) demonstrating AI. Using GTEx data we identified an additional 31 lncRNAs whose expression was associated (false discovery rate [FDR] correction < 0.05) with BMD through TWAS and had a colocalizing eQTL (regional colocalization probability [RCP] > 0.1). The 58 lncRNAs are located in 43 BMD associations. To further support a causal role for the identified lncRNAs, we show that 23 of the 58 lncRNAs are differentially expressed as a function of osteoblast differentiation. Our approach identifies lncRNAs that are potentially responsible for BMD GWAS associations and suggest that lncRNAs play a role in the genetics of osteoporosis. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

An Osteoporosis Susceptibility Allele at 11p15 Regulates SOX6 Expression by Modulating TCF4 Chromatin Binding.

Osteoporosis is an age-related complex disease clinically diagnosed with bone mineral density (BMD). Although several genomewide association studies (GWASs) have discovered multiple noncoding genetic variants at 11p15 influencing osteoporosis risk, the functional mechanisms of these variants remain unknown. Through integrating bioinformatics and functional experiments, a potential functional single-nucleotide polymorphism (SNP; rs1440702) located in an enhancer element was identified and the A allele of rs1440702 acted as an allelic specificities enhancer to increase its distal target gene SOX6 (~600 Kb upstream) expression, which plays a key role in bone formation. We also validated this long-range regulation via conducting chromosome conformation capture (3C) assay. Furthermore, we demonstrated that SNP rs1440702 with a risk allele (rs1440702-A) could increase the activity of the enhancer element by altering the binding affinity of the transcription factor TCF4, resulting in the upregulation expression of SOX6 gene. Collectively, our integrated analyses revealed how the noncoding genetic variants (rs1440702) affect osteoporosis predisposition via long-range gene regulatory mechanisms and identified its target gene SOX6 for downstream biomarker and drug development. © 2022 American Society for Bone and Mineral Research (ASBMR).

Molecular Diagnoses of X-Linked and Other Genetic Hypophosphatemias: Results From a Sponsored Genetic Testing Program.

X-linked hypophosphatemia (XLH), a dominant disorder caused by pathogenic variants in the PHEX gene, affects both sexes of all ages and results in elevated serum fibroblast growth factor 23 (FGF23) and below-normal serum phosphate. In XLH, rickets, osteomalacia, short stature, and lower limb deformity may be present with muscle pain and/or weakness/fatigue, bone pain, joint pain/stiffness, hearing difficulty, enthesopathy, osteoarthritis, and dental abscesses. Invitae and Ultragenyx collaborated to provide a no-charge sponsored testing program using a 13-gene next-generation sequencing panel to confirm clinical XLH or aid diagnosis of suspected XLH/other genetic hypophosphatemia. Individuals aged ≥6 months with clinical XLH or suspected genetic hypophosphatemia were eligible. Of 831 unrelated individuals tested between February 2019 and June 2020 in this cross-sectional study, 519 (62.5%) individuals had a pathogenic or likely pathogenic variant in PHEX (PHEX-positive). Among the 312 PHEX-negative individuals, 38 received molecular diagnoses in other genes, including ALPL, CYP27B1, ENPP1, and FGF23; the remaining 274 did not have a molecular diagnosis. Among 319 patients with a provider-reported clinical diagnosis of XLH, 88.7% (n = 283) had a reportable PHEX variant; 81.5% (n = 260) were PHEX-positive. The most common variant among PHEX-positive individuals was an allele with both the gain of exons 13-15 and c.*231A>G (3'UTR variant) (n = 66/519). Importantly, over 80% of copy number variants would have been missed by traditional microarray analysis. A positive molecular diagnosis in 41 probands (4.9%; 29 PHEX positive, 12 non-PHEX positive) resulted in at least one family member receiving family testing. Additional clinical or family member information resulted in variant(s) of uncertain significance (VUS) reclassification to pathogenic/likely pathogenic (P/LP) in 48 individuals, highlighting the importance of segregation and clinical data. In one of the largest XLH genetic studies to date, 65 novel PHEX variants were identified and a high XLH diagnostic yield demonstrated broad insight into the genetic basis of XLH. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Osteopenia is associated with wasting in pancreatic adenocarcinoma and predicts survival after surgery.

Pancreatic ductal adenocarcinoma (PDAC) is among the deadliest of all common malignancies. Treatment is difficult and often complicated by the presence of cachexia. The clinical portrait of cachexia contributes to the poor prognosis experienced by PDAC patients and worsens therapeutic outcomes. We propose that low bone mineral density is a component of cachexia, which we explore herein through a retrospective review of all patients at our facility that underwent surgery for PDAC between 2011 and 2018 and compared to sex-, age- and comorbidity-matched control individuals. Data were abstracted from the medical record and pre-operative computed tomography scans. Muscle mass and quality were measured at the L3 level and bone mineral density was measured as the radiation attenuation of the lumbar vertebral bodies. Patients with PDAC displayed typical signs of cachexia such as weight loss and radiologically appreciable deterioration of skeletal muscle. Critically, PDAC patients had significantly lower bone mineral density than controls, with 61.2% of PDAC patients categorized as osteopenic compared to 36.8% of controls. PDAC patients classified as osteopenic had significantly reduced survival (1.01 years) compared to patients without osteopenia (2.77 years). The presence of osteopenia was the strongest clinical predictor of 1- and 2-year disease-specific mortality, increasing the risk of death by 107% and 80%, respectively. Osteopenia serves as a test of 2-year mortality with sensitivity of 76% and specificity of 58%. These data therefore identify impaired bone mineral density as a key component of cachexia and predictor of postoperative survival in patients with PDAC. The mechanisms that lead to bone wasting in tumor-bearing hosts deserve further study.

Polygenic Risk Score of Adolescent Idiopathic Scoliosis for Potential Clinical Use.

Adolescent idiopathic scoliosis (AIS) is a common disease causing three-dimensional spinal deformity in as many as 3% of adolescents. Development of a method that can accurately predict the onset and progression of AIS is an immediate need for clinical practice. Because the heritability of AIS is estimated as high as 87.5% in twin studies, prediction of its onset and progression based on genetic data is a promising option. We show the usefulness of polygenic risk score (PRS) for the prediction of onset and progression of AIS. We used AIS genomewide association study (GWAS) data comprising 79,211 subjects in three cohorts and constructed a PRS based on association statistics in a discovery set including 31,999 female subjects. After calibration using a validation data set, we applied the PRS to a test data set. By integrating functional annotations showing heritability enrichment in the selection of variants, the PRS demonstrated an association with AIS susceptibility (p = 3.5 × 10-40 with area under the receiver-operating characteristic [AUROC] = 0.674, sensitivity = 0.644, and specificity = 0.622). The decile with the highest PRS showed an odds ratio of as high as 3.36 (p = 1.4 × 10-10 ) to develop AIS compared with the fifth in decile. The addition of a predictive model with only a single clinical parameter (body mass index) improved predictive ability for development of AIS (AUROC = 0.722, net reclassification improvement [NRI] 0.505 ± 0.054, p = 1.6 × 10-8 ), potentiating clinical use of the prediction model. Furthermore, we found the Cobb angle (CA), the severity measurement of AIS, to be a polygenic trait that showed a significant genetic correlation with AIS susceptibility (rg = 0.6, p = 3.0 × 10-4 ). The AIS PRS demonstrated a significant association with CA. These results indicate a shared polygenic architecture between onset and progression of AIS and the potential usefulness of PRS in clinical settings as a predictor to promote early intervention of AIS and avoid invasive surgery. © 2021 American Society for Bone and Mineral Research (ASBMR).

A Meta-Analysis of the Transferability of Bone Mineral Density Genetic Loci Associations From European to African Ancestry Populations.

Genetic studies of bone mineral density (BMD) largely have been conducted in European populations. We therefore conducted a meta-analysis of six independent African ancestry cohorts to determine whether previously reported BMD loci identified in European populations were transferable to African ancestry populations. We included nearly 5000 individuals with both genetic data and assessments of BMD. Genotype imputation was conducted using the 1000G reference panel. We assessed single-nucleotide polymorphism (SNP) associations with femoral neck and lumbar spine BMD in each cohort separately, then combined results in fixed effects (or random effects if study heterogeneity was high, I2 index >60) inverse variance weighted meta-analyses. In secondary analyses, we conducted locus-based analyses of rare variants using SKAT-O. Mean age ranged from 12 to 68 years. One cohort included only men and another cohort included only women; the proportion of women in the other four cohorts ranged from 52% to 63%. Of 56 BMD loci tested, one locus, 6q25 (C6orf97, p = 8.87 × 10-4 ), was associated with lumbar spine BMD and two loci, 7q21 (SLC25A13, p = 2.84 × 10-4 ) and 7q31 (WNT16, p = 2.96 × 10-5 ), were associated with femoral neck BMD. Effects were in the same direction as previously reported in European ancestry studies and met a Bonferroni-adjusted p value threshold, the criteria for transferability to African ancestry populations. We also found associations that met locus-specific Bonferroni-adjusted p value thresholds in 11q13 (LRP5, p < 2.23 × 10-4 ), 11q14 (DCDC5, p < 5.35 × 10-5 ), and 17p13 (SMG6, p < 6.78 × 10-5 ) that were not tagged by European ancestry index SNPs. Rare single-nucleotide variants in AKAP11 (p = 2.32 × 10-2 ), MBL2 (p = 4.09 × 10-2 ), MEPE (p = 3.15 × 10-2 ), SLC25A13 (p = 3.03 × 10-2 ), STARD3NL (p = 3.35 × 10-2 ), and TNFRSF11A (p = 3.18 × 10-3 ) were also associated with BMD. The majority of known BMD loci were not transferable. Larger genetic studies of BMD in African ancestry populations will be needed to overcome limitations in statistical power and to identify both other loci that are transferable across populations and novel population-specific variants. © 2020 American Society for Bone and Mineral Research (ASBMR).

Post-GWAS Polygenic Risk Score: Utility and Challenges.

Over the past decade, through genome-wide association studies, more than 300 genetic variants have been identified to be associated with either BMD or fracture risk. These genetic variants are common in the general population, but they exert small to modest effects on BMD, suggesting that the utility of any single variant is limited. However, a combination of effect sizes from multiple variants in the form of the polygenic risk score (PRS) can provide a useful indicator of fracture risk beyond that obtained by conventional clinical risk factors. In this perspective, we review the progress of genetics of osteoporosis and approaches for creating PRSs, their uses, and caveats. Recent studies support the idea that the PRS, when integrated into existing fracture prediction models, can help clinicians and patients alike to better assess the fracture risk for an individual, and raise the possibility of precision risk assessment. © 2020 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

The Membrane Transporter OAT7 (SLC22A9) Is Not a Susceptibility Factor for Osteoporosis in Europeans.

Bone production, maintenance, and modeling are a well-balanced process involving mineralization by osteoblasts and resorption by osteoclasts. Sex steroid hormones, including their conjugated forms, contribute majorly to maintaining this balance. Recently, variants in the SLC22A9 gene have been associated with osteoporosis in Korean females. We had recently shown that SLC22A9, encoding organic anion transporter 7 (OAT7), is an uptake transporter of estrone sulfate and identified several genetic variants in Europeans leading to functional consequences in vitro. We therefore hypothesized that SLC22A9 genetic variants may contribute to the pathophysiology of osteoporosis in Europeans. To test this hypothesis, we examined the associations of SLC22A9 variants with bone quality, fractures, and bone turnover markers. We genotyped SLC22A9 variants in 5,701 (2,930 female) subjects (age range, 20-93 years) extracted from the population-based Study of Health in Pomerania (SHIP and SHIP-TREND) covered by the Illumina Infinium HumanExome BeadChip version v1.0 (Exome Chip). Descriptive data (e.g., history of fractures), ultrasonography of the calcaneus, as well as serum concentrations of carboxy-terminal telopeptide of type I collagen, amino-terminal propeptide of type I procollagen, and vitamin D were determined. Comprehensive statistical analyses revealed no association between low-frequency and rare SLC22A9 variants and bone quality, fractures, and bone turnover markers. Our results indicate that single genetic SLC22A9 variants do not have a major impact on osteoporosis risk prediction in Europeans, yet findings need to be replicated in larger-scale studies.

Genetic Burden Contributing to Extremely Low or High Bone Mineral Density in a Senior Male Population From the Osteoporotic Fractures in Men Study (MrOS).

Worldwide, one in five men aged over 50 years will experience osteoporosis or a clinical bone fracture, with a greater fracture-related mortality rate than women. However, the genetic etiology of osteoporosis in men is still poorly understood. We aimed to identify the genetic variants and candidate genes associated with extremely low or high BMD for a better understanding of the biology underlying low bone density that may point to potential therapeutic targets for increasing bone mass. Subjects from the Osteoporotic Fractures in Men Study (MrOS) cohort were evaluated by age and BMI-adjusted total hip BMD. Those with BMD values 3 SDs away from the mean were selected and the remaining individuals whose adjusted BMD ranked at the highest or lowest 100 were included. Men with the lowest adjusted BMD (N = 98) and highest adjusted BMD (N = 110) were chosen for exome sequencing. Controls (N = 82) were men of Northern and Western European descent from the US Utah population of the 1000 Genomes Project. Fisher's exact test was performed to compare low- or high-BMD subjects with controls for single-gene associations. Additionally, sets of candidate genes causative of heritable disorders of connective tissue, including osteogenesis imperfecta (OI) and Ehlers-Danlos syndrome (EDS), were grouped for multigene and mutation burden analyses. No single-gene associations with rare variants were found for either the low BMD group (33 genes) or high BMD group (18 genes). In the group of OI genes, we detected a significant threefold increased accumulation of rare variants in low-BMD subjects compared with controls (p = 0.009). Additionally, genes associated with EDS had a twofold increased frequency in low-BMD subjects compared with controls (p = 0.03). These findings reveal a rare variant burden in OI and EDS disease genes at low BMD, which suggests a potential gene-panel approach to screen for multivariant associations in larger cohorts. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.