RNA-seq in Skeletal Biology

Sodium Alginate Hydrogel Infusion of Bone Marrow Mesenchymal Stem Cell-Derived Extracellular Vesicles and p38α Antagonistic Peptides in Myocardial Infarction Fibrosis Mitigation

BACKGROUND

Myocardial fibrosis is a pathological hallmark of heart failure post infarction, emphasizing the need for innovative treatment strategies. This research assesses the antifibrotic potential of a sodium alginate (SA) hydrogel loaded with extracellular vesicles (EVs) from bone marrow mesenchymal stem cells and PAP (p38α antagonistic peptides), aiming to interfere with fibrosis-inducing pathways in myocardial tissue after infarction.

METHODS

We induced fibrosis in mouse cardiac fibroblasts through hypoxia and disrupted the Mapk14 gene to study its contribution to fibrosis. Mesenchymal stem cell-derived EVs, loaded with PAP, were encapsulated in the SA hydrogel (EVs-PAP@SA). The formulation was tested in vitro for its effect on fibrotic marker expression and cell behavior, and in vivo in a murine model of myocardial infarction for its therapeutic efficacy.

RESULTS

Map k14 silencing showed a decrease in the fibrotic response of cardiac fibroblasts. Treatment with the EVs-PAP@SA hydrogel notably reduced profibrotic signaling, increased cell proliferation and migration, and lowered apoptosis rates. The in vivo treatment with the hydrogel post myocardial infarction significantly diminished myocardial fibrosis and improved cardiac performance.

CONCLUSIONS

The study endorses the SA hydrogel as an effective vehicle for delivering mesenchymal stem cell-derived EVs and PAP to the heart post myocardial infarction, providing a novel approach for modulating myocardial fibrosis and promoting cardiac healing.

Ethyl butyrate inhibits caudal fin regeneration in adult zebrafish by disrupting extracellular matrix remodeling

Wound healing and tissue regeneration are influenced by a variety of factors. Adverse lifestyle habits, such as excessive alcohol consumption, delay wound healing and increase the risk of secondary infections. Ethyl butyrate is a common food additive widely used to enhance the aroma of alcoholic beverages. This additive is generally considered harmless to human health in both industrial and domestic settings. However, the ecotoxicity and its effects on wound healing have not been elucidated. In this study, we used zebrafish as the experimental animal, and the caudal fins were amputated to explore the effects of ethyl butyrate on wound healing and tissue regeneration. The effect of ethyl butyrate on blastema and bone regeneration and its impact on the transcriptional levels of regeneration-related genes and inflammation-related genes were evaluated. RNA-seq was conducted to determine the differentially expressed genes (DEGs) between the treatment and the control groups. KEGG and GO analysis was conducted to explore the functions of DEGs. Significantly enriched GO terms and KEGG pathways were identified to explore the molecular mechanism underlying the inhibition of zebrafish caudal fin regeneration by ethyl butyrate. The results demonstrated that ethyl butyrate significantly inhibited the regeneration of zebrafish caudal fins, including blastema and bone regeneration. Ethyl butyrate exposure significantly downregulated the expression of genes associated with bone and blastema regeneration and inflammation response. KEGG and GO functional analyses revealed that the DEGs were associated with significant enrichment of extracellular matrix-receptor interactions. Ethyl butyrate treatment downregulated the expression of most extracellular matrix-related genes. These findings indicate that ethyl butyrate potentially modulates pathways associated with the structure, adhesion, modification, and degradation of the extracellular matrix, thereby disrupting extracellular matrix remodeling, inhibiting wound inflammation, impairing blastema and bone regeneration and ultimately hindering caudal fin regeneration. In summary, the findings demonstrate that ethyl butyrate disrupts extracellular matrix remodeling and inhibits the regeneration of zebrafish caudal fins. These results provide valuable insights into the rational use of ethyl butyrate and further investigation of wound healing mechanisms.

A New Tissue Engineering Strategy to Promote Tendon-bone Healing: Regulation of Osteogenic and Chondrogenic Differentiation of Tendon-derived Stem Cells

In the field of sports medicine, repair surgery for anterior cruciate ligament (ACL) and rotator cuff (RC) injuries are remarkably common. Despite the availability of relatively effective treatment modalities, outcomes often fall short of expectations. This comprehensive review aims to thoroughly examine current strategies employed to promote tendon-bone healing and analyze pertinent preclinical and clinical research. Amidst ongoing investigations, tendon-derived stem cells (TDSCs), which have comparatively limited prior exploration, have garnered increasing attention in the context of tendon-bone healing, emerging as a promising cell type for regenerative therapies. This review article delves into the potential of combining TDSCs with tissue engineering methods, with ACL reconstruction as the main focus. It comprehensively reviews relevant research on ACL and RC healing to address the issues of graft healing and bone tunnel integration. To optimize tendon-bone healing outcomes, our emphasis lies in not only reconstructing the original microstructure of the tendon-bone interface but also achieving proper bone tunnel integration, encompassing both cartilage and bone formation. In this endeavor, we thoroughly analyze the transcriptional and molecular regulatory variables governing TDSCs differentiation, incorporating a retrospective analysis utilizing single-cell sequencing, with the aim of unearthing relevant signaling pathways and processes. By presenting a novel strategy rooted in TDSCs-driven osteogenic and chondrogenic differentiation for tendon-bone healing, this study paves the way for potential future research avenues and promising therapeutic applications. It is anticipated that the findings herein will contribute to advancing the field of tendon-bone healing and foster the exploration of TDSCs as a viable option for regenerative therapies in the future.

SIRT1-mediated deacetylation of FOXO3 enhances mitophagy and drives hormone resistance in endometrial cancer

BACKGROUND

