Erin L. Heinzen

Genetic variation in IL28B predicts hepatitis C treatment-induced viral clearance

Chronic infection with hepatitis C virus (HCV) affects 170 million people worldwide and is the leading cause of cirrhosis in North America. Although the recommended treatment for chronic infection involves a 48-week course of peginterferon-α-2b (PegIFN-α-2b) or -α-2a (PegIFN-α-2a) combined with ribavirin (RBV), it is well known that many patients will not be cured by treatment, and that patients of European ancestry have a significantly higher probability of being cured than patients of African ancestry. In addition to limited efficacy, treatment is often poorly tolerated because of side effects that prevent some patients from completing therapy. For these reasons, identification of the determinants of response to treatment is a high priority. Here we report that a genetic polymorphism near the IL28B gene, encoding interferon-λ-3 (IFN-λ-3), is associated with an approximately twofold change in response to treatment, both among patients of European ancestry (P = 1.06 × 10-25) and African-Americans (P = 2.06 × 10-3). Because the genotype leading to better response is in substantially greater frequency in European than African populations, this genetic polymorphism also explains approximately half of the difference in response rates between African-Americans and patients of European ancestry.

De novo mutations in epileptic encephalopathies

Epileptic encephalopathies are a devastating group of severe childhood epilepsy disorders for which the cause is often unknown. Here we report a screen for de novo mutations in patients with two classical epileptic encephalopathies: infantile spasms (n = 149) and Lennox–Gastaut syndrome (n = 115). We sequenced the exomes of 264 probands, and their parents, and confirmed 329 de novo mutations. A likelihood analysis showed a significant excess of de novo mutations in the ∼4,000 genes that are the most intolerant to functional genetic variation in the human population (P = 2.9 × 10−3). Among these are GABRB3, with de novo mutations in four patients, and ALG13, with the same de novo mutation in two patients; both genes show clear statistical evidence of association with epileptic encephalopathy. Given the relevant site-specific mutation rates, the probabilities of these outcomes occurring by chance are P = 4.1 × 10−10 and P = 7.8 × 10−12, respectively. Other genes with de novo mutations in this cohort include CACNA1A, CHD2, FLNA, GABRA1, GRIN1, GRIN2B, HNRNPU, IQSEC2, MTOR and NEDD4L. Finally, we show that the de novo mutations observed are enriched in specific gene sets including genes regulated by the fragile X protein (P < 10−8), as has been reported previously for autism spectrum disorders.

HLA-A*3101 and Carbamazepine-Induced Hypersensitivity Reactions in Europeans

BACKGROUND Carbamazepine causes various forms of hypersensitivity reactions, ranging from maculopapular exanthema to severe blistering reactions. The HLA-B*1502 allele has been shown to be strongly correlated with carbamazepine-induced Stevens-Johnson syndrome and toxic epidermal necrolysis (SJS-TEN) in the Han Chinese and other Asian populations but not in European populations. METHODS We performed a genomewide association study of samples obtained from 22 subjects with carbamazepine-induced hypersensitivity syndrome, 43 subjects with carbamazepine-induced maculopapular exanthema, and 3987 control subjects, all of European descent. We tested for an association between disease and HLA alleles through proxy single-nucleotide polymorphisms and imputation, confirming associations by high-resolution sequence-based HLA typing. We replicated the associations in samples from 145 subjects with carbamazepine-induced hypersensitivity reactions. RESULTS The HLA-A*3101 allele, which has a prevalence of 2 to 5% in Northern European populations, was significantly associated with the hypersensitivity syndrome (P=3.5×10(-8)). An independent genomewide association study of samples from subjects with maculopapular exanthema also showed an association with the HLA-A*3101 allele (P=1.1×10(-6)). Follow-up genotyping confirmed the variant as a risk factor for the hypersensitivity syndrome (odds ratio, 12.41; 95% confidence interval [CI], 1.27 to 121.03), maculopapular exanthema (odds ratio, 8.33; 95% CI, 3.59 to 19.36), and SJS-TEN (odds ratio, 25.93; 95% CI, 4.93 to 116.18). CONCLUSIONS The presence of the HLA-A*3101 allele was associated with carbamazepine-induced hypersensitivity reactions among subjects of Northern European ancestry. The presence of the allele increased the risk from 5.0% to 26.0%, whereas its absence reduced the risk from 5.0% to 3.8%. (Funded by the U.K. Department of Health and others.).

