Building a Human Ovarian Antioxidant ceRNA Network "OvAnOx": A Bioinformatic Perspective for Research on Redox-Related Ovarian Functions and Dysfunctions.

The ovary is a major determinant of female reproductive health. Ovarian functions are mainly related to the primordial follicle pool, which is gradually lost with aging. Ovarian aging and reproductive dysfunctions share oxidative stress as a common underlying mechanism. ROS signaling is essential for normal ovarian processes, yet it can contribute to various ovarian disorders when disrupted. Therefore, balance in the redox system is crucial for proper ovarian functions. In the present study, by focusing on mRNAs and ncRNAs described in the ovary and taking into account only validated ncRNA interactions, we built an ovarian antioxidant ceRNA network, named OvAnOx ceRNA, composed of 5 mRNAs (SOD1, SOD2, CAT, PRDX3, GR), 10 miRNAs and 5 lncRNAs (XIST, FGD5-AS1, MALAT1, NEAT1, SNHG1). Our bioinformatic analysis indicated that the components of OvAnOx ceRNA not only contribute to antioxidant defense but are also involved in other ovarian functions. Indeed, antioxidant enzymes encoded by mRNAs of OvAnOx ceRNA operate within a regulatory network that impacts ovarian reserve, follicular dynamics, and oocyte maturation in normal and pathological conditions. The OvAnOx ceRNA network represents a promising tool to unravel the complex dialog between redox potential and ovarian signaling pathways involved in reproductive health, aging, and diseases.

Mitigating off-target effects of small RNAs: conventional approaches, network theory and artificial intelligence.

Three types of highly promising small RNA therapeutics, namely, small interfering RNAs (siRNAs), microRNAs (miRNAs) and the RNA subtype of antisense oligonucleotides (ASOs), offer advantages over small-molecule drugs. These small RNAs can target any gene product, opening up new avenues of effective and safe therapeutic approaches for a wide range of diseases. In preclinical research, synthetic small RNAs play an essential role in the investigation of physiological and pathological pathways as silencers of specific genes, facilitating discovery and validation of drug targets in different conditions. Off-target effects of small RNAs, however, could make it difficult to interpret experimental results in the preclinical phase and may contribute to adverse events of small RNA therapeutics. Out of the two major types of off-target effects we focused on the hybridization-dependent, especially on the miRNA-like off-target effects. Our main aim was to discuss several approaches, including sequence design, chemical modifications and target prediction, to reduce hybridization-dependent off-target effects that should be considered even at the early development phase of small RNA therapy. Because there is no standard way of predicting hybridization-dependent off-target effects, this review provides an overview of all major state-of-the-art computational methods and proposes new approaches, such as the possible inclusion of network theory and artificial intelligence (AI) in the prediction workflows. Case studies and a concise survey of experimental methods for validating in silico predictions are also presented. These methods could contribute to interpret experimental results, to minimize off-target effects and hopefully to avoid off-target-related adverse events of small RNA therapeutics.

Novel insights into sevoflurane-induced developmental neurotoxicity mechanisms.

Aim: This study explores Sevoflurane (Sevo)-induced neurotoxicity mechanisms in neonates through transcriptome sequencing and models.Methods: Seven-day-old mice were exposed to 3% Sevo, and hippocampal tissue was collected for analysis of differentially expressed lncRNAs and mRNAs compared with normal mice. MiR-152-3p was selected, and the interaction between H19, USP30, and miR-152-3p was explored in BV2 microglial cells and mouse hippocampal neurons.Results: Sevo disrupts mitochondrial autophagy via USP30 upregulation, exacerbating neurotoxicity and activating NLRP1 inflammasome-mediated inflammation.Conclusion: Sevo neurotoxicity is mediated through the H19/miR-152-3p/USP30 axis, implicating microglial regulation of neuronal pyroptosis.

[Box: see text].

Integration Analysis of miRNA Circulating Expression Following Cerebellar Transcranial Direct Current Stimulation in Patients with Ischemic Stroke.

The aim of this study was to explore the molecular mechanisms underlying cerebellar transcranial direct current stimulation (ctDCS) as a rehabilitation intervention for patients with ischemic stroke, focusing on the role of microRNAs (miRNAs). Whole-transcriptome sequencing was employed to obtain circulating expression profiles of miRNAs, long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and mRNAs in patients with ischemic stroke before and after 3-week ctDCS. miRanda software was used to predict the target genes of miRNAs, while Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted to identify biological functions and signaling pathways. Subsequently, competing endogenous RNA (ceRNA) regulatory networks comprising circRNA-miRNA-mRNA and lncRNA-miRNA-mRNA interactions were constructed. Key miRNAs in blood samples were validated through quantitative RT-PCR. In total, 43 miRNAs, 807 lncRNAs, 1,111 circRNAs, and 201 mRNAs were differentially expressed after ctDCS compared with before ctDCS. Bioinformatics analyses revealed significant enrichment of target genes regulated by differentially expressed miRNAs across multiple biological pathways. CeRNA regulatory networks implied that several miRNAs were closely related to the ctDCS. Among them, hsa-miR-181a-5p, hsa-miR-224-5p, and hsa-miR-340-3p showed significantly downregulated expression levels as confirmed by qRT-PCR. This study conducted the first-ever assessment of miRNA expression patterns in patients with ischemic stroke undergoing ctDCS. The findings revealed that ctDCS induces alterations in miRNA levels, suggesting their potential utility as therapeutic markers.

