Ion channel mutations are an important cause of rare Mendelian disorders affecting brain, heart, and other tissues. We performed parallel exome sequencing of 237 channel genes in a well characterized human sample, comparing variant profiles of unaffected individuals to those with the most common neuronal excitability disorder, sporadic idiopathic epilepsy. Rare missense variation in known Mendelian disease genes is prevalent in both groups at similar complexity, revealing that even deleterious ion channel mutations confer uncertain risk to an individual depending on the other variants with which they are combined. Our findings indicate that variant discovery via large scale sequencing efforts is only a first step in illuminating the complex allelic architecture underlying personal disease risk. We propose that in silico modeling of channel variation in realistic cell and network models will be crucial to future strategies assessing mutation profile pathogenicity and drug response in individuals with a broad spectrum of excitability disorders.

Background

Next-generation DNA sequencing is opening new avenues for genetic association studies in common diseases that, like deep vein thrombosis (DVT), have a strong genetic predisposition still largely unexplained by currently identified risk variants. In order to develop sequencing and analytical pipelines for the application of next-generation sequencing to complex diseases, we conducted a pilot study sequencing the coding area of 186 hemostatic/proinflammatory genes in 10 Italian cases of idiopathic DVT and 12 healthy controls.

Results

A molecular-barcoding strategy was used to multiplex DNA target capture and sequencing, while retaining individual sequence information. Genomic libraries with barcode sequence-tags were pooled (in pools of 8 or 16 samples) and enriched for target DNA sequences. Sequencing was performed on ABI SOLiD-4 platforms. We produced > 12 gigabases of raw sequence data to sequence at high coverage (average: 42X) the 700-kilobase target area in 22 individuals. A total of 1876 high-quality genetic variants were identified (1778 single nucleotide substitutions and 98 insertions/deletions). Annotation on databases of genetic variation and human disease mutations revealed several novel, potentially deleterious mutations. We tested 576 common variants in a case-control association analysis, carrying the top-5 associations over to replication in up to 719 DVT cases and 719 controls. We also conducted an analysis of the burden of nonsynonymous variants in coagulation factor and anticoagulant genes. We found an excess of rare missense mutations in anticoagulant genes in DVT cases compared to controls and an association for a missense polymorphism of FGA (rs6050; p = 1.9 × 10-5, OR 1.45; 95% CI, 1.22-1.72; after replication in > 1400 individuals).

Conclusions

We implemented a barcode-based strategy to efficiently multiplex sequencing of hundreds of candidate genes in several individuals. In the relatively small dataset of our pilot study we were able to identify bona fide associations with DVT. Our study illustrates the potential of next-generation sequencing for the discovery of genetic variation predisposing to complex diseases.

Clinical cancer genetic susceptibility analysis typically proceeds sequentially beginning with the most likely causative gene. The process is time consuming and the yield is low particularly for families with unusual patterns of cancer. We determined the results of in parallel mutation analysis of a large cancer-associated gene panel. We performed deletion analysis and sequenced the coding regions of 45 genes (8 oncogenes and 37 tumor suppressor or DNA repair genes) in 48 childhood cancer patients who also (1) were diagnosed with a second malignancy under age 30, (2) have a sibling diagnosed with cancer under age 30 and/or (3) have a major congenital anomaly or developmental delay. Deleterious mutations were identified in 6 of 48 (13%) families, 4 of which met the sibling criteria. Mutations were identified in genes previously implicated in both dominant and recessive childhood syndromes including SMARCB1, PMS2, and TP53. No pathogenic deletions were identified. This approach has provided efficient identification of childhood cancer susceptibility mutations and will have greater utility as additional cancer susceptibility genes are identified. Integrating parallel analysis of large gene panels into clinical testing will speed results and increase diagnostic yield. The failure to detect mutations in 87% of families highlights that a number of childhood cancer susceptibility genes remain to be discovered.

Ciliary dysfunction leads to a broad range of overlapping phenotypes, termed collectively as ciliopathies. This grouping is underscored by genetic overlap, where causal genes can also contribute modifying alleles to clinically distinct disorders. Here we show that mutations in TTC21B/IFT139, encoding a retrograde intraflagellar transport (IFT) protein, cause both isolated nephronophthisis (NPHP) and syndromic Jeune Asphyxiating Thoracic Dystrophy (JATD). Moreover, although systematic medical resequencing of a large, clinically diverse ciliopathy cohort and matched controls showed a similar frequency of rare changes, in vivo and in vitro evaluations unmasked a significant enrichment of pathogenic alleles in cases, suggesting that TTC21B contributes pathogenic alleles to ∼5% of ciliopathy patients. Our data illustrate how genetic lesions can be both causally associated with diverse ciliopathies, as well as interact in trans with other disease-causing genes, and highlight how saturated resequencing followed by functional analysis of all variants informs the genetic architecture of disorders.

