miRNA-mRNA Regulatory Network Research Articles

Screening of miRNAs as prognostic biomarkers and their associated hub targets across Hepatocellular carcinoma using survival-based bioinformatics approach

BackgroundThe hepatocellular carcinoma (HCC) incident rate is gradually increasing yearly despite all the research and efforts taken by scientific communities and governing bodies. Approximately 90% of all liver cancer cases belong to HCC. Usually, HCC patients approach the treatment in the late stages of this malignancy which becomes the primary cause of high mortality rate. The knowledge about molecular pathogenesis of HCC is limited and needs more attention from researchers to identify the driver genes and miRNAs, which causes to translate this information into clinical practice. Therefore, the key regulators identification of miRNA-mRNA regulatory network is essential to identify HCC-associated genes. MethodologyWe extracted microRNA (miRNA) and messenger RNA (mRNA) expression datasets of normal and tumor HCC patient samples from UCSC Xena followed by identifying differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs). Univariate and multivariate cox-proportional hazard models were utilized to identify DEMs having significant association with overall survival (OS). Kaplan-Meier (KM) plotter was used to validate the presence of prognostic DEMs. A risk-score model was used to evaluate the effectiveness of KM-plotter validated DEMs combination on risk of samples. Target DEGs of prognostic miRNAs were identified via sources such as miRTargetLink and miRWalk followed by their validation in an external microarray cohort and enrichment analysis. Results562 DEGs and 388 DEMs were identified followed by seven prognostic miRNAs (i.e., miR-19a, miR-19b, miR-30d-5p, miR-424-5p, miR-3677-5p, miR-3913-5p, miR-7705) post univariate, multivariate, risk-score model evaluation and KM-plotter analyses. ANLN, MRO, CPEB3 were their targets and were also validated in GSE84005 dataset. ConclusionsThe findings of this study decipher that most significant miRNAs and their identified target genes have association with apoptosis, inflammation, cell cycle regulation and cancer-related pathways, which appear to contribute to HCC pathogenesis and therefore, the discovery of new targets.

Lnc-PSMA8-1 activated by GEFT promotes rhabdomyosarcoma progression via upregulation of mTOR expression by sponging miR-144-3p

BackgroundGEFT is a key regulator of tumorigenesis in rhabdomyosarcoma (RMS), and overexpression of GEFT is significantly correlated with distant metastasis, lymph node metastasis, and a poor prognosis, yet the underlying molecular mechanism is still poorly understood. This study aimed to investigate and validate the molecular mechanism of GEFT-activated lncRNAs in regulating mTOR expression to promote the progression of RMS.MethodsGEFT-regulated lncRNAs were identified through microarray analysis. The effects of GEFT-regulated lncRNAs on the proliferation, apoptosis, invasion, and migration of RMS cells were confirmed through cell functional experiments. The target miRNAs of GEFT-activated lncRNAs in the regulation of mTOR expression were predicted by bioinformatics analysis combined with quantitative real-time polymerase chain reaction (qRT–PCR) analysis. The expression of lnc-PSMA8-1, miR-144-3p, and mTOR was measured by qRT–PCR in RMS tissue samples and cell lines. The regulatory mechanisms of the lnc-PSMA8-1-miR-144-3p-mTOR signaling axis were verified by RNA-binding protein immunoprecipitation (RIP), a luciferase reporter assay, qRT–PCR analysis, Western blot analysis, and cell functional experiments.ResultsThe microarray-based analysis identified 31 differentially expressed lncRNAs (fold change > 2.0, P < 0.05). Silencing the 4 upregulated lncRNAs (lnc-CEACAM19-1, lnc-VWCE-2, lnc-GPX7-1, and lnc-PSMA8-1) and overexpressing the downregulated lnc-FAM59A-1 inhibited the proliferation, invasion, and migration and induced the apoptosis of RMS cells. Among the factors analyzed, the expression of lnc-PSMA8-1, miR-144-3p, and mTOR in RMS tissue samples and cells was consistent with the correlations among their expression indicated by the lncRNA–miRNA–mRNA regulatory network based on the ceRNA hypothesis. lnc-PSMA8-1 promoted RMS progression by competitively binding to miR-144-3p to regulate mTOR expression.ConclusionOur research demonstrated that lnc-PSMA8-1 was activated by GEFT and that the former positively regulated mTOR expression by sponging miR-144-3p to promote the progression of RMS. Therefore, targeting this network may constitute a potential therapeutic approach for the management of RMS.