The complex interplay between Sirtuin 1 (SIRT1) and FOXO3 in endometrial cancer (EC) remains understudied. This research aims to unravel the interactions of deacetylase SIRT1 and transcription factor FOXO3 in EC, focusing on their impact on mitophagy and hormone resistance.

METHODS

High-throughput sequencing, cell experiments, and bioinformatics tools were employed to investigate the roles and interactions of SIRT1 and FOXO3 in EC. Co-immunoprecipitation (Co-IP) assay was used to assess the interaction between SIRT1 and FOXO3 in RL95-2 cells. Functional assays were used to assess cell viability, proliferation, migration, invasion, apoptosis, and the expression of related genes and proteins. A mouse model of EC was established to evaluate tumor growth and hormone resistance under different interventions. Immunohistochemistry and TUNEL assays were used to assess protein expression and apoptosis in tumor tissues.

RESULTS

High-throughput transcriptome sequencing revealed a close association between SIRT1, FOXO3, and EC development. Co-IP showed a protein-protein interaction between SIRT1 and FOXO3. Overexpression of SIRT1 enhanced FOXO3 deacetylation and activity, promoting BNIP3 transcription and PINK1/Parkin-mediated mitophagy, which in turn promoted cell proliferation, migration, invasion, and inhibited apoptosis in vitro, as well as increased tumor growth and hormone resistance in vivo. These findings highlighted SIRT1 as an upstream regulator and potential therapeutic target in EC.

CONCLUSION

This study reveals a novel molecular mechanism underlying the functional relevance of SIRT1 in regulating mitophagy and hormone resistance through the deacetylation of FOXO3 in EC, thereby providing valuable insights for new therapeutic strategies.

A novel spherical GelMA-HAMA hydrogel encapsulating APET×2 polypeptide and CFIm25-targeting sgRNA for immune microenvironment modulation and nucleus pulposus regeneration in intervertebral discs

METHODS

Single-cell transcriptomics and high-throughput transcriptomics were used to screen factors significantly correlated with intervertebral disc degeneration (IDD). Expression changes of CFIm25 were determined via RT-qPCR and Western blot. NP cells were isolated from mouse intervertebral discs and induced to degrade with TNF-α and IL-1β. CFIm25 was knocked out using CRISPR-Cas9, and CFIm25 knockout and overexpressing nucleus pulposus (NP) cell lines were generated through lentiviral transfection. Proteoglycan expression, protein expression, inflammatory factor expression, cell viability, proliferation, migration, gene expression, and protein expression were analyzed using various assays (alcian blue staining, immunofluorescence, ELISA, CCK-8, EDU labeling, transwell migration, scratch assay, RT-qPCR, Western blot). The GelMA-HAMA hydrogel loaded with APET×2 polypeptide and sgRNA was designed, and its effects on NP regeneration were assessed through in vitro and mouse model experiments. The progression of IDD in mice was evaluated using X-ray, H&E staining, and Safranin O-Fast Green staining. Immunohistochemistry was performed to determine protein expression in NP tissue. Proteomic analysis combined with in vitro and in vivo experiments was conducted to elucidate the mechanisms of hydrogel action.

RESULTS

CFIm25 was upregulated in IDD NP tissue and significantly correlated with disease progression. Inhibition of CFIm25 improved NP cell degeneration, enhanced cell proliferation, and migration. The hydrogel effectively knocked down CFIm25 expression, improved NP cell degeneration, promoted cell proliferation and migration, and mitigated IDD progression in a mouse model. The hydrogel inhibited inflammatory factor expression (IL-6, iNOS, IL-1β, TNF-α) by targeting the p38/NF-κB signaling pathway, increased collagen COLII and proteoglycan Aggrecan expression, and suppressed NP degeneration-related factors (COX-2, MMP-3).

CONCLUSION

The study highlighted the crucial role of CFIm25 in IDD and introduced a promising therapeutic strategy using a porous spherical GelMA-HAMA hydrogel loaded with APET×2 polypeptide and sgRNA. This innovative approach offers new possibilities for treating degenerated intervertebral discs.

Screening of mRNA markers in early bovine tuberculosis blood samples

Bovine tuberculosis (bTB) is a chronic zoonotic disease caused by Mycobacterium bovis. A large number of cattle are infected with bTB every year, resulting in huge economic losses. How to control bTB is an important issue in the current global livestock economy. In this study, the original transcriptome sequences related to this study were obtained from the dataset GSE192537 by searching the Gene Expression Omnibus (GEO) database. Our differential gene analysis showed that there were obvious biological activities related to immune activation and immune regulation in the early stage of bTB. Immune-related biological processes were more active in the early stage of bTB than in the late. There were obvious immune activation and immune cell recruitment in the early stage of bTB. Regulations in immune receptors are associated with pathophysiological processes of the early stage of bTB. A gene module consisting of 236 genes significantly related to the early stage of bTB was obtained by weighted gene co-expression network analysis, and 18 hub genes were further identified as potential biomarkers or therapeutic targets. Finally, by random forest algorithm and logistic regression modeling, FCRL1 was identified as a representative mRNA marker in early bTB blood. FCRL1 has the potential to be a diagnostic biomarker in early bTB.

Understanding blaNDM-1 gene regulation in CRKP infections: toward novel antimicrobial strategies for hospital-acquired pneumonia

BACKGROUND

The escalating challenge of Carbapenem-resistant Klebsiella pneumoniae (CRKP) in hospital-acquired pneumonia (HAP) is closely linked to the blaNDM-1 gene. This study explores the regulatory mechanisms of blaNDM-1 expression and aims to enhance antibacterial tactics to counteract the spread and infection of resistant bacteria.

METHODS

