Targets Of Differentially Expressed miRNAs Research Articles

Identification of serum prognostic marker miRNAs and construction of microRNA-mRNA networks of esophageal cancer.

Esophageal cancer is a common tumor of the digestive system with poor prognosis. This study was to gain a better understanding of the mechanisms involved in esophageal cancer and to identify new prognostic markers. We downloaded the esophageal cancer miRNA expression profile microarray data (GSE113740, GSE112264, GSE122497, GSE113486, and GSE106817) from the GEO database, extracted the esophageal cancer miRNA sequencing data from The Cancer Genome Atlas (TCGA) database, and then used a bioinformatics approach to select common differentially expressed miRNAs (DEMs). Differentially expressed genes (DEGs) were selected by predicting DEM target genes using the miRWalk database and intersecting with differential genes obtained from TCGA database for esophageal cancer. The STRING database was used to obtain protein–protein interaction (PPI) relationships to construct the DEM-DEG network. Furthermore, we selected core genes and core miRNAs associated with esophageal cancer prognosis by performing survival and univariate/multivariate COX analysis on DEMs and DEGs in the network and performed GSEA analysis on core genes alone, and finally the expression of the markers was verified by qPCR in esophageal cancer cell lines Eca109, SKGT-4 and normal esophageal epithelial cells HEEC. Nine DEMs were obtained, of which three were upregulated and six were downregulated, and 326 DEGs were obtained, of which 105 were upregulated and 221 were downregulated. Survival univariate/multivariate COX analysis revealed that five genes, ZBTB16, AQP4, ADCYAP1R1, PDGFD, and VIPR2, and two microRNAs, miR-99a-5p, and miR-508-5p, were related to esophageal cancer prognosis. GSEA analysis showed that the following genes may be involved in esophageal cancer prognosis: ZBTB16 may through the MTOR signaling pathway, AQP4 through the GNRH signaling pathway, ADCYAP1R1 through the PPAR signaling pathway, VIPR2 through the P53 signaling pathway and PDGFD through the PENTOSE-PHOSPHATE signaling pathway.

Identification of microRNA transcriptome throughout the lifespan of yak (Bos grunniens) corpus luteum

The corpus luteum (CL) is a temporary organ that plays a critical role for female fertility by maintaining the estrous cycle. MicroRNA (miRNA) is a class of non-coding RNAs involved in various biological processes. However, there exists limited knowledge of the role of miRNA in yak CL. In this study, we used high-throughput sequencing to study the transcriptome dynamics of miRNA in yak early (eCL), middle (mCL) and late-stage CL (lCL). A total of 6,730 miRNAs were identified, including 5,766 known and 964 novels miRNAs. Three miRNAs, including bta-miR-126-3p, bta-miR-143 and bta-miR-148a, exhibited the highest expressions in yak CLs of all the three stages. Most of the miRNAs were 20–24 nt in length and the peak was at 22 nt. Besides, most miRNAs with different lengths displayed significant uracil preference at the 5'-end. Furthermore, 1,067, 280 and 112 differentially expressed (DE) miRNAs were found in eCL vs. mCL, mCL vs. lCL, and eCL vs. lCL, respectively. Most of the DE miRNAs were down-regulated in the eCL vs. mCL and eCL vs. lCL groups, and up-regulated in the mCL vs. lCL group. A total of 18,904 target genes were identified, with 18,843 annotated. Pathway enrichment analysis of the DE miRNAs target genes illustrated that the most enriched cellular process in each group included pathways in cancer, PI3K-Akt pathway, endocytosis, and focal adhesion. A total of 20 putative target genes in 47 DE miRNAs were identified to be closely associated with the formation, function or regression of CL. Three DE miRNAs, including bta-miR-11972, novel-miR-619 and novel-miR-153, were proved to directly bind to the 3'-UTR of their predicated target mRNAs, including CDK4, HSD17B1 and MAP1LC3C, respectively. Both of these DE miRNAs and their target mRNAs exhibited dynamic expression profiles across the lifespan of yak CL. This study presents a general basis for understanding of the regulation of miRNA on yak CL and also provides a novel genetic resource for future analysis of the gene network during the estrous cycle in the yak.

MicroRNA Omics Analysis of Camellia sinesis Pollen Tubes in Response to Low-Temperature and Nitric Oxide.