Unraveling the regulatory network of miRNA expression in Potato Y virus-infected of Nicotiana benthamiana using integrated small RNA and transcriptome sequencing.

Potato virus Y (PVY) disease is a global problem that causes significant damage to crop quality and yield. As traditional chemical control methods are ineffective against PVY, it is crucial to explore new control strategies. MicroRNAs (miRNAs) play a crucial role in plant and animal defense responses to biotic and abiotic stresses. These endogenous miRNAs act as a link between antiviral gene pathways and host immunity. Several miRNAs target plant immune genes and are involved in the virus infection process. In this study, we conducted small RNA sequencing and transcriptome sequencing on healthy and PVY-infected N. benthamiana tissues (roots, stems, and leaves). Through bioinformatics analysis, we predicted potential targets of differentially expressed miRNAs using the N. benthamiana reference genome and the PVY genome. We then compared the identified differentially expressed mRNAs with the predicted target genes to uncover the complex relationships between miRNAs and their targets. This study successfully constructed a miRNA-mRNA network through the joint analysis of Small RNA sequencing and transcriptome sequencing, which unveiled potential miRNA targets and identified potential binding sites of miRNAs on the PVY genome. This miRNA-mRNA regulatory network suggests the involvement of miRNAs in the virus infection process.

Unravelling microRNA regulation and miRNA-mRNA regulatory networks in osteogenesis driven by 3D nanotopographical cues.

Three-dimensional (3D) culturing of cells is being adopted for developing tissues for various applications such as mechanistic studies, drug testing, tissue regeneration, and animal-free meat. These approaches often involve cost-effective differentiation of stem or progenitor cells. One approach is to exploit architectural cues on a 3D substrate to drive cellular differentiation, which has been shown to be effective in various studies. Although extensive gene expression data from such studies have shown that gene expression patterns might differ, the gene regulatory networks controlling the expression of genes are rarely studied. In this study, we profiled genes and microRNAs (miRNAs) via next-generation sequencing (NGS) in human mesenchymal stem cells (hMSCs) driven toward osteogenesis via architectural cues in 3D matrices (3D conditions) and compared with cells in two-dimensional (2D) culture driven toward osteogenesis via soluble osteoinductive factors (OF conditions). The total number of differentially expressed genes was smaller in 3D compared to OF conditions. A distinct set of genes was observed under these conditions that have been shown to control osteogenic differentiation via different pathways. Small RNA sequencing revealed a core set of miRNAs to be differentially expressed under these conditions, similar to those that have been previously implicated in osteogenesis. We also observed a distinct regulation of miRNAs in these samples that can modulate gene expression, suggesting supplementary gene regulatory networks operative under different stimuli. This study provides insights into studying gene regulatory networks for identifying critical nodes to target for enhanced cellular differentiation and reveal the differences in physical and biochemical cues to drive cell fates.

What happens to the osteoporotic bone mesenchymal stem cells? Evidence from RNA sequencing.

Evidence presented that osteoporosis is closely related to the dysfunction of bone mesenchymal stem cells (BMSCs). But most studies are insufficient to reveal what actually happens to the osteoporotic BMSCs. In this study, BMSCs were harvested from ovariectomized and sham-operated rats. After checking the characteristics of rat models and stem cells, the BMSCs were carried out for RNA sequencing. Part of the findings were verified that seven mRNAs (Abi3bp, Aifm3, Ccl11, Cdkn1c, Chst10, Id2, Vcam1) were significantly up-regulated in osteoporotic BMSCs while seven mRNAs (Cep63, Fgfr3, Myc, Omd, Pou2f1, Smarcal1, Timm10b) were down-regulated. In addition, potential miRNA-mRNA and lncRNA-mRNA regulatory networks were illustrated. The changes in osteoporotic BMSCs covered a large set of biological processes, including cell viability, differentiation, immunoreaction, bone repairment and estrogen defect. This study enriched the pathophysiological mechanisms of BMSCs and osteporosis, as well as provided dozens of attractive RNA targets for further treatment.

The landscape of miRNA-mRNA regulatory network and cellular sources in inflammatory bowel diseases: insights from text mining and single cell RNA sequencing analysis.

