Identification of novel transcription factor-microRNA-mRNA co-regulatory networks in pulmonary large-cell neuroendocrine carcinoma

Application of Machine Learning in Predicting Hepatic Metastasis or Primary Site in Gastroenteropancreatic Neuroendocrine Tumors

Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) account for 80% of gastroenteropancreatic neuroendocrine neoplasms (GEP-NENs). GEP-NETs are well-differentiated tumors, highly heterogeneous in biology and origin, and are often diagnosed at the metastatic stage. Diagnosis is commonly through clinical symptoms, histopathology, and PET-CT imaging, while molecular markers for metastasis and the primary site are unknown. Here, we report the identification of multi-gene signatures for hepatic metastasis and primary sites through analyses on RNA-SEQ datasets of pancreatic and small intestinal NETs tissue samples. Relevant gene features, identified from the normalized RNA-SEQ data using the mRMRe algorithm, were used to develop seven Machine Learning models (LDA, RF, CART, k-NN, SVM, XGBOOST, GBM). Two multi-gene random forest (RF) models classified primary and metastatic samples with 100% accuracy in training and test cohorts and >90% accuracy in an independent validation cohort. Similarly, three multi-gene RF models identified the pancreas or small intestine as the primary site with 100% accuracy in training and test cohorts, and >95% accuracy in an independent cohort. Multi-label models for concurrent prediction of hepatic metastasis and primary site returned >98.42% and >87.42% accuracies on training and test cohorts, respectively. A robust molecular signature to predict liver metastasis or the primary site for GEP-NETs is reported for the first time and could complement the clinical management of GEP-NETs.

Identification of potential biomarkers for idiopathic pulmonary fibrosis and validation of TDO2 as a potential therapeutic target

BACKGROUND

Idiopathic pulmonary fibrosis (IPF) is a progressive interstitial lung disease with a high mortality rate. On this basis, exploring potential therapeutic targets to meet the unmet needs of IPF patients is important.

AIM

To explore novel hub genes for IPF therapy.

METHODS

Here, we used public datasets to identify differentially expressed genes between IPF patients and healthy donors. Potential targets were considered based on multiple bioinformatics analyses, especially the correlation between hub genes and carbon monoxide diffusing capacity of carbon monoxide, forced vital capacity, and patient survival rate. The mRNA levels of the hub genes were determined through quantitative real-time polymerase chain reaction.

RESULTS

We found that TDO2 was upregulated in IPF patients and predicted poor prognosis. Surprisingly, single-cell RNA sequencing data analysis revealed significant enrichment of TDO2 in alveolar fibroblasts, indicating that TDO2 may participate in the regulation of proliferation and survival. Therefore, we verified the upregulated expression of TDO2 in an experimental mouse model of transforming growth factor-β (TGF-β)-induced pulmonary fibrosis. Furthermore, the results showed that a TDO2 inhibitor effectively suppressed TGF-β-induced fibroblast activation. These findings suggest that TDO2 may be a potential target for IPF treatment. Based on transcription factors-microRNA prediction and scRNA-seq analysis, elevated TDO2 promoted the IPF proliferation of fibroblasts and may be involved in the P53 pathway and aggravate ageing and persistent pulmonary fibrosis.

CONCLUSION

We provided new target genes prediction and proposed blocking TGF-β production as a potential treatment for IPF.

Identification of potential microRNA diagnostic panels and uncovering regulatory mechanisms in breast cancer pathogenesis

