FGF2 and FAM201A affect the development of osteonecrosis of the femoral head after femoral neck fracture

Investigation of the association between the Toll-like receptor 1 rs4833095 variation and gastric adenocarcinoma recurrence

BACKGROUND

Toll-like receptors (TLRs) are a family of transmembrane receptors that play key roles in identifying invading pathogens and activating innate immunity. TLR1 has been reported to be associated with the risk of gastric cancer (GC) but that was based on only a simple statistical analysis.

METHODS

We genotyped the TLR1 in 526 GC patients to investigate the association between the variation and gastric cancer survival by the multiplex polymerase chain reaction and sequencing method. The rs4833095 variation (chr4:38798089 [GRCh38. p14], T > C) in the TLR1 gene was genotyped in 526 patients who underwent GC resection. The associations between genotype, survival, and recurrence were investigated. The potential role of TLR1 in stomach cancer was investigated using clinical data from formalin-fixed, paraffin-embedded tissue samples.

RESULTS

Patients with the T/C and C/C genotypes of rs4833095 had a lower risk of recurrence than those with the T/T genotype. Recurrence-free periods were substantially longer in patients with the T/C or C/C genotypes (22.6 and 22.3 months, respectively) than in those with the T/T genotype (20.7 months). Patients with the T/C or C/C genotype, low expression levels of VEGF1, high expression levels of ERBB2 and ERCC1, the absence of cancer nodules, a tumor size of less than 5 cm, and poor differentiation had a considerably reduced risk of recurrence.

CONCLUSIONS

TLR1 rs4833095 was correlated with the postresection prognosis of patients with gastric cancer, suggesting that TLR1 may have a role in the onset or progression of gastric cancer.

SOX9 Inhibits the Progression of Osteonecrosis of the Femoral Head via the Activation of the Wnt/Beta-Catenin Pathway

Background: In this study, we aimed to explore whether the SRY-box transcription factor 9 (SOX9) can play protective roles against the occurrence and development of osteonecrosis of the femoral head (ONFH) by regulating the proliferation, apoptosis, and osteogenic differentiation of human bone marrow stromal cells (hBMSCs) via the Wnt/β-catenin pathway.Methods: We used 1600 mg of glucocorticoid (GC) to induce hBMSCs to establish an ONFH cell model and performed various experiments. Reverse transcription-quantitative polymerase chain reaction and western blotting assays were used to determine the expression levels of SOX9 and osteoblast markers, such as the RUNX family transcription factor 2 (RUNX2), alkaline phosphatase (ALP), osterix, Wnt3a, and β-catenin. An ALP detection kit was used to measure the ALP activity. Flow cytometry and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assays were performed to determine the cell viability.Results: GC treatment decreased the expression levels of RUNX2, ALP, and osterix, suppressed ALP activity, and inhibited SOX9 expression. SOX9 overexpression promoted GC-induced cell proliferation and decreased cell apoptosis. Additionally, hBMSCs were transfected with SOX9-small interfering RNA in GC treatment, and SOX9 knockdown was found to suppress the osteogenic differentiation of cells and decrease their viability.Conclusion: Our results revealed that SOX9 is related to the Wnt/β-catenin pathway in ONFH. Moreover, SOX9 participated in ONFH development by activating the Wnt/β-catenin pathway.

Novel Gene Polymorphisms for Stable Warfarin Dose in a Korean Population: Genome-Wide Association Study

Warfarin has a narrow therapeutic window and high intra- and inter-individual variability. Considering that many published papers on genotype-guided dosing are derived from European populations, the aim of this study was to investigate novel genetic variants associated with the variability of stable warfarin dose in the Korean population with cardiac valve replacement, using the GWAS approach. This retrospective cohort study was performed from January 1982 to December 2020 at the Severance Cardiovascular Hospital of Yonsei University College of Medicine. GWAS was performed to identify associations between genotypes and the warfarin maintenance dose, by comparing the allele frequency of genetic variants between individuals. Then, the extent of genetic and non-genetic factors on the dose variability was determined by multivariable regression analysis. The study enrolled 214 participants, and the most robust signal cluster was detected on chromosome 16 around VKORC1. Followed by VKORC1, three novel variants (NKX2-6 rs310279, FRAS1 rs4386623, and FAM201A rs1890109) showed an association with stable warfarin dose requirement in univariate analysis. The algorithm was constructed by using multivariable analysis that includes genetic and non-genetic factors, and it could explain 58.5% of the variations in stable warfarin doses. In this variability, VKORC1 rs9934438 and FRAS1 rs4386623 accounted for 33.0% and 9.9%, respectively. This GWAS analysis identified the fact that three novel variants (NKX2-6 rs310279, FRAS1 rs4386623, and FAM201A rs1890109) were associated with stable warfarin doses. Additional research is necessary to validate the results and establish personalized treatment strategies for the Korean population.

Construction of ceRNA network based on RNA-seq for identifying prognostic lncRNA biomarkers in Perthes disease

Introduction: Legg-Calvé-Perthes disease or Perthes disease is a condition that occurs in children aged 2 to 15 years, and is characterized by osteonecrosis of the femoral head, which results in physical limitations. Despite ongoing research, the pathogenesis and molecular mechanisms underlying the development of Perthes disease remain unclear. In order to obtain further insights, the expression patterns of long non-coding RNAs (lncRNAs), miRNAs, and mRNAs in a rabbit model of Perthes disease were analyzed in this study by transcriptome sequencing. Methods and results: The results of RNA-seq analyses revealed that 77 lncRNAs, 239 miRNAs, and 1027 mRNAs were differentially expressed in the rabbit model. This finding suggested that multiple genetic pathways are involved in the development of Perthes disease. A weighted gene co-expression network analysis (WGCNA) network was subsequently constructed using the differentially expressed mRNAs (DEmRNAs), and network analysis revealed that the genes associated with angiogenesis and platelet activation were downregulated, which was consistent with the findings of Perthes disease. A competing endogenous RNA (ceRNA) network was additionally constructed using 29 differentially expressed lncRNAs (including HIF3A and LOC103350994), 28 differentially expressed miRNAs (including ocu-miR-574-5p and ocu-miR-324-3p), and 76 DEmRNAs (including ALOX12 and PTGER2). Disscusion: The results obtained herein provide novel perspectives regarding the pathogenesis and molecular mechanisms underlying the development of Perthes disease. The findings of this study can pave the way for the development of effective therapeutic strategies for Perthes disease in future.

