Tissue and Process Specific microRNA–mRNA Co-Expression in Mammalian Development and Malignancy

An association between enrichment and depletion of microRNA (miRNA) binding sites, 3′ UTR length, and mRNA expression has been demonstrated in various developing tissues and tissues from different mature organs; but functional, context-dependent miRNA regulations have yet to be elucidated. Towards that goal, we examined miRNA–mRNA interactions by measuring miRNA and mRNA in the same tissue during development and also in malignant conditions. We identified significant miRNA-mediated biological process categories in developing mouse cerebellum and lung using non-targeted mRNA expression as the negative control. Although miRNAs in general suppress target mRNA messages, many predicted miRNA targets demonstrate a significantly higher level of co-expression than non-target genes in developing cerebellum. This phenomenon is tissue specific since it is not observed in developing lungs. Comparison of mouse cerebellar development and medulloblastoma demonstrates a shared miRNA–mRNA co-expression program for brain-specific neurologic processes such as synaptic transmission and exocytosis, in which miRNA target expression increases with the accumulation of multiple miRNAs in developing cerebellum and decreases with the loss of these miRNAs in brain tumors. These findings demonstrate the context-dependence of miRNA–mRNA co-expression.     

The mRNA–miRNA–lncRNA Regulatory Network and Factors Associated with Prognosis Prediction of Hepatocellular Carcinoma

The aim of this study was to identify novel prognostic mRNA and microRNA (miRNA) biomarkers for hepatocellular carcinoma (HCC) using methods in systems biology. Differentially expressed mRNAs, miRNAs, and long non-coding RNAs (lncRNAs) were compared between HCC tumor tissues and normal liver tissues in The Cancer Genome Atlas (TCGA) database. Subsequently, a prognosis-associated mRNA co-expression network, an mRNA–miRNA regulatory network, and an mRNA–miRNA–lncRNA regulatory network were constructed to identify prognostic biomarkers for HCC through Cox survival analysis. Seven prognosis-associated mRNA co-expression modules were obtained by analyzing these differentially expressed mRNAs. An expression module including 120 mRNAs was significantly correlated with HCC patient survival. Combined with patient survival data, several mRNAs and miRNAs, including CHST4, SLC22A8, STC2, hsa-miR-326, and hsa-miR-21 were identified from the network to predict HCC patient prognosis. Clinical significance was investigated using tissue microarray analysis of samples from 258 patients with HCC. Functional annotation of hsa-miR-326 and hsa-miR-21-5p indicated specific associations with several cancer-related pathways. The present study provides a bioinformatics method for biomarker screening, leading to the identification of an integrated mRNA–miRNA–lncRNA regulatory network and their co-expression patterns in relation to predicting HCC patient survival.     

Integrated miRNA and mRNA Expression Profiling in Inflamed Colon of Patients with Ulcerative Colitis

Background

Ulcerative colitis (UC) is associated with differential colonic expression of genes involved in immune response (e.g. IL8) and barrier integrity (e.g. cadherins). MicroRNAs (miRNAs) are regulators of gene expression and are involved in various immune-related diseases. In this study, we investigated (1) if miRNA expression in UC mucosa is altered and (2) if any of these changes correlate with mucosal mRNA expression. Integration of mRNA and miRNA expression profiling may allow the identification of functional links between dysregulated miRNAs and their target mRNA.

Methodology

Colonic mucosal biopsies were obtained from 17 UC (10 active and 7 inactive) patients and 10 normal controls. Total RNA was used to analyze miRNA and mRNA expression via Affymetrix miRNA 2.0 and Affymetrix Human Gene 1.0ST arrays, respectively. Both miRNA and gene expression profiles were integrated by correlation analysis to identify dysregulated miRNAs with their corresponding predicted target mRNA. Microarray data were validated with qRT-PCR. Regulation of IL8 and CDH11 expression by hsa-miR-200c-3p was determined by luciferase reporter assays.

Results

