Genome-Wide Analyses of Radioresistance-Associated miRNA Expression Profile in Nasopharyngeal Carcinoma Using Next Generation Deep Sequencing

Background

Rapidly growing evidence suggests that microRNAs (miRNAs) are involved in a wide range of cancer malignant behaviours including radioresistance. Therefore, the present study was designed to investigate miRNA expression patterns associated with radioresistance in NPC.

Methods

The differential expression profiles of miRNAs and mRNAs associated with NPC radioresistance were constructed. The predicted target mRNAs of miRNAs and their enriched signaling pathways were analyzed via biological informatical algorithms. Finally, partial miRNAs and pathways-correlated target mRNAs were validated in two NPC radioreisitant cell models.

Results

50 known and 9 novel miRNAs with significant difference were identified, and their target mRNAs were narrowed down to 53 nasopharyngeal-/NPC-specific mRNAs. Subsequent KEGG analyses demonstrated that the 53 mRNAs were enriched in 37 signaling pathways. Further qRT-PCR assays confirmed 3 down-regulated miRNAs (miR-324-3p, miR-93-3p and miR-4501), 3 up-regulated miRNAs (miR-371a-5p, miR-34c-5p and miR-1323) and 2 novel miRNAs. Additionally, corresponding alterations of pathways-correlated target mRNAs were observed including 5 up-regulated mRNAs (ICAM1, WNT2B, MYC, HLA-F and TGF-β1) and 3 down-regulated mRNAs (CDH1, PTENP1 and HSP90AA1).

Conclusions

Our study provides an overview of miRNA expression profile and the interactions between miRNA and their target mRNAs, which will deepen our understanding of the important roles of miRNAs in NPC radioresistance.     

Inferring the perturbed microRNA regulatory networks from gene expression data using a network propagation based method

Background

MicroRNAs (miRNAs) are a class of endogenous small regulatory RNAs. Identifications of the dys-regulated or perturbed miRNAs and their key target genes are important for understanding the regulatory networks associated with the studied cellular processes. Several computational methods have been developed to infer the perturbed miRNA regulatory networks by integrating genome-wide gene expression data and sequence-based miRNA-target predictions. However, most of them only use the expression information of the miRNA direct targets, rarely considering the secondary effects of miRNA perturbation on the global gene regulatory networks.

Results

We proposed a network propagation based method to infer the perturbed miRNAs and their key target genes by integrating gene expressions and global gene regulatory network information. The method used random walk with restart in gene regulatory networks to model the network effects of the miRNA perturbation. Then, it evaluated the significance of the correlation between the network effects of the miRNA perturbation and the gene differential expression levels with a forward searching strategy. Results show that our method outperformed several compared methods in rediscovering the experimentally perturbed miRNAs in cancer cell lines. Then, we applied it on a gene expression dataset of colorectal cancer clinical patient samples and inferred the perturbed miRNA regulatory networks of colorectal cancer, including several known oncogenic or tumor-suppressive miRNAs, such as miR-17, miR-26 and miR-145.

Conclusions

Our network propagation based method takes advantage of the network effect of the miRNA perturbation on its target genes. It is a useful approach to infer the perturbed miRNAs and their key target genes associated with the studied biological processes using gene expression data.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-255) contains supplementary material, which is available to authorized users.     

Microarray-Based Approach Identifies Differentially Expressed MicroRNAs in Porcine Sexually Immature and Mature Testes

Background

MicroRNAs (miRNAs) are short non-coding RNA molecules which are proved to be involved in mammalian spermatogenesis. Their expression and function in the porcine germ cells are not fully understood.

Methodology