Nitric oxide (NO) as a momentous signal molecule participates in plant reproductive development and responds to various abiotic stresses. Here, the inhibitory effects of the NO-dominated signal network on the pollen tube growth of Camellia sinensis under low temperature (LT) were studied by microRNA (miRNA) omics analysis. The results showed that 77 and 71 differentially expressed miRNAs (DEMs) were induced by LT and NO treatment, respectively. Gene ontology (GO) analysis showed that DEM target genes related to microtubules and actin were enriched uniquely under LT treatment, while DEM target genes related to redox process were enriched uniquely under NO treatment. In addition, the target genes of miRNA co-regulated by LT and NO are only located on the cell membrane and cell wall, and most of them are enriched in metal ion binding and/or transport and cell wall organization. Furthermore, DEM and its target genes related to metal ion binding/transport, redox process, actin, cell wall organization and carbohydrate metabolism were identified and quantified by functional analysis and qRT-PCR. In conclusion, miRNA omics analysis provides a complex signal network regulated by NO-mediated miRNA, which changes cell structure and component distribution by adjusting Ca2+ gradient, thus affecting the polar growth of the C. sinensis pollen tube tip under LT.

Mir-139-5p inhibits glioma cell proliferation and progression by targeting GABRA1

Glioma is an extremely aggressive malignant neoplasm of the central nervous system. MicroRNA (miRNA) are known to bind to specific target mRNA to regulate post-transcriptional gene expression and are, therefore, currently regarded as promising biomarkers for glioma diagnosis and prognosis. The aim of the present study was to examine the pathogenesis and potential molecular markers of glioma by comparing the differential expression of miRNA and mRNA between glioma tissue and peritumor brain tissue. We explored the impact of screened core miRNA and mRNA on cell proliferation, invasion, and migration of glioma. An miRNA expression profile dataset (GSE90603) and a transcriptome profile dataset (GSE90598) were downloaded from combined miRNA-mRNA microarray chips in the Gene Expression Omnibus (GEO) database. Overall, 59 differentially expressed miRNAs (DEMs) and 419 differentially expressed genes (DEGs) were identified using the R limma software package. FunRich software was used to predict DEM target genes and miRNA-gene pairs, and Perl software was used to find overlapping genes between DEGs and DEM target genes. There were 129 overlapping genes regulated by nine miRNAs between target genes of the DEMs and DEGs. The Chinese Glioma Genome Atlas(CGGA) was analyzed in order to identify miRNAs with diagnostic and prognostic significance. MiR-139-5p, miR-137, and miR-338-3p were validated to be significantly linked to prognosis in glioma patients. Finally, we validated that miR-139-5p affected glioma malignant biological behavior via targeting gamma-aminobutyric acid A receptor alpha 1(GABRA1) through rescue experiments. Low miR-139-5p expression was correlated with survival probability and World Health Organization (WHO) grade. MiR-139-5p overexpression inhibited cell proliferation, migration, and invasion of glioma in vitro. GABRA1 was identified as a functional downstream target of miR-139-5p. Decreased GABRA1 expression was related to similar biological roles as miR-139-5p overexpression while upregulation of GABRA1 effectively reversed the inhibition effects of miR-139-5p. These results demonstrate a novel axis for miR-139-5p/GABRA1 in glioma progression and provide potential prognostic predictors and therapeutic target for glioma patients.

Differential expression of miRNAs in the presence of B chromosome in the cichlid fish Astatotilapia latifasciata

BackgroundB chromosomes (Bs) are extra elements observed in diverse eukaryotes, including animals, plants and fungi. Although Bs were first identified a century ago and have been studied in hundreds of species, their biology is still enigmatic. Recent advances in omics and big data technologies are revolutionizing the B biology field. These advances allow analyses of DNA, RNA, proteins and the construction of interactive networks for understanding the B composition and behavior in the cell. Several genes have been detected on the B chromosomes, although the interaction of B sequences and the normal genome remains poorly understood.ResultsWe identified 727 miRNA precursors in the A. latifasciata genome, 66% which were novel predicted sequences that had not been identified before. We were able to report the A. latifasciata-specific miRNAs and common miRNAs identified in other fish species. For the samples carrying the B chromosome (B+), we identified 104 differentially expressed (DE) miRNAs that are down or upregulated compared to samples without B chromosome (B−) (p < 0.05). These miRNAs share common targets in the brain, muscle and gonads. These targets were used to construct a protein-protein-miRNA network showing the high interaction between the targets of differentially expressed miRNAs in the B+ chromosome samples. Among the DE-miRNA targets there are protein-coding genes reported for the B chromosome that are present in the protein-protein-miRNA network. Additionally, Gene Ontology (GO) terms related to nuclear matrix organization and response to stimulus are exclusive to DE miRNA targets of B+ samples.ConclusionsThis study is the first to report the connection of B chromosomes and miRNAs in a vertebrate species. We observed that the B chromosome impacts the miRNAs expression in several tissues and these miRNAs target several mRNAs involved with important biological processes.

