Stefania Bortoluzzi

Identification of microRNA expression patterns and definition of a microRNA/mRNA regulatory network in distinct molecular groups of multiple myeloma

To date, little evidence of miRNA expression/deregulation in multiple myeloma has been reported. To characterize miRNA in the context of the major multiple myeloma molecular types, we generated miRNA expression profiles of highly purified malignant plasma cells from 40 primary tumors. Furthermore, transcriptional profiles, available for all patients, were used to investigate the occurrence of miRNA/predicted target mRNA pair anticorrelations, and the miRNA and genome-wide DNA data were integrated in a subset of patients to evaluate the influence of allelic imbalances on miRNA expression. Differential miRNA expression patterns were identified, which were mainly associated with the major IGH translocations; particularly, t(4;14) patients showed specific overexpression of let-7e, miR-125a-5p, and miR-99b belonging to a cluster at 19q13.33. The occurrence of other lesions (ie, 1q gain, 13q and 17p deletions, and hyperdiploidy) was slightly characterized by specific miRNA signatures. Furthermore, the occurrence of several allelic imbalances or loss of heterozygosity was found significantly associated with the altered expression of miRNAs located in the involved regions, such as let-7b at 22q13.31 or miR-140-3p at 16q22. Finally, the integrative analysis based on computational target prediction and miRNA/mRNA profiling defined a network of putative functional miRNA-target regulatory relations supported by expression data.

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.

Novel definition files for human GeneChips based on GeneAnnot

BackgroundImprovements in genome sequence annotation revealed discrepancies in the original probeset/gene assignment in Affymetrix microarray and the existence of differences between annotations and effective alignments of probes and transcription products. In the current generation of Affymetrix human GeneChips, most probesets include probes matching transcripts from more than one gene and probes which do not match any transcribed sequence.ResultsWe developed a novel set of custom Chip Definition Files (CDF) and the corresponding Bioconductor libraries for Affymetrix human GeneChips, based on the information contained in the GeneAnnot database. GeneAnnot-based CDFs are composed of unique custom-probesets, including only probes matching a single gene.ConclusionGeneAnnot-based custom CDFs solve the problem of a reliable reconstruction of expression levels and eliminate the existence of more than one probeset per gene, which often leads to discordant expression signals for the same transcript when gene differential expression is the focus of the analysis. GeneAnnot CDFs are freely distributed and fully compliant with Affymetrix standards and all available software for gene expression analysis. The CDF libraries are available from http://www.xlab.unimo.it/GA_CDF, along with supplementary information (CDF libraries, installation guidelines and R code, CDF statistics, and analysis results).

Transcriptome sequencing and microarray development for the Manila clam, Ruditapes philippinarum: genomic tools for environmental monitoring

BackgroundThe Manila clam, Ruditapes philippinarum, is one of the major aquaculture species in the world and a potential sentinel organism for monitoring the status of marine ecosystems. However, genomic resources for R. philippinarum are still extremely limited. Global analysis of gene expression profiles is increasingly used to evaluate the biological effects of various environmental stressors on aquatic animals under either artificial conditions or in the wild. Here, we report on the development of a transcriptomic platform for global gene expression profiling in the Manila clam.ResultsA normalized cDNA library representing a mixture of adult tissues was sequenced using a ultra high-throughput sequencing technology (Roche 454). A database consisting of 32,606 unique transcripts was constructed, 9,747 (30%) of which could be annotated by similarity. An oligo-DNA microarray platform was designed and applied to profile gene expression of digestive gland and gills. Functional annotation of differentially expressed genes between different tissues was performed by enrichment analysis. Expression of Natural Antisense Transcripts (NAT) analysis was also performed and bi-directional transcription appears a common phenomenon in the R. philippinarum transcriptome. A preliminary study on clam samples collected in a highly polluted area of the Venice Lagoon demonstrated the applicability of genomic tools to environmental monitoring.ConclusionsThe transcriptomic platform developed for the Manila clam confirmed the high level of reproducibility of current microarray technology. Next-generation sequencing provided a good representation of the clam transcriptome. Despite the known limitations in transcript annotation and sequence coverage for non model species, sufficient information was obtained to identify a large set of genes potentially involved in cellular response to environmental stress.

Differential expression of genes coding for ribosomal proteins in different human tissues

MOTIVATION To perform a computational and statistical study on a large set of gene expression data pertaining six adult human tissues (brain, liver, skeletal muscle, ovary, retina and uterus) for analyzing the expression of ribosomal protein genes. RESULTS Unexpectedly, in each of the considered tissues large variations in the expression of ribosomal protein genes were observed. Moreover, when comparing the expression levels of 89 ribosomal protein genes in six different tissues, 13 genes appeared differentially expressed among tissues. AVAILABILITY The expression data of the ribosomal protein genes together with supplementary material (complete transcriptional profiles of the considered human tissues) are freely available at the site GETProfiles (http://telethon.bio.unipd.it/GETProfiles/). CONTACT danieli@bio.unipd.it

MAGIA2: from miRNA and genes expression data integrative analysis to microRNA–transcription factor mixed regulatory circuits (2012 update)

MAGIA2 (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. MAGIA2 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 MAGIA2 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.

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

BackgroundQualitative 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.ResultsMost 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.ConclusionsUsing 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.

Sequencing, de novo annotation and analysis of the first Anguilla anguilla transcriptome: EeelBase opens new perspectives for the study of the critically endangered european eel

BackgroundOnce highly abundant, the European eel (Anguilla anguilla L.; Anguillidae; Teleostei) is considered to be critically endangered and on the verge of extinction, as the stock has declined by 90-99% since the 1980s. Yet, the species is poorly characterized at molecular level with little sequence information available in public databases.ResultsThe first European eel transcriptome was obtained by 454 FLX Titanium sequencing of a normalized cDNA library, produced from a pool of 18 glass eels (juveniles) from the French Atlantic coast and two sites in the Mediterranean coast. Over 310,000 reads were assembled in a total of 19,631 transcribed contigs, with an average length of 531 nucleotides. Overall 36% of the contigs were annotated to known protein/nucleotide sequences and 35 putative miRNA identified.ConclusionsThis study represents the first transcriptome analysis for a critically endangered species. EeelBase, a dedicated database of annotated transcriptome sequences of the European eel is freely available at http://compgen.bio.unipd.it/eeelbase. Considering the multiple factors potentially involved in the decline of the European eel, including anthropogenic factors such as pollution and human-introduced diseases, our results will provide a rich source of data to discover and identify new genes, characterize gene expression, as well as for identification of genetic markers scattered across the genome to be used in various applications.

Computational reconstruction of the human skeletal muscle secretome

In multicellular organisms, secreted proteins play pivotal regulatory roles in intercellular communication. Proteins secreted by skeletal muscle can act locally on muscle cells through autocrine/paracrine loops and on surrounding tissues such as muscle blood vessels, or they can be released into the blood stream, thus producing systemic effects. By a computational approach, we have screened 6255 products of genes expressed in normal human skeletal muscle. Putatively secreted proteins were identified by sequential steps of sieving, through prediction of signal peptide, recognition of transmembrane regions, and analysis of protein annotation. The resulting putative skeletal muscle secretome consists of 319 proteins, including 78 still uncharacterized proteins. This is the first human skeletal muscle secretome produced by computational analysis. Knowledge of proteins secreted by skeletal muscle could stimulate development of novel treatments for different diseases, including muscle atrophy and dystrophy. In addition, better knowledge of the secretion process in skeletal muscle can be useful for future gene therapy approaches. Proteins 2006.

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.