Construction of a ceRNA network of regulated ferroptosis in doxorubicin-induced myocardial injury.

Ferroptosis and long-noncoding RNAs (lncRNAs) play crucial roles in doxorubicin (DOX)-induced myocardial injury (DIMI). Nevertheless, there is no research to construct competing endogenous RNAs (ceRNAs) network between lncRNAs and ferroptosis-related key gene. So our research was designed to screen ferroptosis-related genes from differentially expressed mRNAs in DIMI and construct lncRNAs regulated ferroptosis-related key gene ceRNAs network. The male mice were injected with DOX intraperitoneally to induce myocardial injury, myocardial injury was evaluated by hematoxylin and eosin (HE) staining, and ferroptosis-related protein-glutathione peroxidase 4 (GPx4) protein expression was detected. The differentially expressed lncRNAs and mRNAs were detected by microarray, and the ferroptosis-related genes were screened to construct a protein-protein associations (PPA) network, the highest maximal clique centrality (MCC) score gene were identified by Cytoscape software, miRNAs bound to key genes and lncRNAs bound to miRNAs were predicted; then, the obtained lncRNAs were intersected with differentially expressed lncRNAs detected by microarray. Finally, the lncRNA/miRNA/mRNA ceRNA network of the highest MCC score gene regulating ferroptosis in DIMI was constructed. The expressions of the key components in ceRNA network were detected by qRT-PCR. Compared with the control group, in the DOX group, myocardial enzymes and HE staining showed that myocardium structure was changed, and GPx4 protein expression was decreased. The differentially expressed 10,265 lncRNAs and 6,610 mRNAs in the DOX group were detected via microarray. Among them, 114 ferroptosis-related genes were obtained to construct PPA networks, and Becn1 was identified as the key gene. Finally, the ceRNA network including Becn1, three miRNAs and four lncRNAs was constructed by predicting data of the Starbase database. The relative expressions of these components in ceRNA net were up-regulated and consistent with microarray results. Based on the microarray detection results and bioinformatics analysis, we screened ferroptosis-related gene Becn1 and constructed the lncRNA/miRNA/mRNA ceRNA network of regulated ferroptosis in DIMI.

Ye H ，Li Y ，Li L ，Huang Y ，Wang J ，Gao Q ... -  《PeerJ》

Excavating novel diagnostic and prognostic long non-coding RNAs (lncRNAs) for head and neck squamous cell carcinoma: an integrated bioinformatics analysis of competing endogenous RNAs (ceRNAs) and gene co-expression networks.

Yang L ，Lu P ，Yang X ，Li K ，Chen X ，Qu S ... -  《-》

Bioinformatics identification and integrative analysis of ferroptosis-related key lncRNAs in patients with osteoarthritis.

Ferroptosis and dysregulation of long non-coding RNA (lncRNA) have been described to be strictly relevant to the pathogenesis of osteoarthritis (OA). However, the connection between ferroptosis and lncRNA in OA is poorly appreciated. Herein, we investigated the functional contribution of lncRNA markers correlated with the progression of human OA by comprehensive bioinformatics analysis of a panoramic network of competing endogenous RNA (ceRNA) based on ferroptosis-related genes (FRGs). FRGs-related competing endogenous RNA (ceRNA) networks were generated using differentially expressed genes based on OA-related whole transcriptome data from the Gene Expression Omnibus (GEO) database via starBase, miRTarBase, and miRWalk databases. The pivotal lncRNAs were ascertained by topological features (degree, betweenness, and closeness) and subceRNA networks were re-visualized. The expression difference of pivotal lncRNAs was verified by quantitative real-time polymerase chain reaction (qRT-PCR). The latent molecular mechanisms of the global ceRNA and subceRNA networks were uncovered by the R package clusterProfiler-based enrichment analysis. A total of 98 dysregulated lncRNA-miRNA-mRNA regulatory relationships were attained in the FRGs-related panoramic ceRNA network of OA, covering 26 mRNAs, 20 miRNAs, and 20 lncRNAs. Three lncRNAs (AC011511.5, AL358072.1, and C9orf139) were ascertained as the central lncRNAs in the panoramic ceRNA network. Functional ensemble analysis illustrated that both the panoramic ceRNA network and the subceRNA network were integrally affiliated with the immune-inflammatory response, oxygen homeostasis, and cell death (apoptosis, autophagy, and ferroptosis). Comprehensive bioinformatics analysis of the FRGs-related ceRNA network determined three molecular biomarkers of lncRNAs that might be affiliated with OA progression.

Yang T ，Yang G ，Wang G ，Jia D ，Xiong B ，Lu X ，Li Y ... -  《-》

Identification and characterization of the ferroptosis-related ceRNA network in irreversible pulpitis.

