GenPipes: an open-source framework for distributed and scalable genomic analyses.

BACKGROUND: With the decreasing cost of sequencing and the rapid developments in genomics technologies and protocols, the need for validated bioinformatics software that enables efficient large-scale data processing is growing. FINDINGS: Here we present GenPipes, a flexible Python-based framework that facilitates the development and deployment of multi-step workflows optimized for high-performance computing clusters and the cloud. GenPipes already implements 12 validated and scalable pipelines for various genomics applications, including RNA sequencing, chromatin immunoprecipitation sequencing, DNA sequencing, methylation sequencing, Hi-C, capture Hi-C, metagenomics, and Pacific Biosciences long-read assembly. The software is available under a GPLv3 open source license and is continuously updated to follow recent advances in genomics and bioinformatics. The framework has already been configured on several servers, and a Docker image is also available to facilitate additional installations. CONCLUSIONS: GenPipes offers genomics researchers a simple method to analyze different types of data, customizable to their needs and resources, as well as the flexibility to create their own workflows.

Whole exome sequencing identifies novel predisposing genes in neural tube defects.

BACKGROUND: Neural tube defects (NTD) are among the most common defects affecting 1:1000 births. They are caused by a failure of neural tube closure during development. Their clinical presentation is diverse and dependent on the site and severity of the original defect on the embryonic axis. The etiology of NTD is multifactorial involving environmental factors and genetic variants that remain largely unknown. METHODS: We have conducted a whole exome sequencing (WES) study in five new NTD families and pooled the results with WES data from three NTD families and 43 trios that were previously investigated by our group. We analyzed the data using biased candidate gene and unbiased gene burden approaches. RESULTS: We identified four novel loss-of-function variants in three genes, MTHFR, DLC1, and ITGB1, previously associated with NTD. Notably, DLC1 carried two protein truncating variants in two independent cases. We also demonstrated an enrichment of variants in MYO1E involved in cytoskeletal remodeling. This enrichment reached borderline significance in a replication cohort supporting the association of this new candidate gene to NTD. CONCLUSION: These data provide some key insights into the pathogenic mechanisms of human NTD and demonstrate the power of next-generation sequencing in deciphering the genetics of this complex trait.

Increased missense mutation burden of Fatty Acid metabolism related genes in nunavik inuit population.

BACKGROUND: Nunavik Inuit (northern Quebec, Canada) reside along the arctic coastline where for generations their daily energy intake has mainly been derived from animal fat. Given this particular diet it has been hypothesized that natural selection would lead to population specific allele frequency differences and unique variants in genes related to fatty acid metabolism. A group of genes, namely CPT1A, CPT1B, CPT1C, CPT2, CRAT and CROT, encode for three carnitine acyltransferases that are important for the oxidation of fatty acids, a critical step in their metabolism. METHODS: Exome sequencing and SNP array genotyping were used to examine the genetic variations in the six genes encoding for the carnitine acyltransferases in 113 Nunavik Inuit individuals. RESULTS: Altogether ten missense variants were found in genes CPT1A, CPT1B, CPT1C, CPT2 and CRAT, including three novel variants and one Inuit specific variant CPT1A p.P479L (rs80356779). The latter has the highest frequency (0.955) compared to other Inuit populations. We found that by comparison to Asians or Europeans, the Nunavik Inuit have an increased mutation burden in CPT1A, CPT2 and CRAT; there is also a high level of population differentiation based on carnitine acyltransferase gene variations between Nunavik Inuit and Asians. CONCLUSION: The increased number and frequency of deleterious variants in these fatty acid metabolism genes in Nunavik Inuit may be the result of genetic adaptation to their diet and/or the extremely cold climate. In addition, the identification of these variants may help to understand some of the specific health risks of Nunavik Inuit.

Loss-of-function de novo mutations play an important role in severe human neural tube defects.

BACKGROUND: Neural tube defects (NTDs) are very common and severe birth defects that are caused by failure of neural tube closure and that have a complex aetiology. Anencephaly and spina bifida are severe NTDs that affect reproductive fitness and suggest a role for de novo mutations (DNMs) in their aetiology. METHODS: We used whole-exome sequencing in 43 sporadic cases affected with myelomeningocele or anencephaly and their unaffected parents to identify DNMs in their exomes. RESULTS: We identified 42 coding DNMs in 25 cases, of which 6 were loss of function (LoF) showing a higher rate of LoF DNM in our cohort compared with control cohorts. Notably, we identified two protein-truncating DNMs in two independent cases in SHROOM3, previously associated with NTDs only in animal models. We have demonstrated a significant enrichment of LoF DNMs in this gene in NTDs compared with the gene specific DNM rate and to the DNM rate estimated from control cohorts. We also identified one nonsense DNM in PAX3 and two potentially causative missense DNMs in GRHL3 and PTPRS. CONCLUSIONS: Our study demonstrates an important role of LoF DNMs in the development of NTDs and strongly implicates SHROOM3 in its aetiology.