Genome-wide identification and characterization of perirenal adipose tissue microRNAs in rabbits fed a high-fat diet.

MicroRNAs (miRNAs) are a class of endogenous single-stranded RNA molecules that play an important role in gene regulation in animals by pairing with target gene mRNA. Extensive evidence shows that miRNAs are key players in metabolic regulation and the development of obesity. However, the systemic understanding of miRNAs in the adipogenesis of obese rabbits need further investigation. Here, seven small RNA libraries from rabbits fed either a standard normal diet (SND; n=3) or high-fat diet (HFD; n=4) were constructed and sequenced. Differentially expressed (DE) miRNAs were identified using the edgeR data analysis package from R. Software miRanda and RNAhybrid were used to predict the target genes of miRNAs. To further explore the functions of DE miRNAs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed. A total of 81449996 clean reads were obtained from the seven libraries, of which, 52 known DE miRNAs (24 up-regulated, 28 down-regulated) and 31 novel DE miRNAs (14 up-regulated, 17 down-regulated) were identified. GO enrichment analysis revealed that the DE miRNAs target genes were involved in intermediate filament cytoskeleton organization, intermediate filament-based process, and α-tubulin binding. DE miRNAs were involved in p53 signaling, linoleic acid metabolism, and other adipogenesis-related KEGG pathways. Our study further elucidates the possible functions of DE miRNAs in rabbit adipogenesis, contributing to the understanding of rabbit obesity.

Genome-Wide Identification and Analysis of Nilaparvata lugens microRNAs during Challenge with the Entomopathogenic Fungus Metarhizium anisopliae.

The resistance of the notorious rice pest Nilaparvata lugens to many insecticides has caused significant concerns. Our previous study demonstrated that the fungus Metarhizium anisopliae CQMa421 shows great potential for the control of this pest, but the interactions between them are still unclear. Thus, we further investigated fungal infection-related microRNAs (miRNAs) in N. lugens during M. anisopliae CQMa421 challenge using Illumina sequencing. In this study, we constructed twenty-four small RNA libraries over different time courses (i.e., 4 h, 8 h, 16 h, and 24 h). A total of 478.62 M clean reads were collected, with each sample producing more than 13.37 M reads, after the removal of low-quality reads. We identified 2324 miRNAs and their 11,076 target genes within the twenty-four libraries by bioinformatics analysis. Differentially expressed miRNAs (DEmiRNAs), including 58 (32 upregulated vs. 26 downregulated), 62 (30 upregulated vs. 32 downregulated), 126 (71 upregulated vs. 55 downregulated), and 109 (40 upregulated vs. 69 downregulated) DEmiRNAs were identified at 4 h, 8 h, 16 h, and 24 h post-infection, respectively. We further conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis to predict the functions of all target genes of DEmiRNAs. These DEmiRNAs targets identified during 24 h of infection were primarily involved in energy metabolism, lysine degradation, the FoxO signaling pathway, ubiquitin-mediated proteolysis, the mRNA surveillance pathway, and the MAPK signaling pathway. Taken together, our results provide essential information for further study of the interactions between the entomopathogenic fungus M. anisopliae and N. lugens at the posttranscriptional level.

Bioinformatics Analysis of Dysregulated MicroRNA-Messenger RNA Networks in Small Cell Lung Cancer