In-Silico and In-Vitro Investigation of Key Long Non-coding RNAs Involved in 5-Fluorouracil Resistance in Colorectal Cancer Cells: Analyses Highlighting NEAT1 and MALAT1 as Contributors.

Background Acquired resistance to 5-fluorouracil (5-FU) frequently results in chemotherapy failure and disease recurrence in advanced colorectal cancer (CRC) patients. Research has demonstrated that dysregulation of long non-coding RNAs (lncRNAs) mediates the development of chemotherapy resistance in cancerous cells. The present study aims to identify key lncRNAs associated with 5-FU resistance in CRC using bioinformatic and experimental validation approaches. Methods The Gene Expression Omnibus (GEO) dataset GSE119481, which contains miRNA expression profiles of the parental CRC HCT116 cell line (HCT116/P) and its in-vitro established 5-FU-resistant sub-cell line (HCT116/FUR), was downloaded. Firstly, differentially expressed microRNAs (DEmiRNAs) between the parental and 5-FU resistance cells were identified. LncRNAs and mRNAs were then predicted using online databases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to uncover relevant biological mechanisms and pathways. Networks integrating lncRNAs, miRNAs, and mRNAs interactions were constructed, and topological analyses were used to identify key lncRNAs associated with 5-FU resistance. An in-vitro model of the HCT116/FUR sub-cell line was developed by exposing the HCT116/P cell line to increasing concentrations of 5-FU. Finally, real-time quantitative PCR (RT-qPCR) was performed on total RNA extracted from the HCT116/P cell line and the HCT116/FUR sub-cell line to validate the in-silico predictions of key lncRNAs. Results A total of 32 DEmiRNAs were identified. Enrichment analysis demonstrated that these DEmiRNAs were mainly enriched in several cancer hallmark pathways that regulate cell growth, cell cycle, cell survival, inflammation, immune response, and apoptosis. The predictive analysis identified 237 unique lncRNAs and 123 mRNAs interacting with these DEmiRNAs. The pathway analysis indicated that most of these predicted genes were enriched in the cellular response to starvation, protein polyubiquitination, chromatin remodeling, and negative regulation of gene expression. Topological analyses of the lncRNA-miRNA-mRNA network highlighted the nuclear enriched abundant transcript 1 (NEAT1), metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), and Opa interacting protein 5 antisense RNA 1 (OIP5-AS1) as central lncRNAs. Experimental analysis by RT-qPCR confirmed that the expression levels of NEAT1 and MALAT1 were significantly increased in HCT116/FUR cells compared to HCT116/P cells. However, no significant difference was observed in the OIP5-AS1 expression level between the two cells. Conclusion Our findings specifically highlight MALAT1 and NEAT1 as significant contributors to 5-FU resistance in CRC. These lncRNAs are promising biomarkers for diagnosing and predicting outcomes in CRC.

Comprehensive analysis of cuproptosis-related ceRNA network and immune infiltration in diabetic kidney disease.

Background: Diabetic kidney disease (DKD) is the primary contributor to renal failure and poses a severe threat to human health. Accumulating studies demonstrated that competing endogenous RNA (ceRNA) network is involved in cuproptosis and DKD progression. However, the role of cuproptosis-associated ceRNA network and immune infiltration in DKD remains largely unclear. This study aimed to investigate the cuproptosis-related ceRNA regulation network and immune infiltration in DKD. Methods: The rat model of DKD was induced by combining the nephrectomy of the left kidney, high-fat diet, and streptozotocin. Differentially expressed genes (DEGs), miRNAs (DEMs), and lncRNAs (DELs) between normal and DKD rats were obtained. DEGs were intersected with cuproptosis-related genes (CRGs) to obtain DE-CRGs. LncRNAs and miRNAs were predicted based on the DE-CRGs, and they were intersected with DEMs and DELs, respectively. Subsequently, a cuproptosis-associated lncRNA-miRNA-mRNA network was established in DKD. In addition, the relative proportion of 22 infiltrating immune cell types in each sample was calculated, and the relationship between hub DE-CRGs and immune cells was explored. Results: In total, there were 429 DEGs, 22 DEMs, and 48 DELs between CON and MOD groups. Then, 73 DE-CRGs were obtained, which were significantly enriched in 22 pathways, such as MAPK signaling pathway, IL-17 signaling pathway, and TNF signaling pathway. In addition, a core cuproptosis-related ceRNA network that included one lncRNA (USR0000B2476D), one miRNA (miR-34a-3p), and eight mRNAs (Mmp9, Pik3c3, Prom1, Snta1, Slc51b, Ntrk3, Snca, Egf) was established. In addition, 18 hub DE-CRGs were obtained. CIBERSORT algorithms showed that resting dendritic cells and resting NK cells were more infiltrated whereas regulatory T cells were less infiltrated in DKD rats than in normal rats. Spearman's correlation analysis revealed that hub DE-CRGs showed significant positive or negative correlations with naive B cells, regulatory T cells, resting NK cells, M0 macrophages, resting dendritic cells, and resting mast cells. Conclusion: A ceRNA network was comprehensively constructed, and 18 hub DE-CRGs were obtained, which will provide novel insights into the pathologic mechanism elucidation and targeted therapy development of DKD.