Genic Intolerance to Functional Variation and the Interpretation of Personal Genomes

A central challenge in interpreting personal genomes is determining which mutations most likely influence disease. Although progress has been made in scoring the functional impact of individual mutations, the characteristics of the genes in which those mutations are found remain largely unexplored. For example, genes known to carry few common functional variants in healthy individuals may be judged more likely to cause certain kinds of disease than genes known to carry many such variants. Until now, however, it has not been possible to develop a quantitative assessment of how well genes tolerate functional genetic variation on a genome-wide scale. Here we describe an effort that uses sequence data from 6503 whole exome sequences made available by the NHLBI Exome Sequencing Project (ESP). Specifically, we develop an intolerance scoring system that assesses whether genes have relatively more or less functional genetic variation than expected based on the apparently neutral variation found in the gene. To illustrate the utility of this intolerance score, we show that genes responsible for Mendelian diseases are significantly more intolerant to functional genetic variation than genes that do not cause any known disease, but with striking variation in intolerance among genes causing different classes of genetic disease. We conclude by showing that use of an intolerance ranking system can aid in interpreting personal genomes and identifying pathogenic mutations.

A Genome-Wide Investigation of SNPs and CNVs in Schizophrenia

We report a genome-wide assessment of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) in schizophrenia. We investigated SNPs using 871 patients and 863 controls, following up the top hits in four independent cohorts comprising 1,460 patients and 12,995 controls, all of European origin. We found no genome-wide significant associations, nor could we provide support for any previously reported candidate gene or genome-wide associations. We went on to examine CNVs using a subset of 1,013 cases and 1,084 controls of European ancestry, and a further set of 60 cases and 64 controls of African ancestry. We found that eight cases and zero controls carried deletions greater than 2 Mb, of which two, at 8p22 and 16p13.11-p12.4, are newly reported here. A further evaluation of 1,378 controls identified no deletions greater than 2 Mb, suggesting a high prior probability of disease involvement when such deletions are observed in cases. We also provide further evidence for some smaller, previously reported, schizophrenia-associated CNVs, such as those in NRXN1 and APBA2. We could not provide strong support for the hypothesis that schizophrenia patients have a significantly greater “load” of large (>100 kb), rare CNVs, nor could we find common CNVs that associate with schizophrenia. Finally, we did not provide support for the suggestion that schizophrenia-associated CNVs may preferentially disrupt genes in neurodevelopmental pathways. Collectively, these analyses provide the first integrated study of SNPs and CNVs in schizophrenia and support the emerging view that rare deleterious variants may be more important in schizophrenia predisposition than common polymorphisms. While our analyses do not suggest that implicated CNVs impinge on particular key pathways, we do support the contribution of specific genomic regions in schizophrenia, presumably due to recurrent mutation. On balance, these data suggest that very few schizophrenia patients share identical genomic causation, potentially complicating efforts to personalize treatment regimens.

Tissue-Specific Genetic Control of Splicing: Implications for the Study of Complex Traits

Numerous genome-wide screens for polymorphisms that influence gene expression have provided key insights into the genetic control of transcription. Despite this work, the relevance of specific polymorphisms to in vivo expression and splicing remains unclear. We carried out the first genome-wide screen, to our knowledge, for SNPs that associate with alternative splicing and gene expression in human primary cells, evaluating 93 autopsy-collected cortical brain tissue samples with no defined neuropsychiatric condition and 80 peripheral blood mononucleated cell samples collected from living healthy donors. We identified 23 high confidence associations with total expression and 80 with alternative splicing as reflected by expression levels of specific exons. Fewer than 50% of the implicated SNPs however show effects in both tissue types, reflecting strong evidence for distinct genetic control of splicing and expression in the two tissue types. The data generated here also suggest the possibility that splicing effects may be responsible for up to 13 out of 84 reported genome-wide significant associations with human traits. These results emphasize the importance of establishing a database of polymorphisms affecting splicing and expression in primary tissue types and suggest that splicing effects may be of more phenotypic significance than overall gene expression changes.

Somatic Activation of AKT3 Causes Hemispheric Developmental Brain Malformations