The present study aimed to identify a key module of differentially expressed miRNAs (DE-miRNAs) together with the corresponding differentially expressed mRNAs (DE-mRNAs) within small cell lung cancer (SCLC). Linear models were applied to ascertain the DE-miRNAs and DE-mRNAs in SCLC versus matched non-carcinoma samples obtained from the RNA expression datasets of GSE19945, GSE74190 and GSE6044. The common DE-miRNAs were identified using the Venn plot. Then, 3 databases were used to retrieve the DE-miRNAs target genes, and the intersection was taken for validating the shared target genes. Besides, Cytoscape was utilized for constructing the miRNAmRNA network for SCLC. Finally, a key module of five DE-miRNAs and four hub genes was determined based on the degree. In addition, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted for exploring those hub genes in terms of their functions along with the involved signal transduction pathways. Altogether 106 shared DE-miRNAs were identified, which were used to predict 63 common target genes. In addition, a key module of five DE-miRNAs (hsa-miR-17-5p, hsa-miR-20a-5p, hsa-miR-20b-5p, hsa-miR-93-5p and hsa-miR- 106b-5p) and four hub genes (SOX4, DPYSL2, TGFBR2 and F3) were extracted from the miRNAmRNA network according to their degree. Finally, the hub genes were subjected to GO as well as KEGG analysis, which revealed that cell cycle G1/S phase transition, the extracellular matrix, and cellular senescence signaling pathways exerted vial parts during SCLC progression. A key module of five DE-miRNAs and four hub genes may be potentially used as clinical biomarkers to predict SCLC.

Abstract P737: Human Homologs Of Differentially Expressed Plasma MicroRNAs In Preclinical Murine Models Target Fundamental Mechanisms Of Cerebral Cavernous Malformation

Introduction: There is a current need for sensitive and specific biomarkers of Cerebral Cavernous Malformation (CCM) that can be readily translated from preclinical to human models to accurately diagnose and monitor disease states and response to novel therapeutics. MiRNAs are small non-coding RNAs that influence gene expression and whose levels can be affected by disease states. We hypothesize that there are human homologs of differentially expressed (DE) miRNA in the plasma of CCM murine models that can be identified in CCM patients. We further hypothesize that these miRNAs have gene targets within previously published CCM transcriptomes, mechanistically linking them to CCM disease. Methods: Plasma miRNAs from homozygous and heterozygous Ccm1 and Ccm3 mice, as well as their respective wild type controls were sequenced and analyzed. Putative gene targets of DE miRNAs [p&lt;0.05, false discovery rate (FDR) corrected] were queried in previously published mouse CCM transcriptomes. The human homologs of the DE miRNAs in the plasma of Ccm1 and Ccm3 mouse models were identified and assessed in the plasma of healthy controls (n=13), CCM1 (n=11), and CCM3 (n=11) patients using RT-qPCR. Results: 5 miRNAs in homozygous and 10 in heterozygous for Ccm1 , while 45 in homozygous and 2 in heterozygous for Ccm3 were DE in the plasma of mouse models (p&lt;0.05, FDR corrected), had gene targets within CCM mouse transcriptomes, and have a human homolog. Preliminary results show mmu-miR-375-3p as DE in both Ccm1 +/- and Ccm3 -/- mice. RT-qPCR assays show that plasma relative quantification values of the human homolog hsa-miR-375-3p were higher in CCM3 than in CCM1 patients ( p &lt;0.001) and in healthy controls ( p &lt;0.05). Conclusion: DE plasma miRNAs identified in mouse models with human homologs and putative gene targets mechanistically implicated in CCM disease may be used as candidate biomarkers for specific clinical contexts.

Identification and characterization of sex-biased and differentially expressed miRNAs in gonadal developments of the Chinese mitten crab, Eriocheir sinensis.

MicroRNA (miRNA) is a posttranscriptional downregulator that plays a vital role in a wide variety of biological processes. In this study, we constructed five ovarian and testicular small RNA libraries using two somatic libraries as reference controls for the identification of sex-biased miRNAs and gonadal differentially expressed miRNAs (DEMs) of the Chinese mitten crab, Eriocheir sinensis. A total of 535 known and 243 novel miRNAs were identified, including 312 sex-biased miRNAs and 402 gonadal DEMs. KEGG pathway analysis showed that DEM target genes were statistically enriched in MAPK, Wnt, and GnRH signaling pathway, and so on. A number of the sex-biased miRNAs target genes associated with sex determination/differentiation, such as IAG, Dsx, Dmrt1, and Fem1, while others target the genes related to gonadal development, such as P450s, Wnt, Ef1, and Tra-2c. Dual-luciferase reporter assay in vitro further confirmed that miR-34 and let-7b can downregulate IAG expression, miR-9-5p, let-7d, let-7b, and miR-8915 can downregulate Dsx. Taken together, these data strongly suggest a potential role for the sex-biased miRNAs in sex determination/differentiation and gonadal development in the crab.