Differential Expression of Non-Coding RNAs in Stem Cell Development and Therapeutics of Bone Disorders

Stem cells' self-renewal and multi-lineage differentiation are regulated by a complex network consisting of signaling factors, chromatin regulators, transcription factors, and non-coding RNAs (ncRNAs). Diverse role of ncRNAs in stem cell development and maintenance of bone homeostasis have been discovered recently. The ncRNAs, such as long non-coding RNAs, micro RNAs, circular RNAs, small interfering RNA, Piwi-interacting RNAs, etc., are not translated into proteins but act as essential epigenetic regulators in stem cells' self-renewal and differentiation. Different signaling pathways are monitored efficiently by the differential expression of ncRNAs, which function as regulatory elements in determining the fate of stem cells. In addition, several species of ncRNAs could serve as potential molecular biomarkers in early diagnosis of bone diseases, including osteoporosis, osteoarthritis, and bone cancers, ultimately leading to the development of new therapeutic strategies. This review aims to explore the specific roles of ncRNAs and their effective molecular mechanisms in the growth and development of stem cells, and in the regulation of osteoblast and osteoclast activities. Furthermore, we focus on and explore the association of altered ncRNA expression with stem cells and bone turnover.

Transcriptomic and metabolomic insights into the role of the flgK gene in the pathogenicity of Pseudomonas plecoglossicida to orange-spotted grouper (Epinephelus coioides)

Pseudomonas plecoglossicida is the pathogen responsible for visceral white spot disease in large yellow croaker (Larimichthys crocea) and orange-spotted grouper (Epinephelus coioides). Previously, RNA sequencing showed that P. plecoglossicida flgK gene expression was significantly up-regulated in orange-spotted grouper spleens during infection. To explore the role of flgK in P. plecoglossicida pathogenicity, RNA interference (RNAi) was performed to silence the P. plecoglossicida flgK gene, and the mutant (flgK-RNAi strain) with the best silencing efficiency (89.40%) was chosen for further study. Results showed that flgK gene silencing significantly attenuated P. plecoglossicida motility, adhesion, and biofilm formation. Compared to those fish infected with the wild-type strain of P. plecoglossicida, orange-spotted grouper infected with the flgK-RNAi strain showed a 55% increase in the survival rate and a one-day delay in time of first death, with fewer pathogens in the spleen and fewer white spots on the spleen surface. RNAi of flgK significantly affected the transcriptome and metabolome of the spleen in infected orange-spotted grouper. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the C-type lectin receptor signaling pathway was the most significantly changed immune-related pathway and the mitogen-activated protein kinase (MAPK) signaling pathway was related to multiple immune-related pathways. Furthermore, arginine biosynthesis and glycerophospholipid metabolism were the most significantly changed metabolism-related pathways. These findings suggest that flgK is a virulence gene of P. plecoglossicida. Furthermore, flgK appears to be involved in the regulation of motility, adhesion, and biofilm formation in P. plecoglossicida, as well as in the regulation of inflammatory and immune responses of orange-spotted grouper to P. plecoglossicida infection.

Therapeutic perspectives of exosomes in glucocorticoid-induced osteoarthrosis

Exosomes are widely involved in a variety of physiological and pathological processes. These important roles are also hidden in the physiological processes related to bone. Chondrocytes, osteoblasts, synovial fibroblasts, and bone marrow mesenchymal stem cells produce and secrete exosomes, thereby affecting the biology process of target cells. Furthermore, in the primary pathogenesis of osteoarthrosis induced by steroid hormones, mainly involve glucocorticoid (GC), the exosomes have also widely participated. Therefore, exosomes may also play an important role in glucocorticoid-induced osteoarthrosis and serve as a promising treatment for early intervention of osteoarthrosis in addition to playing a regulatory role in malignant tumors. This review summarizes the previous results on this direction, systematically combs the role and therapeutic potential of exosomes in GC-induced osteoarthrosis, discusses the potential role of exosomes in the treatment and prevention of GC-induced osteoarthrosis, and reveals the current challenges we confronted.

FAM201A, a long noncoding RNA potentially associated with atrial fibrillation identified by ceRNA network analyses and WGCNA

Background

Being the most common arrhythmia in clinic, atrial fibrillation (AF) causes various comorbidities to patients such as heart failure and stroke. LncRNAs were reported involved in pathogenesis of AF, yet, little is known about AF-associated lncRNAs. The present study aims to explore lncRNAs associated with AF susceptibility based on competing endogenous RNA (ceRNA) network analysis and weighted gene co-expression network analysis (WGCNA).

Methods

GSE41177 and GSE79768 datasets were obtained from the Gene Expression Omnibus (GEO) database. Competing endogenous RNA (ceRNA) network analysis was performed using GSE41177. Differentially expressed lncRNAs (DElncRNAs), mRNAs (DEmRNAs) between AF patients and patients with sinus rhythm (SR) were identified from GSE41177 using R software. Then, the ceRNA network was constructed based on DElncRNAs, the predicted target miRNAs and DEmRNAs. Weighted gene co-expression network analysis (WGCNA) was performed using GSE79768 to validate the AF-related lncRNAs identified from GSE41177. LncRNA modules and crucial lncRNAs relevant to AF and were identified.

Results