Identification of key immune-related genes and potential therapeutic targets in immune checkpoint inhibitor-associated myocarditis.

BACKGROUND: Immune checkpoint inhibitors (ICIs) are widely used in cancer treatment; however, the emergence of ICI-associated myocarditis (ICI-MC) presents a severe and potentially fatal complication with poorly understood pathophysiological mechanisms. This study aimed to identify crucial immune-related genes in ICI-MC and uncover potential therapeutic targets using bioinformatics. METHODS: Using the GSE180045 dataset, which includes three groups-Group A: ICI patients without immune adverse events, Group B: ICI patients with non-myocarditis immune adverse events, and Group C: ICI patients with myocarditis-we analyzed differentially expressed genes (DEGs) between ICI-MC samples (Group C) and non-myocarditis controls (Groups A and B). These DEGs were then cross-referenced with 1796 immune-related genes from the immPort database to identify immune-related DEGs. We conducted functional enrichment analyses (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, gene set enrichment analysis), constructed a protein-protein interaction network, and identified hub genes. Validation using the GSE4172 dataset led to the identification of optimal feature genes from the overlap between hub genes and DEGs. Predictions of target MicroRNAs (miRNAs) were made, and a competing endogenous RNA (ceRNA) network was constructed. Target drugs for hub genes were predicted using the Connectivity Map database. RESULTS: We identified 58 DEGs between ICI-MC and controls, which led to the identification of 32 immune-related DEGs after intersection with 1796 immune-related genes. Functional analyses revealed enrichment in cell lysis, CD8+ T-cell receptor, natural killer cell-mediated cytotoxicity, and RAGE signaling. Notably upregulated hub genes included IL7R, PRF1, GNLY, CD3G, NKG7, GZMH, GZMB, KLRB1, KLRK1, and CD247. In the validation dataset, 407 DEGs were uncovered, resulting in the identification of 3 optimal feature genes (KLRB1, NKG7, GZMH). The predicted target miRNAs, lincRNAs, and circRNAs constituted a comprehensive ceRNA network. Among the top 10 drugs with elevated connectivity scores was acetohydroxamic acid, indicating a need for caution in ICI treatment. CONCLUSION: KG7, GZMH, and KLRB1 were identified as pivotal immune-related genes in ICI-MC. Biological enrichments included pathways involved in cell lysis, the CD8+ T-cell receptor pathway, natural killer cell-mediated cytotoxicity, RAGE signaling, and proinflammatory responses. The ceRNA network illuminated the role of critical molecules and underscored the importance of avoiding drugs such as acetohydroxamic acid in ICI treatment. Key message What is already known on this topic  Myocarditis is recognized as a serious ICI-associated toxicity, seemingly infrequent yet often fulminant and lethal. The underlying mechanisms of ICI-associated myocarditis remain not fully understood. Although the significance of T cells and cytotoxic T lymphocyte-associated protein 4 (CTLA-4) is evident, the inciting antigens, the reasons for their recognition, and the mechanisms causing cardiac cell injury are not well characterized. An improved understanding of ICI-associated myocarditis will provide insights into the equilibrium between the immune and cardiovascular systems. What this study adds  Our study further validates the significance of T cells and CTLA-4 in ICI-associated myocarditis. More importantly, we identified three genes-NKG7, GZMH, and KLRB1-essential for the development of ICI-MC and proposed ceRNA networks involving these three key genes. How this study might affect research, practice or policy  The newly discovered key genes and their intricate molecular interactions offer a comprehensive perspective on the mechanisms underlying ICI-MC. Furthermore, our findings advise caution regarding the use of drugs like acetohydroxamic acid during ICI treatment. As our understanding of these regulatory networks deepens, our study provides valuable insights that could inform future therapeutic strategies for ICI-MC.

Bioinformatics analysis of neutrophil-associated hub genes and ceRNA network construction in septic cardiomyopathy.

