RNA-Seq and miRNA-Seq data from Epstein-Barr virus-infected tree shrews reveal a ceRNA network contributing to immune microenvironment regulation

Epstein-Barr virus (EBV) infection in humans is ubiquitous and associated with various diseases. Remodeling of the immune microenvironment is the primary cause of EBV infection and pathogenesis; however, the underlying mechanism has not been fully elucidated. In this study, we used whole-transcriptome RNA-Seq to detect mRNAs, long non-coding RNAs (lncRNA), and microRNA (miRNA) profiles in the control group, 3 days, and 28 days after EBV infection, based on the tree shrew model that we reported previously. First, we estimated the proportion of 22 cell types in each sample using CIBERSORT software and identified 18 high-confidence DElncRNAs related to immune microenvironment regulation after EBV infection. Functional enrichment analysis of these differentially expressed lncRNAs primarily focused on the autophagy, endocytosis, and ferroptosis signalling pathways. Moreover, EBV infection affects miRNA expression patterns, and many miRNAs are silenced. Finally, three competing endogenous RNA regulatory networks were built using lncRNAs that significantly correlated with immune cell types, miRNAs that responded to EBV infection, and potentially targeted the mRNA of the miRNAs. Among them, MRPL42-AS-5 might act as an hsa-miR-296-5p "sponge" and compete with target mRNAs, thus increasing mRNA expression level, which could induce immune cell infiltration through the cellular senescence signalling pathway against EBV infection. Overall, we conducted a complete transcriptomic analysis of EBV infection in vivo for the first time and provided a novel perspective for further investigation of EBV-host interactions.

Comprehensive investigation of differentially expressed ncRNAs, mRNAs, and their ceRNA networks in the regulation of shell color formation in clam, Cyclina sinensis

Noncoding RNAs (ncRNAs) have gained significant attention in recent years due to their crucial roles in various biological processes. However, our understanding of the expression and functions of ncRNAs in Cyclina sinensis, an economically important marine bivalve, remains limited. This study aimed to address this knowledge gap by systematically identifying ncRNAs in the mantles of C. sinensis with purple and white shells. Through our analysis, we identified a differential expression of 1244 mRNAs, 196 lncRNAs, 49 circRNAs, and 23 miRNAs between purple- and white-shell clams. Functional enrichment analysis revealed the involvement of these differentially expressed ncRNAs in biomineralization and pigmentation processes. To gain further insights into the regulatory mechanisms underlying shell color formation, we established competitive endogenous RNA (ceRNA) networks. These networks allowed us to identify targeted differentially expressed miRNAs (DEMis) and genes associated with shell color formation. Based on the ceRNA networks, we obtained an up-down-up lncRNA-miRNA-mRNA network consisting of 13 upregulated lncRNAs and a circRNA-miRNA-mRNA network comprising three upregulated circRNAs (novel_circ_0004787, novel_circ_0001165, novel_circ_0000245). Through these networks, we identified and selected an upregulated novel gene (evm.TU.Hic_asm_7.988) and a downregulated sponge miRNA (hru-miR-1985) as potential contributors to shell color regulation. In summary, the present investigation offers a comprehensive analysis of ncRNA catalogs expressed in the clam mantle of C. sinensis. The findings enhance our comprehension of the molecular mechanisms governing shell coloration and offer new perspectives for selective breeding of C. sinensis in the future.

Identification of hub genes and potential therapeutic mechanisms related to HPV positive head and neck squamous carcinoma based on full transcriptomic detection and ceRNA network construction

Infection with human papillomavirus (HPV) is a major risk factor for head and neck squamous cell carcinoma (HNSCC). The objective of this study is to investigate the gene expression profiles and signaling pathways that are specific to HPV-positive HNSCC (HPV+ HNSCC). Moreover, a competing endogenous RNA (ceRNA) network analysis was utilized to identify the core gene of HPV+ HNSCC and potential targeted therapeutic drugs. Transcriptome sequencing analysis identified 3,253 coding RNAs and 3,903 non-coding RNAs (ncRNAs) that exhibited preferentially expressed in HPV+ HNSCC. Four key signaling pathways were selected through pathway enrichment analysis. By combining ceRNA network and protein-protein interaction (PPI) network topology analysis, RNA Polymerase II Associated Protein 2 (RPAP2), which also exhibited high expression in HPV+ HNSCC based on the TCGA database, was identified as the hub gene. Gene set enrichment analysis (GSEA) results revealed RPAP2's involvement in various signaling pathways, encompassing basal transcription factors, ubiquitin-mediated proteolysis, adherens junction, other glycan degradation, ATP-binding cassette (ABC) transporters, and oglycan biosynthesis. Five potential small molecule targeted drugs (enzastaurin, brequinar, talinolol, phenylbutazone, and afuresertib) were identified using the cMAP database, with enzastaurin showing the highest affinity for RPAP2. Cellular functional experiments confirmed the inhibitory effect of enzastaurin on cell viability of HPV+ HNSCC and RPAP2 expression levels. Additionally, enzastaurin treatment suppressed the expression levels of the top-ranked long non-coding RNA (lncRNA), circular RNA (circRNA), and microRNA (miRNA) in the ceRNA network. This study based on the ceRNA network provides valuable insights into the molecular mechanisms and potential therapeutic strategies for HPV+ HNSCC, and provide theoretical basis for the exploration of HPV+ HNSCC biomarkers and the development of targeted drugs.

A prognostic four-gene signature and a therapeutic strategy for hepatocellular carcinoma: Construction and analysis of a circRNA-mediated competing endogenous RNA network

BACKGROUND

Hepatocellular carcinoma (HCC) has a poor long-term prognosis. The competition of circular RNAs (circRNAs) with endogenous RNA is a novel tool for predicting HCC prognosis. Based on the alterations of circRNA regulatory networks, the analysis of gene modules related to HCC is feasible.

METHODS

Multiple expression datasets and RNA element targeting prediction tools were used to construct a circRNA-microRNA-mRNA network in HCC. Gene function, pathway, and protein interaction analyses were performed for the differentially expressed genes (DEGs) in this regulatory network. In the protein-protein interaction network, hub genes were identified and subjected to regression analysis, producing an optimized four-gene signature for prognostic risk stratification in HCC patients. Anti-HCC drugs were excavated by assessing the DEGs between the low- and high-risk groups. A circRNA-microRNA-hub gene subnetwork was constructed, in which three hallmark genes, KIF4A, CCNA2, and PBK, were subjected to functional enrichment analysis.

RESULTS

A four-gene signature (KIF4A, CCNA2, PBK, and ZWINT) that effectively estimated the overall survival and aided in prognostic risk assessment in the The Cancer Genome Atlas (TCGA) cohort and International Cancer Genome Consortium (ICGC) cohort was developed. CDK inhibitors, PI3K inhibitors, HDAC inhibitors, and EGFR inhibitors were predicted as four potential mechanisms of drug action (MOA) in high-risk HCC patients. Subsequent analysis has revealed that PBK, CCNA2, and KIF4A play a crucial role in regulating the tumor microenvironment by promoting immune cell invasion, regulating microsatellite instability (MSI), and exerting an impact on HCC progression.

CONCLUSIONS

The present study highlights the role of the circRNA-related regulatory network, identifies a four-gene prognostic signature and biomarkers, and further identifies novel therapy for HCC.

Role of circular RNA as competing endogenous RNA in ovarian cancer (Review)

Circular RNA (circRNA), a type of non‑coding RNA, plays a regulatory role in biological processes. The special loop structure of circRNA makes it highly stable and specific in diseased tissues and cells, especially in tumors. Competing endogenous RNAs (ceRNAs) compete for the binding of microRNA (miRNA) at specific binding sites and thus regulate gene expression. ceRNAs play an important role in various diseases and are currently recognized as the most prominent mechanism of action of circRNAs. circRNAs can modulate the proliferation, migration, invasion and apoptosis of tumor cells through the ceRNA mechanism. With further research, circRNAs may serve as novel markers and therapeutic targets for ovarian cancer (OC). In the present review, the research progress of circRNAs as ceRNAs in OC was summarized, focusing on the effects of the circRNA/miRNA/mRNA axis on the biological functions of OC cells through mediating pivotal signaling pathways. The role of circRNAs in the diagnosis, prognostic assessment and treatment of OC was also discussed in the present review.

Circ-10720 as a ceRNA adsorbs microRNA-1238 and modulates ZEB2 to boost NSCLC development by activating EMT

BACKGROUND

Circular RNAs (circRNAs) are critical regulators in the progression of tumors. This experimental design aimed to explore the mechanism of circ-10720 in non-small cell lung cancer (NSCLC).

METHODS

We used RT-qPCR to measure circ-10720 expression in clinical samples and analyzed its relationship with the clinicopathological characteristics of NSCLC patients. The expression levels of microRNA-1238 (miR-1238) and Zinc Finger E-box-binding Homeobox 2 (ZEB2) in clinical samples were detected by RT-qPCR. NSCLC cells were transfected with relevant plasmids or sequences. Circ-10720, miR-1238, and ZEB2 expressions in cells were analyzed via RT-qPCR or western blot. Cell proliferation, apoptosis, migration, and invasion were assessed with CCK-8, flow cytometry, and transwell assay, respectively. The protein expression of ZEB2 and epithelial-mesenchymal transition (EMT)-related markers (E-cadherin, Vimentin, N-cadherin) were detected via western blot. Xenograft assay was used to determine the effect of circ-10720 on NSCLC in vivo. Circ-10720 and ZEB2 expressions in tumors were detected using RT-qPCR or Western blot. Immunohistochemistry was used to evaluate E-cadherin and N-cadherin expression in tumors. Finally, the binding relationship between miR-1238 with circ-10720 or ZEB2 was verified by the bioinformatics website, dual luciferase reporter assay, RNA pull-down assay, and RIP assay.

RESULTS

Circ-10720 was upregulated in NSCLC and correlated with TNM stage of NSCLC patients. MiR-1238 was lowly expressed but ZEB2 was highly expressed in NSCLC. Circ-10720 silencing suppressed the proliferation, metastasis, and EMT of NSCLC cells. Mechanically, circ-10720 was a competitive endogenous RNA (ceRNA) for miR-1238, and ZEB2 was a target of miR-1238. circ-10720-modulated ZEB2 via competitively binding with miR-1238 to control NSCLC progression. In addition, circ-10720 knockdown suppressed tumor growth in vivo.

CONCLUSIONS

Circ-10720 acts as a ceRNA to adsorb miR-1238 and modulate ZEB2 to facilitate the proliferation, migration, invasion, and EMT of NSCLC cells.