Autism spectrum disorders (ASDs) are a heterogeneous group of neuro-developmental disorders. While significant progress has been made in the identification of genes and copy number variants associated with syndromic autism, little is known to date about the etiology of idiopathic non-syndromic autism. Sanger sequencing of 21 known autism susceptibility genes in 339 individuals with high-functioning, idiopathic ASD revealed de novo mutations in at least one of these genes in 6 of 339 probands (1.8%). Additionally, multiple events of oligogenic heterozygosity were seen, affecting 23 of 339 probands (6.8%). Screening of a control population for novel coding variants in CACNA1C, CDKL5, HOXA1, SHANK3, TSC1, TSC2 and UBE3A by the same sequencing technology revealed that controls were carriers of oligogenic heterozygous events at significantly (P < 0.01) lower rate, suggesting oligogenic heterozygosity as a new potential mechanism in the pathogenesis of ASDs.

Despite rapid advances in disease gene identification, the predictive power of the genotype remains limited, in part due to poorly understood effects of second-site modifiers. Here we demonstrate that a polymorphic coding variant of RPGRIP1L (retinitis pigmentosa GTPase regulator-interacting protein-1 like), a ciliary gene mutated in Meckel-Gruber (MKS) and Joubert (JBTS) syndromes, is associated with the development of retinal degeneration in patients with ciliopathies caused by mutations in other genes. As part of our resequencing efforts of the ciliary proteome, we identified several putative loss of function RPGRIP1L mutations, including one common variant, A229T. Multiple genetic lines of evidence showed this allele to be associated with photoreceptor loss in ciliopathies. Moreover, we show that RPGRIP1L interacts biochemically with RPGR, loss of which causes retinal degeneration, and that the 229T-encoded protein significantly compromises this interaction. Our data represent an example of modification of a discrete phenotype of syndromic disease and highlight the importance of a multifaceted approach for the discovery of modifier alleles of intermediate frequency and effect.

Determining the genetic basis of cancer requires comprehensive analyses of large collections of histopathologically well-classified primary tumours. Here we report the results of a collaborative study to discover somatic mutations in 188 human lung adenocarcinomas. DNA sequencing of 623 genes with known or potential relationships to cancer revealed more than 1,000 somatic mutations across the samples. Our analysis identified 26 genes that are mutated at significantly high frequencies and thus are probably involved in carcinogenesis. The frequently mutated genes include tyrosine kinases, among them the EGFR homologue ERBB4; multiple ephrin receptor genes, notably EPHA3; vascular endothelial growth factor receptor KDR; and NTRK genes. These data provide evidence of somatic mutations in primary lung adenocarcinoma for several tumour suppressor genes involved in other cancers—including NF1, APC, RB1 and ATM—and for sequence changes in PTPRD as well as the frequently deleted gene LRP1B. The observed mutational profiles correlate with clinical features, smoking status and DNA repair defects. These results are reinforced by data integration including single nucleotide polymorphism array and gene expression array. Our findings shed further light on several important signalling pathways involved in lung adenocarcinoma, and suggest new molecular targets for treatment.

Pulmonary alveolar proteinosis (PAP) is a rare lung disorder in which surfactant-derived lipoproteins accumulate excessively within pulmonary alveoli, causing severe respiratory distress. The importance of granulocyte/macrophage colony-stimulating factor (GM-CSF) in the pathogenesis of PAP has been confirmed in humans and mice, wherein GM-CSF signaling is required for pulmonary alveolar macrophage catabolism of surfactant. PAP is caused by disruption of GM-CSF signaling in these cells, and is usually caused by neutralizing autoantibodies to GM-CSF or is secondary to other underlying diseases. Rarely, genetic defects in surfactant proteins or the common β chain for the GM-CSF receptor (GM-CSFR) are causal. Using a combination of cellular, molecular, and genomic approaches, we provide the first evidence that PAP can result from a genetic deficiency of the GM-CSFR α chain, encoded in the X-chromosome pseudoautosomal region 1.

Accurately determining the distribution of rare variants is an important goal of human genetics, but resequencing of a sample large enough for this purpose has been unfeasible until now. Here, we applied Sanger sequencing of genomic PCR amplicons to resequence the diabetes-associated genes KCNJ11 and HHEX in 13,715 people (10,422 European Americans and 3,293 African Americans) and validated amplicons potentially harbouring rare variants using 454 pyrosequencing. We observed far more variation (expected variant-site count ∼578) than would have been predicted on the basis of earlier surveys, which could only capture the distribution of common variants. By comparison with earlier estimates based on common variants, our model shows a clear genetic signal of accelerating population growth, suggesting that humanity harbours a myriad of rare, deleterious variants, and that disease risk and the burden of disease in contemporary populations may be heavily influenced by the distribution of rare variants.

To fully catalogue rare genetic variation in humans, many samples need to be examined. In this study, Coventry et al. resequenced two genes, KCNJ11 and HHEX, in 13,715 humans, and concluded that most of the sequence variation arose recently and that variation is greater than expected.