GWAS-informed data integration and non-coding CRISPRi screen illuminate genetic etiology of bone mineral density.

Over 1,100 independent signals have been identified with genome-wide association studies (GWAS) for bone mineral density (BMD), a key risk factor for mortality-increasing fragility fractures; however, the effector gene(s) for most remain unknown. Informed by a variant-to-gene mapping strategy implicating 89 non-coding elements predicted to regulate osteoblast gene expression at BMD GWAS loci, we executed a single-cell CRISPRi screen in human fetal osteoblast 1.19 cells (hFOBs). The BMD relevance of hFOBs was supported by heritability enrichment from cross-cell type stratified LD-score regression involving 98 cell types grouped into 15 tissues. 24 genes showed perturbation in the screen, with four (ARID5B, CC2D1B, EIF4G2, and NCOA3) exhibiting consistent effects upon siRNA knockdown on three measures of osteoblast maturation and mineralization. Lastly, additional heritability enrichments, genetic correlations, and multi-trait fine-mapping revealed that many BMD GWAS signals are pleiotropic and likely mediate their effects via non-bone tissues that warrant attention in future screens.

A perspective on muscle phenotyping in musculoskeletal research.

Musculoskeletal research should synergistically investigate bone and muscle to inform approaches for maintaining mobility and to avoid bone fractures. The relationship between sarcopenia and osteoporosis, integrated in the term 'osteosarcopenia', is underscored by the close association shown between these two conditions in many studies, whereby one entity emerges as a predictor of the other. In a recent workshop of Working Group (WG) 2 of the EU Cooperation in Science and Technology (COST) Action 'Genomics of MusculoSkeletal traits Translational Network' (GEMSTONE) consortium (CA18139), muscle characterization was highlighted as being important, but currently under-recognized in the musculoskeletal field. Here, we summarize the opinions of the Consortium and research questions around translational and clinical musculoskeletal research, discussing muscle phenotyping in human experimental research and in two animal models: zebrafish and mouse.

Sex-specific effects of Fat-1 transgene on bone material properties, size, and shape in mice.

Western diets are becoming increasingly common around the world. Western diets have high omega 6 (ω-6) and omega 3 (ω-3) fatty acids and are linked to bone loss in humans and animals. Dietary fats are not created equal; therefore, it is vital to understand the effects of specific dietary fats on bone. We aimed to determine how altering the endogenous ratios of ω-6:ω-3 fatty acids impacts bone accrual, strength, and fracture toughness. To accomplish this, we used the Fat-1 transgenic mice, which carry a gene responsible for encoding a ω-3 fatty acid desaturase that converts ω-6 to ω-3 fatty acids. Male and female Fat-1 positive mice (Fat-1) and Fat-1 negative littermates (WT) were given either a high-fat diet (HFD) or low-fat diet (LFD) at 4 wk of age for 16 wk. The Fat-1 transgene reduced fracture toughness in males. Additionally, male BMD, measured from DXA, decreased over the diet duration for HFD mice. In males, neither HFD feeding nor the presence of the Fat-1 transgene impacted cortical geometry, trabecular architecture, or whole-bone flexural properties, as detected by main group effects. In females, Fat-1-LFD mice experienced increases in BMD compared to WT-LFD mice; however, cortical area, distal femur trabecular thickness, and cortical stiffness were reduced in Fat-1 mice compared to pooled WT controls. However, reductions in stiffness were caused by a decrease in bone size and were not driven by changes in material properties. Together, these results demonstrate that the endogenous ω-6:ω-3 fatty acid ratio influences bone material properties in a sex-dependent manner. In addition, Fat-1 mediated fatty acid conversion was not able to mitigate the adverse effects of HFD on bone strength and accrual.

Zebrafish as a Model for Osteoporosis: Functional Validations of Genome-Wide Association Studies.

PURPOSE OF REVIEW: GWAS, as a largely correlational analysis, requires in vitro or in vivo validation. Zebrafish (Danio rerio) have many advantages for studying the genetics of human diseases. Since gene editing in zebrafish has been highly valuable for studying embryonic skeletal developmental processes that are prenatally or perinatally lethal in mammalian models, we are reviewing pros and cons of this model. RECENT FINDINGS: The true power for the use of zebrafish is the ease by which the genome can be edited, especially using the CRISPR/Cas9 system. Gene editing, followed by phenotyping, for complex traits such as BMD, is beneficial, but the major physiological differences between the fish and mammals must be considered. Like mammals, zebrafish do have main bone cells; thus, both in vivo stem cell analyses and in vivo imaging are doable. Yet, the "long" bones of fish are peculiar, and their bone cavities do not contain bone marrow. Partial duplication of the zebrafish genome should be taken into account. Overall, small fish toolkit can provide unmatched opportunities for genetic modifications and morphological investigation as a follow-up to human-first discovery.

