RAP: RNA-Seq Analysis Pipeline, a new cloud-based NGS web application.

BACKGROUND: The study of RNA has been dramatically improved by the introduction of Next Generation Sequencing platforms allowing massive and cheap sequencing of selected RNA fractions, also providing information on strand orientation (RNA-Seq). The complexity of transcriptomes and of their regulative pathways make RNA-Seq one of most complex field of NGS applications, addressing several aspects of the expression process (e.g. identification and quantification of expressed genes and transcripts, alternative splicing and polyadenylation, fusion genes and trans-splicing, post-transcriptional events, etc.). METHODS: In order to provide researchers with an effective and friendly resource for analyzing RNA-Seq data, we present here RAP (RNA-Seq Analysis Pipeline), a cloud computing web application implementing a complete but modular analysis workflow. This pipeline integrates both state-of-the-art bioinformatics tools for RNA-Seq analysis and in-house developed scripts to offer to the user a comprehensive strategy for data analysis. RAP is able to perform quality checks (adopting FastQC and NGS QC Toolkit), identify and quantify expressed genes and transcripts (with Tophat, Cufflinks and HTSeq), detect alternative splicing events (using SpliceTrap) and chimeric transcripts (with ChimeraScan). This pipeline is also able to identify splicing junctions and constitutive or alternative polyadenylation sites (implementing custom analysis modules) and call for statistically significant differences in genes and transcripts expression, splicing pattern and polyadenylation site usage (using Cuffdiff2 and DESeq). RESULTS: Through a user friendly web interface, the RAP workflow can be suitably customized by the user and it is automatically executed on our cloud computing environment. This strategy allows to access to bioinformatics tools and computational resources without specific bioinformatics and IT skills. RAP provides a set of tabular and graphical results that can be helpful to browse, filter and export analyzed data, according to the user needs.

SpliceAid-F: a database of human splicing factors and their RNA-binding sites.

A comprehensive knowledge of all the factors involved in splicing, both proteins and RNAs, and of their interaction network is crucial for reaching a better understanding of this process and its functions. A large part of relevant information is buried in the literature or collected in various different databases. By hand-curated screenings of literature and databases, we retrieved experimentally validated data on 71 human RNA-binding splicing regulatory proteins and organized them into a database called 'SpliceAid-F' (http://www.caspur.it/SpliceAidF/). For each splicing factor (SF), the database reports its functional domains, its protein and chemical interactors and its expression data. Furthermore, we collected experimentally validated RNA-SF interactions, including relevant information on the RNA-binding sites, such as the genes where these sites lie, their genomic coordinates, the splicing effects, the experimental procedures used, as well as the corresponding bibliographic references. We also collected information from experiments showing no RNA-SF binding, at least in the assayed conditions. In total, SpliceAid-F contains 4227 interactions, 2590 RNA-binding sites and 1141 'no-binding' sites, including information on cellular contexts and conditions where binding was tested. The data collected in SpliceAid-F can provide significant information to explain an observed splicing pattern as well as the effect of mutations in functional regulatory elements.

MitoZoa 2.0: a database resource and search tools for comparative and evolutionary analyses of mitochondrial genomes in Metazoa.

The MITOchondrial genome database of metaZOAns (MitoZoa) is a public resource for comparative analyses of metazoan mitochondrial genomes (mtDNA) at both the sequence and genomic organizational levels. The main characteristics of the MitoZoa database are the careful revision of mtDNA entry annotations and the possibility of retrieving gene order and non-coding region (NCR) data in appropriate formats. The MitoZoa retrieval system enables basic and complex queries at various taxonomic levels using different search menus. MitoZoa 2.0 has been enhanced in several aspects, including: a re-annotation pipeline to check the correctness of protein-coding gene predictions; a standardized annotation of introns and of precursor ORFs whose functionality is post-transcriptionally recovered by RNA editing or programmed translational frameshifting; updates of taxon-related fields and a BLAST sequence similarity search tool. Database novelties and the definition of standard mtDNA annotation rules, together with the user-friendly retrieval system and the BLAST service, make MitoZoa a valuable resource for comparative and evolutionary analyses as well as a reference database to assist in the annotation of novel mtDNA sequences. MitoZoa is freely accessible at http://www.caspur.it/mitozoa.

WEP: a high-performance analysis pipeline for whole-exome data.

BACKGROUND: The advent of massively parallel sequencing technologies (Next Generation Sequencing, NGS) profoundly modified the landscape of human genetics.In particular, Whole Exome Sequencing (WES) is the NGS branch that focuses on the exonic regions of the eukaryotic genomes; exomes are ideal to help us understanding high-penetrance allelic variation and its relationship to phenotype. A complete WES analysis involves several steps which need to be suitably designed and arranged into an efficient pipeline.Managing a NGS analysis pipeline and its huge amount of produced data requires non trivial IT skills and computational power. RESULTS: Our web resource WEP (Whole-Exome sequencing Pipeline web tool) performs a complete WES pipeline and provides easy access through interface to intermediate and final results. The WEP pipeline is composed of several steps:1) verification of input integrity and quality checks, read trimming and filtering; 2) gapped alignment; 3) BAM conversion, sorting and indexing; 4) duplicates removal; 5) alignment optimization around insertion/deletion (indel) positions; 6) recalibration of quality scores; 7) single nucleotide and deletion/insertion polymorphism (SNP and DIP) variant calling; 8) variant annotation; 9) result storage into custom databases to allow cross-linking and intersections, statistics and much more. In order to overcome the challenge of managing large amount of data and maximize the biological information extracted from them, our tool restricts the number of final results filtering data by customizable thresholds, facilitating the identification of functionally significant variants. Default threshold values are also provided at the analysis computation completion, tuned with the most common literature work published in recent years. CONCLUSIONS: Through our tool a user can perform the whole analysis without knowing the underlying hardware and software architecture, dealing with both paired and single end data. The interface provides an easy and intuitive access for data submission and a user-friendly web interface for annotated variant visualization.Non-IT mastered users can access through WEP to the most updated and tested WES algorithms, tuned to maximize the quality of called variants while minimizing artifacts and false positives.The web tool is available at the following web address: http://www.caspur.it/wep.

