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.

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.

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.

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.

Mouse genome-wide association and systems genetics identifies Lhfp as a regulator of bone mass.

Bone mineral density (BMD) is a strong predictor of osteoporotic fracture. It is also one of the most heritable disease-associated quantitative traits. As a result, there has been considerable effort focused on dissecting its genetic basis. Here, we performed a genome-wide association study (GWAS) in a panel of inbred strains to identify associations influencing BMD. This analysis identified a significant (P = 3.1 x 10-12) BMD locus on Chromosome 3@52.5 Mbp that replicated in two separate inbred strain panels and overlapped a BMD quantitative trait locus (QTL) previously identified in a F2 intercross. The association mapped to a 300 Kbp region containing four genes; Gm2447, Gm20750, Cog6, and Lhfp. Further analysis found that Lipoma HMGIC Fusion Partner (Lhfp) was highly expressed in bone and osteoblasts. Furthermore, its expression was regulated by a local expression QTL (eQTL), which overlapped the BMD association. A co-expression network analysis revealed that Lhfp was strongly connected to genes involved in osteoblast differentiation. To directly evaluate its role in bone, Lhfp deficient mice (Lhfp-/-) were created using CRISPR/Cas9. Consistent with genetic and network predictions, bone marrow stromal cells (BMSCs) from Lhfp-/- mice displayed increased osteogenic differentiation. Lhfp-/- mice also had elevated BMD due to increased cortical bone mass. Lastly, we identified SNPs in human LHFP that were associated (P = 1.2 x 10-5) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass.

Mouse Models and Online Resources for Functional Analysis of Osteoporosis Genome-Wide Association Studies.

Osteoporosis is a complex genetic disease in which the number of loci associated with the bone mineral density, a clinical risk factor for fracture, has increased at an exponential rate in the last decade. The identification of the causative variants and candidate genes underlying these loci has not been able to keep pace with the rate of locus discovery. A large number of tools and data resources have been built around the use of the mouse as model of human genetic disease. Herein, we describe resources available for functional validation of human Genome Wide Association Study (GWAS) loci using mouse models. We specifically focus on large-scale phenotyping efforts focused on bone relevant phenotypes and repositories of genotype-phenotype data that exist for transgenic and mutant mice, which can be readily mined as a first step toward more targeted efforts designed to deeply characterize the role of a gene in bone biology.

Characterization of expression and alternative splicing of the gene cadherin-like and PC esterase domain containing 1 (Cped1).

Cadherin-like and PC-esterase domain containing 1 (CPED1) is an uncharacterized gene with no known function. Human genome wide association studies (GWAS) for bone mineral density (BMD) have repeatedly identified a significant locus on Chromosome 7 that contains the gene CPED1, but it remains unclear if this gene could be causative. While an open reading frame for this gene has been predicted, there has been no systematic exploration of expression or alternate splicing for CPED1 in humans or mice.Using mouse models, we demonstrate that Cped1 is alternately spliced whereby transcripts are generated with exon 3 or exons 16 and 17 removed. In calvarial-derived pre-osteoblasts, Cped1 utilizes the predicted promoter upstream of exon 1 as well as alternate promoters upstream of exon 3 and exon 12.Lastly, we have determined that some transcripts terminate at the end of exon 10 and therefore do not contain the cadherin like and the PC esterase domains.Together, these data suggest that multiple protein products may be produced by this gene, with some products either lacking or containing both the predicted functional domains. Our data provide a framework upon which future functional studies will be built to understand the role of this gene in bone biology.

Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis.

Osteoporosis is a common disease diagnosed primarily by measurement of bone mineral density (BMD). We undertook a genome-wide association study (GWAS) in 142,487 individuals from the UK Biobank to identify loci associated with BMD as estimated by quantitative ultrasound of the heel. We identified 307 conditionally independent single-nucleotide polymorphisms (SNPs) that attained genome-wide significance at 203 loci, explaining approximately 12% of the phenotypic variance. These included 153 previously unreported loci, and several rare variants with large effect sizes. To investigate the underlying mechanisms, we undertook (1) bioinformatic, functional genomic annotation and human osteoblast expression studies; (2) gene-function prediction; (3) skeletal phenotyping of 120 knockout mice with deletions of genes adjacent to lead independent SNPs; and (4) analysis of gene expression in mouse osteoblasts, osteocytes and osteoclasts. The results implicate GPC6 as a novel determinant of BMD, and also identify abnormal skeletal phenotypes in knockout mice associated with a further 100 prioritized genes.

A mutagenesis-derived Lrp5 mouse mutant with abnormal retinal vasculature and low bone mineral density.

PURPOSE: Familial exudative vitreoretinopathy (FEVR) is caused by mutations in the genes encoding low-density lipoprotein receptor-related protein (LRP5) or its interacting partners, namely frizzled class receptor 4 (FZD4) and norrin cystine knot growth factor (NDP). Mouse models for Lrp5, Fzd4, and Ndp have proven to be important for understanding the retinal pathophysiology underlying FEVR and systemic abnormalities related to defective Wnt signaling. Here, we report a new mouse mutant, tvrm111B, which was identified by electroretinogram (ERG) screening of mice generated in the Jackson Laboratory Translational Vision Research Models (TVRM) mutagenesis program. METHODS: ERGs were used to examine outer retinal physiology. The retinal vasculature was examined by in vivo retinal imaging, as well as by histology and immunohistochemistry. The tvrm111B locus was identified by genetic mapping of mice generated in a cross to DBA/2J, and subsequent sequencing analysis. Gene expression was examined by real-time PCR of retinal RNA. Bone mineral density (BMD) was examined by peripheral dual-energy X-ray absorptiometry. RESULTS: The tvrm111B allele is inherited as an autosomal recessive trait. Genetic mapping of the decreased ERG b-wave phenotype of tvrm111B mice localized the mutation to a region on chromosome 19 that included Lrp5. Sequencing of Lrp5 identified the insertion of a cytosine (c.4724_4725insC), which is predicted to cause a frameshift that disrupts the last three of five conserved PPPSPxS motifs in the cytoplasmic domain of LRP5, culminating in a premature termination. In addition to a reduced ERG b-wave, Lrp5tvrm111B homozygotes have low BMD and abnormal features of the retinal vasculature that have been reported previously in Lrp5 mutant mice, including persistent hyaloid vessels, leakage on fluorescein angiography, and an absence of the deep retinal capillary bed. CONCLUSIONS: The phenotype of the Lrp5tvrm111B mutant includes abnormalities of the retinal vasculature and of BMD. This model may be a useful resource to further our understanding of the biological role of LRP5 and to evaluate experimental therapies for FEVR or other conditions associated with LRP5 dysfunction.