mRNA Interactions Research Articles

The Parasite-Derived Peptide, FhHDM-1, Selectively Modulates miRNA Expression in β-Cells to Prevent Apoptotic Pathways Induced by Proinflammatory Cytokines.

We have previously identified a parasite-derived peptide, FhHDM-1, that prevented the progression of diabetes in nonobese diabetic (NOD) mice. Disease prevention was mediated by the activation of the PI3K/Akt pathway to promote β-cell survival and metabolism without inducing proliferation. To determine the molecular mechanisms driving the antidiabetogenic effects of FhHDM-1, miRNA:mRNA interactions and in silico predictions of the gene networks were characterised in β-cells, which were exposed to the proinflammatory cytokines that mediate β-cell destruction in Type 1 diabetes (T1D), in the presence and absence of FhHDM-1. The predicted gene targets of miRNAs differentially regulated by FhHDM-1 mapped to the biological pathways that regulate β-cell biology. Six miRNAs were identified as important nodes in the regulation of PI3K/Akt signaling. Additionally, IGF-2 was identified as a miRNA gene target that mediated the beneficial effects of FhHDM-1 on β-cells. The findings provide a putative mechanism by which FhHDM-1 positively impacts β-cells to permanently prevent diabetes. As β-cell death/dysfunction underlies diabetes development, FhHDM-1 opens new therapeutic avenues.

IGF2BP3 regulates the expression of RRM2 and promotes the progression of rheumatoid arthritis via RRM2/Akt/MMP-9 pathway.

Rheumatoid arthritis (RA) is a common inflammatory and autoimmune disease. Ribonucleotide Reductase Regulatory Subunit M2 (RRM2) is a crucial and a rate-limiting enzyme responsible for deoxynucleotide triphosphate(dNTP) production. We have found a high expression level of RRM2 in patients with RA, but the molecular mechanism of its action remains unclear. We analyzed the expression of hub genes in RA using GSE77298 datasets downloaded from Gene Expression Omnibus database. RRM2 and insulin-like growth factor-2 messenger ribonucleic acid (mRNA)-binding protein 3 (IGF2BP3) gene knockdown was achieved by infection with lentiviruses. The expression of RRM2, IGF2BP3, matrix metalloproteinase (MMP)-1, and MMP-9 were detected via western blotting assay. Cell viability was detected via 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. MeRIP-qRT-PCR was performed to test the interaction of IGF2BP3 and RRM2 mRNA via m6A modification. Cell proliferation was determined by clone formation assay. Migration and invasion assays were performed using transwell Boyden chamber. RRM2 and IGF2BP3 were highly expressed in clinical specimens and tumor necrosis factor alpha (TNF-α) and interleukin (IL)-1β-stimulated synovial cells. RRM2 and IGF2BP3 knockdown inhibited the proliferation, migration, and invasion of MH7A cells. The inhibitory effects of IGF2BP3 knockdown were effectively reversed by simultaneously overexpressing RRM2 in MH7A cells. By analyzing N6-methyladenosine (m6A)2Target database, five m6A regulatory target binding sites for IGF2BP3 were identified in RRM2 mRNA, suggesting a direct relationship between IGF2BP3 and RRM2 mRNA. Additionally, in RRM2 small hairpin (sh)RNA lentivirus-infected cells, the levels of phosphorylated Akt and MMP-9 were significantly decreased compared with control shRNA lentivirus-infected cells. The present study demonstrated that RRM2 promoted the Akt phosphorylation leading to high expression of MMP-9 to promote the migration and invasive capacities of MH7A cells. Overall, IGF2BP promotes the expression of RRM2, and regulates the migration and invasion of MH7A cells via Akt/MMP-9 pathway to promote RA progression.

Regulation of Oocyte mRNA Metabolism: A Key Determinant of Oocyte Developmental Competence.

