MEPE loss-of-function variant associates with decreased bone mineral density and increased fracture risk

Osteosarcopenia, osteoarthritis and frailty: a two-sample Mendelian randomization study

BACKGROUND

Musculoskeletal disease, which has a complicated relationship with frailty, is a common clinical problem among elderly individuals.

AIMS

This study evaluated the potential causal relationships between osteosarcopenia, osteoarthritis and frailty by Mendelian Randomization (MR) analysis.

METHODS

This study employed a two-sample MR approach to investigate the causal relationships among osteosarcopenia, osteoarthritis and frailty. Published summary statistics were used to obtain instrumental variables at the genome-wide significance level.

RESULTS

Among the age groups with osteoporosis, high total bone mineral density (TBMD) (45-60, OR = 0.966, 95% CI 0.940-0.993, P = 0.013) and TBMD (over 60, OR = 0.974, 95% CI 0.954-0.994, P = 0.011) reduced the risk of frailty. Similarly, high forearm BMD (FA-BMD), high ultradistal forearm BMD (UFA-BMD), and high Heel-BMD at different sites also reduced the risk of frailty (OR = 0.966, 95% CI 0.936-0.996, P = 0.028; OR = 0.975, 95% CI 0.953-0.997, P = 0.029; OR = 0.981, 95% CI 0.967-0.995, P = 0.008). Among the characteristics related to sarcopenia, grip strength in the left hand, grip strength in the right hand, appendicular lean mass, and walking pace were all protective factors for frailty (OR = 0.788, 95% CI 0.721-0.862, P < 0.001; OR = 0.800, 95% CI 0.737-0.869, P < 0.001; OR = 0.955, 95% CI 0.937-0.974, P = 0.000; OR = 0.480, 95% CI 0.388-0.593, P < 0.001), with low grip strength in those over 60 years of age significantly positively correlated with frailty (OR = 1.168, 95% CI 1.059-1.289, P = 0.002). The MR results of osteoarthritis and frailty revealed a causal relationship between specific joint sites and frailty, including KOA (OR = 1.086, 95% CI 1.017-1.160, P = 0.014), HOA (OR = 1.028, 95% CI 1.007-1.049, P = 0.009), and KOA/HOA (OR = 1.082, 95% CI 1.053-1.113, P = 0.000), increasing the risk of frailty.

CONCLUSION

Osteosarcopenia, osteoarthritis and frailty exhibit significant causal effects, rendering them risk factors for frailty. Therefore, in clinical practice, patients with osteosarcopenia and osteoarthritis should be required to undergo relevant interventions to reduce the risk of frailty.

Gastroesophageal reflux disease and osteoporosis: A bidirectional Mendelian randomization study

In observational studies, associations between osteoporosis (OP) and gastroesophageal reflux disease (GERD) have been found. We conducted a 2-way, 2-sample Mendelian randomization (MR) analysis to determine whether these associations have a causal relationship. Data on GERD at the summary-level were sourced from extensive genome-wide association studies encompassing 129,080 cases and 473,524 control subjects. Bone mineral density (BMD) served as the phenotypic indicator for OP. BMD metrics were compiled from a cohort of 537,750 individuals, encompassing total body BMD (TB-BMD) and stratified TB-BMD across age groups, along with BMD measurements at 4 anatomical locations: lumbar spine, femoral neck, heel, and ultra-distal forearm. Multiple MR approaches, such as the inverse-variance weighted (IVW) method, MR-Egger regression, and the MR-PRESSO test, were employed, among which findings obtained by IVW method were designated as the primary outcomes. For quality assurance, sensitivity analyses were conducted using the MR-Egger intercept, Cochran Q, and leave-one-out test. There were no significant causal links between genetic inclination towards GERD and reduced BMD levels. Nonetheless, the genetic evidence suggests a causal link between higher BMD levels and lower incidence of GERD [TB-BMD: OR = 0.941, 95% confidence intervals (CI) = 0.910-0.972, P < .001; TB-BMD-1: OR = 0.919, 95% CI = 0.885-0.954, P < .001; TB-BMD-3: OR = 0.945, 95% CI = 0.915-0.977, P = .001; TB-BMD-4: OR = 0.926, 95% CI = 0.896-0.957, P < .001]. Sensitivity analyses corroborate our findings. The MR analysis indicates no significant causal link between genetic inclination towards GERD and OP or reduced BMD within the European demographic. In addition, the study suggests that lower BMD or OP, as predicted by genetics, may contribute to the development of GERD.