Construction and analysis of competitive endogenous RNA networks and prognostic models associated with ovarian cancer based on the exoRBase database

OBJECTIVE

To construct a competitive endogenous RNA (ceRNA) regulatory network in blood exosomes of patients with ovarian cancer (OC) using bioinformatics and explore its pathogenesis.

METHODS

The exoRbase2.0 database was used to download blood exosome gene sequencing data from patients OC and normal controls and the expression profiles of exosomal mRNA, long non-coding RNA (lncRNA), and circular RNA (circRNA) were detected independently using R language for differential expression analysis. TargetScan and miRanda databases were combined for the prediction and differential expression of mRNA-binding microRNAs (miRNA). The miRcode and starBase databases were used to predict miRNAs that bind to differentially expressed lncRNAs and circRNAs repectively. The relevant mRNA, circRNA, lncRNA and their corresponding miRNA prediction data were imported into Cytoscape software for visualization of the ceRNA network. The R language and KEGG Orthology-based Annotation System (KOBAS) were used to execute and illustrate the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Hub genes were identified using The CytoHubba plugin.

RESULTS

Thirty-one differentially expressed mRNAs, 17 differentially expressed lncRNAs, and 24 differentially expressed circRNAs were screened. Cytoscape software was used to construct the ceRNA network with nine mRNA nodes, two lncRNA nodes, eight circRNA nodes, and 51 miRNA nodes. Both GO and KEGG were focused on the Spliceosome pathway, indicating that spliceosomes are closely linked with the development of OC, while heterogenous nuclear ribonucleoprotein K and RNA binding motif protein X-linked genes were the top 10 score Hub genes screened by Cytoscape software, including two lncRNAs, four mRNAs, and four circRNAs. In patients with OC, the expression of eukaryotic translation initiation factor 4 gamma 2 (EIF4G2), SERPINE 1 mRNA binding protein 1 (SERBP1), ribosomal protein L15 (RPL15) and human leukocyte antigen complex P5 (HCP5) was significantly higher whereas that of testis expressed transcript, Y-linked 15 and DEAD-box helicase 3 Y-linked genes was lower compared to normal controls Immunocorrelation scores revealed that SERBP1 was significantly and negatively correlated with endothelial cells and CD4+ T cells and positively correlated with natural killer (NK) cells and macrophages, respectively; RPL15 was significantly positively correlated with macrophages and endothelial cells and negatively correlated with CD8+ T cells and uncharacterized cells, respectively. EIF4G2 was significantly and negatively correlated with endothelial cells and CD4+ T cells, and positively correlated with uncharacterized cells, respectively. Based on the survival data and the significant correlation characteristics derived from the multifactorial Cox analysis (P < 0.05), the survival prediction curves demonstrated that the prognostic factors associated with 3-year survival in patients with OC were The prognostic factors associated with survival were Macrophage, RPL15.

CONCLUSION

This study successfully constructs a ceRNA regulatory network in blood exosomes of OV patients, which provides the specific targets for diagnosis and treatment of OC.

A long non-coding RNA functions as a competitive endogenous RNA to modulate TaNAC018 by acting as a decoy for tae-miR6206

Increasing evidence indicates a strong correlation between the deposition of cuticular waxes and drought tolerance. However, the precise regulatory mechanism remains elusive. Here, we conducted a comprehensive transcriptome analysis of two wheat (Triticum aestivum) near-isogenic lines, the glaucous line G-JM38 rich in cuticular waxes and the non-glaucous line NG-JM31. We identified 85,143 protein-coding mRNAs, 4,485 lncRNAs, and 1,130 miRNAs. Using the lncRNA-miRNA-mRNA network and endogenous target mimic (eTM) prediction, we discovered that lncRNA35557 acted as an eTM for the miRNA tae-miR6206, effectively preventing tae-miR6206 from cleaving the NAC transcription factor gene TaNAC018. This lncRNA-miRNA interaction led to higher transcript abundance for TaNAC018 and enhanced drought-stress tolerance. Additionally, treatment with mannitol and abscisic acid (ABA) each influenced the levels of tae-miR6206, lncRNA35557, and TaNAC018 transcript. The ectopic expression of TaNAC018 in Arabidopsis also improved tolerance toward mannitol and ABA treatment, whereas knocking down TaNAC018 transcript levels via virus-induced gene silencing in wheat rendered seedlings more sensitive to mannitol stress. Our results indicate that lncRNA35557 functions as a competing endogenous RNA to modulate TaNAC018 expression by acting as a decoy target for tae-miR6206 in glaucous wheat, suggesting that non-coding RNA has important roles in the regulatory mechanisms responsible for wheat stress tolerance.

hsa_circ_0087100/hsa-miR-6743-5p affects Th1 cell differentiation by regulating STAT1 in diabetic retinopathy

Objective: To elucidate the role of the competitive endogenous RNA (ceRNA) network in immune infiltration of diabetic retinopathy (DR). Methods: We obtained differentially expressed (DE) circRNAs, miRNAs and mRNAs from the Gene Expression Omnibus database. Then, we identified immune infiltration by CIBERSORT and single-sample gene set enrichment analysis and discovered co-expression genes by weighted gene co-expression network analysis. Furthermore, STAT1-mediated Th1 differentiation was determined in DR cell models, DR patients and DR mouse models. Results: hsa_circ_0087100/hsa-miR-6743-5p/STAT1 was involved in immune infiltration of Th1 cells. Aberrant expression of the ceRNA network and STAT1-mediated Th1 differentiation was thus verified in vitro and in vivo. Conclusion: hsa_circ_0087100/hsa-miR-6743-5p/STAT1 may affect Th1 cell differentiation in DR.

The SLC19A1-AS/miR-1343/WNT11 axis is a novel positive regulatory ceRNA network governing goat granulosa cell proliferation

Long noncoding RNAs (lncRNAs), as competitive endogenous RNAs (ceRNAs), can directly or indirectly affect the proliferation and apoptosis of granulosa cells by regulating microRNA (miRNA) pathways. A ceRNA network of the SLC19A1-AS-miR-1343-WNT11 axis was constructed via comprehensive transcriptome sequencing of ovaries from goats with various fertility levels to further elucidate the function and regulatory mechanism of SLC19A1-AS in modulating miR-1343 and WNT11 during granulosa cell proliferation and apoptosis. Subsequent validation experiments were conducted in vitro using granulosa cells. In these experiments, we performed RNA immunoprecipitation (RIP) and identified SLC19A1-AS as a ceRNA in goat granulosa cells that promoted proliferation. Through bioinformatics prediction, luciferase reporter gene assays, and RNA pulldown assays, we confirmed that SLC19A1-AS acts as a sponge for miR-1343, preventing its binding to WNT11 mRNA and thereby increasing the expression of WNT11. This interaction also influenced the proliferation and apoptosis of granulosa cells. Our study systematically validated the biological function of the lncRNA-miRNA-mRNA ceRNA network in goat ovaries and revealed the potential regulatory mechanism by which SLC19A1-AS functions as a ceRNA in granulosa cells. These findings are expected to provide an important experimental foundation for further elucidating the physiological regulatory network of the ovary and contributing to reproductive health in goats.

The effect Of vascular related CeRNA genes and corresponding imaging biomarkers on survival in lower grade glioma

BACKGROUND & AIMS

To investigate the differential expression of vascular related ceRNA regulatory genes in LGG with different mutations of IDH1 and MGMT, and to verify imaging gene markers that can be closely associated with vascular related ceRNA regulatory genes.

METHOD

Five hundred fifteen patients with LGG were collected from TCGA database. CeRNA network analysis, GO analysis and Cox risk regression were used to find vascular ceRNA regulatory genes and their genetic markers related to survival. The preoperative MRI image data and postoperative tumor tissues of 14 patients with WHO grade III glioma were collected for full transcriptome analysis. The correlation between image characteristics of LGG and survival related vascular ceRNA regulatory genes was compared using nonparametric U test and Pearson correlation coefficient analysis.

RESULTS

Vascular related genes ranked first in the functional enrichment analysis of differentially expressed genes in LGG. EPHA2, ETS1, YAP1 and MEIS1 could significantly affect the survival of patients in each group of LGG. The volume of enhanced region was negatively correlated with IDH1 (r = -0.622, P = 0.009) mutation and TMEM100 (r = -0.535, P = 0.024), and positively correlated with MEIS1 (r = 0.551, P = 0.021), rCBFmax value was negatively correlated with TMEM100 (r = -0.492, P = 0.037).

CONCLUSIONS

Under different IDH1 mutations, lncRNA-dominated vascular-related ceRNA regulatory genes were the first differentially expressed subset of each group, and could be used as an effective risk factor affecting the survival of LGG. The image characteristics of LGG was an ideal image gene marker. It was a reliable imaging biological marker which can truly reflect the pathophysiological characteristics of glioma.

CircRNA expression profiles and regulatory networks in the vitreous humor of people with high myopia

Myopia is a global health and economic issue. Circular RNAs (circRNAs) have been shown to play an important role in the pathogenesis of many ocular diseases. We first evaluated the circRNA profiles and possible roles in vitreous humor samples of individuals with high myopia by a competitive endogenous RNA (ceRNA) array. Vitreous humor samples were collected from 15 high myopic (5 for ceRNA array, and 10 for qPCR) and 15 control eyes (5 for ceRNA array, and 10 for qPCR) with idiopathic epiretinal membrane (ERM) and macular hole (MH). 486 circRNAs (339 upregulated and 147 downregulated) and 264 mRNAs (202 upregulated and 62 downregulated) were differentially expressed between the high myopia and control groups. The expression of hsa_circ_0033079 (hsa-circDicer1), hsa_circ_0029989 (hsa-circNbea), hsa_circ_0019072 (hsa-circPank1) and hsa_circ_0089716 (hsa-circEhmt1) were validated by qPCR. Pearson analysis and multivariate regression analysis showed positive and significant correlations for axial length with hsa-circNbea and hsa-circPank1. KEGG analysis showed that the target genes of circRNAs were enriched in the mTOR, insulin, cAMP, and VEGF signaling pathways. GO analysis indicated that circRNAs mainly targeted transcription, cytoplasm, and protein binding. CircRNA-associated ceRNA network analysis and PPI network analysis identified several critical genes for myopia. The expression of circNbea, circPank1, miR-145-5p, miR-204-5p, Nras, Itpr1 were validated by qPCR in the sclera of form-deprivation myopia (FDM) mice model. CircPank1/miR-145-5p/NRAS and circNbea/miR-204-5p/ITPR1 were identified and may be important in the progression of myopia. Our findings suggest that circRNAs may contribute to the pathogenesis of myopia and may serve as potential biomarkers.