Septic cardiomyopathy (SCM) is a critical sepsis complication characterized by reversible cardiac depression during early septic shock. Neutrophils, integral to innate immunity, can mediate organ damage when abnormal, but their specific role in sepsis-induced myocardial damage remains elusive. Our study focuses on elucidating the role of Neutrophil-Related Genes (NRGs) in SCM, finding early diagnosis and treatment biomarkers. We identified shared differentially expressed genes (DEGs) from datasets GSE79962 and GSE44363 and pinpointed hub DEGs using the cytoHubba plugin in Cytoscape software. The Neutrophil-Related Hub Gene (NRHG) MRC1 was identified via intersecting hub DEGs with NRGs from WGCNA. We validated MRC1's abnormal expression in SCM using our data and external datasets. Furthermore, a neutrophil-related ceRNA network (AC145207.5/ miR-23a-3p/MRC1) was constructed and validated. Our findings reveal MRC1 as a potential NRHG in SCM pathogenesis, offering insights into neutrophil-mediated mechanisms in SCM and providing a novel molecular target for early diagnosis and intervention in SCM.

Source leaves are regulated by sink strengths through non-coding RNAs and alternative polyadenylation in cucumber (Cucumis sativus L.).

BACKGROUND: The yield of major crops is generally limited by sink capacity and source strength. Cucumber is a typical raffinose family oligosaccharides (RFOs)-transporting crop. Non-coding RNAs and alternative polyadenylation (APA) play important roles in the regulation of growth process in plants. However, their roles on the sink‒source regulation have not been demonstrated in RFOs-translocating species. RESULTS: Here, whole-transcriptome sequencing was applied to compare the leaves of cucumber under different sink strength, that is, no fruit-carrying leaves (NFNLs) and fruit-carrying leaves (FNLs) at 12th node from the bottom. The results show that 1101 differentially expressed (DE) mRNAs, 79 DE long non-coding RNAs (lncRNAs) and 23 DE miRNAs were identified, which were enriched in photosynthesis, energy production and conversion, plant hormone signal transduction, starch and carbohydrate metabolism and protein synthesis pathways. Potential co-expression networks like, DE lncRNAs-DE mRNAs/ DE miRNAs-DE mRNAs, and competing endogenous RNA (ceRNA) regulation models (DE lncRNAs-DE miRNAs-DE mRNAs) associated with sink‒source allocation, were constructed. Furthermore, 37 and 48 DE genes, which enriched in MAPK signaling and plant hormone signal transduction pathway, exist differentially APA, and SPS (CsaV3_2G033300), GBSS1 (CsaV3_5G001560), ERS1 (CsaV3_7G029600), PNO1 (CsaV3_3G003950) and Myb (CsaV3_3G022290) may be regulated by both ncRNAs and APA between FNLs and NFNLs, speculating that ncRNAs and APA are involved in the regulation of gene expression of cucumber sink‒source carbon partitioning. CONCLUSIONS: These results reveal a comprehensive network among mRNAs, ncRNAs, and APA in cucumber sink-source relationships. Our findings also provide valuable information for further research on the molecular mechanism of ncRNA and APA to enhance cucumber yield.

Whole-transcriptome sequencing revealed the ceRNA regulatory network during the proliferation and differentiation of goose myoblast.

The Shitou goose, the largest meat-type goose breed, is an ideal model for offering insights into enhancing meat production efficiency through understanding its genetic regulation of muscle development. Here, through whole-transcriptomic analysis of embryonic leg muscles, we identified 847 differentially expressed genes (DEG), 244 differentially expressed lncRNAs (DEL), 37 differentially expressed circRNAs (DEC), and 84 differentially expressed miRNAs (DEM). Gene ontology (GO) analysis highlighted the significant enrichment of differentially expressed RNAs in muscle structure development, actin filament-based processes, and the actin cytoskeleton pathway. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified pathways associated with the FoxO signaling pathway, AMPK signaling pathway, Wnt signaling pathway and calcium signaling pathway. Furthermore, we utilized Miranda, TargetScan, and miRDB to identify regulatory networks that involve interactions between lncRNA-mRNA, circRNA-mRNA, miRNA-mRNA, lncRNA-miRNA-mRNA, and circRNA-miRNA-mRNA, which regulated the growth and development of skeletal muscle. Notably, differentially expressed genes within the ceRNA network were most significantly enriched in the regulation of actin cytoskeletal organization. Additionally, a lncRNA/circRNA-miRNA-mRNA ceRNA network related to muscle growth and development was constructed based on protein-protein interaction (PPI) analysis and hub genes selection using Cytoscape. This further elucidated the regulatory roles of noncoding RNAs (ncRNA) in the formation of muscle fibers in Shitou goose. In summary, this study provides a valuable transcriptional regulatory network for goose muscle development laying the groundwork for further exploration of the molecular regulatory mechanisms underlying the excellent meat production performance of Shitou goose.