Inflammatory Bowel Diseases (IBDs), encompassing Ulcerative Colitis (UC) and Crohn's Disease (CD), are chronic, recurrent inflammatory conditions of the gastrointestinal tract. The microRNA (miRNA) -mRNA regulatory network is pivotal in the initiation and progression of IBDs. Although individual studies provide valuable insights into miRNA mechanisms in IBDs, they often have limited scope due to constraints in population diversity, sample size, sequencing platform variability, batch effects, and potential researcher bias. Our study aimed to construct comprehensive miRNA-mRNA regulatory networks and determine the cellular sources and functions of key miRNAs in IBD pathogenesis. To minimize potential bias from individual studies, we utilized a text mining-based approach on published scientific literature from PubMed and PMC databases to identify miRNAs and mRNAs associated with IBDs and their subtypes. We constructed miRNA-mRNA regulatory networks by integrating both predicted and experimentally validated results from DIANA, Targetscan, PicTar, Miranda, miRDB, and miRTarBase (all of which are databases for miRNA target annotation). The functions of miRNAs were determined through gene enrichment analysis of their target mRNAs. Additionally, we used two large-scale single-cell RNA sequencing datasets to identify the cellular sources of miRNAs and the association of their expression levels with clinical status, molecular and functional alternation in CD and UC. Our analysis systematically summarized IBD-related genes using text-mining methodologies. We constructed three comprehensive miRNA-mRNA regulatory networks specific to IBD, CD, and UC. Through cross-analysis with two large-scale scRNA-seq datasets, we determined the cellular sources of the identified miRNAs. Despite originating from different cell types, hsa-miR-142, hsa-miR-145, and hsa-miR-146a were common to both CD and UC. Notably, hsa-miR-145 was identified as myofibroblast-specific in both CD and UC. Furthermore, we found that higher tissue repair and enhanced glucose and lipid metabolism were associated with hsa-miR-145 in myofibroblasts in both CD and UC contexts. This comprehensive approach revealed common and distinct miRNA-mRNA regulatory networks in CD and UC, identified cell-specific miRNA expressions (notably hsa-miR-145 in myofibroblasts), and linked miRNA expression to functional alterations in IBD. These findings not only enhance our understanding of IBD pathogenesis but also offer promising diagnostic biomarkers and therapeutic targets for clinical practice in managing IBDs.

Exosomal hsa-miR-151a-3p and hsa-miR-877-5p are potential novel biomarkers for predicting bone metastasis in lung cancer.

Exosomal miRNAs (exo-miRNAs) have arisen as novel diagnostic biomarkers for various cancers. However, few reports on exo-miRNAs related to bone metastasis (BM) in lung cancer exist. This study aims to screen out key exo-miRNAs and estimate their prognostic values for predicting BM in lung cancer. The differentially expressed exo-miRNAs between the highly-metastatic (95D) and lowly-metastatic (A549) human lung cancer cell lines were comprehensively analyzed using high-throughput sequencing followed by bioinformatic analyses. 29 candidate exo-miRNAs were identified, and 101 BM-related target genes were predicted. Enrichment analysis revealed that these target genes were mainly involved in regulating transcription and pathways in cancer. An exosomal miRNA-mRNA regulatory network consisting of 7 key miRNAs and 10 hub genes was constructed. Further function analysis indicated that these 10 hub genes were mainly enriched in regulating cancer's apoptosis and central carbon metabolism. The survival analysis indicated that 7 of 10 hub genes were closely related to prognosis. Mutation analysis showed that lung cancer patients presented certain genetic alterations in the 7 real hub genes. GSEA for a single hub gene suggested that 6 of 7 real hub genes had close associations with lung cancer development. Finally, ROC analysis revealed that hsa-miR-151a-3p and hsa-miR-877-5p provided high diagnostic accuracy in discriminating patients with bone metastasis (BM+) from patients without bone metastasis (BM-). These findings provided a comprehensive analysis of exo-miRNAs and target genes in the regulatory network of BM in lung cancer. In particular, hsa-miR-151a-3p and hsa-miR-877-5p may be novel biomarkers for predicting BM in lung cancer.

Identification of Potential miRNA-mRNA Regulatory Network Associated with Growth and Development of Hair Follicles in Forest Musk Deer.