Hemimegalencephaly (HMG) is a developmental brain disorder characterized by an enlarged, malformed cerebral hemisphere, typically causing epilepsy that requires surgical resection. We studied resected HMG tissue to test whether the condition might reflect somatic mutations affecting genes critical to brain development. We found that two out of eight HMG samples showed trisomy of chromosome 1q, which encompasses many genes, including AKT3, a gene known to regulate brain size. A third case showed a known activating mutation in AKT3 (c.49G→A, creating p.E17K) that was not present in the patients blood cells. Remarkably, the E17K mutation in AKT3 is exactly paralogous to E17K mutations in AKT1 and AKT2 recently discovered in somatic overgrowth syndromes. We show that AKT3 is the most abundant AKT paralog in the brain during neurogenesis and that phosphorylated AKT is abundant in cortical progenitor cells. Our data suggest that somatic mutations limited to the brain could represent an important cause of complex neurogenetic disease.

De novo mutations in ATP1A3 cause alternating hemiplegia of childhood

Alternating hemiplegia of childhood (AHC) is a rare, severe neurodevelopmental syndrome characterized by recurrent hemiplegic episodes and distinct neurological manifestations. AHC is usually a sporadic disorder and has unknown etiology. We used exome sequencing of seven patients with AHC and their unaffected parents to identify de novo nonsynonymous mutations in ATP1A3 in all seven individuals. In a subsequent sequence analysis of ATP1A3 in 98 other patients with AHC, we found that ATP1A3 mutations were likely to be responsible for at least 74% of the cases; we also identified one inherited mutation in a case of familial AHC. Notably, most AHC cases are caused by one of seven recurrent ATP1A3 mutations, one of which was observed in 36 patients. Unlike ATP1A3 mutations that cause rapid-onset dystonia-parkinsonism, AHC-causing mutations in this gene caused consistent reductions in ATPase activity without affecting the level of protein expression. This work identifies de novo ATP1A3 mutations as the primary cause of AHC and offers insight into disease pathophysiology by expanding the spectrum of phenotypes associated with mutations in ATP1A3.

Rare deletions at 16p13.11 predispose to a diverse spectrum of sporadic epilepsy syndromes.

Deletions at 16p13.11 are associated with schizophrenia, mental retardation, and most recently idiopathic generalized epilepsy. To evaluate the role of 16p13.11 deletions, as well as other structural variation, in epilepsy disorders, we used genome-wide screens to identify copy number variation in 3812 patients with a diverse spectrum of epilepsy syndromes and in 1299 neurologically-normal controls. Large deletions (> 100 kb) at 16p13.11 were observed in 23 patients, whereas no control had a deletion greater than 16 kb. Patients, even those with identically sized 16p13.11 deletions, presented with highly variable epilepsy phenotypes. For a subset of patients with a 16p13.11 deletion, we show a consistent reduction of expression for included genes, suggesting that haploinsufficiency might contribute to pathogenicity. We also investigated another possible mechanism of pathogenicity by using hybridization-based capture and next-generation sequencing of the homologous chromosome for ten 16p13.11-deletion patients to look for unmasked recessive mutations. Follow-up genotyping of suggestive polymorphisms failed to identify any convincing recessive-acting mutations in the homologous interval corresponding to the deletion. The observation that two of the 16p13.11 deletions were larger than 2 Mb in size led us to screen for other large deletions. We found 12 additional genomic regions harboring deletions > 2 Mb in epilepsy patients, and none in controls. Additional evaluation is needed to characterize the role of these exceedingly large, non-locus-specific deletions in epilepsy. Collectively, these data implicate 16p13.11 and possibly other large deletions as risk factors for a wide range of epilepsy disorders, and they appear to point toward haploinsufficiency as a contributor to the pathogenicity of deletions.

The Characterization of Twenty Sequenced Human Genomes

We present the analysis of twenty human genomes to evaluate the prospects for identifying rare functional variants that contribute to a phenotype of interest. We sequenced at high coverage ten “case” genomes from individuals with severe hemophilia A and ten “control” genomes. We summarize the number of genetic variants emerging from a study of this magnitude, and provide a proof of concept for the identification of rare and highly-penetrant functional variants by confirming that the cause of hemophilia A is easily recognizable in this data set. We also show that the number of novel single nucleotide variants (SNVs) discovered per genome seems to stabilize at about 144,000 new variants per genome, after the first 15 individuals have been sequenced. Finally, we find that, on average, each genome carries 165 homozygous protein-truncating or stop loss variants in genes representing a diverse set of pathways.