Preliminary construction of non-coding RNAs and ceRNA regulatory networks mediated by exosomes in porcine follicular fluid.

BACKGROUND: Follicles are fundamental units of the ovary, regulated intricately during development. Exosomes and ovarian granulosa cells (OGCs) play pivotal roles in follicular development, yet the regulatory mechanisms governing exosomes remain elusive. RESULTS: High-throughput sequencing was employed to evaluate the complete transcript expression profiles of six samples (three porcine ovarian granulosa cells-exosome co-culture samples (GCE) and three porcine ovarian granulosa cells (POGCs) samples). Differential expression analysis revealed 924 lncRNAs, 35 circRNAs, 49 miRNAs, and 9823 mRNAs in the GCE group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated enrichment of differentially expressed transcripts in pathways related to cell proliferation and apoptosis. Furthermore, a ceRNA regulatory network comprising 43 lncRNAs, 6 circRNAs, 11 miRNAs, and 126 mRNAs was constructed based on intergene co-expression correlations. Seven miRNAs associated with cell proliferation and apoptosis regulation were identified within this network, encompassing 92 subnet pairs as candidate genes for further exploration of exosome regulatory mechanisms. Additionally, preliminary verification at the cellular level demonstrated that exosomal miR-200b enhances the viability of POGCs. CONCLUSIONS: Transcriptome analysis unveiled a pivotal candidate ceRNA network potentially implicated in exosome-mediated regulation of granulosa cell proliferation and apoptosis, thereby influencing porcine follicular development. These findings offer insights into the molecular mechanisms of follicular fluid exosome regulation, encompassing both coding and non-coding RNA perspectives.

TRPM2-AS promotes ovarian cancer cell proliferation and inhibits cell apoptosis by upregulating the nearby gene TRPM2 via miR-6764-5p.

Ovarian cancer (OC) becomes a fatal gynecologic malignant cancer in females worldwide. Target therapy is a promising therapeutical choice for patients with OC, and identifying biomarkers and exploring molecular mechanisms are necessary. In this study, the functions and mechanism of long noncoding RNA transient receptor potential cation channel subfamily M member 2 antisense RNA (TRPM2-AS) in OC were explored. TRPM2-AS expression in OC cells was analyzed utilizing reverse transcription quantitative polymerase chain reaction (RT-qPCR). Cell counting kit-8 (CCK-8) and colony forming assays were carried out to explore the influence of TRPM2-AS on OC cell viability and proliferation. Cell apoptosis was detected using TdT-mediated dUTP Nick-End labeling (TUNEL) and flow cytometry analysis. Protein expression of apoptotic markers was subjected to western blotting. RNA pulldown or luciferase reporter assays were applied to explore the interaction between TRPM2-AS and miR-6764-5p or the binding of miR-6764-5p and TRPM2. The results showed that TRPM2-AS is highly expressed in OC cells and was mainly localized in cytoplasm. TRPM2-AS depletion suppressed OC cell viability and proliferation while increasing cell apoptotic rate. TRPM2 displayed a high level in OC cells and was positively regulated by TRPM2-AS. TRPM2-AS interacted with miR-6764-5p and thereby upregulated TRPM2 expression. In addition, TRPM2 overexpression reversed the repressive impact of TRPM2-AS depletion on malignant OC cellular process. In conclusion, TRPM2-AS promotes OC cell viability and proliferation while enhancing cell apoptosis through interaction with miR-6764-5p to regulate TRPM2 level.

Analysis of the miR482 Gene Family in Plants.

MicroRNA482 (miR482) is a conserved microRNA family in plants, playing critical regulatory roles in different biological activities. Though the members of the miR482 gene family have been identified in plants, a systematic study has not been reported yet. In the present research, 140 mature sequences generated by 106 precursors were used for molecular characterization, phylogenetic analysis, and target gene prediction, and the competing endogenous RNA (ceRNA) network mediated by miR482 was summarized. The length of mature sequences ranged from 17 nt to 25 nt, with 22 nt being the most abundant, and the start and end of the mature sequences had a preference for uracil (U). By sequence multiplex comparison, it was found that the mature sequences of 5p were clustered into one group, and others were clustered into the other group. Phylogenetic analysis revealed that the 140 mature sequences were categorized into six groups. Meanwhile, all the precursor sequences formed a stable hairpin structure, and the 106 precursors were divided into five groups. However, the expression of miR482 varied significantly between different species and tissues. In total, 149 target genes were predicted and their functions focused on single-organism process, cellular process, and cell and cell part. The ceRNA network of miR482 in tomato, cotton, and peanut was summarized based on related publications. In conclusion, this research will provide a foundation for further understanding of the miR482 gene family.

The role of lncRNAs related ceRNA regulatory network in multiple hippocampal pathological processes during the development of perioperative neurocognitive disorders.