When comparing active UC patients vs. controls, 51 miRNAs and 1543 gene probe sets gave significantly different signals. In contrast, in inactive UC vs. controls, no significant miRNA expression differences were found while 155 gene probe sets had significantly different signals. We then identified potential target genes of the significantly dysregulated miRNAs and genes in active UC vs. controls and found a highly significant inverse correlation between hsa-miR-200c-3p and IL8, an inflammatory marker, and between hsa-miR-200c-3p and CDH11, a gene related to intestinal epithelial barrier function. We could demonstrate that hsa-miR-200c-3p directly regulates IL8 and CDH11 expression.

Conclusion

Differential expression of immune- and barrier-related genes in inflamed UC mucosa may be influenced by altered expression of miRNAs. Integrated analysis of miRNA and mRNA expression profiles revealed hsa-miR-200c-3p for use of miRNA mimics as therapeutics.     

miRNA在人骨髓来源间充质干细胞软骨诱导分化过程中的表达

目的：研究miRNAs在人骨髓来源间充质干细胞软骨诱导分化过程中的表达情况。方法：以从骨髓中分离培养的MSCs及软骨诱导培养后的细胞为实验对象,利用基因芯片检测miRNAs的表达情况,由SAM分析得到MSCs较其诱导培养细胞中差异表达的miRNAs,再进行生物信息学分析。结果：①分离培养出的MSCs经软骨诱导培养21天后,已具有软骨细胞特性,经芯片检测并SAM分析,软骨诱导培养的细胞较MSCs高表达的miRNAs有6个：hsa-miR-572、hsa-miR-130b、hsa-miR-193b、hsa-miR-28、hsa-miR-152、hsa-miR-560;软骨诱导培养的细胞较MSCs低表达的miRNAs有2个：hsa-miR-424、hsa-miR-122a。②利用TargetScan预测其靶基因,并行生物信息学分析,其中hsa-miR-130b、hsa-miR-193b、hsa-miR-152及hsa-miR-424的预测靶基因中多为参与细胞分化、骨形成、软骨形成及干细胞表型相关的基因。结论：hsa-miR-130b、hsa-miR-193b、hsa-miR-152和hsa-miR-424等对人骨髓来源间充质干细胞的软骨分化起着重要调控作用。     

Comparative microarray analysis of microRNA expression profiles in primary cutaneous malignant melanoma, cutaneous malignant melanoma metastases, and benign melanocytic nevi

Perturbations in microRNA (miRNA) expression profiles have been reported for cutaneous malignant melanoma (CMM) predominantly when examined in cell lines. Despite the rapidly growing number of newly discovered human miRNA sequences, the availability of up-to-date miRNA expression profiles for clinical samples of primary cutaneous malignant melanoma (PCMM), cutaneous malignant melanoma metastases (CMMM), and benign melanocytic nevi (BMN) is limited. Specimens excised from the center of tumors (lesional) from patients with PCMM (n=9), CMMM (n=4), or BMN (n=8) were obtained during surgery. An exploratory microarray analysis was performed by miRNA expression profiling based on Agilent platform screening for 1205 human miRNAs. The results from the microarray analysis were validated by TaqMan quantitative real-time polymerase chain reaction. In addition to several miRNAs previously known to be associated with CMM, 19 unidentified miRNA candidates were found to be dysregulated in CMM patient samples. Among the 19 novel miRNA candidates, the genes hsa-miR-22, hsa-miR-130b, hsa-miR-146b-5p, hsa-miR-223, hsa-miR-301a, hsa-miR-484, hsa-miR-663, hsa-miR-720, hsa-miR-1260, hsa-miR-1274a, hsa-miR-1274b, hsa-miR-3663-3p, hsa-miR-4281, and hsa-miR-4286 were upregulated, and the genes hsa-miR-24-1*, hsa-miR-26a, hsa-miR-4291, hsa-miR-4317, and hsa-miR-4324 were downregulated. The results of this study partially confirm previous CMM miRNA profiling studies identifying miRNAs that are dysregulated in CMM. However, we report several novel miRNA candidates in CMM tumors; these miRNA sequences require further validation and functional analysis to evaluate whether they play a role in the pathogenesis of CMM.     