Genetic analysis of cis-enhancers associated with bone mineral density and periodontitis in the gene SOST

A haplotype block at the sclerostin (SOST) gene correlates with bone mineral density (BMD) and increased periodontitis risk in smokers. Investigating the putative causal variants within this block, our study aimed to elucidate the impact of linked enhancer elements on gene expression and to evaluate their role in transcription factor (TF) binding. Using CRISPR/dCas9 activation (CRISPRa) screening in SaOS-2 cells, we quantified disease-related enhancer activities regulating SOST expression. Additionally, in SaOS-2 cells, we investigated the influence of the candidate TFs CCAAT/enhancer-binding protein beta (CEBPB) on gene expression by antisense (GapmeR) knockdown, followed by RNA sequencing. The periodontitis-linked SNP rs9783823 displayed a significant cis-activating effect (25-fold change in SOST expression), with the C-allele containing a CEBPB binding motif (position weight matrix (PWM) =  0.98, Pcorrected =  7.7 x 10-7). CEBPB knockdown induced genome-wide upregulation but decreased epithelial-mesenchymal transition genes (P =  0.71, AUC =  2.2 x 10-11). This study identifies a robust SOST cis-activating element linked to BMD and periodontitis, carrying CEBPB binding sites, and highlights CEBPB's impact on epithelial-mesenchymal transition.

Regulation of WNT16 in bone may involve upstream enhancers within CPED1

WNT16 stands up as an essential gene for bone homeostasis. Here, we present new evidence of the functional role of a particular region within WNT16. Performing 4 C chromatin conformation analysis in three osteoblast-related cells (the human fetal osteoblast hFOB 1.19 cell line, Saos 2 osteosarcoma cell line and mesenchymal Stem Cells -MSC-), we identify physical interactions between the proximal part of WNT16 intron 2, shown here to be an active promoter in Saos 2 osteosarcoma cells, and several putative regulatory regions within CPED1. Analysis of previously published RNA-seq data from hFOB cells disclosed low expression of a region located downstream of this promoter. Our results suggest a novel regulatory mechanism of WNT16 in bone, mediated by physical interaction with various enhancer regions within CPED1.

Association between bone mineral density and scoliosis: a two-sample mendelian randomization study in european populations

BACKGROUND

Previous studies have shown that bone mineral density (BMD) has a certain impact on scoliosis. However, up to now, there is no clear evidence that there is a causal association between the two. The aim of this study is to investigate whether there is a causal association between BMD at different body positions and scoliosis by two-sample Mendelian randomization (MR).

METHODS

Genetic variants (SNPS) strongly associated with BMD (total body BMD (TB-BMD), lumbar spine BMD (LS-BMD), femoral neck BMD (FN-BMD), heel BMD (HE-BMD), and forearm BMD (FA-BMD)) were extracted from GEFOS and genome-wide association analysis (GWAS) databases SNPs) were used as instrumental variables (IVs). Scoliosis was also selected from the Finnish database as the outcome. Inverse variance weighting (IVW) method was used as the main analysis method, and multiple sensitivity analysis was performed by combining weighted median, MR-Egger, MR Multi-effect residuals and outliers.

RESULTS

IVW results showed that TB-BMD (OR = 0.83, 95%CI: 0.66-1.55 P = 0.13), LS-BMD (OR = 0.72, 95%CI: 0.52-0.99, P = 0.04), FN-BMD (OR = 0.74, 95%CI: 0.50-1.09, P = 0.13), FA-BMD (OR = 0.95,95%CI: 0.70-1.28, P = 0.75), HE-BMD (OR = 0.91, 95%CI: 0.77-1.08, P = 0.29). Sensitivity analyses showed no evidence of pleiotropy or heterogeneity (p > 0.05) (MR-PRESSO and Cochrane). The results were further validated by leave-one-out test and MR-Egger intercept, which confirmed the robustness of the study results.

