Bivariate genome-wide association meta-analysis of pediatric musculoskeletal traits reveals pleiotropic effects at the SREBF1/TOM1L2 locus

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated Cadherin-like and PC-esterase Domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

Bone health index in the assessment of bone health: The Generation R Study

Bone Health Index (BHI) has been proposed as a useful instrument for assessing bone health in children. However, its relationship with fracture risk remains unknown. We aimed to investigate whether BHI is associated with bone mineral density (BMD) and prevalent fracture odds in children from the Generation R Study. We also implemented genome-wide association study (GWAS) and polygenic score (PGS) approaches to improve our understanding of BHI and its potential. In total, 4150 children (49.4 % boys; aged 9.8 years) with genotyped data and bone assessments were included in this study. BMD was measured across the total body (less head following ISCD guidelines) using a GE-Lunar iDXA densitometer; and BHI was determined from the hand DXA scans using BoneXpert®. Fractures were self-reported collected with home questionnaires. The association of BHI with BMD and fractures was evaluated using linear models corrected for age, sex, ethnicity, height, and weight. We observed a positive correlation between BHI and BMD (ρ = 0.32, p-value<0.0001). Further, every SD decrease in BHI was associated with an 11 % increased risk of prevalent fractures (OR:1.11, 95 % CI 1.00-1.24, p-value = 0.05). Our BHI GWAS identified variants (lead SNP rs1404264-A, p-value = 2.61 × 10-14) mapping to the ING3/CPED1/WNT16 locus. Children in the extreme tails of the BMD PGS presented a difference in BHI values of -0.10 standard deviations (95% CI -0.14 to -0.07; p-value<0.0001). On top of the demonstrated epidemiological association of BHI with both BMD and fracture risk, our results reveal a partially shared biological background between BHI and BMD. These findings highlight the potential value of using BHI to screen children at risk of fracture.

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis

The precise roles of chromatin organization at osteoporosis risk loci remain largely elusive. Here, we combined chromatin interaction conformation (Hi-C) profiling and self-transcribing active regulatory region sequencing (STARR-seq) to qualify enhancer activities of prioritized osteoporosis-associated single-nucleotide polymorphisms (SNPs). We identified 319 SNPs with biased allelic enhancer activity effect (baaSNPs) that linked to hundreds of candidate target genes through chromatin interactions across 146 loci. Functional characterizations revealed active epigenetic enrichment for baaSNPs and prevailing osteoporosis-relevant regulatory roles for their chromatin interaction genes. Further motif enrichment and network mapping prioritized several putative, key transcription factors (TFs) controlling osteoporosis binding to baaSNPs. Specifically, we selected one top-ranked TF and deciphered that an intronic baaSNP (rs11202530) could allele-preferentially bind to YY2 to augment PAPSS2 expression through chromatin interactions and promote osteoblast differentiation. Our results underline the roles of TF-mediated enhancer-promoter contacts for osteoporosis, which may help to better understand the intricate molecular regulatory mechanisms underlying osteoporosis risk loci.

Evaluating modularity in the hominine skull related to feeding biomechanics

OBJECTIVES

Modular architecture of traits in complex organisms can be important for morphological evolution at micro- and sometimes macroevolutionary scales as it may influence the tempo and direction of changes to groups of traits that are essential for particular functions, including food acquisition and processing. We tested several distinct hypotheses about craniofacial modularity in the hominine skull in relation to feeding biomechanics.

MATERIALS AND METHODS

First, we formulated hypothesized functional modules for craniofacial traits reflecting specific demands of feeding biomechanics (e.g., masseter leverage/gape or tooth crown mechanics) in Homo sapiens, Pan troglodytes, and Gorilla gorilla. Then, the pattern and strength of modular signal was quantified by the covariance ratio coefficient and compared across groups using covariance ratio effect size. Hierarchical clustering analysis was then conducted to examine whether a priori-defined functional modules correspond to empirically recovered clusters.

RESULTS

There was statistical support for most a priori-defined functional modules in the cranium and half of the functional modules in the mandible. Modularity signal was similar in the cranium and mandible, and across the three taxa. Despite a similar strength of modularity, the empirically recovered clusters do not map perfectly onto our priori functional modules, indicating that further work is needed to refine our hypothesized functional modules.

CONCLUSION

The results suggest that modular structure of traits in association with feeding biomechanics were mostly shared with humans and the two African apes. Thus, conserved patterns of functional modularity may have facilitated evolutionary changes to the skull during human evolution.

Genetic and Gene Expression Resources for Osteoporosis and Bone Biology Research

PURPOSE OF REVIEW

The integration of data from multiple genomic assays from humans and non-human model organisms is an effective approach to identify genes involved in skeletal fragility and fracture risk due to osteoporosis and other conditions. This review summarizes genome-wide genetic variation and gene expression data resources relevant to the discovery of genes contributing to skeletal fragility and fracture risk.

RECENT FINDINGS

Genome-wide association studies (GWAS) of osteoporosis-related traits are summarized, in addition to gene expression in bone tissues in humans and non-human organisms, with a focus on rodent models related to skeletal fragility and fracture risk. Gene discovery approaches using these genomic data resources are described. We also describe the Musculoskeletal Knowledge Portal (MSKKP) that integrates much of the available genomic data relevant to fracture risk. The available genomic resources provide a wealth of knowledge and can be analyzed to identify genes related to fracture risk. Genomic resources that would fill particular scientific gaps are discussed.

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.

Identification of four-gene signature to diagnose osteoarthritis through bioinformatics and machine learning methods

BACKGROUND

Although osteoarthritis (OA) is one of the most prevalent joint disorders, effective biomarkers to diagnose OA are still unavailable. This study aimed to acquire some key synovial biomarkers (hub genes) and analyze their correlation with immune infiltration in OA.

METHODS

Gene expression profiles and clinical characteristics of OA and healthy synovial samples were retrieved from the Gene Expression Omnibus (GEO) database. Hub genes for OA were mined based on a combination of weighted gene co-expression network analysis (WGCNA), the least absolute shrinkage and selection operator (LASSO), support vector machine recursive feature elimination (SVM-RFE), and random forest (RF) algorithms. A diagnostic nomogram model for OA prediction was developed based on the hub genes. Receiver operating characteristic curves (ROC) were performed to confirm the abnormal expression of hub genes in the experimemtal and validation datasets. qRT-PCR using patients' samples were conducted as well. In addition, the infiltration level of 28 immune cells in the expression profile and their relationship with hub genes were analyzed using single-sample GSEA (ssGSEA).

RESULTS

4 hub genes (ZBTB16, TNFSF11, SCRG1 and KDELR3) were obtained by WGCNA, lasso, SVM-RFE, RF algorithms as potential biomarkers for OA. The immune infiltration analyses revealed that hub genes were most correlated with regulatory T cell and natural killer cell.

CONCLUSION

A machine learning model to diagnose OA based on ZBTB16, TNFSF11, SCRG1 and KDELR3 using synovial tissue was constructed, providing theoretical foundation and guideline for diagnostic and treatment targets in OA.

Obesity and head and neck cancer risk: a mendelian randomization study

BACKGROUND

Observational studies have reported controversial results on the association between obesity and head and neck cancer risk. This study aimed to perform a two-sample Mendelian randomization (MR) analysis to assess the causal association between obesity and head and neck cancer risk using publicly available genome-wide association studies (GWAS) summary statistics.

METHODS

Single-nucleotide polymorphisms (SNPs) for obesity [body mass index (BMI), waist-to-hip ratio (WHR), whole body fat mass, lean body mass, and trunk fat mass] and head and neck cancer (total head and neck cancer, oral cavity cancer, oropharyngeal cancer, and oral cavity and oropharyngeal cancer) were retrieved from published GWASs and used as genetic instrumental variables. Five methods including inverse-variance-weighted (IVW), weighted-median, MR-Egger, weighted mode, and MR-PRESSO were used to obtain reliable results, and odds ratio with 95% confidence interval (CI) were calculated. Tests for horizontal pleiotropy, heterogeneity, and sensitivity were performed separately.

RESULTS

