Sequence pre-training-based graph neural network for predicting lncRNA-miRNA associations.

MicroRNAs (miRNAs) silence genes by binding to messenger RNAs, whereas long non-coding RNAs (lncRNAs) act as competitive endogenous RNAs (ceRNAs) that can relieve miRNA silencing effects and upregulate target gene expression. The ceRNA association between lncRNAs and miRNAs has been a research hotspot due to its medical importance, but it is challenging to verify experimentally. In this paper, we propose a novel deep learning scheme, i.e. sequence pre-training-based graph neural network (SPGNN), that combines pre-training and fine-tuning stages to predict lncRNA-miRNA associations from RNA sequences and the existing interactions represented as a graph. First, we utilize a sequence-to-vector technique to generate pre-trained embeddings based on the sequences of all RNAs during the pre-training stage. In the fine-tuning stage, we use Graph Neural Network to learn node representations from the heterogeneous graph constructed using lncRNA-miRNA association information. We evaluate our proposed scheme SPGNN on our newly collected animal lncRNA-miRNA association dataset and demonstrate that combining the $k$-mer technique and Doc2vec model for pre-training with the Simple Graph Convolution Network for fine-tuning is effective in predicting lncRNA-miRNA associations. Our approach outperforms state-of-the-art baselines across various evaluation metrics. We also conduct an ablation study and hyperparameter analysis to verify the effectiveness of each component and parameter of our scheme. The complete code and dataset are available on GitHub: https://github.com/zixwang/SPGNN.

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.

LncRNA LUESCC promotes esophageal squamous cell carcinoma by targeting the miR-6785-5p/NRSN2 axis.

Esophageal squamous cell carcinoma (ESCC) is one of the most prevalent gastrointestinal malignancies with high mortality worldwide. Emerging evidence indicates that long noncoding RNAs (lncRNAs) are involved in human cancers, including ESCC. However, the detailed mechanisms of lncRNAs in the regulation of ESCC progression remain incompletely understood. LUESCC was upregulated in ESCC tissues compared with adjacent normal tissues, which was associated with gender, deep invasion, lymph node metastasis, and poor prognosis of ESCC patients. LUESCC was mainly localized in the cytoplasm of ESCC cells. Knockdown of LUESCC inhibited cell proliferation, colony formation, migration, and invasion in vitro and suppressed tumor growth in vivo. Mechanistic investigation indicated that LUESCC functions as a ceRNA by sponging miR-6785-5p to enhance NRSN2 expression, which is critical for the malignant behaviors of ESCC. Furthermore, ASO targeting LUESCC substantially suppressed ESCC both in vitro and in vivo. Collectively, these data demonstrate that LUESCC may exerts its oncogenic role by sponging miR-6785-5p to promote NRSN2 expression in ESCC, providing a potential diagnostic marker and therapeutic target for ESCC patients.

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.

Integrated analysis of lncRNA, miRNA and mRNA expression profiles reveals regulatory pathways associated with pig testis function.

Long noncoding RNA (lncRNA) and microRNA (miRNA) are known to play pivotal roles in mammalian testicular function and spermatogenesis. However, their impact on porcine male reproduction has yet to be well unraveled. Here, we sequenced and identified lncRNA and miRNA expressed in the testes of Chinese indigenous Banna mini-pig inbred line (BMI) and introduced Western Duroc (DU) and Large White (LW) pigs. By pairwise comparison (BMI vs DU, BMI vs LW, and DU vs LW), we found the gene expression differences in the testes between Chinese local pigs and introduced Western commercial breeds were more striking than those between introduced commercial breeds. Furthermore, we found 1622 co-differentially expressed genes (co-DEGs), 122 co-differentially expressed lncRNAs (co-DELs), 39 co-differentially expressed miRNAs (co-DEMs) in BMI vs introduced commercial breeds (DU and LW). Functional analysis revealed that these co-DEGs and co-DELs/co-DEMs target genes were enriched in male sexual function pathways, including MAPK, AMPK, TGF-ß/Smad, Hippo, NF-kappa B, and PI3K/Akt signaling pathways. Additionally, we established 10,536 lncRNA-mRNA, 11,248 miRNA-mRNA pairs, and 62 ceRNA (lncRNA-miRNA-mRNA) networks. The ssc-miR-1343 had the most interactive factors in the ceRNA network, including 20 mRNAs and 3 lncRNAs, consisting of 56 ceRNA pairs. These factors played extremely important roles in the regulation of testis function as key nodes in the interactive regulatory network. Our results provide insight into the functional roles of lncRNAs and miRNAs in porcine testis and offer a valuable resource for understanding the differences between Chinese indigenous and introduced Western pigs.