CONCLUSION

In conclusion, the present study demonstrates that the causal role of genetic prediction of scoliosis increases with decreasing lumbar BMD. There was no evidence that BMD at the remaining sites has a significant causal effect on scoliosis. Our results suggest that the lumbar spine BMD should be routinely measured in the population at high risk of scoliosis. If osteoporosis occurs, appropriate treatment should be given to reduce the incidence of scoliosis.

CLINICAL TRIAL NUMBER

Not applicable.

Is there a causal relationship between resistin levels and bone mineral density, fracture occurrence? A mendelian randomization study

BACKGROUND

In a great many of observational studies, whether there is a relevance of resistin levels on bone mineral density (BMD) and fracture occurrence has been inconsistently reported, and the causality is unclear.

METHODS

We aim to assess the resistin levels on BMD and fracture occurrence within a Mendelian randomization (MR) analysis. Exposure and outcome data were derived from the Integrative Epidemiology Unit (IEU) Open genome wide association studies (GWAS) database. Screening of instrumental variables (IVs) was performed subject to conditions of relevance, exclusivity, and independence. Inverse variance weighting (IVW) was our primary method for MR analysis based on harmonized data. Weighted median and MR-Egger were chosen to evaluate the robustness of the results of IVW. Simultaneously, heterogeneity and horizontal pleiotropy were also assessed and the direction of potential causality was detected by MR Steiger. Multivariable MR (MVMR) analysis was used to identify whether confounding factors affected the reliability of the results.

RESULTS

After Bonferroni correction, the results showed a suggestively positive causality between resistin levels and total body BMD (TB-BMD) in European populations over the age of 60 [β(95%CI): 0.093(0.021, 0.165), P = 0.011]. The weighted median [β(95%CI): 0.111(0.067, 0.213), P = 0.035] and MR-Egger [β(95%CI): 0.162(0.025, 0.2983), P = 0.040] results demonstrate the robustness of the IVW results. No presence of pleiotropy or heterogeneity was detected between them. MR Steiger supports the causal inference result and MVMR suggests its direct effect.

CONCLUSIONS

In European population older than 60 years, genetically predicted higher levels of resistin were associated with higher TB-BMD. A significant causality between resistin levels on BMD at different sites, fracture in certain parts of the body, and BMD in four different age groups between 0-60 years of age was not found in our study.

Causal Relationship of Obstructive Sleep Apnea with Bone Mineral Density and the Role of BMI

Background

Observational studies have yielded conflicting evidence concerning the relationships between obstructive sleep apnea (OSA) and bone mineral density (BMD). As the exact causal inferences remain inconclusive, we conducted a two-sample Mendelian randomization (MR) to identify the causal associations between OSA and BMD.

Methods

Single-nucleotide polymorphisms associated with OSA were extracted from the FinnGen study. Summary statistics for 10 BMD measured at different age or skeletal sites were obtained from the publicly available IEU GWAS database. Inverse-variance weighted (IVW) method was chosen as the primary analysis, combined with several sensitivity analyses to evaluate the robustness of results. The study design included two-sample MR and network MR.

Results

Our primary MR analysis revealed that genetically predicted OSA was positively linked to increased forearm BMD (β = 0.24, 95% confidence interval [CI]: 0.06-0.41, p = 0.009) and heel BMD (β=0.10, 95% CI = 0.02-0.18, p = 0.018), while no significant causal relationships were observed between OSA and total body BMD, lumbar spine BMD, or femoral neck BMD (all p > 0.05). Network MR suggests that OSA might act as a mediating factor in the effect of BMI on forearm BMD and heel BMD, with a mediated portion estimated at 73% and 84%, respectively.