SEESAW: detecting isoform-level allelic imbalance accounting for inferential uncertainty.

Detecting allelic imbalance at the isoform level requires accounting for inferential uncertainty, caused by multi-mapping of RNA-seq reads. Our proposed method, SEESAW, uses Salmon and Swish to offer analysis at various levels of resolution, including gene, isoform, and aggregating isoforms to groups by transcription start site. The aggregation strategies strengthen the signal for transcripts with high uncertainty. The SEESAW suite of methods is shown to have higher power than other allelic imbalance methods when there is isoform-level allelic imbalance. We also introduce a new test for detecting imbalance that varies across a covariate, such as time.

Functional screens refine height GWAS loci.

Genome-wide association studies (GWASs) have demonstrated the complexity of human height. Baronas et al.1 used a high-throughput CRISPR screen to identify genes that participate in growth plate chondrocyte maturation as a functional follow-up and validation screen to refine loci and establish causality after GWASs.

Proceedings of the Post-Genome Analysis for Musculoskeletal Biology Workshop.

PURPOSE OF THE REVIEW: Herein, we report on the proceedings of the workshop entitled "Post-Genome analysis for musculoskeletal biology" that was held in July of 2022 in Safed, Galilee, Israel. Supported by the Israel Science Foundation, the goal of this workshop was to bring together established investigators and their trainees who were interested in understanding the etiology of musculoskeletal disease, from Israel and from around the world. RECENT FINDINGS: Presentations at this workshop spanned the spectrum from basic science to clinical studies. A major emphasis of the discussion centered on genetic studies in humans, and the limitations and advantages of such studies. The power of coupling studies using human data with functional follow-up studies in pre-clinical models such as mice, rats, and zebrafish was discussed in depth. The advantages and limitations of mice and zebrafish for faithfully modelling aspects of human disease were debated, specifically in the context of age-related diseases such as osteoporosis, osteoarthritis, adult-onset auto-immune disease, and osteosarcopenia. There remain significant gaps in our understanding of the nature and etiology of human musculoskeletal disease. While therapies and medications exist, much work is still needed to find safe and effective interventions for all patients suffering from diseases associated with age-related deterioration of musculoskeletal tissues. The potential of forward and reverse genetic studies has not been exhausted for diseases of muscles, joints, and bones.

Genetics of Intervertebral Disc Degeneration.

PURPOSE OF REVIEW: Intervertebral disc degeneration is a contributor to chronic back pain. While a part of the natural aging process, early or rapid intervertebral disc degeneration is highly heritable. In this review, we summarize recent progress towards unraveling the genetics associated with this degenerative process. RECENT FINDINGS: Use of large cohorts of patient data to conduct genome-wide association studies (GWAS) for intervertebral disc disease, and to lesser extent for aspects of this process, such as disc height, has resulted in a large increase in our understanding of the genetic etiology. Genetic correlation suggests that intervertebral disc disease is pleiotropic with risk factors for other diseases such as osteoporosis. The use of Mendelian Randomization is slowly establishing what are the causal relationships between intervertebral disc disease and factors previously correlated with this disease. The results from these human genetic studies highlight the complex nature of this disease and have the potential to lead to improved clinical management of intervertebral disc disease. Much additional work should now be focused on characterizing the causative relationship various co-morbid conditions have with intervertebral disc degeneration and on finding interventions to slow or halt this disease.

Whole-genome RNA sequencing identifies distinct transcriptomic profiles in impingement cartilage between patients with femoroacetabular impingement and hip osteoarthritis.

