Identifying MicroRNA and mRNA Expression Profiles in Embryonic Stem Cells Derived from Parthenogenetic, Androgenetic and Fertilized Blastocysts

MicroRNAs (miRNAs) are a class of highly conserved small non-coding RNA molecules that play a pivotal role in several cellular functions. In this study, miRNA and messenger RNA (mRNA) profiles were examined by Illumina microarray in mouse embryonic stem cells (ESCs) derived from parthenogenetic, androgenetic, and fertilized blastocysts. The global analysis of miRNA-mRNA target pairs provided insight into the role of miRNAs in gene expression. Results showed that a total of 125 miRNAs and 2394 mRNAs were differentially expressed between androgenetic ESCs (aESCs) and fertilized ESCs (fESCs), a total of 42 miRNAs and 87 mRNAs were differentially expressed between parthenogenetic ESCs (pESCs) and fESCs, and a total of 99 miRNAs and 1788 mRNAs were differentially expressed between aESCs and pESCs. In addition, a total of 575, 5 and 376 miRNA-mRNA target pairs were observed in aESCs vs. fESCs, pESCs vs. fESCs, and aESCs vs. pESCs, respectively. Furthermore, 15 known imprinted genes and 16 putative uniparentally expressed miRNAs with high expression levels were confirmed by both microarray and real-time RT-PCR. Finally, transfection of miRNA inhibitors was performed to validate the regulatory relationship between putative maternally expressed miRNAs and target mRNAs. Inhibition of miR-880 increased the expression of Peg3, Dyrk1b, and Prrg2 mRNA, inhibition of miR-363 increased the expression of Nfat5 and Soat1 mRNA, and inhibition of miR-883b-5p increased Nfat5, Tacstd2, and Ppapdc1 mRNA. These results warrant a functional study to fully understand the underlying regulation of genomic imprinting in early embryo development.     

Next-Generation Sequencing Analysis Reveals Differential Expression Profiles of MiRNA-mRNA Target Pairs in KSHV-Infected Cells

Kaposi’s sarcoma associated herpesvirus (KSHV) causes several tumors, including primary effusion lymphoma (PEL) and Kaposi’s sarcoma (KS). Cellular and viral microRNAs (miRNAs) have been shown to play important roles in regulating gene expression. A better knowledge of the miRNA-mediated pathways affected by KSHV infection is therefore important for understanding viral infection and tumor pathogenesis. In this study, we used deep sequencing to analyze miRNA and cellular mRNA expression in a cell line with latent KSHV infection (SLKK) as compared to the uninfected SLK line. This approach revealed 153 differentially expressed human miRNAs, eight of which were independently confirmed by qRT-PCR. KSHV infection led to the dysregulation of ~15% of the human miRNA pool and most of these cellular miRNAs were down-regulated, including nearly all members of the 14q32 miRNA cluster, a genomic locus linked to cancer and that is deleted in a number of PEL cell lines. Furthermore, we identified 48 miRNAs that were associated with a total of 1,117 predicted or experimentally validated target mRNAs; of these mRNAs, a majority (73%) were inversely correlated to expression changes of their respective miRNAs, suggesting miRNA-mediated silencing mechanisms were involved in a number of these alterations. Several dysregulated miRNA-mRNA pairs may facilitate KSHV infection or tumor formation, such as up-regulated miR-708-5p, associated with a decrease in pro-apoptotic caspase-2 and leukemia inhibitory factor LIF, or down-regulated miR-409-5p, associated with an increase in the p53-inhibitor MDM2. Transfection of miRNA mimics provided further evidence that changes in miRNAs are driving some observed mRNA changes. Using filtered datasets, we also identified several canonical pathways that were significantly enriched in differentially expressed miRNA-mRNA pairs, such as the epithelial-to-mesenchymal transition and the interleukin-8 signaling pathways. Overall, our data provide a more detailed understanding of KSHV latency and guide further studies of the biological significance of these changes.     

microRNA expression profiles and the potential competing endogenous RNA networks in NELL-1-induced human adipose-derived stem cell osteogenic differentiation