In summary, 18 DElncRNAs and 350 DEmRNAs were found between AF patients and SR patients. A total of 5 lncRNAs, 10 miRNAs, and 21 mRNAs were contained in the final ceRNA network. Taking into consideration both the ceRNA theory and inference scores from the comparative toxicogenomics database (CTD) database, the ceRNA axis FAM201A-miR-33a-3p-RAC3 was identified as mostly relevant to AF susceptibility. FAM201A (Gene significance, GS = − 0.62; Module membership, MM = 0.75) was also proved in the blue module, which was identified most highly relevant with AF by WGCNA.

Conclusions

These results demonstrated that decreased expression of FAM201A might be associated with susceptibility of AF. Working as the ceRNA to regulate RAC3 might be one function of FAM201A in AF susceptibility, which requires further exploration in future research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01232-w.

Positive feedback loop of lncRNA FAM201A/miR‑146a‑5p/POU2F1 regulates IL‑1β‑induced chondrocyte injury in vitro

Numerous studies have previously demonstrated that long non-coding RNAs (lncRNAs) serve an important regulatory role in osteoarthritis (OA). In particular, the lncRNA family with sequence similarity 201 member A (FAM201A) was previously found to be downregulated in necrotic femoral head samples. However, the role of FAM201A in IL-1β-induced chondrocyte injury remains unclear. It was hypothesized that FAM201A may exert a protective effect on IL-1β-induced chondrocyte injury in OA by sponging microRNAs (miRNAs/miRs). The purpose of the present study was to explore the role and molecular mechanism of FAM201A in IL-1β-induced chondrocyte injury. A model of OA was established by stimulation C-28/I2 cell with IL-1β in vitro. The expression levels of FAM201A following IL-1β-induced chondrocyte injury were detected via reverse transcription-quantitative PCR. Luciferase reporter assay was used to assess the possible associations among FAM201A, miR-146a-5p and POU class 2 homeobox 1 (POU2F1). Chromatin immunoprecipitation assay was performed to analyze the interaction between POU2F1 and miR-146a-5p. ELISA, TUNEL and western blotting were performed to measure the level of inflammation, lactate dehydrogenase release, apoptosis and the expression of apoptosis-related proteins (Bcl-2, Bax, cleaved caspase 3 and cleaved caspase 9), respectively. The expression levels of FAM201A were found to be downregulated following IL-1β-induced chondrocyte injury. Overexpression of FAM201A exerted a protective effect against IL-1β-induced chondrocyte injury. In addition, FAM201A could upregulate the expression levels of POU2F1 by sponging miR-146a-5p. Further experiments revealed that POU2F1 could bind to the promoter region of FAM201A and subsequently regulate the expression levels of POU2F1, indicating a role for the FAM201A/miR-146a-5p/POU2F1 positive feedback loop in IL-1β-induced chondrocyte injury. The present study revealed the protective effects of the FAM201A/miR-146a-5p/POU2F1 positive feedback loop on IL-1β-induced chondrocyte injury and provided a potential therapeutic target for OA.

FAM201A knockdown inhibits proliferation and invasion of lung adenocarcinoma cells by regulating miR-7515/GLO1 axis

Lung adenocarcinoma (LUAD) is the most important histological type of lung cancer. We aimed to identify the role of long noncoding RNA family with sequence similarity 201-member A (FAM201A) in the occurrence and development of LUAD. The expressions of FAM201A in LUAD tissues and cells were determined via reverse transcription-quantitative polymerase chain reaction. The effects of FAM201A knockdown on LUAD cell malignant phenotypes were examined by cell counting kit-8, 5-ethynyl-2'-deoxyuridine, flow cytometry, transwell assay and wound healing assay. The underlying mechanism by which FAM201A regulated LUAD progression was also studied. Nude mice LUAD xenograft model was constructed, to explore the in vivo effect of FAM201A. Our results showed that the FAM201A expression in LUAD tissues and cell lines was notably higher than normal tissues and cells. Downregulation of FAM201A suppressed the cell proliferation, migration and invasion and promoted the cell apoptosis in LUAD cells. While, FAM201A overexpression showed tumorigenesis effect on LUAD cells. Moreover, we demonstrated that FAM201A affected LUAD progression via targeting miR-7515 to promote GLO1 expression. FAM201A downregulation also suppressed LUAD development in vivo experiment. Our results indicated that FAM201A was an oncogene in LUAD and might be a novel therapeutic target for LUAD.

ATF6 aggravates angiogenesis-osteogenesis coupling during ankylosing spondylitis by mediating FGF2 expression in chondrocytes

Although angiogenesis-osteogenesis coupling is important in ankylosing spondylitis (AS), therapeutic agents targeting the vasculature remain elusive. Here, we identified activating transcription factor 6 (ATF6) as an important regulator of angiogenesis in the pathogenesis of AS. First, we found that ATF6 and fibroblast growth factor 2 (FGF2) levels were higher in SKG mice and in cartilage of pateints with AS1. The proangiogenic activity of human chondrocytes was enhanced by the activation of the ATF6-FGF2 axis following 7 days of stimulation with inflammatory factors, e.g., tumor necrosis factor alpha (TNF-α), interferon-γ (IFN-γ) or interleukin-17 (IL-17). Mechanistically, ATF6 interacted with the FGF2 promotor and promoted its transcription. Treatment with the ATF6 inhibitor Ceapin-A7 inhibited angiogenesis in vitro and angiogenesis-osteogenesis coupling in vivo. ATF6 may aggravate angiogenesis-osteogenesis coupling during AS by mediating FGF2 transcription in chondrocytes, implying that ATF6 represents a promising therapeutic target for AS.

Eight-gene signature predicts recurrence in lung adenocarcinoma

BACKGROUND

Recurrence significantly influences the survival in patients with lung adenocarcinoma (LUAD). However, there are less gene signatures that predict recurrence risk of LUAD.

OBJECTIVE

We performed this study to construct a model to predict risk of recurrence in LUAD.

METHODS