Identification of miRNA-Target Gene Pairs in the Parietal and Frontal Lobes of the Brain in Patients with Alzheimer's Disease Using Bioinformatic Analyses.

Alzheimer's disease (AD) is a growing health concern worldwide. MicroRNAs (miRNAs) have been extensively studied in many diseases, including AD. To identify differentially expressed miRNAs (DEmiRNAs) and genes specific to AD, we used bioinformatic analyses to investigate candidate miRNA-mRNA pairs involved in the pathogenesis of AD. We focused on differentially expressed genes (DEGs) that are targets of DEmiRNAs. The GEO2R tool and the HISAT2-DESeq2 software were used to identify DEmiRNAs and DEGs. Bioinformatic tools available online, such as TAM and the Database for Annotation, Visualization and Integrated Discovery (DAVID), were used to perform functional annotation and enrichment analysis. Targets of miRNAs were predicted using the miRTarBase. The Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape, which are available online, were utilized to construct protein-protein interaction (PPI) networks and identify hub genes. Furthermore, transcription factors (TFs) encoded by the DEGs were predicted using the TransmiR database and TF-miRNA-mRNA networks were constructed. Finally, the expression profile of a hub gene in peripheral blood mononuclear cells was compared between healthy individuals and AD patients. We identified 26 correlated miRNA-mRNA pairs. In the parietal lobe, miRNA-mRNA pairs involved in protein folding were enriched, and in the frontal lobe, miRNA-mRNA pairs involved in synaptic transmission, abnormal protein degradation, and apoptosis were enriched. In addition, HSP90AB1 in peripheral blood mononuclear cells was found to be significantly downregulated in AD patients, and this was consistent with its expression profile in the parietal lobe of AD patients. Our results provide brain region-specific changes in miRNA-mRNA associations in AD patients, further our understanding of potential underlying molecular mechanisms of AD, and reveal promising diagnostic and therapeutic targets for AD.

Identification of microRNA-16-5p and microRNA-21-5p in feces as potential noninvasive biomarkers for inflammatory bowel disease.

Background: Inflammatory bowel disease (IBD) is a chronic idiopathic gastrointestinal disease. Increasing evidence suggests that microRNAs (miRNAs) may participate in the pathophysiology of IBD.Methods: A miRCURY™ LNA Array and in situ hybridization were employed to screen for differentially expressed miRNAs (DEMs) in fecal specimens from 41 IBD patients (22 ulcerative colitis (UC), 19 Crohn’s disease (CD)) and 23 healthy controls (HC). RT-qPCR was performed to confirm the findings. The DEMs target genes and corresponding biological functions were predicted by bioinformatics analysis.Results: Compared with HC, miR-16-5p in the feces was up-regulated both in UC and CD patients (p < 0.01), while miR-21-5p was up-regulated only in UC patients (p < 0.01). TargetScan 7.2, miRWalk, and miRDB were used to predict 216 public target genes of miR-16-5p and miR-21-5p, and six hub genes (PIK3R1, GRB2, SUZ12, NTRK2, Smurf2, and WWP1) were analyzed using the STRING database and Cytoscape. All the hub genes promote the occurrence and development of IBD-related colorectal cancer.Conclusions: The elevated levels of miR-16-5p and miR-21-5p in feces of IBD patients have to guide significance for the noninvasive clinical diagnosis of IBD and have a warning effect on the occurrence of IBD-related colorectal cancer.

Integrative Analysis Extracts a Core ceRNA Network of the Fetal Hippocampus With Down Syndrome.

Accumulating evidence suggests that circular RNAs (circRNAs)—miRNA–mRNA ceRNA regulatory network—may play an important role in neurological disorders, such as Alzheimer’s disease (AD). Interestingly, neuropathological changes that closely resemble AD have been found in nearly all Down syndrome (DS) cases > 35 years. However, few studies have reported circRNA transcriptional profiling in DS cases, which is caused by a chromosomal aberration of trisomy 21. Here, we characterized the expression profiles of circRNAs in the fetal hippocampus of DS patients (n = 8) and controls (n = 6) by using microarray. MiRNA, mRNA expression profiling of DS from our previous study and scRNA-seq data describing normal fetal hippocampus development (GEO) were also integrated into the analysis. The similarity between circRNAs/genes with traits/cell-types was calculated by weighted correlation network analysis (WGCNA). miRanda and miRWalk2 were used to predict ceRNA network interactions. We identified a total of 7,078 significantly differentially expressed (DE) circRNAs, including 2,637 upregulated and 4,441 downregulated genes, respectively. WGCNA obtained 15 hub circRNAs and 6 modules with cell type–specific expression patterns among scRNA-seq data. Finally, a core ceRNA network was constructed by 14 hub circRNAs, 17 DE miRNA targets and 245 DE mRNA targets with a cell type–specific expression pattern annotation. Known functional molecules in DS or neurodegeneration (e.g., miR-138, OLIG1, and TPM2) were also included in this network. Our findings are the first to delineate the landscape of circRNAs in DS and the first to effectively integrate ceRNA regulation with scRNA-seq data. These data may provide a valuable resource for further research on the molecular mechanisms or therapeutic targets underlying DS neuropathy.