Unveiling the role of YARS1 in bladder cancer: A prognostic biomarker and therapeutic target

YARS is responsible for catalysing the binding of tyrosine to its cognate tRNA and plays a crucial role in basic biosynthesis. However, its biological functions in bladder cancer remains to be proven. We analysed variations in YARS1 expression and survival in bladder cancer using multiple data sets, including TCGA-BLCA, GSE13507 and bladder cancer-specific tissue microarrays. Furthermore, we explored the biological functions of YARS1 using transcriptome data. Our findings revealed a noteworthy correlation between YARS1 and immune infiltration in bladder cancer, as determined using the XCELL algorithm and single-cell analysis. In addition, we employed the TIDE algorithm to evaluate the responsiveness of different cohorts to immune checkpoint therapy. We investigated the regulatory associations between YARS1 and various aspects of bladder cancer, including senescence, ferroptosis and stemness. Finally, we established a ceRNA network that is directly linked to the overall prognosis, YARS1 can serve as a prognostic biomarker for bladder cancer; its interaction with MYC has implications for bladder cancer cell senescence, ferroptosis and stemness. Moreover, the identified ceRNA network has potential as a therapeutic target in bladder cancer.

Competitive endogenous RNA networks: Decoding the role of long non-coding RNAs and circular RNAs in colorectal cancer chemoresistance

Colorectal cancer (CRC) is recognized as one of the most common gastrointestinal malignancies across the globe. Despite significant progress in designing novel treatments for CRC, there is a pressing need for more effective therapeutic approaches. Unfortunately, many patients undergoing chemotherapy develop drug resistance, posing a significant challenge for cancer treatment. Non-coding RNAs (ncRNAs) have been found to play crucial roles in CRC development and its response to chemotherapy. However, there are still gaps in our understanding of interactions among various ncRNAs, such as long non-coding RNAs (lncRNAs), circular RNAs (circRNAs) and microRNAs (miRNAs). These ncRNAs can act as either oncogenes or tumour suppressors, affecting numerous biological functions in different cancers including CRC. A class of ncRNA molecules known as competitive endogenous RNAs (ceRNAs) has emerged as a key player in various cellular processes. These molecules form networks through lncRNA/miRNA/mRNA and circRNA/miRNA/mRNA interactions. In CRC, dysregulation of ceRNA networks has been observed across various cellular processes, including proliferation, apoptosis and angiogenesis. These dysregulations are believed to play a significant role in the progression of CRC and, in certain instances, may contribute to the development of chemoresistance. Enriching our knowledge of these dysregulations holds promise for advancing the field of diagnostic and therapeutic modalities for CRC. In this review, we discuss lncRNA- and circRNA-associated ceRNA networks implicated in the emergence and advancement of drug resistance in colorectal carcinogenesis.

Mining key circRNA-associated-ceRNA networks for milk fat metabolism in cows with varying milk fat percentages

BACKGROUND

Cow milk fat is an essential indicator for evaluating and measuring milk quality and cow performance. Growing research has identified the molecular functions of circular RNAs (circRNAs) necessary for mammary gland development and lactation in mammals.

METHOD

The present study analyzed circRNA expression profiling data in mammary epithelial cells (MECs) from cows with highly variable milk fat percentage (MFP) using differential expression analysis and weighted gene co-expression network analysis (WGCNA).

RESULTS

A total of 309 differentially expressed circRNAs (DE-circRNAs) were identified in the high and low MFP groups. WGCNA analysis revealed that the pink module was significantly associated with MFP (r = - 0.85, P = 0.007). Parental genes of circRNAs in this module were enriched mainly in lipid metabolism-related signaling pathways, such as focal adhesion, ECM-receptor interaction, adherens junction and AMPK. Finally, six DE-circRNAs were screened from the pink module: circ_0010571, circ_0007797, circ_0002746, circ_0003052, circ_0004319, and circ_0012840. Among them, circ_0002746, circ_0003052, circ_0004319, and circ_0012840 had circular structures and were highly expressed in mammary tissues. Subcellular localization revealed that these four DE-circRNAs may play a regulatory role in the mammary glands of dairy cows, mainly as competitive endogenous RNAs (ceRNAs). Seven hub target genes (GNB1, GNG2, PLCB1, PLCG1, ATP6V0C, NDUFS4, and PIGH) were obtained by constructing the regulatory network of their ceRNAs and then analyzed by CytoHubba and MCODE plugins in Cytoscape. Functional enrichment analysis revealed that these genes are crucial and most probable ceRNA regulators in milk fat metabolism.

CONCLUSIONS

Our study identified several vital circRNAs and ceRNAs affecting milk fat synthesis, providing new research ideas and a theoretical basis for cow lactation, milk quality, and breed improvement.

Integrated analysis of circRNA- related ceRNA network targeting neuroinflammation in medial temporal lobe epilepsy

BACKGROUND

medial temporal lobe epilepsy (mTLE) is among the most common types of temporal lobe epilepsy (TLE) ，it is generally resistant to drug treatment, which significantly impacts the quality of life and treatment. Research on novel therapeutic approaches for mTLE has become a current focus. Our study aims to construct and analyze a competing endogenous RNA (ceRNA) network that targets neuroinflammation using publicly available data, which may offer a novel therapeutic approach for mTLE.

METHODS

we utilized the R package to analyze GSE186334 downloaded from Gene Expression Omnibus database, subsequently constructing and identifying hub network within the ceRNA network using public databases. Lastly, we validated the expressions and interactions of some nodes within the hub ceRNA network in Sombati cell model.

RESULTS

our transcriptome analysis identified 649 differentially expressed (DE) mRNAs (273 up-regulated, 376 down-regulated) and 36 DE circRNAs (11 up-regulated, 25 down-regulated) among mTLE patients. A total of 23 candidate DE mRNAs associated with neuroinflammation were screened, and two ceRNA networks were constructed. A hub network was further screened which included 3 mRNAs, 22 miRNAs, and 11 circRNAs. Finally, we confirmed the hsa-miR-149-5p is crucial in the regulatory effect of hsa_circ_0005145 on IL - 1α in the hub network.

CONCLUSIONS

In summary, our study identified a hub ceRNA network and validated a potential circRNA-miRNA-mRNA axis targeting neuroinflammation. The results of our research may serve as a potential therapeutic target for mTLE.

Analysis of ceRNA Network and Identification of Potential Treatment Target and Biomarkers of Endothelial Cell Injury in Sepsis

Background: Sepsis is a complex clinical syndrome caused by a dysregulated host immune response to infection. This study aimed to identify a competing endogenous RNA (ceRNA) network that can greatly contribute to understanding the pathophysiological process of sepsis and determining sepsis biomarkers. Methods: The GSE100159, GSE65682, GSE167363, and GSE94717 datasets were obtained from the Gene Expression Omnibus (GEO) database. Weighted gene coexpression network analysis was performed to find modules possibly involved in sepsis. A long noncoding RNA-microRNA-messenger RNA (lncRNA-miRNA-mRNA) network was constructed based on the findings. Single-cell analysis was performed. Human umbilical vein endothelial cells were treated with lipopolysaccharide (LPS) to create an in vitro model of sepsis for network verification. Reverse transcription-polymerase chain reaction, fluorescence in situ hybridization, and luciferase reporter genes were used to verify the bioinformatic analysis. Result: By integrating data from three GEO datasets, we successfully constructed a ceRNA network containing 18 lncRNAs, 7 miRNAs, and 94 mRNAs based on the ceRNA hypothesis. The lncRNA ZFAS1 was found to be highly expressed in LPS-stimulated endothelial cells and may thus play a role in endothelial cell injury. Univariate and multivariate Cox analyses showed that only SLC26A6 was an independent predictor of prognosis in sepsis. Overall, our findings indicated that the ZFAS1/hsa-miR-449c-5p/SLC26A6 ceRNA regulatory axis may play a role in the progression of sepsis. Conclusion: The sepsis ceRNA network, especially the ZFAS1/hsa-miR-449c-5p/SLC26A6 regulatory axis, is expected to reveal potential biomarkers and therapeutic targets for sepsis management.

Whole-transcriptome sequencing reveals a melanin-related ceRNA regulatory network in the breast muscle of Xichuan black-bone chicken

The economic losses incurred due to reduced muscle pigmentation highlight the crucial role of melanin-based coloration in the meat of black-bone chickens. Melanogenesis in the breast muscle of black-bone chickens is currently poorly understood in terms of molecular mechanisms. This study employed whole-transcriptome sequencing to analyze black and white breast muscle samples from black-bone chickens, leading to the identification of 367 differentially expressed (DE) mRNAs, 48 DElncRNAs, 104 DEcircRNAs, and 112 DEmiRNAs involved in melanin deposition. Based on these findings, a competitive endogenous RNA (ceRNA) network was developed to better understand the complex mechanisms of melanin deposition. Furthermore, our analysis revealed key DEmRNAs (TYR, DCT, EDNRB, MLPH and OCA2) regulated by DEmiRNAs (gga-miR-140-5p, gga-miR-1682, gga-miR-3529, gga-miR-499-3p, novel-m0012-3p, gga-miR-200b-5p, gga-miR-203a, gga-miR-6651-5p, gga-miR-7455-3p, gga-miR-31-5p, miR-140-x, miR-455-x, novel-m0065-3p, gga-miR-29b-1-5p, miR-455-y, novel-m0085-3p, and gga-miR-196-1-3p). These DEmiRNAs competitively interacted with DElncRNAs including MSTRG.2609.2, MSTRG.4185.1, LOC112530666, LOC112533366, LOC771030, LOC107054724, LOC121107411, LOC100859072, LOC101750037, LOC121108550, LOC121109224, LOC121110876, and LOC101749016, as well as DEcircRNAs, such as novel_circ_000158, novel_circ_000623, novel_001518, and novel_circ_003596. The findings from this study provide insight into the mechanisms that regulate lncRNA, circRNA, miRNA, and mRNA expression in chicken melanin deposition.