In this study, sRNA libraries and mRNA libraries of HFs of FMD were constructed and sequenced using an Illumina HiSeq 2500, and the expression profiles of miRNAs and genes in the HFs of FMD were obtained at the anagen and catagen stages. In total, 565 differentially expressed unigenes (DEGs) were identified, 90 of which were upregulated and 475 of which were downregulated. In the BP category of GO enrichment, the DEGs were enriched in the processes related to HF development and differentiation, including the hair cycle regulation and processes, HF development, skin epidermis development, regulation of HF development, skin development, the Wnt signaling pathway, and the BMP signaling pathway. Through KEGG analysis it was found that DEGs were significantly enriched in pathways associated with HF development and growth. A total of 186 differentially expressed miRNAs (DEmiRNAs) were screened (p < 0.05) in the HFs of FMD at the anagen stage vs. the catagen stage, 33 of which were upregulated and 153 of which were downregulated. Through DEmiRNA-mRNA association analysis, we found DEmiRNAs and target genes that mainly play regulatory roles in HF development and growth. The enrichment analysis of DEmiRNA target genes revealed similarities with the enrichment results of DEGs associated with HF development. Notably, both sets of genes were enriched in key pathways such as the Notch signaling pathway, melanogenesis, the cAMP signaling pathway, and cGMP-PKG. To validate our findings, we selected 11 DEGs and 11 DEmiRNAs for experimental verification using RT-qPCR. The results of the experimental validation were consistent with the RNA-Seq results.

Identifying genes for regulating osteogenic differentiation of human periodontal ligament stem cells in inflammatory environments by bioinformatics analysis.

Periodontitis is an immuno-inflammatory disease caused by dental plaque biofilms and inflammations. The regeneration of bone tissue in inflammatory environment is of great significance for the treatment of periodontal disease, but the specific molecular mechanism of bone formation in periodontitis still needs further exploration. The objective of this study was to identify key osteogenesis-related genes (ORGs) in periodontitis. We used two datasets from the Gene Expression Omnibus (GEO) database to find differentially expressed mRNAs and miRNAs, further performed Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Then we predicted the downstream genes of the differentially expressed miRNAs (DEMs) by the TargetScan database and established a miRNA-mRNA regulatory network. Finally, the osteogenic mechanism of periodontitis was explored through quantitative real-time PCR (qRT-PCR) by inducing inflammatory environment and osteogenic differentiation of hPDLSCs. Through differential expression analysis and prediction of downstream target genes of DEMs, we created a miRNA-mRNA regulatory network consisting of 29 DEMs and 11 differentially expressed osteogenesis-related genes (DEORGs). In addition, the qRT-PCR results demonstrated that BTBD3, PLAT, AKAP12, SGK1, and GLCE expression levels were significantly upregulated, while those of TIMP3, ZCCHC14, LIN7A, DNAH6, NNT, and ITGA6 were downregulated under the dual effects of inflammatory stimulation and osteogenic induction. DEORGs might be important factors in the osteogenic phase of periodontitis, and the miRNA-mRNA network may shed light on the clarification of the role and mechanism of osteogenesis in periodontitis and contribute to the development of novel therapeutic strategies.

Characterization of MicroRNA expression profiles in the ovarian tissue of goats during the sexual maturity period