Identification and Validation of Potential miRNAs, as Biomarkers for Sepsis and Associated Lung Injury: A Network-Based Approach.

Sepsis is a dysregulated immune response disease affecting millions worldwide. Delayed diagnosis, poor prognosis, and disease heterogeneity make its treatment ineffective. miRNAs are imposingly involved in personalized medicine such as therapeutics, due to their high sensitivity and accuracy. Our study aimed to reveal the biomarkers that may be involved in the dysregulated immune response in sepsis and lung injury using a computational approach and in vivo validation studies. A sepsis miRNA Gene Expression Omnibus (GEO) dataset based on the former analysis of blood samples was used to identify differentially expressed miRNAs (DEMs) and associated hub genes. Sepsis-associated genes from the Comparative Toxicogenomics Database (CTD) that overlapped with identified DEM targets were utilized for network construction. In total, 317 genes were found to be regulated by 10 DEMs (three upregulated, namely miR-4634, miR-4638-5p, and miR-4769-5p, and seven downregulated, namely miR-4299, miR-451a, miR181a-2-3p, miR-16-5p, miR-5704, miR-144-3p, and miR-1290). Overall hub genes (HIP1, GJC1, MDM4, IL6R, and ERC1) and for miR-16-5p (SYNRG, TNRC6B, and LAMTOR3) were identified based on centrality measures (degree, betweenness, and closeness). In vivo validation of miRNAs in lung tissue showed significantly downregulated expression of miR-16-5p corroborating with our computational findings, whereas expression of miR-181a-2-3p and miR-451a were found to be upregulated in contrast to the computational approach. In conclusion, the differential expression pattern of miRNAs and hub genes reported in this study may help to unravel many unexplored regulatory pathways, leading to the identification of critical molecular targets for increased prognosis, diagnosis, and drug efficacy in sepsis and associated organ injuries.

Application of microRNA profiling to understand sevoflurane-induced adverse effects on developing monkey brain

We have described that prolonged sevoflurane exposure at a clinically-relevant concentration of 2.5 % caused neuronal cell death in the developing monkey brain. Postnatal day 5 or 6 rhesus monkeys (n = 3) were exposed to 2.5 % sevoflurane for 8 h. Monkeys kept at environmental conditions in the procedure room served as controls (n = 3). Brain tissues were harvested four hours after sevoflurane exposure for histological analysis, and RNA or protein extraction. MicroRNA (miRNA) profiling on the frontal cortex of monkey brains was performed using next-generation sequencing. 417 miRNAs were identified in the frontal cortex, where most neuronal cell death was observed. 7 miRNAs were differentially expressed in frontal cortex after sevoflurane exposure. Five of these were expressed at significantly lower levels than controls and the other two miRNAs were expressed significantly higher. These differentially expressed miRNAs (DEMs) were then loaded to the Ingenuity Pathway Analysis database for pathway analysis, in which targeting information was available for 5 DEMs. The 5 DEMs target 2,919 mRNAs which are involved in pathways that contribute to various cellular functions. Of note, 78 genes that are related to axon guidance signaling were targeted, suggesting that development of the neural circuit may be affected after sevoflurane exposure. Such changes may have long-term effects on brain development and function. These findings are supplementary to our previous observations and provide more evidence for better understanding the adverse effects of sevoflurane on the developing brain after an 8 -h exposure.

Genome-wide identification of low phosphorus responsive microRNAs in two soybean genotypes by high-throughput sequencing.