Femoroacetabular impingement (FAI) has a strong clinical association with the development of hip osteoarthritis (OA); however, the pathobiological mechanisms underlying the transition from focal impingement to global joint degeneration remain poorly understood. The purpose of this study is to use whole-genome RNA sequencing to identify and subsequently validate differentially expressed genes (DEGs) in femoral head articular cartilage samples from patients with FAI and hip OA secondary to FAI. Thirty-seven patients were included in the study with whole-genome RNA sequencing performed on 10 gender-matched patients in the FAI and OA cohorts and the remaining specimens were used for validation analyses. We identified a total of 3531 DEGs between the FAI and OA cohorts with multiple targets for genes implicated in canonical OA pathways. Quantitative reverse transcription-polymerase chain reaction validation confirmed increased expression of FGF18 and WNT16 in the FAI samples, while there was increased expression of MMP13 and ADAMTS4 in the OA samples. Expression levels of FGF18 and WNT16 were also higher in FAI samples with mild cartilage damage compared to FAI samples with severe cartilage damage or OA cartilage. Our study further expands the knowledge regarding distinct genetic reprogramming in the cartilage between FAI and hip OA patients. We independently validated the results of the sequencing analysis and found increased expression of anabolic markers in patients with FAI and minimal histologic cartilage damage, suggesting that anabolic signaling may be increased in early FAI with a transition to catabolic and inflammatory gene expression as FAI progresses towards more severe hip OA. Clinical significance:Cam-type FAI has a strong clinical association with hip OA; however, the cellular pathophysiology of disease progression remains poorly understood. Several previous studies have demonstrated increased expression of inflammatory markers in FAI cartilage samples, suggesting the involvement of these inflammatory pathways in the disease progression. Our study further expands the knowledge regarding distinct genetic reprogramming in the cartilage between FAI and hip OA patients. In addition to differences in inflammatory gene expression, we also identified differential expression in multiple pathways involved in hip OA progression.

Osteoporosis: An Update on Screening, Diagnosis, Evaluation, and Treatment.

Osteoporosis screening, diagnosis, and treatment have gained much attention in the health care community over the past 2 decades. During this time, creation of multispecialty awareness programs (eg, "Own the Bone," American Orthopedic Association; "Capture the Fracture," International Osteoporosis Foundation) and improvements in diagnostic protocols have been evident. Significant advances in technology have elucidated elements of genetic predisposition for decreased bone mineral density in the aging population. Additionally, several novel drug therapies have entered the market and provide more options for primary care and osteoporosis specialists to medically manage patients at risk for fragility fractures. Despite this, adherence to osteoporosis screening and treatment protocols has been surprisingly low by health care practitioners, including orthopedic surgeons. Continued awareness and education of this skeletal disorder is crucial to effectively care for our aging population. [Orthopedics. 2023;46(1):e20-e26.].

Transcriptome-wide association study and eQTL colocalization identify potentially causal genes responsible for human bone mineral density GWAS associations.

Genome-wide association studies (GWASs) for bone mineral density (BMD) in humans have identified over 1100 associations to date. However, identifying causal genes implicated by such studies has been challenging. Recent advances in the development of transcriptome reference datasets and computational approaches such as transcriptome-wide association studies (TWASs) and expression quantitative trait loci (eQTL) colocalization have proven to be informative in identifying putatively causal genes underlying GWAS associations. Here, we used TWAS/eQTL colocalization in conjunction with transcriptomic data from the Genotype-Tissue Expression (GTEx) project to identify potentially causal genes for the largest BMD GWAS performed to date. Using this approach, we identified 512 genes as significant using both TWAS and eQTL colocalization. This set of genes was enriched for regulators of BMD and members of bone relevant biological processes. To investigate the significance of our findings, we selected PPP6R3, the gene with the strongest support from our analysis which was not previously implicated in the regulation of BMD, for further investigation. We observed that Ppp6r3 deletion in mice decreased BMD. In this work, we provide an updated resource of putatively causal BMD genes and demonstrate that PPP6R3 is a putatively causal BMD GWAS gene. These data increase our understanding of the genetics of BMD and provide further evidence for the utility of combined TWAS/colocalization approaches in untangling the genetics of complex traits.

Understanding the Transcriptomic Landscape to Drive New Innovations in Musculoskeletal Regenerative Medicine.

PURPOSE OF REVIEW: RNA-sequencing (RNA-seq) is a novel and highly sought-after tool in the field of musculoskeletal regenerative medicine. The technology is being used to better understand pathological processes, as well as elucidate mechanisms governing development and regeneration. It has allowed in-depth characterization of stem cell populations and discovery of molecular mechanisms that regulate stem cell development, maintenance, and differentiation in a way that was not possible with previous technology. This review introduces RNA-seq technology and how it has paved the way for advances in musculoskeletal regenerative medicine. RECENT FINDINGS: Recent studies in regenerative medicine have utilized RNA-seq to decipher mechanisms of pathophysiology and identify novel targets for regenerative medicine. The technology has also advanced stem cell biology through in-depth characterization of stem cells, identifying differentiation trajectories and optimizing cell culture conditions. It has also provided new knowledge that has led to improved growth factor use and scaffold design for musculoskeletal regenerative medicine. This article reviews recent studies utilizing RNA-seq in the field of musculoskeletal regenerative medicine. It demonstrates how transcriptomic analysis can be used to provide insights that can aid in formulating a regenerative strategy.