We employed a miRNA microarray containing 1260 unique miRNA probes to evaluate the miRNA expression patterns between sexually immature (60-day) and mature (180-day) pig testes. One hundred and twenty nine miRNAs representing 164 reporter miRNAs were expressed differently (p<0.1). Fifty one miRNAs were significantly up-regulated and 78 miRNAs were down-regulated in mature testes. Nine of these differentially expressed miRNAs were validated using quantitative RT-PCR assay. Totally 15919 putative miRNA-target sites were detected by using RNA22 method to align 445 NCBI pig cDNA sequences with these 129 differentially expressed miRNAs, and seven putative target genes involved in spermatogenesis including DAZL, RNF4 gene were simply confirmed by quantitative RT-PCR.

Conclusions

Overall, the results of this study indicated specific miRNAs expression in porcine testes and suggested that miRNAs had a role in regulating spermatogenesis.     

SNP Regulation of microRNA Expression and Subsequent Colon Cancer Risk

Introduction

MicroRNAs (miRNAs) regulate messenger RNAs (mRNAs) and as such have been implicated in a variety of diseases, including cancer. MiRNAs regulate mRNAs through binding of the miRNA 5’ seed sequence (~7–8 nucleotides) to the mRNA 3’ UTRs; polymorphisms in these regions have the potential to alter miRNA-mRNA target associations. SNPs in miRNA genes as well as miRNA-target genes have been proposed to influence cancer risk through altered miRNA expression levels.

Methods

MiRNA-SNPs and miRNA-target gene-SNPs were identified through the literature. We used SNPs from Genome-Wide Association Study (GWAS) data that were matched to individuals with miRNA expression data generated from an Agilent platform for colon tumor and non-tumor paired tissues. These samples were used to evaluate 327 miRNA-SNP pairs for associations between SNPs and miRNA expression levels as well as for SNP associations with colon cancer.

Results

Twenty-two miRNAs expressed in non-tumor tissue were significantly different by genotype and 21 SNPs were associated with altered tumor/non-tumor differential miRNA expression across genotypes. Two miRNAs were associated with SNP genotype for both non-tumor and tumor/non-tumor differential expression. Of the 41 miRNAs significantly associated with SNPs all but seven were significantly differentially expressed in colon tumor tissue. Two of the 41 SNPs significantly associated with miRNA expression levels were associated with colon cancer risk: rs8176318 (BRCA1), OR_AA 1.31 95% CI 1.01, 1.78, and rs8905 (PRKAR1A), OR_GG 2.31 95% CI 1.11, 4.77.

Conclusion

Of the 327 SNPs identified in the literature as being important because of their potential regulation of miRNA expression levels, 12.5% had statistically significantly associations with miRNA expression. However, only two of these SNPs were significantly associated with colon cancer.     

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.     

Cell-Free MicroRNA Expression Profiles in Malignant Effusion Associated with Patient Survival in Non-Small Cell Lung Cancer

Objective

MicroRNAs (miRNAs) expression is altered in cancer cells, and miRNAs could serve as diagnostic and prognostic biomarker for cancer patients. This study was designed to analyze circulating miRNAs expression in the malignant pleural effusion (MPE) and their association with patient survival in non-small cell lung cancer (NSCLC).

Methods

Pleural effusion from 184 patients with NSCLC and MPE were collected. MiRNA microarray and bioinformatics interpretation were used to evaluate miRNA expression profiles in 10 NSCLC patients with different survival prognosis. Associations were validated in 184 patients (randomly classified into training and validation set with equal number in each group) using quantitative RT-PCR. Risk scores were formulated based on the expression signature of miRNAs. Clinical data, such as patient survival, were collected for correlation analysis.

Results

Thirty-three miRNAs were found to be altered more than two-fold by microarray in malignant effusions between longer-survival and shorter-survival groups, and levels of five miRNAs (miRNA-93, miRNA-100, miRNA-134, miRNA-151 and miRNA-345) were significantly associated with overall survival. High expression of miR-100 and low expression of miRNA-93, miRNA-134, miRNA-151 and miRNA-345 were associated with poor survival in both the training and validation cohort. Patients with high risk scores had overall poor survival compared to the patients with low risk scores. Risk score was an independent predictor of patient survival.

Conclusions

