Search

Abstract Colorectal cancer (CRC) is a wide-spread gastrointestinal cancer that is associated with augmented morbidity and mortality, and we do not yet have a deep understanding of its epidemiology and carcinogenicity. The transcriptome can reveal the complexity and heterogeneity of tumors and uncover new biomarkers or treatment options. In this study, we identified messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), round RNAs (circRNAs), and microRNAs (miRNAs) using whole-transcriptome sequencing and generated competing endogenous RNA (ceRNA) modulatory axes. We conducted whole transcriptome sequencing on 10 CRC and para-cancer (CRCP) samples and discovered 2465 differentially expressed (DE) mRNAs (DEmRNAs), 77 DE miRNAs (DEmiRNAs). 2852 DE lncRNAs (DElncRNAs) and 1477 DE circRNAs (DEcircRNAs). In addition, utilizing co-DE analysis, we generated the ceRNA axis. Subsequently, we employed the ceRNA axis to identify essential genes and corresponding associations with lncRNAs, circRNAs, and miRNAs in CRC. ceRNA regulatory network including mRNA-miRNA-lncRNA and mRNA-miRNA-circRNA. These modulatory axes potentially modulate the positive regulation of smooth muscle contraction, melanosome, plasma membrane, integral plasma membrane component and so on. Finally, the results of RNA sequencing (RNA-SEQ) were combined with the TCGA and GEO databases, and the DEGs strongly correlated with the TCGA-COAD overall survival (OS) as estimated by univariate cox and logarithmic rank analyses were cross-analyzed, and the co-upregulated DEGs were screened. Among the many DEs, KPNA2 was chosen for additional analysis. Using invitro experimentations, western blot, CCK8, EdU and other experiments were performed to verify the results. We found siRNA-based KPNA2 depletion reduces bladder cancer cells’ viability, migratory, and proliferative activities, which showed that the DEmRNA profiles were comparable to the sequencing information, confirming that the sequencing data were very reliable. These evidences highlight the ceRNA regulatory mechanisms in CRC and will aid future research into the molecular mechanisms behind colorectal cancer prevention and treatment.

Abstract Circular RNA (CircRNA)–microRNA (miRNA) interaction (CMI) is an important model for the regulation of biological processes by non-coding RNA (ncRNA), which provides a new perspective for the study of human complex diseases. However, the existing CMI prediction models mainly rely on the nearest neighbor structure in the biological network, ignoring the molecular network topology, so it is difficult to improve the prediction performance. In this paper, we proposed a new CMI prediction method, BEROLECMI, which uses molecular sequence attributes, molecular self-similarity, and biological network topology to define the specific role feature representation for molecules to infer the new CMI. BEROLECMI effectively makes up for the lack of network topology in the CMI prediction model and achieves the highest prediction performance in three commonly used data sets. In the case study, 14 of the 15 pairs of unknown CMIs were correctly predicted.

Abstract Purpose Multiple myeloma (MM) is the second most common hematologic malignancy, and there is no cure for this disease. This study aimed to explore the prognostic value of long noncoding RNAs (lncRNAs) in MM and to reveal related immune and chemotherapy resistance mechanisms. Methods In this study, lncRNA profiles from the Multiple Myeloma Research Foundation (MMRF) and Gene Expression Omnibus (GEO) databases were analyzed to identify lncRNAs linked to MM patient survival. A risk assessment model stratified patients into high- and low-risk groups, and survival was evaluated. Additionally, a triple-ceRNA (lncRNA–miRNA–mRNA) network was constructed, and functional analysis was performed. The research also involved immune function analysis and chemotherapy drug sensitivity assessment using oncoPredict and the GDSC dataset. Results We identified 422 lncRNAs significantly associated with overall survival in MM patients and ultimately focused on the 6 with the highest prognostic value. These lncRNAs were used to develop a risk score formula that stratified patients into high- and low-risk groups. Kaplan–Meier analysis revealed shorter survival in high-risk patients. We integrated this lncRNA signature with clinical parameters to construct a nomogram for predicting MM prognosis. Additionally, a triple-ceRNA network was constructed to reveal potential miRNA targets, coding genes related to these lncRNAs and significantly enriched pathways. Immune checkpoint gene expression and immune cell composition were also analyzed in relation to the lncRNA risk score. Finally, using the oncoPredict tool, we observed that high-risk patients exhibited decreased sensitivity to key MM chemotherapeutics, suggesting that lncRNA profiles are linked to chemotherapy resistance.

This manuscript summarizes information on the diverse range of RNA molecules and their role as competing endogenous RNAs (ceRNAs). Moreover, it provides an overview of ceRNA regulatory networks and their applications in ruminant biology. Knowledge of co-expression networks has increased with microarrays, RNA-seq, and scRNA-seq characterizing molecular mediators across various biological scales, using sequences from numerous blood and tissue samples. By synthesizing existing knowledge, this study summarizes interactions between coding and non-coding RNAs through microRNA response elements (MREs), elucidating large-scale regulatory networks throughout the transcriptome that influence the expression and activities of various ceRNAs. Identification of non-coding RNAs with important regulatory functions will revolutionize understanding of RNA biology, shifting from an mRNA-centric model to a complex network of RNA crosstalk. The ceRNA networks offer a more comprehensive and arguably more realistic perspective compared to protein–protein interaction (PPI) networks and weighted gene co-expression networks (WGCN). These ceRNA regulatory networks can describe potential molecular regulatory mechanisms related to functional and economically important traits in ruminants, plus contribute to disease and pathology research, by elucidating pathogenesis and potential drug effects in disease and cancer models. Furthermore, they can provide insights into farm animal biology, e.g., reproductive traits in goats and sheep, regulation of fat metabolism in beef cattle, heat stress responses, and lactation regulation in dairy cattle, fertility and muscle characteristics in buffalo, and resistance to high-salt and water-deprivation conditions in camels. In conclusion, ceRNA and associated regulatory networks should promote a new understanding of molecular mechanisms and identify candidate genes and metabolic-signaling pathways in ruminants.

Abstract Background Lipid metabolism plays a pivotal role in asthma pathogenesis. However, a comprehensive analysis of the importance of lipid metabolism-related genes (LMRGs) in regulating the immune microenvironment in asthma remains lacking. The transcriptome matrix was downloaded from the Gene Expression Omnibus (GEO) dataset. Differentially expressed analysis and weighted gene coexpression network analysis (WGCNA) were conducted on the GSE74986 dataset to select hub LMRGs, and gene set enrichment analysis (GSEA) was conducted to explore their biological functions. The CIBERSORT algorithm was used to determine immune infiltration in the asthma and control groups, and the correlation of diagnostic biomarkers and immune cells was performed via Spearman correlation analysis. Subsequently, a competitive endogenous RNA (ceRNA) network was constructed to investigate the hidden molecular mechanism of asthma. The expression levels of the hub genes were further validated in the GSE143192 dataset, and RT‒qPCR and immunofluorescence were performed to verify the reliability of the results in the OVA asthma model. Lastly, the ceRNA network was confirmed by qRT-PCR and RNAi experiments in the characteristic cytokine (IL-13)-induced asthma cellular model. Results ASAH1, ACER3 and SGPP1 were identified as hub LMRGs and were mainly involved in protein secretion, mTORC1 signaling, and fatty acid metabolism. We found more infiltration of CD8+ T cells, activated NK cells, and monocytes and less M0 macrophage infiltration in the asthma group than in the healthy control group. In addition, ASAH1, ACER3, and SGPP1 were negatively correlated with CD8+ T cells and activated NK cells, but positively correlated with M0 macrophages. Within the ceRNA network, SNHG9-hsa-miR-615-3p-ACER3, hsa-miR-212-5p and hsa-miR-5682 may play crucial roles in asthma pathogenesis. The low expression of ASAH1 and SGPP1 in asthma was also validated in the GSE74075 dataset. After SNHG9 knockdown, miR-615-3p expression was significantly upregulated, while that of ACER3 was significantly downregulated. Conclusion ASAH1, ACER3 and SGPP1 might be diagnostic biomarkers for asthma, and are associated with increased immune system activation. In addition, SNHG9-hsa-miR-615-3p-ACER3 may be viewed as effective therapeutic targets for asthma. Our findings might provide a novel perspective for future research on asthma.

Small cell lung cancer (SCLC) is characterized by early metastatic dissemination and rapid emergence of chemoresistance resulting in a very dismal prognosis. SCLC tumors are characterized by nearly universal biallelic inactivation of TP53 and RB1 genes and are classified into four molecular subtypes based on the expression of lineage-specific transcription factors. The integration of information encoded by the coding and non-coding genome has significantly improved our understanding of the contribution of various non-coding RNAs (ncRNAs), such as long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs), in the pathogenesis of SCLC. This has led to the concept of competing endogenous RNAs (ceRNAs) where the presence of the same miRNA response elements in one or more coding and ncRNAs may result in them competing for the same miRNA. Several studies have looked at the role of lncRNAs and circRNAs as ceRNAs by constructing ceRNA regulatory networks (ceRNETs). In this review, we discuss the role of ceRNETs in regulating various cancer hallmarks, including cell proliferation, invasion, migration, EMT, apoptosis, and chemoresistance of SCLC cells. We also discuss the potential of lncRNAs and circRNAs as biomarkers for diagnosis, prognosis, and predicting chemoresistance of SCLC.

In recent years, burgeoning research has underscored the pivotal role of non-coding RNA in orchestrating the growth, development, and pathogenesis of various diseases across organisms. However, despite these advances, our understanding of the specific contributions of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) to lens development remains notably limited. Clarifying the intricate gene regulatory networks is imperative for unraveling the molecular underpinnings of lens-related disorders. In this study, we aimed to address this gap by conducting a comprehensive analysis of the expression profiles of messenger RNAs (mRNAs), lncRNAs, and circRNAs at critical developmental time points of the mouse lens, encompassing both embryonic (E10.5, E12.5, and E16.5) and postnatal stages (P0.5, P10.5, and P60). Leveraging RNA-sequencing technology, we identified key transcripts pivotal to lens development. Our analysis revealed differentially expressed (DE) mRNAs, lncRNAs, and circRNAs across various developmental stages. Particularly noteworthy, there were 1831 co-differentially expressed (CO-DE) mRNAs, 150 CO-DE lncRNAs, and 13 CO-DE circRNAs identified during embryonic stages. Gene Ontology (GO) enrichment analysis unveiled associations primarily related to lens development, DNA conformational changes, and angiogenesis among DE mRNAs and lncRNAs. Furthermore, employing protein–protein interaction networks, mRNA–lncRNA co-expression networks, and circRNA–microRNA–mRNA networks, we predicted candidate key molecules implicated in lens development. Our findings underscore the pivotal roles of lncRNAs and circRNAs in this process, offering fresh insights into the pathogenesis of lens-related disorders and paving the way for future exploration in this field.

