Prioritization of Genes Relevant to Bone Fragility Through the Unbiased Integration of Aging Mouse Bone Transcriptomics and Human GWAS Analyses

TMEM175 plays a crucial role in osteoblast differentiation by regulating lysosomal function and autophagy

Bone provides structural support, enables movement, protects internal organs, regulates calcium and phosphorus levels, and contains bone marrow essential for hematopoiesis. Osteoblasts are specialized cells responsible for bone formation through the secretion of extracellular matrix components. Transmembrane protein 175 (TMEM175), which functions as an endosomal/lysosomal K+ channel and a lysosomal H+ channel, regulates lysosomal function and autophagy. Despite the recognized importance of lysosomes and autophagy in osteoblast differentiation, the specific role of TMEM175 in osteoblast differentiation has not been revealed. In this study, we investigated whether TMEM175 is associated with human bone mineral density and fracture and examined the role of TMEM175 in osteoblast differentiation. In analyses of single nucleotide polymorphisms of pore ion channel genes using the mouse2human database, a significant correlation between TMEM175 single nucleotide polymorphisms and human bone mineral density and fracture was identified. TMEM175 expression levels were found to increase during osteoblast differentiation from bone chip-derived mesenchymal stem cells (BMSCs). Knockdown of TMEM175 in BMSCs suppressed osteoblast differentiation, as evidenced by decreased matrix mineralization and lower expression levels of osteoblast marker genes. Further analysis indicated that TMEM175 deficiency leads to lysosomal dysfunction and partially impairs autophagic clearance during osteoblast differentiation. Moreover, the TMEM175 inhibitor 4-aminopyridine decreased osteoblast differentiation of BMSCs. Taken together, this study reveals that TMEM175 plays an important role in osteoblast differentiation by regulating lysosomal function and autophagic clearance.

An integrated multi-omics analysis reveals osteokines involved in global regulation

Bone secretory proteins, termed osteokines, regulate bone metabolism and whole-body homeostasis. However, fundamental questions as to what the bona fide osteokines and their cellular sources are and how they are regulated remain unclear. In this study, we analyzed bone and extraskeletal tissues, osteoblast (OB) conditioned media, bone marrow supernatant (BMS), and serum, for basal osteokines and those responsive to aging and mechanical loading/unloading. We identified 375 candidate osteokines and their changes in response to aging and mechanical dynamics by integrating data from RNA-seq, scRNA-seq, and proteomic approaches. Furthermore, we analyzed their cellular sources in the bone and inter-organ communication facilitated by them (bone-brain, liver, and aorta). Notably, we discovered that senescent OBs secrete fatty-acid-binding protein 3 to propagate senescence toward vascular smooth muscle cells (VSMCs). Taken together, we identified previously unknown candidate osteokines and established a dynamic regulatory network among them, thus providing valuable resources to further investigate their systemic roles.

Aging impairs the osteocytic regulation of collagen integrity and bone quality

Poor bone quality is a major factor in skeletal fragility in elderly individuals. The molecular mechanisms that establish and maintain bone quality, independent of bone mass, are unknown but are thought to be primarily determined by osteocytes. We hypothesize that the age-related decline in bone quality results from the suppression of osteocyte perilacunar/canalicular remodeling (PLR), which maintains bone material properties. We examined bones from young and aged mice with osteocyte-intrinsic repression of TGFβ signaling (TβRIIocy-/-) that suppresses PLR. The control aged bone displayed decreased TGFβ signaling and PLR, but aging did not worsen the existing PLR suppression in male TβRIIocy-/- bone. This relationship impacted the behavior of collagen material at the nanoscale and tissue scale in macromechanical tests. The effects of age on bone mass, density, and mineral material behavior were independent of osteocytic TGFβ. We determined that the decline in bone quality with age arises from the loss of osteocyte function and the loss of TGFβ-dependent maintenance of collagen integrity.

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.

Whole transcriptome analysis in offspring whose fathers were exposed to a developmental insult: a novel avian model

