VarGenius executes cohort-level DNA-seq variant calling and annotation and allows to manage the resulting data through a PostgreSQL database.

BACKGROUND: Targeted resequencing has become the most used and cost-effective approach for identifying causative mutations of Mendelian diseases both for diagnostics and research purposes. Due to very rapid technological progress, NGS laboratories are expanding their capabilities to address the increasing number of analyses. Several open source tools are available to build a generic variant calling pipeline, but a tool able to simultaneously execute multiple analyses, organize, and categorize the samples is still missing. RESULTS: Here we describe VarGenius, a Linux based command line software able to execute customizable pipelines for the analysis of multiple targeted resequencing data using parallel computing. VarGenius provides a database to store the output of the analysis (calling quality statistics, variant annotations, internal allelic variant frequencies) and sample information (personal data, genotypes, phenotypes). VarGenius can also perform the "joint analysis" of hundreds of samples with a single command, drastically reducing the time for the configuration and execution of the analysis. VarGenius executes the standard pipeline of the Genome Analysis Tool-Kit (GATK) best practices (GBP) for germinal variant calling, annotates the variants using Annovar, and generates a user-friendly output displaying the results through a web page. VarGenius has been tested on a parallel computing cluster with 52 machines with 120GB of RAM each. Under this configuration, a 50 M whole exome sequencing (WES) analysis for a family was executed in about 7 h (trio or quartet); a joint analysis of 30 WES in about 24 h and the parallel analysis of 34 single samples from a 1 M panel in about 2 h. CONCLUSIONS: We developed VarGenius, a "master" tool that faces the increasing demand of heterogeneous NGS analyses and allows maximum flexibility for downstream analyses. It paves the way to a different kind of analysis, centered on cohorts rather than on singleton. Patient and variant information are stored into the database and any output file can be accessed programmatically. VarGenius can be used for routine analyses by biomedical researchers with basic Linux skills providing additional flexibility for computational biologists to develop their own algorithms for the comparison and analysis of data. The software is freely available at: https://github.com/frankMusacchia/VarGenius.

PML-RARA-associated cooperating mutations belong to a transcriptional network that is deregulated in myeloid leukemias.

It has been shown that individual acute myeloid leukemia (AML) patients are characterized by one of few initiating DNA mutations and 5-10 cooperating mutations not yet defined among hundreds identified by massive sequencing of AML genomes. We report an in vivo insertional-mutagenesis screen for genes cooperating with one AML initiating mutations (PML-RARA, oncogene of acute promyelocytic leukemia, APL), which allowed identification of hundreds of genetic cooperators. The cooperators are mutated at low frequency in APL or AML patients but are always abnormally expressed in a cohort of 182 APLs and AMLs analyzed. These deregulations appear non-randomly distributed and present in all samples, regardless of their associated genomic mutations. Reverse-engineering approaches showed that these cooperators belong to a single transcriptional gene network, enriched in genes mutated in AMLs, where perturbation of single genes modifies expression of others. Their gene-ontology analysis showed enrichment of genes directly involved in cell proliferation control. Therefore, the pool of PML-RARA cooperating mutations appears large and heterogeneous, but functionally equivalent and deregulated in the majority of APLs and AMLs. Our data suggest that the high heterogeneity of DNA mutations in APLs and AMLs can be reduced to patterns of gene expression deregulation of a single 'mutated' gene network.

Expanding the phenotype of RTTN variations: a new family with primary microcephaly, severe growth failure, brain malformations and dermatitis.

Primary autosomal recessive microcephaly (MCPH) is a developmental disorder characterized by prenatal onset of abnormal brain growth. MCPH occurs both alone and as part of a broad range of neurodevelopmental syndromes with or without cortical malformations and growth retardation. Here we report a consanguineous Moroccan family with two siblings affected by severe primary microcephaly, failure to thrive, congenital dermatitis and severe developmental delay. Brain magnetic resonance imaging showed lissencephaly of frontal lobes and periventricular heterotopia of the gray matter. We performed both Comparative Genomic Hybridization array and whole exome sequencing (WES) analyses of the kindred. No quantitative defects were detected. However, WES identified a new homozygous missense variation in the penultimate nucleotide of exon 23 of RTTN gene (c.2953A>G;pArg985Gly). cDNA sequencing revealed two abnormal spliced products, one lacking only exon 23 and the other lacking exons 22 and 23 (out-of-frame). RTTN is a protein involved in cilia structure and function. Homozygous mutations in RTTN gene have been described in bilateral diffuse isolated polymicrogyria and, more recently, in microcephalic primordial dwarfism (PD). We found a novel homozygous mutation in RTTN associated with microcephalic PD as well as complex brain malformations and congenital dermatitis, thus expanding the phenotypic spectrum of both RTTN-associated diseases and ciliary dysfunction.

Genetic association of ARHGAP21 gene variant with mandibular prognathism.

Mandibular prognathism (MP) is a recognizable phenotype associated with dentoskeletal class III malocclusion. MP is a complex genetic trait, although familial recurrence also suggests the contribution of single inherited variations. To date, the genetic causes of MP have been investigated using linkage analysis or association studies in pooled families. Here for the first time, next-generation sequencing was used to study a single family with a large number of MP-affected members and to identify MP-related candidate genes. A 6-generation kindred with MP segregating as an autosomal dominant character was recruited. To identify family members affected by MP, a standard cephalometric procedure was used. In 5 MP subjects separated by the largest number of meioses, whole-exome sequencing was performed. Five promising missense gene variants (BMP3, ANXA2, FLNB, HOXA2, and ARHGAP21) associated with MP were selected and genotyped in most other family members. In this family, MP seemed to consist of 2 distinct genetic branches. Interestingly, the Gly1121Ser variant in the ARHGAP21 gene was found to be shared by all MP individuals in the larger branch of the family with nearly complete penetrance. This variant is rare in the Caucasian population (frequency 0.00034) and is predicted as damaging by all bioinformatic algorithms. ARHGAP21 protein strengthens cell-cell adhesions and may be regulated by bone morphogenetic factors, thus influencing mandibular growth. Further studies in both animal models and human patients are required to clarify the significance of this association.