The regulation of mRNA transcription and translation is uncoupled during oogenesis. The reason for this uncoupling is two-fold. Chromatin is only accessible to the transcriptional machinery during the growth phase as it condenses prior to resumption of meiosis to ensure faithful segregation of chromosomes during meiotic maturation. Thus, transcription rates are high during this time period in order to produce all of the transcripts needed for meiosis, fertilization, and embryo cleavage until the newly formed embryonic genome becomes transcriptionally active. To ensure appropriate timing of key developmental milestones including chromatin condensation, resumption of meiosis, segregation of chromosomes, and polar body extrusion, the translation of protein from transcripts synthesized during oocyte growth must be temporally regulated. This is achieved by the regulation of mRNA interaction with RNA binding proteins and shortening and lengthening of the poly(A) tail. This chapter details the essential factors that regulate the dynamic changes in mRNA synthesis, storage, translation, and degradation during oocyte growth and maturation.

Implications for miR-339-5p regulation of trophoblast proliferation and migration in placentas associated with porcine intrauterine growth retardation using integrated transcriptome sequencing analysis

Placental dysfunction is considered as one of the main etiologies of fetal intrauterine growth retardation (IUGR). MicroRNAs (miRNAs) have been demonstrated to be a vital epigenetic modification involved in regulating the placental function and pregnancy outcomes in mammals. However, the mechanisms underlying placenta-specific miRNAs involved in the occurrence and development of pig IUGR remain unclear. In this work, we compared the placental morphologies of piglets with IUGR and normal birth weight (NBW) by using histomorphological analysis and performed a miRNA–mRNA integrative analysis of the gene expression profiles of IUGR and NBW placentas through RNA sequencing. We also investigated the role of differentially expressed ssc-miR-339-5p/GRIK3 through an in vitro experiment on porcine trophoblast cells (PTr2). IUGR piglets had significantly lower birth weight, placental weight, placental efficiency, and placental villus and capillary densities compared with the NBW piglets (P < 0.05). A total of 81 differentially expressed miRNAs and 726 differentially expressed genes in the placentas were screened out between the IUGR and NBW groups. The miRNA–mRNA interaction networks revealed the key core miRNA (ssc-miR-339-5p) and its corresponding target genes. Subsequently, we found that upregulation of ssc-miR-339-5p significantly inhibited the migration and proliferation of PTr2 cells (P < 0.05). The dual-luciferase reporter system showed that GRIK3 was the target gene of ssc-miR-339-5p, and the transcription level of GRIK3 may be negatively regulated by ssc-miR-339-5p. Additionally, overexpression of ssc-miR-339-5p significantly increased (P < 0.05) the mRNA expression levels of genes involved in the cytokine–cytokine receptor interaction pathway. These results indicate that ssc-miR-339-5p may affect the migration and proliferation of trophoblast cells by regulating the expression of GRIK3 and altering the placental inflammatory response, resulting in a suboptimal morphology and function of the placenta and the development of pig IUGR.

Genome-Wide Identification of Cotton MicroRNAs Predicted for Targeting Cotton Leaf Curl Kokhran Virus-Lucknow

Cotton leaf curl Kokhran virus (CLCuKoV) (genus, Begomovirus; family, Geminiviridae) is one of several plant virus pathogens of cotton (Gossypium hirsutum L.) that cause cotton leaf curl disease in Pakistan. Begomoviruses are transmitted by the whitefly Bemisia tabaci cryptic species group and cause economic losses in cotton and other crops worldwide. The CLCuKoV strain, referred to as CLCuKoV-Bur, emerged in the vicinity of Burewala, Pakistan, and was the primary causal virus associated with the second CLCuD epidemic in Pakistan. The monopartite ssDNA genome of (2.7 Kb) contains six open reading frames that encode four predicted proteins. RNA interference (RNAi)-mediated antiviral immunity is a sequence-specific biological process in plants and animals that has evolved to combat virus infection. The objective of this study was to design cotton locus-derived microRNA (ghr-miRNA) molecules to target strains of CLCuKoV, with CLCuKoV-Lu, as a typical CLCuD-begomovirus genome, predicted by four algorithms, miRanda, RNA22, psRNATarget, and RNA hybrid. Mature ghr-miRNA sequences (n = 80) from upland cotton (2n = 4x = 52) were selected from miRBase and aligned with available CLCuKoV-Lu genome sequences. Among the 80 cotton locus-derived ghr-miRNAs analyzed, ghr-miR2950 was identified as the most optimal, effective ghr-miRNA for targeting the CLCuKoV-Lu genome (nucleotide 82 onward), respectively, based on stringent criteria. The miRNA targeting relies on the base pairing of miRNA–mRNA targets. Conservation and potential base pairing of binding sites with the ghr-miR2950 were validated by multiple sequence alignment with all available CLCuKoV sequences. A regulatory interaction network was constructed to evaluate potential miRNA–mRNA interactions with the predicted targets. The efficacy of miRNA targeting of CLCuKoV was evaluated in silico by RNAi-mediated mRNA cleavage. This predicted targets for the development of CLCuD-resistant cotton plants.