ASPicDB: a database of annotated transcript and protein variants generated by alternative splicing.

Alternative splicing is emerging as a major mechanism for the expansion of the transcriptome and proteome diversity, particularly in human and other vertebrates. However, the proportion of alternative transcripts and proteins actually endowed with functional activity is currently highly debated. We present here a new release of ASPicDB which now provides a unique annotation resource of human protein variants generated by alternative splicing. A total of 256,939 protein variants from 17,191 multi-exon genes have been extensively annotated through state of the art machine learning tools providing information of the protein type (globular and transmembrane), localization, presence of PFAM domains, signal peptides, GPI-anchor propeptides, transmembrane and coiled-coil segments. Furthermore, full-length variants can be now specifically selected based on the annotation of CAGE-tags and polyA signal and/or polyA sites, marking transcription initiation and termination sites, respectively. The retrieval can be carried out at gene, transcript, exon, protein or splice site level allowing the selection of data sets fulfilling one or more features settled by the user. The retrieval interface also enables the selection of protein variants showing specific differences in the annotated features. ASPicDB is available at http://www.caspur.it/ASPicDB/.

Transcriptional regulation of miR-224 upregulated in human HCCs by NFκB inflammatory pathways.

BACKGROUND & AIMS: miR-224 is up-regulated in human HCCs as compared to both paired peri-tumoral cirrhotic tissues and cirrhotic livers without HCC. Here, we have cloned the miR-224 regulatory region and characterized its transcriptional regulation by the NFκB-dependent inflammatory pathways. METHODS: Mature miRNA expression was evaluated by a 2 step stem-loop real-time RT-PCR. The recruitment of polymerase II and transcription factors on the pre-miR-224 promoter has been assessed by ChIPSeq and ChIP. RESULTS: We found miR-224 levels strongly up-regulated in both peri-tumoral cirrhotic livers and HCC samples as compared to normal livers. In silico analysis of the putative miR-224 promoter revealed multiple NFκB sites. We showed that LTα and TNFα activate transcription from the miR-224 promoter and of endogenous miR-224 expression in HCC cell lines, whereas the expression of miR-224 target API5 was reduced. Exogenously expressed p65/RelA activates the miR-224 promoter and a dominant negative form of IκBα (IκBSR) represses it. ChIP analysis showed that p65/NFκB is recruited on the miR-224 promoter and that its binding sharply increases after exposure to LPS, TNFα, and LTα. Altogether these findings link the inflammatory signals to NFκB-mediated activation of miR-224 expression. An antago-miR specific for miR-224 blocked LPS and LTα stimulated HCC cells migration and invasion. Conversely, the IKK inhibitor BMS-345541 blocks pre-miR-224-induced cellular migration and invasion. CONCLUSIONS: Our results identify p65/NFκB as a direct transcriptional regulator of miR-224 expression and link miR-224 up-regulation with the activation of the LPS, LTα, and TNFα inflammatory pathways and cell migration/invasion in HCC.

VDR primary targets by genome-wide transcriptional profiling.

There is growing evidence that 1α,25-dihydroxyvitamin D3 (1α,25(OH)2D3) plays a role in breast cancer prevention and survival. It elicits a variety of antitumor activities like controlling cellular differentiation, proliferation and angiogenesis. Most of its biological effects are exerted via its nuclear receptor which acts as a transcriptional regulator. Here, we carried out a genome-wide investigation of the primary transcriptional targets of 1α,25(OH)2D3 in breast epithelial cancer cells using RNA-Seq technology. We identified early transcriptional targets of 1α,25(OH)2D3 involved in adhesion, growth regulation, angiogenesis, actin cytoskeleton regulation, hexose transport, inflammation and immunomodulation, apoptosis, endocytosis and signaling. Furthermore, we found several transcription factors to be regulated by 1α,25(OH)2D3 that subsequently amplify and diversify the transcriptional output driven by 1α,25(OH)2D3 leading finally to a growth arrest of the cells. Moreover, we could show that 1α,25(OH)2D3 elevates the trimethylation of histone H3 lysine 4 at several target gene promoters. Our present transcriptomic analysis of differential expression after 1α,25(OH)2D3 treatment provides a resource of primary 1α,25(OH)2D3 targets that might drive the antiproliferative action in breast cancer epithelial cells.