RNA-seq data from 426 patients with LUAD were downloaded from The Cancer Genome Atlas (TCGA) and were randomly assigned into the training (n= 213) and validation set (n= 213). Differentially expressed genes (DEGs) between recurrent and non-recurrent tumors in the training set were identified. Recurrence-associated DEGs were selected using multivariate Cox regression analysis. The recurrence risk model that identifies patients at low and high risk for recurrence was constructed, followed by the validation of its performance in the validation set and a microarray dataset.

RESULTS

In total, 378 DEGs, including 20 recurrence-associated DEGs, were identified between the recurrent and non-recurrent tumors in the training set. The signatures of 8 genes (including AZGP1, INPP5J, MYBPH, SPIB, GUCA2A, HTR1B, SLC15A1 and TNFSF11) were used to construct the prognostic model to assess the risk of recurrence. This model indicated that patients with high risk scores had shorter recurrence-free survival time compared with patients with low risk scores. ROC curve analysis of this model showed it had high predictive accuracy (AUC > 0.8) to predict LUAD recurrence in the TCGA cohort (the training and validation sets) and GSE50081 dataset. This prognostic model showed high predictive power and performance in predicting recurrence in LUAD.

CONCLUSION

We concluded that this model might be of great value for evaluating the risk of recurrence of LUAD in clinics.

Non-coding RNA Identification in Osteonecrosis of the Femoral Head Using Competitive Endogenous RNA Network Analysis

OBJECTIVE

To investigate the regulatory network of long non-coding RNA (lncRNA) as competing endogenous RNAs (ceRNAs) in osteonecrosis of the femoral head (ONFH).

METHODS

The gene expression profile GSE74089 of ONFH and microRNA (miRNA) expression profile of GSE89587 were obtained from the Gene Expression Omnibus (GEO) database. The GSE74089 contained four ONFH samples and four controls. The GSE89587 included 10 ONFH samples and 10 control samples. The differentially expressed lncRNAs (DE-lncRNAs) and DE-mRNAs between ONFH group and control group were identified from GSE74089 using the limma package based on criteria of adjusted P value <0.05 and |log fold change (FC)| ≥2. The DEmiRNAs between ONFH group and control group were screened from GSE89587 on the basis of adjusted P value <0.05. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway for DE-mRNAs were analyzed using DAVID 6.7 and GSEA 3.0, respectively. Coexpressed lncRNA-mRNA pairs were identified by corr.test method in R based on the criteria of adjusted P value <0.01 and |r| ≥ 0.9. A ceRNA network was constructed and visualized using cytoscape 3.7.0 by integrating the DE-lncRNA, DE-miRNA, and DEmRNA data. The key mRNAs and lncRNAs in the ceRNA network were further validated in an independent dataset of GSE123568.

RESULTS

Based on our analysis, a total of 28 DE-lncRNAs, 1403 DE-mRNAs, and 134 DE-miRNAs were identified, respectively. The DE-mRNAs were significantly enriched in the function of "skeletal system development," "collagen fibril organization," "blood vessel development," and "regulation of nervous system development." Besides, 72 KEGG pathways, including eight active pathways and 64 suppressed pathways were identified, including which immune pathway was the most significantly activated one and which ribosome-related function was the most suppressed. A co-expression network including 161 DE-mRNAs and 16 DE-lncRNAs was built. Highly connected nodes were identified among lncRNAs such as H19, C20orf203, LINC00355, SFTA3, CRNDE, CASC2, LINC00494, C9orf163, C10orf91, and LINC00301. The ceRNA network indicated that lncRNA H19 functioned as a ceRNA of hsa-miR-519b-3p and hsa-miR-296-5p in ANKH and ECHDC1 regulation; lncRNA C9orf163 functioned as a ceRNA of hsa-miR-424-5p in CCNT1 regulation. The expression trends of ANKH, CCNT1, and C9orf163 were successfully validated in independent dataset of GSE123568.

CONCLUSION

The ceRNAs of lncRNA H19- hsa-miR-519b-3p/hsa-miR-296-5p-ANKH and lncRNA c9orf163- hsa-miR-424-5p-CCNT1 might play important roles in ONFH development. Our research provided an understanding of the important role of lncRNA-related ceRNAs in ONFH.

Integrated analysis of transcriptome and proteome to explore the genes related to steroid-induced femoral head necrosis

PURPOSE

Femoral head necrosis (FHN) is a common disease of hip. However, the pathogenesis of FHN is not well understood. This study attempted to explore the potentially important genes and proteins involved in FHN.

METHODS

We integrated the transcriptomic and proteomic methods to quantitatively screen the differentially expressed genes (DEGs) and proteins (DEPs) between Control and FHN groups. Gene ontology (GO) terms and KEGG pathway enrichment analysis were used to assess the roles of DEGs and DEPs. qRT-PCR and western blot were performed to verify the key genes/proteins in FHN. CCK-8 assay was performed to measure cell viability. The protein expression of Bax and Bcl-2 were used to evaluate cell apoptosis.

RESULTS

Transcriptome and proteome studies indicated 758 DEGs and 1097 DEPs between Control and FHN groups, respectively. Cell division, extracellular exosome, and serine-type endopeptidase activity were the most common terms in biological process (BP), cellular component (CC), and molecular function (MF) enrichment, respectively. DEPs were mainly enriched in cellular process, cell, and binding for BP, CC, and MF categories, respectively. DEGs were mainly involved in PI3K-Akt pathway and DEPs were mainly focused in glycolysis/gluconeogenesis pathway. Notably, 14 down-regulated and 22 up-regulated genes/proteins were detected at both the transcript and protein level. LRG1, SERPINE2, STMN1, COL14A1, SLC37A2, and MMP2 were determined as the key genes/proteins in FHN. SERPINE2/STMN1 overexpression increased viability and decreased apoptosis of dexamethasone-treated MC3T3-E1 cells.