Background: Perioperative neurocognitive disorders (PND) refer to neurocognitive abnormalities during perioperative period, which are a great challenge for elderly patients and associated with increased morbidity and mortality. Our studies showed that long non-coding RNAs (lncRNAs) regulate mitochondrial function and aging-related pathologies in the aged hippocampus after anesthesia, and lncRNAs are associated with multiple neurodegenerations. However, the regulatory role of lncRNAs in PND-related pathological processes remains unclear. Methods: A total of 18-month mice were assigned to control and surgery (PND) groups, mice in PND group received sevoflurane anesthesia and laparotomy. Cognitive function was assessed with fear conditioning test. Hippocampal RNAs were isolated for sequencing, lncRNA and microRNA libraries were constructed, mRNAs were identified, Gene Ontology (GO) analysis were performed, and lncRNA-microRNA-mRNA networks were established. qPCR was performed for gene expression verification. Results: A total of 312 differentially expressed (DE) lncRNAs, 340 DE-Transcripts of Uncertain Coding Potential (TUCPs), and 2,003 DEmRNAs were identified in the hippocampus between groups. The lncRNA-microRNA-mRNA competing endogenous RNA (ceRNA) network was constructed with 29 DElncRNAs, 90 microRNAs, 493 DEmRNAs, 148 lncRNA-microRNA interaction pairs, 794 microRNA-mRNA interaction pairs, and 110 lncRNA-mRNA co-expression pairs. 795 GO terms were obtained. Based on the frequencies of involved pathological processes, BP terms were divided into eight categories: neurological system alternation, neuronal development, metabolism alternation, immunity and neuroinflammation, apoptosis and autophagy, cellular communication, molecular modification, and behavior changes. LncRNA-microRNA-mRNA ceRNA networks in these pathological categories were constructed, and involved pathways and targeted genes were revealed. The top relevant lncRNAs in these ceRNA networks included RP23-65G6.4, RP24-396L14.1, RP23-251I16.2, XLOC_113622, RP24-496E14.1, etc., and the top relevant mRNAs in these ceRNA networks included Dlg4 (synaptic function), Avp (lipophagy), Islr2 (synaptic function), Hcrt (regulation of awake behavior), Tnc (neurotransmitter uptake). Conclusion: In summary, we have constructed the lncRNA-associated ceRNA network during PND development in mice, explored the role of lncRNAs in multiple pathological processes in the mouse hippocampus, and provided insights into the potential mechanisms and therapeutic gene targets for PND.

Core competing endogenous RNA network based on mRNA and non-coding RNA expression profiles in chicken fatty liver.

Fatty liver disease is a common metabolic disease in chickens. This disease can lead to a decrease in egg production and increase the risk of death in chickens. Long non-coding RNAs (lncRNAs) are involved in fatty liver formation by directly targeting genes or regulating gene expression by competitively binding microRNAs. However, a large proportion of competing endogenous RNA (ceRNA) networks in fatty liver diseases are still unclear. The total of 300 Jingxing-Huang chickens were used for fatty liver model construction. Then, differentially expressed (DE) genes (DEGs) identified through whole-transcriptome sequencing from four chickens with fatty liver and four chickens without fatty liver were chosen from the F1 generation. A total of 953 DEGs were identified between the fatty liver group and the control group, including 26 DE micro (mi)RNAs and 56 DE lncRNAs. Differential expression heatmaps and volcano plots were obtained after clustering expression analysis. Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that these DEGs were involved in many biological processes and signaling pathways related to fatty acid metabolism and lipid synthesis. Furthermore, cytoscape was used to construct a ceRNA network of the DE miRNAs, DE mRNAs, and DE lncRNAs. Eleven DE lncRNAs, seven DE miRNAs, and 13 DE mRNAs were found to be associated with the pathogenesis of fatty liver disease. An lncRNA-miRNA-mRNA ceRNA network was constructed to elucidate the mechanisms of fatty liver diseases, and the ENSGALT00000079786-miR-140/miR-143/miR-1a/miR-22/miR-375 network was identified. These results provide a valuable resource for further elucidating the posttranscriptional regulatory mechanisms of chicken liver and adipose fat development or deposition.

Revealing the novel metabolism-related genes in the ossification of the ligamentum flavum based on whole transcriptomic data.