De novo mutations in moderate or severe intellectual disability.

Genetics is believed to have an important role in intellectual disability (ID). Recent studies have emphasized the involvement of de novo mutations (DNMs) in ID but the extent to which they contribute to its pathogenesis and the identity of the corresponding genes remain largely unknown. Here, we report a screen for DNMs in subjects with moderate or severe ID. We sequenced the exomes of 41 probands and their parents, and confirmed 81 DNMs affecting the coding sequence or consensus splice sites (1.98 DNMs/proband). We observed a significant excess of de novo single nucleotide substitutions and loss-of-function mutations in these cases compared to control subjects, suggesting that at least a subset of these variations are pathogenic. A total of 12 likely pathogenic DNMs were identified in genes previously associated with ID (ARID1B, CHD2, FOXG1, GABRB3, GATAD2B, GRIN2B, MBD5, MED13L, SETBP1, TBR1, TCF4, WDR45), resulting in a diagnostic yield of ∼29%. We also identified 12 possibly pathogenic DNMs in genes (HNRNPU, WAC, RYR2, SET, EGR1, MYH10, EIF2C1, COL4A3BP, CHMP2A, PPP1CB, VPS4A, PPP2R2B) that have not previously been causally linked to ID. Interestingly, no case was explained by inherited mutations. Protein network analysis indicated that the products of many of these known and candidate genes interact with each other or with products of other ID-associated genes further supporting their involvement in ID. We conclude that DNMs represent a major cause of moderate or severe ID.

Increased exonic de novo mutation rate in individuals with schizophrenia.

Schizophrenia is a severe psychiatric disorder that profoundly affects cognitive, behavioral and emotional processes. The wide spectrum of symptoms and clinical variability in schizophrenia suggest a complex genetic etiology, which is consistent with the numerous loci thus far identified by linkage, copy number variation and association studies. Although schizophrenia heritability may be as high as ∼80%, the genes responsible for much of this heritability remain to be identified. Here we sequenced the exomes of 14 schizophrenia probands and their parents. We identified 15 de novo mutations (DNMs) in eight probands, which is significantly more than expected considering the previously reported DNM rate. In addition, 4 of the 15 identified DNMs are nonsense mutations, which is more than what is expected by chance. Our study supports the notion that DNMs may account for some of the heritability reported for schizophrenia while providing a list of genes possibly involved in disease pathogenesis.

A population genetic approach to mapping neurological disorder genes using deep resequencing.

Deep resequencing of functional regions in human genomes is key to identifying potentially causal rare variants for complex disorders. Here, we present the results from a large-sample resequencing (n  =  285 patients) study of candidate genes coupled with population genetics and statistical methods to identify rare variants associated with Autism Spectrum Disorder and Schizophrenia. Three genes, MAP1A, GRIN2B, and CACNA1F, were consistently identified by different methods as having significant excess of rare missense mutations in either one or both disease cohorts. In a broader context, we also found that the overall site frequency spectrum of variation in these cases is best explained by population models of both selection and complex demography rather than neutral models or models accounting for complex demography alone. Mutations in the three disease-associated genes explained much of the difference in the overall site frequency spectrum among the cases versus controls. This study demonstrates that genes associated with complex disorders can be mapped using resequencing and analytical methods with sample sizes far smaller than those required by genome-wide association studies. Additionally, our findings support the hypothesis that rare mutations account for a proportion of the phenotypic variance of these complex disorders.

Excess of de novo deleterious mutations in genes associated with glutamatergic systems in nonsyndromic intellectual disability.

Little is known about the genetics of nonsyndromic intellectual disability (NSID). We hypothesized that de novo mutations (DNMs) in synaptic genes explain an important fraction of sporadic NSID cases. In order to investigate this possibility, we sequenced 197 genes encoding glutamate receptors and a large subset of their known interacting proteins in 95 sporadic cases of NSID. We found 11 DNMs, including ten potentially deleterious mutations (three nonsense, two splicing, one frameshift, four missense) and one neutral mutation (silent) in eight different genes. Calculation of point-substitution DNM rates per functional and neutral site showed significant excess of functional DNMs compared to neutral ones. De novo truncating and/or splicing mutations in SYNGAP1, STXBP1, and SHANK3 were found in six patients and are likely to be pathogenic. De novo missense mutations were found in KIF1A, GRIN1, CACNG2, and EPB41L1. Functional studies showed that all these missense mutations affect protein function in cell culture systems, suggesting that they may be pathogenic. Sequencing these four genes in 50 additional sporadic cases of NSID identified a second DNM in GRIN1 (c.1679_1681dup/p.Ser560dup). This mutation also affects protein function, consistent with structural predictions. None of these mutations or any other DNMs were identified in these genes in 285 healthy controls. This study highlights the importance of the glutamate receptor complexes in NSID and further supports the role of DNMs in this disorder.