Role of circular RNA/miRNA axes in the pathophysiology of polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is a disorder resulted from interactions between genetic and environmental factors. Based on the importance of epigenetic factors in the pathoetiology of PCOS, the current review focused on identification of circular RNAs (circRNAs) that are involved in PCOS through acting as molecular sponges for microRNAs (miRNAs). The literature search led to identification of circ_0043533/miR-1179, circ_0030018/miR-136, circ_FURIN/miR-423-5p, circ-FURIN/miR-195-5p, circ_0043532/miR-182, circ_RANBP9/miR-136-5p, circRHBG/miR-515-5p, circMTO1/miR-320b, circASPH/miR-375, circPSMC3/miR-296-3p, circLDLR/miR-1294, circPUM1/miR-760, and hsa_circ_0118530/miR-136 as molecular axes contributing to the pathogenesis of PCOS. To set the stage for future research on the role of the ceRNA network in PCOS, in-silico analyses were performed using miRWalk, miRNet, and miRDIP databases. miRWalk identified 80 genes regulated by 5 miRNAs, miRNet revealed 6449 circRNAs potentially controlling 11 miRNAs, and miRDIP identified 11 miRNAs associated with 35 human pathways. These targets can be used in the treatment options, design of personalized medicine and prediction of prognosis of PCOS.

Immune dysfunction mediated by the competitive endogenous RNA network in fetal side placental tissue of polycystic ovary syndrome

Polycystic ovary syndrome (PCOS), a common endocrine and metabolic disorder affecting women in their reproductive years. Emerging evidence suggests that the maternal-fetal immune system is crucial for proper pregnancy. However, whether immune function is altered at the end of pregnancy in PCOS women and the underlying molecular mechanisms is currently unexplored. Herein, the basic maternal immune system was investigated (n = 136 in the control group; n = 103 in the PCOS group), and whole-transcriptome sequencing was carried out to quantify the mRNAs, miRNAs, and lncRNAs expression levels in fetal side placental tissue of women with PCOS. GO, KEGG, and GSEA analysis were employed for functional enrichment analysis. The process of identifying hub genes was conducted utilizing the protein-protein interaction network. CIBERSORT and Connectivity Map were deployed to determine immune cell infiltration and predict potential drugs, respectively. A network of mRNA-miRNA-lncRNA was constructed and then validated by qRT-PCR. First, red blood cell count, neutrophil count, lymphocyte count, hypersensitive C-reactive protein, and procalcitonin were significantly elevated, while placental growth factor was hindered in PCOS women. We identified 308 DEmRNAs, 77 DEmiRNAs, and 332 DElncRNAs in PCOS samples. Functional enrichment analysis revealed that there were significant changes observed in terms of the immune system, especially the chemokine pathway. Eight genes, including FOS, JUN, EGR1, CXCL10, CXCR1, CXCR2, CXCL11, and CXCL8, were considered as hub genes. Furthermore, the degree of infiltration of neutrophils was dramatically decreased in PCOS tissues. In total, 57 ceRNA events were finally obtained, and immune-related ceRNA networks were validated. Some potential drug candidates, such as enalapril and RS-100329, could have a function in PCOS therapy. This study represents the inaugural attempt to evaluate the immune system at the end of pregnancy and placental ceRNA networks in PCOS, indicating alterations in the chemokine pathway, which may impact fetal and placental growth, and provides new therapy targets.

Maternal nicotine exposure promotes hippocampal CeRNA-mediated excitotoxicity and social barriers in adolescent offspring mice

Nicotine, an addictive component of cigarettes, causes cognitive defects, particularly when exposure occurs early in life. However, the exact mechanism through which nicotine causes toxicity and alters synaptic plasticity is still not fully understood. The aim of the current study is to examine how non-coding developmental regulatory RNA impacts the hippocampus of mice offspring whose mothers were exposed to nicotine. Female C57BL/6J mice were given nicotine water from one week before pregnancy until end of lactation. Hippocampal tissue from offspring at 20 days post-birth was used for LncRNA and mRNA microarray analysis. Differential expression of LncRNAs and mRNAs associated with neuronal development were screened and validated, and the CeRNA pathway mediating neuronal synaptic plasticity GM13530/miR-7119-3p/mef2c was predicted using LncBase Predicted v.2. Using protein immunoblotting, Golgi staining and behavioral tests, our findings revealed that nicotine exposure in offspring mice increased hippocampal NMDAR receptor, activated receptor-dependent calcium channels, enhanced the formation of NMDAR/nNOS/PSD95 ternary complexes, increased NO synthesis, mediated p38 activation, induced neuronal excitability toxicity. Furthermore, an epigenetic CeRNA regulatory mechanism was identified, which suppresses Mef2c-mediated synaptic plasticity and leads to modifications in the learning and social behavior of the offspring during adolescence. This study uncovers the way in which maternal nicotine exposure results in neurotoxicity in offspring.

Discovering therapeutic possibilities for polycystic ovary syndrome by targeting XIST and its associated ceRNA network through the analysis of transcriptome data

Long non-coding RNA (lncRNA) regulates many physiological processes by acting as competitive endogenous RNA (ceRNA). The dysregulation of lncRNA X-inactive specific transcript (XIST) has been shown in various human disorders. However, its role in the pathogenesis of polycystic ovary syndrome (PCOS) is yet to be explored. This study aimed to explore the underlying mechanism of XIST in the pathogenesis of PCOS, specifically through dataset functional analysis. GEO PCOS datasets including RNA-seq, microarray, and miRNA-seq in granulosa cells (GCs) and blood, were examined and comprehensively analyzed. Enrichment analysis, ROC curve constructions, lncRNA-miRNA-mRNA interaction network analyses, and qRT-PCR validation were performed followed by a series of drug signature screenings. Our results revealed significant dysregulation in the expression of 1131 mRNAs, 30 miRNAs, and XIST in GCs of PCOS patients compared to healthy individuals. Of the120 XIST-correlated upregulated genes, 25 were enriched in inflammation-related pathways. Additionally, 5 miRNAs were identified as negative regulators of XIST-correlated genes. Accordingly, a ceRNA network containing XIST-miRNAs-mRNAs interactions was constructed. Furthermore, 6 genes, including AQP9, ETS2, PLAU, PLEK, SOCS3, and TNFRSF1B served as both GCs and blood-based biomarkers. By analyzing the number of interactions among XIST, miRNAs, and mRNAs, we pinpointed ETS2 as the pivotal gene within the ceRNA network. Our findings reveal a novel XIST- hsa-miR-146a-5p, hsa-miR-144-3p, and hsa-miR-1271-5p-ETS2 axis that comprehensively elucidates the XIST-associated mechanism underlying PCOS onset. qRT-PCR analysis further confirmed the, overexpression of both XIST and ETS2 . Furthermore, our results demonstrated that XIST and ETS2 were correlated with some assisted reproductive technologies outcomes. Finally, we identified two novel compounds including, methotrexate/folate and threonine using drug-gene interaction databases for PCOS management. These findings provide novel insights into the molecular etiology, diagnosis, and potential therapeutic interventions for PCOS.

Comprehensive exosomal microRNA profile and construction of competing endogenous RNA network in autism spectrum disorder: A pilot study

Exosomes have been demonstrated to exert momentous roles in autism spectrum disorder (ASD). However, few studies have reported a correlation between exosomal microRNAs (miRNAs) and ASD. To date, our understanding of crucial competing endogenous RNA (ceRNA) networks in ASD remains limited. Herein, the exosomal miRNA profile in the peripheral blood of children with ASD and healthy controls was investigated and the level of immune cell infiltration in ASD was evaluated to determine the distribution of immune cell subtypes. Exosomes were isolated from the peripheral blood of ten children with ASD and ten healthy controls, and further identified using transmission electron microscopy and western blot analysis. RNA sequencing was conducted to investigate exosomal miRNA profiles in patients with ASD. The mRNA and circular RNA (circRNA) expression profiles were acquired from the Gene Expression Omnibus (GEO) database. Differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs), and circRNAs (DEcircRNAs) were identified and ceRNA regulatory networks were constructed. Furthermore, the immune cell infiltration levels in patients with ASD were evaluated. Exosomes were spherical, approximately 100 nm in size, and were confirmed via western blot analysis using exosome-associated markers CD9, CD63, and CD81. Thirty-five DEmRNAs, 63 DEmiRNAs, and 494 DEcircRNAs were identified in patients with ASD. CeRNA regulatory networks, including 6 DEmRNAs, 14 DEmiRNAs, and 86 DEcircRNAs, were established. Correlation analysis indicated that leucine-rich glioma inactivated protein 1 (LGI1) expression was significantly positively correlated with the content of CD8+ T cells. Our findings may be conducive to offering novel insights into this disease and providing further evidence of transcriptomic abnormalities in ASD.

Construction of LncRNA-mediated CeRNA network for investigating the immune pathogenesis of myocardial infarction

BACKGROUND

Myocardial infarction (MI) is a cardiovascular disease that seriously threatens human health. However, an immune-related competitive endogenous RNA (ceRNA) network has not been reported in MI.

METHODS

The GSE66360, GSE19339, GSE97320, GSE61741, and GSE168281 datasets were acquired from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) and differentially expressed miRNAs (DEmiRNAs) from MI patients and healthy controls were screened and an immune-related ceRNA network was constructed. Furthermore, the key long noncoding RNAs(lncRNAs) highly related to the immune mechanism of MI were identified utilizing the random walk with restart algorithm. Finally, the expression of the hub genes was further verified in the GSE66360, GSE19339, and GSE97320 datasets, and quantitativereal-time polymerase chain reaction (qRT-PCR) was performed for the MI patients and healthy controls.

RESULTS

A total of 184 differentially expressed immune-related genes(DE-IRGs) and 432 DE-miRNAs were obtained, and an immune-related ceRNA network comprising 1421 lncRNAs, 61 DE-miRNAs, and 139 DE-IRGs was constructed. According to the order of stress, betweenness, and closeness, NEAT1, KCNQ1OT1, and XIST were identified as key lncRNAs. Moreover, random walk with restart analysis also suggested that NEAT1, KCNQ1OT1, and XIST are key lncRNAs. Subsequently, a ceRNA network of 10 hub genes and 3 lncRNAs was constructed. Finally, we found that the expression of FCER1G and TYROBP significantly differed between MI patients and control individuals in the GSE66360, GSE19339, and GSE97320 datasets. qRT-PCR revealed that the expression of NEAT1, KCNQ1OT1, XIST, FCER1G, and TYROBP was significantly elevated in MI tissue samples compared to healthy control tissue samples.

CONCLUSION

NEAT1, KCNQ1OT1, XIST, FCER1G, and TYROBP are involved in MI and can be used as molecular biomarkers for the screening and diagnosis of MI. Furthermore, the immune system plays an essential role in the onset and progression of MI.

TGF-β-activated circRYK drives glioblastoma progression by increasing VLDLR mRNA expression and stability in a ceRNA- and RBP-dependent manner

BACKGROUND

The TGF-β signalling pathway is intricately associated with the progression of glioblastoma (GBM). The objective of this study was to examine the role of circRNAs in the TGF-β signalling pathway.

METHODS