BackgroundThe ovary is an important reproductive organ in mammals, and its development directly affects the sexual maturity and reproductive capacity of individuals. MicroRNAs (miRNAs) are recognized as regulators of reproductive physiological processes in various animals and have been shown to regulate ovarian development through typical targeting and translational repression. However, little is known about the regulatory role of miRNAs in ovarian tissue development during sexual maturity in goats. To comprehensively profile the different physiological stages of sexual maturation in goats, we performed small-RNA sequencing of ovarian tissue samples collected at four specific time points (1 day after birth (D1), 2 months old (M2), 4 months old (M4), and 6 months old (M6)). In addition, we used ELISAs to measure serum levels of reproductive hormones to study their temporal changes.ResultsThe results showed that serum levels of gonadotropin-releasing hormone, follicle-stimulating hormone, luteinizing hormone, oestradiol, progesterone, oxytocin, and prolactin were lower in goats at the D1 stage than at the other three developmental stages (P < 0.05). The secretion patterns of these seven hormones show a similar trend, with hormone levels reaching their peaks at 4 months of age. A total of 667 miRNAs were detected in 20 libraries, and 254 differentially expressed miRNAs and 3 groups of miRNA clusters that had unique expression patterns were identified (|log2-fold change|> 1, FDR < 0.05) in the 6 comparison groups. RT‒qPCR was employed to confirm that the expression pattern of the 15 selected miRNAs was consistent with the Illumina sequencing results. Gene ontology analyses revealed significant enrichment of GO terms such as cell proliferation regulation, epithelial cell development, and amino acid transport, as well as important signaling pathways including the MAPK signaling pathway, the PI3K-Akt signaling pathway, and the oestrogen signaling pathway. Further miRNA‒mRNA regulation network analysis revealed that 8 differentially expressed miRNAs (chi-miR-1343, chi-miR-328-3p, chi-miR-877-3p, chi-miR-296-3p, chi-miR-128-5p, chi-miR-331-3p, chi-miR-342-5p and chi-miR-34a) have important regulatory roles in ovarian cell proliferation, hormone secretion and metabolism-related biological processes.ConclusionsOverall, our study investigated the changes in serum hormone and miRNA levels in the ovaries. These data provide a valuable resource for understanding the molecular regulatory mechanisms of miRNAs in ovarian tissue during the sexual maturity period in goats.

Prognostic Value and the Immune Microenvironment-associated Role of SMAD4 in Pancreatic Adenocarcinoma.

SMAD4 is a well-known cancer suppressor gene that regulates cell proliferation, differentiation, autoimmunity, pluripotency, and immune responses in pancreatic adenocarcinoma (PAAD). We herein investigated a novel involvement of SMAD4 within the PAAD microenvironment. Transcriptome data, derived from The Cancer Genome Atlas and Genotype-Tissue Expression Using TIMER 2.0, CIBERSORT, and ImmuCellAI, were used to identify the immune cell infiltration pattern of PAAD. We then knocked-down SMAD4 in the PANC-1 cell line and acquired RNA-seq data through the Illumina microarray technology. Kyoto Encyclopedia of Genes and Genomes was utilized along with Gene Ontology enrichment analyses, and protein-protein interaction network analysis to screen for genes that were differentially expressed. We constructed a miRNA-mRNA regulatory network and analyzed SMAD4 copy number variation (CNV) data. In PAAD, decreased levels of SMAD4 expression were found to be connected to an unfavorable prognosis. There was a significantly higher infiltration level of DC cells, CD8+ T cells, TgD, Tc, and Tex cells and a lower level of B cells and Th2 in the SMAD4-high group. The expression of SMAD4 and immune cell infiltration including CD8+ T cells, myeloid dendritic cells, neutrophils, and macrophages are significantly positively correlated. DEGs were found enriched in the hypoxia response pathway. The six hypoxia-related genes exhibited a significant correlation with immune cell infiltration and survival rates. SMAD4 CNV levels were associated with MSI, stemness, infiltration of immune cells, and survival rates. A statistically significant correlation was found between SMAD4 expression and immune cell infiltration. SMAD4 could mediate hypoxia response in pancreatic cancer. The CNV levels of SMAD4 were associated with prognosis. SMAD4 has potential as a prognostic biomarker and provides a new orientation for the immunotherapy of PAAD.

Screening and Bioinformatics Analysis of Differential Genes in Autism Spectrum Disorder Based on GEO Database.

The prevalence of autism spectrum disorder (ASD) in children has been increasing year by year, which has seriously affected the quality of life of children. There are many theories about the cause of ASDs, with some studies suggesting that it may be related to gene expression levels or inflammation and immune system dysfunction. But the exact mechanism is not fully understood. profile of gene expression The protein interaction network (PPI) of differentially expressed genes was created using the STRING web tool and GSE77103, which was chosen from the gene expression omnibus (GEO) database. Using the CytoHubba plugin of Cytoscape program, the hub genes were examined. The hub gene regulatory network for miRNA-mRNA was then built. We identified 551 differentially expressed genes(DEGs) in 8 children with ASD and normal children. In addition, we screened out 10 hub genes (MX1, ISG15, IRF7, DDX58, IFIT1, BCL2L1, HPGDS, CTSD, PTGS2 and CD68) that were most associated with the development of ASDs. Then, microRNAs (miRNAs) closely related to hub genes (such as has-miR-27a-5p) were screened, and the miRNA-mRNA regulatory network was constructed. In this study, a total of 10 hub genes were identified, including MX1, ISG15, IRF7, DDX58, IFIT1, BCL2L1, HPGDS, CTSD, PTGS2 and CD68, which are closely related to ASD. These genes may play a key role in the occurrence and progression of ASD. In addition, we also revealed some miRNAs that regulate the hub genes of ASD. These results may deepen our understanding of ASD and provide potential biomarkers and targets for future treatment of patients with ASD.