KP and CRKP strains were isolated from HAP patients' blood samples. Transcriptomic sequencing (RNA-seq) identified significant upregulation of blaNDM-1 gene expression in CRKP strains. Bioinformatics analysis revealed blaNDM-1 gene involvement in beta-lactam resistance pathways. CRISPR-Cas9 was used to delete the blaNDM-1 gene, restoring sensitivity. In vitro and in vivo experiments demonstrated enhanced efficacy with Imipenem and Thanatin or Subatan combination therapy.

RESULTS

KP and CRKP strains were isolated with significant upregulation of blaNDM-1 in CRKP strains identified by RNA-seq. The Beta-lactam resistance pathway was implicated in bioinformatics analysis. Knockout of blaNDM-1 reinstated sensitivity in CRKP strains. Further, co-treatment with Imipenem, Thanatin, or Subactam markedly improved antimicrobial effectiveness.

CONCLUSION

Silencing blaNDM-1 in CRKP strains from HAP patients weakens their Carbapenem resistance and optimizes antibacterial strategies. These results provide new theoretical insights and practical methods for treating resistant bacterial infections.

Single-cell RNA sequencing in osteoarthritis

Osteoarthritis is a progressive and heterogeneous joint disease with complex pathogenesis. The various phenotypes associated with each patient suggest that better subgrouping of tissues associated with genotypes in different phases of osteoarthritis may provide new insights into the onset and progression of the disease. Recently, single-cell RNA sequencing was used to describe osteoarthritis pathogenesis on a high-resolution view surpassing traditional technologies. Herein, this review summarizes the microstructural changes in articular cartilage, meniscus, synovium and subchondral bone that are mainly due to crosstalk amongst chondrocytes, osteoblasts, fibroblasts and endothelial cells during osteoarthritis progression. Next, we focus on the promising targets discovered by single-cell RNA sequencing and its potential applications in target drugs and tissue engineering. Additionally, the limited amount of research on the evaluation of bone-related biomaterials is reviewed. Based on the pre-clinical findings, we elaborate on the potential clinical values of single-cell RNA sequencing for the therapeutic strategies of osteoarthritis. Finally, a perspective on the future development of patient-centred medicine for osteoarthritis therapy combining other single-cell multi-omics technologies is discussed. This review will provide new insights into osteoarthritis pathogenesis on a cellular level and the field of applications of single-cell RNA sequencing in personalized therapeutics for osteoarthritis in the future.

Immunotoxicity and transcriptome analysis of zebrafish embryos exposure to Nitazoxanide

Nitazoxanide (NTZ) is a broad-spectrum immunomodulatory drug, and little information is about the immunotoxicity of aquatic organisms induced by NTZ. In the present study, reduced body length and decreased yolk sac absorption in the NTZ-treated group were observed. Meanwhile, the number of innate immune cells and adaptive immune cells was substantially reduced upon NTZ exposure, and the migration and retention of macrophages and neutrophils in the injured area were inhibited. Following NTZ stimulation, oxidative stress levels in the zebrafish increased obviously. Mechanistically, RNA-seq, a high-throughput method, was performed to analyze the global expression of differentially expressed genes (DEGs) in zebrafish embryos treated with NTZ. 531 DEGs were identified by comparative transcriptome analysis, including 121 up-regulated and 420 down-regulated genes in zebrafish embryos after NTZ exposure. The transcriptome sequences were further subjected to the Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) and analysis, showing phototransduction and metabolic pathway, respectively, and were most enriched. In addition, some immune-related genes were inhibited after NTZ exposure. RNA-seq results confirmed by qRT-PCR were used to verify the expression of the 6 selected genes. The other immune-related genes such as two pro-inflammatory cytokines (IL-1β, tnfα) and two chemokines (CXCL8b.3, CXCL-c1c) were further confirmed and were differentially regulated after NTZ exposure. In summary, NTZ exposure could lead to immunotoxicity and increased ROS in zebrafish embryos, this study provides valuable information for future elucidating the molecular mechanism of exogenous stimuli-induced immunotoxicity in aquatic ecosystems.

Identification of Differentially Expressed Genes and Molecular Pathways Involved in Osteoclastogenesis Using RNA-seq

Osteoporosis is a disease that is characterised by reduced bone mineral density (BMD) and can be exacerbated by the excessive bone resorption of osteoclasts (OCs). Bioinformatic methods, including functional enrichment and network analysis, can provide information about the underlying molecular mechanisms that participate in the progression of osteoporosis. In this study, we harvested human OC-like cells differentiated in culture and their precursor peripheral blood mononuclear cells (PBMCs) and characterised the transcriptome of the two cell types using RNA-sequencing in order to identify differentially expressed genes. Differential gene expression analysis was performed in RStudio using the edgeR package. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to identify enriched GO terms and signalling pathways, with inter-connected regions characterised using protein-protein interaction analysis. In this study, we identified 3201 differentially expressed genes using a 5% false discovery rate; 1834 genes were upregulated, whereas 1367 genes were downregulated. We confirmed a significant upregulation of several well-established OC genes including CTSK, DCSTAMP, ACP5, MMP9, ITGB3, and ATP6V0D2. The GO analysis suggested that upregulated genes are involved in cell division, cell migration, and cell adhesion, while the KEGG pathway analysis highlighted oxidative phosphorylation, glycolysis and gluconeogenesis, lysosome, and focal adhesion pathways. This study provides new information about changes in gene expression and highlights key biological pathways involved in osteoclastogenesis.

Polydopamine-coated 3D-printed β-tricalcium phosphate scaffolds to promote the adhesion and osteogenesis of BMSCs for bone-defect repair: mRNA transcriptomic sequencing analysis