Expression patterns of miRNAs are systematically altered in MPE of patient with NSCLC. The five miRNA signature from the effusion may serve as a predictor for the overall survival of patients with lung cancers.     

MiR-107 and MiR-185 Can Induce Cell Cycle Arrest in Human Non Small Cell Lung Cancer Cell Lines

Background

MicroRNAs (miRNAs) are short single stranded noncoding RNAs that suppress gene expression through either translational repression or degradation of target mRNAs. The annealing between messenger RNAs and 5′ seed region of miRNAs is believed to be essential for the specific suppression of target gene expression. One miRNA can have several hundred different targets in a cell. Rapidly accumulating evidence suggests that many miRNAs are involved in cell cycle regulation and consequentially play critical roles in carcinogenesis.

Methodology/Principal Findings

Introduction of synthetic miR-107 or miR-185 suppressed growth of the human non-small cell lung cancer cell lines. Flow cytometry analysis revealed these miRNAs induce a G1 cell cycle arrest in H1299 cells and the suppression of cell cycle progression is stronger than that by Let-7 miRNA. By the gene expression analyses with oligonucleotide microarrays, we find hundreds of genes are affected by transfection of these miRNAs. Using miRNA-target prediction analyses and the array data, we listed up a set of likely targets of miR-107 and miR-185 for G1 cell cycle arrest and validate a subset of them using real-time RT-PCR and immunoblotting for CDK6.

Conclusions/Significance

We identified new cell cycle regulating miRNAs, miR-107 and miR-185, localized in frequently altered chromosomal regions in human lung cancers. Especially for miR-107, a large number of down-regulated genes are annotated with the gene ontology term ‘cell cycle’. Our results suggest that these miRNAs may contribute to regulate cell cycle in human malignant tumors.     

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.     

Plant microRNA-Target Interaction Identification Model Based on the Integration of Prediction Tools and Support Vector Machine

Background

Confident identification of microRNA-target interactions is significant for studying the function of microRNA (miRNA). Although some computational miRNA target prediction methods have been proposed for plants, results of various methods tend to be inconsistent and usually lead to more false positive. To address these issues, we developed an integrated model for identifying plant miRNA–target interactions.

Results

Three online miRNA target prediction toolkits and machine learning algorithms were integrated to identify and analyze Arabidopsis thaliana miRNA-target interactions. Principle component analysis (PCA) feature extraction and self-training technology were introduced to improve the performance. Results showed that the proposed model outperformed the previously existing methods. The results were validated by using degradome sequencing supported Arabidopsis thaliana miRNA-target interactions. The proposed model constructed on Arabidopsis thaliana was run over Oryza sativa and Vitis vinifera to demonstrate that our model is effective for other plant species.

Conclusions

The integrated model of online predictors and local PCA-SVM classifier gained credible and high quality miRNA-target interactions. The supervised learning algorithm of PCA-SVM classifier was employed in plant miRNA target identification for the first time. Its performance can be substantially improved if more experimentally proved training samples are provided.     

miR-29b,miR-205 and miR-221 Enhance Chemosensitivity to Gemcitabine in HuH28 Human Cholangiocarcinoma Cells

Background and Aims

Cholangiocarcinoma (CCA) is highly resistant to chemotherapy, including gemcitabine (Gem) treatment. MicroRNAs (miRNAs) are endogenous, non-coding, short RNAs that can regulate multiple genes expression. Some miRNAs play important roles in the chemosensitivity of tumors. Here, we examined the relationship between miRNA expression and the sensitivity of CCA cells to Gem.

Methods

Microarray analysis was used to determine the miRNA expression profiles of two CCA cell lines, HuH28 and HuCCT1. To determine the effect of candidate miRNAs on Gem sensitivity, expression of each candidate miRNA was modified via either transfection of a miRNA mimic or transfection of an anti-oligonucleotide. Ontology-based programs were used to identify potential target genes of candidate miRNAs that were confirmed to affect the Gem sensitivity of CCA cells.

Results

