Paolo D'Onorio De Meo

The PMDB Protein Model Database

The Protein Model Database (PMDB) is a public resource aimed at storing manually built 3D models of proteins. The database is designed to provide access to models published in the scientific literature, together with validating experimental data. It is a relational database and it currently contains >74 000 models for ∼240 proteins. The system is accessible at and allows predictors to submit models along with related supporting evidence and users to download them through a simple and intuitive interface. Users can navigate in the database and retrieve models referring to the same target protein or to different regions of the same protein. Each model is assigned a unique identifier that allows interested users to directly access the data.

Developmental factor IRF6 exhibits tumor suppressor activity in squamous cell carcinomas

The transcription factor interferon regulatory factor 6 (IRF6) regulates craniofacial development and epidermal proliferation. We recently showed that IRF6 is a component of a regulatory feedback loop that controls the proliferative potential of epidermal cells. IRF6 is transcriptionally activated by p63 and induces its proteasome-mediated down-regulation, thereby limiting keratinocyte proliferative potential. We hypothesized that IRF6 may also be involved in skin carcinogenesis. Hence, we analyzed IRF6 expression in a large series of squamous cell carcinomas (SCCs) and found a strong down-regulation of IRF6 that correlated with tumor invasive and differentiation status. IRF6 down-regulation in SCC cell lines and primary tumors correlates with methylation on a CpG dinucleotide island located in its promoter region. To identify the molecular mechanisms regulating IRF6 potential tumor suppressive activity, we performed a genome-wide analysis by combining ChIP sequencing for IRF6 binding sites and gene expression profiling in primary human keratinocytes after siRNA-mediated IRF6 depletion. We observed dysregulation of cell cycle-related genes and genes involved in differentiation, cell adhesion, and cell–cell contact. Many of these genes were direct IRF6 targets. We also performed in vitro invasion assays showing that IRF6 down-regulation promotes invasive behavior and that reintroduction of IRF6 into SCC cells strongly inhibits cell growth. These results indicate a function for IRF6 in suppression of tumorigenesis in stratified epithelia.

MitoZoa: A curated mitochondrial genome database of metazoans for comparative genomics studies

MitoZoa is a relational database collecting curated metazoan entries of complete or nearly complete mitochondrial genomes (mtDNA), specifically designed to assist comparative studies of mitochondrial genome-level features in a given taxon or in congeneric species of Metazoa. The principal novelties of MitoZoa are extensive corrections/improvements of the mtDNA annotations and the possibility of easily searching for data on: (1) gene order, a genomic feature useful as phylogenetic marker; (2) sequence, size and location of non-coding regions, likely containing the regulatory signals for mtDNA replication and transcription; (3) mt features/sequences of congeneric species, where saturation phenomena in nucleotide substitutions and gene order changes are expected to be absent or at least minimal. In addition, MitoZoa allows the exploration of basic mt features such as molecule topology, genetic code, gene content, and compositional parameters of the entire genome. Finally, in order to facilitate downstream analyses of retrieved data, MitoZoa entry lists can be visualized and downloaded in a tabular format, while sequences and gene order data are provided in FASTA and FASTA-like formats, respectively. The MitoZoa database is available at http://www.caspur.it/mitozoa.

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/.

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

BackgroundThe 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.ResultsOur 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.ConclusionsThrough 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

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.

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

BackgroundThe 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.).Moreover, the huge volume of data generated by NGS platforms introduces unprecedented computational and technological challenges to efficiently analyze and store sequence data and results.MethodsIn 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).ResultsThrough 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.

ASPIC: a web resource for alternative splicing prediction and transcript isoforms characterization

Alternative splicing (AS) is now emerging as a major mechanism contributing to the expansion of the transcriptome and proteome complexity of multicellular organisms. The fact that a single gene locus may give rise to multiple mRNAs and protein isoforms, showing both major and subtle structural variations, is an exceptionally versatile tool in the optimization of the coding capacity of the eukaryotic genome. The huge and continuously increasing number of genome and transcript sequences provides an essential information source for the computational detection of genes AS pattern. However, much of this information is not optimally or comprehensively used in gene annotation by current genome annotation pipelines. We present here a web resource implementing the ASPIC algorithm which we developed previously for the investigation of AS of user submitted genes, based on comparative analysis of available transcript and genome data from a variety of species. The ASPIC web resource provides graphical and tabular views of the splicing patterns of all full-length mRNA isoforms compatible with the detected splice sites of genes under investigation as well as relevant structural and functional annotation. The ASPIC web resource-available at http://www.caspur.it/ASPIC/--is dynamically interconnected with the Ensembl and Unigene databases and also implements an upload facility.

Che-1-induced inhibition of mTOR pathway enables stress-induced autophagy

Mammalian target of rapamycin (mTOR) is a key protein kinase that regulates cell growth, metabolism, and autophagy to maintain cellular homeostasis. Its activity is inhibited by adverse conditions, including nutrient limitation, hypoxia, and DNA damage. In this study, we demonstrate that Che‐1, a RNA polymerase II‐binding protein activated by the DNA damage response, inhibits mTOR activity in response to stress conditions. We found that, under stress, Che‐1 induces the expression of two important mTOR inhibitors, Redd1 and Deptor, and that this activity is required for sustaining stress‐induced autophagy. Strikingly, Che‐1 expression correlates with the progression of multiple myeloma and is required for cell growth and survival, a malignancy characterized by high autophagy response.