Cellular bioactivity and tissue regeneration can be affected by coatings on tissue-engineered scaffolds. Using mussel-inspired polydopamine (PDA) is a convenient and effective approach to surface modification. Therefore, 3D-printed β-tricalcium phosphate (β-TCP) scaffolds were coated with PDA in this study. The effects of the scaffolds on the adhesion and osteogenic differentiation of seeded bone marrow mesenchymal stem cells (BMSCs) in vitro and on new-bone formation in vivo were investigated. The potential mechanisms and related differential genes were assessed using mRNA sequencing. It was seen that PDA coating increased the surface roughness of the 3D-printed β-TCP scaffolds. Furthermore, it prompted the adhesion and osteogenic differentiation of seeded BMSCs. mRNA sequencing analysis revealed that PDA coating might affect the osteogenic differentiation of BMSCs through the calcium signaling pathway, Wnt signaling pathway, TGF-beta signaling pathway, etc. Moreover, the expression of osteogenesis-related genes, such as R-spondin 1 and chemokine c-c-motif ligand 2, was increased. Finally, both the 3D-printed β-TCP scaffolds and PDA-coated scaffolds could significantly accelerate the formation of new bone in critical-size calvarial defects in rats compared with the control group; and the new bone formation was obviously higher in the PDA-coated scaffolds than in β-TCP scaffolds. In summary, 3D-printed β-TCP scaffolds with a PDA coating can improve the physicochemical characteristics and cellular bioactivity of the scaffold surface for bone regeneration. Potential differential genes were identified, which can be used as a foundation for further research.

Integration of metabolomics and transcriptomics provides insights into enhanced osteogenesis in Ano5Cys360Tyr knock-in mouse model

INTRODUCTION