HMGB1-mediated restriction of EPO signaling contributes to anemia of inflammation.

Anemia of inflammation, also known as anemia of chronic disease, is refractory to erythropoietin (EPO) treatment, but the mechanisms underlying the EPO refractory state are unclear. Here, we demonstrate that high mobility group box-1 protein (HMGB1), a damage-associated molecular pattern molecule recently implicated in anemia development during sepsis, leads to reduced expansion and increased death of EPO-sensitive erythroid precursors in human models of erythropoiesis. HMGB1 significantly attenuates EPO-mediated phosphorylation of the Janus kinase 2/STAT5 and mTOR signaling pathways. Genetic ablation of receptor for advanced glycation end products, the only known HMGB1 receptor expressed by erythroid precursors, does not rescue the deleterious effects of HMGB1 on EPO signaling, either in human or murine precursors. Furthermore, surface plasmon resonance studies highlight the ability of HMGB1 to interfere with the binding between EPO and the EPOR. Administration of a monoclonal anti-HMGB1 antibody after sepsis onset in mice partially restores EPO signaling in vivo. Thus, HMGB1-mediated restriction of EPO signaling contributes to the chronic phase of anemia of inflammation.

A computational approach for identification of core modules from a co-expression network and GWAS data.

This protocol describes the application of the "omnigenic" model of the genetic architecture of complex traits to identify novel "core" genes influencing a disease-associated phenotype. Core genes are hypothesized to directly regulate disease and may serve as therapeutic targets. This protocol leverages GWAS data, a co-expression network, and publicly available data, including the GTEx database and the International Mouse Phenotyping Consortium Database, to identify modules enriched for genes with "core-like" characteristics. For complete details on the use and execution of this protocol, please refer to Sabik et al. (2020).

Genome-wide meta-analysis of muscle weakness identifies 15 susceptibility loci in older men and women.

Low muscle strength is an important heritable indicator of poor health linked to morbidity and mortality in older people. In a genome-wide association study meta-analysis of 256,523 Europeans aged 60 years and over from 22 cohorts we identify 15 loci associated with muscle weakness (European Working Group on Sarcopenia in Older People definition: n = 48,596 cases, 18.9% of total), including 12 loci not implicated in previous analyses of continuous measures of grip strength. Loci include genes reportedly involved in autoimmune disease (HLA-DQA1 p = 4 × 10-17), arthritis (GDF5 p = 4 × 10-13), cell cycle control and cancer protection, regulation of transcription, and others involved in the development and maintenance of the musculoskeletal system. Using Mendelian randomization we report possible overlapping causal pathways, including diabetes susceptibility, haematological parameters, and the immune system. We conclude that muscle weakness in older adults has distinct mechanisms from continuous strength, including several pathways considered to be hallmarks of ageing.

Isolation and Culture of Neonatal Mouse Calvarial Osteoblasts.

This chapter describes the isolation and culture of neonatal mouse calvarial osteoblasts. This primary cell population is obtained by sequential enzymatic digestion of the calvarial bone matrix and is capable of differentiating in vitro into mature osteoblasts that deposit a collagen extracellular matrix and form mineralized bone nodules. Maturation of the cultures can be monitored by gene expression analyses and staining for the presence of alkaline phosphatase or matrix mineralization. This culture system, therefore, provides a powerful model in which to test how various experimental conditions, such as the manipulation of gene expression, may affect osteoblast maturation and/or function.

Identification of a Core Module for Bone Mineral Density through the Integration of a Co-expression Network and GWAS Data.

The "omnigenic" model of the genetic architecture of complex traits proposed two categories of causal genes: core and peripheral. Core genes are hypothesized to directly regulate disease and may serve as therapeutic targets. Using a cell-type- and time-point-specific gene co-expression network for mineralizing osteoblasts, we identify a co-expression module enriched for genes implicated by bone mineral density (BMD) genome-wide association studies (GWASs), correlated with in vitro osteoblast mineralization and associated with skeletal phenotypes in human monogenic disease and mouse knockouts. Four genes from this module (B4GALNT3, CADM1, DOCK9, and GPR133) are located within the BMD GWAS loci with colocalizing expression quantitative trait loci (eQTL) and exhibit altered BMD in mouse knockouts, suggesting that they are causal genetic drivers of BMD in humans. Our network-based approach identifies a "core" module for BMD and provides a resource for expanding our understanding of the genetics of bone mass.