Backgrounds: The ossification of the ligamentum flavum (OLF) is one of the major causes of thoracic myelopathy. Previous studies indicated there might be a potential link between metabolic disorder and pathogenesis of OLF. The aim of this study was to determine the potential role of metabolic disorder in the pathogenesis of OLF using the strict bioinformatic workflow for metabolism-related genes and experimental validation. Methods: A series of bioinformatic approaches based on metabolism-related genes were conducted to compare the metabolism score between OLF tissues and normal ligamentum flavum (LF) tissues using the single sample gene set enrichment analysis. The OLF-related and metabolism-related differentially expressed genes (OMDEGs) were screened out, and the biological functions of OMDEGs were explored, including the Gene Ontology enrichment analysis, Kyoto Encyclopedia of Genes and Genomes enrichment analysis, and protein-protein interaction. The competing endogenous RNA (ceRNA) network based on pairs of miRNA-hub OMDEGs was constructed. The correlation analysis was conducted to explore the potential relationship between metabolic disorder and immunity abnormality in OLF. In the end, the cell experiments were performed to validate the roles of GBE1 and TNF-α in the osteogenic differentiation of LF cells. Results: There was a significant difference of metabolism score between OLF tissues and normal LF tissues. Forty-nine OMDEGs were screened out and their biological functions were determined. The ceRNA network containing three hub OMDEGs and five differentially expressed miRNAs (DEmiRNAs) was built. The correlation analysis between hub OMDEGs and OLF-related infiltrating immune cells indicated that metabolic disorder might contribute to the OLF via altering the local immune status of LF tissues. The cell experiments determined the important roles of GBE1 expression and TNF-α in the osteogenic differentiation of LF cells. Conclusions: This research, for the first time, preliminarily illustrated the vital role of metabolic disorder in the pathogenesis of OLF using strict bioinformatic algorithms and experimental validation for metabolism-related genes, which could provide new insights for investigating disease mechanism and screening effective therapeutic targets of OLF in the future.

Exploration of ferroptosis and necroptosis-related genes and potential molecular mechanisms in psoriasis and atherosclerosis.

Objective: Ferroptosis and necroptosis are two recently identified forms of non-apoptotic cell death. Their dysregulation plays a critical role in the development and progression of Psoriasis (PsD) and Atherosclerosis (AS). This study explores shared Ferroptosis and necroptosis-related genes and elucidates their molecular mechanisms in PsD and AS through the analysis of public databases. Methods: Data sets for PsD (GSE30999) and AS (GSE28829) were retrieved from the GEO database. Differential gene expression (DEG) and weighted gene co-expression network analysis (WGCNA) were performed. Machine learning algorithms identified candidate biomarkers, whose diagnostic values were assessed using Receiver Operating Characteristic (ROC) curve analysis. Additionally, the expression levels of these biomarkers in cell models of AS and PsD were quantitatively measured using Western Blot (WB) and real-time quantitative PCR (RT-qPCR). Furthermore, CIBERSORT evaluated immune cell infiltration in PsD and AS tissues, highlighting the correlation between characteristic genes and immune cells. Predictive analysis for candidate drugs targeting characteristic genes was conducted using the DGIdb database, and an lncRNA-miRNA-mRNA network related to these genes was constructed. Results: We identified 44 differentially expressed ferroptosis-related genes (DE-FRGs) and 30 differentially expressed necroptosis-related genes (DE-NRGs). GO and KEGG enrichment analyses revealed significant enrichment of these genes in immune-related and inflammatory pathways, especially in NOD-like receptor and TNF signaling pathways. Two ferroptosis-related genes (NAMPT, ZFP36) and eight necroptosis-related genes (C7, CARD6, CASP1, CTSD, HMOX1, NOD2, PYCARD, TNFRSF21) showed high sensitivity and specificity in ROC curve analysis. These findings were corroborated in external validation datasets and cell models. Immune infiltration analysis revealed increased levels of T cells gamma delta, Macrophages M0, and Macrophages M2 in PsD and AS samples. Additionally, we identified 43 drugs targeting 5 characteristic genes. Notably, the XIST-miR-93-5p-ZFP36/HMOX1 and NEAT1-miR-93-5p-ZFP36/HMOX1 pathways have been identified as promising RNA regulatory pathways in AS and PsD. Conclusion: The two ferroptosis-related genes (NAMPT, ZFP36) and eight necroptosis-related genes (C7, CARD6, CASP1, CTSD, HMOX1, NOD2, PYCARD, TNFRSF21) are potential key biomarkers for PsD and AS. These genes significantly influence the pathogenesis of PsD and AS by modulating macrophage activity, participating in immune regulation, and mediating inflammatory responses.

Insights into the role of non-coding RNAs in the development of insecticide resistance in insects.