Identification of key regulatory molecules in the early development stage of Alzheimer's disease.

Alzheimer's disease (AD) is one of the most common neurodegenerative diseases, the incidence of which increases with age, and the pathological changes in the brain are irreversible. Recent studies have highlighted the essential role of long noncoding RNAs (lncRNAs) in AD by acting as competing endogenous RNAs (ceRNAs). Our aim was to construct lncRNA-associated ceRNA regulatory networks composed of potential biomarkers for the early stage of AD. AD related datasets come from AlzData and GEO databases. The R package 'Limma' identifies differentially expressed genes (DEGs), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases for functional enrichment analysis. Protein-protein interactions (PPIs) in DEGs were constructed in the STRING database, and Cytoscape software identified DEGs. Convergent functional genomics (CFG) analysis of differentially expressed hub genes (referred to as early-DEGs) in the brain before the development of AD pathology. The AlzData database analyses the expression levels of early-DEGs in different nerve cells. The lncRNA-miRNA-mRNA regulatory network was established according to the ceRNA hypothesis. We identified four lncRNAs (XIST, NEAT1, KCNQ1OT1 and HCG18) and four miRNAs (hsa-let-7c-5p, hsa-miR-107, hsa-miR-129-2-3p and hsa-miR-214-3p) were preliminarily identified as potential biomarkers for early AD, competitively regulating Atp6v0b, Atp6v1e1 Atp6v1f and Syt1. This study indicates that NEAT1, XIST, HCG18 and KCNQ1OT1 act as ceRNAs in competitive binding with miRNAs to regulate the expression of Atp6v0b, Atp6v1e1, Atp6v1f and Syt1 before the occurrence of pathological changes in AD.

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.

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.

Hepatocyte-derived exosomes deliver the lncRNA CYTOR to hepatic stellate cells and promote liver fibrosis.

Liver fibrosis is characterized by the activation and transformation of hepatic stellate cells (HSCs) induced by various injury factors. The degree of liver fibrosis can be significantly improved, but persistent injury factors present a significant therapeutic challenge. Hepatocytes are the most important parenchymal cell type in the liver. In this study, we explored the molecular mechanisms by which damaged liver cells activate HSCs through extracellular vesicles. We established a coculture model of LO2 and LX2 and validated its exosomal transmission activity. Subsequently, differentially expressed long noncoding RNAs (lncRNAs) were screened through RNA sequencing and their mechanisms of action as competing endogenous RNAs (ceRNAs) further confirmed using biological methods, such as FISH and luciferase assays. Damaged liver cells induced activation of LX2 and upregulation of liver fibrosis-related markers. Exosomes extracted and identified from the supernatant fraction contained differentially expressed lncRNA cytoskeleton regulator RNA (CYTOR) that competed with microRNA-125 (miR-125) for binding to glial cell line-derived neurotrophic factor (GDNF) in HSCs, in turn, promoting LX2 activation. MiR-125 could target and regulate both CYTOR and GDNF and vice versa, as verified using the luciferase assay. In an in vivo model, damaged liver extracellular vesicles induced the formation of liver fibrosis. Notably, downregulation of CYTOR within extracellular vesicles effectively inhibited liver fibrosis. The lncRNA CYTOR in exosomes of damaged liver cells is upregulated and modulates the expression of downstream GDNF through activity as a ceRNA, providing an effective mechanism for activation of HSCs.