Genetically predicted BMI was negatively associated with the risk of total head and neck cancer, which was significant in the IVW [OR (95%CI), 0.990 (0.984-0.996), P = 0.0005], weighted-median [OR (95%CI), 0.984 (0.975-0.993), P = 0.0009], and MR-PRESSO [OR (95%CI), 0.990 (0.984-0.995), P = 0.0004] analyses, but suggestive significant in the MR-Egger [OR (95%CI), 0.9980 (0.9968-0.9991), P < 0.001] and weighted mode [OR (95%CI), 0.9980 (0.9968-0.9991), P < 0.001] analyses. Similar, genetically predicted BMI adjust for smoking may also be negatively associated with the risk of total head and neck cancer (P < 0.05). Genetically predicted BMI may be negatively related to the risk of oral cavity cancer, oropharyngeal cancer, and oral cavity and oropharyngeal cancer (P < 0.05), but no causal association was observed for BMI adjust for smoking (P > 0.05). In addition, no causal associations were observed for other exposures and outcomes (all P > 0.05).

CONCLUSION

This MR analysis supported the causal association of BMI-related obesity with decreased risk of total head and neck cancer. However, the effect estimates from the MR analysis were close to 1, suggesting a slight protective effect of BMI-related obesity on head and neck cancer risk.

Multiple Mechanisms Explain Genetic Effects at the CPED1-WNT16 Bone Mineral Density Locus

PURPOSE OF REVIEW

Chromosome region 7q31.31, also known as the CPED1-WNT16 locus, is robustly associated with BMD and fracture risk. The aim of the review is to highlight experimental studies examining the function of genes at the CPED1-WNT16 locus.

RECENT FINDINGS

Genes that reside at the CPED1-WNT16 locus include WNT16, FAM3C, ING3, CPED1, and TSPAN12. Experimental studies in mice strongly support the notion that Wnt16 is necessary for bone mass and strength. In addition, roles for Fam3c and Ing3 in regulating bone morphology in vivo and/or osteoblast differentiation in vitro have been identified. Finally, a role for wnt16 in dually influencing bone and muscle morphogenesis in zebrafish has recently been discovered, which has brought forth new questions related to whether the influence of WNT16 in muscle may conspire with its influence in bone to alter BMD and fracture risk.

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.

Transcriptome-wide association study identifies novel genes associated with bone mineral density and lean body mass in children

OBJECTIVE

To identify novel candidate genes whose expression is associated with bone mineral density (BMD) and body lean mass (LM) in children.

METHODS