Conclusion

Our findings provide support for a causal relationship between genetically predicted OSA and increased forearm BMD and heel BMD. Furthermore, our results suggest that OSA may play a role in mediating the influence of BMI on BMD.

Association of coffee intake with bone mineral density: a Mendelian randomization study

Background

In observational studies, the relationship between coffee intake and bone mineral density (BMD) is contradictory. However, residual confounding tends to bias the results of these studies. Therefore, we used a two-sample Mendelian randomization (MR) approach to further investigate the potential causal relationship between the two.

Methods

Genetic instrumental variables (IVs) associated with coffee intake were derived from genome-wide association studies (GWAS) of the Food Frequency Questionnaire (FFQ) in 428,860 British individuals and matched using phenotypes in PhenoScanner. Summarized data on BMD were obtained from 537,750 participants, including total body BMD (TB-BMD), TB-BMD in five age brackets ≥60, 45-60, 30-45, 15-30, and 0-15 years, and BMD in four body sites: the lumbar spine, the femoral neck, the heel, and the ultradistal forearm. We used inverse variance weighting (IVW) methods as the primary analytical method for causal inference. In addition, several sensitivity analyses (MR-Egger, Weighted median, MR-PRESSO, Cochran's Q test, and Leave-one-out test) were used to test the robustness of the results.

Results

After Bonferroni correction, Coffee intake has a potential positive correlation with total body BMD (effect estimate [Beta]: 0.198, 95% confidence interval [Cl]: 0.05-0.35, P=0.008). In subgroup analyses, coffee intake was potentially positively associated with TB-BMD (45-60, 30-45 years) (Beta: 0.408, 95% Cl: 0.12-0.69, P=0.005; Beta: 0.486, 95% Cl: 0.12-0.85, P=0.010). In addition, a significant positive correlation with heel BMD was also observed (Beta: 0.173, 95% Cl: 0.08-0.27, P=0.002). The results of the sensitivity analysis were generally consistent.

Conclusion

The results of the present study provide genetic evidence for the idea that coffee intake is beneficial for bone density. Further studies are needed to reveal the biological mechanisms and offer solid support for clinical guidelines on osteoporosis prevention.

The osteoporosis susceptibility SNP rs188303909 at 2q14.2 regulates EN1 expression by modulating DNA methylation and E2F6 binding

EN1 encodes a homeodomain-containing transcription factor and is a determinant of bone density and fracture. Previous powerful genome-wide association studies (GWASs) have identified multiple single-nucleotide polymorphisms (SNPs) near EN1 at 2q14.2 locus for osteoporosis, but the causal SNPs and functional mechanisms underlying these associations are poorly understood. The target genes regulated by the transcription factor EN1 are also unclear. In this study, we identified rs188303909, a functional CpG-SNP, as a causal SNP for osteoporosis at 2q14.2 through the integration of functional and epigenomic analyses. Functional experiments demonstrated that unmethylated rs188303909 acted as a strong allele-specific distal enhancer to regulate EN1 expression by modifying the binding of transcription factor E2F6, but rs188303909 methylation attenuated the active effect of E2F6 on EN1 expression. Importantly, transcription factor EN1 could differentially bind osteoporosis GWAS lead SNPs rs4869739-T and rs4355801-G to upregulate CCDC170 and COLEC10 expression, thus promoting bone formation. Our study provided a mechanistic insight into expression regulation of the osteoporosis susceptibility gene EN1, which could be a potential therapeutic target for osteoporosis precision medicine. KEY MESSAGES: CpG-SNP rs188303909 is a causal SNP at the osteoporosis susceptibility locus 2q14.2. Rs188303909 distally regulates EN1 expression by modulating DNA methylation and E2F6 binding. EN1 upregulates CCDC170 and COLEC10 expression through osteoporosis GWAS lead SNPs rs4869739 and rs4355801.

Converging evidence from exome sequencing and common variants implicates target genes for osteoporosis

In this study, we leveraged the combined evidence of rare coding variants and common alleles to identify therapeutic targets for osteoporosis. We undertook a large-scale multiancestry exome-wide association study for estimated bone mineral density, which showed that the burden of rare coding alleles in 19 genes was associated with estimated bone mineral density (P < 3.6 × 10-7). These genes were highly enriched for a set of known causal genes for osteoporosis (65-fold; P = 2.5 × 10-5). Exome-wide significant genes had 96-fold increased odds of being the top ranked effector gene at a given GWAS locus (P = 1.8 × 10-10). By integrating proteomics Mendelian randomization evidence, we prioritized CD109 (cluster of differentiation 109) as a gene for which heterozygous loss of function is associated with higher bone density. CRISPR-Cas9 editing of CD109 in SaOS-2 osteoblast-like cell lines showed that partial CD109 knockdown led to increased mineralization. This study demonstrates that the convergence of common and rare variants, proteomics and CRISPR can highlight new bone biology to guide therapeutic development.

Metabolome-Wide Mendelian Randomization Assessing the Causal Relationship Between Blood Metabolites and Bone Mineral Density

Mounting evidence has supported osteoporosis (OP) as a metabolic disorder. Recent metabolomics studies have discovered numerous metabolites related to bone mineral density (BMD). However, the causal effects of metabolites on BMD at distinct sites remained underexplored. Leveraging genome-wide association datasets, we conducted two-sample Mendelian randomization (MR) analyses to investigate the causal relationship between 486 blood metabolites and bone mineral density at five skeletal sites including heel (H), total body (TB), lumbar spine (LS), femoral neck (FN), and ultra-distal forearm (FA). Sensitivity analyses were performed to test the presence of the heterogeneity and the pleiotropy. To exclude the influences of reverse causation, genetic correlation, and linkage disequilibrium (LD), we further performed reverse MR, linkage disequilibrium regression score (LDSC), and colocalization analyses. In the primary MR analyses, 22, 10, 3, 7, and 2 metabolite associations were established respectively for H-BMD, TB-BMD, LS-BMD, FN-BMD, and FA-BMD at the nominal significance level (IVW, P < 0.05) and passing sensitivity analyses. Among these, one metabolite, androsterone sulfate showed a strong effect on four out of five BMD phenotypes (Odds ratio [OR] for H-BMD = 1.045 [1.020, 1.071]; Odds ratio [OR] for TB-BMD = 1.061 [1.017, 1.107]; Odds ratio [OR] for LS-BMD = 1.088 [1.023, 1.159]; Odds ratio [OR] for FN-BMD = 1.114 [1.054, 1.177]). Reverse MR analysis provided no evidence for the causal effects of BMD measurements on these metabolites. Colocalization analysis have found that several metabolite associations might be driven by shared genetic variants such as mannose for TB-BMD. This study identified some metabolites causally related to BMD at distinct sites and several key metabolic pathways, which shed light on predictive biomarkers and drug targets for OP.

The Michigan Genomics Initiative: A biobank linking genotypes and electronic clinical records in Michigan Medicine patients

Biobanks of linked clinical patient histories and biological samples are an efficient strategy to generate large cohorts for modern genetics research. Biobank recruitment varies by factors such as geographic catchment and sampling strategy, which affect biobank demographics and research utility. Here, we describe the Michigan Genomics Initiative (MGI), a single-health-system biobank currently consisting of >91,000 participants recruited primarily during surgical encounters at Michigan Medicine. The surgical enrollment results in a biobank enriched for many diseases and ideally suited for a disease genetics cohort. Compared with the much larger population-based UK Biobank, MGI has higher prevalence for nearly all diagnosis-code-based phenotypes and larger absolute case counts for many phenotypes. Genome-wide association study (GWAS) results replicate known findings, thereby validating the genetic and clinical data. Our results illustrate that opportunistic biobank sampling within single health systems provides a unique and complementary resource for exploring the genetics of complex diseases.

Genome-wide screen of otosclerosis in population biobanks: 27 loci and shared associations with skeletal structure

Otosclerosis is one of the most common causes of conductive hearing loss, affecting 0.3% of the population. It typically presents in adulthood and half of the patients have a positive family history. The pathophysiology of otosclerosis is poorly understood. A previous genome-wide association study (GWAS) identified a single association locus in an intronic region of RELN. Here, we report a meta-analysis of GWAS studies of otosclerosis in three population-based biobanks comprising 3504 cases and 861,198 controls. We identify 23 novel risk loci (p < 5 × 10-8) and report an association in RELN and three previously reported candidate gene or linkage regions (TGFB1, MEPE, and OTSC7). We demonstrate developmental stage-dependent immunostaining patterns of MEPE and RUNX2 in mouse otic capsules. In most association loci, the nearest protein-coding genes are implicated in bone remodelling, mineralization or severe skeletal disorders. We highlight multiple genes involved in transforming growth factor beta signalling for follow-up studies.

Kidney tonifying traditional Chinese medicine: Potential implications for the prevention and treatment of osteoporosis

The aging global population is increasingly affected by osteoporosis (OP), which is one of the most significant threats to the elderly. Moreover, its prevention and treatment situations have become increasingly severe. Therefore, it is imperative to develop alternatives or complementary drugs for preventing and treating osteoporosis. Kidney tonifying traditional Chinese medicine (KTTCM) has been used for the treatment of osteoporosis for a long time. Pharmacological studies have shown that kidney tonifying traditional Chinese medicine can promote osteoblasts, inhibit osteoclasts, and regulate the level of estrogen and plays vital roles in stimulating osteogenesis, restraining adipogenesis of marrow mesenchymal stem cells (MSCs), regulating the metabolism of calcium and phosphorus, and inhibiting oxidative stress. These effects are mediated by OPG/RANKL/RANK, BMP/Smads, MAPKs, and Wnt/β-catenin systems. To develop a safe, synergistic, effective, and homogenized TCM formula with robust scientific evidence to provide faster and more economical alternatives, the anti-osteoporosis ingredients and pharmacological mechanisms of kidney tonifying traditional Chinese medicine are recapitulated from the perspective of molecular and cell biology, and the safety and toxicity of kidney tonifying traditional Chinese medicine have also been reviewed in this paper.

Polygenic Liability to Depression Is Associated With Multiple Medical Conditions in the Electronic Health Record: Phenome-wide Association Study of 46,782 Individuals

BACKGROUND

Major depressive disorder (MDD) is a leading cause of disease-associated disability, with much of the increased burden due to psychiatric and medical comorbidity. This comorbidity partly reflects common genetic influences across conditions. Integrating molecular-genetic tools with health records enables tests of association with the broad range of physiological and clinical phenotypes. However, standard phenome-wide association studies analyze associations with individual genetic variants. For polygenic traits such as MDD, aggregate measures of genetic risk may yield greater insight into associations across the clinical phenome.

METHODS

We tested for associations between a genome-wide polygenic risk score for MDD and medical and psychiatric traits in a phenome-wide association study of 46,782 unrelated, European-ancestry participants from the Michigan Genomics Initiative.

RESULTS

The MDD polygenic risk score was associated with 211 traits from 15 medical and psychiatric disease categories at the phenome-wide significance threshold. After excluding patients with depression, continued associations were observed with respiratory, digestive, neurological, and genitourinary conditions; neoplasms; and mental disorders. Associations with tobacco use disorder, respiratory conditions, and genitourinary conditions persisted after accounting for genetic overlap between depression and other psychiatric traits. Temporal analyses of time-at-first-diagnosis indicated that depression disproportionately preceded chronic pain and substance-related disorders, while asthma disproportionately preceded depression.

CONCLUSIONS

The present results can inform the biological links between depression and both mental and systemic diseases. Although MDD polygenic risk scores cannot currently forecast health outcomes with precision at the individual level, as molecular-genetic discoveries for depression increase, these tools may augment risk prediction for medical and psychiatric conditions.

Biobanking as a Tool for Genomic Research: From Allele Frequencies to Cross-Ancestry Association Studies

In recent years, great advances have been made in the field of collection, storage, and analysis of biological samples. Large collections of samples, biobanks, have been established in many countries. Biobanks typically collect large amounts of biological samples and associated clinical information; the largest collections include over a million samples. In this review, we summarize the main directions in which biobanks aid medical genetics and genomic research, from providing reference allele frequency information to allowing large-scale cross-ancestry meta-analyses. The largest biobanks greatly vary in the size of the collection, and the amount of available phenotype and genotype data. Nevertheless, all of them are extensively used in genomics, providing a rich resource for genome-wide association analysis, genetic epidemiology, and statistical research into the structure, function, and evolution of the human genome. Recently, multiple research efforts were based on trans-biobank data integration, which increases sample size and allows for the identification of robust genetic associations. We provide prominent examples of such data integration and discuss important caveats which have to be taken into account in trans-biobank research.

The HUNT study: A population-based cohort for genetic research

The Trøndelag Health Study (HUNT) is a population-based cohort of ∼229,000 individuals recruited in four waves beginning in 1984 in Trøndelag County, Norway. Approximately 88,000 of these individuals have available genetic data from array genotyping. HUNT participants were recruited during four community-based recruitment waves and provided information on health-related behaviors, self-reported diagnoses, family history of disease, and underwent physical examinations. Linkage via the Norwegian personal identification number integrates digitized health care information from doctor visits and national health registries including death, cancer and prescription registries. Genome-wide association studies of HUNT participants have provided insights into the mechanism of cardiovascular, metabolic, osteoporotic, and liver-related diseases, among others. Unique features of this cohort that facilitate research include nearly 40 years of longitudinal follow-up in a motivated and well-educated population, family data, comprehensive phenotyping, and broad availability of DNA, RNA, urine, fecal, plasma, and serum samples.

Osteoporosis epidemiology using international cohorts

PURPOSE OF REVIEW

To provide an update on the most important new cohort studies within osteoporosis and their bearing on clinical management and directions for future research.

RECENT FINDINGS

We identified a collection of new observational cohort studies - including new reports from already established large cohorts - and intervention studies providing new insights into osteoporosis pathophysiology, risk finding, intervention, and treatment barriers.

SUMMARY

Recent cohort studies in osteoporosis highlight the importance of timely identification and treatment of people who are at high risk of suffering osteoporotic fractures. Physical performance is a strong indicator of fracture risk and one that is tightly linked to a number of chronic conditions, not least inflammatory conditions like rheumatoid arthritis. Advances in case finding may involve opportunistic screening for low bone mineral density and vertebral fractures of radiology images obtained for other purposes, polygenic risk scores, and routinely collected medication and comorbidity information.

CNest: A novel copy number association discovery method uncovers 862 new associations from 200,629 whole-exome sequence datasets in the UK Biobank

Copy number variation (CNV) is known to influence human traits, having a rich history of research into common and rare genetic disease, and although CNV is accepted as an important class of genomic variation, progress on copy-number-based genome-wide association studies (GWASs) from next-generation sequencing (NGS) data has been limited. Here we present a novel method for large-scale copy number analysis from NGS data generating robust copy number estimates and allowing copy number GWASs (CN-GWASs) to be performed genome-wide in discovery mode. We provide a detailed analysis in the UK Biobank resource and a specifically designed software package. We use these methods to perform CN-GWAS analysis across 78 human traits, discovering over 800 genetic associations that are likely to contribute strongly to trait distributions. Finally, we compare CNV and SNP association signals across the same traits and samples, defining specific CNV association classes.

Genetic variation in WNT16 and its association with bone mineral density, fractures and osteoporosis in children with bone fragility

Several genome-wide association studies (GWAS), GWAS meta-analyses, and mouse studies have demonstrated that wingless-related integration site 16 (WNT16) gene is associated with bone mineral density (BMD), cortical bone thickness, bone strength and fracture risk. Practically no data exist regarding the significance of WNT16 in childhood-onset osteoporosis and related fractures. We hypothesized that pathogenic variants and genetic variations in WNT16 could explain skeletal fragility in affected children. We screened the WNT16 gene by Sanger sequencing in three pediatric cohorts: 35 with primary osteoporosis, 59 with multiple fractures, and in 95 healthy controls. Altogether, we identified 12 variants in WNT16. Of them one was a rare 5′UTR variant rs1386898215 in genome aggregate and medical trans-omic databases (GnomAD, TOPMED; minor allele frequency (MAF) 0.00 and 0.000008, respectively). One variant rs1554366753, overrepresented in children with osteoporosis (MAF = 0.06 vs healthy controls MAF = 0.01), was significantly associated with lower BMD. This variant was found associated with increased WNT16 gene expression at mRNA level in fibroblast cultures. None of the other identified variants were rare (MAF < 0.001) or deemed pathogenic by predictor programs. WNT16 may play a role in childhood osteoporosis but genetic WNT16 variation is not a common cause of skeletal fragility in childhood.

•

No pathogenic WNT16 variants were found associated with pediatric osteoporosis or fracture-prone patients

•

Altogether, twelve WNT16 variants were found in pediatric osteoporosis or fracture-prone patients

•

The genetic variation rs1554366753 in the WNT16 gene is associated with bone mineral density and primary osteoporosis

A timeseries analysis of the fracture callus extracellular matrix proteome during bone fracture healing

While most bones fully self-heal, certain diseases require bone allograft to assist with fracture healing. Bone allografts offer promise as treatments for such fractures due to their osteogenic properties. However, current bone allografts made of decellularized bone extracellular matrix (ECM) have high failure rates, and thus grafts which improve fracture healing outcomes are needed. Understanding specific changes to the ECM proteome during normal fracture healing would enable the identification of key proteins that could be used enhance osteogenicity of bone allograft. Here, we performed a timeseries analysis of the fracture callus in mice to investigate proteomic and mineralization changes to the ECM at key stages of fracture healing. We found that changes to the ECM proteome largely coincide with the distinct phases of fracture healing. Basement membrane proteins (AGRN, COL4, LAMA), cartilage proteins (COL2A1, ACAN), and collagen crosslinking enzymes (LOXL, PLOD, ITIH) were initially upregulated, followed by bone specific proteoglycans and collagens (IBSP, COL1A1). Various tissue proteases (MMP2, 9, 13, 14; CTSK, CTSG, ELANE) were expressed at different levels throughout fracture healing. These changes coordinated with mineralization of the fracture callus, which increased steeply during the initial stages of healing. Interestingly the later timepoint was characterized by a response to wound healing and high expression of clotting factors (F2, 7, 9, 10). We identified ELANE and ITIH2 as tissue remodeling enzymes having no prior known involvement with fracture healing. This data can be further mined to identify regenerative proteins for enhanced bone graft design.

snpXplorer: a web application to explore human SNP-associations and annotate SNP-sets

Genetic association studies are frequently used to study the genetic basis of numerous human phenotypes. However, the rapid interrogation of how well a certain genomic region associates across traits as well as the interpretation of genetic associations is often complex and requires the integration of multiple sources of annotation, which involves advanced bioinformatic skills. We developed snpXplorer, an easy-to-use web-server application for exploring Single Nucleotide Polymorphisms (SNP) association statistics and to functionally annotate sets of SNPs. snpXplorer can superimpose association statistics from multiple studies, and displays regional information including SNP associations, structural variations, recombination rates, eQTL, linkage disequilibrium patterns, genes and gene-expressions per tissue. By overlaying multiple GWAS studies, snpXplorer can be used to compare levels of association across different traits, which may help the interpretation of variant consequences. Given a list of SNPs, snpXplorer can also be used to perform variant-to-gene mapping and gene-set enrichment analysis to identify molecular pathways that are overrepresented in the list of input SNPs. snpXplorer is freely available at https://snpxplorer.net. Source code, documentation, example files and tutorial videos are available within the Help section of snpXplorer and at https://github.com/TesiNicco/snpXplorer.