Background: Hepatocellular carcinoma (HCC) continues to endanger human health worldwide. Regulatory networks of competing endogenous RNAs (ceRNAs) play important roles in HCC. TP53 is the second most often altered gene in HCC and has a significant role in regulating target genes such as miRNAs and lncRNAs. Methods: Data from patients with TP53 mutation were collected through the cBioPortal database and differential analysis was performed to screen RNAs related to TP53 mutation. The lncRNA-miRNA–mRNA relationship was predicted by the miRcode, miRDB, and TargetScan databases. The ceRNA networks were screened and visualized by Cytoscape. Core ceRNA networks were generated by differential analysis, coexpression analysis, prognostic analysis and subcellular localization. Finally, methylation, mutation, PPI, GSEA, immunity and drug sensitivity analyses of MEX3A were performed to determine the role of MEX3A in HCC. Results: We identified 1508 DEmRNAs, 85 DEmiRNAs and 931 DElncRNAs and obtained a ceRNA network including 28 lncRNAs, 4 miRNAs and 31 mRNAs. Twenty hub DERNAs in the TP53-altered-related ceRNA network were screened out by Cytoscape and the core ceRNA network (LINC00491/TCL6-hsa-miR-139-5p-MEX3A) was obtained by multiple analyses. In addition, we discovered that the methylation level of MEX3A was decreased and the mutation frequency was raised in HCC. Furthermore, elevated MEX3A expression was associated with alterations in the HCC immunological microenvironment. Conclusion: We successfully constructed a reciprocal ceRNA network, which could provide new ideas for exploring HCC mechanisms and therapeutic approaches.

Abstract Background Breast cancer (BC) is a prevalent malignancy with complex etiology and varied clinical behavior. Long non-coding RNAs (lncRNAs) have emerged as key regulators in cancer progression, including BC. Among these, lncRNA TDRKH-AS1 has been implicated in several cancers, but its role in BC remains unclear. Methods We conducted a comprehensive investigation to elucidate the role of TDRKH-AS1 in BC. Clinical samples were collected from BC patients, and BC cell lines were cultured. Bioinformatics analysis using the starBase database was carried out to assess TDRKH-AS1 expression levels in BC tissue samples. Functional experiments, including knockdown, colony formation, CCK-8, Transwell, and wound-healing assays, were conducted to determine the role of TDRKH-AS1 in BC cell proliferation and invasion. Luciferase reporter and RIP assays were used to examine the interactions between TDRKH-AS1 and miR-134-5p. In addition, the downstream target gene of miR-134-5p, cAMP response element-binding protein 1 (CREB1), was identified and studied using various methods, including RT-qPCR, immunoprecipitation, and rescue experiments. In vivo experiments using mouse tumor xenograft models were conducted to examine the role of TDRKH-AS1 in BC tumorigenesis. Results TDRKH-AS1 was found to be significantly upregulated in BC tissues and cell lines. High TDRKH-AS1 expression correlated with advanced BC stages and worse patient outcomes. Knockdown of TDRKH-AS1 led to decreased BC cell proliferation and invasion. Mechanistically, TDRKH-AS1 acted as a sponge for miR-134-5p, thereby reducing the inhibitory effects of miR-134-5p on CREB1 expression. Overexpression of CREB1 partially rescued the effects of TDRKH-AS1 knockdown in BC cells. In vivo studies further confirmed the tumor-promoting role of TDRKH-AS1 in BC. Conclusions Our study unveiled a novel regulatory axis involving TDRKH-AS1, miR-134-5p, and CREB1 in BC progression. TDRKH-AS1 functioned as an oncogenic lncRNA by promoting BC cell proliferation and invasion through modulation of the miR-134-5p/CREB1 axis. These findings highlighted TDRKH-AS1 as a potential diagnostic biomarker and therapeutic target for BC treatment.

IntroductionLow-temperature stress negatively impacts rice yield, posing a significant risk to food security. While previous studies have explored the physiological and linear gene expression alterations in rice under low-temperature conditions, the changes in competing endogenous RNA (ceRNA) networks remain largely unexamined.MethodsWe conducted RNA sequencing on two japonica rice varieties with differing cold-tolerance capabilities to establish ceRNA networks. This enabled us to investigate the transcriptional regulatory network and molecular mechanisms that rice employs in response to low-temperature stress.ResultsWe identified 364 differentially expressed circular RNAs (circRNAs), 224 differentially expressed microRNAs (miRNAs), and 12,183 differentially expressed messenger RNAs (mRNAs). WRKY family was the most prominent transcription factor family involved in cold tolerance. Based on the expression patterns and targeted relationships of these differentially expressed RNAs, we discerned five potential ceRNA networks related to low-temperature stress in rice: osa-miR166j-5p from the miR166 family was associated with cold tolerance; osa-miR528-3p and osa-miR156j-3p were linked to stress response; and osa-miR156j-3p was involved in the antioxidant system. In addition, Os03g0152000 in the antioxidant system, as well as Os12g0491800 and Os05g0381400, correlated with the corresponding stress response and circRNAs in the network. A gene sequence difference analysis and phenotypic validation of Os11g0685700 (OsWRKY61) within the WRKY family suggested its potential role in regulating cold tolerance in rice.Discussion and conclusionWe identified Os11g0685700 (OsWRKY61) as a promising candidate gene for enhancing cold tolerance in japonica rice. The candidate miRNAs, mRNAs, and circRNAs uncovered in this study are valuable targets for researchers and breeders. Our findings will facilitate the development of cold-tolerant rice varieties from multiple angles and provide critical directions for future research into the functions of cold-tolerance-related miRNAs, mRNAs, and circRNAs in rice.

IntroductionMenstrual blood-derived stem cells (MenSCs) are vital in treating many degenerative and traumatic disorders. However, the underlying molecular mechanisms remain obscure in MenSCs-treating spinal cord injury (SCI) rats.MethodsMenSCs were adopted into the injured sites of rat spinal cords at day 7 post surgery and the tissues were harvested for total RNA sequencing analysis at day 21 after surgery to investigate the expression patterns of RNAs. The differentially expressed genes (DEGs) were analyzed with volcano and heatmap plot. DEGs were sequentially analyzed by weighted gene co-expression network, functional enrichment, and competitive endogenous RNAs (ceRNA) network analysis. Next, expression of selected miRNAs, lncRNAs, circRNAs and mRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR). Bioinformatics packages and extra databases were enrolled to scoop the genes functions and their interaction relationships.ResultsA total of 89 lncRNAs, 65 circRNAs, 120 miRNAs and 422 mRNAs were significantly upregulated and 65 lncRNAs, 72 circRNAs, 74 miRNAs, and 190 mRNAs were significantly downregulated in the MenSCs treated rats compared to SCI ones. Current investigation revealed that MenSCs treatment improve the recovery of the injured rats and the most significantly involved pathways in SCI regeneration were cell adhesion molecules, nature killer cell mediated cytotoxicity, primary immunodeficiency, chemokine signaling pathway, T cell receptor signaling pathway and B cell receptor signaling pathway. Moreover, the lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA ceRNA network of SCI was constructed. Finally, the protein-protein interaction (PPI) network was constructed using the top 100 DE mRNAs. The constructed PPI network included 47 nodes and 70 edges.DiscussionIn summary, the above results revealed the expression profile and potential functions of differentially expressed (DE) RNAs in the injured spinal cords of rats in the MenSCs-treated and SCI groups, and this study may provide new clues to understand the mechanisms of MenSCs in treating SCI.

The prognosis of patients with advanced renal cell carcinoma (RCC) has improved with newer therapies, including molecular-targeted therapies and immuno-oncology agents. Despite these therapeutic advances, many patients with metastatic disease remain uncured. Inhibition of glycogen synthase kinase-3β (GSK-3β) is a promising new therapeutic strategy for RCC; however, the precise regulatory mechanism has not yet been fully elucidated. MicroRNAs (miRNAs) act as post-translational regulators of target genes, and we investigated the potential regulation of miRNAs on GSK-3β in RCC. We selected nine candidate miRNAs from three databases that could potentially regulate GSK-3β. Among these, hsa-miR-4465 (miR-4465) was downregulated in RCC cell lines and renal cancer tissues. Furthermore, luciferase assays revealed that miR-4465 directly interacted with the 3′ untranslated region of GSK-3β, and Western blot analysis showed that overexpression of miR-4465 significantly decreased GSK-3β protein expression. Functional assays showed that miR-4465 overexpression significantly suppressed cell invasion of A498 and Caki-1 cells; however, cell proliferation and migration were suppressed only in Caki-1 and A498 cells, respectively, with no effect on cell cycle and apoptosis. In conclusion, miR-4465 regulates GSK-3β expression but does not consistently affect RCC cell function as a single molecule. Further comprehensive investigation of regulatory networks is required in this field.

Ulcerative colitis (UC), an inflammatory bowel disease (IBD), may increase the risk of colorectal cancer (CRC) by activating chronic proinflammatory pathways. The goal of this study was to find serum prediction biomarkers in UC to CRC development by combining low-density miRNA microarray and biocomputational approaches. The UC and CRC miRNA expression profiles were compared by low-density miRNA microarray, finding five upregulated miRNAs specific to UC progression to CRC (hsa-let-7d-5p, hsa-miR-16-5p, hsa-miR-145-5p, hsa-miR-223-5p, and hsa-miR-331-3p). The circRNA/miRNA/mRNA competitive endogenous RNA (ceRNA) network analysis showed that the candidate miRNAs were connected to well-known colitis-associated CRC ACVR2A, SOCS1, IGF2BP1, FAM126A, and CCDC85C mRNAs, and circ-SHPRH circRNA. SST and SCARA5 genes regulated by hsa-let-7d-5p, hsa-miR-145-5p, and hsa-miR-331-3p were linked to a poor survival prognosis in a CRC patient dataset from The Cancer Genome Atlas (TCGA). Lastly, our mRNA and miRNA candidates were validated by comparing their expression to differentially expressed mRNAs and miRNAs from colitis-associated CRC tissue databases. A high level of hsa-miR-331-3p and a parallel reduction in SOCS1 mRNA were found in tissue and serum. We propose hsa-miR-331-3p and possibly hsa-let-7d-5p as novel serum biomarkers for predicting UC progression to CRC. More clinical sample analysis is required for further validation.