Resequencing of 29 candidate genes in patients with familial and sporadic amyotrophic lateral sclerosis.

OBJECTIVE: To identify novel disease-causing genes for amyotrophic lateral sclerosis (ALS). DESIGN, SETTING, AND PATIENTS: We carried out a systematic mutation screening of the entire coding regions of 29 candidate genes encoding critically important proteins for proper differentiation and development of corticospinal motor neurons in 190 patients with familial and sporadic ALS. MAIN OUTCOME MEASURES: We focused our analysis on coding variants and evaluated the distribution of nonsynonymous and synonymous variants in our cohort of patients with ALS. RESULTS: We identified 40 novel nonsynonymous variants and showed a significant excess of unique nonsynonymous variants in our cohort of patients with ALS, which suggests the presence of ALS-predisposing mutations. CONCLUSIONS: Using a multifaceted approach based on the functional prediction of missense variants, the conservation of the altered amino acid, and the cosegregation of the variants identified in familial cases, we identified several promising novel genes for ALS such as LUM and CRYM. We have also highlighted the analytical challenges of large-scale sequencing screens to detect disease-causing variants.

Direct measure of the de novo mutation rate in autism and schizophrenia cohorts.

The role of de novo mutations (DNMs) in common diseases remains largely unknown. Nonetheless, the rate of de novo deleterious mutations and the strength of selection against de novo mutations are critical to understanding the genetic architecture of a disease. Discovery of high-impact DNMs requires substantial high-resolution interrogation of partial or complete genomes of families via resequencing. We hypothesized that deleterious DNMs may play a role in cases of autism spectrum disorders (ASD) and schizophrenia (SCZ), two etiologically heterogeneous disorders with significantly reduced reproductive fitness. We present a direct measure of the de novo mutation rate (µ) and selective constraints from DNMs estimated from a deep resequencing data set generated from a large cohort of ASD and SCZ cases (n = 285) and population control individuals (n = 285) with available parental DNA. A survey of ∼430 Mb of DNA from 401 synapse-expressed genes across all cases and 25 Mb of DNA in controls found 28 candidate DNMs, 13 of which were cell line artifacts. Our calculated direct neutral mutation rate (1.36 × 10(-8)) is similar to previous indirect estimates, but we observed a significant excess of potentially deleterious DNMs in ASD and SCZ individuals. Our results emphasize the importance of DNMs as genetic mechanisms in ASD and SCZ and the limitations of using DNA from archived cell lines to identify functional variants.

De novo truncating mutation in Kinesin 17 associated with schizophrenia.

BACKGROUND: Schizophrenia (SCZ) is one of the most disabling psychiatric disorders. It is thought to be due to a complex interplay between polygenic and various environmental risk factors, although recent reports on genomic copy number variations suggest that a fraction of the cases could result from variably penetrant de novo variants. The gene encoding the synaptic motor protein kinesin 17 (KIF17) involved in glutamatergic synapse is a candidate gene for SCZ. METHODS: As part of our Synapse to Disease project, we resequenced KIF17 in a cohort of individuals with sporadic SCZ (188 subjects). Additional populations included autism spectrum disorder (142 subjects), nonsyndromic mental retardation (95 subjects), and control subjects (568 subjects). Functional validation of the human mutation was done in developing zebrafish. RESULTS: Here we report the identification of a de novo nonsense truncating mutation in one patient with SCZ, in kinesin 17, a synaptic motor protein. No de novo or truncating KIF17 mutations were found in the additional samples. We further validated the pathogenic nature of this mutation by knocking down its expression in zebrafish embryos, which resulted in a developmental defect. CONCLUSIONS: Together our findings suggest that disruption of KIF17, although rare, could result in a schizophrenia phenotype and emphasize the possible involvement of rare de novo mutations in this disorder.

De novo STXBP1 mutations in mental retardation and nonsyndromic epilepsy.

We sequenced genes coding for components of the SNARE complex (STX1A, VAMP2, SNAP25) and their regulatory proteins (STXBP1/Munc18-1, SYT1), which are essential for neurotransmission, in 95 patients with idiopathic mental retardation. We identified de novo mutations in STXBP1 (nonsense, p.R388X; splicing, c.169+1G>A) in two patients with severe mental retardation and nonsyndromic epilepsy. Reverse transcriptase polymerase chain reaction and sequencing showed that the splicing mutation creates a stop codon downstream of exon-3. No de novo or deleterious mutations in STXBP1 were found in 190 control subjects, or in 142 autistic patients. These results suggest that STXBP1 disruption is associated with autosomal dominant mental retardation and nonsyndromic epilepsy.