A Bioinformatic Approach to Utilize a Patient's Antibody-Secreting Cells against Staphylococcus aureus to Detect Challenging Musculoskeletal Infections.

Noninvasive diagnostics for Staphylococcus aureus musculoskeletal infections (MSKI) remain challenging. Abs from newly activated, pathogen-specific plasmablasts in human blood, which emerge during an ongoing infection, can be used for diagnosing and tracking treatment response in diabetic foot infections. Using multianalyte immunoassays on medium enriched for newly synthesized Abs (MENSA) from Ab-secreting cells, we assessed anti-S. aureus IgG responses in 101 MSKI patients (63 culture-confirmed S. aureus, 38 S. aureus-negative) and 52 healthy controls. MENSA IgG levels were assessed for their ability to identify the presence and type of S. aureus MSKI using machine learning and multivariate receiver operating characteristic curves. Eleven S. aureus-infected patients were presented with prosthetic joint infections, 15 with fracture-related infections, 5 with native joint septic arthritis, 15 with diabetic foot infections, and 17 with suspected orthopedic infections in the soft tissue. Anti-S. aureus MENSA IgG levels in patients with non-S. aureus infections and healthy controls were 4-fold (***p = 0.0002) and 8-fold (****p < 0.0001) lower, respectively, compared with those with culture-confirmed S. aureus infections. Comparison of MENSA IgG responses among S. aureus culture-positive patients revealed Ags predictive of active MSKI (IsdB, SCIN, Gmd) and Ags predictive of MSKI type (IsdB, IsdH, Amd, Hla). When combined, IsdB, IsdH, Gmd, Amd, SCIN, and Hla were highly discriminatory of S. aureus MSKI (area under the ROC curve = 0.89 [95% confidence interval 0.82-0.93, p < 0.01]). Collectively, these results demonstrate the feasibility of a bioinformatic approach to use a patient's active immune proteome against S. aureus to diagnose challenging MSKI.

Genetic analysis of osteoblast activity identifies Zbtb40 as a regulator of osteoblast activity and bone mass.

Osteoporosis is a genetic disease characterized by progressive reductions in bone mineral density (BMD) leading to an increased risk of fracture. Over the last decade, genome-wide association studies (GWASs) have identified over 1000 associations for BMD. However, as a phenotype BMD is challenging as bone is a multicellular tissue affected by both local and systemic physiology. Here, we focused on a single component of BMD, osteoblast-mediated bone formation in mice, and identified associations influencing osteoblast activity on mouse Chromosomes (Chrs) 1, 4, and 17. The locus on Chr. 4 was in an intergenic region between Wnt4 and Zbtb40, homologous to a locus for BMD in humans. We tested both Wnt4 and Zbtb40 for a role in osteoblast activity and BMD. Knockdown of Zbtb40, but not Wnt4, in osteoblasts drastically reduced mineralization. Additionally, loss-of-function mouse models for both genes exhibited reduced BMD. Our results highlight that investigating the genetic basis of in vitro osteoblast mineralization can be used to identify genes impacting bone formation and BMD.

Genetic Dissection of Femoral and Tibial Microarchitecture.

Our understanding of the genetic control of bone strength has relied mainly on estimates of bone mineral density. Here we have mapped genetic factors that influence femoral and tibial microarchitecture using high-resolution x-ray computed tomography (8-µm isotropic voxels) across a family of 61 BXD strains of mice, roughly 10 isogenic cases per strain and balanced by sex. We computed heritabilities for 25 cortical and trabecular traits. Males and females have well-matched heritabilities, ranging from 0.25 to 0.75. We mapped 16 genetic loci most of which were detected only in females. There is also a bias in favor of loci that control cortical rather than trabecular bone. To evaluate candidate genes, we combined well-established gene ontologies with bone transcriptome data to compute bone-enrichment scores for all protein-coding genes. We aligned candidates with those of human genome-wide association studies. A subset of 50 strong candidates fell into three categories: (1) experimentally validated genes already known to modulate bone function (Adamts4, Ddr2, Darc, Adam12, Fkbp10, E2f6, Adam17, Grem2, Ifi204); (2) candidates without any experimentally validated function in bone (eg, Greb1, Ifi202b), but linked to skeletal phenotypes in human cohorts; and (3) candidates that have high bone-enrichment scores, but for which there is not yet any functional link to bone biology or skeletal system disease (including Ifi202b, Ly9, Ifi205, Mgmt, F2rl1, Iqgap2). Our results highlight contrasting genetic architecture between sexes and among major bone compartments. The alignment of murine and human data facilitates function analysis and should prove of value for preclinical testing of molecular control of bone structure. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.