MiR&moRe2: A Bioinformatics Tool to Characterize microRNAs and microRNA-Offset RNAs from Small RNA-Seq Data.

MicroRNA-offset RNAs (moRNAs) are microRNA-like small RNAs generated by microRNA precursors. To date, little is known about moRNAs and bioinformatics tools to inspect their expression are still missing. We developed miR&moRe2, the first bioinformatics method to consistently characterize microRNAs, moRNAs, and their isoforms from small RNA sequencing data. To illustrate miR&moRe2 discovery power, we applied it to several published datasets. MoRNAs identified by miR&moRe2 were in agreement with previous research findings. Moreover, we observed that moRNAs and new microRNAs predicted by miR&moRe2 were downregulated upon the silencing of the microRNA-biogenesis pathway. Further, in a sizeable dataset of human blood cell populations, tens of novel miRNAs and moRNAs were discovered, some of them with significantly varied expression levels among the cell types. Results demonstrate that miR&moRe2 is a valid tool for a comprehensive study of small RNAs generated from microRNA precursors and could help to investigate their biogenesis and function.

Expanding the repertoire of miRNAs and miRNA-offset RNAs expressed in multiple myeloma by small RNA deep sequencing.

Microarray analysis of the multiple myeloma (MM) miRNome has unraveled the differential expression of miRNAs in cytogenetic subgroups, their involvement in the tumor biology and their effectiveness in prognostic models. Herein, the small RNA transcriptional landscape in MM has been investigated exploiting the possibilities offered by small RNA-seq, including accurate quantification of known mature species, discovery and characterization of isomiRs, and miRNA-offset RNAs (moRNAs). Matched small RNA-seq and miRNA GeneChip® microarray expression profiles were obtained in a representative panel of 30 primary MM tumors, fully characterized for genomic aberrations and mutations. RNA-seq and microarray gave concordant estimations of known species. Enhanced analysis of RNA-seq data with the miR&moRe pipeline led to the characterization of 655 known and 17 new mature miRNAs and of 74 moRNAs expressed in the considered cohort, 5 of which (moR-150-3p, moR-24-2-5p, moR-421-5p, moR-21-5p, and moR-6724-5p) at high level. Ectopic expression of miR-135a-3p in t(4;14) patients, upregulation of moR-150-3p and moR-21-5p in t(14;16)/t(14;20) samples, and of moR-6724-1-5p in patients overexpressing CCND1 were uncovered and validated by qRT-PCR. Overall, RNA-seq offered a more complete overview of small non-coding RNA in MM tumors, indicating specific moRNAs that demand further investigations to explore their role in MM biology.

Disentangling the microRNA regulatory milieu in multiple myeloma: integrative genomics analysis outlines mixed miRNA-TF circuits and pathway-derived networks modulated in t(4;14) patients.

The identification of overexpressed miRNAs in multiple myeloma (MM) has progressively added a further level of complexity to MM biology. miRNA and gene expression profiles of two large representative MM datasets, available from retrospective and prospective series and encompassing a total of 249 patients at diagnosis, were analyzed by means of in silico integrative genomics methods, based on MAGIA2 and Micrographite computational procedures. We first identified relevant miRNA/transcription factors/target gene regulation circuits in the disease and linked them to biological processes. Members of the miR-99b/let-7e/miR-125a cluster, or of its paralog, upregulated in t(4;14), were connected with the specific transcription factors PBX1 and CEBPA and several target genes. These results were validated in two additional independent plasma cell tumor datasets. Then, we reconstructed a non-redundant miRNA-gene regulatory network in MM, linking miRNAs, such as let-7g, miR-19a, mirR-20a, mir-21, miR-29 family, miR-34 family, miR-125b, miR-155, miR-221 to pathways associated with MM subtypes, in particular the ErbB, the Hippo, and the Acute myeloid leukemia associated pathways.

Small RNA Sequencing Uncovers New miRNAs and moRNAs Differentially Expressed in Normal and Primary Myelofibrosis CD34+ Cells.