Identification of necroptosis-related genes in ankylosing spondylitis by bioinformatics and experimental validation.

The pathogenesis of ankylosing spondylitis (AS) remains unclear, and while recent studies have implicated necroptosis in various autoimmune diseases, an investigation of its relationship with AS has not been reported. In this study, we utilized the Gene Expression Omnibus database to compare gene expressions between AS patients and healthy controls, identifying 18 differentially expressed necroptosis-related genes (DENRGs), with 8 upregulated and 10 downregulated. Through the application of three machine learning algorithms-least absolute shrinkage and selection operation, support vector machine-recursive feature elimination and random forest-two hub genes, FASLG and TARDBP, were pinpointed. These genes demonstrated high specificity and sensitivity for AS diagnosis, as evidenced by receiver operating characteristic curve analysis. These findings were further supported by external datasets and cellular experiments, which confirmed the downregulation of FASLG and upregulation of TARDBP in AS patients. Immune cell infiltration analysis suggested that CD4+ T cells, CD8+ T cells, NK cells and neutrophils may be associated with the development of AS. Notably, in the group with high FASLG expression, there was a significant infiltration of CD8+ T cells, memory-activated CD4+ T cells and resting NK cells, with relatively less infiltration of memory-resting CD4+ T cells and neutrophils. Conversely, in the group with high TARDBP expression, there was enhanced infiltration of naïve CD4+ T cells and M0 macrophages, with a reduced presence of memory-resting CD4+ T cells. In summary, FASLG and TARDBP may contribute to AS pathogenesis by regulating the immune microenvironment and immune-related signalling pathways. These findings offer new insights into the molecular mechanisms of AS and suggest potential new targets for therapeutic strategies.

Improved therapeutic effects on vascular intimal hyperplasia by mesenchymal stem cells expressing MIR155HG that function as a ceRNA for microRNA-205.

Coronary artery bypass grafting (CABG) is an effective treatment for coronary heart disease, with vascular transplantation as the key procedure. Intimal hyperplasia (IH) gradually leads to vascular stenosis, seriously affecting the curative effect of CABG. Mesenchymal stem cells (MSCs) were used to alleviate IH, but the effect was not satisfactory. This work aimed to investigate whether lncRNA MIR155HG could improve the efficacy of MSCs in the treatment of IH and to elucidate the role of the competing endogenous RNA (ceRNA). The effect of MIR155HG on MSCs function was investigated, while the proteins involved were assessed. IH was detected by HE and Van Gieson staining. miRNAs as the target of lncRNA were selected by bioinformatics analysis. qRT-PCR and dual-luciferase reporter assay were performed to verify the binding sites of lncRNA-miRNA. The apoptosis, Elisa and tube formation assay revealed the effect of ceRNA on the endothelial protection of MIR155HG-MSCs. We observed that MIR155HG improved the effect of MSCs on IH by promoting viability and migration. MIR155HG worked as a sponge for miR-205. MIR155HG/miR-205 significantly improved the function of MSCs, avoiding apoptosis and inducing angiogenesis. The improved therapeutic effects of MSCs on IH might be due to the ceRNA role of MIR155HG/miR-205.

Long noncoding RNA Glis2 regulates podocyte mitochondrial dysfunction and apoptosis in diabetic nephropathy via sponging miR-328-5p.