Introduction: IgA nephropathy (IgAN) is a prevalent worldwide glomerular disease with a complex pathophysiology that has significant economic implications. Despite the lack of successful research, this study aims to discover the potential competing endogenous RNA (ceRNA) network of autophagy-associated genes in IgAN and examine their correlation with immune cell infiltration. Methods: Autophagy-related hub genes were discovered by assessing the GSE116626 dataset and constructing a protein-protein interaction network. Nephroseq v5 analysis engine was used to analyze correlations between hub genes and proteinuria, glomerular filtration rate (GFR), and serum creatinine levels. Then, a ceRNA network construction and the CIBERSORT tool for immune cell infiltration analysis were also performed. Additionally, the differentially expressed autophagy-related genes were used to predict potential targeted medications for IgAN. Results: Overall, 1,396 differentially expressed genes were identified in IgAN along with 25 autophagy-related differentially expressed messenger RNAs. Enrichment analysis revealed significant involvement of autophagy and apoptosis in biological processes. Next, we evaluated the top hub nodes based on their highest degrees. The ability of IgAN discrimination was confirmed in the GSE35487 and GSE37460 datasets by validating the five hub genes: SIRT1, FOS, CCL2, CDKN1A, and MYC. In the Nephroseq v5 analysis engine, the clinical correlation of the five hub genes was confirmed. Furthermore, the ceRNA network identified 18 circular RNAs and 2 microRNAs associated with hub autophagy-related genes in IgAN. Our investigation identified hsa-miR-32-3p and hsa-let-7i-5p as having elevated expression levels and substantial diagnostic value. Finally, four distinctively infiltrated immune cells were found to be associated with the hub autophagy-related genes, and 67 drugs were identified as potential therapeutic options for IgAN. Conclusion: This study sheds light on a novel ceRNA regulatory network mechanism associated with autophagy in IgAN development.

Abstract Background Abdominal fat deposition depends on both the proliferation of preadipocytes and their maturation into adipocytes, which is a well-orchestrated multistep process involving many regulatory molecules. Circular RNAs (circRNAs) have emergingly been implicated in mammalian adipogenesis. However, circRNA-mediated regulation in chicken adipogenesis remains unclear. Our previous circRNA sequencing data identified a differentially expressed novel circRNA, 8:27,886,180|27,889,657, during the adipogenic differentiation of chicken abdominal preadipocytes. This study aimed to investigate the regulatory role of circDOCK7 in the proliferation and adipogenic differentiation of chicken abdominal preadipocytes, and explore its molecular mechanisms of competing endogenous RNA underlying chicken adipogenesis. Results Our results showed that 8:27,886,180|27,889,657 is an exonic circRNA derived from the head-to-tail splicing of exons 19–22 of the dedicator of cytokinesis 7 (DOCK7) gene, abbreviated as circDOCK7. CircDOCK7 is mainly distributed in the cytoplasm of chicken abdominal preadipocytes and is stable because of its RNase R resistance and longer half-life. CircDOCK7 is significantly upregulated in the abdominal fat tissues of fat chickens compared to lean chickens, and its expression gradually increases during the proliferation and adipogenic differentiation of chicken abdominal preadipocytes. Functionally, the gain- and loss-of-function experiments showed that circDOCK7 promoted proliferation, G0/G1- to S-phase progression, and glucose uptake capacity of chicken abdominal preadipocytes, in parallel with adipogenic differentiation characterized by remarkably increased intracellular lipid droplet accumulation and triglyceride and acetyl coenzyme A content in differentiated chicken abdominal preadipocytes. Mechanistically, a pull-down assay and a dual-luciferase reporter assay confirmed that circDOCK7 interacted with gga-miR-301b-3p, which was identified as an inhibitor of chicken abdominal adipogenesis. Moreover, the ACSL1 gene was demonstrated to be a direct target of gga-miR-301b-3p. Chicken ACSL1 protein is localized in the endoplasmic reticulum and mitochondria of chicken abdominal preadipocytes and acts as an adipogenesis accelerator. Rescue experiments showed that circDOCK7 could counteract the inhibitory effects of gga-miR-301b-3p on ACSL1 mRNA abundance as well as the proliferation and adipogenic differentiation of chicken abdominal preadipocytes. Conclusions CircDOCK7 serves as a miRNA sponge that directly sequesters gga-miR-301b-3p away from the ACSL1 gene, thus augmenting adipogenesis in chickens. These findings may elucidate a new regulatory mechanism underlying abdominal fat deposition in chickens.

Abstract Background The differential expression of circular RNAs (circRNAs) in individuals with very severe chronic obstructive pulmonary disease (COPD) and healthy individuals was screened using microarray technology. The related functions and mechanisms were analyzed using bioinformatic methods to explore the potential of target circRNAs as biomarkers of COPD and provide insights for future pathogenesis. Patients and methods Thirty patients with very severe COPD and thirty healthy controls were diagnosed at The Second People’s Hospital of Hefei from September 2021 to September 2022. The differential expression of circRNAs was compared and analyzed using a gene microarray and verified using quantitative real-time polymerase chain reaction (qRT-PCR) technology. Results A total of 90 upregulated and 29 downregulated circRNAs were screened in patients with very severe COPD and compared with those in healthy controls. qRT-PCR analysis showed that hsa_circ_0062683 of patients with very severe COPD was significantly upregulated, and hsa_circ_0089763 and hsa_circ_0008882 were significantly downregulated. By constructing the circRNA-miRNA interaction network, it was found that hsa-miR-612, hsa-miR-593-5p, hsa-miR-765, and hsa-miR-103a-2-5p are the miRNAs regulated by more differentially expressed circRNAs (DEcircRNAs). DEcircRNAs may participate in the development of COPD through hypoxia or regulation of various immune cells. Conclusion Plasma circRNAs may play a helpful role in the diagnosis and assessment of COPD and be valuable disease biomarkers.

Abstract Background Gastric cancer (GC) carries significant morbidity and mortality globally. An increasing number of studies have confirmed that circular RNA (circRNA) is tightly associated with the carcinogenesis and development of GC, especially acting as a competing endogenous RNA for miRNAs. Objective Our study aimed to construct the circRNA-miRNA-mRNA regulatory network and analyze the function and prognostic significance of the network using bioinformatics tools. Methods We first downloaded the GC expression profile from the Gene Expression Omnibus database and identified differentially expressed genes and differentially expressed circRNAs. Then, we predicted the miRNA-mRNA interaction pairs and constructed the circRNA-miRNA-mRNA regulatory network. Next, we established a protein-protein interaction network and analyzed the function of these networks. Finally, we primarily validated our results by comparison with The Cancer Genome Atlas cohort and by performing qRT-PCR. Results We screened the top 15 hub genes and 3 core modules. Functional analysis showed that in the upregulated circRNA network, 15 hub genes were correlated with extracellular matrix organization and interaction. The function of downregulated circRNAs converged on physiological functions, such as protein processing, energy metabolism and gastric acid secretion. We ascertained 3 prognostic and immune infiltration-related genes, COL12A1, COL5A2, and THBS1, and built a nomogram for clinical application. We validated the expression level and diagnostic performance of key prognostic differentially expressed genes. Conclusions In conclusion, we constructed two circRNA-miRNA-mRNA regulatory networks and identified 3 prognostic and screening biomarkers, COL12A1, COL5A2, and THBS1. The ceRNA network and these genes could play important roles in GC development, diagnosis and prognosis.

Abstract Background Retinal neovascularization (RNV) is a leading cause of blindness worldwide. Long non-coding RNA (lncRNA) and competing endogenous RNA (ceRNA) regulatory networks play vital roles in angiogenesis. The RNA-binding protein galectin-1 (Gal-1) participates in pathological RNV in oxygen-induced retinopathy mouse models. However, the molecular associations between Gal-1 and lncRNAs remain unclear. Herein, we aimed to explore the potential mechanism of action of Gal-1 as an RNA-binding protein. Results A comprehensive network of Gal-1, ceRNAs, and neovascularization-related genes was constructed based on transcriptome chip data and bioinformatics analysis of human retinal microvascular endothelial cells (HRMECs). We also conducted functional enrichment and pathway enrichment analyses. Fourteen lncRNAs, twenty-nine miRNAs, and eleven differentially expressed angiogenic genes were included in the Gal-1/ceRNA network. Additionally, the expression of six lncRNAs and eleven differentially expressed angiogenic genes were validated by qPCR in HRMECs with or without siLGALS1. Several hub genes, such as NRIR, ZFPM2-AS1, LINC0121, apelin, claudin-5, and C-X-C motif chemokine ligand 10, were found to potentially interact with Gal-1 via the ceRNA axis. Furthermore, Gal-1 may be involved in regulating biological processes related to chemotaxis, chemokine-mediated signaling, the immune response, and the inflammatory response. Conclusions The Gal-1/ceRNA axis identified in this study may play a vital role in RNV. This study provides a foundation for the continued exploration of therapeutic targets and biomarkers associated with RNV.

Bladder cancer is a common disease associated with high rates of morbidity and mortality. Although immunotherapy approaches such as adoptive T-cell therapy and immune checkpoint blockade have been investigated for the treatment of bladder cancer, their off-target effects and ability to affect only single targets have led to clinical outcomes that are far from satisfactory. Therefore, it is important to identify novel targets that can effectively control tumor growth and metastasis. It is well known that long noncoding RNAs (lncRNAs) are powerful regulators of gene expression. Increasing evidence has shown that dysregulated lncRNAs in bladder cancer are involved in cancer cell proliferation, migration, invasion, apoptosis, and epithelial-mesenchymal transition (EMT). In this review, we focus on the roles and underlying mechanisms of lncRNA-mediated competing endogenous RNA (ceRNA) networks in the regulation of bladder cancer progression. In addition, we discuss the potential of targeting lncRNA-mediated ceRNA networks to overcome cancer treatment resistance and its association with clinicopathological features and outcomes in bladder cancer patients. We hope this review will stimulate research to develop more effective therapeutic approaches for bladder cancer treatment.