MicroRNAs (miRNAs) have been reported to be correlated with various stress responses in soybean, but only a few miRNAs have been demonstrated to respond to low phosphorus (LP) stress. To unravel the response mechanisms of miRNAs to low-P stress, the roots of two representative soybean genotypes with different P efficiency, Nannong94-156 (a LP-tolerant genotype) and Bogao (a LP-sensitive genotype), were used for the construction of RNA sequencing (RNA-seq) libraries under low/normal-P treatment by high-throughput sequencing. In total, 603 existing miRNAs and 1699 novel miRNAs belonging to 248 and 1582 families in all samples were identified, respectively. Among these miRNAs, 777 miRNAs were differentially expressed (DE) across different P levels and genotypes. Furthermore, putative targets of DE miRNAs were predicted, and these miRNAs mainly targeted ERF (ethylene responsive factor), auxin response factors (ARF), zinc finger protein, MYB, and NAC domain transcription factors. Gene ontology (GO) analysis showed that targets of DE miRNAs were significantly enriched in binding, metabolic processes, biological regulation, response to stress, and phosphorus metabolic processes. In addition, the expression profiles of chosen P-responsive miRNAs and target genes were validated by quantitative real-time PCR (qRT-PCR). Our study focused on genome-wide miRNA identification in two representative soybean genotypes under low-P stress. Overall, the DE miRNAs across different P levels and genotypes and their putative target genes will provide useful information for further study of miRNAs mediating low-P response and facilitate improvements in soybean breeding.

Identification of microRNAs enriched in exosomes in human pericardial fluid of patients with atrial fibrillation based on bioinformatic analysis.

BackgroundAtrial fibrillation (AF) is related to structural and electrical atria remodeling. Atrial fibrosis development and progression is characteristic of structural remodeling and is taken as the AF perpetuation substrate. Increasing evidence has confirmed that microRNAs (miRNAs) are associated with AF, including cardiac fibrosis.MethodsPericardial fluid (PF) samples were collected from nine adult patients who had congenital heart disease with persistent AF or sinus rhythm (SR) undergoing surgery. Abnormally expressed miRNAs were acquired, and P<0.05 and fold change >2 were taken as the thresholds of differentially expressed miRNAs (DE-miRNAs). The predicted target genes were obtained by miRTarBase. The Database for Annotation, Visualization and Integrated Discovery was used to annotate functions and analyze pathway abundance for latent targets of DE-miRNAs. STRING database was applied to construct a protein–protein interplay (PPI) network, and Cytoscape software was used to visualize the miRNA-hub gene-Kyoto Encyclopedia of Genes and Genomes (KEGG) network. DE-miRNA expressions were evaluated by quantitative polymerase chain reaction (qPCR).ResultsFifty-five exosomal DE-miRNAs were found between the AF and SR samples; these included 24 miRNAs that were upregulated and 31 that were downregulated. For the top 3 downregulated miRNAs (miR-382-3p, miR-3126-5p, and miR-450a-2-3p) 283 predicted target genes were identified, and were implicated in cardiac fibrosis-related pathways, including the hypoxia-inducible factor-1 (HIF1), mitogen-activated protein kinase (MAPK), and adrenergic and insulin pathways. The top 10 hub genes in the PPI network, including mitogen-activated protein kinase 1 (MAPK1) and AKT serine/threonine kinase 1 (AKT1), were identified as hub genes. By establishing the miRNA-hub gene-KEGG network, we observed that these hub genes, which were regulated by miR-382-3p, miR-3126-5p, and miR-450a-2-3p, were involved in many KEGG pathways associated with cardiac fibrosis, such as the AKT1/glycogen synthase kinsase-3β (GSK-3β) and transforming growth factor-β (TGF-β)/MAPK1 pathways.ConclusionsThe findings of the present study suggest that miR-382-3p, miR-450a-2-3p, and miR-3126-5p contained in exosomes in human PF are pivotal in the progression of AF. The results of qPCR showed that miR-382-3p was consistent with our sequencing data, which indicates its potential value as a therapeutic target for AF.

Identification of hub genes in unstable atherosclerotic plaque by conjoint analysis of bioinformatics