Podocyte apoptosis exerts a crucial role in the pathogenesis of DN. Recently, long noncoding RNAs (lncRNAs) have been gradually identified to be functional in a variety of different mechanisms associated with podocyte apoptosis. This study aimed to investigate whether lncRNA Glis2 could regulate podocyte apoptosis in DN and uncover the underlying mechanism. The apoptosis rate was detected by flow cytometry. Mitochondrial membrane potential (ΔΨM) was measured using JC-1 staining. Mitochondrial morphology was detected by MitoTracker Deep Red staining. Then, the histopathological and ultrastructure changes of renal tissues in diabetic mice were observed using periodic acid-Schiff (PAS) staining and transmission electron microscopy. We found that lncRNA Glis2 was significantly downregulated in high-glucose cultured podocytes and renal tissues of db/db mice. LncRNA Glis2 overexpression was found to alleviate podocyte mitochondrial dysfunction and apoptosis. The direct interaction between lncRNA Glis2 and miR-328-5p was confirmed by dual luciferase reporter assay. Furthermore, lncRNA Glis2 overexpression alleviated podocyte apoptosis in diabetic mice. Taken together, this study demonstrated that lncRNA Glis2, acting as a competing endogenous RNA (ceRNA) of miRNA-328-5p, regulated Sirt1-mediated mitochondrial dysfunction and podocyte apoptosis in DN.

Comprehensive analysis of the expression profiles of mRNA, lncRNA, circRNA, and miRNA in primary hair follicles of coarse sheep fetal skin.

BACKGROUND: The Qinghai Tibetan sheep, a local breed renowned for its long hair, has experienced significant deterioration in wool characteristics due to the absence of systematic breeding practices. Therefore, it is imperative to investigate the molecular mechanisms underlying follicle development in order to genetically enhance wool-related traits and safeguard the sustainable utilization of valuable germplasm resources. However, our understanding of the regulatory roles played by coding and non-coding RNAs in hair follicle development remains largely elusive. RESULTS: A total of 20,874 mRNAs, 25,831 circRNAs, 4087 lncRNAs, and 794 miRNAs were annotated. Among them, we identified 58 DE lncRNAs, 325 DE circRNAs, 924 DE mRNAs, and 228 DE miRNAs during the development of medullary primary hair follicle development. GO and KEGG functional enrichment analyses revealed that the JAK-STAT, TGF-ß, Hedgehog, PPAR, cGMP-PKG signaling pathway play crucial roles in regulating fibroblast and epithelial development during skin and hair follicle induction. Furthermore, the interactive network analysis additionally identified several crucial mRNA, circRNA, and lncRNA molecules associated with the process of primary hair follicle development. Ultimately, by investigating DEmir's role in the ceRNA regulatory network mechanism, we identified 113 circRNA-miRNA pairs and 14 miRNA-mRNA pairs, including IGF2BP1-miR-23-x-novel-circ-01998-MSTRG.7111.3, DPT-miR-370-y-novel-circ-005802-MSTRG.14857.1 and TSPEAR-oar-miR-370-3p-novel-circ-005802- MSTRG.10527.1. CONCLUSIONS: Our study offers novel insights into the distinct expression patterns of various transcription types during hair follicle morphogenesis, establishing a solid foundation for unraveling the molecular mechanisms that drive hair development and providing a scientific basis for selectively breeding desirable wool-related traits in this specific breed.

Identification of different myofiber types in pigs muscles and construction of regulatory networks.