Studies have indicated that Nel-like molecule-1 (NELL-1) was an osteoblast-specific cytokine and some specific microRNAs (miRNAs) could serve as competing endogenous RNA (ceRNA) to partake in osteogenic differentiation of human adipose-derived stem cells (hASCs). The aim of this study was to explore the potential functional mechanisms of recombinant human NELL-1 protein (rhNELL-1) during hASCs osteogenic differentiation. rhNELL-1 was added to osteogenic medium to activate osteogenic differentiation of hASCs. High-throughput RNA sequencing (RNA-Seq) was performed and validated by real-time quantitative polymerase chain reaction. Gene ontology functional annotation and Kyoto Encyclopedia of Genes and Genomes pathway analysis were performed to detect the functions of differentially expressed miRNAs and genes. Coding-noncoding gene co-expression network and ceRNA networks were constructed to predict the potential regulatory role of miRNAs. A total of 1010 differentially expressed miRNAs and 1762 differentially expressed messenger RNAs (mRNAs) were detected. miRNA-370-3p, bone morphogenetic protein 2 (BMP2), and parathyroid hormone like hormone (PTHLH) were differentially expressed during NELL-1-induced osteogenesis. Bioinformatic analyses demonstrated that these differentially expressed miRNAs and mRNAs enriched in Rap1 signaling pathway, PI3K-Akt signaling pathway, p53 signaling pathway, Glucagon signaling pathway, and hypoxia-inducible factor-1 signaling pathway, which were important pathways related to osteogenic differentiation. In addition, miRNA-370-3p and has-miR-485-5p were predicted to interact with circ0001543, circ0002405, and ENST00000570267 in ceRNA networks. Based on the gain or loss of functional experiments by transfection, the results showed that miR-370-3p was a key regulator in osteogenic differentiation by targeting BMP2 and disturbing the expression of PTHLH, and participated in NELL-1-stimulated osteogenesis. The present study provided the primary data and evidence for further exploration on the roles of miRNAs and ceRNAs during NELL-1-induced ossification of hASCs.     

The RNA expression signature of the HepG2 cell line as determined by the integrated analysis of miRNA and mRNA expression profiles

Understanding miRNAs' regulatory networks and target genes could facilitate the development of therapies for human diseases such as cancer. Although much useful gene expression profiling data for tumor cell lines is available, microarray data for miRNAs and mRNAs in the human HepG2 cell line have only been compared with that of other cell lines separately. The relationship between miRNAs and mRNAs in integrated expression profiles for HepG2 cells is still unknown. To explore the miRNA–mRNA correlations in hepatocellular carcinoma (HCC) cells, we performed miRNA and mRNA expression profiling in HepG2 cells and normal liver HL-7702 cells at the genome scale using next-generation sequencing technology. We identified 193 miRNAs that are differentially expressed in these two cell lines. Of these, 89 miRNAs were down-regulated in HepG2 cells compared with HL-7702 cells, while 104 miRNAs were up-regulated. We also observed 3035 mRNAs that are significantly dys-regulated in HepG2 cells. We then performed an integrated analysis of the expression data for differentially expressed miRNAs and mRNAs and found several miRNA–mRNA pairs that are significantly correlated in HepG2 cells. Further analysis suggested that these differentially expressed genes were enriched in four tumorigenesis-related signaling pathways, namely, ErbB, JAK–STAT, mTOR, and WNT, which until now had not been fully reported. Our results could be helpful in understanding the mechanisms of HCC occurrence and development.     

Integrated analysis of microRNA expression and mRNA transcriptome in lungs of avian influenza virus infected broilers

ABSTRACT: BACKGROUND: Avian influenza virus (AIV) outbreaks are worldwide threats to both poultry and humans. Our previous study suggested microRNAs (miRNAs) play significant roles in the regulation of host response to AIV infection in layer chickens. The objective of this study was to test the hypothesis if genetic background play essential role in the miRNA regulation of AIV infection in chickens and if miRNAs that were differentially expressed in layer with AIV infection would be modulated the same way in broiler chickens. Furthermore, by integrating with parallel mRNA expression profiling, potential molecular mechanisms of host response to AIV infection can be further exploited. RESULTS: Total RNA isolated from the lungs of non-infected and low pathogenic H5N3 infected broilers at four days post-infection were used for both miRNA deep sequencing and mRNA microarray analyses. A total of 2.6M and 3.3M filtered high quality reads were obtained from infected and non-infected chickens by Solexa GA-I Sequencer, respectively. A total of 271 miRNAs in miRBase 16.0 were identified and one potential novel miRNA was discovered. There were 121 miRNAs differentially expressed at the 5% false discovery rate by Fisher's exact test. More miRNAs were highly expressed in infected lungs (108) than in non-infected lungs (13), which was opposite to the findings in layer chickens. This result suggested that a different regulatory mechanism of host response to AIV infection mediated by miRNAs might exist in broiler chickens. Analysis using the chicken 44K Agilent microarray indicated that 508 mRNAs (347 down-regulated) were differentially expressed following AIV infection. CONCLUSION: A comprehensive analysis combining both miRNA and targeted mRNA gene expression suggests that gga-miR-34a, 122-1, 122-2, 146a, 155, 206, 1719, 1594, 1599 and 451, and MX1, IL-8, IRF-7, TNFRS19 are strong candidate miRNAs or genes involved in regulating the host response to AIV infection in the lungs of broiler chickens. Further miRNA or gene specific knock-down assay is warranted to elucidate underlying mechanism of AIV infection regulation in the chicken.     