CONCLUSIONS

Our study investigated some pivotal regulatory genes/proteins in the pathogenesis of FHN, providing novel insight into the genes/proteins involved in FHN.

Identifying the Potential Differentially Expressed miRNAs and mRNAs in Osteonecrosis of the Femoral Head Based on Integrated Analysis

Purpose

Osteonecrosis of the femoral head is a common disease of the hip that leads to severe pain or joint disability. We aimed to identify potential differentially expressed miRNAs and mRNAs in osteonecrosis of the femoral head.

Methods

The data of miRNA and mRNA were firstly downloaded from the database. Secondly, the regulatory network of miRNAs-mRNAs was constructed, followed by function annotation of mRNAs. Thirdly, an in vitro experiment was applied to validate the expression of miRNAs and targeted mRNAs. Finally, GSE123568 dataset was used for electronic validation and diagnostic analysis of targeted mRNAs.

Results

Several regulatory interaction pairs between miRNA and mRNAs were identified, such as hsa-miR-378c-WNT3A/DACT1/CSF1, hsa-let-7a-5p-RCAN2/IL9R, hsa-miR-28-5p-RELA, hsa-miR-3200-5p-RELN, and hsa-miR-532-5p-CLDN18/CLDN10. Interestingly, CLDN10, CLDN18, CSF1, DACT1, IL9R, RCAN2, RELN, and WNT3A had the diagnostic value for osteonecrosis of the femoral head. Wnt signaling pathway (involved WNT3A), chemokine signaling pathway (involved RELA), focal adhesion and ECM-receptor interaction (involved RELN), cell adhesion molecules (CAMs) (involved CLDN18 and CLDN10), cytokine-cytokine receptor interaction, and hematopoietic cell lineage (involved CSF1 and IL9R) were identified.

Conclusion

The identified differentially expressed miRNAs and mRNAs may be involved in the pathology of osteonecrosis of the femoral head.

Bone marrow mesenchymal stem cell-derived exosomes promote plasminogen activator inhibitor 1 expression in vascular cells in the local microenvironment during rabbit osteonecrosis of the femoral head

Background

Nontraumatic osteonecrosis of the femoral head (NONFH) is a highly disabling orthopedic disease in young individuals. Plasminogen activator inhibitor 1 (PAI-1) has been reported to be positively associated with NONFH. We aimed to investigate the dysregulating PAI-1 in bone marrow mesenchymal stem cells (BMMSCs) and vascular cells in rabbit steroid-induced NONFH.

Methods

To verify the hypothesis that BMMSCs could promote thrombus formation in a paracrine manner, we collected exosomes from glucocorticoid-treated BMMSCs (GB-Exo) to determine their regulatory effects on vascular cells. microRNA sequencing was conducted to find potential regulators in GB-Exo. Utilizing gain-of-function and knockdown approaches, we testified the regulatory effect of microRNA in exosomes.

Results

The expression of PAI-1 was significantly increased in the local microenvironment of the femoral head in the ONFH model. GB-Exo promoted PAI-1 expression in vascular smooth muscle cells and vascular endothelial cells. We also revealed that miR-451-5p in GB-Exo plays a crucial role for the elevated PAI-1. Moreover, we identified miR-133b-3p and tested its role as a potential inhibitor of PAI-1.

Conclusions

This study provided considerable evidence for BMMSC exosomal miR-mediated upregulation of the fibrinolytic regulator PAI-1 in vascular cells. The disruption of coagulation and low fibrinolysis in the femoral head will eventually lead to a disturbance in the microcirculation of NONFH. We believe that our findings could be of great significance for guiding clinical trials in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-020-01991-2.

Long Non-Coding RNAs as Strategic Molecules to Augment the Radiation Therapy in Esophageal Squamous Cell Carcinoma

Intrinsic resistance to ionizing radiation is the major impediment in the treatment and clinical management of esophageal squamous cell carcinoma (ESCC), leading to tumor relapse and poor prognosis. Although several biological and molecular mechanisms are responsible for resistance to radiotherapy in ESCC, the molecule(s) involved in predicting radiotherapy response and prognosis are still lacking, thus requiring a detailed understanding. Recent studies have demonstrated an imperative correlation amongst several long non-coding RNAs and their involvement in complex cellular networks like DNA damage and repair, cell cycle, apoptosis, proliferation, and epithelial-mesenchymal transition. Additionally, accumulating evidence has suggested abnormal expression of lncRNAs in malignant tumor cells before and after radiotherapy effects in tumor cells' sensitivity. Thus, lncRNAs indeed represent unique molecules that can influence tumor cell susceptibility for various clinical interventions. On this note, herein, we have summarized the current status of lncRNAs in augmenting resistance/sensitivity in ESCC against radiotherapy. In addition, we have also discussed various strategies to increase the radiosensitivity in ESCC cells under clinical settings.

A three-mRNA status risk score has greater predictive ability compared with a lncRNA-based risk score for predicting prognosis in patients with hepatocellular carcinoma

Hepatocellular carcinoma (HCC) represents the fifth most common cause of cancer-associated mortality in men, and the seventh in women, worldwide. The aim of the present study was to identify a reliable and robust RNA-based risk score for the survival prediction of patients with hepatocellular carcinoma (HCC). Gene expression data from HCC and healthy control samples were obtained from The Cancer Genome Atlas to screen differentially expressed mRNAs and long non-coding RNAs (lncRNAs). Univariate and multivariate Cox proportional-hazards regression models and the LASSO algorithm for the Cox proportional-hazards model (LASSO Cox-PH model) were used to identify the prognostic mRNAs and lncRNAs among differentially expressed mRNAs (DEMs) and differentially expressed lncRNAs (DELs), respectively. Prognostic risk scores were generated based on the expression level or status of the prognostic lncRNAs and mRNAs, and the predictive abilities of these RNAs in TCGA and validation datasets were compared. Functional enrichment analyses were also performed. The results revealed a total of 154 downregulated and 625 upregulated mRNAs and 18 upregulated lncRNAs between tumor and control samples in TCGA dataset. A three-mRNA and a five-lncRNA expression signatures were identified using the LASSO Cox-PH model. Three-mRNA and five-lncRNA expression and status risk scores were generated. Using likelihood ratio P-values and area under the curve values from TCGA and the validation datasets, the three-mRNA status risk score was more accurate compared with the other risk scores in predicting the mortality of patients with HCC. The three identified mRNAs, including hepatitis A virus cellular receptor 1, MYCN proto-oncogene BHLH transcription factor and stratifin, were associated with the cell cycle and oocyte maturation pathways. Therefore, a three-mRNA status risk score may be valuable and robust for risk stratification of patients with HCC. The three-mRNA status risk score exhibited greater prognostic value compared with the lncRNA-based risk score.