The current situation of hepatocellular carcinoma (HCC) management is challenging due to its high incidence, mortality, recurrence and metastasis. Recent advances in gene genetic and expression regulation have unveiled the significant role of non-coding RNA (ncRNA) in various cancers. This led to the formulation of the competing endogenous RNA (ceRNA) hypothesis, which posits that both coding RNA and ncRNA, containing miRNA response elements (MRE), can share the same miRNA sequence. This results in a competitive network between ncRNAs, such as lncRNA and mRNA, allowing them to regulate each other. Extensive research has highlighted the crucial role of the ceRNA network in HCC development, impacting various cellular processes including proliferation, metastasis, cell death, angiogenesis, tumor microenvironment, organismal immunity, and chemotherapy resistance. Additionally, the ceRNA network, mediated by lncRNA or circRNA, offers potential in early diagnosis and prevention of HCC. Consequently, ceRNAs are emerging as therapeutic targets for HCC. The complexity of these gene networks aligns with the multi-target approach of traditional Chinese medicine (TCM), presenting a novel perspective for TCM in combating HCC. Research is beginning to show that TCM compounds and prescriptions can affect HCC progression through the ceRNA network, inhibiting proliferation and metastasis, and inducing apoptosis. Currently, the lncRNAs TUG1, NEAT1, and CCAT1, along with their associated ceRNA networks, are among the most promising ncRNAs for HCC research. However, this field is still in its infancy, necessitating advanced technology and extensive basic research to fully understand the ceRNA network mechanisms of TCM in HCC treatment.

Accumulating evidence showed that competing endogenous RNA (ceRNA) mechanism plays a pivotal role in salt sensitivity of blood pressure (SSBP). We constructed a ceRNA network based on SSBP-related differently expressed lncRNAs (2), mRNAs (73) and miRNAs (18). Bioinformatic analyses were utilized to analyze network and found network genes participate in biological pathways related to SSBP pathogenesis such as regulation of nitric oxide biosynthetic process (GO:0045,428) and cellular response to cytokine stimulus (GO:0071,345). Fourteen candidate ceRNA pathways were selected from network to perform qRT-PCR validation and found nine RNAs (KCNQ1OT1, SLC8A1-AS1, IL1B, BCL2L11, KCNJ15, CX3CR1, KLF2, hsa-miR-362–5p and hsa-miR-423–5p) differently expressed between salt-sensitive (SS) and salt-resistant (SR) groups (P

Background: Acute myeloid leukemia (AML) is characterized by a high recurrence and mortality rate. Cuproptosis is involved in cell death regulation in in a variety of solid tumors. Long non-coding RNAs that regulate cuproptosis genes in the pathogenesis of acute leukemia have yet to be explored. Methods: First, cuproptosis genes with distinct expression levels were discovered by contrasting AML with normal samples from the TCGA and GTEx cohorts. Pearson correlation and univariate Cox-regression analysis were performed to identify cuproptosis-associated lncRNAs with significant prognostic values. Then the least absolute shrinkage and selection operator (LASSO) Cox regression was utilized to establish a multi-gene signature to predict AML prognosis. Next, Kaplan-Meier estimator, receiver operating characteristic curve, and a nomogram were performed to evaluate the predictive capacity of the risk signature. Functional enrichment analyses were employed to assess their function. Moreover, qRT-PCR testing of lncRNA expression in AML samples was conducted. The competing endogenous RNA (ceRNA) network was constructed to find the target genes. Results: A risk model based on the signature of three cuproptosis-associated lncRNAs was developed. The results showed that the model possessed excellent prognostic potential. The nomogram raised the accuracy in predicting AML survival. In addition, functional enrichment analyses demonstrated an enrichment of inflammatory and immune-related pathways. Moreover, correlations between the risk signature and clinicopathological variables, tumor mutational burden, RNA stemness score, immune profile, and drug sensitivity were observed. Furthermore, we discovered that TRAF3IP2-AS1 may function as a ceRNA to regulate cuproptosis and ferroptosis gene expression. Conclusion: The risk signature established in this study could serve as a reliable biosignature for AML prognosis. And the findings presented here may facilitate research on cuproptosis in AML.

As the critical components of tumor microenvironment, cancer-associated fibroblasts (CAFs) support the development of various type of cancers, including laryngeal squamous cell carcinoma (LSCC), but the detailed molecular mechanisms by which cancer-associated fibroblasts interact with LSCC cells to facilitate its progression have not been fully uncovered. In the present study, by analyzing the contents from normal fibroblasts (NFs) and cancer-associated fibroblasts-derived extracellular vesicles (EVs) with Real-Time qPCR analysis, we found that the tumor-initiating LncRNA TUC338 was significantly upregulated in the cancer-associated fibroblasts-derived extracellular vesicles, compared to the normal fibroblasts-secreted extracellular vesicles. Further experiments confirmed that cancer-associated fibroblasts-derived extracellular vesicles promoted cell proliferation, colony formation abilities, epithelial-mesenchymal transition (EMT) and tumorigenesis of LSCC cells via delivering LncRNA TUC338. The mechanical experiments verified that LncRNA TUC338 was stabilized by METTL3/YTHDF1-mediated N6-methyladenosine (m6A) modifications, and elevated LncRNA TUC338 sponged miR-8485 to upregulate chromobox homolog 2 (CBX2) in LSCC cells in a competing endogenous RNA mechanisms-dependent manner. Moreover, our rescue experiments evidenced that cancer-associated fibroblasts-derived LncRNA TUC338-containing extracellular vesicles-induced supportive effects in LSCC aggressiveness were all abrogated by overexpressing miR-8485 and silencing CBX2. Collectively, this study is the first to identify a novel m6A/LncRNA TUC338/miR-8485/CBX2 axis in CAFs-EVs-mediated LSCC development, and to show its potential as a diagnostic biomarker for LSCC.

Objective: Circular RNAs (circRNAs) have been shown to participate in various cancers via sponging miRNAs (microRNAs). However, their role in lung adenocarcinoma (LUAD) remains elusive.Methods: The transcriptome data and corresponding clinical information of lung adenocarcinoma samples were extracted from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed genes (DEgenes) were identified and further used to constructed a circRNA-associated competing endogenous RNA (ceRNA) network. Real-Time qPCR analysis was conducted to examine gene expression at transcriptional level. The regulatory mechanisms of circRNA-miRNA-gene were validated by dual-luciferase reporter array and RNA pull-down assay. Cell growth, migration and invasion were evaluated by CCK-8 assay, colony formation assay and transwell assay, respectively.Results: Based on public microarray data, we systematically constructed a circRNA-associated ceRNA network including 11 DEcircRNAs, 8 DEmiRNAs and 49 DEgenes. Among the ceRNA network, we found that circ-0002727 was a key regulatory and was further confirmed to be upregulated in LUAD cancer cells. Subsequently, we found that silencing of circ-0002727 significantly suppressed the LUAD cell proliferation, migration and invasion in vitro. Mechanistically, we showed that circ-0002727 could competitively bind miR-144-3p to enhance the KIF14 expression in LUAD cells. Rescue assays indicated that circ-0002727 could regulate LUAD cell proliferation through modulating miR-144-3p/KIF14 pathway. Besides, KIF14 expression level was positively correlated with TNM stage and metastasis, and patients with high KIF14 expression suffered poor prognosis.Conclusion: Taken together, our study revealed that circ-0002727 could act as a ceRNA to regulate LUAD progression via modulating miR-144-3p/KIF14 pathway, providing a potential therapeutic target for LUAD.

Abstract Background Abnormal miRNA and mRNA expression and dysregulated immune microenvironment have been found to frequently induce the progression of hepatocellular carcinoma (HCC) in recent reports. In particular, the immune-related competing endogenous RNAs (ceRNA) mechanism plays a crucial role in HCC progression. However, the underlying mechanisms remain unclear. Methods Differentially expressed immune-related genes were obtained from the Immport, GEO, and TCGA databases. The mRNA and protein expression levels in HCC tissues and adjacent normal tissues were confirmed, and we further investigated the methylation levels of these biomarkers to explore their function. Then, the TIMER and TISCH databases were used to assess the relationship between immune infiltration and hub genes. Survival analysis and univariate and multivariate Cox models were used to evaluate the association between hub genes and HCC diagnosis. Hub gene expression was experimentally validated in six HCC cell lines and 15 HCC samples using qRT-PCR and immunohistochemistry. The hub genes were uploaded to DSigDB for drug prediction enrichment analysis. Results We identified that patients with abnormal miRNAs (hsa-miR-125b-5p and hsa-miR-21-5p) and their targeted genes (NTF3, PSMD14, CD320, and SORT1) had a worse prognosis. Methylation analysis of miRNA-targeted genes suggested that alteration of methylation levels is also a factor in the induction of tumorigenesis. We also found that the development of HCC progression caused by miRNA-mRNA interactions may be closely correlated with the infiltration of immunocytes. Moreover, the GSEA, GO, and KEGG analysis suggested that several common immune-related biological processes and pathways were related to miRNA-targeted genes. The results of qRT-PCR, immunohistochemistry, and western blotting were consistent with our bioinformatics results, suggesting that abnormal miRNAs and their targeted genes may affect HCC progression. Conclusions Briefly, our study systematically describes the mechanisms of miRNA-mRNA interactions in HCC and predicts promising biomarkers that are associated with immune filtration for HCC progression.

Abstract Background Our previous studies have confirmed that miR-154-5p can regulate pRb expression, and thus, play a tumor suppressor role in HPV16 E7-induced cervical cancer. However, its upstream molecules have not been elucidated in the progression of cervical cancer. This study aimed to explore the role of the miR-154-5p upstream molecule, hsa_circ_0000276 in cervical cancer development and its possible mechanisms of action. Methods We detected differences in whole transcriptome expression profiles of cervical squamous carcinoma and tissues adjacent to cervical cancer tissues from patients using microarray technology to predict circular RNAs (circRNAs) with binding sites to miR-154-5p. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to detect the expression of hsa_circ_0000276 (which had the strongest binding capacity to miR-154 and was selected as the target molecule) in cervical cancer tissues, followed by in vitro functional assays. Downstream microRNAs (miRNAs) and mRNAs of hsa_circ_0000276 were identified using transcriptome microarray data and databases, while the protein–protein interaction networks were obtained using STRING. A competing endogenous RNA (ceRNA) network centered on hsa_circ_0000276 was constructed using Cytoscape and GO and KEGG databases. Abnormal expression and prognosis of critical downstream molecules were analyzed using gene databases and molecular experiments. qRT-PCR and western blot analysis was performed to verify the expression of candidate genes. Results We identified 4,001 differentially expressed circRNAs between HPV16-positive cervical squamous carcinoma and benign cervical tissues and 760 circRNAs targeting miR-154-5p, including hsa_circ_0000276. hsa_circ_0000276 and miR-154-5p directly bound, and hsa_circ_0000276 was upregulated, in cervical precancerous lesions and cervical cancer tissues and cells. Silencing hsa_circ_0000276 inhibited G1/S transition and cell proliferation and promoted apoptosis in SiHa and CaSki cells. Bioinformatics analysis showed that the hsa_circ_0000276 ceRNA network included 17 miRNAs and seven mRNAs, and downstream molecules of hsa_circ_0000276 were upregulated in cervical cancer tissues. These downstream molecules were associated with a poor prognosis and affected cervical cancer-associated immune infiltration. Of these, expression of CD47, LDHA, PDIA3, and SLC16A1 was downregulated in sh_hsa_circ_0000276 cells. Conclusions Our findings show that hsa_circ_0000276 exerts cancer-promoting effects in cervical cancer and is an underlying biomarker for cervical squamous cell carcinoma.

Abstract Background Stroke attributable to atrial fibrillation (AF related stroke, AFST) accounts for 13 ~ 26% of ischemic stroke. It has been found that AFST patients have a higher risk of disability and mortality than those without AF. Additionally, it’s still a great challenge to treat AFST patients because its exact mechanism at the molecular level remains unclear. Thus, it’s vital to investigate the mechanism of AFST and search for molecular targets of treatment. Long non-coding RNAs (lncRNAs) are related to the pathogenesis of various diseases. However, the role of lncRNAs in AFST remains unclear. In this study, AFST-related lncRNAs are explored using competing endogenous RNA (ceRNA) network analysis and weighted gene co-expression network analysis (WGCNA). Methods GSE66724 and GSE58294 datasets were downloaded from GEO database. After data preprocessing and probe reannotation, differentially expressed lncRNAs (DELs) and differentially expressed mRNAs (DEMs) between AFST and AF samples were explored. Then, functional enrichment analysis and protein-protein interaction (PPI) network analysis of the DEMs were performed. At the meantime, ceRNA network analysis and WGCNA were performed to identify hub lncRNAs. The hub lncRNAs identified both by ceRNA network analysis and WGCNA were further validated by Comparative Toxicogenomics Database (CTD). Results In all, 19 DELs and 317 DEMs were identified between the AFST and AF samples. Functional enrichment analysis suggested that the DEMs associated with AFST were mainly enriched in the activation of the immune response. Two lncRNAs which overlapped between the three lncRNAs identified by the ceRNA network analysis and the 28 lncRNAs identified by the WGCNA were screened as hub lncRNAs for further validation. Finally, lncRNA GAS6-AS1 turned out to be associated with AFST by CTD validation. Conclusion These findings suggested that low expression of GAS6-AS1 might exert an essential role in AFST through downregulating its downstream target mRNAs GOLGA8A and BACH2, and GAS6-AS1 might be a potential target for AFST therapy.

Low-temperature chilling is a major abiotic stress leading to reduced rice yield and is a significant environmental threat to food security. Low-temperature chilling studies have focused on physiological changes or coding genes. However, the competitive endogenous RNA mechanism in rice at low temperatures has not been reported. Therefore, in this study, antioxidant physiological indices were combined with whole-transcriptome data through weighted correlation network analysis, which found that the gene modules had the highest correlation with the key antioxidant enzymes superoxide dismutase and peroxidase. The hub genes of the superoxide dismutase-related module included the UDP-glucosyltransferase family protein, sesquiterpene synthase and indole-3-glycerophosphatase gene. The hub genes of the peroxidase-related module included the WRKY transcription factor, abscisic acid signal transduction pathway-related gene plasma membrane hydrogen-ATPase and receptor-like kinase. Therefore, we selected the modular hub genes and significantly enriched the metabolic pathway genes to construct the key competitive endogenous RNA networks, resulting in three competitive endogenous RNA networks of seven long non-coding RNAs regulating three co-expressed messenger RNAs via four microRNAs. Finally, the negative regulatory function of the WRKY transcription factor OsWRKY61 was determined via subcellular localization and validation of the physiological indices in the mutant.

Objective: The mechanisms of chondrocytes ferroptosis in osteoarthritis (OA) have not yet been fully elucidated. This study aimed to identify key ferroptosis related genes (FRGs) involved in chondrocytes ferroptosis. Methods: LASSO, SVM-RFE, and receiver operating characteristic curve (ROC) were performed to screen key differentially expressed FRGs (DEFRGs). Functional analyses were conducted using GO, and KEGG analyses. Unsupervised clustering analysis was used to identify ferroptosis related patterns. The CeRNA network was constructed to predict the upstream miRNAs and lncRNAs. Finally, we validated the role of EGFR in chondrocytes ferroptosis using in vivo and in vitro experiments. Results: A total of 42 DEFRGs were identified between OA and normal cartilages. GO and KEGG analyses indicated that these DEFRGs were significantly engaged in ferroptosis related biological processes and pathways, such as cellular response to oxidative stress, positive regulation of programmed cell death, MAPK and PI3K-Akt signaling pathways. Moreover, four key DEFRGs, including ACSF2, AURKA, EGFR, and KLHL24, were considered as potential biomarkers of OA. Moreover, two distinct ferroptosis related patterns were determined, and a total of 882 differentially expressed genes were identified which might participate in extracellular matrix degradation and inflammatory response. In addition, the CeRNA network showed that EGFR could be competitively regulated by 3 lncRNAs and 4 miRNAs. Significantly, the expression of EGFR was downregulated in human OA cartilages, OA mouse model, and erastin induced chondrocytes. EGFR inhibition could induce the occurrence of chondrocytes ferroptosis and ECM degradation which could be reversed by the addition of Ferrostatin-1. Conclusion: Our study has identified ACSF2, AURKA, EGFR, and KLHL24 as ferroptosis-related biomarkers in OA. Furthermore, we have conducted a preliminary investigation into the role of EGFR in regulating chondrocytes ferroptosis. These findings offer novel insights into the molecular mechanisms underlying OA.

Long noncoding RNAs (lncRNAs) play significant roles in the regulation of mRNA expression or in shaping the competing endogenous RNA (ceRNA) network by targeting miRNA. The insect gut is one of the most important tissues due to direct contact with external pathogens and functions in the immune defense against pathogen infection through the innate immune system and symbionts, but there are limited observations on the role of the lncRNA-involved ceRNA network of the Toll/Imd pathway and correlation analysis between this network and bacterial microbiota in the Altica viridicyanea gut. In this research, we constructed and sequenced six RNA sequencing libraries using normal and antibiotic-reared samples, generating a total of 17,193 lncRNAs and 26,361 mRNAs from massive clean data by quality control and bioinformatic analysis. Furthermore, a set of 8,539 differentially expressed lncRNAs (DELs) and 13,263 differentially expressed mRNAs (DEMs), of which related to various immune signaling pathways, such as the Toll/Imd, JAK/STAT, NF-κB, and PI3K-Akt signaling pathways, were obtained between the two experimental groups in A. viridicyanea. In addition, numerous GO and KEGG enrichment analyses were used to annotate the DELs and their target genes. Moreover, six Toll family members and nineteen signal genes from the Toll/Imd signaling pathway were identified and characterized using online tools, and phylogenetic analyses of the above genes proved their classification. Next, a lncRNA-miRNA-mRNA network of the Toll/Imd pathway was built, and it contained different numbers of DEMs in this pathway and related DELs based on prediction and annotation. In addition, qRT-PCR validation and sequencing data were conducted to show the expression patterns of the above DELs and DEMs related to the Toll/Imd signaling pathway. Finally, the correlated investigations between DELs or DEMs of the Toll/Imd signaling pathway and most changes in the gut bacterial microbiota revealed significantly positive or negative relationships between them. The present findings provide essential evidence for innate immune ceRNAs in the beetle gut and uncover new potential relationships between innate immune pathways and the gut bacterial microbiota in insects.

Background: Cardiac fibrosis has gradually gained significance in the field of cardiovascular disease; however, its specific pathogenesis remains unclear. This study aims to establish the regulatory networks based on whole-transcriptome RNA sequencing analyses and reveal the underlying mechanisms of cardiac fibrosis. Methods: An experimental model of myocardial fibrosis was induced using the chronic intermittent hypoxia (CIH) method. Expression profiles of long non-coding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) were acquired from right atrial tissue samples of rats. Differentially expressed RNAs (DERs) were identified, and functional enrichment analysis was performed. Moreover, a protein-protein interaction (PPI) network and competitive endogenous RNA (ceRNA) regulatory network that are related to cardiac fibrosis were constructed, and the relevant regulatory factors and functional pathways were identified. Finally, the crucial regulators were validated using qRT-PCR. Results: DERs, including 268 lncRNAs, 20 miRNAs, and 436 mRNAs, were screened. Further, 18 relevant biological processes, such as "chromosome segregation, " and 6 KEGG signaling pathways, such as "cell cycle, " were significantly enriched. The regulatory relationship of miRNA–mRNA–KEGG pathways showed eight overlapping disease pathways, including "pathways in cancer." In addition, crucial regulatory factors, such as Arnt2, WNT2B, GNG7, LOC100909750, Cyp1a1, E2F1, BIRC5, and LPAR4, were identified and verified to be closely related to cardiac fibrosis. Conclusion: This study identified the crucial regulators and related functional pathways in cardiac fibrosis by integrating the whole transcriptome analysis in rats, which might provide novel insights into the pathogenesis of cardiac fibrosis.