Analysis of deep sequencing exosome-microRNA expression profile derived from CP-II reveals potential role of gga-miRNA-451 in inflammation

Mycoplasma gallisepticum (MG) can cause chronic respiratory disease (CRD) in chickens. While several studies have reported the inflammatory functions of microRNAs during MG infection, the mechanism by which exosomal miRNAs regulate MG-induced inflammation remains to be elucidated. The expression of exosome-microRNA derived from MG-infected chicken type II pneumocytes (CP-II) was screened, and the target genes and function of differentially expressed miRNAs (DEGs) were predicted. To verify the role of exosomal gga-miR-451, Western blot, ELISA and RT-qPCR were used in this study. The results showed that a total of 722 miRNAs were identified from the two exosomal small RNA (sRNA) libraries, and 30 miRNAs (9 up-regulated and 21 down-regulated) were significantly differentially expressed. The target miRNAs were significantly enriched in the treatment group, such as cell cycle, Toll-like receptor signalling pathway and MAPK signalling pathway. The results have also confirmed that gga-miR-451-absent exosomes derived from MG-infected CP-II cells increased inflammatory cytokine production in chicken fibroblast cells (DF-1), and wild-type CP-II cell–derived exosomes displayed protective effects. Collectively, our work suggests that exosomes from MG-infected CP-II cells alter the dynamics of the DF-1 cells, and may contribute to pathology of the MG infection via exosomal gga-miR-451 targeting YWHAZ involving in inflammation.     

ceRNA network construction and comparison of gastric cancer with or without Helicobacter pylori infection

Gastric cancer (GC) is a lethal disease, and among its variety of etiological factors, Helicobacter pylori (H. pylori) infection is the strongest risk factor. However, the genetic and molecular mechanisms underlying H. pylori-related GC need further elucidation. We investigated the competing endogenous RNA (ceRNA) network differences between H. pylori (+) and H. pylori (−) GC. The long noncoding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) expression data from 32 adjacent noncancerous samples and 18 H. pylori (+) and 141 H. pylori (−) stomach adenocarcinoma samples were downloaded from the TCGA database. After construction of lncRNA–miRNA–mRNA ceRNA networks of H. pylori (+) and H. pylori (−) GC, Panther and Kobas databases were used to analyze the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Finally, survival analysis was used to discover the key genes. In H. pylori (+) GC, we identified a total of 1,419 lncRNAs, 82 miRNAs, and 2,501 mRNAs with differentially expressed profiles. In H. pylori (−) GC, 2,225 lncRNAs, 130 miRNAs, and 3,146 mRNAs were differentially expressed. Furthermore, three unique pathways (cytokine–cytokine receptor interaction, HIF-1 signaling pathway, and Wnt signaling pathway) were enriched in H. pylori (+) GC. According to the overall survival analysis, three lncRNAs (AP002478.1, LINC00111, and LINC00313) and two mRNAs (MYB and COL1A1) functioned as prognostic biomarkers for patients with H. pylori (+) GC. In conclusion, our study has identified the differences in ceRNA regulatory networks between H. pylori (+) and H. pylori (−) GC and provides a rich candidate reservoir for future studies.     

Expression and localization of microRNAs in perinatal rat pancreas: role of miR-21 in regulation of cholesterol metabolism

Objective

To investigate the expression of pancreatic microRNAs (miRNAs) during the period of perinatal beta-cell expansion and maturation in rats, determine the localization of these miRNAs and perform a pathway analysis with predicted target mRNAs expressed in perinatal pancreas.

Research Design and Methods

RNA was extracted from whole pancreas at embryonic day 20 (E20), on the day of birth (P0) and two days after birth (P2) and hybridized to miRNA microarrays. Differentially expressed miRNAs were verified by northern blotting and their pancreatic localization determined by in situ hybridization. Pathway analysis was done using regulated sets of mRNAs predicted as targets of the miRNAs. Possible target genes were tested using reporter-gene analysis in INS-1E cells.