Mesenchymal stem cells-derived and siRNAs-encapsulated exosomes inhibit osteonecrosis of the femoral head

Osteonecrosis of the femoral head (ONFH) is a progressive, obstinate and disabling disease. At present, the treatment of ONFH is still a global medical problem. We aim to explore the role of bone mesenchymal stem cells (BMSCs)‐derived and siRNAs‐encapsulated exosomes (siRNAs‐encapsulated BMSCexos) in ONFH. We first isolated BMSCexos and screened siRNAs of 6 ONFH‐related genes for siRNAs‐encapsulated BMSCexo. The expression of these 6 ONFH‐related genes in dexamethasone (DXM)‐treated MC3T3‐E1 cell, cell model of ONFH, was detected by RT‐qPCR and Western blot analysis. And then, we performed CCK‐8 assay, angiogenesis assay and HE staining analysis to test the promotion role of the siRNAs‐encapsulated BMSCexo for angiogenesis during ONFH repair. The results suggest that the obtained particles were BMSCexos. The screened effective siRNAs could effectively knock down their expression in VECs. Moreover, siRNAs‐encapsulated BMSCexo could effectively knock down the expression of these genes in VECs. In addition, siRNAs‐encapsulated BMSCexo promote angiogenesis during ONFH repair. In conclusion, we found siRNAs‐encapsulated BMSCexos could promote ONFH repair by angiogenesis, and indicated exosome as the new siRNA carrier is of great significance to improve the efficiency of RNAi.

Comparative Transcriptomic Analysis of the Pituitary Gland between Cattle Breeds Differing in Growth: Yunling Cattle and Leiqiong Cattle

The hypothalamic-pituitary-thyroid (HPT) axis hormones regulate the growth and development of ruminants, and the pituitary gland plays a decisive role in this process. In order to identify pivotal genes in the pituitary gland that could affect the growth of cattle by regulating the secretion of hormones, we detected the content of six HPT hormones related to growth in the plasma of two cattle breeds (Yunling and Leiqiong cattle, both also known as the zebu cattle) with great differences in growth and compared the transcriptome data of their pituitary glands. Our study found that the contents of GH, IGF, TSH, thyroxine, triiodothyronine, and insulin were significantly different between the two breeds, which was the main cause of the difference in growth; 175 genes were identified as differentially expressed genes (DEGs). Functional association analyses revealed that DEGs were mainly involved in the process of transcription and signal transduction. Combining the enrichment analysis and protein interaction analysis, eight DEGs were predicted to control the growth of cattle by affecting the expression of growth-related hormones in the pituitary gland. In summary, our results suggested that SLC38A1, SLC38A3, DGKH, GNB4, GNAQ, ESR1, NPY, and GAL are candidates in the pituitary gland for regulating the growth of Yunling and Leiqiong cattle by regulating the secretion of growth-related hormones. This study may help researchers further understand the growth mechanisms and improve the artificial selection of zebu cattle.

An optimal prognostic model based on gene expression for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most prevalent type of RCC; however, prognostic prediction tools for ccRCC are scant. Developing mRNA or long non-coding RNA (lncRNA)-based risk assessment tools may improve the prognosis in patients with ccRCC. RNA-sequencing and prognostic data from patients with ccRCC were downloaded from The Cancer Genome Atlas and the European Bioinformatics Institute Array database at the National Center for Biotechnology Information. Differentially expressed (DE) RNAs (DERs) and prognostic DERs were screened between less favorable and favorable prognoses using the limma package in R 3.4.1, and analyzed using univariate and multivariate Cox regression analyses, respectively. Risk score models were constructed using optimal combinations of DEmRNAs and DElncRNAs identified using the Least Absolute Shrinkage And Selection Operator Cox regression model of the penalized package. Associations between risk score models and overall survival time were evaluated. Independent prognostic clinical factors were screened using univariate and multivariate Cox regression analyses, and nomogram models were constructed. Gene Ontology biological processes and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were conducted using the clusterProfiler package in R3.4.1. A total of 451 DERs were identified, including 404 mRNAs and 47 lncRNAs, between less favorable and favorable prognoses, and 269 DERs, including 233 mRNAs and 36 lncRNAs, were identified as independent prognostic factors. Optimal combinations including 10 DEmRNAs or 10 DElncRNAs were screened using four risk score models based on the status or expression levels of the 10 DEmRNAs or 10 DElncRNAs. The model based on the expression levels of the 10 DEmRNAs had the highest prognostic power. These prognostic DEmRNAs may be involved in biological processes associated with the inflammatory response, complement and coagulation cascades and neuroactive ligand-receptor interaction pathways. The present validated risk assessment tool based on the expression levels of these 10 DEmRNAs may help to identify patients with ccRCC at a high risk of mortality. These 10 DEmRNAs in optimal combinations may serve as prognostic biomarkers and help to elucidate the pathogenesis of ccRCC.

An eight-mRNA signature outperforms the lncRNA-based signature in predicting prognosis of patients with glioblastoma

BACKGROUND

