InteractomeSeq: a web server for the identification and profiling of domains and epitopes from phage display and next generation sequencing data.

High-Throughput Sequencing technologies are transforming many research fields, including the analysis of phage display libraries. The phage display technology coupled with deep sequencing was introduced more than a decade ago and holds the potential to circumvent the traditional laborious picking and testing of individual phage rescued clones. However, from a bioinformatics point of view, the analysis of this kind of data was always performed by adapting tools designed for other purposes, thus not considering the noise background typical of the 'interactome sequencing' approach and the heterogeneity of the data. InteractomeSeq is a web server allowing data analysis of protein domains ('domainome') or epitopes ('epitome') from either Eukaryotic or Prokaryotic genomic phage libraries generated and selected by following an Interactome sequencing approach. InteractomeSeq allows users to upload raw sequencing data and to obtain an accurate characterization of domainome/epitome profiles after setting the parameters required to tune the analysis. The release of this tool is relevant for the scientific and clinical community, because InteractomeSeq will fill an existing gap in the field of large-scale biomarkers profiling, reverse vaccinology, and structural/functional studies, thus contributing essential information for gene annotation or antigen identification. InteractomeSeq is freely available at https://InteractomeSeq.ba.itb.cnr.it/.

Combined microRNA and mRNA expression analysis in pediatric multiple sclerosis: an integrated approach to uncover novel pathogenic mechanisms of the disease.

Multiple sclerosis (MS) is a complex disease of the CNS that usually affects young adults, although 3-5% of cases are diagnosed in childhood and adolescence (hence called pediatric MS, PedMS). Genetic predisposition, among other factors, seems to contribute to the risk of the onset, in pediatric as in adult ages, but few studies have investigated the genetic 'environmentally naïve' load of PedMS. The main goal of this study was to identify circulating markers (miRNAs), target genes (mRNAs) and functional pathways associated with PedMS; we also verified the impact of miRNAs on clinical features, i.e. disability and cognitive performances. The investigation was performed in 19 PedMS and 20 pediatric controls (PCs) using a High-Throughput Next-generation Sequencing (HT-NGS) approach followed by an integrated bioinformatics/biostatistics analysis. Twelve miRNAs were significantly upregulated (let-7a-5p, let-7b-5p, miR-25-3p, miR-125a-5p, miR-942-5p, miR-221-3p, miR-652-3p, miR-182-5p, miR-185-5p, miR-181a-5p, miR-320a, miR-99b-5p) and 1 miRNA was downregulated (miR-148b-3p) in PedMS compared with PCs. The interactions between the significant miRNAs and their targets uncovered predicted genes (i.e. TNFSF13B, TLR2, BACH2, KLF4) related to immunological functions, as well as genes involved in autophagy-related processes (i.e. ATG16L1, SORT1, LAMP2) and ATPase activity (i.e. ABCA1, GPX3). No significant molecular profiles were associated with any PedMS demographic/clinical features. Both miRNAs and mRNA expressions predicted the phenotypes (PedMS-PC) with an accuracy of 92% and 91%, respectively. In our view, this original strategy of contemporary miRNA/mRNA analysis may help to shed light in the genetic background of the disease, suggesting further molecular investigations in novel pathogenic mechanisms.

ITSoneDB: a comprehensive collection of eukaryotic ribosomal RNA Internal Transcribed Spacer 1 (ITS1) sequences.

