Stefano Castellana

Congruency in the prediction of pathogenic missense mutations: state-of-the-art web-based tools

A remarkable degree of genetic variation has been found in the protein-encoding regions of DNA through deep sequencing of samples obtained from thousands of subjects from several populations. Approximately half of the 20 000 single nucleotide polymorphisms present, even in normal healthy subjects, are nonsynonymous amino acid substitutions that could potentially affect protein function. The greatest challenges currently facing investigators are data interpretation and the development of strategies to identify the few gene-coding variants that actually cause or confer susceptibility to disease. A confusing array of options is available to address this problem. Unfortunately, the overall accuracy of these tools at ultraconserved positions is low, and predictions generated by current computational tools may mislead researchers involved in downstream experimental and clinical studies. First, we have presented an updated review of these tools and their primary functionalities, focusing on those that are naturally prone to analyze massive variant sets, to infer some interesting similarities among their results. Additionally, we have evaluated the prediction congruency for real whole-exome sequencing data in a proof-of-concept study on some of these web-based tools.

Mutations in B9D1 and MKS1 cause mild Joubert syndrome: expanding the genetic overlap with the lethal ciliopathy Meckel syndrome

Joubert syndrome is a clinically and genetically heterogeneous ciliopathy characterized by a typical cerebellar and brainstem malformation (the “molar tooth sign”), and variable multiorgan involvement. To date, 24 genes have been found mutated in Joubert syndrome, of which 13 also cause Meckel syndrome, a lethal ciliopathy with kidney, liver and skeletal involvement. Here we describe four patients with mild Joubert phenotypes who carry pathogenic mutations in either MKS1 or B9D1, two genes previously implicated only in Meckel syndrome.

MitImpact: an exhaustive collection of pre-computed pathogenicity predictions of human mitochondrial non-synonymous variants.

Mitochondrial DNA carries a tiny, but fundamental portion of the eukaryotic genetic code. As its nuclear counterpart, it is susceptible to point mutations. Their level of pathogenicity has been assessed for the newly discovered mutations only, leaving some degree of uncertainty on the potential impact of the unknown mutations. Here we present Mitochondrial mutation Impact (MitImpact), a queryable lightweight web interface to a reasoned collection of structurally and evolutionary annotated pathogenicity predictions, obtained by assembling pre‐computed with on‐the‐fly‐computed sets of pathogenicity estimations, for all the possible mitochondrial missense variants. It presents itself as a resource for fast and reliable evaluation of gene‐specific susceptibility of unknown and verified amino acid changes. MitImpact is freely available at http://bioinformatics.css‐mendel.it/ (tools section). ©2014 Wiley Periodicals, Inc.

Infantile and childhood onset PLA2G6-associated neurodegeneration in a large North African cohort.

Mutations in the PLA2G6 gene are causative of PLA2G6‐associated neurodegeneration (PLAN), a spectrum of neurodegenerative conditions including infantile, childhood and adult onset forms.

A solid quality-control analysis of AB SOLiD short-read sequencing data

Next generation sequencers have greatly improved our ability to mine polymorphisms and mutations out of entire (or portions of) genomes. The reliability of their outputs, though, showed to be very related to the sequencing chemistry and to deeply affect the quality of the downstream analyses. We focus here on the two-base color code chemistry of AB SOLiD sequencers and propose a comprehensive quality control methodological and software pipeline. We used existing and custom tools to detect and purge short-reads of some common flaws due to sequencing errors and chemical hitches. We apply them to a cohort of SOLiD 4 runs and measure their joint efficacy in terms of the resulting ability to detect the greatest possible number of true variants.

High-confidence assessment of functional impact of human mitochondrial non-synonymous genome variations by APOGEE

24,189 are all the possible non-synonymous amino acid changes potentially affecting the human mitochondrial DNA. Only a tiny subset was functionally evaluated with certainty so far, while the pathogenicity of the vast majority was only assessed in-silico by software predictors. Since these tools proved to be rather incongruent, we have designed and implemented APOGEE, a machine-learning algorithm that outperforms all existing prediction methods in estimating the harmfulness of mitochondrial non-synonymous genome variations. We provide a detailed description of the underlying algorithm, of the selected and manually curated training and test sets of variants, as well as of its classification ability.

Stepwise analysis of MIR9 loci identifies miR-9-5p to be involved in Oestrogen regulated pathways in breast cancer patients