In our research, we used transcriptome analysis to search for circRNAs that were activated by TGF-β. After confirming the expression pattern of the selected circRYK, we carried out in vitro and in vivo cell function assays. The underlying mechanisms were analysed via RNA pull-down, luciferase reporter, and RNA immunoprecipitation assays.

RESULTS

CircRYK expression was markedly elevated in GBM, and this phenotype was strongly associated with a poor prognosis. Functionally, circRYK promotes epithelial-mesenchymal transition and GSC maintenance in GBM. Mechanistically, circRYK sponges miR-330-5p and promotes the expression of the oncogene VLDLR. In addition, circRYK could enhance the stability of VLDLR mRNA via the RNA-binding protein HuR.

CONCLUSION

Our findings show that TGF-β promotes epithelial-mesenchymal transition and GSC maintenance in GBM through the circRYK-VLDLR axis, which may provide a new therapeutic target for the treatment of GBM.

Circular RNA circ-SLC7A5 Functions as a Competing Endogenous RNA to Impact Cell Biological Behaviors in Esophageal Squamous Cell Carcinoma (ESCC)

BACKGROUND

Circular RNAs (circRNAs) have profound effects on establishment and pathogenesis of esophageal squamous cell carcinoma (ESCC). Here, we defined whether circRNA solute carrier family 7 member 5 (circ-SLC7A5, also called hsa_circ_0040796) is causally involved in the pathogenesis of ESCC.

METHODS

Circ-SLC7A5, microRNA (miR)-874-3p and coronin-1C (CORO1C) expression levels were gauged by qRT-PCR or immunoblotting. Cell functional phenotypes were tested by colony formation, EdU, flow cytometry, transwell and wound-healing assays. RNA immunoprecipitation (RIP) and dual-luciferase reporter assays were applied to ascertained circ-SLC7A5/miR-874-3p and miR-874-3p/CORO1C relationships.

RESULTS

Circ-SLC7A5 was highly expressed in human ESCC. Circ-SLC7A5 depletion impaired cell growth, migration, invasiveness, and promoted apoptosis. Circ-SLC7A5 knockdown diminished ESCC cell tumorigenicity. Mechanistically, circ-SLC7A5 contained a binding site for miR-874-3p. Also, miR-874-3p was responsible for circ-SLC7A5's function in ESCC cells. CORO1C was a direct miR-874-3p target. Circ-SLC7A5 functioned as a competing endogenous RNA (ceRNA) to control CORO1C by competing for shared miR-874-3p. Furthermore, CORO1C knockdown phenocopied miR-874-3p overexpression in impacting the biological behaviors of ESCC cells.

CONCLUSION

These findings identify circ-SLC7A5 as a crucial modulator of ESCC cells and establish a novel circ-SLC7A5/miR-874-3p/CORO1C ceRNA network in ESCC.

Effect of runx2b deficiency in intermuscular bones on the regulatory network of lncRNA-miRNA-mRNA

Intermuscular bones (IBs) are mineralized spicules that negatively impact the quality and value of fish products. Runx2b is a crucial modulator in promoting bone formation through regulating osteoblast differentiation. Previous studies suggested that loss of runx2b gene completely inhibited IBs formation in zebrafish. However, how the whole transcriptome, including mRNA and non-coding RNA (ncRNA), affects the IBs development in runx2b-/- zebrafish are not known. The aim of this study was to identify the regulatory networks of differentially expressed (DE) lncRNAs, miRNAs, and mRNAs in zebrafish with and without IBs (runx2b+/+ fish and runx2b-/- fish) utilizing high-throughput sequencing techniques. All together there are 1051 mRNAs, 456 lncRNAs, and 18 miRNAs differentially expressed were found between these two strains. The analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) has highlighted significant pathways linked to the development of IBs, specifically the TGF-beta and Wnt signaling pathways, and a number of genes concentrated on these two signaling pathways related to the formation of IBs. Further, 1989 competing endogenous RNA (ceRNA) networks were created according to the correlation among mRNAs, miRNAs and lncRNAs. The ceRNA networks results revealed 52 ceRNA pairs related to the IBs formation, consisting of 52 mRNAs, 37 lncRNAs, and 6 miRNAs. Of these, we found that dre-miR-2189 was the key element of ceRNA pairs, interacting with 19 mRNAs and 11 lncRNAs, and MSTRG.13175.1 could regulate sp7 expression by interacting with dre-miR-2189 to function in osteogenic differentiation. Subsequent experiments at the cellular level also revealed the interaction mechanism. The outcomes indicated a crucial role of miRNAs and lncRNAs in the development of fish IBs, which offer new views into the functions of ncRNAs involved in IBs formation.

METTL3 as a novel diagnosis and treatment biomarker and its association with glycolysis, cuproptosis and ceRNA in oesophageal carcinoma

METTL3 has been shown to be involved in regulating a variety of biological processes. However, the relationship between METTL3 expression and glycolysis, cuproptosis-related genes and the ceRNA network in oesophageal carcinoma (ESCA) remains unclear. ESCA expression profiles from databases were obtained, and target genes were identified using differential analysis and visualization. Immunohistochemistry (IHC) staining assessed METTL3 expression differences. Functional enrichment analysis using GO, KEGG and GSEA was conducted on the co-expression profile of METTL3. Cell experiments were performed to assess the effect of METTL3 interference on tumour cells. Correlation and differential analyses were carried out to assess the relationship between METTL3 with glycolysis and cuproptosis. qRT-PCR was used to validate the effects of METTL3 interference on glycolysis-related genes. Online tools were utilized to screen and construct ceRNA networks based on the ceRNA theory. METTL3 expression was significantly higher in ESCA compared to the controls. The IHC results were consistent with the above results. Enrichment analysis revealed that METTL3 is involved in multiple pathways associated with tumour development. Significant correlations were observed between METTL3 and glycolysis-related genes and cuproptosis-related gene. Experiments confirmed that interfered with METTL3 significantly inhibited glucose uptake and lactate production in tumour cells, and affected the expression of glycolytic-related genes. Finally, two potential ceRNA networks were successfully predicted and constructed. Our study establishes the association between METTL3 overexpression and ESCA progression. Additionally, we propose potential links between METTL3 and glycolysis, cuproptosis and ceRNA, presenting a novel targeted therapy strategy for ESCA.

Bovine enhancer-regulated circSGCB acts as a ceRNA to regulate skeletal muscle development via enhancing KLF3 expression

Skeletal muscle growth and development in livestock and poultry play a pivotal role in determining the quality and yield of meat production. However, the mechanisms of myogenesis are remained unclear due to it finely regulated by a complex network of biological macromolecules. In this study, leveraging previous sequencing data, we investigated a differentially expressed circular RNA (circSGCB) present in fetal and adult muscle tissues among various ruminant species, including cattle, goat, and sheep. Our analysis revealed that circSGCB is a single exon circRNA, potentially regulated by an adjacent bovine enhancer. Functional analysis through loss-of-function tests demonstrated that circSGCB exerts inhibitory effects on bovine myoblast proliferation while promoting myocytes generation. Furthermore, we discovered that circSGCB primarily localizes to the cytoplasm, where it functions as a molecular sponge by binding to bta-miR-27a-3p. This interaction releases the mRNAs of KLF3 gene and further activates downstream functional pathways. In vivo, studies provided evidence that up-regulation of KLF3 contributes to muscle regeneration. These findings collectively suggest that circSGCB operates via a competing endogenous RNA (ceRNA) mechanism to regulate KLF3, thereby influencing myogenesis in ruminants and highlights it may as potential molecular targets for enhancing meat production in livestock and poultry industries.

Screening of hair follicle telogen-associated circRNAs in sheep and construction of their ceRNA network

Sheep breeds with hair-shedding traits have many advantages over non-shedding sheep breeds, not only because of reduced shearing labor and feeding management costs but also because it reduces in vitro parasites and improves adaptability to summer heat stress. The wool of Dorper sheep naturally sheds in spring due to the periodic growth of hair follicles. CircRNAs primarily regulate the morphogenesis of hair follicles through the ceRNA mechanism. In this study, five 2-year-old Dorper ewes with extreme hair-shedding phenotype (S) and three Dorper ewes with non-shedding (N) phenotype were selected for subsequent analyses. For RNA extraction, skin tissues were collected on 27th September 2019 (S1, N1), 3rd January 2020 (S2, N2), and 17th March 2020 (S3, N3), which were then subjected to RNA-seq. RNA-seq technology revealed 20,185 novel circRNAs in the hair follicles of Dorper sheep. Among them, 1450 circRNAs were differentially expressed (DE). Clustering heatmap and expression pattern analyses were performed on DE circRNAs, which indicated 78 circRNAs with T pattern (Telogen, highly expressed in telogen), and the source genes for candidate circRNAs were further screened by functional enrichment analysis, which identified 13 crucial genes enriched in pathways associated with hair follicle development. Additionally, a ceRNA regulatory network comprising 4 circRNAs, 11 miRNAs, and 13 target genes was constructed. Overall, this study screened circRNAs that may be associated with the telogen phase of hair follicles in sheep, providing a relevant theoretical basis for wool shedding in sheep and for breeding Dorper sheep with automatic wool shedding.

Long noncoding RNA XIST inhibition promotes Leydig cell apoptosis by acting as a competing endogenous RNA for microRNA-145a-5p that targets SIRT1 in late-onset hypogonadism

Leydig cell (LCs) apoptosis is responsible for decreased serum testosterone levels during late-onset hypogonadism (LOH). Our study was designed to illustrate the regulatory effect of lncRNA XIST on LCs and to clarify its molecular mechanism of action in LOH. The Leydig cells (TM3) was treated by 300 μM H2O2 for 8 h to establish Leydig cell oxidative stress model in vitro. The expression levels of lncRNA XIST in the testicular tissues of patients with LOH were measured using fluorescence in situ hybridization (FISH). The interaction between lncRNA XIST/SIRT1 and miR-145a-5p was assessed using starBase and dual-luciferase reporter gene assays. Apoptotic cells and Caspase3 activity were determined by flow cytometry (FCM) assay. Testosterone concentration was determined by ELISA. Moreover, histological assessment of testicles in mice was performed by using HE staining and the TUNEL assay was used to determine apoptosis. We found that the lncRNA XIST was downregulated in the testicular tissues of LOH patients and mice and in H2O2-induced TM3 cells. XIST siRNA significantly promoted apoptosis, enhanced Caspase3 activity and reduced testosterone levels in H2O2-stimulated TM3 cells. Further studies showed that the miR-145a-5p inhibitor reversed the effect of XIST-siRNA on H2O2-induced Leydig cell apoptosis. MiR-145a-5p negatively regulated SIRT1 expression, and SIRT1-siRNA reversed the effects of the miR-145a-5p inhibitor on H2O2 stimulated TM3 cells. The in vivo experiments indicated that silencing of the lncRNA XIST aggravated LOH symptoms in mice. Inhibition of lncRNA XIST induces Leydig cell apoptosis through the miR-145a-5p/SIRT1 axis in the progression of LOH.