Myeloproliferative neoplasms (MPN) are chronic myeloid cancers thought to arise at the level of CD34+ hematopoietic stem/progenitor cells. They include essential thrombocythemia (ET), polycythemia vera (PV) and primary myelofibrosis (PMF). All can progress to acute leukemia, but PMF carries the worst prognosis. Increasing evidences indicate that deregulation of microRNAs (miRNAs) might plays an important role in hematologic malignancies, including MPN. To attain deeper knowledge of short RNAs (sRNAs) expression pattern in CD34+ cells and of their possible role in mediating post-transcriptional regulation in PMF, we sequenced with Illumina HiSeq2000 technology CD34+ cells from healthy subjects and PMF patients. We detected the expression of 784 known miRNAs, with a prevalence of miRNA up-regulation in PMF samples, and discovered 34 new miRNAs and 99 new miRNA-offset RNAs (moRNAs), in CD34+ cells. Thirty-seven small RNAs were differentially expressed in PMF patients compared with healthy subjects, according to microRNA sequencing data. Five miRNAs (miR-10b-5p, miR-19b-3p, miR-29a-3p, miR-379-5p, and miR-543) were deregulated also in PMF granulocytes. Moreover, 3'-moR-128-2 resulted consistently downregulated in PMF according to RNA-seq and qRT-PCR data both in CD34+ cells and granulocytes. Target predictions of these validated small RNAs de-regulated in PMF and functional enrichment analyses highlighted many interesting pathways involved in tumor development and progression, such as signaling by FGFR and DAP12 and Oncogene Induced Senescence. As a whole, data obtained in this study deepened the knowledge of miRNAs and moRNAs altered expression in PMF CD34+ cells and allowed to identify and validate a specific small RNA profile that distinguishes PMF granulocytes from those of normal subjects. We thus provided new information regarding the possible role of miRNAs and, specifically, of new moRNAs in this disease.

A constitutive active MAPK/ERK pathway due to BRAFV600E positively regulates AHR pathway in PTC.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor mediating the toxicity and tumor-promoting properties of dioxin. AHR has been reported to be overexpressed and constitutively active in a variety of solid tumors, but few data are currently available concerning its role in thyroid cancer. In this study we quantitatively explored a series of 51 paired-normal and papillary thyroid carcinoma (PTC) tissues for AHR-related genes. We identified an increased AHR expression/activity in PTC, independently from its nuclear dimerization partner and repressor but strictly related to a constitutive active MAPK/ERK pathway. The AHR up-regulation followed by an increased expression of AHR target genes was confirmed by a meta-analysis of published microarray data, suggesting a ligand-independent active AHR pathway in PTC. In-vitro studies using a PTC-derived cell line (BCPAP) and HEK293 cells showed that BRAFV600E may directly modulate AHR localization, induce AHR expression and activity in an exogenous ligand-independent manner. The AHR pathway might represent a potential novel therapeutic target for PTC in the clinical practice.

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 (K562shFHC) comparing it with K562 transduced with scrambled RNA (K562shRNA). 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.

miRNA-mRNA integrative analysis in primary myelofibrosis CD34+ cells: role of miR-155/JARID2 axis in abnormal megakaryopoiesis.

Primary myelofibrosis (PMF) is a myeloproliferative neoplasm characterized by megakaryocyte (MK) hyperplasia, bone marrow fibrosis, and abnormal stem cell trafficking. PMF may be associated with somatic mutations in JAK2, MPL, or CALR. Previous studies have shown that abnormal MKs play a central role in the pathophysiology of PMF. In this work, we studied both gene and microRNA (miRNA) expression profiles in CD34(+) cells from PMF patients. We identified several biomarkers and putative molecular targets such as FGR, LCN2, and OLFM4. By means of miRNA-gene expression integrative analysis, we found different regulatory networks involved in the dysregulation of transcriptional control and chromatin remodeling. In particular, we identified a network gathering several miRNAs with oncogenic potential (eg, miR-155-5p) and targeted genes whose abnormal function has been previously associated with myeloid neoplasms, including JARID2, NR4A3, CDC42, and HMGB3. Because the validation of miRNA-target interactions unveiled JARID2/miR-155-5p as the strongest relationship in the network, we studied the function of this axis in normal and PMF CD34(+) cells. We showed that JARID2 downregulation mediated by miR-155-5p overexpression leads to increased in vitro formation of CD41(+) MK precursors. These findings suggest that overexpression of miR-155-5p and the resulting downregulation of JARID2 may contribute to MK hyperplasia in PMF.

