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.

Biallelic mutations in DYNC2LI1 are a rare cause of Ellis-van Creveld syndrome.

Ellis-van Creveld syndrome (EvC) is a chondral and ectodermal dysplasia caused by biallelic mutations in the EVC, EVC2 and WDR35 genes. A proportion of cases with clinical diagnosis of EvC, however, do not carry mutations in these genes. To identify the genetic cause of EvC in a cohort of mutation-negative patients, exome sequencing was undertaken in a family with 3 affected members, and mutation scanning of a panel of clinically and functionally relevant genes was performed in 24 additional subjects with features fitting/overlapping EvC. Compound heterozygosity for the c.2T>C (p.Met1?) and c.662C>T (p.Thr221Ile) variants in DYNC2LI1, which encodes a component of the intraflagellar transport-related dynein-2 complex previously found mutated in other short-rib thoracic dysplasias, was identified in the 3 affected members of the first family. Targeted resequencing detected compound heterozygosity for the same missense variant and a truncating change (p.Val141*) in 2 siblings with EvC from a second family, while a newborn with a more severe phenotype carried 2 DYNC2LI1 truncating variants. Our findings indicate that DYNC2LI1 mutations are associated with a wider clinical spectrum than previously appreciated, including EvC, with the severity of the phenotype likely depending on the extent of defective DYNC2LI1 function.

Mutations in the IRBIT domain of ITPR1 are a frequent cause of autosomal dominant nonprogressive congenital ataxia.

Congenital ataxias are nonprogressive neurological disorders characterized by neonatal hypotonia, developmental delay and ataxia, variably associated with intellectual disability and other neurological or extraneurological features. We performed trio-based whole-exome sequencing of 12 families with congenital cerebellar and/or vermis atrophy in parallel with targeted next-generation sequencing of known ataxia genes (CACNA1A, ITPR1, KCNC3, ATP2B3 and GRM1) in 12 additional patients with a similar phenotype. Novel pathological mutations of ITPR1 (inositol 1,4,5-trisphosphate receptor, type 1) were found in seven patients from four families (4/24, ∼16.8%) all localized in the IRBIT (inositol triphosphate receptor binding protein) domain which plays an essential role in the regulation of neuronal plasticity and development. Our study expands the mutational spectrum of ITPR1-related congenital ataxia and indicates that ITPR1 gene screening should be implemented in this subgroup of ataxias.

Analysis of co-receptor usage of circulating viral and proviral HIV genome quasispecies by ultra-deep pyrosequencing in patients who are candidates for CCR5 antagonist treatment.

UDPS combined with genotypic algorithms for prediction of HIV-1 co-receptor usage may provide quantitative data about the tropism of each variant present in the viral quasispecies. The aim of the present study was to assess co-receptor usage by ultra-deep pyrosequencing (UDPS), in comparison with the reference phenotypic test (Trofile), in patients who are candidates for CCR5 antagonist treatment, in both circulating and proviral HIV-1. Seventeen patients who were tested by Trofile were enrolled. UDPS of the V3 loop region was carried out on both plasma RNA and proviral DNA. Genotypic prediction of co-receptor usage was established by position-specific score matrices (PSSM) and confirmed, in discordant cases, with geno2pheno. Genetic heterogeneity of the RNA and DNA quasispecies was assessed as well. A total of 196,729 V3 sequences were considered (mean coverage per site, 6346). Concordance between phenotypic test and UDPS with PSSM was 0.82. Geno2pheno results were in line with those obtained with PSSM. Proviral quasispecies were more heterogeneous than those found in circulating HIV. In most patients eligible for CCR5 antagonist treatment, X4 variants were detected in proviral DNA, ranging from 1.0% to 52.7%. UDPS combined with genotypic algorithms for co-receptor usage prediction highlighted the presence of minority variants, with a discordant tropism with respect to the predominant population, in both circulating viral and proviral HIV. In most patients treated with Maraviroc the virological response was independent of the presence of X4 in proviral DNA. The clinical impact of minority X4 variants present in patients who are candidates for anti-CCR5 antagonists remains a crucial point to be addressed.

Comparison of real-time PCR methods for measurement of HIV-1 proviral DNA.

The use of HIV-1 DNA quantitation in cellular reservoirs to predict disease progression and treatment outcome in infected patients is hampered by the lack of standardization among the available methods. In the present study, real-time PCR methods used commonly for HIV-1 proviral DNA evaluation were compared, showing strong differences in the results, probably as a consequence of genome variability in the target regions. Standardization of HIV-1 proviral DNA quantitation assays is needed for use in clinical management of patients with HIV-1.

Role of the linker domain and the 203-214 N-terminal residues in the human topoisomerase I DNA complex dynamics.

The influence of the N-terminal residues 203-214 and the linker domain on motions in the human topoisomerase I-DNA complex has been investigated by comparing the molecular dynamics simulations of the system with (topo70) or without (topo58/6.3) these regions. Topo58/6.3 is found to fluctuate more than topo70, indicating that the presence of the N-terminal residues and the linker domain dampen the core and C-terminal fluctuations. The simulations also show that residues 203-207 and the linker domain participate in a network of correlated movements with key regions of the enzyme, involved in the human topoisomerase I catalytic cycle, providing a structural-dynamical explanation for the better DNA relaxation activity of topo70 when compared to topo58/6.3. The data have been examined in relation to a wealth of biochemical, site-directed mutagenesis and crystallographic data on human topoisomerase I. The simulations finally show the occurrence of a network of direct and water mediated hydrogen bonds in the proximity of the active site, and the presence of a water molecule in the appropriate position to accept a proton from the catalytic Tyr-723 residue, suggesting that water molecules have an important role in the stabilization and function of this enzyme.