Abstract Objective Osteoporosis is a progressive systemic skeletal disorder. Multiple profiling studies have contributed to characterizing biomarkers and therapeutic targets for osteoporosis. However, due to the limitation of the platform of miRNA sequencing, only a part of miRNA can be sequenced based on one platform. Materials and methods In this study, we performed miRNA sequencing in osteoporosis bone samples based on a novel platform Illumina Hiseq 2500. Bioinformatics analysis was performed to construct osteoporosis-related competing endogenous RNA (ceRNA) networks. Gene interference and osteogenic induction were used to examine the effect of identified ceRNA networks on osteogenic differentiation of human bone marrow-derived mesenchymal stem cells (HBMSCs). Results miR-1303 was lowly expressed, while cochlin (COCH) and KCNMA1-AS1 were highly expressed in the osteoporosis subjects. COCH knockdown improved the osteogenic differentiation of HBMSCs. Meanwhile, COCH inhibition compensated for the suppression of osteogenic differentiation of HBMSCs by miR-1303 knockdown. Further, KCNMA1-AS1 knockdown promoted osteogenic differentiation of HBMSCs through downregulating COCH by sponging miR-1303. Conclusions Our findings suggest that the KCNMA1-AS1/miR-1303/COCH axis is a promising biomarker and therapeutic target for osteoporosis.