miRNome of Italian Large White pig subcutaneous fat tissue: new miRNAs, isomiRs and moRNAs.

Small RNAs, such as micro-RNAs (miRNAs), are decisive regulators of gene expression, and they could determine adipose tissue traits. A better knowledge of porcine fat genomics is relevant given that the pig is a biomedical model for metabolic and cardiovascular human pathologies. Adipose tissue is particularly important for the meat industry. We explored the miRNome of two adult Italian Large White pig backfat samples by Illumina RNA-Seq. Using custom bioinformatic methods, the expressed miRNAs were identified and quantified and the nucleotide sequence variability of miRNA isoforms were analysed. We detected 222 known miRNAs, 68 new miRNAs and 17 miRNA-offset RNAs (moRNAs) expressed from known hairpins, and 312 new miRNAs expressed from 253 new hairpins. Porcine transcripts targeted by the most expressed miRNAs were predicted, showing that these miRNAs may have an impact on Wnt, insulin signalling and axon guidance pathways. The expression of five small RNAs, including moRNA ssc-5'-moR-21 and a miRNA from a new hairpin, was validated by a qRT-PCR assay, thus confirming the robustness of our results. The depicted miRNome complexity suggests that quantitative and qualitative features of miRNAs and non-canonical products of their precursors are worthy of further investigation to clarify their roles in the adipose tissue biology.

An integrative framework identifies alternative splicing events in colorectal cancer development.

Alternative splicing (AS) is a common mechanism which creates diverse RNA isoforms from a single gene, potentially increasing protein variety. Growing evidence suggests that this mechanism is closely related to cancer progression. In this study, whole transcriptome analysis was performed with GeneChip Human exon 1.0 ST Array from 80 samples comprising 23 normal colon mucosa, 30 primary colorectal cancer and 27 liver metastatic specimens from 46 patients, to identify AS events in colorectal cancer progression. Differentially expressed genes and exons were estimated and AS events were reconstructed by combining exon-level analyses with AltAnalyze algorithms and transcript-level estimations (MMBGX probabilistic method). The number of AS genes in the transition from normal colon mucosa to primary tumor was the most abundant, but fell considerably in the next transition to liver metastasis. 206 genes with probable AS events in colon cancer development and progression were identified, that are involved in processes and pathways relevant to tumor biology, as cell-cell and cell-matrix interactions. Several AS events in VCL, CALD1, B3GNT6 and CTHRC1 genes, differentially expressed during tumor development were validated, at RNA and at protein level. Taken together, these results demonstrate that cancer-specific AS is common in early phases of colorectal cancer natural history.

Human miRNome profiling in colorectal cancer and liver metastasis development.

Qualitative alterations or abnormal expression of microRNAs (miRNAs) in colorectal cancer has mainly been demonstrated in primary tumors. The miRNA expression profiles in 78 samples from 46 patients were analyzed to identify changes in miRNA expression level among normal colon mucosa, primary tumor and liver metastasis samples. Using this dataset, we describe the interplay of miRNA groups in regulating pathways that are important for tumor development. Here we describe in details the contents and quality controls for the miRNA expression and clinical data associated with the study published by Pizzini and colleagues in the BMC Genomics in 2013 (Pizzini et al., 2013). Data are deposited in GEO database as GSE35834 series.

Impact of microRNAs on regulatory networks and pathways in human colorectal carcinogenesis and development of metastasis.