The role of ferroptosis in irreversible pulpitis (IP) remains unclear. The competing endogenous RNA (ceRNA) theory that has been widely investigated is rarely used studied in IP. Hub lncRNAs selected from a ceRNA network may provide a novel hypothesis for the interaction of ferroptosis and IP. Differentially expressed genes (DEGs) were intersected with 484 ferroptosis markers to identify differentially expressed ferroptosis-related genes (DE-FRGs). Functional analysis and protein-protein interaction (PPI) networks were constructed to reveal the functions of DE-FRGs. Then, coexpression analyses were conducted between DE-FRGs and DElncRNAs to define ferroptosis-related DElncRNAs (FR-DElncRNAs). Predictions of DE-FRG- and FR-DElncRNA-related miRNAs were obtained, and members of both groups were selected. Additionally, two ceRNA networks consisting of FR-DElncRNAs, miRNAs and DE-FRGs from upregulated and downregulated groups were built. Finally, the hub lncRNAs of the ceRNA networks were used for immuno-infiltration analysis and qPCR verification. According to the results of PCA and clustering analysis, 5 inflamed and 5 healthy pulp tissue samples were selected for analysis. The intersection of DEGs with 484 ferroptosis marker genes identified 72 DE-FRGs. The response to stimulus, cellular process, signaling, localization, and biological regulation pathways related to DE-FRGs were enriched. In total, 161 downregulated and 40 upregulated FR-DElncRNAs were chosen by coexpression analysis for further investigation. The MultimiR package and starBase were used to predict miRNAs of DE-FRGs and FR-DElncRNAs, respectively. The upregulated ceRNA network contained 2 FR-DElncRNAs (↑), 19 miRNAs (↓) and 22 DE-FRGs (↑). The downregulated network contained 44 FR-DElncRNAs (↓), 251 miRNAs (↑) and 10 DE-FRGs (↓). Six hub lncRNAs were identified based on the MCC method (LUCAT1 and AC106897.1 ↑; LINC00943, AL583810.1, AC068888.1, and AC125257.1↓). In addition, strong relationships between hub lncRNAs and immune cells were shown by immune infiltration analysis. Finally, validated by qPCR assays of the pulp tissue of IP patients, the expression levels in clinical samples were consistent with the microarray data. Two ceRNA networks were comprehensively constructed, and 6 hub lncRNAs were identified. These genes provide novel insights into the relationship between ferroptosis and IP. Intriguingly, the LINC00943/hsa-miR-29a-3p/PDK4 axis was deemed to be the key node in this network.

Xie Q ，Yu H ，Liu Z ，Zhou B ，Fang F ，Qiu W ，Wu H ... -  《Frontiers in Immunology》

Comprehensive analysis of lncRNA-miRNA-mRNA networks during osteogenic differentiation of bone marrow mesenchymal stem cells.

Long non-coding RNA (lncRNA) plays crucial role in osteogenic differentiation of bone marrow mesenchymal stem cells (BMMSCs), involving in regulation of competing endogenous RNA (ceRNA) mechanisms and conduction of signaling pathways. However, its mechanisms are poorly understood. This study aimed to investigate lncRNAs, miRNAs and mRNAs expression profiles in rat BMMSCs (rBMMSCs) osteogenic differentiation, screen the potential key lncRNA-miRNA-mRNA networks, explore the putative functions and identify the key molecules, as the basis of studying potential mechanism of rBMMSCs osteogenic differentiation driven by lncRNA, providing molecular targets for the management of bone defect. High-throughput RNA sequencing (RNA-seq) was used to determine lncRNAs, miRNAs, and mRNAs expression profiles at 14-day rBMMSCs osteogenesis. The pivotal lncRNA-miRNA and miRNA-mRNA networks were predicted from sequencing data and bioinformatic analysis, and the results were exported by Cytoscape 3.9.0 software. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used for functional exploration. Real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was performed to validate lncRNAs, miRNAs and mRNAs. rBMMSCs were identified, and the osteogenic and adipogenic differentiation ability were detected. A total of 8634 lncRNAs were detected by RNA-seq, and 1524 differential expressed lncRNAs, of which 812 up-regulated and 712 down-regulated in osteo-inductive groups compared with control groups. 30 up-regulated and 61 down-regulated miRNAs, 91 miRNAs were differentially expressed in total. 2453 differentially expressed mRNAs including 1272 up-expressed and 1181 down-expressed were detected. 10 up-regulated lncRNAs were chosen to predict 21 down-regulated miRNAs and 650 up-regulated mRNAs. 49 lncRNA-miRNA and 1515 miRNA-mRNA interactive networks were constructed. GO analysis showed the most important enrichment in cell component and molecular function were "cytoplasm" and "protein binding", respectively. Biological process related to osteogenic differentiation such as "cell proliferation", "wound healing", "cell migration", "osteoblast differentiation", "extracellular matrix organization" and "response to hypoxia" were enriched. KEGG analysis showed differentially expressed genes were mainly enriched in "PI3K-Akt signaling pathway", "Signaling pathway regulating pluripotency of stem cells", "cGMP-PKG signaling pathway", "Axon guidance" and "Calcium signaling pathway". qRT-PCR verified that lncRNA Tug1, lncRNA AABR07011996.1, rno-miR-93-5p, rno-miR-322-5p, Sgk1 and Fzd4 were consistent with the sequencing results, and 4 lncRNA-miRNA-mRNA networks based on validations were constructed, and enrichment pathways were closely related to "PI3K-Akt signaling pathway", "Signaling pathway regulating pluripotency of stem cells" and "Wnt signaling pathway". lncRNAs, miRNAs and mRNAs expression profiles provide clues for future studies on their roles for BMMSCs osteogenic differentiation. Furthermore, lncRNA-miRNA-mRNA networks give more information on potential new mechanisms and targets for management on bone defect.

Liu J ，Yao Y ，Huang J ，Sun H ，Pu Y ，Tian M ，Zheng M ，He H ，Li Z ... -  《BMC GENOMICS》