A tissue-specific transcriptome-wide association study (TWAS) was conducted utilizing a large-scale genome-wide association study (GWAS) dataset associated with BMD and LM and involving 10,414 participants. The measurement of BMD and LM phenotypes was made based on total-body dual-energy X-ray absorptiometry (TB-DXA) scans. TWAS was conducted by using FUSION software. Reference panels for muscle skeleton (MS), peripheral blood (NBL) and whole blood (YBL) were used for TWAS analysis. Functional enrichment and protein-protein interaction (PPI) analyses of the genes identified by TWAS were performed by using the online tool Metascape ( http://metascape.org ).

RESULTS

For BMD, we identified 174 genes with P < 0.05, such as IKZF1 (P = 1.46 × 10^-9) and CHKB (P = 8.31 × 10^-7). For LM, we identified 208 genes with P < 0.05, such as COPS5 (P = 3.03 × 10^-12) and MRPS33 (P = 5.45 × 10^-10). Gene ontology (GO) enrichment analysis of the BMD-associated genes revealed 200 GO terms, such as protein catabolic process (Log P = -5.09) and steroid hormone-mediated signaling pathway (Log P = -3.13). GO enrichment analysis of the LM-associated genes detected 287 GO terms, such as the apoptotic signaling pathway (Log P = -8.08) and lipid storage (Log P = -3.55).

CONCLUSION

This study identified several candidate genes for BMD and LM in children, providing novel clues to the genetic mechanisms underlying the development of childhood BMD and LM.

Causal effects of B vitamins and homocysteine on obesity and musculoskeletal diseases: A Mendelian randomization study

Objectives

Although homocysteine (Hcy) increases the risk of cardiovascular diseases, its effects on obesity and musculoskeletal diseases remain unclear. We performed a Mendelian randomization study to estimate the associations between Hcy and B vitamin concentrations and their effects on obesity and musculoskeletal-relevant diseases in the general population.

Methods

We selected independent single nucleotide polymorphisms of Hcy (n = 44,147), vitamin B12 (n = 45,576), vitamin B6 (n = 1864), and folate (n = 37,465) at the genome-wide significance level as instruments and applied them to the studies of summary-level data for fat and musculoskeletal phenotypes from the UK Biobank study (n = 331,117), the FinnGen consortium (n = 218,792), and other consortia. Two-sample Mendelian randomization (MR) approaches were utilized in this study. The inverse variance weighting (IVW) was adopted as the main analysis. MR-PRESSO, MR-Egger, the weighted median estimate, bidirectional MR, and multivariable MR were performed as sensitivity methods.

Results

Higher Hcy concentrations were robustly associated with an increased risk of knee osteoarthritis [odds ratio (OR) 1.119; 95% confidence interval (CI) 1.032-1.214; P = 0.007], hospital-diagnosed osteoarthritis (OR 1.178; 95% CI 1.012-1.37; P = 0.034), osteoporosis with pathological fracture (OR 1.597; 95% CI 1.036-2.46; P = 0.034), and soft tissue disorder (OR 1.069; 95% CI 1.001-1.141; P = 0.045) via an inverse variance weighting method and other MR approaches. Higher vitamin B12 levels were robustly associated with decreased body fat percentage and its subtypes (all P < 0.05). Bidirectional analyses showed no reverse causation. Multivariable MR analyses and other sensitivity analyses showed directionally similar results.

Conclusions

There exist significant causal effects of vitamin B12 in the serum and Hcy in the blood on fat and musculoskeletal diseases, respectively. These findings may have an important insight into the pathogenesis of obesity and musculoskeletal diseases and other possible future therapies.

Simultaneous detection of novel genes and SNPs by adaptive p-value combination

Combining SNP p-values from GWAS summary data is a promising strategy for detecting novel genetic factors. Existing statistical methods for the p-value-based SNP-set testing confront two challenges. First, the statistical power of different methods depends on unknown patterns of genetic effects that could drastically vary over different SNP sets. Second, they do not identify which SNPs primarily contribute to the global association of the whole set. We propose a new signal-adaptive analysis pipeline to address these challenges using the omnibus thresholding Fisher’s method (oTFisher). The oTFisher remains robustly powerful over various patterns of genetic effects. Its adaptive thresholding can be applied to estimate important SNPs contributing to the overall significance of the given SNP set. We develop efficient calculation algorithms to control the type I error rate, which accounts for the linkage disequilibrium among SNPs. Extensive simulations show that the oTFisher has robustly high power and provides a higher balanced accuracy in screening SNPs than the traditional Bonferroni and FDR procedures. We applied the oTFisher to study the genetic association of genes and haplotype blocks of the bone density-related traits using the summary data of the Genetic Factors for Osteoporosis Consortium. The oTFisher identified more novel and literature-reported genetic factors than existing p-value combination methods. Relevant computation has been implemented into the R package TFisher to support similar data analysis.

wnt16 regulates spine and muscle morphogenesis through parallel signals from notochord and dermomyotome

Bone and muscle are coupled through developmental, mechanical, paracrine, and autocrine signals. Genetic variants at the CPED1-WNT16 locus are dually associated with bone- and muscle-related traits. While Wnt16 is necessary for bone mass and strength, this fails to explain pleiotropy at this locus. Here, we show wnt16 is required for spine and muscle morphogenesis in zebrafish. In embryos, wnt16 is expressed in dermomyotome and developing notochord, and contributes to larval myotome morphology and notochord elongation. Later, wnt16 is expressed at the ventral midline of the notochord sheath, and contributes to spine mineralization and osteoblast recruitment. Morphological changes in wnt16 mutant larvae are mirrored in adults, indicating that wnt16 impacts bone and muscle morphology throughout the lifespan. Finally, we show that wnt16 is a gene of major effect on lean mass at the CPED1-WNT16 locus. Our findings indicate that Wnt16 is secreted in structures adjacent to developing bone (notochord) and muscle (dermomyotome) where it affects the morphogenesis of each tissue, thereby rendering wnt16 expression into dual effects on bone and muscle morphology. This work expands our understanding of wnt16 in musculoskeletal development and supports the potential for variants to act through WNT16 to influence bone and muscle via parallel morphogenetic processes.

In humans, genetic variants (DNA sequences that vary amongst individuals) have been identified that appear to influence two tissues, bone and skeletal muscle. However, how single genes and genetic variants exert dual influence on both tissues is not well understood. In this study, we found that the wnt16 gene is necessary for specifying the size and shape of both muscle and bone during development in zebrafish. We also disentangled how wnt16 affects both tissues: distinct cellular populations adjacent to muscle and bone secrete Wnt16, where it acts as a signal guiding the size and shape of each tissue. This is important because in humans, genetic variants near the WNT16 gene have effects on both bone- and muscle-related traits. This study expands our understanding of the role of WNT16 in bone and muscle development, and helps to explain how genetic variants near WNT16 affect traits for both tissues. Moreover, WNT16 is actively being explored as a target for osteoporosis therapies, thus our study could have implications with regard to the potential of targeting WNT16 to treat bone and muscle simultaneously.

Discovery of exercise-related genes and pathway analysis based on comparative genomes of Mongolian originated Abaga and Wushen horse

The Mongolian horses have excellent endurance and stress resistance to adapt to the cold and harsh plateau conditions. Intraspecific genetic diversity is mainly embodied in various genetic advantages of different branches of the Mongolian horse. Since people pay progressive attention to the athletic performance of horse, we expect to guide the exercise-oriented breeding of horses through genomics research. We obtained the clean data of 630,535,376,400 bp through the entire genome second-generation sequencing for the whole blood of four Abaga horses and ten Wushen horses. Based on the data analysis of single nucleotide polymorphism, we severally detected that 479 and 943 positively selected genes, particularly exercise related, were mainly enriched on equine chromosome 4 in Abaga horses and Wushen horses, which implied that chromosome 4 may be associated with the evolution of the Mongolian horse and athletic performance. Four hundred and forty genes of positive selection were enriched in 12 exercise-related pathways and narrowed in 21 exercise-related genes in Abaga horse, which were distinguished from Wushen horse. So, we speculated that the Abaga horse may have oriented genes for the motorial mechanism and 21 exercise-related genes also provided a molecular genetic basis for exercise-directed breeding of the Mongolian horse.

Identifying Causes of Fracture Beyond Bone Mineral Density: Evidence From Human Genetics

New therapies may help to prevent osteoporotic fractures other than through increasing bone mineral density (BMD). Because fracture risk has an important genetic component, we aim to identify loci increasing fracture risk that do not decrease BMD, using a recently-proposed structural equation model adapted to remove genetic influences of BMD on fracture risk. We used summary statistics of the largest genome-wide association studies (GWASs) for BMD and for fracture in these analyses. We next estimated the genetic correlation between the non-BMD or BMD-related genetic effects and other clinical risk factors for fracture. Last, based on white British participants in the UK Biobank, we conducted genetic risk score analyses to assess whether the aggregated genetic effects conferred increased major osteoporotic fracture risk. We found that only three loci affecting fracture risk exhibited genetic effects not mediated by BMD: SOST, CPED1-WNT16, and RSPO3, while these three loci simultaneously conferred BMD-related effects. No strong genetic associations between non-BMD or BMD-related effects and 16 clinical risk factors were observed. However, non-BMD effects might be genetic correlated with hip bone size. In the UK Biobank, a 1 standard deviation (1-SD) increase in the non-BMD genetic risk score conferred an odds ratio of 1.17 for incident major osteoporotic fracture, compared to 1.29 by a BMD-related genetic risk score. Our study suggests that the majority of common genetic predisposition toward fracture risk acts upon BMD. Although non-BMD genetic effects may exist, they are not strongly correlated with most traditional clinical risk factors. Risk loci harboring non-BMD genetic effects may influence other perspectives of bone quality, or confer effects that existing GWASs fail to capture, but they demonstrate weaker impact on fracture risk than BMD-related genetic effects. These findings suggest that most successful drug development programs for osteoporosis should focus on pathways identified through BMD-associated loci. © 2022 American Society for Bone and Mineral Research (ASBMR).

Ras and Rab Interactor 3: From Cellular Mechanisms to Human Diseases

Ras and Rab interactor 3 (RIN3) functions as a Guanine nucleotide Exchange Factor (GEF) for some members of the Rab family of small GTPase. By promoting the activation of Rab5, RIN3 plays an important role in regulating endocytosis and endocytic trafficking. In addition, RIN3 activates Ras, another small GTPase, that controls multiple signaling pathways to regulate cellular function. Increasing evidence suggests that dysregulation of RIN3 activity may contribute to the pathogenesis of several disease conditions ranging from Paget’s Disease of the Bone (PDB), Alzheimer’s Disease (AD), Chronic Obstructive Pulmonary Disease (COPD) and to obesity. Recent genome-wide association studies (GWAS) identified variants in the RIN3 gene to be linked with these disease conditions. Interestingly, some variants appear to be missense mutations in the functional domains of the RIN3 protein while most variants are located in the noncoding regions of the RIN3 gene, potentially altering its gene expression. However, neither the protein structure of RIN3 nor its exact function(s) (except for its GEF activity) has been fully defined. Furthermore, how the polymorphisms/variants contribute to disease pathogenesis remain to be understood. Herein, we examine, and review published studies in an attempt to provide a better understanding of the physiological function of RIN3; More importantly, we construct a framework linking the polymorphisms/variants of RIN3 to altered cell signaling and endocytic traffic, and to potential disease mechanism(s).

Functional Validation of Osteoporosis Genetic Findings Using Small Fish Models

The advancement of human genomics has revolutionized our understanding of the genetic architecture of many skeletal diseases, including osteoporosis. However, interpreting results from human association studies remains a challenge, since index variants often reside in non-coding regions of the genome and do not possess an obvious regulatory function. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary, such as the one offered by animal models. These models enable us to identify causal mechanisms, clarify the underlying biology, and apply interventions. Over the past several decades, small teleost fishes, mostly zebrafish and medaka, have emerged as powerful systems for modeling the genetics of human diseases. Due to their amenability to genetic intervention and the highly conserved genetic and physiological features, fish have become indispensable for skeletal genomic studies. The goal of this review is to summarize the evidence supporting the utility of Zebrafish (Danio rerio) for accelerating our understanding of human skeletal genomics and outlining the remaining gaps in knowledge. We provide an overview of zebrafish skeletal morphophysiology and gene homology, shedding light on the advantages of human skeletal genomic exploration and validation. Knowledge of the biology underlying osteoporosis through animal models will lead to the translation into new, better and more effective therapeutic approaches.

The genetic overlap between osteoporosis and craniosynostosis

Osteoporosis is the most prevalent bone condition in the ageing population. This systemic disease is characterized by microarchitectural deterioration of bone, leading to increased fracture risk. In the past 15 years, genome-wide association studies (GWAS), have pinpointed hundreds of loci associated with bone mineral density (BMD), helping elucidate the underlying molecular mechanisms and genetic architecture of fracture risk. However, the challenge remains in pinpointing causative genes driving GWAS signals as a pivotal step to drawing the translational therapeutic roadmap. Recently, a skull BMD-GWAS uncovered an intriguing intersection with craniosynostosis, a congenital anomaly due to premature suture fusion in the skull. Here, we recapitulate the genetic contribution to both osteoporosis and craniosynostosis, describing the biological underpinnings of this overlap and using zebrafish models to leverage the functional investigation of genes associated with skull development and systemic skeletal homeostasis.

Identification of PDXDC1 as a novel pleiotropic susceptibility locus shared between lumbar spine bone mineral density and birth weight

An increasing number of epidemiological studies have suggested that birth weight (BW) may be a determinant of bone health later in life, although the underlying genetic mechanism remains unclear. Here, we applied a pleiotropic conditional false discovery rate (cFDR) approach to the genome-wide association study (GWAS) summary statistics for lumbar spine bone mineral density (LS BMD) and BW, aiming to identify novel susceptibility variants shared between these two traits. We detected 5 novel potential pleiotropic loci which are located at or near 7 different genes (NTAN1, PDXDC1, CACNA1G, JAG1, FAT1P1, CCDC170, ESR1), among which PDXDC1 and FAT1P1 have not previously been linked to these phenotypes. To partially validate the findings, we demonstrated that the expression of PDXDC1 was dramatically reduced in ovariectomized (OVX) mice in comparison with sham-operated (SHAM) mice in both the growth plate and trabecula bone. Furthermore, immunohistochemistry assay with serial sections showed that both osteoclasts and osteoblasts express PDXDC1, supporting its potential role in bone metabolism. In conclusion, our study provides insights into some shared genetic mechanisms for BMD and BW as well as a novel potential therapeutic target for the prevention of OP in the early stages of the disease development. KEY MESSAGES : We investigated pleiotropy-informed enrichment between LS BMD and BW. We identified genetic variants related to both LS BMD and BW by utilizing a cFDR approach. PDXDC1 is a novel pleiotropic gene which may be related to both LS BMD and BW. Elevated expression of PDXDC1 is related to higher BMD and lower ratio n-6/n-3 PUFA indicating a bone protective effect of PDXDC1.

Integrative analysis of genomic and epigenomic data reveals underlying super-enhancer-mediated microRNA regulatory network for human bone mineral density

Bone mineral density (BMD) is a highly heritable complex trait and is a key indicator for diagnosis and treatment for osteoporosis. In the last decade, numerous susceptibility loci for BMD and fracture have been identified by genome wide association studies (GWAS); however, fine mapping of these loci is challengeable. Here, we proposed a new long-range fine-mapping approach that combined super-enhancers (SEs) and microRNAs (miRNAs) data, which were two important factors in control of cell identity and specific differentiation, with the GWAS summary datasets in cell-type-restricted way. Genome-wide SE-based analysis found that the BMD-related variants were significantly enriched in the osteoblast SE regions, indicative of potential long-range effects of such SNPs. With the SNP-mapped SEs (mSEs), 13 accessible long-range mSE-interacted miRNAs (mSE-miRNAs) were identified by integrating osteoblast Hi-C and ATAC-seq data, including three known bone-related miRNAs (miR-132-3p, miR-212-3p and miR-125b-5p). The putative targets of the two newly identified mSE-miRNAs (miR-548aj-3p and miR-190a-3p) were found largely enriched in osteogenic-related pathway and processes, suggesting that these mSE-miRNAs could be functional in the regulation of osteoblast differentiation. Furthermore, we identified 54 genes with the long-range 'mSE-miRNA' approach, and 24 of them were previously reported to be related to skeletal development. Besides, enrichment analysis found that these genes were specifically enriched in the post-transcriptional regulation and bone formation processes. This study provided a new insight into the approach of fine-mapping of GWAS loci. A tool was provided for the genome-wide SE-based analysis and the detection of long-range osteoblast-restricted mSE-miRNAs (https://github.com/Zheng-Lab-Westlake/Osteo-Fine-Mapp-SNP2SE2miRNA).

Identification of novel pleiotropic gene for bone mineral density and lean mass using the cFDR method

Bone mineral density (BMD) and whole-body lean mass (WBLM) are two important phenotypes of osteoporosis and sarcopenia. Previous studies have shown that BMD and lean mass were phenotypically and genetically correlated. To identify the novel common genetic factors shared between BMD and WBLM, we performed the conditional false discovery rate (cFDR) analysis using summary data of the genome-wide association study of femoral neck BMD (n = 53,236) and WBLM (n = 38,292) from the Genetic Factors for Osteoporosis Consortium (GEFOS). We identified eight pleiotropic Single Nucleotide Polymorphism (SNPs) (PLCL1 rs11684176 and rs2880389, JAZF1 rs198, ADAMTSL3 rs10906982, RFTN2/MARS2 rs7340470, SH3GL3 rs1896797, ST7L rs10776755, ANKRD44/SF3B1 rs11888760) significantly associated with femoral neck BMD and WBLM (ccFDR < 0.05). Bayesian fine-mapping analysis showed that rs11888760, rs198, and rs1896797 were the possible functional variants in the ANKRD44/SF3B1, JAZF1i, and SH3GL3 loci, respectively. Functional annotation suggested that rs11888760 was likely to comprise a DNA regulatory element and linked to the expression of RFTN2 and PLCL1. PLCL1 showed differential expression in laryngeal posterior cricoarytenoid muscle between rats of 6 months and 30 months of age. Our findings, together with PLCL1's potential functional relevance to bone and skeletal muscle function, suggested that rs11888760 was the possible pleiotropic functional variants appearing to coregulate both bone and muscle metabolism through regulating the expression of PLCL1. The findings enhanced our knowledge of genetic associations between BMD and lean mass and provide a rationale for subsequent functional studies of the implicated genes in the pathophysiology of diseases, such as osteoporosis and sarcopenia.

Relationship between Serum Nutritional Factors and Bone Mineral Density: A Mendelian Randomization Study

PURPOSE

Multiple risk factors have been implicated in the development of osteoporosis. This study examined potential associations between serum nutritional factors and bone mineral density (BMD).

METHODS

Six nutritional factors were selected as exposures. Outcomes included total body BMD (n = 66 945); BMD at the forearm (FA), femoral neck (FN) and lumbar spine (LS) (n = 8143, n = 32 735, and n = 28 498, respectively); estimated heel BMD (HL eBMD) (n = 394 929); and HL eBMD stratified by sex (n = 206 496). A 2-sample Mendelian randomization approach was adopted to estimate the association between serum nutritional factors and BMD. The threshold for adjusted P value was 1.39 × 10-3.

RESULTS

Serum calcium levels were inversely associated with LS BMD (effect = -0.55; 95% CI, -0.86 to -0.24; P = 0.001), whereas serum selenium levels were positively correlated with HL eBMD (effect = 0.22; 95% CI, 0.10 to 0.33; P = 1.70 × 10-4). Regarding nominal significance, there was a positive association between serum selenium levels and FA BMD. Nominally significant results were also obtained for serum retinol as well as vitamin E levels and HL eBMD. Moreover, sex-specific effects of serum retinol and vitamin E levels on BMD were observed in men.

CONCLUSION

Serum calcium and selenium levels influence BMD at specific skeletal sites. This implies that these nutritional factors play crucial roles in bone metabolism.

Enhanced Identification of Novel Potential Variants for Appendicular Lean Mass by Leveraging Pleiotropy With Bone Mineral Density

Strong relationships have been found between appendicular lean mass (ALM) and bone mineral density (BMD). It may be due to a shared genetic basis, termed pleiotropy. By leveraging the pleiotropy with BMD, the aim of this study was to detect more potential genetic variants for ALM. Using the conditional false discovery rate (cFDR) methodology, a combined analysis of the summary statistics of two large independent genome wide association studies (GWAS) of ALM (n = 73,420) and BMD (n = 10,414) was conducted. Strong pleiotropic enrichment and 26 novel potential pleiotropic SNPs were found for ALM and BMD. We identified 156 SNPs for ALM (cFDR <0.05), of which 74 were replicates of previous GWASs and 82 were novel SNPs potentially-associated with ALM. Eleven genes annotated by 31 novel SNPs (13 pleiotropic and 18 ALM specific) were partially validated in a gene expression assay. Functional enrichment analysis indicated that genes corresponding to the novel potential SNPs were enriched in GO terms and/or KEGG pathways that played important roles in muscle development and/or BMD metabolism (adjP <0.05). In protein–protein interaction analysis, rich interactions were demonstrated among the proteins produced by the corresponding genes. In conclusion, the present study, as in other recent studies we have conducted, demonstrated superior efficiency and reliability of the cFDR methodology for enhanced detection of trait-associated genetic variants. Our findings shed novel insight into the genetic variability of ALM in addition to the shared genetic basis underlying ALM and BMD.

Muscle-Bone Crosstalk in Chronic Kidney Disease: The Potential Modulatory Effects of Exercise

Chronic kidney disease (CKD) is a prevalent worldwide public burden that increasingly compromises overall health as the disease progresses. Two of the most negatively affected tissues are bone and skeletal muscle, with CKD negatively impacting their structure, function and activity, impairing the quality of life of these patients and contributing to morbidity and mortality. Whereas skeletal health in this population has conventionally been associated with bone and mineral disorders, sarcopenia has been observed to impact skeletal muscle health in CKD. Indeed, bone and muscle tissues are linked anatomically and physiologically, and together regulate functional and metabolic mechanisms. With the initial crosstalk between the skeleton and muscle proposed to explain bone formation through muscle contraction, it is now understood that this communication occurs through the interaction of myokines and osteokines, with the skeletal muscle secretome playing a pivotal role in the regulation of bone activity. Regular exercise has been reported to be beneficial to overall health. Also, the positive regulatory effect that exercise has been proposed to have on bone and muscle anatomical, functional, and metabolic activity has led to the proposal of regular physical exercise as a therapeutic strategy for muscle and bone-related disorders. The detection of bone- and muscle-derived cytokine secretion following physical exercise has strengthened the idea of a cross communication between these organs. Hence, this review presents an overview of the impact of CKD in bone and skeletal muscle, and narrates how these tissues intrinsically communicate with each other, with focus on the potential effect of exercise in the modulation of this intercommunication.

Variants in NEB and RIF1 genes on chr2q23 are associated with skeletal muscle index in Koreans: genome-wide association study

Although skeletal muscle plays a crucial role in metabolism and influences aging and chronic diseases, little is known about the genetic variations with skeletal muscle, especially in the Asian population. We performed a genome-wide association study in 2,046 participants drawn from a population-based study. Appendicular skeletal muscle mass was estimated based on appendicular lean soft tissue measured with a multi-frequency bioelectrical impedance analyzer and divided by height squared to derive the skeletal muscle index (SMI). After conducting quality control and imputing the genotypes, we analyzed 6,391,983 autosomal SNPs. A genome-wide significant association was found for the intronic variant rs138684936 in the NEB and RIF1 genes (β = 0.217, p = 6.83 × 10^–9). These two genes are next to each other and are partially overlapped on chr2q23. We conducted extensive functional annotations to gain insight into the directional biological implication of significant genetic variants. A gene-based analysis identified the significant TNFSF9 gene and confirmed the suggestive association of the NEB gene. Pathway analyses showed the significant association of regulation of multicellular organism growth gene-set and the suggestive associations of pathways related to skeletal system development or skeleton morphogenesis with SMI. In conclusion, we identified a new genetic locus on chromosome 2 for SMI with genome-wide significance. These results enhance the biological understanding of skeletal muscle mass and provide specific leads for functional experiments.

Identification of pleiotropic loci underlying hip bone mineral density and trunk lean mass

Bone mineral density (BMD) and lean body mass (LBM) not only have a considerable heritability each, but also are genetically correlated. However, common genetic determinants shared by both traits are largely unknown. In the present study, we performed a bivariate genome-wide association study (GWAS) meta-analysis of hip BMD and trunk lean mass (TLM) in 11,335 subjects from 6 samples, and performed replication in estimated heel BMD and TLM in 215,234 UK Biobank (UKB) participants. We identified 2 loci that nearly attained the genome-wide significance (GWS, p < 5.0 × 10-8) level in the discovery GWAS meta-analysis and that were successfully replicated in the UKB sample: 11p15.2 (lead SNP rs12800228, discovery p = 2.88 × 10-7, replication p = 1.95 × 10-4) and 18q21.32 (rs489693, discovery p = 1.67 × 10-7, replication p = 1.17 × 10-3). The above 2 pleiotropic loci may play a pleiotropic role for hip BMD and TLM development. So our findings provide useful insights that further enhance our understanding of genetic interplay between BMD and LBM.

Deletion of SREBF1, a Functional Bone-Muscle Pleiotropic Gene, Alters Bone Density and Lipid Signaling in Zebrafish

Through a genome-wide analysis of bone mineral density (BMD) and muscle mass, identification of a signaling pattern on 17p11.2 recognized the presence of sterol regulatory element-binding factor 1 (SREBF1), a gene responsible for the regulation of lipid homeostasis. In conjunction with lipid-based metabolic functions, SREBF1 also codes for the protein, SREBP-1, a transcription factor known for its role in adipocyte differentiation. We conducted a quantitative correlational study. We established a zebrafish (ZF) SREBF1 knockout (KO) model and used a targeted customized lipidomics approach to analyze the extent of SREBF1 capabilities. For lipidomics profiling, we isolated the dorsal muscles of wild type (WT) and KO fishes, and we performed liquid chromatography-tandem mass spectrometry screening assays of these samples. In our analysis, we profiled 48 lipid mediators (LMs) derived from various essential polyunsaturated fatty acids to determine potential targets regulated by SREBF1, and we found that the levels of 11,12 epoxyeicosatrienoic acid (11,12-EET) were negatively associated with the number of SREBF1 alleles (P = 0.006 for a linear model). We also compared gene expression between KO and WT ZF by genome-wide RNA-sequencing. Significantly enriched pathways included fatty acid elongation, linoleic acid metabolism, arachidonic acid metabolism, adipocytokine signaling, and DNA replication. We discovered trends indicating that BMD in adult fish was significantly lower in the KO than in the WT population (P < 0.03). These studies reinforce the importance of lipidomics investigation by detailing how the KO of SREBF1 affects both BMD and lipid-signaling mediators, thus confirming the importance of SREBF1 for musculoskeletal homeostasis.

Three pleiotropic loci associated with bone mineral density and lean body mass

Both bone mineral density (BMD) and lean body mass (LBM) are important physiological measures with strong genetic determination. Besides, BMD and LBM might have common genetic factors. Aiming to identify pleiotropic genomic loci underlying BMD and LBM, we performed bivariate genome-wide association study meta-analyses of femoral neck bone mineral density and LBM at arms and legs, and replicated in the large-scale UK Biobank cohort sample. Combining the results from discovery meta-analysis and replication sample, we identified three genomic loci at the genome-wide significance level (p < 5.0 × 10-8): 2p23.2 (lead SNP rs4477866, discovery p = 3.47 × 10-8, replication p = 1.03 × 10-4), 16q12.2 (rs1421085, discovery p = 2.04 × 10-9, replication p = 6.47 × 10-14) and 18q21.32 (rs11152213, discovery p = 3.47 × 10-8, replication p = 6.69 × 10-6). Our findings not only provide useful insights into lean mass and bone mass development, but also enhance our understanding of the potential genetic correlation between BMD and LBM.

No Evidence for Passive Gene-Environment Correlation or the Influence of Genetic Risk for Psychiatric Disorders on Adult Body Composition via the Adoption Design

The relationship between genetic and environmental risk is complex and for many traits, estimates of genetic effects may be inflated by passive gene-environment correlation. This arises because biological offspring inherit both their genotypes and rearing environment from their parents. We tested for passive gene-environment correlation in adult body composition traits using the 'natural experiment' of childhood adoption, which removes passive gene-environment correlation within families. Specifically, we compared 6165 adoptees with propensity score matched non-adoptees in the UK Biobank. We also tested whether passive gene-environment correlation inflates the association between psychiatric genetic risk and body composition. We found no evidence for inflation of heritability or polygenic scores in non-adoptees compared to adoptees for a range of body composition traits. Furthermore, polygenic risk scores for anorexia nervosa, attention-deficit/hyperactivity disorder and schizophrenia did not differ in their influence on body composition traits in adoptees and non-adoptees. These findings suggest that passive gene-environment correlation does not inflate genetic effects for body composition, or the influence of psychiatric disorder genetic risk on body composition. Our design does not look at passive gene-environment correlation in childhood, and does not test for 'pure' environmental effects or the effects of active and evocative gene-environment correlations, where child genetics directly influences home environment. However, these findings suggest that genetic influences identified for body composition in this adult sample are direct, and not confounded by the family environment provided by biological relatives.

A powerful method for pleiotropic analysis under composite null hypothesis identifies novel shared loci between Type 2 Diabetes and Prostate Cancer

There is increasing evidence that pleiotropy, the association of multiple traits with the same genetic variants/loci, is a very common phenomenon. Cross-phenotype association tests are often used to jointly analyze multiple traits from a genome-wide association study (GWAS). The underlying methods, however, are often designed to test the global null hypothesis that there is no association of a genetic variant with any of the traits, the rejection of which does not implicate pleiotropy. In this article, we propose a new statistical approach, PLACO, for specifically detecting pleiotropic loci between two traits by considering an underlying composite null hypothesis that a variant is associated with none or only one of the traits. We propose testing the null hypothesis based on the product of the Z-statistics of the genetic variants across two studies and derive a null distribution of the test statistic in the form of a mixture distribution that allows for fractions of variants to be associated with none or only one of the traits. We borrow approaches from the statistical literature on mediation analysis that allow asymptotic approximation of the null distribution avoiding estimation of nuisance parameters related to mixture proportions and variance components. Simulation studies demonstrate that the proposed method can maintain type I error and can achieve major power gain over alternative simpler methods that are typically used for testing pleiotropy. PLACO allows correlation in summary statistics between studies that may arise due to sharing of controls between disease traits. Application of PLACO to publicly available summary data from two large case-control GWAS of Type 2 Diabetes and of Prostate Cancer implicated a number of novel shared genetic regions: 3q23 (ZBTB38), 6q25.3 (RGS17), 9p22.1 (HAUS6), 9p13.3 (UBAP2), 11p11.2 (RAPSN), 14q12 (AKAP6), 15q15 (KNL1) and 18q23 (ZNF236).

We propose a new approach PLACO that uses aggregate-level genotype-phenotype association statistics—commonly referred to as GWAS summary statistics—to identify genetic variants that influence risk of two traits or diseases. It allows correlation in summary statistics between studies that may arise due to sharing of controls between disease traits. We demonstrate that PLACO can achieve major power gain over alternative methods that are typically used. We applied PLACO to Type 2 Diabetes and Prostate Cancer summary data from two large case-control studies. Many previous studies have reported an inverse association of these two chronic diseases suggesting shared risk factors; however, shared genetic mechanisms underlying this association is poorly understood. PLACO identified a number of novel shared genetic regions that are not detected by individual trait analysis. Many of the loci implicated by PLACO increase risk for one disease while decreasing risk for the other. PLACO can similarly be used on other traits to shed light on shared genetic risk factors.

The genetic architecture of appendicular lean mass characterized by association analysis in the UK Biobank study

Appendicular lean mass (ALM) is a heritable trait associated with loss of lean muscle mass and strength, or sarcopenia, but its genetic determinants are largely unknown. Here we conducted a genome-wide association study (GWAS) with 450,243 UK Biobank participants to uncover its genetic architecture. A total of 1059 conditionally independent variants from 799 loci were identified at the genome-wide significance level (p < 5 × 10-9), all of which were also significant at p < 5 × 10-5 in both sexes. These variants explained ~15.5% of the phenotypic variance, accounting for more than one quarter of the total ~50% GWAS-attributable heritability. There was no difference in genetic effect between sexes or among different age strata. Heritability was enriched in certain functional categories, such as conserved and coding regions, and in tissues related to the musculoskeletal system. Polygenic risk score prediction well distinguished participants with high and low ALM. The findings are important not only for lean mass but also for other complex diseases, such as type 2 diabetes, as ALM is shown to be a protective factor for type 2 diabetes.

Helicobacter pylori Related Diseases and Osteoporotic Fractures (Narrative Review)

Osteoporosis (OP) and osteoporotic fractures (OFs) are common multifactorial and heterogenic disorders of increasing incidence. Helicobacter pylori (H.p.) colonizes the stomach approximately in half of the world's population, causes gastroduodenal diseases and is prevalent in numerous extra-digestive diseases known to be associated with OP/OF. The studies regarding relationship between H.p. infection (HPI) and OP/OFs are inconsistent. The current review summarizes the relevant literature on the potential role of HPI in OP, falls and OFs and highlights the reasons for controversies in the publications. In the first section, after a brief overview of HPI biological features, we analyze the studies evaluating the association of HPI and bone status. The second part includes data on the prevalence of OP/OFs in HPI-induced gastroduodenal diseases (peptic ulcer, chronic/atrophic gastritis and cancer) and the effects of acid-suppressive drugs. In the next section, we discuss the possible contribution of HPI-associated extra-digestive diseases and medications to OP/OF, focusing on conditions affecting both bone homeostasis and predisposing to falls. In the last section, we describe clinical implications of accumulated data on HPI as a co-factor of OP/OF and present a feasible five-step algorithm for OP/OF risk assessment and management in regard to HPI, emphasizing the importance of an integrative (but differentiated) holistic approach. Increased awareness about the consequences of HPI linked to OP/OF can aid early detection and management. Further research on the HPI-OP/OF relationship is needed to close current knowledge gaps and improve clinical management of both OP/OF and HPI-related disorders.

Genetic Pleiotropy of Bone-Related Phenotypes: Insights from Osteoporosis

PURPOSE OF REVIEW

We summarize recent evidence on the shared genetics within and outside the musculoskeletal system (mostly related to bone density and osteoporosis).

RECENT FINDINGS

Osteoporosis is determined by an interplay between multiple genetic and environmental factors. Significant progress has been made regarding its genetic background revealing a number of robustly validated loci and respective pathways. However, pleiotropic factors affecting bone and other tissues are not well understood. The analytical methods proposed to test for potential associations between genetic variants and multiple phenotypes can be applied to bone-related data. A number of recent genetic studies have shown evidence of pleiotropy between bone density and other different phenotypes (traits, conditions, or diseases), within and outside the musculoskeletal system. Power benefits of combining correlated phenotypes, as well as unbiased discovery, make these studies promising. Studies in humans are supported by evidence from animal models. Drug development and repurposing should benefit from the pleiotropic approach. We believe that future studies should take into account shared genetics between the bone and related traits.

Bone-Muscle Mutual Interactions

PURPOSE OF REVIEW

The purpose of this review is to describe the current state of our thinking regarding bone-muscle interactions beyond the mechanical perspective.

RECENT FINDINGS

Recent and prior evidence has begun to dissect many of the molecular mechanisms that bone and muscle use to communicate with each other and to modify each other's function. Several signaling factors produced by muscle and bone have emerged as potential mediators of these biochemical/molecular interactions. These include muscle factors such as myostatin, Irisin, BAIBA, IL-6, and the IGF family and the bone factors FGF-23, Wnt1 and Wnt3a, PGE2, FGF9, RANKL, osteocalcin, and sclerostin. The identification of these signaling molecules and their underlying mechanisms offers the very real and exciting possibility that new pharmaceutical approaches can be developed that will permit the simultaneous treatments of diseases that often occur in combination, such as osteoporosis and sarcopenia.

Muscle, Bone, and Fat Crosstalk: the Biological Role of Myokines, Osteokines, and Adipokines

PURPOSE OF REVIEW

Skeletal muscle and bone are connected anatomically and physiologically, and play a crucial role in human locomotion and metabolism. Historically, the coupling between muscle and bone has been viewed in light of mechanotransduction, which dictates that the mechanical forces applied to muscle are transmitted to the skeleton to initiate bone formation. However, these organs also communicate through the endocrine system, orchestrated by a family of cytokines namely myokines (derived from myocytes) and osteokines (derived from bone cells). A third player in this biochemical crosstalk is adipose tissue and the secretion of adipokines (derived from adipocytes). In this review, we discuss the bidirectional effects of myokines and osteokines on muscle and bone metabolism, and the impact of adipokines on both of these secretory organs.

RECENT FINDINGS

Several myokines, notably, IL6, irisin, IGF-1, BDNF, myostatin, and FGF2 exert anabolic/catabolic effects on bone, while the osteokines osteocalcin and sclerostin have shown to induce muscle anabolism and catabolism, respectively. Adipokines, such as leptin, resistin, adiponectin, and TNFα (released from adipose tissue), can also modulate muscle and bone metabolism. Contrarily, exercise-mediated release of lipolytic myokines (IL6, irisin, and LIF) stimulates thermogenesis by promoting the browning of adipocytes. Myokines, osteokines, and adipokines exert autocrine/paracrine effects locally as well as through the endocrine system, to regulate muscle, bone, and fat metabolism. Reductions in physical activity and increases in energy intake, both linked with aging, leads to adipocyte hypertrophy and the recruitment of immunological cells (macrophages). In turn, this releases pro-inflammatory adipokines which induces chronic low-grade inflammation (LGI), a key player in the pathology of several diseases. However, exercise-induced stimulation of bioactive cytokines, through muscle-bone-fat crosstalk, increases muscle anabolism, bone formation, mitochondrial biogenesis, glucose utilization, and fatty acid oxidation, and attenuates chronic LGI.

Postmenopausal osteoporosis is a musculoskeletal disease with a common genetic trait which responds to strength training: a translational intervention study

Background:

Clinical evidence suggests that body muscle mass is positively associated with bone mass, of significance for the elderly population at risk of osteoporosis (OP). Furthermore, muscle and bone interact mechanically and functionally, via local interactions as well as remotely via secreted components. Thus, it was of interest to compare muscle transcriptomes in postmenopausal OP and healthy women, and study effects of strength training on the muscle transcriptome, muscle stress proteins and bone mineral density (BMD).

Methods:

Skeletal muscle histological and genetic properties were compared in postmenopausal healthy (n = 18) and OP (n = 17) women before and after heavy-load strength training for 13–15 weeks. The cohorts were of similar age and body mass index without interfering diseases, medication or difference in lifestyle factors. Muscle biopsies obtained before and after intervention were studied histologically, and stress proteins and transcriptomes analyzed.

Results:

The OP women showed distinct muscle transcription profiles when compared with healthy women and had higher levels of the stress proteins HSP70 and α-β-crystalline. A set of 12 muscle transcripts, including ACSS3, FZD4, GNAI1 and IGF1, were differentially expressed before and after intervention (false discovery rate ⩽0.10, p ⩽0.001), and their corresponding bone transcripts were associated with BMD. Experimental data underline and describe the functionality of these genes in bone biology. OP women had 8% (p <0.01) higher proportion of type I fibres, but muscle fibre cross-sectional area did not differ. Muscle strength increased in both groups (p <0.01).

Conclusions:

Postmenopausal healthy and OP women have distinct muscle transcriptomes [messenger ribonucleic acids (mRNAs) and microRNAs] that are modulated by strength training, translating into key protein alterations and muscle fibre changes. The function of common skeletal muscle and bone genes in postmenopausal OP is suggestive of a shared disease trait.

Identification of PIEZO1 polymorphisms for human bone mineral density

Bone mineral density (BMD) is a key indicator for diagnosis and treatment for osteoporosis; the reduction of BMD could increase the risk of osteoporotic fracture. It was very recently found that Piezo1 mediated mechanically evoked responses in bone and further participated in bone formation in mice. Here, we performed cross phenotype meta-analysis for human BMD at lumbar spine (LS), femoral neck (FN), distal radius/forearm (FA) and heel and screened out 14 top SNPs for PIEZO1, these SNPs were overlapped with putative enhancers, DNase-I hypersensitive sites and active promoter flanking regions. We found that the signal of the best SNP rs62048221 was mainly from heel ultrasound estimated BMD (-0.02 SD per T allele, P = 8.50E-09), where calcaneus supported most of the mechanical force of body when standing, walking and doing physical exercises. Each copy of the effect allele T of SNP rs62048221 was associated with a decrease of 0.0035 g/cm² BMD (P = 4.6E-27, SE = 0.0003) in UK Biobank data within 477,760 samples. SNP rs62048221 was located at the enhancer region (HEDD enhancer ID 2331049) of gene PIEZO1, site-directed ChIP assays in human mesenchymal stem cells (hMSCs) showed significant enrichment of H3K4me1 and H3K27ac in this region, luciferase assays showed that rs62048221 could significantly affect the activity of the enhancer where it resides. Our results first suggested that SNP rs62048221 might mediate the PIEZO1 expression level via modulating the activity of cis-regulatory elements and then further affect the BMD.

Methods for multiple outcome meta-analysis of gene-expression data

Meta-analysis is a valuable tool for the synthesis of evidence across a wide range study types including high-throughput experiments such as genome-wide association studies (GWAS) and gene expression studies. There are situations though, in which we have multiple outcomes or multiple treatments, in which the multivariate meta-analysis framework which performs a joint modeling of the different quantities of interest may offer important advantages, such as increasing statistical power and allowing performing global tests. In this work we adapted the multivariate meta-analysis method and applied it in gene expression data. With this method we can test for pleiotropic effects, that is, for genes that influence both outcomes or discover genes that have a change in expression not detectable in the univariate method. We tested this method on data regarding inflammatory bowel disease (IBD), with its two main forms, Crohn’s disease (CD) and Ulcerative colitis (UC), sharing many clinical manifestations, but differing in the location and extent of inflammation and in complications. The Stata code is given in the Appendix and it is available at: www.compgen.org/tools/multivariate-microarrays.•

Multivariate meta-analysis method for gene expression data.

•

Discover genes with pleiotropic effects.

•

Differentially Expressed Genes (DEGs) identification in complex traits.

Replication of FTO Gene associated with lean mass in a Meta-Analysis of Genome-Wide Association Studies

Sarcopenia is characterized by low skeletal muscle, a complex trait with high heritability. With the dramatically increasing prevalence of obesity, obesity and sarcopenia occur simultaneously, a condition known as sarcopenic obesity. Fat mass and obesity-associated (FTO) gene is a candidate gene of obesity. To identify associations between lean mass and FTO gene, we performed a genome-wide association study (GWAS) of lean mass index (LMI) in 2207 unrelated Caucasian subjects and replicated major findings in two replication samples including 6,004 unrelated Caucasian and 38,292 unrelated Caucasian. We found 29 single nucleotide polymorphisms (SNPs) in FTO significantly associated with sarcopenia (combined p-values ranging from 5.92 × 10⁻¹² to 1.69 × 10⁻⁹). Potential biological functions of SNPs were analyzed by HaploReg v4.1, RegulomeDB, GTEx, IMPC and STRING. Our results provide suggestive evidence that FTO gene is associated with lean mass.

Association of 3p27.1 Variants with Whole Body Lean Mass Identified by a Genome-wide Association Study

Whole body lean mass (WBLM) is a heritable trait predicting sarcopenia. To identify genomic locus underlying WBLM, we performed a genome-wide association study of fat-adjusted WBLM in the Framingham Heart Study (FHS, N = 6,004), and replicated in the Kansas City Osteoporosis Study (KCOS, N = 2,207). We identified a novel locus 3p27.1 that was associated with WBLM (lead SNP rs3732593 P = 7.19 × 10^-8) in the discovery FHS sample, and the lead SNP was successfully replicated in the KCOS sample (one-sided P = 0.04). Bioinformatics analysis found that this SNP and its adjacent SNPs had the function of regulating enhancer activity in skeletal muscle myoblasts cells, further confirming the regulation of WBLM by this locus. Our finding provides new insight into the genetics of WBLM and enhance our understanding of sarcopenia.

A road map for understanding molecular and genetic determinants of osteoporosis

Osteoporosis is a highly prevalent disorder characterized by low bone mineral density and an increased risk of fracture, termed osteoporotic fracture. Notably, bone mineral density, osteoporosis and osteoporotic fracture are highly heritable; however, determining the genetic architecture, and especially the underlying genomic and molecular mechanisms, of osteoporosis in vivo in humans is still challenging. In addition to susceptibility loci identified in genome-wide association studies, advances in various omics technologies, including genomics, transcriptomics, epigenomics, proteomics and metabolomics, have all been applied to dissect the pathogenesis of osteoporosis. However, each technology individually cannot capture the entire view of the disease pathology and thus fails to comprehensively identify the underlying pathological molecular mechanisms, especially the regulatory and signalling mechanisms. A change to the status quo calls for integrative multi-omics and inter-omics analyses with approaches in 'systems genetics and genomics'. In this Review, we highlight findings from genome-wide association studies and studies using various omics technologies individually to identify mechanisms of osteoporosis. Furthermore, we summarize current studies of data integration to understand, diagnose and inform the treatment of osteoporosis. The integration of multiple technologies will provide a road map to illuminate the complex pathogenesis of osteoporosis, especially from molecular functional aspects, in vivo in humans.

AGS and NIA Bench-to Bedside Conference Summary: Osteoporosis and Soft Tissue (Muscle and Fat) Disorders

This report summarizes the presentations and recommendations of the eleventh annual American Geriatrics Society and National Institute on Aging research conference, "Osteoporosis and Soft Tissue (Muscle/Fat) Disorders," on March 11-12, 2019, in Bethesda, Maryland. Falls, fractures, and sarcopenia have a major impact on health in older adults, and they are interconnected by known risk factors. The link between osteoporosis, which is common in older adults, and the risk of falls is well known. Sarcopenia, the age-related decline in skeletal muscle mass and function, is also associated with an increased risk of falls and fractures because it reduces strength and leads to functional limitations. In addition to increasing the risk of falls, sarcopenia and osteoporosis can lead to frailty, reduced quality of life, morbidity, and mortality. The conference highlighted the impact of bone and soft tissue disorders on quality of life, morbidity, and mortality in older adults. Presenters described factors that contribute to these disorders; health disparities experienced by various subpopulations; and promising biological, pharmacologic, and behavioral interventions to prevent or treat these disorders. The workshop identified many research gaps and questions along with research recommendations that have the potential to enhance the prospect of healthy aging and improved quality of life for older adults. J Am Geriatr Soc 68:31-38, 2019.

Genomics of body fat percentage may contribute to sex bias in anorexia nervosa

Anorexia nervosa (AN) occurs nine times more often in females than in males. Although environmental factors likely play a role, the reasons for this imbalanced sex ratio remain unresolved. AN displays high genetic correlations with anthropometric and metabolic traits. Given sex differences in body composition, we investigated the possible metabolic underpinnings of female propensity for AN. We conducted sex-specific GWAS in a healthy and medication-free subsample of the UK Biobank (n = 155,961), identifying 77 genome-wide significant loci associated with body fat percentage (BF%) and 174 with fat-free mass (FFM). Partitioned heritability analysis showed an enrichment for central nervous tissue-associated genes for BF%, which was more prominent in females than males. Genetic correlations of BF% and FFM with the largest GWAS of AN by the Psychiatric Genomics Consortium were estimated to explore shared genomics. The genetic correlations of BF%male and BF%female with AN differed significantly from each other (p < .0001, δ = -0.17), suggesting that the female preponderance in AN may, in part, be explained by sex-specific anthropometric and metabolic genetic factors increasing liability to AN.

Osteosarcopenia: A case of geroscience

Many older persons lose their mobility and independence due to multiple diseases occurring simultaneously. Geroscience is aimed at developing innovative approaches to better identify relationships among the biological processes of aging. Osteoporosis and sarcopenia are two of the most prevalent chronic diseases in older people, with both conditions sharing overlapping risk factors and pathogenesis. When occurring together, these diseases form a geriatric syndrome termed "osteosarcopenia," which increases the risk of frailty, hospitalizations, and death. Findings from basic and clinical sciences aiming to understand osteosarcopenia have provided evidence of this syndrome as a case of geroscience. Genetic, endocrine, and mechanical stimuli, in addition to fat infiltration, sedentarism, and nutritional deficiencies, affect muscle and bone homeostasis to characterize this syndrome. However, research is in its infancy regarding accurate diagnostic markers and effective treatments with dual effects on muscle and bone. To date, resistance exercise remains the most promising strategy to increase muscle and bone mass, while sufficient quantities of protein, vitamin D, calcium, and creatine may preserve these tissues with aging. More recent findings, from rodent models, suggest treating ectopic fat in muscle and bone marrow as a possible avenue to curb osteosarcopenia, although this needs testing in human clinical trials.

Genetics of Bone and Muscle Interactions in Humans

PURPOSE OF REVIEW

To summarize the evidence from recent studies on the shared genetics between bone and muscle in humans.

RECENT FINDINGS

Genome-wide association studies (GWAS) have successfully identified a multitude of loci influencing the variability of different bone or muscle parameters, with multiple loci overlapping between the traits. In addition, joint analyses of multiple correlated musculoskeletal traits (i.e., multivariate GWAS) have underscored several genes with possible pleiotropic effects on both bone and muscle including MEF2C and SREBF1. Notably, several of the proposed pleiotropic genes have been validated using human cells or animal models. It is clear that the study of pleiotropy may provide novel insights into disease pathophysiology potentially leading to the identification of new treatment strategies that simultaneously prevent or treat both osteoporosis and sarcopenia. However, the role of muscle factors (myokines) that stimulate bone metabolism, as well as osteokines that affect muscles, is in its earliest stage of understanding.

Genome-scale Capture C promoter interactions implicate effector genes at GWAS loci for bone mineral density

Osteoporosis is a devastating disease with an essential genetic component. GWAS have discovered genetic signals robustly associated with bone mineral density (BMD), but not the precise localization of effector genes. Here, we carry out physical and direct variant to gene mapping in human mesenchymal progenitor cell-derived osteoblasts employing a massively parallel, high resolution Capture C based method in order to simultaneously characterize the genome-wide interactions of all human promoters. By intersecting our Capture C and ATAC-seq data, we observe consistent contacts between candidate causal variants and putative target gene promoters in open chromatin for ~ 17% of the 273 BMD loci investigated. Knockdown of two novel implicated genes, ING3 at ‘CPED1-WNT16’ and EPDR1 at ‘STARD3NL’, inhibits osteoblastogenesis, while promoting adipogenesis. This approach therefore aids target discovery in osteoporosis, here on the example of two relevant genes involved in the fate determination of mesenchymal progenitors, and can be applied to other common genetic diseases.

GWAS have identified numerous genetic loci for bone mineral density (BMD) and fracture risk. Here, the authors map these variants to putative target genes using ATAC-seq and Capture C of human osteoblasts and confirm ING3 and EPDR1 as BMD genes in in vitro osteoblast differentiation experiments.

Use it or lose it to age: A review of bone and muscle communication

Until recently, it was assumed that the only interaction between muscle and bone is mechanical, that the muscle acts as a pulley and the bone as a lever to move the organism. A relatively new concept is that muscle, especially contracted muscle, acts as a secretory organ, regulating metabolism. An even newer concept is that bone, especially the osteocytes in bone, act as endocrine cells targeting other organs such as kidney and more recently, muscle. These two new concepts logically led to the third concept: that muscle and bone communicate via soluble factors. Crosstalk occurs through muscle factors such as myostatin, irisin, and a muscle metabolite, β-aminoisobutyric acid, BAIBA, and through bone factors such as osteocalcin, transforming growth factor beta, TGFβ, Prostaglandin E2, PGE₂ and Wnts. Some of these factors have positive and some negative effects on the opposing tissue. One feature both bone and muscle have in common is that their tissues are mechanically loaded and many of their secreted factors are regulated by load. This mechanical loading, also known as exercise, has beneficial effects on many systems leading to the hypothesis that muscle and bone factors can be responsible for the beneficial effects of exercise. Many of the characteristics of aging and diseases associated with aging such as sarcopenia and osteoporosis and neurological conditions such as Alzheimer's disease and dementia, are delayed by exercise. This beneficial effect has been ascribed to increased blood flow increasing oxygen and nutrients, but could also be due to the secretome of the musculoskeletal system as outlined in this review.

Analysis of RAMP3 gene polymorphism with body composition and bone density in young and elderly women

Background and aim

The Receptor Activity Modifying Proteins (RAMPs) are a group of accessory proteins, of which there are three in humans, that interact with a number of G-protein coupled receptors (GPCR) and play various roles in regulation of endocrine signaling. Studies in RAMP3 knockout (KO) mice reveal an age related phenotype with altered metabolic regulation and high bone mass. To translate these findings into a clinically relevant perspective, we investigated the association between RAMP3 gene variants, body composition and bone phenotypes in two population-based cohorts of Swedish women.

Methods

Five single nucleotide polymorphisms (SNP) in the vicinity of the RAMP3 gene were genotyped in the PEAK-25 cohort (n = 1061; 25 years) and OPRA (n = 1044; 75 years). Bone mineral density (BMD), fat mass and lean mass (total body; regional) were measured by DXA at baseline, 5 and 10 year follow-up.

Results

BMD did not differ with RAMP3 genotype in either cohort, although fracture risk was increased in the elderly women (OR 2.695 [95% CI 1.514–4.801]). Fat mass tended to be higher with RAMP3 SNPs; although only in elderly women. In the young women, changes in BMI and fat mass between ages 25–35 differed by genotype (p = 0.001; p < 0.001).

Conclusion

Variation in RAMP3 may contribute to age-related changes in body composition and risk of fracture.