BACKGROUND: Qualitative alterations or abnormal expression of microRNAs (miRNAs) in colon cancer have mainly been demonstrated in primary tumors. Poorly overlapping sets of oncomiRs, tumor suppressor miRNAs and metastamiRs have been linked with distinct stages in the progression of colorectal cancer. To identify changes in both miRNA and gene expression levels among normal colon mucosa, primary tumor and liver metastasis samples, and to classify miRNAs into functional networks, in this work miRNA and gene expression profiles in 158 samples from 46 patients were analysed. RESULTS: Most changes in miRNA and gene expression levels had already manifested in the primary tumors while these levels were almost stably maintained in the subsequent primary tumor-to-metastasis transition. In addition, comparing normal tissue, tumor and metastasis, we did not observe general impairment or any rise in miRNA biogenesis. While only few mRNAs were found to be differentially expressed between primary colorectal carcinoma and liver metastases, miRNA expression profiles can classify primary tumors and metastases well, including differential expression of miR-10b, miR-210 and miR-708. Of 82 miRNAs that were modulated during tumor progression, 22 were involved in EMT. qRT-PCR confirmed the down-regulation of miR-150 and miR-10b in both primary tumor and metastasis compared to normal mucosa and of miR-146a in metastases compared to primary tumor. The upregulation of miR-201 in metastasis compared both with normal and primary tumour was also confirmed. A preliminary survival analysis considering differentially expressed miRNAs suggested a possible link between miR-10b expression in metastasis and patient survival. By integrating miRNA and target gene expression data, we identified a combination of interconnected miRNAs, which are organized into sub-networks, including several regulatory relationships with differentially expressed genes. Key regulatory interactions were validated experimentally. Specific mixed circuits involving miRNAs and transcription factors were identified and deserve further investigation. The suppressor activity of miR-182 on ENTPD5 gene was identified for the first time and confirmed in an independent set of samples. CONCLUSIONS: Using a large dataset of CRC miRNA and gene expression profiles, we describe the interplay of miRNA groups in regulating gene expression, which in turn affects modulated pathways that are important for tumor development.

miR-142-3p prevents macrophage differentiation during cancer-induced myelopoiesis.

Tumor progression is accompanied by an altered myelopoiesis causing the accumulation of immunosuppressive cells. Here, we showed that miR-142-3p downregulation promoted macrophage differentiation and determined the acquisition of their immunosuppressive function in tumor. Tumor-released cytokines signaling through gp130, the common subunit of the interleukin-6 cytokine receptor family, induced the LAP∗ isoform of C/EBPß transcription factor, promoting macrophage generation. miR-142-3p downregulated gp130 by canonical binding to its messenger RNA (mRNA) 3' UTR and repressed C/EBPß LAP∗ by noncanonical binding to its 5' mRNA coding sequence. Enforced miR expression impaired macrophage differentiation both in vitro and in vivo. Mice constitutively expressing miR-142-3p in the bone marrow showed a marked increase in survival following immunotherapy with tumor-specific T lymphocytes. By modulating a specific miR in bone marrow precursors, we thus demonstrated the feasibility of altering tumor-induced macrophage differentiation as a potent tool to improve the efficacy of cancer immunotherapy.

MAGIA²: from miRNA and genes expression data integrative analysis to microRNA-transcription factor mixed regulatory circuits (2012 update).