Abstract Purpose Steroid-induced osteonecrosis of the femoral head (SONFH) was a refractory orthopedic hip joint disease in the young and middle-aged people, but the pathogenesis of SONFH remained unclear. We aimed to identify the potential genes and screen potential therapeutic compounds for SONFH. Methods The microarray was obtained for blood tissue from the GEO database, and then it identifies differentially expressed genes (DEGs). The DEGs were analyzed to obtain the differences in immune cell infiltration. The gene functional enrichment analysis of SONFH was analyzed. The PPI of DEGs was identified through the STRING database, and the cluster modules and hub genes were ascertained using MCODE and CytoHubba, and the ROC curve of hub genes was analyzed, and the tissues distribution of hub genes was understood by the HPA, Bgee and BioGPS databases. The hub genes and target miRNAs and corresponding upstream lncRNAs were predicted by TargetScan, miRDB and ENCORI database. Subsequently, we used CMap, DGIdb and L1000FWD databases to identify several potential therapeutic molecular compounds for SONFH. Finally, the AutoDockTools Vina, PyMOL and Discovery Studio were employed for molecular docking analyses between compounds and hub genes. Results The microarray dataset GSE123568 was obtained related to SONFH. There were 372 DEGs including 197 upregulated genes and 175 downregulated genes by adjusted P value 1. Several significant GSEA enrichment analysis and biological processes and KEGG pathway associated with SONFH were identified, which were significantly related to cytoskeleton organization, nucleobase-containing compound catabolic process, NOD-like receptor signaling pathway, MAPK signaling pathway, FoxO signaling pathway, neutrophil-mediated immunity, neutrophil degranulation and neutrophil activation involved in immune response. Activated T cells CD4 memory, B cells naïve, B cells memory, T cells CD8 and T cells gamma delta might be involved in the occurrence and development of SONFH. Three cluster modules were identified in the PPI network, and eleven hub genes including FPR2, LILRB2, MNDA, CCR1, IRF8, TYROBP, TLR1, HCK, TLR8, TLR2 and CCR2 were identified by Cytohubba, which were differed in bone marrow, adipose tissue and blood, and which had good diagnostic performance in SONFH. We identified IRF8 and 10 target miRNAs that was utilized including Targetsan, miRDB and ENCORI databases and 8 corresponding upstream lncRNAs that was revealed by ENCORI database. IRF8 was detected with consistent expression by qRT-PCR. Based on the CMap, DGIdb and L1000FWD databases, the 11 small molecular compounds that were most strongly therapeutic correlated with SONFH were estradiol, genistein, domperidone, lovastatin, myricetin, fenbufen, rosiglitazone, sirolimus, phenformin, vorinostat and vinblastine. All of 11 small molecules had good binding affinity with the IRF8 in molecular docking. Conclusion The occurrence of SONFH was associated with a “multi-target” and “multi-pathway” pattern, especially related to immunity, and IRF8 and its noncoding RNA were closely related to the development of SONFH. The CMap, DGIdb and L1000FWD databases could be effectively used in a systematic manner to predict potential drugs for the prevention and treatment of SONFH. However, additional clinical and experimental research is warranted.

Tao Zuo,1,* Ping Jiang,2,* Jun Fu,1 Yongjian Zhang1 1Department of Thoracic Surgery, The Central Hospital of Wuhan, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, People’s Republic of China; 2Department of Ophthalmology, Zhongnan Hospital of Wuhan University, Wuhan, People’s Republic of China*These authors contributed equally to this workCorrespondence: Tao Zuo, Department of Thoracic Surgery, The Central Hospital of Wuhan, Tongji Medical College, Huazhong University of Science and Technology, No. 26, Shengli Street, Jiang’an District, Wuhan, People’s Republic of China, Tel +86 15002786691, Email zuotaowuhan@126.comPurpose: Gefitinib resistance limits the therapeutic efficacy of gefitinib to lung adenocarcinoma (LUAD). The goal of this research is to learn more about the lncRNA AFAP1-AS1 and how it functions in gefitinib-resistant LUAD cells.Methods: RT-qPCR was performed to test the expression of AFAP1-AS1, miR-653-5p and AGR2 in LUAD tissues with acquired resistance to gefitinib or not as well as in gefitinib-resistant LUAD cells. Cell proliferation, invasion and apoptosis were measured by CCK8 assays, transwell invasion assays and flow cytometry, respectively. Luciferase reporter assay showed that miR-653-5p and AFAP1-AS1 or AGR2 interactions.Results: In gefitinib-resistant LUAD cells and tissues, AFAP1-AS1 was overexpressed. Meanwhile, silencing AFAP1-AS1 reduced proliferation and migration while increasing apoptosis and gefitinib sensitivity. Mechanically, AFAP1-AS1 sequestered the miR-653-5p and blocked the inhibition of miR-653-5p to AGR2 and stepwise upregulated AGR2 overexpression in LUAD gefitinib resistant cells, resulting gefitinib resistance in LUAD.Conclusion: AFAP1-AS1 promotes gefitinib-resistance LUAD cells through a previously unrecognized miR-653-5p/AGR2 axis, suggesting targeting AFAP1-AS1/miR-653-5p/AGR2 axis might be a promising way for LUAD intervention.Keywords: lncRNA, gefitinib resistance, competing endogenous RNA, lung adenocarcinoma

Objective To construct a regulatory network of long noncoding RNA (lncRNA)-related competing endogenous RNAs (ceRNA) in low folate combined with methotrexate (MTX) induced neural tube defects (NTDs) mouse model, and to explore the potential lncRNAs regulation in low folate-induced NTDs. Methods Forty SPF-grade 7~8-week-old healthy female C57BL/6J mice (body weight 18~22 g) and 20 healthy male C57BL/6J mice (20~22 g) were selected to this experiment. After a folate deficiency induced fetal mouse NTD model was successfully established, the embryonic brain tissues were extracted for total RNA and RNA high-throughput sequencing (RNA-seq) was performed. The results of RNA-seq were used to carry out co-expression network construction, GO function enrichment of mRNA in co-expression network and KEGG pathway enrichment analysis, construction of lncRNA-related ceRNA regulation network and other bioinformatics analysis. And the relative expression of lncRNA Xist, Sema4f and Padi2 in the fetal mouse and cell samples was verified by qPCR. Results A total of 47 differential lncRNAs (Log2FC >0.585 or Log2FC < -0.585, FDR < 0.05) were detected by RNA-seq, of which 38 were up-regulated and 9 were down-regulated. The co-expression network between lncRNA and mRNA was constructed. The results of KEGG/GO analysis on mRNAs in the co-expression network indicated that tube development (P < 0.001), collagen metabolic process (P=0.002), system development (P=0.002) and other significant enrichments were highly relevant to NTDs. Besides, KEGG pathway analysis displayed that the PI3K-Akt pathway was closely related to NTDs; lncRNA Xist, LOC102633899, 1700086L19Rik and other targeting differential genes Sema4f, Padi2, Lin28a were found in the lncRNA-related ceRNA network. qPCR confirmed that the above results in the fetal mice and cell samples under folic acid deficiency, with consistent corresponding sequences, and obviously down-regulated (P < 0.05). Conclusion In folic acid deficiency-induced mouse NTD model, lncRNA Xist, and its targeting genes Sema4f and Padi2 are obviously down-regulated.