Although the effects of paternal exposure to insults on the offspring received limited attention in the past, it is currently gaining interest especially after understanding the mechanisms which may mediate such exposure effects. In the current study, the well-controlled avian model (Fayoumi) was utilized to investigate the effects of paternal exposure to the developmental insult, chlorpyrifos on the offspring's gene expression via mRNA and small RNA sequencing. Numerous mRNA gene expression changes were detected in the offspring after paternal exposure to the developmental insult, especially in genes related to neurogenesis, learning and memory. qPCR analysis of several genes, that were significantly changed in mRNA sequencing, confirmed the results obtained in mRNA sequencing. On the other hand, small RNA sequencing did not identify significant microRNA genes expression changes in the offspring after paternal exposure to the developmental insult. The effects of the paternal exposure were more pronounced in the female offspring compared to the male offspring. The results identified expression alterations in major genes (some of which were pertinent to the functional changes observed in other forms of early developmental exposure) after paternal insult exposure and provided a direction for future studies involving the most affected genes.

Single-Cell Integration of BMD GWAS Results Prioritize Candidate Genes Influencing Age-Related Bone Loss

The regulation of bone mineral density (BMD) is highly influenced by genetics and age. Although genome-wide association studies (GWAS) for BMD have uncovered many genes through their proximity to associated variants (variant nearest-neighbor [VNN] genes), the cell-specific mechanisms of each VNN gene remain unclear. This is primarily due to the inability to prioritize these genes by cell type and age-related expression. Using age-related transcriptomics, we found that the expression of many VNN genes was upregulated in the bone and marrow from aged mice. Candidate genes from GWAS were investigated using single-cell RNA-sequencing (scRNA-seq) datasets to enrich for cell-specific expression signatures. VNN candidate genes are highly enriched in osteo-lineage cells, osteocytes, hypertrophic chondrocytes, and Lepr+ mesenchymal stem cells. These data were used to generate a "blueprint" for Cre-loxp mouse line selection for functional validation of candidate genes and further investigation of their role in BMD maintenance throughout aging. In VNN-gene-enriched cells, Sparc, encoding the extracellular matrix (ECM) protein osteonectin, was robustly expressed. This, along with expression of numerous other ECM genes, indicates that many VNN genes likely have roles in ECM deposition by osteoblasts. Overall, we provide data supporting streamlined translation of GWAS candidate genes to potential novel therapeutic targets for the treatment of osteoporosis. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Examining craniofacial variation among crispant and mutant zebrafish models of human skeletal diseases

Genetic diseases affecting the skeletal system present with a wide range of symptoms that make diagnosis and treatment difficult. Genome-wide association and sequencing studies have identified genes linked to human skeletal diseases. Gene editing of zebrafish models allows researchers to further examine the link between genotype and phenotype, with the long-term goal of improving diagnosis and treatment. While current automated tools enable rapid and in-depth phenotyping of the axial skeleton, characterizing the effects of mutations on the craniofacial skeleton has been more challenging. The objective of this study was to evaluate a semi-automated screening tool can be used to quantify craniofacial variations in zebrafish models using four genes that have been associated with human skeletal diseases (meox1, plod2, sost, and wnt16) as test cases. We used traditional landmarks to ground truth our dataset and pseudolandmarks to quantify variation across the 3D cranial skeleton between the groups (somatic crispant, germline mutant, and control fish). The proposed pipeline identified variation between the crispant or mutant fish and control fish for four genes. Variation in phenotypes parallel human craniofacial symptoms for two of the four genes tested. This study demonstrates the potential as well as the limitations of our pipeline as a screening tool to examine multi-dimensional phenotypes associated with the zebrafish craniofacial skeleton.

Bone-cartilage crosstalk informed by aging mouse bone transcriptomics and human osteoarthritis genome-wide association studies