Results

Nine miRNAs were differentially expressed perinatally, seven were confirmed to be regulated at the level of the mature miRNA. The localization studies showed endocrine localization of six of these miRNAs (miR-21, -23a, -29a, -125b-5p, -376b-3p and -451), and all were expressed in exocrine cells at one time point at least. Pathways involving metabolic processes, terpenoid and sterol metabolism were selectively affected by concomitant regulation by miRNAs and mRNAs, and Srebf1 was validated as a target of miR-21.

Conclusions

The findings suggest that miRNAs are involved in the functional maturation of pancreatic exocrine and endocrine tissue following birth. Pathway analysis of target genes identify changes in sterol metabolism around birth as being selectively affected by differential miRNA expression during this period.     

Integrated miRNA-mRNA Analysis Revealing the Potential Roles of miRNAs in Chordomas

Introduction

Emerging evidence suggests that microRNAs (miRNAs) are crucially involved in tumorigenesis and that paired expression profiles of miRNAs and mRNAs can be used to identify functional miRNA-target relationships with high precision. However, no studies have applied integrated analysis to miRNA and mRNA profiles in chordomas. The purpose of this study was to provide insights into the pathogenesis of chordomas by using this integrated analysis method.

Methods

Differentially expressed miRNAs and mRNAs of chordomas (n = 3) and notochord tissues (n = 3) were analyzed by using microarrays with hierarchical clustering analysis. Subsequently, the target genes of the differentially expressed miRNAs were predicted and overlapped with the differentially expressed mRNAs. Then, GO and pathway analyses were performed for the intersecting genes.

Results

The microarray analysis indicated that 33 miRNAs and 2,791 mRNAs were significantly dysregulated between the two groups. Among the 2,791 mRNAs, 911 overlapped with putative miRNA target genes. A pathway analysis showed that the MAPK pathway was consistently enriched in the chordoma tissue and that miR-149-3p, miR-663a, miR-1908, miR-2861 and miR-3185 likely play important roles in the regulation of MAPK pathways. Furthermore, the Notch signaling pathway and the loss of the calcification or ossification capacity of the notochord may also be involved in chordoma pathogenesis.

Conclusion

This study provides an integrated dataset of the miRNA and mRNA profiles in chordomas, and the results demonstrate that not only the MAPK pathway and its related miRNAs but also the Notch pathway may be involved in chordoma development. The occurrence of chordoma may be associated with dysfunctional calcification or ossification of the notochord.     

Novel transcriptome data analysis implicates circulating microRNAs in epigenetic inheritance in mammals

Experimental evidence supports a role of mobile small non-coding RNAs in mediating soma to germline hereditary information transfer in epigenetic inheritance in plants and worms. Similar evidence in mammals has not been reported so far. In this bioinformatic analysis, differentially expressed microRNAs (miRNAs) or mRNAs reported previously in genome level expression profiling studies related to or relevant in epigenetic inheritance in mammals were examined for circulating miRNA association. The reported sets of differentially expressed miRNAs or miRNAs that are known to target the reported sets of differentially expressed genes, in that order, showed enrichment of circulating miRNAs across environmental factors, tissues, life cycle stages, generations, genders and species. Circulating miRNAs commonly representing the expression profiles enriched various epigenetic processes. These results provide bioinformatic evidence for a role of circulating miRNAs in epigenetic inheritance in mammals.     

The novel circular RNA CircMef2c is positively associated with muscle growth in Nile tilapia

Muscle growth in teleosts is a complex biological process orchestrated by numerous protein-coding genes and non-coding RNAs. A few recent studies suggest that circRNAs are involved in teleost myogenesis, but the molecular networks involved remain poorly understood. In this study, an integrative omics approach was used to determine myogenic circRNAs in Nile tilapia by quantifying and comparing the expression profile of mRNAs, miRNAs, and circRNAs in fast muscle from full-sib fish with distinct growth rates. There were 1947 mRNAs, 9 miRNAs, and 4 circRNAs differentially expressed between fast- and slow-growing individuals. These miRNAs can regulate myogenic genes and have binding sites for the novel circRNA circMef2c. Our data indicate that circMef2c may interact with three miRNAs and 65 differentially expressed mRNAs to form multiple competing endogenous RNA networks that regulate growth, thus providing novel insights into the role of circRNAs in the regulation of muscle growth in teleosts.