HuCCT1 cells were more sensitive to Gem than were HuH28 cells, and 18 miRNAs were differentially expressed whose ratios over ± 2log2 between HuH28 and HuCCT1. Among these 18 miRNAs, ectopic overexpression of each of three downregulated miRNAs in HuH28 (miR-29b, miR-205, miR-221) restored Gem sensitivity to HuH28. Suppression of one upregulated miRNA in HuH28, miR-125a-5p, inhibited HuH28 cell proliferation independently to Gem treatment. Selective siRNA-mediated downregulation of either of two software-predicted targets, PIK3R1 (target of miR-29b and miR-221) or MMP-2 (target of miR-29b), also conferred Gem sensitivity to HuH28.

Conclusions

miRNA expression profiling was used to identify key miRNAs that regulate Gem sensitivity in CCA cells, and software that predicts miRNA targets was used to identify promising target genes for anti-tumor therapies.     

Age and sex differences in kidney microRNA expression during the life span of F344 rats

Background

Growing evidence suggests that epigenetic mechanisms of gene regulation may play a role in susceptibilities to specific toxicities and adverse drug reactions. MiRNAs in particular have been shown to be important regulators in cancer and other diseases and show promise as predictive biomarkers for diagnosis and prognosis. In this study, we characterized the global kidney miRNA expression profile in untreated male and female F344 rats throughout the life span. These findings were correlated with sex-specific susceptibilities to adverse renal events, such as male-biased renal fibrosis and inflammation in old age.

Methods

Kidney miRNA expression was examined in F344 rats at 2, 5, 6, 8, 15, 21, 78, and 104 weeks of age in both sexes using Agilent miRNA microarrays. Differential expression was determined using filtering criteria of ≥1.5 fold change and ANOVA or pairwise t-test (FDR <5%) to determine significant age and sex effects, respectively. Pathway analysis software was used to investigate the possible roles of these target genes in age- and sex-specific differences.

Results

Three hundred eleven miRNAs were found to be expressed in at least one age and sex. Filtering criteria revealed 174 differentially expressed miRNAs in the kidney; 173 and 34 miRNAs exhibiting age and sex effects, respectively. Principal component analysis revealed age effects predominated over sex effects, with 2-week miRNA expression being much different from other ages. No significant sexually dimorphic miRNA expression was observed from 5 to 8 weeks, while the most differential expression (13 miRNAs) was observed at 21 weeks. Potential target genes of these differentially expressed miRNAs were identified.

Conclusions

The expression of 56% of detected renal miRNAs was found to vary significantly with age and/or sex during the life span of F344 rats. Pathway analysis suggested that 2-week-expressed miRNAs may be related to organ and cellular development and proliferation pathways. Male-biased miRNA expression at older ages correlated with male-biased renal fibrosis and mononuclear cell infiltration. These miRNAs showed high representation in renal inflammation and nephritis pathways, and included miR-214, miR-130b, miR-150, miR-223, miR-142-5p, miR-185, and miR-296*. Analysis of kidney miRNA expression throughout the rat life span will improve the use of current and future renal biomarkers and inform our assessments of kidney injury and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-014-0019-1) contains supplementary material, which is available to authorized users.     

Comparison of Microarray Platforms for Measuring Differential MicroRNA Expression in Paired Normal/Cancer Colon Tissues

Background

Microarray technology applied to microRNA (miRNA) profiling is a promising tool in many research fields; nevertheless, independent studies characterizing the same pathology have often reported poorly overlapping results. miRNA analysis methods have only recently been systematically compared but only in few cases using clinical samples.

Methodology/Principal Findings