GCH1-regulated miRNAs are potential targets for microglial activation in neuropathic pain

Neuropathic pain (NP) is a chronic pain directly caused by injury or disease of the somatosensory nervous system. Previous studies suggest that GTP cyclohydrolase I (GCH1) may play a pivotal role in microglial activation, which has been shown to be essential for NP. However, its underlying mechanisms in microglial activation remain unclear. A wide range of microRNAs (miRNAs) have been found to be involved in microglial activation-induced NP. To identify the miRNAs regulated by GCH1 and predict their functions in the progression of microglial activation, we analyzed the miRNA expression profiles of GCH1-knockdown (KD) BV2 microglial cells. Small RNA-sequencing analysis revealed 13 differentially expressed (DE) miRNAs in GCH1-KD cells. The target genes of DE miRNAs mainly participate in PI3K-Akt signaling pathway, peroxisome and ferroptosis. The miRNA–mRNA regulatory network analysis showed that GCH1, MAP4K5 and YWHAB acted as hub genes. qRT-PCR results further verified the expression levels of mmu-miR-1a-3p, mmu-miR-133a-3p, mmu-miR-7a-5p and mmu-miR-10a-5p in GCH1-KD cells, which were consistent with the sequencing data. In addition, our data indicated that overexpression of mmu-miR-133a-3p alleviated the pro-inflammatory cytokines IL-1β and IL-6 production induced by lipopolysaccharide (LPS), indicating that mmu-miR-133a-3p has a negative effect on microglial activation. Taken together, our findings suggest that many miRNAs regulated by GCH1 may be involved in microglial activation, which may provide new potential targets for GCH1 in the pathogenesis of NP.     

Systematic study of cis-antisense miRNAs in animal species reveals miR-3661 to target PPP2CA in human cells

MicroRNAs (miRNAs) suppress targeting gene expression through blocking translation or triggering mRNA degradation and, in general, act in trans, through a partially complementary interaction with the 3′ untranslated region (3′ UTR) or coding regions of a target gene. Although it has been reported previously that some miRNAs suppress their target genes on the opposite strand with a fully complementary sequence (i.e., natural antisense miRNAs that act in cis), there is no report to systematically study such cis-antisense miRNAs in different animal species. Here we report that cis-antisense miRNAs do exist in different animal species: 48 in Caenorhabditis elegans, 17 in Drosophila, 36 in Mus musculus, and 52 in Homo sapiens using a systematical bioinformatics approach. We show that most of these cis-antisense miRNAs can efficiently reduce the expression levels of their target genes in human cells. We further investigate hsa-miR-3661, one of the predicted cis-antisense miRNAs, in detail and demonstrate that this miRNA directly targets the coding sequence of PPP2CA located on the opposite DNA strand and inhibits the PPP2CA expression. Taken together, these results indicate that cis-antisense miRNAs are conservative and functional in animal species including humans.     

A Quest for miRNA Bio-Marker: A Track Back Approach from Gingivo Buccal Cancer to Two Different Types of Precancers

Deregulation of miRNA expression may contribute to tumorigenesis and other patho-physiology associated with cancer. Using TLDA, expression of 762 miRNAs was checked in 18 pairs of gingivo buccal cancer-adjacent control tissues. Expression of significantly deregulated miRNAs was further validated in cancer and examined in two types of precancer (leukoplakia and lichen planus) tissues by primer-specific TaqMan assays. Biological implications of these miRNAs were assessed bioinformatically. Expression of hsa-miR-1293, hsa-miR-31, hsa-miR-31* and hsa-miR-7 were significantly up-regulated and those of hsa-miR-206, hsa-miR-204 and hsa-miR-133a were significantly down-regulated in all cancer samples. Expression of only hsa-miR-31 was significantly up-regulated in leukoplakia but none in lichen planus samples. Analysis of expression heterogeneity divided 18 cancer samples into clusters of 13 and 5 samples and revealed that expression of 30 miRNAs (including the above-mentioned 7 miRNAs), was significantly deregulated in the cluster of 13 samples. From database mining and pathway analysis it was observed that these miRNAs can significantly target many of the genes present in different cancer related pathways such as “proteoglycans in cancer”, PI3K-AKT etc. which play important roles in expression of different molecular features of cancer. Expression of hsa-miR-31 was significantly up-regulated in both cancer and leukoplakia tissues and, thus, may be one of the molecular markers of leukoplakia which may progress to gingivo-buccal cancer.     