Circular RNA RHBDD1 regulates tumorigenicity and ferroptosis in colorectal cancer by mediating the ELAVL1/SCD mRNA interaction.

Circular RNAs (circRNAs) are covalently closed noncoding RNA molecules that play multiple roles in tumorigenesis and metastasis. Ferroptosis is an iron-dependent, regulated form of cell death and has emerged as a promising target for cancer treatment. However, whether and how circRNAs regulate ferroptotic cell death in colorectal cancer (CRC) remains largely unknown. Three circRNA microarrays were used to screen differentially expressed circRNAs in CRC tissues. A series of functional experiments were conducted to investigate the effects of circRNA on CRC cell proliferation, migration and ferroptosis. We found that hsa_circ_0058495 (circRHBDD1), a novel circRNA, was significantly upregulated in colorectal cancer tissues and cells. The expression levels of circRHBDD1 in serum samples were strongly associated with the advancement of CRC. Silencing of circRHBDD1 remarkably suppressed the proliferation and migration of CRC cells in vitro. Moreover, the depletion of circRHBDD1 notably increased ferroptotic cell death and enhanced RSL3-induced ferroptosis in CRC cells. Mechanistically, circRHBDD1 upregulated the expression of stearoyl-CoA desaturase (SCD), a ferroptosis suppressor mediating lipid remodelling, by enhancing the ELAVL1/SCD mRNA interaction. Finally, circRHBDD1 knockdown repressed the tumorigenesis and ferroptosis of CRC cells in vivo. In conclusion, circRHBDD1 facilitates tumour progression and obstructs ferroptosis in CRC by regulating SCD expression in an ELAVL1-dependent manner.

Establishing mRNA and microRNA interactions driving disease heterogeneity in amyotrophic lateral sclerosis patient survival.

Amyotrophic lateral sclerosis is a fatal neurodegenerative disease, associated with the degeneration of both upper and lower motor neurons of the motor cortex, brainstem and spinal cord. Death in most patients results from respiratory failure within 3-4 years from symptom onset. However, due to disease heterogeneity some individuals survive only months from symptom onset while others live for several years. Identifying specific biomarkers that aid in establishing disease prognosis, particularly in terms of predicting disease progression, will help our understanding of amyotrophic lateral sclerosis pathophysiology and could be used to monitor a patient's response to drugs and therapeutic agents. Transcriptomic profiling technologies are continually evolving, enabling us to identify key gene changes in biological processes associated with disease. MicroRNAs are small non-coding RNAs typically associated with regulating gene expression, by degrading mRNA or reducing levels of gene expression. Being able to associate gene expression changes with corresponding microRNA changes would help to distinguish a more complex biomarker signature enabling us to address key challenges associated with complex diseases such as amyotrophic lateral sclerosis. The present study aimed to investigate the transcriptomic profile (mRNA and microRNA) of lymphoblastoid cell lines from amyotrophic lateral sclerosis patients to identify key signatures that are distinguishable in those patients who suffered a short disease duration (<12 months) (n = 22) compared with those that had a longer disease duration (>6 years) (n = 20). Transcriptional profiling of microRNA-mRNA interactions from lymphoblastoid cell lines in amyotrophic lateral sclerosis patients revealed differential expression of genes involved in cell cycle, DNA damage and RNA processing in patients with longer survival from disease onset compared with those with short survival. Understanding these particular microRNA-mRNA interactions and the pathways in which they are involved may help to distinguish potential therapeutic targets that could exert neuroprotective effects to prolong the life expectancy of amyotrophic lateral sclerosis patients.

Interaction of mRNA with the C-Terminal Domain of PCID2, a Subunit of the TREX-2 Complex, Is Required for Its Export from the Nucleus to the Cytoplasm in Drosophila melanogaster