We investigated the inter-platform reproducibility of four miRNA microarray platforms (Agilent, Exiqon, Illumina, and Miltenyi), comparing nine paired tumor/normal colon tissues. The most concordant and selected discordant miRNAs were further studied by quantitative RT-PCR. Globally, a poor overlap among differentially expressed miRNAs identified by each platform was found. Nevertheless, for eight miRNAs high agreement in differential expression among the four platforms and comparability to qRT-PCR was observed. Furthermore, most of the miRNA sets identified by each platform are coherently enriched in data from the other platforms and the great majority of colon cancer associated miRNA sets derived from the literature were validated in our data, independently from the platform. Computational integration of miRNA and gene expression profiles suggested that anti-correlated predicted target genes of differentially expressed miRNAs are commonly enriched in cancer-related pathways and in genes involved in glycolysis and nutrient transport.

Conclusions

Technical and analytical challenges in measuring miRNAs still remain and further research is required in order to increase consistency between different microarray-based methodologies. However, a better inter-platform agreement was found by looking at miRNA sets instead of single miRNAs and through a miRNAs – gene expression integration approach.     

MicroRNA-5p and -3p co-expression and cross-targeting in colon cancer cells

Background

Two mature miRNA species may be generated from the 5’ and 3’ arms of a pre-miRNA precursor. In most cases, only one species remains while the complementary species is degraded. However, co-existence of miRNA-5p and -3p species is increasingly being reported. In this work, we aimed to systematically investigate co-expression of miRNA-5p/3p in colon cancer cells in a genome-wide analysis, and to examine cross-targeting of the dysregulated miRNAs and 5p/3p species.

Results

Four colon cancer cell lines were examined relative to two normal colon tissues. Of the 1,190 miRNAs analyzed, 92 and 36 were found to be up- or down-regulated, respectively, in cancer cells. Nineteen co-expressed miRNA-5p/3p pairs were further identified suggesting frequent 5p/3p co-accumulation in colon cancer cells. Of these, 14 pairs were co-up-regulated and 3 pairs were co-down-regulated indicating concerted 5p/3p dysregulation. Nine dysregulated miRNA pairs fell into three miRNA gene families, namely let-7, mir-8/200 and mir-17, which showed frequent cross-targeting in the metastasis process. Focusing on the let-7d-5p/3p pair, the respectively targeted IGF1R and KRAS were shown to be in a reverse relationship with expression of the respective miRNA, which was confirmed in transient transfection assays using let-7d mimic or inhibitor. Targeting of KRAS by let-7d was previous reported; targeting of IGF1R by let-7d-5p was confirmed in luciferase assays in this study. The findings of let-7d-5p/3p and multiple other miRNAs targeting IGF1R, KRAS and other metastasis-related factors suggest that 5p/3p miRNAs contribute to cross-targeting of multiple cancer-associated factors and processes possibly to evade functional abolishment when any one of the crucial factors are inactivated.

Conclusions

miRNA-5p/3p species are frequently co-expressed and are coordinately regulated in colon cancer cells. In cancer cells, multiple cross-targeting by the miRNAs, including the co-existing 5p/3p species, frequently occurs in an apparent safe-proof scheme of miRNA regulation of important tumorigenesis processes. Further systematic analysis of co-existing miRNA-5p/3p pairs in clinical tissues is important in elucidating 5p/3p contributions to cancer pathogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-014-0095-x) contains supplementary material, which is available to authorized users.     

MicroRNA expression profiling of the fifth-instar posterior silk gland of Bombyx mori

Background

The growth and development of the posterior silk gland and the biosynthesis of the silk core protein at the fifth larval instar stage of Bombyx mori are of paramount importance for silk production.

Results

Here, aided by next-generation sequencing and microarry assay, we profile 1,229 microRNAs (miRNAs), including 728 novel miRNAs and 110 miRNA/miRNA* duplexes, of the posterior silk gland at the fifth larval instar. Target gene prediction yields 14,222 unique target genes from 1,195 miRNAs. Functional categorization classifies the targets into complex pathways that include both cellular and metabolic processes, especially protein synthesis and processing.

Conclusion

The enrichment of target genes in the ribosome-related pathway indicates that miRNAs may directly regulate translation. Our findings pave a way for further functional elucidation of these miRNAs and their targets in silk production.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-410) contains supplementary material, which is available to authorized users.