Identification and validation of m6A RNA regulatory network in pulpitis

BackgroundN6-methyladenosine (m6A) RNA modification regulators play an important role in many human diseases, and its abnormal expression can lead to the occurrence and development of diseases. However, their significance in pulpitis remains largely unknown. Here, we sought to identify and validate the m6A RNA regulatory network in pulpitis.MethodsGene expression data for m6A regulators in human pulpitis and normal pulp tissues from public GEO databases were analyzed. Bioinformatics analysis including Gene ontology (GO) functional, and Kyoto encyclopedia of genes and genomes (KEGG) pathway analyses were performed by R package, and Cytoscape software was used to study the role of m6A miRNA-mRNA regulatory network in pulpitis. Quantitative real-time PCR (qRT-PCR) was performed to validate the expression of key m6A regulators in collected human pulpitis specimens.ResultsDifferential genes between pulpitis and normal groups were found from the GEO database, and further analysis found that there were significant differences in the m6A modification-related genes ALKBH5, METTL14, METTL3, METTL16, RBM15B and YTHDF1. And their interaction relationships and hub genes were determined. The hub m6A regulator targets were enriched in immune cells differentiation, glutamatergic synapse, ephrin receptor binding and osteoclast differentiation in pulpitis. Validation by qRT-PCR showed that the expression of methylases METTL14 and METTL3 was decreased, thus these two genes may play a key role in pulpitis.ConclusionOur study identified and validated the m6A RNA regulatory network in pulpitis. These findings will provide valuable resource to guide the mechanistic and therapeutic analysis of the role of key m6A modulators in pulpitis.

Construction of miRNA-mRNA regulatory network indicates potential biomarkers for primary open-angle glaucoma

BackgroundTrabecular meshwork (TM) dysfunction-induced elevation of intraocular pressure has been identified as the main risk factor of irreversible optic nerve injury in Primary open‑angle glaucoma (POAG). Increasing evidences suggest that microRNA (miRNA) plays a vital role in the pathogenesis of POAG. This study aims to construct a miRNA-mRNA regulatory network and identify biomarkers for POAG.MethodsmiRNAs and mRNAs expression profiling of TM samples from controls and POAG patients were assessed through microarray analysis. Target genes of differentially expressed miRNAs (DEmiRNAs) were predicted by miEAA and miRNet. Then GO and KEGG pathway analysis of differentially expressed mRNAs (DEmRNAs) were performed. PPI of top 30 hub genes was identified and miRNA-mRNA network was established by STRING database and Cytoscape software. GSE27276 and GSE105269 datasets were used to verify the expression of hub genes and to predict potential biomarkers in TM and aqueous humor (AH) for POAG, respectively. Finally, GSEA analysis was conducted to estimate the main signaling pathway of POAG pathogenesis.ResultsA total of 29 up-regulated and 7 down-regulated miRNAs, 923 up-regulated and 887 down-regulated mRNAs were identified in TM of POAG compared with controls. Target genes and DEmRNAs were mainly enriched in nitric oxide biosynthetic process, vasopressin-regulated water reabsorption, and so on. Through miRNA-mRNA network construction, top 30 hub genes were regulated by 24 DEmiRNAs. 8 genes were aberrantly expressed in dataset GSE27276. 3 genes (CREB1, CAPZA2, SLC2A3) and 2 miRNAs (miR-106b-5p, miR-15a-5p) were identified as potential biomarkers for POAG in TM and AH, respectively. GSEA analysis revealed that these 3 genes modulated POAG through different pathways.ConclusionIn this study, construction of miRNA-mRNA network and identification of biomarkers provide a novel insight into the pathogenesis, early diagnosis and treatment for POAG.

Identification of miRNA-mRNA Regulatory Networks Associated with Diabetic Retinopathy using Bioinformatics Analysis.