Integrated analysis of ceRNA-miRNA changes in paraquat-induced pulmonary epithelial-mesenchymal transition via high-throughput sequencing

Recent studies have shown that epithelial-mesenchymal transition (EMT) plays an important role in paraquat (PQ)-induced tissue fibrosis, which is the main cause of death in patients with PQ poisoning. However, no effective treatment for pulmonary interstitial fibrosis caused by PQ poisoning exists. It is of great significance for us to find new therapeutic targets through bioinformatics in PQ-induced EMT. We conducted transcriptome sequencing to determine the expression profiles of 1210 messenger RNAs (mRNAs), 558 long noncoding RNAs, 28 microRNAs (miRNAs), including 18 known-miRNAs, 10 novel-miRNAs and 154 circular RNAs in the PQ-exposed EMT group mice. Using gene ontology and Kyoto Encyclopaedia of Genes and Genomes analyses, we identified the pathways associated with signal transduction, cancers, endocrine systems and immune systems were involved in PQ-induced EMT. Furthermore, we constructed long noncoding RNA-miRNA-mRNA interrelated networks and found that upregulated genes included Il22ra2, Mdm4, Slc35e2 and Angptl4, and downregulated genes included RGS2, Gabpb2, Acvr1, Prkd3, Sp100, Tlr12, Syt15 and Camk2d. Thirteen new potential competitive endogenous RNA targets were also identified for further treatment of PQ-induced pulmonary tissue fibrosis. Through further study of the pathway and networks, we may identify new molecular targets in PQ-induced pulmonary EMT.

Integrated analysis identifies cuproptosis-related gene DLAT and its competing endogenous RNAs network to predict the prognosis of pancreatic adenocarcinoma patients

Pancreatic adenocarcinoma (PAAD) is a highly malignant tumor with poor prognosis. However, the relationship between cuproptosis-related genes (CRGs) and its competing endogenous RNA (ceRNA) network with the prognosis of PAAD patients remains unclear. To investigate this relationship, we calculated the difference in CRGs between PAAD tissues and normal tissues using the 'limma' R package. Additionally, we employed least absolute shrinkage and selection operator (LASSO) Cox regression analysis to construct a prognostic signature for CRGs. Survival analysis of patients with PAAD was performed using Kaplan-Meier analysis. Furthermore, we used bioinformatics tools to screen for CRGs-related MicroRNA (miRNA) and lncRNAs. To validate these findings, we conducted real-time quantitative polymerase chain reaction (RT-qPCR), CCK-8, colony formation, and Transwell assays to assess the effect of DLAT in vitro. Our results revealed a cuproptosis-related prognostic signature consisting of 3 prognostic genes (DLAT, LIAS, and LIPT1). Notably, patients with a high-risk score for the CRGs signature exhibited poor prognosis in terms of overall survival (OS) (P < .05). The receiver operating characteristic (ROC) curve was used to evaluate the prognostic signature of CRGs. The results showed that the 1-year, 3-year, and 5-year area under the curve values for predicting OS were 0.62, 0.66, and 0.79, respectively. Additionally, the CRGs-related ceRNA network revealed the regulatory axis of LINC00857/has-miR-1179/DLAT in PAAD. In vitro experiments demonstrated that knockdown of LINC00857 and DLAT inhibited the growth and invasion of PAAD cells. This study identified a CRG-related prognostic signature consisting of 3 biomarkers (DLAT, LIAS, and LIPT1) for PAAD. Furthermore, ceRNA network analysis suggested the involvement of the LINC00857/has-miR-1179/DLAT axis in the development of PAAD. Overall, this study provides theoretical support for the investigation of diagnostic and prognostic biomarkers as well as potential therapeutic targets in PAAD.

Hijacking and rewiring of host CircRNA/miRNA/mRNA competitive endogenous RNA (ceRNA) regulatory networks by oncoviruses during development of viral cancers

A significant portion of human cancers are caused by oncoviruses (12%-25%). Oncoviruses employ various strategies to promote their replication and induce tumourigenesis in host cells, one of which involves modifying the gene expression patterns of the host cells, leading to the rewiring of genes and resulting in significant changes in cellular processes and signalling pathways. In recent studies, a specific mode of gene regulation known as circular RNA (circRNA)-mediated competing endogenous RNA (ceRNA) networks has emerged as a key player in this context. CircRNAs, a class of non-coding RNA molecules, can interact with other RNA molecules, such as mRNAs and microRNAs (miRNAs), through a process known as ceRNA crosstalk. This interaction occurs when circRNAs, acting as sponges, sequester miRNAs, thereby preventing them from binding to their target mRNAs and modulating their expression. By rewiring the host cell genome, oncoviruses have the ability to manipulate the expression and activity of circRNAs, thereby influencing the ceRNA networks that can profoundly impact cellular processes such as cell proliferation, differentiation, apoptosis, and immune responses. This review focuses on a comprehensive evaluation of the latest findings on the involvement of virus-induced reprogramming of host circRNA-mediated ceRNA networks in the development and pathophysiology of human viral cancers, including cervical cancer, gastric cancer, nasopharyngeal carcinoma, Kaposi's sarcoma, hepatocellular carcinoma, and diffuse large B cell lymphoma. Understanding these mechanisms can improve our knowledge of how oncoviruses contribute to human tumourigenesis and identify potential targets for developing optimised therapies and diagnostic tools for viral cancers.

Expression profiles of host miRNAs and circRNAs and ceRNA network during Toxoplasma gondii lytic cycle

Toxoplasma gondii is an opportunistic protozoan parasite that is highly prevalent in the human population and can lead to adverse health consequences in immunocompromised patients and pregnant women. Noncoding RNAs, such as microRNAs (miRNAs) and circular RNAs (circRNAs), play important regulatory roles in the pathogenesis of many infections. However, the differentially expressed (DE) miRNAs and circRNAs implicated in the host cell response during the lytic cycle of T. gondii are unknown. In this study, we profiled the expression of miRNAs and circRNAs in human foreskin fibroblasts (HFFs) at different time points after T. gondii infection using RNA sequencing (RNA-seq). We identified a total of 7, 7, 27, 45, 70, 148, 203, and 217 DEmiRNAs and 276, 355, 782, 1863, 1738, 6336, 1229, and 1680 DEcircRNAs at 1.5, 3, 6, 9, 12, 24, 36, and 48 h post infection (hpi), respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses revealed that the DE transcripts were enriched in immune response, apoptosis, signal transduction, and metabolism-related pathways. These findings provide new insight into the involvement of miRNAs and circRNAs in the host response to T. gondii infection.

An in-silico Approach for Recognition of Long non-coding RNA-Associated Competing Endogenous RNA Axes in Prostate Cancer

PURPOSE

Prostate cancer is among the most central sources of cancer-related mortalities. In order to find novel candidates for therapeutic strategies in this kind of cancer, we developed an in-silico method for identification of competing endogenous RNA network.

METHODS

According to the microarray data analyses between prostate tumor and normal specimens, we attained 1312 differentially expressed (DE)mRNAs, including 778 down-regulated DEmRNAs (such as CXCL13 and BMP5) and 584 up-regulated DEmRNAs (such as OR51E2 and LUZP2), 39 DElncRNAs, including 10 down-regulated DElncRNAs (such as UBXN10-AS1 and FENDRR) and 29 up-regulated DElncRNAs (such as PCA3 and LINC00992) and 10 DEmiRNAs, including 2 down-regulated DEmiRNAs (such as MIR675 and MIR1908) and 8 up-regulated DEmiRNAs (such as MIR6773 and MIR4683).

RESULTS

We constructed the ceRNA network between these transcripts. We also evaluated the related signaling pathways and the significance of these RNAs in prediction of survival of patients with prostate cancer.

CONCLUSION

This study provides novel candidates for construction of specific treatment routes for prostate cancer.

Identification of circRNA-associated ceRNA networks in the longissimus dorsi of yak under different feeding systems

BACKGROUND

Yaks (Bos grunniens), prized for their ability to thrive in high-altitude environments, are indispensable livestock in the plateau region. Modifying their feeding systems holds significant promise for improving their growth and meat quality. Tenderness, a key determinant of yak meat quality and consumer appeal, is demonstrably influenced by dietary regimen. Indoor feeding regimes have been shown to enhance tenderness by lowering shear stress and optimizing pH values. CircRNAs, well-known modulators of circulatory function, also play a crucial role in skeletal muscle development across various animal species. However, their functional significance in yak skeletal muscle remains largely unexplored.

RESULTS

In this study, we identified a total of 5,534 circRNAs within the longissimus dorsi muscle, and we found 51 differentially expressed circRNAs (20 up-regulated and 31 down-regulated) between the two feeding groups. Constructing a comprehensive ceRNA network illuminated intricate regulatory mechanisms, with PGP and circRNA_0617 converging on bta-miR-2285q, mirrored by KLF15/circRNA_0345/bta-miR-20b and CTSF/circRNA_0348/bta-miR-146a. These findings shed light on the potential of circRNAs to influence yak muscle development and meat quality, offering valuable insights for future research.

CONCLUSIONS

This investigation unraveled a complex interaction network between circRNAs、mRNAs and miRNAs in yak skeletal muscle. We further elucidated the target genes regulated by these target genes within the network, offering valuable insights into the potential regulatory mechanisms governing muscle development and meat quality-related traits in yaks.

Identification of functional genes in liver fibrosis based on bioinformatics analysis of a lncRNA-mediated ceRNA network

BACKGROUND

Liver fibrosis is a major global healths problem; nevertheless, its molecular mechanism are not completely clear. This study aimed to build a lncRNA-miRNA-mRNA network, identify potentially related lncRNAs, and explore the pathogenesis of liver fibrosis.

MATERIALS AND METHODS

We used the Gene Expression Omnibus databases and bioinformatics analysis to identify differentially expressed genes (DEGs) between liver fibrosis and normal tissues. The ceRNA network was constructed according to the interactions between DElncRNA, miRNA, and DEmRNA. Then, these DEGs were identified using functional enrichment analysis, and a protein-protein interaction (PPI) network was established. The critical lncRNAs were verified using the quantitative real-time polymerase chain reaction (qRT-PCR).