A holistic understanding of environmental communities is the new challenge of metagenomics. Accordingly, the amplicon-based or metabarcoding approach, largely applied to investigate bacterial microbiomes, is moving to the eukaryotic world too. Indeed, the analysis of metabarcoding data may provide a comprehensive assessment of both bacterial and eukaryotic composition in a variety of environments, including human body. In this respect, whereas hypervariable regions of the 16S rRNA are the de facto standard barcode for bacteria, the Internal Transcribed Spacer 1 (ITS1) of ribosomal RNA gene cluster has shown a high potential in discriminating eukaryotes at deep taxonomic levels. As metabarcoding data analysis rely on the availability of a well-curated barcode reference resource, a comprehensive collection of ITS1 sequences supplied with robust taxonomies, is highly needed. To address this issue, we created ITSoneDB (available at http://itsonedb.cloud.ba.infn.it/) which in its current version hosts 985 240 ITS1 sequences spanning over 134 000 eukaryotic species. Each ITS1 is mapped on the NCBI reference taxonomy with its start and end positions precisely annotated. ITSoneDB has been developed in agreement to the FAIR guidelines by enabling the users to query and download its content through a simple web-interface and access relevant metadata by cross-linking to European Nucleotide Archive.

WoPPER: Web server for Position Related data analysis of gene Expression in Prokaryotes.

The structural and conformational organization of chromosomes is crucial for gene expression regulation in eukaryotes and prokaryotes as well. Up to date, gene expression data generated using either microarray or RNA-sequencing are available for many bacterial genomes. However, differential gene expression is usually investigated with methods considering each gene independently, thus not taking into account the physical localization of genes along a bacterial chromosome. Here, we present WoPPER, a web tool integrating gene expression and genomic annotations to identify differentially expressed chromosomal regions in bacteria. RNA-sequencing or microarray-based gene expression data are provided as input, along with gene annotations. The user can select genomic annotations from an internal database including 2780 bacterial strains, or provide custom genomic annotations. The analysis produces as output the lists of positionally related genes showing a coordinated trend of differential expression. Graphical representations, including a circular plot of the analyzed chromosome, allow intuitive browsing of the results. The analysis procedure is based on our previously published R-package PREDA. The release of this tool is timely and relevant for the scientific community, as WoPPER will fill an existing gap in prokaryotic gene expression data analysis and visualization tools. WoPPER is open to all users and can be reached at the following URL: https://WoPPER.ba.itb.cnr.it.

Arena-Idb: a platform to build human non-coding RNA interaction networks.

BACKGROUND: High throughput technologies have provided the scientific community an unprecedented opportunity for large-scale analysis of genomes. Non-coding RNAs (ncRNAs), for a long time believed to be non-functional, are emerging as one of the most important and large family of gene regulators and key elements for genome maintenance. Functional studies have been able to assign to ncRNAs a wide spectrum of functions in primary biological processes, and for this reason they are assuming a growing importance as a potential new family of cancer therapeutic targets. Nevertheless, the number of functionally characterized ncRNAs is still too poor if compared to the number of new discovered ncRNAs. Thus platforms able to merge information from available resources addressing data integration issues are necessary and still insufficient to elucidate ncRNAs biological roles. RESULTS: In this paper, we describe a platform called Arena-Idb for the retrieval of comprehensive and non-redundant annotated ncRNAs interactions. Arena-Idb provides a framework for network reconstruction of ncRNA heterogeneous interactions (i.e., with other type of molecules) and relationships with human diseases which guide the integration of data, extracted from different sources, via mapping of entities and minimization of ambiguity. CONCLUSIONS: Arena-Idb provides a schema and a visualization system to integrate ncRNA interactions that assists in discovering ncRNA functions through the extraction of heterogeneous interaction networks. The Arena-Idb is available at http://arenaidb.ba.itb.cnr.it.

Investigating the Role of MicroRNA and Transcription Factor Co-regulatory Networks in Multiple Sclerosis Pathogenesis.

MicroRNAs (miRNAs) and transcription factors (TFs) play key roles in complex multifactorial diseases like multiple sclerosis (MS). Starting from the miRNomic profile previously associated with a cohort of pediatric MS (PedMS) patients, we applied a combined molecular and computational approach in order to verify published data in patients with adult-onset MS (AOMS). Six out of the 13 selected miRNAs (miR-320a, miR-125a-5p, miR-652-3p, miR-185-5p, miR-942-5p, miR-25-3p) were significantly upregulated in PedMS and AOMS patients, suggesting that they may be considered circulating biomarkers distinctive of the disease independently from age. A computational and unbiased miRNA-based screening of target genes not necessarily associated to MS was then performed in order to provide an extensive view of the genetic mechanisms underlying the disease. A comprehensive MS-specific miRNA-TF co-regulatory network was hypothesized; among others, SP1, RELA, NF-κB, TP53, AR, MYC, HDAC1, and STAT3 regulated the transcription of 61 targets. Interestingly, NF-κB and STAT3 cooperatively regulate the expression of immune response genes and control the cross-talk between inflammatory and immune cells. Further functional analysis will be performed on the identified critical hubs. Above all, in our view, this approach supports the need of multidisciplinary strategies for shedding light into the pathogenesis of MS.

TRIM8 restores p53 tumour suppressor function by blunting N-MYC activity in chemo-resistant tumours.

BACKGROUND: TRIM8 plays a key role in controlling the p53 molecular switch that sustains the transcriptional activation of cell cycle arrest genes and response to chemotherapeutic drugs. The mechanisms that regulate TRIM8, especially in cancers like clear cell Renal Cell Carcinoma (ccRCC) and colorectal cancer (CRC) where it is low expressed, are still unknown. However, recent studies suggest the potential involvement of some microRNAs belonging to miR-17-92 and its paralogous clusters, which could include TRIM8 in a more complex pathway. METHODS: We used RCC and CRC cell models for in-vitro experiments, and ccRCC patients and xenograft transplanted mice for in vivo assessments. To measure microRNAs levels we performed RT-qPCR, while steady-states of TRIM8, p53, p21 and N-MYC were quantified at protein level by Western Blotting as well as at transcript level by RT-qPCR. Luciferase reporter assays were performed to assess the interaction between TRIM8 and specific miRNAs, and the potential effects of this interaction on TRIM8 expression. Moreover, we treated our cell models with conventional chemotherapeutic drugs or tyrosine kinase inhibitors, and measured their response in terms of cell proliferation by MTT and colony suppression assays. RESULTS: We showed that TRIM8 is a target of miR-17-5p and miR-106b-5p, whose expression is promoted by N-MYC, and that alterations of their levels affect cell proliferation, acting on the TRIM8 transcripts stability, as confirmed in ccRCC patients and cell lines. In addition, reducing the levels of miR-17-5p/miR-106b-5p, we increased the chemo-sensitivity of RCC/CRC-derived cells to anti-tumour drugs used in the clinic. Intriguingly, this occurs, on one hand, by recovering the p53 tumour suppressor activity in a TRIM8-dependent fashion and, on the other hand, by promoting the transcription of miR-34a that turns off the oncogenic action of N-MYC. This ultimately leads to cell proliferation reduction or block, observed also in colon cancer xenografts overexpressing TRIM8. CONCLUSIONS: In this paper we provided evidence that TRIM8 and its regulators miR-17-5p and miR-106b-5 participate to a feedback loop controlling cell proliferation through the reciprocal modulation of p53, miR-34a and N-MYC. Our experiments pointed out that this axis is pivotal in defining drug responsiveness of cancers such ccRCC and CRC.

A fuzzy method for RNA-Seq differential expression analysis in presence of multireads.

BACKGROUND: When the reads obtained from high-throughput RNA sequencing are mapped against a reference database, a significant proportion of them - known as multireads - can map to more than one reference sequence. These multireads originate from gene duplications, repetitive regions or overlapping genes. Removing the multireads from the mapping results, in RNA-Seq analyses, causes an underestimation of the read counts, while estimating the real read count can lead to false positives during the detection of differentially expressed sequences. RESULTS: We present an innovative approach to deal with multireads and evaluate differential expression events, entirely based on fuzzy set theory. Since multireads cause uncertainty in the estimation of read counts during gene expression computation, they can also influence the reliability of differential expression analysis results, by producing false positives. Our method manages the uncertainty in gene expression estimation by defining the fuzzy read counts and evaluates the possibility of a gene to be differentially expressed with three fuzzy concepts: over-expression, same-expression and under-expression. The output of the method is a list of differentially expressed genes enriched with information about the uncertainty of the results due to the multiread presence. We have tested the method on RNA-Seq data designed for case-control studies and we have compared the obtained results with other existing tools for read count estimation and differential expression analysis. CONCLUSIONS: The management of multireads with the use of fuzzy sets allows to obtain a list of differential expression events which takes in account the uncertainty in the results caused by the presence of multireads. Such additional information can be used by the biologists when they have to select the most relevant differential expression events to validate with laboratory assays. Our method can be used to compute reliable differential expression events and to highlight possible false positives in the lists of differentially expressed genes computed with other tools.

Meta-Analysis of Differential Connectivity in Gene Co-Expression Networks in Multiple Sclerosis.

Differential gene expression analyses to investigate multiple sclerosis (MS) molecular pathogenesis cannot detect genes harboring genetic and/or epigenetic modifications that change the gene functions without affecting their expression. Differential co-expression network approaches may capture changes in functional interactions resulting from these alterations. We re-analyzed 595 mRNA arrays from publicly available datasets by studying changes in gene co-expression networks in MS and in response to interferon (IFN)-ß treatment. Interestingly, MS networks show a reduced connectivity relative to the healthy condition, and the treatment activates the transcription of genes and increases their connectivity in MS patients. Importantly, the analysis of changes in gene connectivity in MS patients provides new evidence of association for genes already implicated in MS by single-nucleotide polymorphism studies and that do not show differential expression. This is the case of amiloride-sensitive cation channel 1 neuronal (ACCN1) that shows a reduced number of interacting partners in MS networks, and it is known for its role in synaptic transmission and central nervous system (CNS) development. Furthermore, our study confirms a deregulation of the vitamin D system: among the transcription factors that potentially regulate the deregulated genes, we find TCF3 and SP1 that are both involved in vitamin D3-induced p27Kip1 expression. Unveiling differential network properties allows us to gain systems-level insights into disease mechanisms and may suggest putative targets for the treatment.

Clusterin transcript variants expression in thyroid tumor: a potential marker of malignancy?

BACKGROUND: Clusterin (CLU) is a ubiquitous multifunctional factor involved in neoplastic transformation. The CLU transcript variants and protein forms play a crucial role in balancing cells proliferation and death. METHODS: We investigated the regulation of CLU transcript variants expression in an in vivo model system consisting of both neoplastic tissues and fine needle aspiration biopsy (FNAB) samples isolated from patients undergoing thyroidectomy. RESULTS: The immunohistochemical analyses showed an overall CLU up-regulation in papillary carcinoma. A specific CLU2 transcript variant increase was registered using qPCR in papillary carcinomas while CLU1 decreased. In addition, the analysis of CLU transcripts expression level showed an increase of the CLU2 transcript in the TIR 3 patients with histologically confirmed thyroid cancer. CONCLUSIONS: Our results suggest the existence of a specific alteration of CLU2:CLU1 ratio towards CLU2, thus providing the first circumstantial evidence for the potential use of CLU transcript variants as effective biomarkers for a more accurate assessment of the so called "indeterminate" thyroid nodules.

Reference databases for taxonomic assignment in metagenomics.

Metagenomics is providing an unprecedented access to the environmental microbial diversity. The amplicon-based metagenomics approach involves the PCR-targeted sequencing of a genetic locus fitting different features. Namely, it must be ubiquitous in the taxonomic range of interest, variable enough to discriminate between different species but flanked by highly conserved sequences, and of suitable size to be sequenced through next-generation platforms. The internal transcribed spacers 1 and 2 (ITS1 and ITS2) of the ribosomal DNA operon and one or more hyper-variable regions of 16S ribosomal RNA gene are typically used to identify fungal and bacterial species, respectively. In this context, reliable reference databases and taxonomies are crucial to assign amplicon sequence reads to the correct phylogenetic ranks. Several resources provide consistent phylogenetic classification of publicly available 16S ribosomal DNA sequences, whereas the state of ribosomal internal transcribed spacers reference databases is notably less advanced. In this review, we aim to give an overview of existing reference resources for both types of markers, highlighting strengths and possible shortcomings of their use for metagenomics purposes. Moreover, we present a new database, ITSoneDB, of well annotated and phylogenetically classified ITS1 sequences to be used as a reference collection in metagenomic studies of environmental fungal communities. ITSoneDB is available for download and browsing at http://itsonedb.ba.itb.cnr.it/.

A platform independent RNA-Seq protocol for the detection of transcriptome complexity.

BACKGROUND: Recent studies have demonstrated an unexpected complexity of transcription in eukaryotes. The majority of the genome is transcribed and only a little fraction of these transcripts is annotated as protein coding genes and their splice variants. Indeed, most transcripts are the result of antisense, overlapping and non-coding RNA expression. In this frame, one of the key aims of high throughput transcriptome sequencing is the detection of all RNA species present in the cell and the first crucial step for RNA-seq users is represented by the choice of the strategy for cDNA library construction. The protocols developed so far provide the utilization of the entire library for a single sequencing run with a specific platform. RESULTS: We set up a unique protocol to generate and amplify a strand-specific cDNA library representative of all RNA species that may be implemented with all major platforms currently available on the market (Roche 454, Illumina, ABI/SOLiD). Our method is reproducible, fast, easy-to-perform and even allows to start from low input total RNA. Furthermore, we provide a suitable bioinformatics tool for the analysis of the sequences produced following this protocol. CONCLUSION: We tested the efficiency of our strategy, showing that our method is platform-independent, thus allowing the simultaneous analysis of the same sample with different NGS technologies, and providing an accurate quantitative and qualitative portrait of complex whole transcriptomes.

BEAT: Bioinformatics Exon Array Tool to store, analyze and visualize Affymetrix GeneChip Human Exon Array data from disease experiments.

BACKGROUND: It is known from recent studies that more than 90% of human multi-exon genes are subject to Alternative Splicing (AS), a key molecular mechanism in which multiple transcripts may be generated from a single gene. It is widely recognized that a breakdown in AS mechanisms plays an important role in cellular differentiation and pathologies. Polymerase Chain Reactions, microarrays and sequencing technologies have been applied to the study of transcript diversity arising from alternative expression. Last generation Affymetrix GeneChip Human Exon 1.0 ST Arrays offer a more detailed view of the gene expression profile providing information on the AS patterns. The exon array technology, with more than five million data points, can detect approximately one million exons, and it allows performing analyses at both gene and exon level. In this paper we describe BEAT, an integrated user-friendly bioinformatics framework to store, analyze and visualize exon arrays datasets. It combines a data warehouse approach with some rigorous statistical methods for assessing the AS of genes involved in diseases. Meta statistics are proposed as a novel approach to explore the analysis results. BEAT is available at http://beat.ba.itb.cnr.it. RESULTS: BEAT is a web tool which allows uploading and analyzing exon array datasets using standard statistical methods and an easy-to-use graphical web front-end. BEAT has been tested on a dataset with 173 samples and tuned using new datasets of exon array experiments from 28 colorectal cancer and 26 renal cell cancer samples produced at the Medical Genetics Unit of IRCCS Casa Sollievo della Sofferenza.To highlight all possible AS events, alternative names, accession Ids, Gene Ontology terms and biochemical pathways annotations are integrated with exon and gene level expression plots. The user can customize the results choosing custom thresholds for the statistical parameters and exploiting the available clinical data of the samples for a multivariate AS analysis. CONCLUSIONS: Despite exon array chips being widely used for transcriptomics studies, there is a lack of analysis tools offering advanced statistical features and requiring no programming knowledge. BEAT provides a user-friendly platform for a comprehensive study of AS events in human diseases, displaying the analysis results with easily interpretable and interactive tables and graphics.

Bioinformatics approaches for genomics and post genomics applications of next-generation sequencing.

Technical advances such as the development of molecular cloning, Sanger sequencing, PCR and oligonucleotide microarrays are key to our current capacity to sequence, annotate and study complete organismal genomes. Recent years have seen the development of a variety of so-called 'next-generation' sequencing platforms, with several others anticipated to become available shortly. The previously unimaginable scale and economy of these methods, coupled with their enthusiastic uptake by the scientific community and the potential for further improvements in accuracy and read length, suggest that these technologies are destined to make a huge and ongoing impact upon genomic and post-genomic biology. However, like the analysis of microarray data and the assembly and annotation of complete genome sequences from conventional sequencing data, the management and analysis of next-generation sequencing data requires (and indeed has already driven) the development of informatics tools able to assemble, map, and interpret huge quantities of relatively or extremely short nucleotide sequence data. Here we provide a broad overview of bioinformatics approaches that have been introduced for several genomics and functional genomics applications of next-generation sequencing.

UTRdb and UTRsite (RELEASE 2010): a collection of sequences and regulatory motifs of the untranslated regions of eukaryotic mRNAs.

The 5' and 3' untranslated regions of eukaryotic mRNAs (UTRs) play crucial roles in the post-transcriptional regulation of gene expression through the modulation of nucleo-cytoplasmic mRNA transport, translation efficiency, subcellular localization and message stability. UTRdb is a curated database of 5' and 3' untranslated sequences of eukaryotic mRNAs, derived from several sources of primary data. Experimentally validated functional motifs are annotated and also collated as the UTRsite database where more specific information on the functional motifs and cross-links to interacting regulatory protein are provided. In the current update, the UTR entries have been organized in a gene-centric structure to better visualize and retrieve 5' and 3'UTR variants generated by alternative initiation and termination of transcription and alternative splicing. Experimentally validated miRNA targets and conserved sequence elements are also annotated. The integration of UTRdb with genomic data has allowed the implementation of an efficient annotation system and a powerful retrieval resource for the selection and extraction of specific UTR subsets. All internet resources implemented for retrieval and functional analysis of 5' and 3' untranslated regions of eukaryotic mRNAs are accessible at http://utrdb.ba.itb.cnr.it/.

Transmembrane Protein TMEM230, a Target of Glioblastoma Therapy.

Glioblastomas (GBM) are the most aggressive tumors originating in the brain. Histopathologic features include circuitous, disorganized, and highly permeable blood vessels with intermittent blood flow. These features contribute to the inability to direct therapeutic agents to tumor cells. Known targets for anti-angiogenic therapies provide minimal or no effect in overall survival of 12-15 months following diagnosis. Identification of novel targets therefore remains an important goal for effective treatment of highly vascularized tumors such as GBM. We previously demonstrated in zebrafish that a balanced level of expression of the transmembrane protein TMEM230/C20ORF30 was required to maintain normal blood vessel structural integrity and promote proper vessel network formation. To investigate whether TMEM230 has a role in the pathogenesis of GBM, we analyzed its prognostic value in patient tumor gene expression datasets and performed cell functional analysis. TMEM230 was found necessary for growth of U87-MG cells, a model of human GBM. Downregulation of TMEM230 resulted in loss of U87 migration, substratum adhesion, and re-passaging capacity. Conditioned media from U87 expressing endogenous TMEM230 induced sprouting and tubule-like structure formation of HUVECs. Moreover, TMEM230 promoted vascular mimicry-like behavior of U87 cells. Gene expression analysis of 702 patients identified that TMEM230 expression levels distinguished high from low grade gliomas. Transcriptomic analysis of patients with gliomas revealed molecular pathways consistent with properties observed in U87 cell assays. Within low grade gliomas, elevated TMEM230 expression levels correlated with reduced overall survival independent from tumor subtype. Highest level of TMEM230 correlated with glioblastoma and ATP-dependent microtubule kinesin motor activity, providing a direction for future therapeutic intervention. Our studies support that TMEM230 has both glial tumor and endothelial cell intracellular and extracellular functions. Elevated levels of TMEM230 promote glial tumor cell migration, extracellular scaffold remodeling, and hypervascularization and abnormal formation of blood vessels. Downregulation of TMEM230 expression may inhibit both low grade glioma and glioblastoma tumor progression and promote normalization of abnormally formed blood vessels. TMEM230 therefore is both a promising anticancer and antiangiogenic therapeutic target for inhibiting GBM tumor cells and tumor-driven angiogenesis.

Plant miRNAs Reduce Cancer Cell Proliferation by Targeting MALAT1 and NEAT1: A Beneficial Cross-Kingdom Interaction.

MicroRNAs (miRNAs) are ubiquitous regulators of gene expression, evolutionarily conserved in plants and mammals. In recent years, although a growing number of papers debate the role of plant miRNAs on human gene expression, the molecular mechanisms through which this effect is achieved are still not completely elucidated. Some evidence suggest that this interaction might be sequence specific, and in this work, we investigated this possibility by transcriptomic and bioinformatics approaches. Plant and human miRNA sequences from primary databases were collected and compared for their similarities (global or local alignments). Out of 2,588 human miRNAs, 1,606 showed a perfect match of their seed sequence with the 5' end of 3,172 plant miRNAs. Further selections were applied based on the role of the human target genes or of the miRNA in cell cycle regulation (as an oncogene, tumor suppressor, or a biomarker for prognosis, or diagnosis in cancer). Based on these criteria, 20 human miRNAs were selected as potential functional analogous of 7 plant miRNAs, which were in turn transfected in different cell lines to evaluate their effect on cell proliferation. A significant decrease was observed in colorectal carcinoma HCT116 cell line. RNA-Seq demonstrated that 446 genes were differentially expressed 72 h after transfection. Noteworthy, we demonstrated that the plant mtr-miR-5754 and gma-miR4995 directly target the tumor-associated long non-coding RNA metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) and nuclear paraspeckle assembly transcript 1 (NEAT1) in a sequence-specific manner. In conclusion, according to other recent discoveries, our study strengthens and expands the hypothesis that plant miRNAs can have a regulatory effect in mammals by targeting both protein-coding and non-coding RNA, thus suggesting new biotechnological applications.

Integrated Analysis of microRNA and mRNA Expression Profiles: An Attempt to Disentangle the Complex Interaction Network in Attention Deficit Hyperactivity Disorder.

Attention Deficit Hyperactivity Disorder (ADHD) is a childhood-onset neurodevelopmental disorder, whose etiology and pathogenesis are still largely unknown. In order to uncover novel regulatory networks and molecular pathways possibly related to ADHD, we performed an integrated miRNA and mRNA expression profiling analysis in peripheral blood samples of children with ADHD and age-matched typically developing (TD) children. The expression levels of 13 miRNAs were evaluated with microfluidic qPCR, and differentially expressed (DE) mRNAs were detected on an Illumina HiSeq 2500 genome analyzer. The miRNA targetome was identified using an integrated approach of validated and predicted interaction data extracted from seven different bioinformatic tools. Gene Ontology (GO) and pathway enrichment analyses were carried out. Results showed that six miRNAs (miR-652-3p, miR-942-5p, let-7b-5p, miR-181a-5p, miR-320a, and miR-148b-3p) and 560 genes were significantly DE in children with ADHD compared to TD subjects. After correction for multiple testing, only three miRNAs (miR-652-3p, miR-148b-3p, and miR-942-5p) remained significant. Genes known to be associated with ADHD (e.g., B4GALT2, SLC6A9 TLE1, ANK3, TRIO, TAF1, and SYNE1) were confirmed to be significantly DE in our study. Integrated miRNA and mRNA expression data identified critical key hubs involved in ADHD. Finally, the GO and pathway enrichment analyses of all DE genes showed their deep involvement in immune functions, reinforcing the hypothesis that an immune imbalance might contribute to the ADHD etiology. Despite the relatively small sample size, in this study we were able to build a complex miRNA-target interaction network in children with ADHD that might help in deciphering the disease pathogenesis. Validation in larger samples should be performed in order to possibly suggest novel therapeutic strategies for treating this complex disease.

p53FamTaG: a database resource of human p53, p63 and p73 direct target genes combining in silico prediction and microarray data.

BACKGROUND: The p53 gene family consists of the three genes p53, p63 and p73, which have polyhedral non-overlapping functions in pivotal cellular processes such as DNA synthesis and repair, growth arrest, apoptosis, genome stability, angiogenesis, development and differentiation. These genes encode sequence-specific nuclear transcription factors that recognise the same responsive element (RE) in their target genes. Their inactivation or aberrant expression may determine tumour progression or developmental disease. The discovery of several protein isoforms with antagonistic roles, which are produced by the expression of different promoters and alternative splicing, widened the complexity of the scenario of the transcriptional network of the p53 family members. Therefore, the identification of the genes transactivated by p53 family members is crucial to understand the specific role for each gene in cell cycle regulation. We have combined a genome-wide computational search of p53 family REs and microarray analysis to identify new direct target genes. The huge amount of biological data produced has generated a critical need for bioinformatic tools able to manage and integrate such data and facilitate their retrieval and analysis. DESCRIPTION: We have developed the p53FamTaG database (p53 FAMily TArget Genes), a modular relational database, which contains p53 family direct target genes selected in the human genome searching for the presence of the REs and the expression profile of these target genes obtained by microarray experiments. p53FamTaG database also contains annotations of publicly available databases and links to other experimental data. The genome-wide computational search of the REs was performed using PatSearch, a pattern-matching program implemented in the DNAfan tool. These data were integrated with the microarray results we produced from the overexpression of different isoforms of p53, p63 and p73 stably transfected in isogenic cell lines, allowing the comparative study of the transcriptional activity of all the proteins in the same cellular background.p53FamTaG database is available free at http://www2.ba.itb.cnr.it/p53FamTaG/ CONCLUSION: p53FamTaG represents a unique integrated resource of human direct p53 family target genes that is extensively annotated and provides the users with an efficient query/retrieval system which displays the results of our microarray experiments and allows the export of RE sequences. The database was developed for supporting and integrating high-throughput in silico and experimental analyses and represents an important reference source of knowledge for research groups involved in the field of oncogenesis, apoptosis and cell cycle regulation.

Selection of relevant genes in cancer diagnosis based on their prediction accuracy.

MOTIVATIONS: One of the main problems in cancer diagnosis by using DNA microarray data is selecting genes relevant for the pathology by analyzing their expression profiles in tissues in two different phenotypical conditions. The question we pose is the following: how do we measure the relevance of a single gene in a given pathology? METHODS: A gene is relevant for a particular disease if we are able to correctly predict the occurrence of the pathology in new patients on the basis of its expression level only. In other words, a gene is informative for the disease if its expression levels are useful for training a classifier able to generalize, that is, able to correctly predict the status of new patients. In this paper we present a selection bias free, statistically well founded method for finding relevant genes on the basis of their classification ability. RESULTS: We applied the method on a colon cancer data set and produced a list of relevant genes, ranked on the basis of their prediction accuracy. We found, out of more than 6500 available genes, 54 overexpressed in normal tissues and 77 overexpressed in tumor tissues having prediction accuracy greater than 70% with p-value