Early diagnosis of breast cancer (BC), as the most common cancer among women, increases the survival rate and effectiveness of treatment. MicroRNAs (miRNAs) control various cell behaviors, and their dysregulation is widely involved in pathophysiological processes such as BC development and progress. In this study, we aimed to identify potential miRNA biomarkers for early diagnosis of BC. We also proposed a consensus-based strategy to analyze the miRNA expression data to gain a deeper insight into the regulatory roles of miRNAs in BC initiation. Two microarray datasets (GSE106817 and GSE113486) were analyzed to explore the differentially expressed miRNAs (DEMs) in serum of BC patients and healthy controls. Utilizing multiple bioinformatics tools, six serum-based miRNA biomarkers (miR-92a-3p, miR-23b-3p, miR-191-5p, miR-141-3p, miR-590-5p and miR-190a-5p) were identified for BC diagnosis. We applied our consensus and integration approach to construct a comprehensive BC-specific miRNA-TF co-regulatory network. Using different combination of these miRNA biomarkers, two novel diagnostic models, consisting of miR-92a-3p, miR-23b-3p, miR-191-5p (model 1) and miR-92a-3p, miR-23b-3p, miR-141-3p, and miR-590-5p (model 2), were obtained from bioinformatics analysis. Validation analysis was carried out for the considered models on two microarray datasets (GSE73002 and GSE41922). The model based on similar network topology features, comprising miR-92a-3p, miR-23b-3p and miR-191-5p was the most promising model in the diagnosis of BC patients from healthy controls with 0.89 sensitivity, 0.96 specificity and area under the curve (AUC) of 0.98. These findings elucidate the regulatory mechanisms underlying BC and represent novel biomarkers for early BC diagnosis.

KLF9 regulates miR-338-3p/NRCAM axis to block the progression of osteosarcoma cells

Background: MiR-338-3p is revealed to serve as a tumor suppressor in several carcinomas. Whereas, the effect of miR-338-3p in the progression of osteosarcoma has not been explored. The aim of this paper was to analyze the functional influences of miR-338-3p on osteosarcoma progression and the potential mechanism.

Methods: The expression of genes and miRNAs in osteosarcoma cells was assessed via western blotting or quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Osteosarcoma cellular proliferation was explored by MTT and EdU incorporation assay. Osteosarcoma cellular migratory and invasive capacity was explored by wound-healing and transwell assay. Bioinformatics approaches were adopted to predict target genes. The relationships between miR-338-3p and neuron‑glial‑related cell adhesion (NRCAM), between kruppel-like factor 9 (KLF9) and miR-338-3p were verified by dual-luciferase reporter assay.

Results: We found that miR-338-3p was reduced in osteosarcoma and that higher expression of miR-338-3p suppressed proliferative, invasive and migratory ability of osteosarcoma cells. Furthermore, the result showed that overexpression of NRCAM could reduce the anti-tumor role of miR-338-3p in osteosarcoma cells. In addition, we found that overexpression of KLF9 could enhance the expression level of miR-338-3p in osteosarcoma cells.

Conclusion: The KLF9/miR-338-3p/NRCAM axis played a significant role in regulating osteosarcoma progression, which may become a promising therapeutic method for osteosarcoma.

Identifying lncRNA- and Transcription Factor-Associated Regulatory Networks in the Cortex of Rats With Deep Hypothermic Circulatory Arrest

Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) are involved in the mechanism underlying cerebral dysfunction after deep hypothermic circulatory arrest (DHCA), although the exact details have not been elucidated. To explore the expression profiles of lncRNAs and miRNAs in DHCA cerebral injury, we determined the lncRNA, miRNA and mRNA expression profiles in the cerebral cortex of DHCA and sham rats. First, a rat model of DHCA was established, and high-throughput sequencing was performed to analyze the differentially expressed RNAs (DERNAs). Then, the principal functions of the significantly deregulated genes were identified using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Expression networks (lncRNAs-miRNAs-mRNAs and transcription factors (TFs)-miRNAs-mRNAs) were also established. Finally, the expression of DERNAs was confirmed by quantitative real-time PCR (RT-qPCR). We identified 89 lncRNAs, 45 miRNAs and 59 mRNAs between the DHCA and sham groups and constructed a comprehensive competitive endogenous RNAs (ceRNAs) network. A TF-miRNA-mRNA regulatory network was also established. Finally, we predicted that Lcorl-miR-200a-3p-Ttr, BRD4-Ccl2 and Ep300-miR-200b-3p-Tmem72 may participate in the pathogenesis of DHCA cerebral injury.