RESULTS

The ceRNA network was composed of three lncRNAs, five miRNAs, and 93 mRNAs. Gene Ontology functional enrichment analysis revealed significant enhancement in cell components, molecular function, and biological process. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed pathways associated with transcriptional misregulation in cancer, including the Rap1 signaling pathway, proteoglycans in cancer, mineral absorption, HTLV-l infection, and central carbon metabolism in cancer. According to the PPI network and the GSE84044 database, seven hub genes associated with liver fibrosis were identified. In addition, qRT-PCR revealed that lncRNA AC100861 (lncRNA TNFRSF10A-DT) was explicitly decreased in liver fibrosis tissues and activated hepatic stellate cells.

CONCLUSIONS

In summary, this study preliminarily found that lncRNA TNFRSF10A-DT may be a biomarker for the diagnosis and outcome of liver fibrosis. We uncovered a novel lncRNA-mediated ceRNA regulatory mechanism in the pathogenesis of liver fibrosis.

ZNF460-mediated circRPPH1 promotes TNBC progression through ITGA5-induced FAK/PI3K/AKT activation in a ceRNA manner

BACKGROUND

Circular RNAs are highly stable regulatory RNAs that have been increasingly associated with tumorigenesis and progression. However, the role of many circRNAs in triple-negative breast cancer (TNBC) and the related mechanisms have not been elucidated.

METHODS

In this study, we screened circRNAs with significant expression differences in the RNA sequencing datasets of TNBC and normal breast tissues and then detected the expression level of circRPPH1 by qRT‒PCR. The biological role of circRPPH1 in TNBC was then verified by in vivo and in vitro experiments. Mechanistically, we verified the regulatory effects between circRPPH1 and ZNF460 and between circRPPH1 and miR-326 by chromatin immunoprecipitation (ChIP), fluorescence in situ hybridization assay, dual luciferase reporter gene assay and RNA pull-down assay. In addition, to determine the expression of associated proteins, we performed immunohistochemistry, immunofluorescence, and western blotting.

RESULTS

The upregulation of circRPPH1 in TNBC was positively linked with a poor prognosis. Additionally, both in vivo and in vitro, circRPPH1 promoted the biologically malignant behavior of TNBC cells. Additionally, circRPPH1 may function as a molecular sponge for miR-326 to control integrin subunit alpha 5 (ITGA5) expression and activate the focal adhesion kinase (FAK)/PI3K/AKT pathway.

CONCLUSION

Our research showed that ZNF460 could promote circRPPH1 expression and that the circRPPH1/miR-326/ITGA5 axis could activate the FAK/PI3K/AKT pathway to promote the progression of TNBC. Therefore, circRPPH1 can be used as a therapeutic or diagnostic target for TNBC.

Cuproptosis-related ceRNA axis triggers cell proliferation and cell cycle through CBX2 in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) has high morbidity and mortality. Despite substantial advances in treatment, the prognosis of patients with LUAD remains unfavorable. The ceRNA axis has been reported to play an important role in the pathogenesis of LUAD. In addition, cuproptosis is considered an important factor in tumorigenesis. The expression of CBX2 has been associated with the development of multiple tumors, including LUAD. However, the precise molecular mechanisms through which the cuproptosis-related ceRNA network regulates CBX2 remain unclear.

METHODS

The DEGs between tumor and normal samples of LUAD were identified in TCGA database. The "ConsensusClusterPlus" R package was used to perform consensus clustering based on the mRNA expression matrix and cuproptosis-related gene expression profile. Then, LASSO-COX regression analysis was performed to identify potential prognostic biomarkers associated with cuproptosis, and the ceRNA network was constructed. Finally, the mechanisms of ceRNA in LUAD was studied by cell experiments.

RESULTS

In this study, the AC144450.1/miR-424-5p axis was found to promote the progression of LUAD by acting on CBX2. The expression of AC144450.1 and miR-424-5p can be altered to regulate CBX2 and is correlated with cell proliferation and cell cycle of LUAD. Mechanistically, AC144450.1 affects the expression of CBX2 by acting as the ceRNA of miR-424-5p. In addition, a cuproptosis-related model were constructed in this study to predict the prognosis of LUAD.

CONCLUSIONS

This study is the first to demonstrate that the AC144450.1/miR-424-5p/CBX2 axis is involved in LUAD progression and may serve as a novel target for its diagnosis and treatment.

Integrated analysis reveals ceRNA network of cardiac remodeling by SGLT2 inhibitor in middle-aged hypertensive rats

Sodium-glucose cotransporter 2 inhibitors (SGLT2i) represent an innovative class of antidiabetic agents that have demonstrated promise in mitigating cardiac remodeling. However, the transcriptional regulatory mechanisms underpinning their impact on blood pressure and the reversal of hypertension-induced cardiac remodeling remain largely unexplored. Given this context, our study concentrated on comparing the cardiac expression profiles of lncRNAs and mRNAs between Wistar-Kyoto (WKY) rats and spontaneously hypertensive rats (SHR). To validate our results, we performed blood pressure measurements, tissue staining, and qRT-PCR. The treatment led to a significant reduction in systolic blood pressure and improved cardiac remodeling by reducing myocardial fibrosis and regulating the inflammatory response. Our examination disclosed that ventricular tissue mRNA, regulated by hypertension, was primarily concentrated in the complement and coagulation cascades and cytokine-cytokine receptor interactions. Compared with SHR, the SGLT2i treatment group was associated with myocardial contraction. Investigation into the lncRNA-mRNA regulatory network and competing endogenous RNA (ceRNA) network suggested that the potential roles of these differentially expressed (DE) lncRNAs and mRNAs were tied to processes such as collagen fibril organization, inflammatory response, and extracellular matrix (ECM) modifications. We found that the expression of Col3a1, C1qa, and lncRNA NONRATT007139.2 were altered in the SHR group and that SGLT2i treatment reversed these changes. Our results suggest that dapagliflozin effectively reverses hypertension-induced myocardial remodeling through a lncRNA-mRNA transcriptional regulatory network, with immune cell-mediated ECM deposition as a potential regulatory target. This underlines the potentiality of SGLT2i and genes related to immunity as promising targets for the treatment of hypertension.

RNA sequencing reveals differential long noncoding RNA expression profiles in bacterial and viral meningitis in children

BACKGROUND

We aimed to investigate the involvement of long non-coding RNA (lncRNA) in bacterial and viral meningitis in children.

METHODS

The peripheral blood of five bacterial meningitis patients, five viral meningitis samples, and five healthy individuals were collected for RNA sequencing. Then, the differentially expressed lncRNA and mRNA were detected in bacterial meningitis vs. controls, viral meningitis vs. healthy samples, and bacterial vs. viral meningitis patients. Besides, co-expression and the competing endogenous RNA (ceRNA) networks were constructed. Receiver operating characteristic curve (ROC) analysis was performed.

RESULTS

Compared with the control group, 2 lncRNAs and 32 mRNAs were identified in bacterial meningitis patients, and 115 lncRNAs and 54 mRNAs were detected in viral meningitis. Compared with bacterial meningitis, 165 lncRNAs and 765 mRNAs were identified in viral meningitis. 2 lncRNAs and 31 mRNAs were specific to bacterial meningitis, and 115 lncRNAs and 53 mRNAs were specific to viral meningitis. The function enrichment results indicated that these mRNAs were involved in innate immune response, inflammatory response, and immune system process. A total of 8 and 1401 co-expression relationships were respectively found in bacterial and viral meningitis groups. The ceRNA networks contained 1 lncRNA-mRNA pair and 4 miRNA-mRNA pairs in viral meningitis group. GPR68 and KIF5C, identified in bacterial meningitis co-expression analysis, had an area under the curve (AUC) of 1.00, while the AUC of OR52K2 and CCR5 is 0.883 and 0.698, respectively.

CONCLUSIONS

Our research is the first to profile the lncRNAs in bacterial and viral meningitis in children and may provide new insight into understanding meningitis regulatory mechanisms.

Identification and Characterization of a ceRNA Regulatory Network Involving LINC00482 and PRRC2B in Peripheral Blood Mononuclear Cells: Implications for COPD Pathogenesis and Diagnosis

Purpose

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death worldwide, characterized by intense lung infiltrations of immune cells (macrophages and monocytes). While existing studies have highlighted the crucial role of the competitive endogenous RNA (ceRNA) regulatory network in COPD development, the complexity and characteristics of the ceRNA network in monocytes remain unexplored.

Methods

We downloaded messenger RNA (mRNA), microRNA (miRNA), and long noncoding RNA (lncRNA) microarray data from GSE146560, GSE102915, and GSE71220 in the Gene Expression Omnibus (GEO) database. This data was used to identify differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs), and lncRNAs (DElncRNAs). Predicted miRNAs that bind to DElncRNAs were intersected with DEmiRNAs, forming a set of intersecting miRNAs. This set was then used to predict potential binding mRNAs, intersected with DEmRNAs, and underwent functional enrichment analysis using R software and the STRING database. The resulting triple regulatory network and hub genes were constructed using Cytoscape. Comparative Toxicomics Database (CTD) was utilized for disease correlation predictions, and ROC curve analysis assessed diagnostic accuracy.

Results

Our study identified 5 lncRNAs, 4 miRNAs, and 149 mRNAs as differentially expressed. A lncRNA-miRNA-mRNA regulatory network was constructed, and hub genes were selected through hub analysis. Enrichment analysis highlighted terms related to cell movement and gene expression regulation. We established a LINC00482-has-miR-6088-PRRC2B ceRNA network with diagnostic relevance for COPD. ROC analysis demonstrated the diagnostic value of these genes. Moreover, a positive correlation between LINC00482 and PRRC2B expression was observed in COPD PBMCs. The CTD database indicated their involvement in inflammatory responses.

Conclusion

In summary, our study not only identified pivotal hub genes in peripheral blood mononuclear cells (PBMCs) of COPD but also constructed a ceRNA regulatory network. This contributes to understanding the pathophysiological processes of COPD through bioinformatics analysis, expanding our knowledge of COPD, and providing a foundation for potential diagnostic and therapeutic targets for COPD.

Construction of a New Immune-Related Competing Endogenous RNA Network with Prognostic Value in Lung Adenocarcinoma