The prognosis of the glioblastoma (GBM) is dismal. This study aims to select an optimal RNA signature for prognostic prediction of GBM patients.

METHODS

For the training set, the long non-coding RNA (lncRNA) and mRNA expression profiles of 151 patients were downloaded from the TCGA. Differentially expressed mRNAs (DEGs) and lncRNAs (DE-lncRNAs) were identified between good prognosis and bad prognosis patients. Optimal prognostic mRNAs and lncRNAs were selected respectively, by using univariate Cox proportional-hazards (PH) regression model and LASSO Cox-PH model. Subsequently, four prognostic scoring models were built based on expression levels or expression status of the selected prognostic lncRNAs or mRNAs, separately. Each prognostic model was applied to the training set and an independent validation set. Function analysis was used to uncover the biological roles of these prognostic DEGs between different risk groups classified by the mRNA-based signature.

RESULTS

We obtained 261 DEGs and 33 DE-lncRNAs between good prognosis and bad prognosis patients. A panel of eight mRNAs and a combination of ten lncRNAs were determined as predictive RNAs by LASSO Cox-PH model. Among the four prognostic scoring models using the eight-mRNA signature or the ten-lncRNA signature, the one based on the expression levels of the eight mRNAs showed the greatest predictive power. The DEGs between different risk groups using the eight prognostic mRNAs were functionally involved in calcium signaling pathway, neuroactive ligand-receptor interaction pathway, and Wnt signaling pathway.

CONCLUSION

The eight-mRNA signature has greater prognostic value than the ten-lncRNA-based signature for GBM patients based on bioinformatics analysis.

Identification of potential miRNA biomarkers for traumatic osteonecrosis of femoral head

Traumatic osteonecrosis of femoral head (TONFH) is a common orthopedic disease caused by physical injury in hip. However, the unclear pathogenesis mechanism of TONFH and lacking of simple noninvasive early diagnosis method cause the necessity of hip replacement for most patients with TONFH. In this study, we aimed to identify circulating microRNAs (miRNAs) by integrated bioinformatics analyses as potential biomarker of TONFH. mRNA expression profiles were downloaded from the Gene Expression Omnibus database. Then we combined two miRNA screen methods: Weighted gene co-expression network analysis and fold change based differentially expressed miRNAs analysis. As a result, we identified 14 key miRNAs as potential biomarkers for TONFH. Besides, 302 target genes of these miRNAs were obtained and the miRNA-mRNA interaction network was constructed. Furthermore, the results of Kyoto Encyclopedia of Gene and Genome pathway analysis, Gene Ontology function analysis, protein-protein interaction (PPI) network analysis and PPI network module analysis showed close correlation between these 14 key miRNAs and TONFH. Then we established receiver operating characteristic curves and identified 6-miRNA signature with highly diagnosis value including miR-93-5p (area under the curve [AUC] = 0.93), miR-1324 (AUC = 0.92), miR-4666a-3p (AUC = 0.92), miR-5011-3p (AUC = 0.92), and miR-320a (AUC = 0.89), miR-185-5p (AUC = 0.89). Finally, the results of quantitative real-time polymerase chain reaction confirmed the significantly higher expression of miR-93-5p and miR-320a in the serum of patients with ONFH. These circulating miRNAs could serve as candidate early diagnosis markers and potential treatment targets of TONFH.

BMSC-derived exosomes carrying microRNA-122-5p promote proliferation of osteoblasts in osteonecrosis of the femoral head

Mesenchymal stem cells (MSCs) with multipotential differentiation capacity can differentiate into bone cells under specific conditions and can be used to treat osteonecrosis (ON) of the femoral head (ONFH) through cell transplantation. The current study aims to explore the role of bone marrow (BM) MSCs (BMSCs)-derived exosomes carrying microRNA-122-5p (miR-122-5p) in ONFH rabbit models.First, rabbit models with ONFH were established. ONFH-related miRNAs were screened using the Gene Expression Omnibus (GEO) database. A gain-of-function study was performed to investigate the effect of miR-122-5p on osteoblasts and BMSCs and effects of exosomes carrying miR-122-5p on ONFH. Co-culture experiments for osteoblasts and BMSCs were performed to examine the role of exosomal miR-122-5p in osteoblast proliferation and osteogenesis. The target relationship between miR-122-5p and Sprouty2 (SPRY2) was tested.MiR-122, significantly decreased in ONFH in the GSE89587 expression profile, was screened. MiR-122-5p negatively regulated SPRY2 and elevated the activity of receptor tyrosine kinase (RTK), thereby promoting the proliferation and differentiation of osteoblasts. In vivo experiments indicated that bone mineral density (BMD), trabecular bone volume (TBV), and mean trabecular plate thickness (MTPT) of femoral head were increased after over-expressing miR-122-5p in exosomes. Significant healing of necrotic femoral head was also observed.Exosomes carrying over-expressed miR-122-5p attenuated ONFH development by down-regulating SPRY2 via the RTK/Ras/mitogen-activated protein kinase (MAPK) signaling pathway. Findings in the present study may provide miR-122-5p as a novel biomarker for ONFH treatment.

Identification of differentially expressed genes in hip cartilage with femoral head necrosis, based on genome‑wide expression profiles