Pest control heavily relies on chemical pesticides has been going on for decades. However, the indiscriminate use of chemical pesticides often results in the development of resistance in pests. Almost all pests have developed some degree of resistance to pesticides. Research showed that the mechanisms of insecticide resistance in insects encompass metabolic resistance, behavioral resistance, penetration resistance and target-site resistance. Research on the these mechanisms has been mainly focused on the cis-regulatory or trans-regulatory for the insecticide resistance-related genes, with less attention paid to non-coding RNAs (ncRNAs), such as microRNA (miRNA), long non-coding RNA (lncRNA), and circular RNA (circRNA). There has been increased studies focus on understanding how these ncRNAs are involved in post-transcriptional regulation of insecticide resistance-related genes. Besides, the formatted endogenous RNA (ceRNA) regulatory networks (lncRNA/circRNA-miRNA-mRNA) has been identified as a key player in governing insect resistance formation. This review delves into the functions and underlying mechanisms of miRNA, lncRNA, and circRNA in regulating insect resistance. ncRNAs orchestrate insect resistance by modulating the expression of detoxification enzyme genes, insecticide target genes, as well as receptor genes, effectively regulating both target-site, metabolic and penetration resistance in insects. It also explores the regulatory mechanisms of ceRNA networks in the development of resistance. By enhancing our understanding of the mechanisms of ncRNAs in insecticide resistance, it will not only provide valuable insights into the new mechanisms of insecticide resistance but also help to enrich new directions in ncRNAs gene regulation research.

Identification of ceRNA networks in type H and L vascular endothelial cells through integrated bioinformatics methods. / 通过整合的生物信息学方法鉴定H和Låž‹è¡€ç®¡å† çš®ç»†èƒžä¸­çš„ceRNA网络（英文）.

OBJECTIVES: Type H blood vessels are a subtype of bone-specific microvessels (CD31hiEmcnhi) that play an important regulatory role in the coupling of angiogenesis and osteogenesis. Despite reports on the distinct roles of type H and L vessels under physiological and pathological bone conditions, their genetic differences remain to be elucidated. This study aims to construct a competitive endogenous RNA (ceRNA) network of key gene for differencial expression (DE) in type H and L vascular endothelial cells (ECs) through integrated bioinformatic methods. METHODS: We downloaded relevant raw data from the ArrayExpress and the Gene Expression Omnibus (GEO) database and used the Limma R-Bioconductor package to screen for DE lncRNAs, DE miRNAs, and DE mRNAs between type H and L vascular ECs. A total ceRNA network was constructed based on their interactions, followed by refinement using protein-protein interaction (PPI) networks to select upregulated and downregulated key genes. Enrichment analysis was performed on these key genes. Random validation was conducted using flow cytometry and real-time RT-PCR. RESULTS: A total of 1 761 DE mRNAs, 187 DE lncRNAs, and 159 DE miRNAs were identified, and a comprehensive ceRNA network was constructed based on their interactions. Six upregulated (Itga5, Kdr, Tjp1, Pecam1, Cdh5, and Ptk2) and 2 downregulated (Csf1r and Il10) key genes were selected via PPI network to construct a subnetwork of ceRNAs related to these key genes. Upregulated key genes were mainly enriched in negative regulation of angiogenesis and vascular apoptosis. Results from flow cytometry and real-time RT-PCR were consistent with bioinformatics analysis. CONCLUSIONS: This study proposes a ceRNA network associated with upregulated and downregulated type H and L vascular ECs based on selected key genes, providing new insights into the regulatory mechanisms of type H and L vascular ECs in bone metabolism.

A novel hsa_circ_0006903 circular RNA promotes tumor development and dendritic cells activated expression in infantile hemangioma.

Background: Increasing reports revealed that circular RNAs (circRNAs) and immune cells infiltration were related with tumor development. However, its role in infantile hemangioma (IH) is unknown. We will explore a novel hsa_circ_0006903-based ceRNA network and investigate the landscape of dendritic cells activated expression in IH. Material and methods: Differentially expressed circRNAs (DECs) were identified from Gene Expression Omnibus (GEO) database. Regulatory networks and functional enrichment analysis were constructed. CIBERSORT was used to characterize immune cells composition. qRT-PCR was performed to detect the expression of hsa_circ_0006903 in cell lines. Then, the role of hsa_circ_0006903 in IH were validated in vitro using transwell assay. Immunofluorescence was applied to the colocalization of CD11b for dendritic cells activated as a biomarker in IH tissues. Results: Using GEO database, a total of 67 DECs were screened out in IH. Hsa_circ_0006903 was the most significant DECs. Then, a novel hsa_circ_0006903 circular RNA-ceRNA network was constructed. Mechanistically, functional enrichment analysis showed that the p53 signaling pathway played the most important roles, and hsa_circ_0006903/miR-6721-5p/CACNA2D2 and hsa_circ_0006903/miR-4786-3p/ATP13A4 axis were identified. CACNA2D2, ATP13A4, and P53 were significantly downregulated in IH cell lines. We validated that dendritic cell activated was significantly overexpressed. Moreover, CD11b as a biomarker of dendritic cells activated were tested in IH tissues. Finally, hsa_circ_0006903 was significantly overexpressed, and hsa_circ_0006903 promoted infantile hemangioma cell proliferation, invasion, and migration in vitro. Conclusion: Overall, our study revealed that a novel hsa_circ_0006903 promoted tumor progression, and indicated a potential biomarker CD11b of dendritic cells activated in IH.