Following the transcription step, the newly synthesized mRNA is exported from the nucleus to the cytoplasm and further to the translation site. The TREX-2 complex is involved in the step of mRNA export from the nucleus to the cytoplasm. This complex in Drosophila melanogaster consists of four proteins: Xmas-2, PCID2, ENY2, and Sem1p. In our work, we have shown that deletion of the C-terminal sequence of PCID2 leads to a decrease in the interaction of the protein with RNA and to impaired mRNA export from the nucleus to the cytoplasm in D. melanogaster.

Transcriptome Analysis of mRNA and lncRNA Related to Muscle Growth and Development in Gannan Yak and Jeryak

The production performance of Jeryak, resulting from the F1 generation of the cross between Gannan yak and Jersey cattle, exhibits a significantly superior outcome compared with that of Gannan yak. Therefore, we used an RNA-seq approach to identify differentially expressed mRNAs (DEMs) and differentially expressed lncRNAs (DELs) influencing muscle growth and development in Gannan yaks and Jeryaks. A total of 304 differentially expressed lncRNAs and 1819 differentially expressed mRNAs were identified based on the screening criteria of |log 2 FC| &gt; 1 and FDR &lt; 0.05. Among these, 132 lncRNAs and 1081 mRNAs were found to be down-regulated, while 172 lncRNAs and 738 mRNAs were up-regulated. GO and KEGG analyses showed that the identified DELs and DEMs were enriched in the entries of pathways associated with muscle growth and development. On this basis, we constructed an lncRNA–mRNA interaction network. Interestingly, two candidate DELs (MSTRG.16260.9 and MSTRG.22127.1) had targeting relationships with 16 (MYC, IGFBP5, IGFBP2, MYH4, FGF6, etc.) genes related to muscle growth and development. These results could provide a basis for further studies on the roles of lncRNAs and mRNAs in muscle growth in Gannan yaks and Jeryak breeds.

Identification of Candidate mRNA Genes and Their Potential MicroRNA Targets in Lung Cancer Induced by Smoking Tobacco.