Diabetic retinopathy (DR) is a major complication of diabetes and a leading cause of visual loss. This study aimed to explore biomarkers for DR that may provide additional reference to DR pathogenesis and development. The differentially expressed genes (DEGs) between the DR and control samples in the GSE53257 dataset were identified. Logistics analyses were performed to identify DR-associated miRNAs and genes, and correlation analysis was performed to determine the correlation between them in GSE160306. A total of 114 DEGs in DR were identified in GSE53257. Three genes, including ATP5A1 (down), DAUFV2 (down), and OXA1L (down), were differentially expressed between DR and control samples in GSE160306. Univariate logistics analysis identified that ATP5A1 (OR=0.007, p = 1.40E-02), NDUFV2 (OR = 0.003, p = 6.40E-03), and OXA1L (OR = 0.093, p = 3.08E-02) were DR-associated genes. ATP5A1 and OXA1L were regulated by multiple miRNAs, of which hsa-let- 7b-5p (OR = 26.071, p = 4.40E-03) and hsa-miR-31-5p (OR = 4.188, p = 5.09E-02) were related to DR. ATP5A1 and OXA1L were closely correlated with each other in DR. The hsa-miR-31-5p-ATP5A1 and hsa-let-7b-5p-OXA1L axes might play novel and important roles in the pathogenesis and development of DR.

Identifying Fibroblast Growth Factor Receptor 3 as a Mediator of Periosteal Osteochondral Differentiation through the Construction of microRNA-Based Interaction Networks.

Human periosteum-derived progenitor cells (hPDCs) have the ability to differentiate towards both the chondrogenic and osteogenic lineages. This coordinated and complex osteochondrogenic differentiation process permits endochondral ossification and is essential in bone development and repair. We have previously shown that humanised cultures of hPDCs enhance their osteochondrogenic potentials in vitro and in vivo; however, the underlying mechanisms are largely unknown. This study aimed to identify novel regulators of hPDC osteochondrogenic differentiation through the construction of miRNA-mRNA regulatory networks derived from hPDCs cultured in human serum or foetal bovine serum as an alternative in silico strategy to serum characterisation. Sixteen differentially expressed miRNAs (DEMis) were identified in the humanised culture. In silico analysis of the DEMis with TargetScan allowed for the identification of 1503 potential miRNA target genes. Upon comparison with a paired RNAseq dataset, a 4.5% overlap was observed (122 genes). A protein-protein interaction network created with STRING interestingly identified FGFR3 as a key network node, which was further predicted using multiple pathway analyses. Functional analysis revealed that hPDCs with the activating mutation FGFR3N540K displayed increased expressions of chondrogenic gene markers when cultured under chondrogenic conditions in vitro and displayed enhanced endochondral bone formation in vivo. A further histological analysis uncovered known downstream mediators involved in FGFR3 signalling and endochondral ossification to be upregulated in hPDC FGFR3N540K-seeded implants. This combinational approach of miRNA-mRNA-protein network analysis with in vitro and in vivo characterisation has permitted the identification of FGFR3 as a novel mediator of hPDC biology. Furthermore, this miRNA-based workflow may also allow for the identification of drug targets, which may be of relevance in instances of delayed fracture repair.

Analysis of fatty acid composition and miRNA profile in adipose tissue after application of Allium mongolicum Regel extract in diet of sheep