Necrosis of the femoral head (NFH), a severe orthopedic disease in adults, involves the collapse of the femoral head. The pathophysiological mechanisms underlying NFH are yet to be fully investigated. The aim of the present study was to identify potentially important genes and signaling pathways involved in NFH and investigate their molecular mechanisms. Gene expression profiles of patients with NFH and healthy controls were compared using the Gene Expression Omnibus (GEO) database repository of the National Center of Biotechnology Information. GSE74089 from the GEO database included 4 patients with NFH and 4 healthy individuals. A total of 1,191 differentially expressed genes (DEGs) were identified between the patients with NFH and controls, including 743 upregulated and 448 downregulated DEGs. Then, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that upregulated DEGs were mainly involved in the phosphoinositide 3-kinase/protein kinase B signaling pathway, focal adhesion and extracellular matrix-receptor interactions. Additionally, protein-protein interaction (PPI) analysis identified the most central DEGs as vascular endothelial growth factor A, Jun proto-oncogene, cyclin D1, fibroblast growth factor 2, HECT domain and ankyrin repeat-containing E3 ubiquitin protein ligase 1, protein kinase Cα, bone morphogenetic protein 2 and prostaglandin-endoperoxide synthase 2. PPI analysis also identified guanine nucleotide-binding protein, γ13 as the most commonly downregulated gene based on different centrality. The results of the present study may provide novel insight into the genes and associated pathways involved in NFH, and aid the identification of novel therapeutic targets and biomarkers in the treatment of NFH.

Metal trabecular bone reconstruction system better improves clinical efficacy and biomechanical repair of osteonecrosis of the femoral head than free vascularized fibular graft: A case-control study

In this study, we aim to compare and analyze the biomechanical repair and clinical efficacy of osteonecrosis of the femoral head (ONFH) with the use of metal trabecular bone reconstruction system and free vascularized fibular graft. The study enrolled 66 adult patients from medical records of nontraumatic ARCO 2A-3B stage ONFH. A simple ONFH model without surgical treatment was established in 13 cases, 29 cases were treated with metal trabecular bone reconstruction system, and 24 cases were treated with free vascularized fibular graft. Computer-recognized and extracted femur outlines were imported, and three-dimensional reconstructions were performed. The stress concentration and stress peak value were analyzed, and the Harris score, visual analog scale pain score, and operation status of the above patients were compared. Finally, quality of life assessment was performed using SF-36 scale. Metal trabecular bone reconstruction system provided less operation time, blood loss, and the total length of postoperative hospital stay than free vascularized fibular graft. Metal trabecular bone reconstruction system promoted bone reconstruction, increased bone mineral density and Harris score. The total clinical effective rate of young patients (20-40 years) was higher than that of older patients (41-60 years). Metal trabecular bone reconstruction system provided higher physical component summary, mental component summary, and role/social component summary than free vascularized fibular graft. This study demonstrates that both metal trabecular bone reconstruction system and free vascularized fibular graft can prevent or delay the progression of ONFH, while metal trabecular bone reconstruction system is a better choice because of better short-term clinical efficacy.

Long noncoding RNAs: a missing link in osteoporosis

Osteoporosis is a systemic disease that results in loss of bone density and increased fracture risk, particularly in the vertebrae and the hip. This condition and associated morbidity and mortality increase with population ageing. Long noncoding (lnc) RNAs are transcripts longer than 200 nucleotides that are not translated into proteins, but play important regulatory roles in transcriptional and post-transcriptional regulation. Their contribution to disease onset and development is increasingly recognized. Herein, we present an integrative revision on the studies that implicate lncRNAs in osteoporosis and that support their potential use as therapeutic tools. Firstly, current evidence on lncRNAs involvement in cellular and molecular mechanisms linked to osteoporosis and its major complication, fragility fractures, is reviewed. We analyze evidence of their roles in osteogenesis, osteoclastogenesis, and bone fracture healing events from human and animal model studies. Secondly, the potential of lncRNAs alterations at genetic and transcriptomic level are discussed as osteoporosis risk factors and as new circulating biomarkers for diagnosis. Finally, we conclude debating the possibilities, persisting difficulties, and future prospects of using lncRNAs in the treatment of osteoporosis.

Long Noncoding RNA FAM201A Mediates the Radiosensitivity of Esophageal Squamous Cell Cancer by Regulating ATM and mTOR Expression via miR-101

Background: The aim of the present study was to identify the potential long non-coding (lnc.)-RNA and its associated molecular mechanisms involved in the regulation of the radiosensitivity of esophageal squamous cell cancer (ESCC) in order to assess whether it could be a biomarker for the prediction of the response to radiotherapy and prognosis in patients with ESCC.

Methods: Microarrays and bioinformatics analysis were utilized to screen the potential lncRNAs associated with radiosensitivity in radiosensitive (n = 3) and radioresistant (n = 3) ESCC tumor tissues. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was performed in 35 ESCC tumor tissues (20 radiosensitive and 15 radioresistant tissues, respectively) to validate the lncRNA that contributed the most to the radiosensitivity of ESCC (named the candidate lncRNA). MTT, flow cytometry, and western blot assays were conducted to assess the effect of the candidate lncRNA on radiosensitivity in vitro in ECA109/ECA109R ESCC cells. A mouse xenograft model was established to confirm the function of the candidate lncRNA in the radiosensitivity of ESCC in vivo. The putative downstream target genes regulated by the candidate lncRNA were predicted using Starbase 2.0 software and the TargetScan database. The interactions between the candidate lncRNA and the putative downstream target genes were examined by Luciferase reporter assay, and were confirmed by PCR.

Results: A total of 113 aberrantly expressed lncRNAs were identified by microarray analysis, of which family with sequence similarity 201-member A (FAM201A) was identified as the lncRNA that contributed the most to the radiosensitivity of ESCC. FAM201A was upregulated in radioresistant ESCC tumor tissues and had a poorer short-term response to radiotherapy resulting in inferior overall survival. FAM201A knockdown enhanced the radiosensitivity of ECA109/ECA109R cells by upregulating ataxia telangiectasia mutated (ATM) and mammalian target of rapamycin (mTOR) expression via the negative regulation of miR-101 expression. The mouse xenograft model demonstrated that FAM201A knockdown improved the radiosensitivity of ESCC.

Conclusion: The lncRNA FAM201A, which mediated the radiosensitivity of ESCC by regulating ATM and mTOR expression via miR-101 in the present study, may be a potential biomarker for predicting radiosensitivity and patient prognosis, and may be a therapeutic target for enhancing cancer radiosensitivity in ESCC.