Smoking is considered the single highest risk factor for lung cancer and has been suggested to be associated with accelerated somatic mutations in respiratory mucosa that lead to the development of lung cancer. MicroRNAs serve as modulators in smoking-induced mRNA gene expression changes in the human airway epithelium and are linked to the development of lung cancer. The thermodynamics in the microRNA (miRNA)-mRNA interactions may be affected in tobacco smokers, consequently, leading to phenotypic variations in lung cancer patients. Therefore, this study aimed to investigate the impact of smoking tobacco on somatic mutations in mRNA genes and assess their potential impact on miRNA-mRNA interactions in lung cancers. The clinically significant pathogenic variants in mRNA genes in the dataset in lung cancer cases linked to smoking tobacco (n = 330) were obtained from the Cancer Atlas database (TCGA, http://cancergenome.nih.gov/) and used to assess the potential role of tobacco consumption in driving the genetic alterations in proto-oncogenes associated with lung cancer. The analysis of the miRNA interaction with the top five altered mRNA proto-oncogenes in lung cancer cases due to tobacco consumption was performed using the target prediction function in the miRDP program (Database version 5.2.3.1, https://mirdb.org/). We identified the top five mRNA proto-oncogenes enriched with simple somatic mutations (SSM) in lung cancer were TP53, EGFR, KRAS, FAT4, and KMT2D. Interestingly, we observed the highest incidence of SSM in the Tumor Protein p53 (TP53) gene at 63.64%. Similarly, the SSM incidence was 23.94% in the Epidermal Growth Factor Receptor (EGFR), 22.12% in the Kirsten Rat Sarcoma Viral Oncogene Homolog (KRAS), 18.48% in the FAT Atypical Cadherin 4 (FAT4), and 14.24% in the Lysine (K)-Specific Methyltransferase 2D (KMT2D) genes. Subsequently, we used a bioinformatics approach to assess the effect of miRNA-mRNA interactions in lung cancer among the top five SSM-enriched mRNA proto-oncogenes. Among the top 20 identified and selected miRNAs, we observed 18 unique microRNAs that bind specifically to TP53, KRAS, and FAT4 genes and 17 and 19 microRNAs that exclusively bind with the EGFR and KMT2D genes, respectively. Our study found that the top five SSM-enriched mRNA proto-oncogenes in lung cancers among tobacco smokers were TP53, EGFR, KRAS, FAT4, and KMT2D. Further, our results provide an important insight into the involvement of the intricate network of mRNA-miRNA interactions in the development of lung cancer.

GraphTar: applying word2vec and graph neural networks to miRNA target prediction

BackgroundMicroRNAs (miRNAs) are short, non-coding RNA molecules that regulate gene expression by binding to specific mRNAs, inhibiting their translation. They play a critical role in regulating various biological processes and are implicated in many diseases, including cardiovascular, oncological, gastrointestinal diseases, and viral infections. Computational methods that can identify potential miRNA–mRNA interactions from raw data use one-dimensional miRNA–mRNA duplex representations and simple sequence encoding techniques, which may limit their performance.ResultsWe have developed GraphTar, a new target prediction method that uses a novel graph-based representation to reflect the spatial structure of the miRNA–mRNA duplex. Unlike existing approaches, we use the word2vec method to accurately encode RNA sequence information. In conjunction with the novel encoding method, we use a graph neural network classifier that can accurately predict miRNA–mRNA interactions based on graph representation learning. As part of a comparative study, we evaluate three different node embedding approaches within the GraphTar framework and compare them with other state-of-the-art target prediction methods. The results show that the proposed method achieves similar performance to the best methods in the field and outperforms them on one of the datasets.ConclusionsIn this study, a novel miRNA target prediction approach called GraphTar is introduced. Results show that GraphTar is as effective as existing methods and even outperforms them in some cases, opening new avenues for further research. However, the expansion of available datasets is critical for advancing the field towards real-world applications.

Identification and analysis of differential miRNA-mRNA interactions in coronary heart disease: an experimental screening approach.

This aim of this study is to screen the differential molecules of kidney deficiency and blood stasis (KDBS) syndrome in coronary heart disease by high-throughput sequencing. In addition, the study aims to verify the alterations in the expression levels of miR-4685-3p and its regulated downstream, namely, C1QC, C4, and C5, using quantitative polymerase chain reaction (qPCR) and enzyme-linked immunosorbent assay (ELISA), and to determine whether the complement and coagulation cascade pathway is the specific pathogenic pathway. Patients diagnosed with unstable angina pectoris with KDBS syndrome, patients with non-kidney deficiency blood stasis (NKDBS) syndrome, and a Normal group were recruited. The clinical symptoms of each group were further analyzed. Illumina's NextSeq 2000 sequencing platform and FastQC software were used for RNA sequencing and quality control. DESeq software was used for differential gene expression (DGE) analysis. qPCR and ELISA verification were performed on DGE analysis. The DGE profiles of 77 miRNA and 331 mRNA were selected. The GO enrichment analysis comprised 43 biological processes, 49 cell components, and 42 molecular functions. The KEGG enrichment results included 40 KEGG pathways. The PCR results showed that, compared with the Normal group, the miR-4685-3p levels decreased in the CHD_KDBS group (P = 0.001), and were found to be lower than those observed in the CHD_NKDBS group. The downstream mRNA C1 regulated by miR-4685-3p showed an increasing trend in the CHD_KDBS group, which was higher than that in the Normal group (P = 0.0019). The mRNA C4 and C5 in the CHD_KDBS group showed an upward trend, but the difference was not statistically significant. ELISA was utilized for the detection of proteins associated with the complement and coagulation cascade pathway. It was found that the expression level of C1 was significantly upregulated in the CHD_KDBS group compared with the Normal group (P < 0.0001), which was seen to be higher than that in the CHD_NKDBS group (P < 0.0001). The expression levels of C4 and C5 in the CHD_KDBS group were significantly lower than the Normal group, and were lower than that in the CHD_NKDBS group (P < 0.0001). The occurrence of CHD_KDBS might be related to the activation of the complement and coagulation cascade pathway, which is demonstrated by the observed decrease in miR-4685-3p and the subsequent upregulation of its downstream C1QC. In addition, the expression levels of complement C4 and C5 were found to be decreased, which provided a research basis for the prevention and treatment of this disease.

The miRNA–mRNA regulatory networks of the response to NaHCO3 stress in industrial hemp (Cannabis sativa L.)

BackgroundIndustrial hemp is an important industrial crop and has strong resistance to saline-alkaline stress. However, research on the industrial hemp response to NaHCO3 stress is limited. Therefore, the response mechanisms of industrial hemp under NaHCO3 stress were analysed through miRNA–mRNA regulatory networks.ResultsSeedlings of two salt–alkali tolerant and sensitive varieties were cultured in a solution containing 100 mM NaHCO3 and randomly sampled at 0, 6, 12, and 24 h. With prolonged NaHCO3 stress, the seedlings gradually withered, and the contents of jasmonic acid, lignin, trehalose, soluble protein, peroxidase, and superoxide dismutase in the roots increased significantly. The abscisic acid content decreased and then gradually increased. Overall, 18,215 mRNAs and 74 miRNAs were identified as differentially expressed under NaHCO3 stress. The network showed that 230 miRNA–mRNA interactions involved 16 miRNAs and 179 mRNAs, including some key hub novel mRNAs of these crucial pathways. Carbon metabolism, starch, sucrose metabolism, plant hormone signal transduction, and the spliceosome (SPL) were crucial pathways in industrial hemp's response to NaHCO3 stress.ConclusionsIt is speculated that industrial hemp can regulate SPL pathway by upregulating miRNAs such as novel_miR_179 and novel_miR_75, thus affecting starch and sucrose metabolism, plant hormone signal transduction and carbon metabolism and improving key physiological indices such as jasmonic acid content, trehalose content, and peroxidase and superoxide dismutase activities under NaHCO3 stress.

LDLR and PCSK9 3´UTR variants and their putative effects on microRNA molecular interactions in familial hypercholesterolemia: a computational approach.

Familial hypercholesterolemia (FH) is caused by pathogenic variants in low-density lipoprotein (LDL) receptor (LDLR) or its associated genes, including apolipoprotein B (APOB), proprotein convertase subtilisin/kexin type 9 (PCSK9), and LDLR adaptor protein 1 (LDLRAP1). However, approximately 40% of the FH patients clinically diagnosed (based on FH phenotypes) may not carry a causal variant in a FH-related gene. Variants located at 3' untranslated region (UTR) of FH-related genes could elucidate mechanisms involved in FH pathogenesis. This study used a computational approach to assess the effects of 3'UTR variants in FH-related genes on miRNAs molecular interactions and to explore the association of these variants with molecular diagnosis of FH. Exons and regulatory regions of FH-related genes were sequenced in 83 FH patients using an exon-target gene sequencing strategy. In silico prediction tools were used to study the effects of 3´UTR variants on interactions between miRNAs and target mRNAs. Pathogenic variants in FH-related genes (molecular diagnosis) were detected in 44.6% FH patients. Among 59 3'UTR variants identified, LDLR rs5742911 and PCSK9 rs17111557 were associated with molecular diagnosis of FH, whereas LDLR rs7258146 and rs7254521 and LDLRAP1 rs397860393 had an opposite effect (p < 0.05). 3´UTR variants in LDLR (rs5742911, rs7258146, rs7254521) and PCSK9 (rs17111557) disrupt interactions with several miRNAs, and more stable bindings were found with LDLR (miR-4435, miR-509-3 and miR-502) and PCSK9 (miR-4796). LDLR and PCSK9 3´UTR variants disturb miRNA:mRNA interactions that could affect gene expression and are potentially associated with molecular diagnosis of FH.

Seasonal patterns of miRNA and mRNA expression profiles in the testes of plateau zokors (Eospalax baileyi)

The gonads of seasonal breeding animals undergo periodic annual changes in morphology, physiological hormones, and gene expression levels. To clarify the regulatory mechanism of miRNAs in the seasonal testicular development and spermatogenesis of plateau zokors, the miRNA expression profiles in their testicles during breeding and non-breeding seasons were analyzed. In total, 447 miRNAs, including 366, 81, and 167 known, novel, and differentially expressed (DE) miRNAs, respectively, were determined in the testes. Compared to the non-breeding season, 90 DE miRNAs were upregulated and 77 DE miRNAs were downregulated during the breeding season. By analysing the miRNA and mRNA expression profiles, we predicted 2096 significant target mRNAs. According to the miRNA–mRNA interaction network, target mRNAs with DE miRNAs were related to testicular development and spermatogenesis. GO indicated that target mRNAs were enriched in spermatogenesis, cell differentiation, multicellular biological development, and flagellated sperm movement and were associated with regulating testicular development and spermatogenesis. KEGG suggested that target mRNAs were enriched in lipid and fructose metabolism and provided energy and material for spermatogenesis. The target mRNA of rno-miR-24-3p was determined to be Polyubiquitin-B (UBB). Our results provide a reference for revealing the mechanism by which miRNAs regulate testicular development and spermatogenesis in plateau zokors, which has important implications for understanding the regulation of seasonal reproduction in animals.

HGCLMDA: Predicting mRNA-Drug Sensitivity Associations via Hypergraph Contrastive Learning.

The identification of drug sensitivity to mRNA interactions is crucial for drug development and disease treatment, but traditional experimental methods for verifying mRNA-drug sensitivity associations are labor-intensive and time-consuming. In this study, we present a hypergraph contrastive learning approach, HGCLMDA, to predict potential mRNA-drug sensitivity associations. HGCLMDA integrates a graph convolutional network-based method with a hypergraph convolutional network to mine high-order relationships between mRNA-drug association pairs. The proposed cross-view contrastive learning architecture improves the model's learning ability, and the inner product is used to obtain the mRNA-drug sensitivity association score. Our experiments on three mRNA-drug sensitivity association data sets show that HGCLMDA outperforms traditional graph convolutional network-based methods, graph augmentation-based contrastive learning methods, and state-of-the-art association prediction methods. The visualization experiment demonstrates the strong discrimination ability of the mRNA and drug embeddings learned by HGCLMDA, and experiments on sparse data sets showcase the performance and robustness of the method. In-depth analysis of hypergraph structures reveals a crucial role that hypergraphs play in enhancing the performance of models. The case study highlights the potential of HGCLMDA as a valuable tool for predicting mRNA-drug sensitivity associations. The interpretive analysis reveals that HGCLMDA effectively models the similarity between mRNA-mRNA and drug-drug interactions.

Translation factor and RNA binding protein mRNA interactomes support broader RNA regulons for posttranscriptional control

The regulation of translation provides a rapid and direct mechanism to modulate the cellular proteome. In eukaryotes, an established model for the recruitment of ribosomes to mRNA depends upon a set of conserved translation initiation factors. Nevertheless, how cells orchestrate and define the selection of individual mRNAs for translation, as opposed to other potential cytosolic fates, is poorly understood. We have previously found significant variation in the interaction between individual mRNAs and an array of translation initiation factors. Indeed, mRNAs can be separated into different classes based upon these interactions to provide a framework for understanding different modes of translation initiation. Here, we extend this approach to include new mRNA interaction profiles for additional proteins involved in shaping the cytoplasmic fate of mRNAs. This work defines a set of seven mRNA clusters, based on their interaction profiles with 12 factors involved in translation and/or RNA binding. The mRNA clusters share both physical and functional characteristics to provide a rationale for the interaction profiles. Moreover, a comparison with mRNA interaction profiles from a host of RNA binding proteins suggests that there are defined patterns in the interactions of functionally related mRNAs. Therefore, this work defines global cytoplasmic mRNA binding modules that likely coordinate the synthesis of functionally related proteins.

Integrated analysis of high‑throughput sequencing reveals the regulatory potential of hsa_circ_0035431 in HNSCC.

Circular RNAs (circRNAs) are molecular sponges that are involved in regulation of multiple types of cancer. The present study aimed to screen and explore the key circRNA/microRNA (miRNA or miR)/mRNA interactions in head and neck squamous cell carcinoma (HNSCC) using bioinformatics. A total of six pairs of cancerous and adjacent healthy tissue were obtained from patients with HNSCC and genome-wide transcriptional sequencing was performed. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed on differentially expressed genes (DEGs). Moreover, expression levels of DEGs were verified in HNSCC cells and tissues using reverse transcription-quantitative (RT-q)PCR. A molecular regulatory network consisting of three circRNAs, seven miRNAs and seven mRNAs was constructed, resulting in identification of two signaling axes, hsa_circ_0035431/hsa-miR-940/fucosyltransferase 6 (FUT6) and hsa_circ_0035431/hsa-miR-940/cingulin-like 1 (CGNL1). FUT6 and CGNL1 were downregulated in HNSCC compared with adjacent healthy tissue and the expression levels of these genes were associated with tumor stage. Low FUT6 and CGNL1 expression levels were associated with lower overall survival rate and progression-free intervals in HNSCC. RT-qPCR demonstrated that hsa_circ_0035431, FUT6 and CGNL1 were downregulated in HNSCC cells and tissue and hsa-miR-940 was upregulated. Notably, these results were consistent with those obtained using high-throughput sequencing. In conclusion, hsa_circ_0035431 may participate in regulation of FUT6 and CGNL1 expression by sponging hsa-miR-940, thus, impacting the occurrence, development and prognosis of HNSCC.

CircRNA/lncRNA–miRNA–mRNA network and gene landscape in calcific aortic valve disease

BackgroundCalcific aortic valve disease (CAVD) is a common valve disease with an increasing incidence, but no effective drugs as of yet. With the development of sequencing technology, non-coding RNAs have been found to play roles in many diseases as well as CAVD, but no circRNA/lncRNA–miRNA–mRNA interaction axis has been established. Moreover, valve interstitial cells (VICs) and valvular endothelial cells (VECs) play important roles in CAVD, and CAVD differed between leaflet phenotypes and genders. This work aims to explore the mechanism of circRNA/lncRNA–miRNA–mRNA network in CAVD, and perform subgroup analysis on the important characteristics of CAVD, such as key cells, leaflet phenotypes and genders.ResultsWe identified 158 differentially expressed circRNAs (DEcircRNAs), 397 DElncRNAs, 45 DEmiRNAs and 167 DEmRNAs, and constructed a hsa-circ-0073813/hsa-circ-0027587–hsa-miR-525-5p–SPP1/HMOX1/CD28 network in CAVD after qRT-PCR verification. Additionally, 17 differentially expressed genes (DEGs) in VICs, 9 DEGs in VECs, 7 DEGs between different leaflet phenotypes and 24 DEGs between different genders were identified. Enrichment analysis suggested the potentially important pathways in inflammation and fibro-calcification during the pathogenesis of CAVD, and immune cell patterns in CAVD suggest that M0 macrophages and memory B cells memory were significantly increased, and many genes in immune cells were also differently expressed.ConclusionsThe circRNA/lncRNA–miRNA–mRNA interaction axis constructed in this work and the DEGs identified between different characteristics of CAVD provide a direction for a deeper understanding of CAVD and provide possible diagnostic markers and treatment targets for CAVD in the future.

Administration of mRNA-Nanomedicine-Augmented Calvarial Defect Healing via Endochondral Ossification.

Large-area craniofacial defects remain a challenge for orthopaedists, hastening the need to develop a facile and safe tissue engineering strategy; osteoconductive material and a combination of optimal growth factors and microenvironment should be considered. Faced with the unmet need, we propose that abundant cytokines and chemokines can be secreted from the bone defect, provoking the infiltration of endogenous stem cells to assist bone regeneration. We can provide a potent mRNA medicine cocktail to promptly initiate the formation of bone templates, osteogenesis, and subsequent bone matrix deposition via endochondral ossification, which may retard rapid fibroblast infiltration and prevent the formation of atrophic non-union. We explored the mutual interaction of BMP2 and TGFβ3 mRNA, both potent chondrogenic factors, on inducing endochondral ossification; examined the influence of in vitro the transcribed polyA tail length on mRNA stability; prepared mRNA nanomedicine using a PEGylated polyaspartamide block copolymer loaded in a gelatin sponge and grafted in a critical-sized calvarial defect; and evaluated bone regeneration using histological and μCT examination. The BMP2 and TGFβ3 composite mRNA nanomedicine resulted in over 10-fold new bone volume (BV) regeneration in 8 weeks than the BMP2 mRNA nanomedicine administration alone, demonstrating that the TGFβ3 mRNA nanomedicine synergistically enhances the bone's formation capability, which is induced by BMP2 mRNA nanomedicine. Our data demonstrated that mRNA-medicine-mediated endochondral ossification provides an alternative cell-free tissue engineering methodology for guiding craniofacial defect healing.