AimsUnstable atherosclerotic plaque is the main pathological basis of acute coronary syndrome, which is the leading cause of death and disability worldwide. Therefore, we combined multiple bioinformatics tools to identify key genes related to unstable plaque. Main methodsGSE94605 contained 7 plasma sample pools of 175 healthy and 6 sample pools of 150 unstable angina pectoris (UAP) patients, and detected with miRNA array while GSE60993 collected peripheral blood from 7 normal and 9 UAP, and detected with mRNA array. GSE120521 collected carotid plaques from 4 patients and dissected in stable and unstable regions, then detected with RNA-seq. Differentially expressed miRNAs (DEMs) and genes (DEGs) in UAP were re-analyzed. Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) and Protein-protein interaction (PPI) network were applied on top 10 up-regulated or down-regulated DEMs targets, and whole DEGs. MiRNAs-mRNAs network was constructed with these DEMs and DEGs, and the expression profile of genes within the network was finally validated in GSE120521. Key findingsTotally, 263 up-regulated and 201 down-regulated DEMs were identified in GSE94605, and 78 up-regulated and 29 down-regulated DEGs were identified in GSE60993. Subsequently, a miRNAs-mRNAs network was constructed with 6 up-regulated miRNAs targeted to 12 down-regulated genes, and 4 down-regulated miRNAs targeted to 8 up-regulated genes. Finally, MORF4L2, RAB3IL1 and MMP9 within the network were considered as hub genes in unstable plaque progression after being validated in GSE120521. SignificanceThese 3 genes may provide new targets for diagnosis and therapy of unstable atherosclerotic plaque.

Characterization of Differentially Expressed miRNAs and Their Predicted Target Transcripts during Smoltification and Adaptation to Seawater in Head Kidney of Atlantic Salmon.

Smoltification and early seawater phase are critical developmental periods with physiological and biochemical changes in Atlantic salmon that facilitates survival in saltwater. MicroRNAs (miRNAs) are known to have important roles in development, but whether any miRNAs are involved in regulation of gene expression during smoltification and the adaption to seawater is largely unknown. Here, small RNA sequencing of materials from head kidney before, during smoltification and post seawater transfer were used to study expression dynamics of miRNAs, while microarray analysis was applied to study mRNA expression dynamics. Comparing all timepoints, 71 miRNAs and 2709 mRNAs were identified as differentially expressed (DE). Hierarchical clustering analysis of the DE miRNAs showed three major clusters with characteristic expression changes. Eighty-one DE mRNAs revealed negatively correlated expression patterns to DE miRNAs in Cluster I and III. Furthermore, 42 of these mRNAs were predicted as DE miRNA targets. Gene enrichment analysis of negatively correlated target genes showed they were enriched in gene ontology groups hormone biosynthesis, stress management, immune response, and ion transport. The results strongly indicate that post-transcriptional regulation of gene expression by miRNAs is important in smoltification and sea water adaption, and this study identifies several putative miRNA-target pairs for further functional studies.

Construction of Potential miRNA-mRNA Regulatory Network in COPD Plasma by Bioinformatics Analysis.

BackgroundChronic obstructive pulmonary disease (COPD) has become a major cause of morbidity and mortality worldwide. Increasing evidence indicates that aberrantly expressed microRNAs (miRNAs) are involved in the pathogenesis of COPD. However, an integrative exploration of miRNA–mRNA regulatory network in COPD plasma remains lacking.MethodsThe microarray datasets GSE24709, GSE61741, and GSE31568 were downloaded from the GEO database and analyzed using GEO2R tool to identify differentially expressed miRNAs (DEMs) between COPD and normal plasma. The consistently changing miRNAs in the three datasets were screened out as candidate DEMs. Potential upstream transcription factors and downstream target genes of candidate DEMs were predicted by FunRich and miRNet, respectively. Next, GO annotation and KEGG pathway enrichment analysis for target genes were performed using DAVID. Then, PPI and DEM-hub gene network were constructed using the STRING database and Cytoscape software. Finally, GSE56768 was used to evaluate the hub gene expressions.ResultsA total of nine (six upregulated and three downregulated) DEMs were screened out in the above three datasets. SP1 was predicted to potentially regulate most of the downregulated DEMs, while YY1 and E2F1 could regulate both upregulated and downregulated DEMs. 1139 target genes were then predicted, including 596 upregulated DEM target genes and 543 downregulated DEM target genes. Target genes of DEMs were mainly enriched in PI3K/Akt signaling pathway, mTOR signaling pathway, and autophagy. Through the DEM-hub gene network construction, most of the hub genes were found to be potentially modulated by miR-497-5p, miR-130b-5p, and miR-126-5p. Among the top 12 hub genes, MYC and FOXO1 expressions were consistent with that in the GSE56768 dataset.ConclusionIn the study, potential miRNA–mRNA regulatory network was firstly constructed in COPD plasma, which may provide a new insight into the pathogenesis and treatment of COPD.