BACKGROUND: Skeletal muscle is composed of muscle fibers with different physiological characteristics, which plays an important role in regulating skeletal muscle metabolism, movement and body homeostasis. The type of skeletal muscle fiber directly affects meat quality. However, the transcriptome and gene interactions between different types of muscle fibers are not well understood. RESULTS: In this paper, we selected 180-days-old Large White pigs and found that longissimus dorsi (LD) muscle was dominated by fast-fermenting myofibrils and soleus (SOL) muscle was dominated by slow-oxidizing myofibrils by frozen sections and related mRNA and protein assays. Here, we selected LD muscle and SOL muscle for transcriptomic sequencing, and identified 312 differentially expressed mRNA (DEmRs), 30 differentially expressed miRNA (DEmiRs), 183 differentially expressed lncRNA (DElRs), and 3417 differentially expressed circRNA (DEcRs). The ceRNA network included ssc-miR-378, ssc-miR-378b-3p, ssc-miR-24-3p, XR_308817, XR_308823, SMIM8, MAVS and FOS as multiple core nodes that play important roles in muscle development. Moreover, we found that different members of the miR-10 family expressed differently in oxidized and glycolytic muscle fibers, among which miR-10a-5p was highly expressed in glycolytic muscle fibers (LD) and could target MYBPH gene mRNA. Therefore, we speculate that miR-10a-5p may be involved in the transformation of muscle fiber types by targeting the MYHBP gene. In addition, PPI analysis of differentially expressed mRNA genes showed that ACTC1, ACTG2 and ACTN2 gene had the highest node degree, suggesting that this gene may play a key role in the regulatory network of muscle fiber type determination. CONCLUSIONS: We can conclude that these genes play a key role in regulating muscle fiber type transformation. Our study provides transcriptomic profiles and ceRNA interaction networks for different muscle fiber types in pigs, providing reference for the transformation of pig muscle fiber types and the improvement of meat quality.

Transcriptome reveals the roles and potential mechanisms of lncRNAs in the regulation of albendazole resistance in Haemonchus contortus.

BACKGROUND: Haemonchus contortus (H. contortus) is the most common parasitic nematode in ruminants and is prevalent worldwide. H. contortus resistance to albendazole (ABZ) hinders the efficacy of anthelmintic drugs, but little is known about the molecular mechanisms that regulate this of drug resistance. Recent research has demonstrated that long noncoding RNAs (lncRNAs) can exert significant influence as pivotal regulators of the emergence of drug resistance. RESULTS: In this study, transcriptome sequencing was conducted on both albendazole-sensitive (ABZ-sensitive) and albendazole-resistant (ABZ-resistant) H. contortus strains, with three biological replicates for each group. The analysis of lncRNA in the transcriptomic data revealed that there were 276 differentially expressed lncRNA (DElncRNA) between strains with ABZ-sensitive and ABZ-resistant according to the criteria of |log2Foldchange|≥ 1 and FDR < 0.05. Notably, MSTRG.12969.2 and MSTRG.9827.1 exhibited the most significant upregulation and downregulation, respectively, in the resistant strains. The potential roles of the DElncRNAs included catalytic activity, stimulus response, regulation of drug metabolism, and modulation of the immune response. Moreover, we investigated the interactions between DElncRNAs and other RNAs, specifically MSTRG.12741.1, MSTRG.11848.1, MSTRG.5895.1, and MSTRG.14070.1, involved in regulating drug stimulation through cis/trans/antisense/lncRNA‒miRNA-mRNA interaction networks. This regulation leads to a decrease (or increase) in the expression of relevant genes, consequently enhancing the resistance of H. contortus to albendazole. Furthermore, through comprehensive analysis of competitive endogenous RNAs (ceRNAs) involved in drug resistance-related pathways, such as the mTOR signalling pathway and ABC transporter signalling pathway, the relevance of the MSTRG.2499.1-novel-m0062-3p-HCON_00099610 interaction was identified to mainly involve the regulation of catalytic activity, metabolism, ubiquitination and transcriptional regulation of gene promoters. Additionally, quantitative real-time polymerase chain reaction (qRT-PCR) validation indicated that the transcription profiles of six DElncRNAs and six DEmRNAs were consistent with those obtained by RNA-seq. CONCLUSIONS: The results of the present study allowed us to better understand the changes in the lncRNA expression profile of ABZ-resistant H. contortus. In total, these results suggest that the lncRNAs MSTRG.963.1, MSTRG.12741.1, MSTRG.11848.1 and MSTRG.2499.1 play important roles in the development of ABZ resistance and can serve as promising biomarkers for further study.