MAGIA(2) (http://gencomp.bio.unipd.it/magia2) is an update, extension and evolution of the MAGIA web tool. It is dedicated to the integrated analysis of in silico target prediction, microRNA (miRNA) and gene expression data for the reconstruction of post-transcriptional regulatory networks. miRNAs are fundamental post-transcriptional regulators of several key biological and pathological processes. As miRNAs act prevalently through target degradation, their expression profiles are expected to be inversely correlated to those of the target genes. Low specificity of target prediction algorithms makes integration approaches an interesting solution for target prediction refinement. MAGIA(2) performs this integrative approach supporting different association measures, multiple organisms and almost all target predictions algorithms. Nevertheless, miRNAs activity should be viewed as part of a more complex scenario where regulatory elements and their interactors generate a highly connected network and where gene expression profiles are the result of different levels of regulation. The updated MAGIA(2) tries to dissect this complexity by reconstructing mixed regulatory circuits involving either miRNA or transcription factor (TF) as regulators. Two types of circuits are identified: (i) a TF that regulates both a miRNA and its target and (ii) a miRNA that regulates both a TF and its target.

Improving biomarker list stability by integration of biological knowledge in the learning process.

BACKGROUND: The identification of robust lists of molecular biomarkers related to a disease is a fundamental step for early diagnosis and treatment. However, methodologies for biomarker discovery using microarray data often provide results with limited overlap. It has been suggested that one reason for these inconsistencies may be that in complex diseases, such as cancer, multiple genes belonging to one or more physiological pathways are associated with the outcomes. Thus, a possible approach to improve list stability is to integrate biological information from genomic databases in the learning process; however, a comprehensive assessment based on different types of biological information is still lacking in the literature. In this work we have compared the effect of using different biological information in the learning process like functional annotations, protein-protein interactions and expression correlation among genes. RESULTS: Biological knowledge has been codified by means of gene similarity matrices and expression data linearly transformed in such a way that the more similar two features are, the more closely they are mapped. Two semantic similarity matrices, based on Biological Process and Molecular Function Gene Ontology annotation, and geodesic distance applied on protein-protein interaction networks, are the best performers in improving list stability maintaining almost equal prediction accuracy. CONCLUSIONS: The performed analysis supports the idea that when some features are strongly correlated to each other, for example because are close in the protein-protein interaction network, then they might have similar importance and are equally relevant for the task at hand. Obtained results can be a starting point for additional experiments on combining similarity matrices in order to obtain even more stable lists of biomarkers. The implementation of the classification algorithm is available at the link: http://www.math.unipd.it/~dasan/biomarkers.html.

Characterization and discovery of novel miRNAs and moRNAs in JAK2V617F-mutated SET2 cells.

To gain insights into a possible role of microRNAs in myeloproliferative neoplasms, we performed short RNA massive sequencing and extensive bioinformatic analysis in the JAK2V617F-mutated SET2 cell line. Overall, 652 known mature miRNAs were detected, of which 21 were highly expressed, thus being responsible of most of miRNA-mediated gene repression. microRNA putative targets were enriched in specific signaling pathways, providing information about cell activities under massive posttranscriptional regulation. The majority of miRNAs were mixtures of sequence variants, called isomiRs, mainly because of alternative, noncanonical processing of hairpin precursors. We also identified 78 novel miRNAs (miRNA*) derived from known hairpin precursors. Both major and minor (*) forms of miRNAs were expressed concurrently from half of expressed hairpins, highlighting the relevance of miRNA* and the complexity of strand selection bias regulation. Finally, we discovered that SET2 cells express a number of miRNA-offset RNAs (moRNAs), short RNAs derived from genomic regions flanking mature miRNAs. We provide novel data about the possible origin of moRNAs, although their functional role remains to be elucidated. Overall, this study shed light on the complexity of microRNA-mediated gene regulation in SET2 cells and represents the basis for future studies in JAK2V617F-mutated cellular models.

Impact of host genes and strand selection on miRNA and miRNA* expression.

Dysregulation of miRNAs expression plays a critical role in the pathogenesis of genetic, multifactorial disorders and in human cancers. We exploited sequence, genomic and expression information to investigate two main aspects of post-transcriptional regulation in miRNA biogenesis, namely strand selection regulation and expression relationships between intragenic miRNAs and host genes. We considered miRNAs expression profiles, measured in five sizeable microarray datasets, including samples from different normal cell types and tissues, as well as different tumours and disease states. First, the study of expression profiles of "sister" miRNA pairs (miRNA/miRNA*, 5' and 3' strands of the same hairpin precursor) showed that the strand selection is highly regulated since it shows tissue-/cell-/condition-specific modulation. We used information about the direction and the strength of the strand selection bias to perform an unsupervised cluster analysis for the sample classification evidencing that is able to distinguish among different tissues, and sometimes between normal and malignant cells. Then, considering a minimum expression threshold, in few miRNA pairs only one mature miRNA is always present in all considered cell types, whereas the majority of pairs were concurrently expressed in some cell types and alternatively in others. In a significant fraction of concurrently expressed pairs, the major and the minor forms found at comparable levels may contribute to post-transcriptional gene silencing, possibly in a coordinate way. In the second part of the study, the behaved tendency to co-expression of intragenic miRNAs and their "host" mRNA genes was confuted by expression profiles examination, suggesting that the expression profile of a given host gene can hardly be a good estimator of co-transcribed miRNA(s) for post-transcriptional regulatory networks inference. Our results point out the regulatory importance of post-transcriptional phases of miRNAs biogenesis, reinforcing the role of such layer of miRNA biogenesis in miRNA-based regulation of cell activities.

MAGIA, a web-based tool for miRNA and Genes Integrated Analysis.

MAGIA (miRNA and genes integrated analysis) is a novel web tool for the integrative analysis of target predictions, miRNA and gene expression data. MAGIA is divided into two parts: the query section allows the user to retrieve and browse updated miRNA target predictions computed with a number of different algorithms (PITA, miRanda and Target Scan) and Boolean combinations thereof. The analysis section comprises a multistep procedure for (i) direct integration through different functional measures (parametric and non-parametric correlation indexes, a variational Bayesian model, mutual information and a meta-analysis approach based on P-value combination) of mRNA and miRNA expression data, (ii) construction of bipartite regulatory network of the best miRNA and mRNA putative interactions and (iii) retrieval of information available in several public databases of genes, miRNAs and diseases and via scientific literature text-mining. MAGIA is freely available for Academic users at http://gencomp.bio.unipd.it/magia.

A-MADMAN: annotation-based microarray data meta-analysis tool.

BACKGROUND: Publicly available datasets of microarray gene expression signals represent an unprecedented opportunity for extracting genomic relevant information and validating biological hypotheses. However, the exploitation of this exceptionally rich mine of information is still hampered by the lack of appropriate computational tools, able to overcome the critical issues raised by meta-analysis. RESULTS: This work presents A-MADMAN, an open source web application which allows the retrieval, annotation, organization and meta-analysis of gene expression datasets obtained from Gene Expression Omnibus. A-MADMAN addresses and resolves several open issues in the meta-analysis of gene expression data. CONCLUSION: A-MADMAN allows i) the batch retrieval from Gene Expression Omnibus and the local organization of raw data files and of any related meta-information, ii) the re-annotation of samples to fix incomplete, or otherwise inadequate, metadata and to create user-defined batches of data, iii) the integrative analysis of data obtained from different Affymetrix platforms through custom chip definition files and meta-normalization. Software and documentation are available on-line at http://compgen.bio.unipd.it/bioinfo/amadman/.

Motif discovery in promoters of genes co-localized and co-expressed during myeloid cells differentiation.

Genes co-expressed may be under similar promoter-based and/or position-based regulation. Although data on expression, position and function of human genes are available, their true integration still represents a challenge for computational biology, hampering the identification of regulatory mechanisms. We carried out an integrative analysis of genomic position, functional annotation and promoters of genes expressed in myeloid cells. Promoter analysis was conducted by a novel multi-step method for discovering putative regulatory elements, i.e. over-represented motifs, in a selected set of promoters, as compared with a background model. The combination of transcriptional, structural and functional data allowed the identification of sets of promoters pertaining to groups of genes co-expressed and co-localized in regions of the human genome. The application of motif discovery to 26 groups of genes co-expressed in myeloid cells differentiation and co-localized in the genome showed that there are more over-represented motifs in promoters of co-expressed and co-localized genes than in promoters of simply co-expressed genes (CEG). Motifs, which are similar to the binding sequences of known transcription factors, non-uniformly distributed along promoter sequences and/or occurring in highly co-expressed subset of genes were identified. Co-expressed and co-localized gene sets were grouped in two co-expressed genomic meta-regions, putatively representing functional domains of a high-level expression regulation.