Abstract Background Among the most lethal cancers, pancreatic adenocarcinoma (PAAD) is an essential component of digestive system malignancies that still lacks effective diagnosis and treatment methods. As exosomes and competing endogenous RNA (ceRNA) regulatory networks in tumors go deeper, we expect to construct a ceRNA regulatory network derived from blood exosomes of PAAD patients by bioinformatics methods and develop a survival prediction model based on it. Methods Blood exosome sequencing data of PAAD patients and normal controls were downloaded from the exoRbase database, and the expression profiles of exosomal mRNA, lncRNA, and circRNA were differentially analyzed by R. The related mRNA, circRNA, lncRNA, and their corresponding miRNA prediction data were imported into Cytoscape software to visualize the ceRNA network. Then, we conducted GO and KEGG enrichment analysis of mRNA in the ceRNA network. Genes that express differently in pancreatic cancer tissues compared with normal tissues and associate with survival (P

Background: Cervical cancer (CC) is a common cancer in female, which is associated with problems like poor prognosis. Circular RNA (circRNA) is a kind of competing endogenous RNA (ceRNA) that has an important role in regulating microRNA (miRNA) in many cancers. The regulatory mechanisms of CC immune microenvironment and the transcriptome level remain to be fully explored. Methods: In this study, we constructed the ceRNA network through the interaction data and expression matrix of circRNA, miRNA and mRNA. Meanwhile, based on the gene expression matrix, CIBERSORT algorithm was used to reveal contents of tumor-infiltrating immune cells (TIICs). Then, we screened prognostic markers based on ceRNA network and immune infiltration and constructed two nomograms. In order to find immunological differences between the high- and low-risk CC samples, we examined multiple immune checkpoints and predicted the effect of PD-L1 ICI immunotherapy. In addition, the sensitive therapeutics for high-risk patients were screened, and the potential agents with anti-CC activity were predicted by Connective Map (CMap). Results: We mapped a ceRNA network including 5 circRNAs, 17 miRNAs and 129 mRNAs. From the mRNA nodes of the network six genes and two kind of cells were identified as prognostic makers for CC. Among them, there was a significant positive correlation between CD8+ T cells and SNX10 gene. The results of TIDE and single sample GSEA (ssGSEA) showed that T cells CD8 do play a key role in inhibiting tumor progression. Further, our study screened 24 drugs that were more sensitive to high-risk CC patients and several potential therapeutic agents for reference. Conclusions: Our study identified several circRNA-miRNA-mRNA regulatory axes and six prognostic genes based on the ceRNA network. In addition, through TIIC, survival analysis and a series of immunological analyses, T cells were proved to be good prognostic markers, besides play an important role in the immune process. Finally, we screened 24 potentially more effective drugs and multiple potential drug compounds for high- and low-risk patients.

Abstract Background Chronic lymphocytic leukemia (CLL) is a heterogeneous B-cell malignancy that lacks specific biomarkers and drug targets. Competing endogenous RNAs (ceRNAs) play vital roles in oncogenesis and tumor progression by sponging microRNAs (miRNAs). Nevertheless, the regulatory mechanisms of survival-related ceRNA networks in CLL remain to be uncovered. Methods We included 865 de novo CLL patients to investigate RNA expression profiles and Illumina sequencing was performed on four CLL patients, two CLL cell lines and six healthy donors in our center. According to univariate Cox regression, LASSO regression as well as multivariate Cox regression analyses, we established a novel risk score model in CLL patients. Immune signatures were compared between the low- and high-risk groups with CIBERSORT and ESTIMATE program. Afterwards, we analyzed the relationship between differentially expressed miRNAs (DEmiRNAs) and IGHV mutational status, p53 mutation status and del17p. Based on the survival analyses and differentially expressed RNAs with targeting relationships, the lncRNA/circRNA-miRNA-mRNA ceRNA networks were constructed. In addition, the circRNA circ_0002078/miR-185-3p/TCF7L1 axis was verified and their interrelations were delineated by dual-luciferase reporter gene assay. Results Totally, 57 differentially expressed mRNAs (DEmRNAs) and 335 DEmiRNAs were identified between CLL patient specimens and normal B cells. A novel risk score model consisting of HTN3, IL3RA and NCK1 was established and validated. The concordance indexes of the model were 0.825, 0.719 and 0.773 in the training, test and total sets, respectively. The high-risk group was related to del(13q14) as well as shorter overall survival (OS). Moreover, we identified DEmiRNAs that related to cytogenetic abnormality of CLL patients, which revealed that miR-324-3p was associated with IGHV mutation, p53 mutation and del17p. The survival-related lncRNA/circRNA-miRNA-mRNA ceRNA networks were constructed to further facilitate the development of potential predictive biomarkers. Besides, the expression of circ_0002078 and TCF7L1 were significantly elevated and miR-185-3p was obviously decreased in CLL patients. Circ_0002078 regulated TCF7L1 expression by competing with TCF7L1 for miR-185-3p. Conclusions The comprehensive analyses of RNA expression profiles provide pioneering insights into the molecular mechanisms of CLL. The novel risk score model and survival-related ceRNA networks promote the development of prognostic biomarkers and potential therapeutic vulnerabilities for CLL.

Abstract Background Increasing studies have indicated that noncoding RNA (ncRNA)-mediated competing endogenous RNA (ceRNA) network serves as a significant role in cancer progression, but the underlying regulatory mechanisms of which in gastric cancer (GC) remain largely unclear. Methods Based on Gene Expression Omnibus and The Cancer Genome Atlas datasets, potential biomarkers for GC were screened and validated by machine learning. Then, upstream regulatory ncRNA of potential biomarkers was identified to construct a novel ceRNA network in GC through means of stepwise reverse prediction and validation. Ultimately, tumor immune cell infiltration analysis was performed based on the EPIC algorithm. Results A total of 188 differentially expressed genes (DEGs) were screened, and three candidate diagnostic biomarkers (FAP, PSAPL1, and SERPINH1) for GC were identified and validated. Subsequently, H19 and miR-378a-5p were identified as upstream regulatory ncRNAs that could potentially bind SERPINH1 in GC. Moreover, Immune infiltration analysis revealed that each component in the ceRNA network (H19/miR-378a-5p/SERPINH1) was significantly correlated with the infiltration abundances of diverse tumor-infiltrating immune cells. Conclusions H19 may regulate the immune cell infiltration in carcinogenesis of GC through miR-378a-5p/SERPINH1 signaling.

Abstract Circular RNAs (circRNAs) are important transcriptional regulatory RNA molecule that can regulate the transcription of downstream genes by competitive binding of miRNAs or coding proteins or by blocking mRNAs translation. Numerous studies have shown that circRNAs are extensively involved in cell proliferation, differentiation and apoptosis, gene transcription and signal transduction. Fat deposition and muscle development have important effects on beef traits. CircRNAs are involved in regulating bovine fat and muscle cells and are differentially expressed in the tissues composed of these cells, suggesting that circRNAs play an important role in regulating bovine fat formation and muscle development. This review describes differential expression of circRNAs in bovine fat and muscle tissues, research progress in understanding how circRNAs regulate the proliferation and differentiation of bovine fat and muscle cells through competing endogenous RNAs networks, and provide a reference for the subsequent research on the molecular mechanism of circRNAs in regulating fat deposition and muscle development in cattle.

Abstract Background Intervertebral disc degeneration (IDD) has become a serious public health problem, the mechanism of which is complex and still unclear. We aimed to construct a ceRNA network related to IDD to explore its pathogenesis. Methods We downloaded the GSE67566, GSE63492, GSE116726 and GSE124272 datasets from GEO database, and obtained the differentially expressed RNAs. Then, we constructed a ceRNA network and the KEGG and GO enrichment analysis were performed. Finally, we performed immune cell infiltration analysis on the GSE124272 dataset and analysed the correlation between immune cell abundance and hub genes expression levels. Results The ceRNA network included three down-regulated circRNAs: hsa_circ_0074817, hsa_circ_0002702, hsa_circ_0003600, three up-regulated miRNAs: hsa-miR-4741, hsa-miR-3158-5p, hsa-miR-508-5p, and 57 down-regulated mRNAs, including six hub genes: IGF1, CHEK1, CCNB1, OIP5, BIRC5, AR. GO and KEGG analysis revealed that the network is involved in various biological functions. Immune infiltration analysis showed that IDD was closely related to immune cell infiltration, and hub genes could further affect the development of IDD by affecting immune cell infiltration. Conclusion This study identified the hsa_circ_0074817-hsa-miR-508-5p-IGF1/CHEK1/CCNB1, the hsa_circ_0003600-hsa-miR-4741-BIRC5/OIP5/AR and the hsa_circ_0002702-hsa-miR-3158-5p-IGF1/AR as important regulatory axis of IDD, which will help us gain further insight into the pathogenesis of IDD and determine potential therapeutic targets.

Abstract Background Liver cirrhosis is a well-known risk factor for hepatocellular carcinoma (HCC). However, some HCC cases can also originate from non-cirrhotic livers. The aim of this study was to identify key circular RNAs (circRNAs) associated with the tumorigenesis of non-cirrhotic liver disease. Methods The differently expressed circRNAs between non-cirrhotic and cirrhotic HCCs were assessed with use of high-throughput circRNAs sequencing and validated with quantitative reverse transcription polymerase chain reaction (qRT-PCR). Potential biological functions of these dysregulated circRNAs were predicted with use of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. A circRNA-miRNA-mRNA regulation network was constructed as achieved with use of miRanda software and visualized using Cytoscape software. Biological functions of the four most prominent dysregulated circRNAs identified were confirmed by in vitro experiments. Moreover, possible translations of these dysregulated circRNAs were also predicted. Results A total of 393 dysregulated circRNAs were identified between non-cirrhotic and cirrhotic HCC, including 213 that were significantly up-regulated and 180 significantly down-regulated circRNAs. Expression levels of the six most prominent dysregulated circRNAs were further validated using qRT-PCR. Many tumor related miRNAs were involved in the circRNA-miRNA-mRNA networks, including miR-182-5p, miR-561-3p, miR-125a-5p, miR-145, miR-23b-3p and miR-30e-3p, and downstream mRNAs of dysregulated circRNAs were significantly related with biological processes involved in the progression of tumors, including proliferation, migration, differentiation, and focal adhesion. Results from the in vitro experiments demonstrated that the most prominent dysregulated circRNAs exerted notable effects upon the proliferation and migration of HCC cells. Finally, we also identified 19 dysregulated circRNAs having potential for the coding of functional peptides. Conclusion The results of this present study indicate that circRNAs may play important roles in tumorigenesis of non-cirrhotic HCC. Such findings provide some novel insights and pave the way for the development of future studies directed at investigating the initiation and treatment of HCC.

Background Exploring the regulatory network of competing endogenous RNAs (ceRNAs) as hallmarks for breast cancer development has great significance and could provide therapeutic targets. An mRNA signature predictive of prognosis and therapy response in BRCA carriers was developed according to circular RNA homeodomain-interacting protein kinase 3 (circHIPK3)-based ceRNA network. Method We constructed a circHIPK3-based ceRNA network based on GSE173766 dataset and identified potential mRNAs that were associated with BRCA mutation patients within this ceRNA network. A total of 11 prognostic mRNAs and a risk model were identified and developed by univariate Cox regression analysis and the LASSO regression analysis as well as stepAIC method. Genomic landscape was treated by mutect2 and fisher. Immune characteristics was analyzed by ESTIMATE, MCP-counter. TIDE analysis was conducted to predict immunotherapy. The clinical treatment outcomes of BRCA mutation patients were assessed using a nomogram. The proliferation, migration and invasion in breast cancer cell lines were examined using CCK8 assay and transwell assay. Result We found 241 mRNAs within the circHIPK3-based ceRNA network. An 11 mRNA-based signature was identified for prognostic model construction. High risk patients exhibited dismal prognosis, low response to immunotherapy, less immune cell infiltration and tumor mutation burden (TMB). High-risk patients were sensitive to six anti-tumor drugs, while low-risk patient were sensitive to 47 drugs. The risk score was the most effective on evaluating patients’ survival. The robustness and good prediction performance were validated in The Cancer Genome Atlas (TCGA) dataset and immunotherapy datasets, respectively. In addition, circHIPK3 mRNA level was upregulated, and promoted cell viability, migration and invasion in breast cancer cell lines. Conclusion The current study could improve the understanding of mRNAs in relation to BRCA mutation and pave the way to develop mRNA-based therapeutic targets for breast cancer patients with BRCA mutation.

Objective: Using bioinformatics analyses, this study aimed to identify lncRNAs related to the immune status of acute myeloid leukemia (AML) patients and ascertain the potential impact in immunity-related competing endogenous RNA (ceRNA) networks on AML prognosis.Methods: AML-related RNA-seq FPKM data, AML-related miRNA expression microarray data, and gene sets associated with immunity-related pathways were, respectively, obtained from the TCGA, GEO, and ImmReg databases. An immunity-related ceRNA network was then constructed according to the predicted interactions between AML-related mRNAs, lncRNAs, and miRNAs. After performing LASSO and multivariate Cox regression analyses, lncRNAs in the ceRNA network were used to establish an AML prognostic model. According to mutual regulatory relationships and consistent trends of expression among candidate ceRNAs, two ceRNA subnetworks related to the AML prognostic model were determined. Finally, the correlation between the expression levels of mRNAs, lncRNAs, and miRNAs in each ceRNA subnetwork and immune cell infiltration (assessed by combining the ESTIMATE and CIBERSORT methods and ssGSEA) was analyzed.Results: A total of 424 immunity-related differentially expressed (IR-DE) mRNAs (IR-DEmRNAs), 191 IR-DElncRNAs, and 69 IR-DEmiRNAs were obtained, and a ceRNA network of 20 IR-DElncRNAs, 6 IR-DEmRNAs, and 3 IR-DEmiRNAs was established. Univariate Cox regression analysis was conducted on 20 IR-DElncRNAs, and 7 of these were identified to be significantly correlated with the overall survival (OS) time in AML patients. Then, two IR-DElncRNAs (MEG3 and HCP5) were screened as independent OS-related factors by LASSO and multivariable Cox regression analyses, and a prognostic model was constructed to evaluate the survival risk in AML patients. Survival analyses indicated that the OS of patients was often poor in the high-risk group. Additionally, from this model, two ceRNA regulatory pathways, namely, MEG3/miR-125a-5p/SEMA4C and HCP5/miR-125b-5p/IL6R, which were potentially involved in the immune regulation of AML prognosis were identified.Conclusion: lncRNAs HCP5 and MEG3 may act as key ceRNAs in the pathogenesis in AML by regulating immune cell representation as part of the regulatory lncRNA-miRNA-mRNA axes. The candidate mRNAs, lncRNAs, and miRNAs included in the ceRNA network identified here may serve as useful prognostic biomarkers and immunotherapeutic targets for AML.

Yan Lin,1,2,&ast; Gaofeng Zhu,3,&ast; Xiaoyong Li,4,&ast; Huaxiao Yu,5 Yuhang Luo,5 Jiaqiong Lin,6 Renyuan Li,7 Zena Huang2 1Department of Nephrology, The Third Affiliated Hospital of Guangzhou Medical University, Guangzhou, People’s Republic of China; 2Yunkang School of Medicine and Health, Nanfang College, Guangzhou, People’s Republic of China; 3Department of General Medicine, Guangdong Provincial People’s Hospital, Guangdong Academy of Medical Sciences, Guangzhou, People’s Republic of China; 4Department of General Surgery, The Third Affiliated Hospital of Guangzhou Medical University, Guangzhou, People’s Republic of China; 5The Third Clinical School, Guangzhou Medical University, Guangzhou, People’s Republic of China; 6Affiliated Dongguan Maternal and Child Healthcare Hospital, Southern Medical University, Guangzhou, People’s Republic of China; 7Department of Endocrinology and Metabolism, The Third Affiliated Hospital of Guangzhou Medical University, Guangzhou, People’s Republic of China&ast;These authors contributed equally to this workCorrespondence: Zena Huang, Yunkang School of Medicine and Health, Nanfang College, Guangzhou, People’s Republic of China, Tel +86-13570466614, Email huangzn@nfu.edu.cn Renyuan Li, Department of Endocrinology and Metabolism, The Third Affiliated Hospital of Guangzhou Medical University, Guangzhou, People’s Republic of China, Tel +86-18926146852, Email lirenyuan@gzhmu.edu.cnBackground: Icariin presents protective effect in several kidney diseases. However, the role of icariin in contrast-induced acute kidney injury (CIAKI) is still unclear. This study aimed to investigate the effect of icariin in CIAKI, as well as exploring the underlying mechanism from the aspect of interaction between protein-coding genes and non-coding RNAs.Methods: The effect of icariin was evaluated in both in vivo and in vitro CIAKI models. Rat kidneys were collected for genome-wide sequencing. The differentially expressed genes (DEGs) were screened and visualized by R software. The function annotation of DEGs was analyzed by Metascape. By Cytoscape software, the competing endogenous RNA (ceRNA) network was constructed, and hub genes were selected. Expressions of hub genes were validated by PCR. Association of hub genes in the ceRNA network and renal function was also examined.Results: Icariin protected against CIAKI in both in vivo and in vitro models. Based on DEGs in icariin pretreated CIAKI rats, lncRNA- and circRNA-associated ceRNA networks were constructed, respectively. Function annotation showed the ceRNA networks were enriched in ERK1 and ERK2 cascade, MAPK signaling and NF-κB signaling. Further, two circRNAs, six lncRNAs, four miRNAs and nine mRNAs were selected as hub genes of the ceRNA network. Among them, eight mRNAs (Acot1, Cbwd1, Ly6i, Map3k14, Mettl2b, Nyap1, Set and Utp20) were negatively correlated with renal function, while one mRNA (Tmem44) was positively correlated with renal function.Conclusion: Icariin presented a protective effect against CIAKI. The ceRNA network, involving Acot1, Cbwd1, Ly6i, Map3k14, Mettl2, Nyap1, Set, Tmem44 and Utp20, might partially contribute to the underlying mechanism of icariin protection by regulation of ERK1 and ERK2 cascade, MAPK signaling and NF-κB signaling.Keywords: icariin, contrast medium, acute kidney injury, competing endogenous RNA, non-coding RNA

Abstract LINC00894 plays an important role in cancer cell proliferation and invasion in breast and kidney cancer. However, its role in thyroid cancer proliferation and metastasis remains unclear. In this study, data on LINC00894 expression in thyroid cancer tissues were obtained from GEPIA2. miRNA expression in thyroid cancer tissues was obtained from starBase 3.0 and OncomiR. Cell proliferation was evaluated using CCK-8, and Transwell chambers were used for the migration and invasion assays. LINC00894 and let-7e-5p expressions in thyroid cancer cells were measured using qRT–PCR. Meanwhile, TIA-1 expression in thyroid cancer cells was analyzed via western blotting. We found that LINC00894 expression was markedly reduced in thyroid cancer tissues and cells, and low expression of LINC00894 was associated with poor prognosis in thyroid cancer. LINC00894 overexpression inhibited the proliferation, migration, and invasion of CAL-62 and TPC-1 cells. Additionally, let-7e-5p expression was substantially enhanced in CAL-62 and TPC-1 cells. LINC00894 overexpression promoted TIA-1 expression by acting as a sponge of let-7e-5p. Finally, let-7e-5p weakened the function of LINC00894 in thyroid cancer cells via reduction in TIA-1 levels. In conclusion, our data suggest that increased LINC00894 expression reduces the oncogenic properties of thyroid cancer cells by sponging let-7e-5p to promote TIA-1 expression.

Highlights 1. An IDD-related ceRNA network involving four circRNAs, three miRNAs and 11 mRNAs was constructed. 2. The expression levels of miR-185-5p, miR-486-5p, FOXO1, SMAD2 and TGFB1 were consistent in different databases. 3. Our study identified IDD-related circRNA–miRNA–mRNA interaction axes, including circRNA_100086–miR-509-3p–MAPK1.

Xing-Zhi Wang,1,2 Shuo Li,3 Yun Liu,3 Gui-Yun Cui,2 Fu-Ling Yan3 1School of Medicine, Southeast University, Nanjing, Jiangsu, 210009, People’s Republic of China; 2Department of Neurology, The Affiliated Hospital of Xuzhou Medical University, Xuzhou, Jiangsu, 221006, People’s Republic of China; 3Department of Neurology, Affiliated ZhongDa Hospital, School of Medicine, Southeast University, Nanjing, Jiangsu, 210009, People’s Republic of ChinaCorrespondence: Fu-Ling Yan, Tel +861 391 396 9651, Fax +8602583262251, Email yanfuling218@163.comBackground and Purpose: Accumulating evidence suggests that circular RNAs (circRNAs) are involved in immune and inflammatory processes after acute ischemic stroke (AIS). However, the roles of circRNA-mediated competing endogenous RNA (ceRNA) in modulating immune inflammation of AIS have not yet been determined. This study aimed to construct a circRNA-mediated immune-related ceRNA network and identify novel circRNAs in AIS.Methods: Microarray data were downloaded from the GEO database and further analysed by R software. Then, we constructed a circRNA-mediated ceRNA network based on interaction information from the bioinformatics database. A topological property analysis of the ceRNA network was conducted to screen novel circRNAs. Finally, we further applied quantitative real-time polymerase chain reaction (qRT-PCR) to two independent sets.Results: We constructed an AIS immune-related ceRNA (AISIRC) network containing immune-related genes (IRGs), miRNAs, and circRNAs. Additionally, we extracted the subnetwork from the AISIRC network and screened six immune-related circRNAs. After identification and validation, we finally confirmed that plasma levels of circPTP4A2 and circTLK2 were significantly increased in AIS patients compared with both healthy control subjects (HCs) and transient ischemic attack (TIA) patients. Logistic regression and receiver-operating characteristic (ROC) curve analyses demonstrated that these two circRNAs may function as predictive and discriminative biomarkers for AIS. We also confirmed that plasma levels of circPTP4A2 were elevated in TIA patients compared with HCs and might be an independent risk factor for predicting TIA. Longitudinal analysis of circRNA expression up to 90 days after AIS indicated that the ability of circPTP4A2 and circTLK2 to monitor AIS dynamics was highly desirable.Conclusion: In summary, the circRNA-mediated immune-related ceRNA network was successfully constructed, and two circulating circRNAs (circPTP4A2 and circTLK2) improved sensitivity for the diagnosis of AIS and could be considered diagnostic biomarkers.Keywords: circular RNA, acute ischemic stroke, immune mechanism, competing endogenous RNA, biomarker

ObjectiveLong non-coding RNAs (lncRNAs) play pivotal roles in the transcriptional regulation of atrial fibrillation (AF) by acting as competing endogenous RNAs (ceRNAs). In the present study, the expression levels of lncRNAs of sinus rhythm (SR) patients and AF patients were investigated with transcriptomics technology, and the lncRNA-miRNA-mRNA network based on the ceRNA theory in AF was elaborated.MethodsLeft atrial appendage (LAA) tissues were obtained from patients with valvular heart disease during cardiac surgery, and they were divided into SR and AF groups. The expression characterizations of differentially expressed (DE) lncRNAs in the two groups were revealed by high-throughput sequencing methods. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed, and the lncRNA-miRNA-mRNA-mediated ceRNA network was constructed.ResultsA total of differentially expressed 82 lncRNAs, 18 miRNAs, and 495 mRNAs in human atrial appendage tissues were targeted. Compared to SR patients, the following changes were found in AF patients: 32 upregulated and 50 downregulated lncRNAs; 7 upregulated and 11 downregulated miRNAs; and 408 upregulated and 87 downregulated mRNAs. A lncRNA-miRNA-mRNA network was constructed, which included 44 lncRNAs, 18 miRNAs, and 347 mRNAs. qRT-PCR was performed to verify these findings. GO and KEGG analyses suggested that inflammatory response, chemokine signaling pathway, and other biological processes play important roles in the pathogenesis of AF. Network analysis based on the ceRNA theory identified that lncRNA XR_001750763.2 and Toll-like receptor 2 (TLR2) compete for binding to miR-302b-3p. In AF patients, lncRNA XR_001750763.2 and TLR2 were upregulated, and miR-302b-3p was downregulated.ConclusionWe identified a lncRNA XR_001750763.2/miR-302b-3p/TLR2 network based on the ceRNA theory in AF. The present study shed light on the physiological functions of lncRNAs and provided information for exploring potential treatments for AF.