Construction of LncRNA-Related ceRNA Networks in Longissimus Dorsi Muscle of Jinfen White Pigs at Different Developmental Stages.

The development of skeletal muscle in pigs might determine the quality of pork. In recent years, long non-coding RNAs (lncRNAs) have been found to play an important role in skeletal muscle growth and development. In this study, we investigated the whole transcriptome of the longissimus dorsi muscle (LDM) of Jinfen White pigs at three developmental stages (1, 90, and 180 days) and performed a comprehensive analysis of lncRNAs, mRNAs, and micro-RNAs (miRNAs), aiming to find the key regulators and interaction networks in Jinfen White pigs. A total of 2638 differentially expressed mRNAs (DE mRNAs) and 982 differentially expressed lncRNAs (DE lncRNAs) were identified. Compared with JFW_1d, there were 497 up-regulated and 698 down-regulated DE mRNAs and 212 up-regulated and 286 down-regulated DE lncRNAs in JFW_90d, respectively. In JFW_180d, there were 613 up-regulated and 895 down-regulated DE mRNAs and 184 up-regulated and 131 down-regulated DE lncRNAs compared with JFW_1d. There were 615 up-regulated and 477 down-regulated DE mRNAs and 254 up-regulated and 355 down-regulated DE lncRNAs in JFW_180d compared with JFW_90d. Compared with mRNA, lncRNA has fewer exons, fewer ORFs, and a shorter length. We performed GO and KEGG pathway functional enrichment analysis for DE mRNAs and the potential target genes of DE lncRNAs. As a result, several pathways are involved in muscle growth and development, such as the PI3K-Akt, MAPK, hedgehog, and hippo signaling pathways. These are among the pathways through which mRNA and lncRNAs function. As part of this study, bioinformatic screening was used to identify miRNAs and DE lncRNAs that could act as ceRNAs. Finally, we constructed an lncRNA-miRNA-mRNA regulation network containing 26 mRNAs, 7 miRNAs, and 17 lncRNAs; qRT-PCR was used to verify the key genes in these networks. Among these, XLOC_022984/miR-127/ENAH and XLOC_016847/miR-486/NRF1 may function as key ceRNA networks. In this study, we obtained transcriptomic profiles from the LDM of Jinfen White pigs at three developmental stages and screened out lncRNA-miRNA-mRNA regulatory networks that may provide crucial information for the further exploration of the molecular mechanisms during skeletal muscle development.

Circular RNA circESYT2 serves as a microRNA-665 sponge to promote the progression of hepatocellular carcinoma through ENO2.

Circular RNAs (circRNAs) have emerged as crucial regulators in tumor progression, yet their specific role in hepatocellular carcinoma (HCC) remains largely uncharacterized. In this study, we utilized high-transcriptome sequencing to identify the upregulation of circESYT2 (hsa_circ_002142) in HCC tissues. Functional experiments carried out in vivo and in vitro revealed that circESYT2 played a significant role in maintaining the growth and metastatic behaviors of HCC. Through integrative analysis, we identified enolase 2 (ENO2) as a potential target regulated by circESYT2 through the competitive endogenous RNA sponge mechanism. Additional gain- or loss-of-function experiments indicated that overexpression of circESYT2 led to a tumor-promoting effect, which could be reversed by transfection of microRNA-665 (miR-665) mimic or ENO2 knockdown in HCC cells. Furthermore, the direct interaction between miR-665 and circESYT2 and between miR-665 and ENO2 was confirmed using RNA immunoprecipitation, FISH, RNA pull-down, and dual-luciferase reporter assays, highlighting the involvement of the circESYT2/miR-665/ENO2 axis in promoting HCC progression. These findings shed light on the molecular characteristics of circESYT2 in HCC tissues and suggest its potential as a biomarker or therapeutic target for HCC treatment.