H-Ferritin-Regulated MicroRNAs Modulate Gene Expression in K562 Cells

In a previous study, we showed that the silencing of the heavy subunit (FHC) offerritin, the central iron storage molecule in the cell, is accompanied by a modification in global gene expression. In this work, we explored whether different FHC amounts might modulate miRNA expression levels in K562 cells and studied the impact of miRNAs in gene expression profile modifications. To this aim, we performed a miRNA-mRNA integrative analysis in K562 silenced for FHC (K562^shFHC) comparing it with K562 transduced with scrambled RNA (K562^shRNA). Four miRNAs, namely hsa-let-7g, hsa-let-7f, hsa-let-7i and hsa-miR-125b, were significantly up-regulated in silenced cells. The remarkable down-regulation of these miRNAs, following FHC expression rescue, supports a specific relation between FHC silencing and miRNA-modulation. The integration of target predictions with miRNA and gene expression profiles led to the identification of a regulatory network which includes the miRNAs up-regulated by FHC silencing, as well as91 down-regulated putative target genes. These genes were further classified in 9 networks; the highest scoring network, “Cell Death and Survival, Hematological System Development and Function, Hematopoiesis”, is composed by 18 focus molecules including RAF1 and ERK1/2. We confirmed that, following FHC silencing, ERK1/2 phosphorylation is severely impaired and that RAF1 mRNA is significantly down-regulated. Taken all together, our data indicate that, in our experimental model, FHC silencing may affect RAF1/pERK1/2 levels through the modulation of a specific set of miRNAs and add new insights in to the relationship among iron homeostasis and miRNAs.     

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.     

Functional Profiling of Precursor MicroRNAs Identifies MicroRNAs Essential for Glioma Proliferation

Cancer initiation and progression involve microRNAs that can function like tumor suppressors and oncogenes. The functional significance of most miRNAs is currently unknown. To determine systematically which microRNAs are essential for glioma growth, we screened a precursor microRNA library in three human glioblastoma and one astroglial cell line model systems. The most prominent and consistent cell proliferation–reducing hits were validated in secondary screening with an additional apoptosis endpoint. The functional screening data were integrated in the miRNA expression data to find underexpressed true functional tumor suppressor miRNAs. In addition, we used miRNA-target gene predictions and combined siRNA functional screening data to find the most probable miRNA-target gene pairs with a similar functional effect on proliferation. Nine novel functional miRNAs (hsa-miR-129, -136, -145, -155, -181b, -342-5p, -342-3p, -376a/b) in GBM cell lines were validated for their importance in glioma cell growth, and similar effects for six target genes (ROCK1, RHOA, MET, CSF1R, EIF2AK1, FGF7) of these miRNAs were shown functionally. The clinical significance of the functional hits was validated in miRNA expression data from the TCGA glioblastoma multiforme (GBM) tumor cohort. Five tumor suppressor miRNAs (hsa-miR-136, -145, -342, -129, -376a) showed significant underexpression in clinical GBM tumor samples from the TCGA GBM cohort further supporting the role of these miRNAs in vivo. Most importantly, higher hsa-miR-145 expression in GBM tumors yielded significantly better survival (p<0.005) in a subset of patients thus validating it as a genuine tumor suppressor miRNA. This systematic functional profiling provides important new knowledge about functionally relevant miRNAs in GBM biology and may offer new targets for treating glioma.