Extensive regulation of microRNAs (miRNAs) that belong to a type of endogenous small noncoding RNAs on gene expression in sheep has been recognised. Meanwhile, utilisation of water extract of Allium mongolicum Regel (WEA) in diet of sheep has been proved to regulate fatty acid composition and flavour in mutton. However, the regulatory mechanism remains unclear. We conducted miRNA and mRNA sequencing in this study. Briefly, 30 male small-tailed Han sheep (4.5 months old) were randomly assigned to two groups. A basal diet with 3.4 g·sheep−1·d−1 WEA supplementation (WEA group) or with no WEA supplementation (CON group) was fed to sheep, respectively. Subcutaneous and visceral adipose tissue samples were collected during slaughter after 75-day feeding experiment. We obtained 30 differentially expressed miRNAs (DE-miRNAs) after application of WEA in diet of sheep in subcutaneous adipose tissue. Likewise, 15 DE-miRNAs in visceral adipose tissue were identified. A total of 746 differentially expressed mRNAs (DE-mRNAs) with 90 lipid-related DE-mRNAs in subcutaneous adipose tissue and 569 DE-mRNAs with 73 lipid-related DE-mRNAs in visceral adipose tissue were identified after WEA supplementation. The intersection genes from predicted target genes of DE-miRNAs by databases and actual DE-mRNAs by mRNA sequencing were enriched in several lipid-related pathways both in subcutaneous adipose and in visceral adipose tissue. We constructed differential miRNAs-mRNAs networks by Cytoscape software. Our results will help to explain regulatory mechanism of fatty acid composition in mutton by dietary WEA and will provide a theoretical basis for rational use of the plant extract in diet of sheep. HIGHLIGHTS Utilisation of Allium mongolicum Regel water extract led to alterations of endogenous miRNA and gene expression profiles in subcutaneous and visceral adipose tissue of sheep. Allium mongolicum Regel water extract improved nutritional value and flavour of mutton by changing fatty acid composition in adipose tissue. Establishment of miRNA-mRNA regulatory network with regard to manipulation of fatty acid composition in mutton by Allium mongolicum Regel water extract, and provided a theoretical basis for rational utilisation of this plant extract in diet of sheep.

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.

Integrative bioinformatics analysis of miRNA and mRNA expression profiles identified some potential biomarkers for breast cancer

BackgroundBreast cancer is a common cause of cancer death among women with a complex and heterogeneous picture in histological, molecular and clinical features. The aim of this study was to identify hub gene and their target microRNAs in related pathways for breast cancer.MethodsWe selected screening methods for differentially expressed mRNAs and miRNAs using expression profile data of breast cancer from the cancer genome atlas. Using some databases for annotation, the functional and pathway enrichment for differential expression genes was performed. We selected genes and miRNAs with differential expression pattern. Then we determined target genes for differential expression miRNAs (DEMIs) and intersection between them was selected as differentially expressed miRNA–target genes for breast cancer. In the next step, we constructed miRNA–mRNA regulatory network and protein–protein interaction (PPI) network for more information.ResultsTop 10 DEMIs were identified from miRNA profile. Then, we selected 354 genes as target gene for 10 DEMIs. The miRNA–mRNA and PPI network were constructed, and 10 hub genes and 5 miRNAs identified that some of them are new for breast cancer. Also, miRNA–target genes with differential expressions in this study were all mainly involved in signaling pathways and developmental process.ConclusionThis study identified some candidate biomarkers for breast cancer that they have a potential role in pathways related to breast. These findings can be used for research, early diagnosis and therapeutic goals.

Epididymal segment-specific miRNA and mRNA regulatory network at the single cell level

ABSTRACTSpermatozoa released from the testis cannot fertilize an egg before becoming mature and motile in the epididymis. Based on three bulk and one single-cell RNA-seq (scRNA-seq) data series, we compared mRNA or miRNA expression between epididymal segment-specific samples and the other samples. Hereby, we identified 570 differentially expressed mRNAs (DE-mRNAs) and 23 differentially expressed miRNAs (DE-miRNAs) in the caput, 175 DE-mRNAs and 15 DE-miRNAs in the corpus, 946 DE-mRNAs and 12 DE-miRNAs in the cauda. In accordance with respective DE-miRNAs, we predicted upstream transcription factors (TFs) and downstream target genes. Subsequently, we intersected target genes of respective DE-miRNAs with corresponding DE-mRNAs, thereby obtaining 127 upregulated genes in the caput and 92 upregulated genes in cauda. Enriched upregulated pathways included cell motility-related pathways for the caput, smooth muscle-related pathways for the corpus, and immune-associated pathways for the cauda. Protein–protein interaction (PPI) network was constructed to extract key module for the caput and cauda, followed by identifying hub genes through cytohubba. Epididymis tissues from six mice were applied to validate hub genes expression using qRT-PCR, and 7 of the 10 genes displayed identical expression trends in mice caput/cauda. These hub genes were found to be predominantly distributed in spermatozoa using scRNA-seq data. In addition, target genes of DE-miRNAs were intersected with genes in the PPI network for each segment. Subsequently, the miRNA and mRNA regulatory networks for the caput and cauda were constructed. Conclusively, we uncover segment-specific miRNA-mRNA regulatory network, upstream TFs, and downstream pathways of the human epididymis, warranting further investigation into epididymal segment-specific functions.