LINC02086 inhibits ferroptosis and promotes malignant phenotypes of pancreatic cancer via miR-342-3p/CA9 axis.

Long noncoding RNAs (lncRNAs) play important roles in modulating the tumorigenesis and progression of malignant tumors. LINC02086 is a newly identified oncogene associated with tumorigenesis, but its role in pancreatic cancer (PC) has not been fully elucidated. In this study we examined the expression levels of LINC02086, miR-342-3p, and CA9 in PC. The relationship of ferroptosis with these factors was analyzed by detecting the expression levels of Fe2+, reactive oxygen species (ROS), and ferroptosis marker proteins. The expression of these genes was altered to observe their effects on cell proliferation, migration, and invasion ability. Bioinformatics was used to predict target genes, and the binding relationship was verified luciferase reporter assay. Finally, the function of LINC02086 was evaluated in vivo. The findings suggest that LINC02086 is highly expressed in PC tissues and cell lines and is correlated with a poor prognosis. In vitro experiments demonstrated that LINC02086 knockdown promoted ferroptosis in PC cells to suppress their malignant phenotype. LINC02086 acts as a competitive endogenous RNA that adsorbed miR-342-3p. miR-342-3p hinders the malignant progression of PC by promoting ferroptosis. In addition, miR-342-3p targets CA9 and affects its function. Further mechanistic studies revealed that LINC02086 inhibits ferroptosis and promotes PC progression by acting as a sponge for miR-342-3p to upregulate CA9 expression. In vivo experiments further confirmed this mechanism. Taken together, LINC02086 upregulates CA9 expression by competitively binding with miR-342-3p, thereby inhibiting ferroptosis in PC cells and promoting their malignant phenotype. The results of our study provide new insights into how LINC02086 contributes to the progression of PC.

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.

Identification and characterization of circular RNAs involved in the fertility stability of cotton CMS-D2 restorer line under heat stress.

BACKGROUND: As a vital type of noncoding RNAs, circular RNAs (circRNAs) play important roles in plant growth and development and stress response. However, little is known about the biological roles of circRNAs in regulating the stability of male fertility restoration for cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) cotton under high-temperature (HT) stress. RESULTS: In this study, RNA-sequencing and bioinformatics analysis were performed on pollen grains of isonuclear alloplasmic near-isogenic restorer lines NH [N(Rf1rf1)] and SH [S(Rf1rf1)] with obvious differences in fertility stability under HT stress at two environments. A total of 967 circRNAs were identified, with 250 differentially expressed under HT stress. We confirmed the back-splicing sites of eight selected circRNAs using divergent primers and Sanger sequencing. Tissue-specific expression patterns of five differentially expressed circRNAs (DECs) were also verified by RT-PCR and qRT-PCR. Functional enrichment and metabolic pathway analysis revealed that the parental genes of DECs were significantly enriched in fertility-related biological processes such as pollen tube guidance and cell wall organization, as well as the Pentose and glucuronate interconversions, Steroid biosynthesis, and N-Glycan biosynthesis pathways. Moreover, we also constructed a putative circRNA-mediated competing endogenous RNA (ceRNA) network consisting of 21 DECs, eight predicted circRNA-binding miRNAs, and their corresponding 22 mRNA targets, especially the two ceRNA modules circRNA346-miR159a-MYB33 and circRNA484-miR319e-MYB33, which might play important biological roles in regulating pollen fertility stability of cotton CMS-D2 restorer line under HT stress. CONCLUSIONS: Through systematic analysis of the abundance, characteristics and expression patterns of circRNAs, as well as the potential functions of their parent genes, our findings suggested that circRNAs and their mediated ceRNA networks acted vital biological roles in cotton pollen development, and might be also essential regulators for fertility stability of CMS-D2 restorer line under heat stress. This study will open a new door for further unlocking complex regulatory mechanisms underpinning the fertility restoration stability for CMS-D2 in cotton.