Gnathodiaphyseal dysplasia (GDD; OMIM#166260) is a rare autosomal dominant disorder characterized by diaphyseal sclerosis of tubular bones and cemento-osseous lesions in mandibles. GDD is caused by point mutations in the ANO5 gene. However, the mechanisms underlying GDD have not been disclosed. We previously generated the first knock-in mouse model for GDD expressing a human mutation (p.Cys360Tyr) in ANO5 and homozygous Ano5 knock-in (Ano5^KI/KI ) mice exhibited representative traits of human GDD especially including enhanced osteogenesis.

METHODS

Metabolomics and transcriptomics analyses were conducted for wildtype (Ano5^+/+ ) and Ano5^KI/KI mature mouse calvarial osteoblasts (mCOBs) grown in osteogenic cultures for 14 days to identify differential intracellular metabolites and genes involved in GDD. Subsequently, related differential genes were validated by qRT-PCR. Cell proliferation was confirmed by CCK8 assay and calcium content in mineral nodules was detected using SEM-EDS.

RESULTS

Metabolomics identified 42 differential metabolites that are primarily involved in amino acid and pyrimidine metabolism, and endocrine and other factor-regulated calcium reabsorption. Concomitantly, transcriptomic analysis revealed 407 differentially expressed genes in Ano5^KI/KI osteoblasts compared with wildtype. Gene ontology and pathway analysis indicated that Ano5^Cys360Tyr mutation considerably promoted cell cycle progression and perturbed calcium signaling pathway, which were confirmed by validated experiments. qRT-PCR and CCK-8 assays manifested that proliferation of Ano5^KI/KI mCOBs was enhanced and the expression of cell cycle regulating genes (Mki67, Ccnb1, and Ccna2) was increased. In addition, SEM-EDS demonstrated that Ano5^KI/KI mCOBs developed higher calcium contents in mineral nodules than Ano5^+/+ mCOBs, while some calcium-related genes (Cacna1, Slc8a1, and Cyp27b1) were significantly up-regulated. Furthermore, osteocalcin which has been proved to be an osteoblast-derived metabolic hormone was upregulated in Ano5^KI/KI osteoblast cultures.

DISCUSSION

Our data demonstrated that the Ano5^Cys360Tyr mutation could affect the metabolism of osteoblasts, leading to unwonted calcium homeostasis and cellular proliferation that can contribute to the underlying pathogenesis of GDD disorders.

An optimized FACS-free single-nucleus RNA sequencing (snRNA-seq) method for plant science research

Recently, single-cell RNA sequencing (scRNA-seq) provides unprecedented power for accurately understanding gene expression regulatory mechanisms. However, scRNA-seq studies have limitations in plants, due to difficulty in protoplast isolation that requires enzymatic digestion of the cell walls from various plant tissues. Therefore, to overcome this problem, we developed a nuclei isolation approach that does not rely on Fluorescence Activated Cell Sorting (FACS). We validated the robustness of the FACS-free single-nucleus RNA sequencing (snRNA-seq) methodology in mature Arabidopsis plant tissue by comparing it to scRNA-seq results based on protoplasts extracted from the same batch of leaf materials. Sequencing results demonstrated the high quality of snRNA-seq data, as well as its utility in cell type classification and marker gene identification. This approach also showed several advantages, including the ability to use frozen samples, taking less suspension preparation time, and reducing biased cellular coverage and dissociation-induced transcriptional artifacts. Surprisingly, snRNA-seq detected two epidermal pavement cell clusters, while scRNA-seq only had one. Furthermore, we hypothesized that these two epidermal cells represent the top and lower epidermis based on differences in expression patterns of cluster-specific expressed genes. In summary, this study has advanced the application of snRNA-seq in Arabidopsis leaves and confirmed the advantages of snRNA-seq in plant research.

3D bioprinting of in situ vascularized tissue engineered bone for repairing large segmental bone defects

Large bone defects remain an unsolved clinical challenge because of the lack of effective vascularization in newly formed bone tissue. 3D bioprinting is a fabrication technology with the potential to create vascularized bone grafts with biological activity for repairing bone defects. In this study, vascular endothelial cells laden with thermosensitive bio-ink were bioprinted in situ on the inner surfaces of interconnected tubular channels of bone mesenchymal stem cell-laden 3D-bioprinted scaffolds. Endothelial cells exhibited a more uniform distribution and greater seeding efficiency throughout the channels. In vitro, the in situ bioprinted endothelial cells can form a vascular network through proliferation and migration. The in situ vascularized tissue-engineered bone also resulted in a coupling effect between angiogenesis and osteogenesis. Moreover, RNA sequencing analysis revealed that the expression of genes related to osteogenesis and angiogenesis is upregulated in biological processes. The in vivo 3D-bioprinted in situ vascularized scaffolds exhibited excellent performance in promoting new bone formation in rat calvarial critical-sized defect models. Consequently, in situ vascularized tissue-engineered bones constructed using 3D bioprinting technology have a potential of being used as bone grafts for repairing large bone defects, with a possible clinical application in the future.

•

3D bioprinting was used to fabricate in situ vascularized tissue engineered bone.

•

In situ bioprinted endothelial cells exhibited uniform distribution and greater seeding efficiency.

•

3D-bioprinted scaffold produced coupling between angiogenesis and osteogenesis.

Molecular Identification of Spatially Distinct Anabolic Responses to Mechanical Loading in Murine Cortical Bone

Osteoporosis affects over 200 million women worldwide, one-third of whom are predicted to suffer from an osteoporotic fracture in their lifetime. The most promising anabolic drugs involve administration of expensive antibodies. Because mechanical loading stimulates bone formation, our current data, using a mouse model, replicates the anabolic effects of loading in humans and may identify novel pathways amenable to oral treatment. Murine tibial compression produces axially varying deformations along the cortical bone, inducing highest strains at the mid-diaphysis and lowest at the metaphyseal shell. To test the hypothesis that load-induced transcriptomic responses at different axial locations of cortical bone would vary as a function of strain magnitude, we loaded the left tibias of 10-week-old female C57Bl/6 mice in vivo in compression, with contralateral limbs as controls. Animals were euthanized at 1, 3, or 24 hours post-loading or loaded for 1 week (n = 4-5/group). Bone marrow and cancellous bone were removed, cortical bone was segmented into the metaphyseal shell, proximal diaphysis, and mid-diaphysis, and load-induced differential gene expression and enriched biological processes were examined for the three segments. At each time point, the mid-diaphysis (highest strain) had the greatest transcriptomic response. Similarly, biological processes regulating bone formation and turnover increased earlier and to the greatest extent at the mid-diaphysis. Higher strain induced greater levels of osteoblast and osteocyte genes, whereas expression was lower in osteoclasts. Among the top differentially expressed genes at 24-hours post-loading, 17 had known functions in bone biology, of which 12 were present only in osteoblasts, 3 exclusively in osteoclasts, and 2 were present in both cell types. Based on these results, we conclude that murine tibial loading induces spatially unique transcriptomic responses correlating with strain magnitude in cortical bone. © 2022 American Society for Bone and Mineral Research (ASBMR).

Transcriptome profiling of Toona ciliata young stems in response to Hypsipyla robusta Moore

Toona ciliata is a traditional woody plant that can be used as a medicinal material in China. The extracts of its roots, stems, leaves, and flowers all have a wide range of bioactive compounds. However, T. ciliata has been facing an unresolved pest problem caused by Hypsipyla robusta Moore (HRM), which seriously affects its growth and development. In this study, the expression level of TcMYB3 gene reached the maximum (28-fold) at 12 h and transcriptome sequencing of young stems eaten by HRM for 0, 3, 12, and 21 h were performed. A large number of differentially expressed genes (DEGs) were identified including jointly up-regulated genes (263) and down-regulated genes (378). JA synthesis and signaling transduction, terpene biosynthesis, and MAPKs signaling pathway were analyzed in depth and found that TcOPR3, TcJAR1, TcJAZs, and TcTPS9 genes possessed anti-insect potential. Moreover, MYB and ERF transcription factor (TF) families were significantly strengthened to the point that they may participate in induced defense mechanisms in T. ciliata. These data not only provide insights into the molecular mechanisms in resistance of T. ciliata to HRM but also helps to explore the new biocontrol strategies against insects in eco-friendly woody plants.

Vitamin K2 Improves Osteogenic Differentiation by Inhibiting STAT1 via the Bcl-6 and IL-6/JAK in C3H10 T1/2 Clone 8 Cells

Osteogenic activity of vitamin K2 (VK2), a small molecular nutrient, has been suggested. However, the underlying mechanisms have not been fully elucidated. Therefore, this study aimed to explore the mechanisms by which VK2 promotes osteogenic differentiation. The effects of VK2 on osteogenic differentiation indicators were determined in C3H10 T1/2 clone 8 cells. The RNA-seq analysis was used to explore the hypothesis that VK2 promotes osteogenic differentiation. Small interfering RNA (siRNA) assay and plasmid transfection assay were used to determine the potential role of VK2 in the modulation of Bcl-6/STAT axis and IL-6/JAK/STAT signaling pathway. VK2 significantly increased alkaline phosphatase (ALP) activity, ALP, osteocalcin (OCN), and RUNX2 abundance, and RUNX2 protein expression. RNA-seq analysis showed that there were 314 differentially expressed genes (DEGs) upregulated and 1348 DEGs downregulated by VK2. PPI analysis determined the top 10 hub genes upregulated or downregulated by VK2. Overexpression of Bcl-6 increased osteogenic differentiation and decreased expression of STAT1. Administration with VK2 restored the inhibition by siBcl-6 in osteogenic differentiation. Knockdown of IL-6 decreased the mRNA levels of genes associated with the JAK/STAT signaling pathway, and increased markers of osteoblast differentiation. Furthermore, treatment with VK2 improved inhibition in osteogenic differentiation and decreased enhancement of JAK/STAT signaling pathway related genes by overexpression of IL-6. Our study suggests that VK2 could improve osteogenic differentiation via the Bcl-6/STAT axis and IL-6/JAK/STAT signaling pathway.

Transcriptomics analysis of cashmere fineness functional genes

Liaoning cashmere goat (LCG) is a famous cashmere goat breed in China. Cashmere fineness, as an important index to evaluate cashmere quality, is also one of the problems to be improved for Liaoning cashmere goats. Transcriptome studies all mRNA transcribed by a specific tissue or cell in a certain period. It is a key link in the study of gene expression regulation. It plays an important role in the analysis of biological growth and disease. Transcriptome is spatio-temporal specific, that is, gene expression varies in different tissues or at different times. Three coarser and three fine LCG skin samples were sequenced by RNA-seq technology, and a total of 427 differentially expressed genes were obtained, including 291 up-regulated genes and 136 down-regulated genes. In the experiment, we screened out 16 genes that had significant differences in the expression of coarse and fine cashmere of Liaoning cashmere goats, so it was inferred that these 16 genes might have regulatory effects on cashmere fineness. Moreover, GO gene set enrichment analysis revealed that differential genes mainly consist of immune response, MHC protein complex, Heme binding and other pathways. KEGG analysis showed that transplant-versus-host disease and allograft rejection were the main pathways of differential genes.

Transcriptome Analysis Revealed the Symbiosis Niche of 3D Scaffolds to Accelerate Bone Defect Healing

Three dimension (3D) printed scaffolds have been shown to be superior in promoting tissue repair, but the cell-level specific regulatory network activated by 3D printing scaffolds with different material components to form a symbiosis niche have not been systematically revealed. Here, three typical 3D printed scaffolds, including natural polymer hydrogel (gelatin-methacryloyl, GelMA), synthetic polymer material (polycaprolactone, PCL), and bioceramic (β-tricalcium phosphate, β-TCP), are fabricated to explore the regulating effect of the symbiotic microenvironment during bone healing. Enrichment analysis show that hydrogel promotes tissue regeneration and reconstruction by improving blood vessel generation by enhancing oxygen transport and red blood cell development. The PCL scaffold regulates cell proliferation and differentiation by promoting cellular senescence, cell cycle and deoxyribonucleic acid (DNA) replication pathways, accelerating the process of endochondral ossification, and the formation of callus. The β-TCP scaffold can specifically enhance the expression of osteoclast differentiation and extracellular space pathway genes to promote the differentiation of osteoclasts and promote the process of bone remodeling. In these processes, specific biomaterial properties can be used to guide cell behavior and regulate molecular network in the symbiotic microenvironment to reduce the barriers of regeneration and repair.

Advances in single-cell sequencing and its application to musculoskeletal system research

In recent years, single-cell sequencing (SCS) technologies have continued to advance with improved operating procedures and reduced cost, leading to increasing practical adoption among researchers. These emerging technologies have superior abilities to analyse cell heterogeneity at a single-cell level, which have elevated multi-omics research to a higher level. In some fields of research, application of SCS has enabled many valuable discoveries, and musculoskeletal system offers typical examples. This article reviews some major scientific issues and recent advances in musculoskeletal system. In addition, combined with SCS technologies, the research of cell or tissue heterogeneity in limb development and various musculoskeletal system clinical diseases also provides new possibilities for treatment strategies. Finally, this article discusses the challenges and future development potential of SCS and recommends the direction of future applications of SCS to musculoskeletal medicine.

Mineralized Cartilage and Bone-Like Tissues in Chondrichthyans Offer Potential Insights Into the Evolution and Development of Mineralized Tissues in the Vertebrate Endoskeleton

The impregnation of biominerals into the extracellular matrix of living organisms, a process termed biomineralization, gives rise to diverse mineralized (or calcified) tissues in vertebrates. Preservation of mineralized tissues in the fossil record has provided insights into the evolutionary history of vertebrates and their skeletons. However, current understanding of the vertebrate skeleton and of the processes underlying its formation is biased towards biomedical models such as the tetrapods mouse and chick. Chondrichthyans (sharks, skates, rays, and chimaeras) and osteichthyans are the only vertebrate groups with extant (living) representatives that have a mineralized skeleton, but the basal phylogenetic position of chondrichthyans could potentially offer unique insights into skeletal evolution. For example, bone is a vertebrate novelty, but the internal supporting skeleton (endoskeleton) of extant chondrichthyans is commonly described as lacking bone. The molecular and developmental basis for this assertion is yet to be tested. Subperichondral tissues in the endoskeleton of some chondrichthyans display mineralization patterns and histological and molecular features of bone, thereby challenging the notion that extant chondrichthyans lack endoskeletal bone. Additionally, the chondrichthyan endoskeleton demonstrates some unique features and others that are potentially homologous with other vertebrates, including a polygonal mineralization pattern, a trabecular mineralization pattern, and an unconstricted perichordal sheath. Because of the basal phylogenetic position of chondrichthyans among all other extant vertebrates with a mineralized skeleton, developmental and molecular studies of chondrichthyans are critical to flesh out the evolution of vertebrate skeletal tissues, but only a handful of such studies have been carried out to date. This review discusses morphological and molecular features of chondrichthyan endoskeletal tissues and cell types, ultimately emphasizing how comparative embryology and transcriptomics can reveal homology of mineralized skeletal tissues (and their cell types) between chondrichthyans and other vertebrates.

Perspective of the GEMSTONE Consortium on Current and Future Approaches to Functional Validation for Skeletal Genetic Disease Using Cellular, Molecular and Animal-Modeling Techniques

The availability of large human datasets for genome-wide association studies (GWAS) and the advancement of sequencing technologies have boosted the identification of genetic variants in complex and rare diseases in the skeletal field. Yet, interpreting results from human association studies remains a challenge. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary. Multiple unknowns exist for putative causal genes, including cellular localization of the molecular function. Intermediate traits ("endophenotypes"), e.g. molecular quantitative trait loci (molQTLs), are needed to identify mechanisms of underlying associations. Furthermore, index variants often reside in non-coding regions of the genome, therefore challenging for interpretation. Knowledge of non-coding variance (e.g. ncRNAs), repetitive sequences, and regulatory interactions between enhancers and their target genes is central for understanding causal genes in skeletal conditions. Animal models with deep skeletal phenotyping and cell culture models have already facilitated fine mapping of some association signals, elucidated gene mechanisms, and revealed disease-relevant biology. However, to accelerate research towards bridging the current gap between association and causality in skeletal diseases, alternative in vivo platforms need to be used and developed in parallel with the current -omics and traditional in vivo resources. Therefore, we argue that as a field we need to establish resource-sharing standards to collectively address complex research questions. These standards will promote data integration from various -omics technologies and functional dissection of human complex traits. In this mission statement, we review the current available resources and as a group propose a consensus to facilitate resource sharing using existing and future resources. Such coordination efforts will maximize the acquisition of knowledge from different approaches and thus reduce redundancy and duplication of resources. These measures will help to understand the pathogenesis of osteoporosis and other skeletal diseases towards defining new and more efficient therapeutic targets.

Genetic variability affects the skeletal response to immobilization in founder strains of the diversity outbred mouse population

Mechanical unloading decreases bone volume and strength. In humans and mice, bone mineral density is highly heritable, and in mice the response to changes in loading varies with genetic background. Thus, genetic variability may affect the response of bone to unloading. As a first step to identify genes involved in bones' response to unloading, we evaluated the effects of unloading in eight inbred mouse strains: C57BL/6J, PWK/PhJ, WSB/EiJ, A/J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HlLtJ, and CAST/EiJ. C57BL/6J and NOD/ShiLtJ mice had the greatest unloading-induced loss of diaphyseal cortical bone volume and strength. NZO/HlLtJ mice had the greatest metaphyseal trabecular bone loss, and C57BL/6J, WSB/EiJ, NOD/ShiLtJ, and CAST/EiJ mice had the greatest epiphyseal trabecular bone loss. Bone loss in the epiphyses displayed the highest heritability. With immobilization, mineral:matrix was reduced, and carbonate:phosphate and crystallinity were increased. A/J mice displayed the greatest unloading-induced loss of mineral:matrix. Changes in gene expression in response to unloading were greatest in NOD/ShiLtJ and CAST/EiJ mice. The most upregulated genes in response to unloading were associated with increased collagen synthesis and extracellular matrix formation. Our results demonstrate a strong differential response to unloading as a function of strain. Diversity outbred (DO) mice are a high-resolution mapping population derived from these eight inbred founder strains. These results suggest DO mice will be highly suited for examining the genetic basis of the skeletal response to unloading.

•

Mouse strain affects bone's response to immobilization.

•

Magnitude of bone loss from immobilization is heritable.

•

Bone transcriptomic response to immobilization is influenced by genetic variation.

The Aqueous Extract of Eucommia Leaves Promotes Proliferation, Differentiation, and Mineralization of Osteoblast-Like MC3T3-E1 Cells

Eucommia leaves are dry leaves of Eucommia ulmoides which have long been considered as a functional health food for the treatment of hypertension, hypercholesterolemia, fatty liver, and osteoporosis. With the recent development of Chinese medicine, Eucommia leaves are widely used for tonifying the kidneys and strengthening bone. However, the specific molecular mechanism of Eucommia leaves for strengthening bone remains largely unknown. Osteoblasts are the main functional cells of bone formation; thus, it is essential to study the effect of Eucommia leaves on osteoblasts to better understand their mechanism of action. In the present study, we prepared an aqueous extract of Eucommia leaves (ELAE) and determined its content by high-performance liquid chromatography (HPLC). The effects of ELAE on MC3T3-E1 cells were investigated by CCK-8 assay, alkaline phosphatase (ALP), and Alizarin red S staining assays, combined with RNA sequencing (RNA-seq) and qRT-PCR validation. We demonstrated that ELAE had a significant promoting effect on the proliferation of MC3T3-E1 cells and significantly enhanced extracellular matrix synthesis and mineralization, which were achieved by regulating various functional genes and related signaling pathways. ELAE significantly increased the expression level of genes promoting cell proliferation, such as Rpl10a, Adnp, Pex1, Inpp4a, Frat2, and Pcdhga1, and reduced the expression level of genes inhibiting cell proliferation, such as Npm1, Eif3e, Cbx3, Psmc6, Fgf7, Fxr1, Ddx3x, Mbnl1, and Cdc27. In addition, ELAE increased the expression level of gene markers in osteoblasts, such as Col5a2, Ubap2l, Dkk3, Foxm1, Col16a1, Col12a1, Usp7, Col4a6, Runx2, Sox4, and Bmp4. Taken together, our results suggest that ELAE could promote osteoblast proliferation, differentiation, and mineralization and prevent osteoblast apoptosis. These findings not only increase our understanding of ELAE on the regulation of bone development but also provide a possible strategy to further study the prevention and treatment of osteogenic related diseases by ELAE.

The enzymatic hydrolysates from deer sinew promote MC3T3-E1 cell proliferation and extracellular matrix synthesis by regulating multiple functional genes

Background

Deer Sinew serves as a medicinal food, and has been used for treating skeletal diseases, especially bone diseases in a long history. Thus, it could become an alternative option for the prevention and therapeutic remedy of bone-related diseases. In our previous study, we established an optimal extraction process of the enzymatic hydrolysates from Chinese Sika deer sinews (DSEH), and we demonstrated that DSEH significantly promoted the proliferation of MC3T3-E1 cells (an osteoblast-like cell line) with a certain dose-effect relationship. However, the precise molecular mechanism of deer sinew in regulating bone strength is still largely unknown. The aim of this study was to explore the underlying molecular mechanism of DSEH on MC3T3-E1 cells proliferation and extracellular matrix synthesis.

Methods

Preparation and quality control were performed as previously described. The effect of DSEH at different administrated concentrations on cell proliferation was measured using both CCK-8 and MTT assays, and the capacity of DSEH on extracellular matrix synthesis was detected by Alizarin red staining and quantification. The gene expression pattern change of MC3T3-E1 cells under the treatment of DSEH was investigated by RNA-seq analysis accompanied with validation methods.

Results

We demonstrated that DSEH promoted MC3T3-E1 cell proliferation and extracellular matrix synthesis by regulating multiple functional genes. DSEH significantly increased the expression levels of genes that promoted cell proliferation such as Gstp1, Timp1, Serpine1, Cyr61, Crlf1, Thbs1, Ctgf, P4ha2, Sod3 and Nqo1. However, DSEH significantly decreased the expression levels of genes that inhibited cell proliferation such as Mt1, Cdc20, Gas1, Nrp2, Cmtm3, Dlk2, Sema3a, Rbm25 and Hspb6. Furthermore, DSEH mildly increased the expression levels of osteoblast gene markers.

Conclusions

Our findings suggest that DSEH facilitate MC3T3-E1 cell proliferation and extracellular matrix synthesis to consolidate bone formation and stability, but prevent MC3T3-E1 cells from oxidative stress-induced damage, apoptosis and further differentiation. These findings deepened the current understanding of DSEH on regulating bone development, and provided theoretical support for the discovery of optional prevention and treatment for bone-related diseases.

New Targets and Emergent Therapies for Osteoporosis

The 11 existing FDA-approved osteoporosis drug treatments include hormone replacement therapy, 2 SERMs (raloxifene and bazedoxifene), 5 inhibitors of bone-resorbing osteoclasts (4 bisphosphonates and anti-RANKL denosumab), 2 parathyroid hormone analogues (teriparatide and abaloparatide), and 1 WNT signaling enhancer (romosozumab). These therapies are effective and provide multiple options for patients and physicians. As the genomic revolution continues, potential novel targets for future drug development are identified. This review takes a wide perspective to describe potentially rewarding topics to explore, including knowledge of genes and pathways involved in bone cell metabolism, the utility of animal models, targeting drugs to bone, and ongoing advances in drug design and delivery.

High Fidelity of Mouse Models Mimicking Human Genetic Skeletal Disorders

The 2019 International Skeletal Dysplasia Society nosology update lists 441 genes for which mutations result in rare human skeletal disorders. These genes code for enzymes (33%), scaffolding proteins (18%), signal transduction proteins (16%), transcription factors (14%), cilia proteins (8%), extracellular matrix proteins (5%), and membrane transporters (4%). Skeletal disorders include aggrecanopathies, channelopathies, ciliopathies, cohesinopathies, laminopathies, linkeropathies, lysosomal storage diseases, protein-folding and RNA splicing defects, and ribosomopathies. With the goal of evaluating the ability of mouse models to mimic these human genetic skeletal disorders, a PubMed literature search identified 260 genes for which mutant mice were examined for skeletal phenotypes. These mouse models included spontaneous and ENU-induced mutants, global and conditional gene knockouts, and transgenic mice with gene over-expression or specific base-pair substitutions. The human X-linked gene ARSE and small nuclear RNA U4ATAC, a component of the minor spliceosome, do not have mouse homologs. Mouse skeletal phenotypes mimicking human skeletal disorders were observed in 249 of the 260 genes (96%) for which comparisons are possible. A supplemental table in spreadsheet format provides PubMed weblinks to representative publications of mutant mouse skeletal phenotypes. Mutations in 11 mouse genes (Ccn6, Cyp2r1, Flna, Galns, Gna13, Lemd3, Manba, Mnx1, Nsd1, Plod1, Smarcal1) do not result in similar skeletal phenotypes observed with mutations of the homologous human genes. These discrepancies can result from failure of mouse models to mimic the exact human gene mutations. There are no obvious commonalities among these 11 genes. Body BMD and/or radiologic dysmorphology phenotypes were successfully identified for 28 genes by the International Mouse Phenotyping Consortium (IMPC). Forward genetics using ENU mouse mutagenesis successfully identified 37 nosology gene phenotypes. Since many human genetic disorders involve hypomorphic, gain-of-function, dominant-negative and intronic mutations, future studies will undoubtedly utilize CRISPR/Cas9 technology to examine transgenic mice having genes modified to exactly mimic variant human sequences. Mutant mice will increasingly be employed for drug development studies designed to treat human genetic skeletal disorders.

SIGNIFICANCE

Great progress is being made identifying mutant genes responsible for human rare genetic skeletal disorders and mouse models for genes affecting bone mass, architecture, mineralization and strength. This review organizes data for 441 human genetic bone disorders with regard to heredity, gene function, molecular pathways, and fidelity of relevant mouse models to mimic the human skeletal disorders. PubMed weblinks to citations of 249 successful mouse models are provided.