miR-9 was initially identified as an epigenetically regulated miRNA in tumours, but inconsistent findings have been reported so far. We analysed the expression of miR-9-5p, miR-9-3p, pri-miRs and MIR9 promoters methylation status in 131 breast cancer cases and 12 normal breast tissues (NBTs). The expression of both mature miRs was increased in tumours as compared to NBTs (P < 0.001) and negatively correlated with ER protein expression (P = 0.005 and P = 0.003, for miR-9-3p and miR-9-5p respectively). In addition, miR-9-5p showed a significant negative correlation with PgR (P = 0.002). Consistently, miR-9-5p and miR-9 3p were differentially expressed in the breast cancer subgroups identified by ER and PgR expression and HER2 amplification. No significant correlation between promoter methylation and pri-miRNAs expressions was found either in tumours or in NBTs. In the Luminal breast cancer subtype the expression of miR-9-5p was associated with a worse prognosis in both univariable and multivariable analyses. Ingenuity Pathway Analysis exploring the putative interactions among miR-9-5p/miR-9-3p, ER and PgR upstream and downstream regulators suggested a regulatory loop by which miR-9-5p but not miR-9-3p is induced by steroid hormone receptor and acts within hormone-receptor regulated pathways.

Identification of p53-target genes in Danio rerio.

To orchestrate the genomic response to cellular stress signals, p53 recognizes and binds to DNA containing specific and well-characterized p53-responsive elements (REs). Differences in RE sequences can strongly affect the p53 transactivation capacity and occur even between closely related species. Therefore, the identification and characterization of a species-specific p53 Binding sistes (BS) consensus sequence and of the associated target genes may help to provide new insights into the evolution of the p53 regulatory networks across different species. Although p53 functions were studied in a wide range of species, little is known about the p53-mediated transcriptional signature in Danio rerio. Here, we designed and biochemically validated a computational approach to identify novel p53 target genes in Danio rerio genome. Screening all the Danio rerio genome by pattern-matching-based analysis, we found p53 RE-like patterns proximal to 979 annotated Danio rerio genes. Prioritization analysis identified a subset of 134 candidate pattern-related genes, 31 of which have been investigated in further biochemical assays. Our study identified runx1, axin1, traf4a, hspa8, col4a5, necab2, and dnajc9 genes as novel direct p53 targets and 12 additional p53-controlled genes in Danio rerio genome. The proposed combinatorial approach resulted to be highly sensitive and robust for identifying new p53 target genes also in additional animal species.

A Broad Overview of Computational Methods for Predicting the Pathophysiological Effects of Non-synonymous Variants.

Next-generation sequencing has provided extraordinary opportunities to investigate the massive human genetic variability. It helped identifying several kinds of genomic mismatches from the wild-type reference genome sequences and to explain the onset of several pathogenic phenotypes and diseases susceptibility. In this context, distinguishing pathogenic from functionally neutral amino acid changes turns out to be a task as useful as complex, expensive, and time-consuming.Here, we present an exhaustive and up-to-dated survey of the algorithms and software packages conceived for the estimation of the putative pathogenicity of mutations, along with a description of the most popular mutation datasets that these tools used as training sets. Finally, we present and describe software for the prediction of cancer-related mutations.

The Biological Clock and the Molecular Basis of Lysosomal Storage Diseases

The lysosomal storage disorders encompass nearly fifty diseases provoked by lack or deficiency of enzymes essential for the breakdown of complex molecules and hallmarked by accumulation in the lysosomes of metabolic residues. Histochemistry and cytochemistry studies evidenced patterns of circadian variation of the lysosomal marker enzymes, suggesting that lysosomal function oscillates rhythmically during the 24-h day. The circadian rhythmicity of cellular processes is driven by the biological clock ticking through transcriptional/translational feedback loops hardwired by circadian genes and proteins. Malfunction of the molecular clockwork may provoke severe deregulation of downstream gene expression regulating a complex array of cellular functions leading to anatomical and functional changes. In this review we highlight that all the genes mutated in lysosomal storage disorders encode circadian transcripts suggesting a direct participation of the biological clock in the pathophysiological mechanisms underlying cellular and tissue derangements hallmarking these hereditary diseases. The 24-h periodicity of oscillation of gene transcription and translation could lead in physiological conditions to circadian rhythmicity of fluctuation of enzyme levels and activity, so that gene transfer could be envisaged to reproduce 24-h periodicity of variation of enzymatic dynamics and circadian rhythmicity could have an impact on the schedule of enzyme replacement therapy.