Subchondral bone participates in crosstalk with articular cartilage to maintain joint homeostasis, and disruption of either tissue results in overall joint degeneration. Among the subchondral bone changes observed in osteoarthritis (OA), subchondral bone plate (SBP) thickening has a time-dependent relationship with cartilage degeneration and has recently been shown to be regulated by osteocytes. Here, we evaluate the effect of age on SBP thickness and cartilage degeneration in aging mice. We find that SBP thickness significantly increases by 18-months of age, corresponding temporally with increased cartilage degeneration. To identify factors in subchondral bone that may participate in bone cartilage crosstalk or OA, we leveraged mouse transcriptomic data from one joint tissue compartment - osteocyte-enriched bone - to search for enrichment with human OA in UK Biobank and Arthritis Research UK Osteoarthritis Genetics (arcOGEN) GWAS using the mouse2human (M2H, www.mouse2human.org) strategy. Genes differentially expressed in aging mouse bone are significantly enriched for human OA, showing joint site-specific (knee vs. hip) relationships, exhibit temporal associations with age, and unique gene clusters are implicated in each type of OA. Application of M2H identifies genes with known and unknown functions in osteocytes and OA development that are clinically associated with human OA. Altogether, this work prioritizes genes with a potential role in bone/cartilage crosstalk for further mechanistic study based on their association with human OA in GWAS.

Differentially Expressed Genes and Molecular Susceptibility to Human Age-Related Diseases

Mainstream transcriptome profiling of susceptibility versus resistance to age-related diseases (ARDs) is focused on differentially expressed genes (DEGs) specific to gender, age, and pathogeneses. This approach fits in well with predictive, preventive, personalized, participatory medicine and helps understand how, why, when, and what ARDs one can develop depending on their genetic background. Within this mainstream paradigm, we wanted to find out whether the known ARD-linked DEGs available in PubMed can reveal a molecular marker that will serve the purpose in anyone's any tissue at any time. We sequenced the periaqueductal gray (PAG) transcriptome of tame versus aggressive rats, identified rat-behavior-related DEGs, and compared them with their known homologous animal ARD-linked DEGs. This analysis yielded statistically significant correlations between behavior-related and ARD-susceptibility-related fold changes (log₂ values) in the expression of these DEG homologs. We found principal components, PC1 and PC2, corresponding to the half-sum and the half-difference of these log₂ values, respectively. With the DEGs linked to ARD susceptibility and ARD resistance in humans used as controls, we verified these principal components. This yielded only one statistically significant common molecular marker for ARDs: an excess of Fcγ receptor IIb suppressing immune cell hyperactivation.

High Bone Mass Disorders: New Insights From Connecting the Clinic and the Bench

Monogenic high bone mass (HBM) disorders are characterized by an increased amount of bone in general, or at specific sites in the skeleton. Here, we describe 59 HBM disorders with 50 known disease-causing genes from the literature, and we provide an overview of the signaling pathways and mechanisms involved in the pathogenesis of these disorders. Based on this, we classify the known HBM genes into HBM (sub)groups according to uniform Gene Ontology (GO) terminology. This classification system may aid in hypothesis generation, for both wet lab experimental design and clinical genetic screening strategies. We discuss how functional genomics can shape discovery of novel HBM genes and/or mechanisms in the future, through implementation of omics assessments in existing and future model systems. Finally, we address strategies to improve gene identification in unsolved HBM cases and highlight the importance for cross-laboratory collaborations encompassing multidisciplinary efforts to transfer knowledge generated at the bench to the clinic. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Pathways Controlling Formation and Maintenance of the Osteocyte Dendrite Network

PURPOSE OF REVIEW

The purpose of this review is to discuss the molecular mechanisms involved in osteocyte dendrite formation, summarize the similarities between osteocytic and neuronal projections, and highlight the importance of osteocyte dendrite maintenance in human skeletal disease.

RECENT FINDINGS

It is suggested that there is a causal relationship between the loss of osteocyte dendrites and the increased osteocyte apoptosis during conditions including aging, microdamage, and skeletal disease. A few mechanisms are proposed to control dendrite formation and outgrowth, such as via the regulation of actin polymerization dynamics. This review addresses the impact of osteocyte dendrites in bone health and disease. Recent advances in multi-omics, in vivo and in vitro models, and microscopy-based imaging have provided novel approaches to reveal the underlying mechanisms that regulate dendrite development. Future therapeutic approaches are needed to target the process of osteocyte dendrite formation.