Tumor microenvironment has significant influence on the gene expression of tumor tissues and on the clinical outcomes in lung adenocarcinoma. Infiltrating immune and stromal cells not only perturb the tumor signal in molecular studies, but also play crucial roles in cancer biology. The competing endogenous RNAs (ceRNAs) are useful to explain the post-transcriptional layer regulated by gene translation and play an important role in the occurrence and progression of lung adenocarcinoma. Therefore, identifying novel molecular markers by constructing ceRNA associated with immune infiltration is of great significance to guide the treatment of lung adenocarcinoma in the future. According to the immune and stromal scores of lung adenocarcinoma samples in The Cancer Genome Atlas (TCGA) database calculated by the ESTIMATE algorithm, we identified differentially expressed lncRNAs, miRNAs and mRNAs associated with immune infiltration, including 60 dysregulated lncRNAs, 38 dysregulated mRNAs, and 29 dysregulated miRNAs. Based on the PPI network and Cox regression analysis, 5 mRNAs including CNR2, P2RY12, ZNF831, RSPO1, and F2 were identified to be related to immune infiltration and prognosis in lung adenocarcinoma, and their differential expression, prognosis and correlation with immune cells were verified. Next, through target binding prediction, pearson correlation analysis and expression analysis, a novel immune-related ceRNA network containing 6 lncRNAs, 4 miRNAs, and 3 mRNAs was finally constructed. The present study constructed a novel immune-associated lncRNA-miRNA-mRNA ceRNA network, which deepens our understanding on the molecular network mechanism of lung adenocarcinoma and provides potential prognostic markers and novel therapeutic targets for the patients with lung adenocarcinoma.

The identification of regulatory ceRNA network involved in Drosophila Toll immune responses

Non-coding RNAs play important roles in the innate immunity of Drosophila, with various lncRNAs and miRNAs identified to maintain Drosophila innate immune homeostasis by regulating protein functions. However, it remains unclear whether interactions between lncRNAs and miRNAs give rise to a ceRNA network. In our previous study, we observed the highest differential expression levels of lncRNA-CR11538, lncRNA-CR33942, and lncRNA-CR46018 in wild-type flies after Gram-positive bacterial infection, prompting us to investigate their role in the regulation of Drosophila Toll immune response through RNA-seq analysis. Herein, our comprehensive bioinformatics analysis revealed that lncRNA-CR11538, lncRNA-CR33942, and lncRNA-CR46018 are involved in defense mechanisms and stimulus response. Moreover, lncRNA-CR11538 and lncRNA-CR46018 can also participate in the metabolic recovery processes following Gram-positive bacterial infection. Subsequently, we employed GSEA screening and RT-qPCR to identify seven miRNAs (miR-957, miR-1015, miR-982, miR-993, miR-1007, miR-193, and miR-978) that may be regulated by these three lncRNAs. Furthermore, we predicted the potential target genes in the Toll signaling pathway for these miRNAs and their interaction with the three lncRNAs using TargetScan and miRanda software and preliminary verification. As a result, we established a potential ceRNA regulatory network for Toll immune responses in Drosophila, comprising three lncRNAs and seven miRNAs. This study provides evidence of a ceRNA regulatory network in Drosophila Toll immune responses and offers novel insights into understanding the regulatory networks involved in the innate immunity of other animals.

LncRNA pol-lnc78 as a ceRNA regulates antibacterial responses via suppression of pol-miR-n199-3p-mediated SARM down-regulation in Paralichthys olivaceus

Long non-coding RNAs (lncRNAs) function as key modulators in mammalian immunity, particularly due to their involvement in lncRNA-mediated competitive endogenous RNA (ceRNA) crosstalk. Despite their recognized significance in mammals, research on lncRNAs in lower vertebrates remains limited. In the present study, we characterized the first immune-related lncRNA (pol-lnc78) in the teleost Japanese flounder ( Paralichthys olivaceus). Results indicated that pol-lnc78 acted as a ceRNA for pol-miR-n199-3p to target the sterile alpha and armadillo motif-containing protein (SARM), the fifth discovered member of the Toll/interleukin 1 (IL-1) receptor (TIR) adaptor family. This ceRNA network regulated the antibacterial responses of flounder via the Toll-like receptor (TLR) signaling pathway. Specifically, SARM acted as a negative regulator and exacerbated bacterial infection by inhibiting the expression of inflammatory cytokines IL-1β and tumor necrosis factor-α (TNF-α). Pol-miR-n199-3p reduced SARM expression by specifically interacting with the 3' untranslated region (UTR), thereby promoting SARM-dependent inflammatory cytokine expression and protecting the host against bacterial dissemination. Furthermore, pol-lnc78 sponged pol-miR-n199-3p to ameliorate the inhibition of SARM expression. During infection, the negative regulators pol-lnc78 and SARM were significantly down-regulated, while pol-miR-n199-3p was significantly up-regulated, thus favoring host antibacterial defense. These findings provide novel insights into the mechanisms underlying fish immunity and open new horizons to better understand ceRNA crosstalk in lower vertebrates.

A circRNA ceRNA network involved in cognitive dysfunction after chronic cerebral hypoperfusion

Chronic Cerebral Hypoperfusion (CCH) is associated with cognitive dysfunction, the underlying mechanisms of which remain elusive, hindering the development of effective therapeutic approaches. In this study, we employed an established CCH animal model to delve into neuropathological alterations like oxidative stress, inflammation, neurotransmitter synthesis deficits, and other morphological alterations. Our findings revealed that while the number of neurons remained unchanged, there was a significant reduction in neuronal fibers post-CCH, as evidenced by microtubule-associated protein 2 (MAP2) staining. Moreover, myelin basic protein (MBP) staining showed exacerbated demyelination of neuronal fibers. Furthermore, we observed increased neuroinflammation, proliferation, and activation of astrocytes and microglia, as well as synaptic loss and microglial-mediated synapse engulfment post-CCH. Utilizing RNA sequencing, differential expression analysis displayed alterations in both mRNAs and circRNAs. Following gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, both showed significant enrichment in immunological and inflammation-related terms and pathways. Importantly, the differentially expressed circular RNAs (DE circRNAs) exhibited a notable coexpression pattern with DE mRNAs. The ternary circRNA-miRNA-mRNA competing endogenous RNAs (ceRNA) network was constructed, and subsequent analysis reiterated the significance of neuroimmunological and neuroinflammatory dysfunction in CCH-induced neuropathological changes and cognitive dysfunction. This study underscores the potential role of circRNAs in these processes, suggesting them as promising therapeutic targets to mitigate the detrimental effects of CCH.

SPC25 as a novel therapeutic and prognostic biomarker and its association with glycolysis, ferroptosis and ceRNA in lung adenocarcinoma

OBJECTIVE

Spindle pole body component 25 (SPC25) is an important cyclin involved in chromosome segregation and spindle dynamics regulation during mitosis. However, the role of SPC25 in lung adenocarcinoma (LAUD) is unclear.

MATERIALS AND METHODS

The differential expression of SPC25 in tumor samples and normal samples was analyzed using TIMER, TCGA, GEO databases, and the correlation between its expression and clinicopathological features and prognosis in LUAD patients. Biological pathways that may be enriched by SPC25 were analyzed using GSEA. In vitro cell experiments were used to evaluate the effect of knocking down SPC25 expression on LUAD cells. Correlation analysis and differential analysis were used to assess the association of SPC25 expression with genes related to cell cycle, glycolysis, and ferroptosis. A ceRNA network involving SPC25 was constructed using multiple database analyses.

RESULTS

SPC25 was highly expressed in LUAD, and its expression level could guide staging and predict prognosis. GSEA found that high expression of SPC25 involved multiple cell cycles and glycolytic pathways. Knocking down SPC25 expression significantly affected the proliferation, migration and apoptosis of LUAD cells. Abnormal SPC25 expression levels can affect cell cycle progression, glycolytic ability and ferroptosis regulation. A ceRNA network containing SPC25, SNHG15/hsa-miR-451a/SPC25, was successfully predicted and constructed.

CONCLUSIONS

Our findings reveal the association of up-regulation of SPC25 in LUAD and its expression with clinical features, prognosis prediction, proliferation migration, cell cycle, glycolysis, ferroptosis, and ceRNA networks. Our results indicate that SPC25 can be used as a biomarker in LUAD therapy and a target for therapeutic intervention.

Tumor-Derived Exosomal LINC01480 Upregulates VCAM1 Expression by Acting as a Competitive Endogenous RNA of miR-204-5p to Promote Gastric Cancer Progression

Exosomes are a type of cell-derived vesicles that range in size from 30 to 100 nm. They are widely present in various organisms and participate in diverse biological processes, playing crucial roles in tumorigenesis and progression. This study aimed to investigate whether LINC01480 in tumor-derived exosomes is involved in the molecular mechanism of gastric cancer by competitively upregulating the VCAM1 expression through binding miR-204-5p. The study analyzed transcriptome data related to gastric cancer from the cancer genome atlas database and constructed a risk-scoring model for epithelial-mesenchymal transition (EMT)-related lncRNAs to identify eight EMT-related lncRNAs associated with prognosis. EMT-related mRNAs positively correlated with LINC01480 were screened in the ExoRBase database. In vitro cell experiments showed that exosomal LINC01480 can promote the proliferation, migration, invasion, and EMT of gastric cancer cells by upregulating VCAM1 expression through competitive binding with miR-204-5p. In vivo experiments on nude mice showed that exosomal LINC01480 promotes the development of gastric cancer. These results suggest that exosomal LINC01480 could be a potential therapeutic target for gastric cancer.

Transcriptomic analysis identifies diagnostic genes in polycystic ovary syndrome and periodontitis

PURPOSE

To investigate underlying co-mechanisms of PCOS and periodontitis through transcriptomic approach.

METHODS

PCOS and periodontitis gene expression data were downloaded from the GEO database to identify differentially expressed genes. GO and KEGG pathway enrichment analysis and random forest algorithm were used to screen hub genes. GSEA analyzed the functions of hub genes. Correlations between hub genes and immune infiltration in two diseases were examined, constructing a TF-ceRNA regulatory network. Clinical samples were gathered from PCOS and periodontitis patients and RT-qPCR was performed to verify the connection.

RESULTS

There were 1661 DEGs in PCOS and 701 DEGs in periodontitis. 66 intersected genes were involved and were enriched in immune and inflammation-related biological pathways. 40 common genes were selected from the PPI network. RF algorithm demonstrated that ACSL5, NLRP12, CCRL2, and CEACAM3 were hub genes, and GSEA results revealed their close relationship with antigen processing and presentation, and chemokine signaling pathway. RT-qPCR results confirmed the upregulated gene expression in both PCOS and periodontitis.

CONCLUSION

The 4 hub genes ACSL5, NLRP12, CCRL2, and CEACAM3 may be diagnostic genes for PCOS and periodontitis. The created ceRNA network could provide a molecular basis for future studies on the association between PCOS and periodontitis.