Evan A. Boyle

Saturation editing of genomic regions by multiplex homology-directed repair

Saturation mutagenesis—coupled to an appropriate biological assay—represents a fundamental means of achieving a high-resolution understanding of regulatory and protein-coding nucleic acid sequences of interest. However, mutagenized sequences introduced in trans on episomes or via random or “safe-harbour” integration fail to capture the native context of the endogenous chromosomal locus. This shortcoming markedly limits the interpretability of the resulting measurements of mutational impact. Here, we couple CRISPR/Cas9 RNA-guided cleavage with multiplex homology-directed repair using a complex library of donor templates to demonstrate saturation editing of genomic regions. In exon 18 of BRCA1, we replace a six-base-pair (bp) genomic region with all possible hexamers, or the full exon with all possible single nucleotide variants (SNVs), and measure strong effects on transcript abundance attributable to nonsense-mediated decay and exonic splicing elements. We similarly perform saturation genome editing of a well-conserved coding region of an essential gene, DBR1, and measure relative effects on growth that correlate with functional impact. Measurement of the functional consequences of large numbers of mutations with saturation genome editing will potentially facilitate high-resolution functional dissection of both cis-regulatory elements and trans-acting factors, as well as the interpretation of variants of uncertain significance observed in clinical sequencing.

Recurrent de novo mutations implicate novel genes underlying simplex autism risk

Autism spectrum disorder (ASD) has a strong but complex genetic component. Here we report on the resequencing of 64 candidate neurodevelopmental disorder risk genes in 5,979 individuals: 3,486 probands and 2,493 unaffected siblings. We find a strong burden of de novo point mutations for these genes and specifically implicate nine genes. These include CHD2 and SYNGAP1, genes previously reported in related disorders, and novel genes TRIP12 and PAX5. We also show that mutation carriers generally have lower IQs and enrichment for seizures. These data begin to distinguish genetically distinct subtypes of autism important for etiological classification and future therapeutics.

A germline homozygous mutation in the base-excision repair gene NTHL1 causes adenomatous polyposis and colorectal cancer

The genetic cause underlying the development of multiple colonic adenomas, the premalignant precursors of colorectal cancer (CRC), frequently remains unresolved in patients with adenomatous polyposis. Here we applied whole-exome sequencing to 51 individuals with multiple colonic adenomas from 48 families. In seven affected individuals from three unrelated families, we identified a homozygous germline nonsense mutation in the base-excision repair (BER) gene NTHL1. This mutation was exclusively found in a heterozygous state in controls (minor allele frequency of 0.0036; n = 2,329). All three families showed recessive inheritance of the adenomatous polyposis phenotype and progression to CRC in at least one member. All three affected women developed an endometrial malignancy or premalignancy. Genetic analysis of three carcinomas and five adenomas from different affected individuals showed a non-hypermutated profile enriched for cytosine-to-thymine transitions. We conclude that a homozygous loss-of-function germline mutation in the NTHL1 gene predisposes to a new subtype of BER-associated adenomatous polyposis and CRC.

Mammalian target of rapamycin pathway mutations cause hemimegalencephaly and focal cortical dysplasia

Focal malformations of cortical development, including focal cortical dysplasia (FCD) and hemimegalencephaly (HME), are important causes of intractable childhood epilepsy. Using targeted and exome sequencing on DNA from resected brain samples and nonbrain samples from 53 patients with FCD or HME, we identified pathogenic germline and mosaic mutations in multiple PI3K/AKT pathway genes in 9 patients, and a likely pathogenic variant in 1 additional patient. Our data confirm the association of DEPDC5 with sporadic FCD but also implicate this gene for the first time in HME. Our findings suggest that modulation of the mammalian target of rapamycin pathway may hold promise for malformation‐associated epilepsy. Ann Neurol 2015;77:720–725

Detection of human adaptation during the past 2000 years

Identifying genes under recent selection Evolutionary analyses aim to identify recent genetic changes that are likely to have been subject to selection. Field et al. present a method to identify such changes, the singleton density score, which they applied to over 3000 human genomes. Over the past ∼100 generations (2000 to 3000 years), Europeans are likely to have experienced selection for genetic variants, including those that affect skin and hair pigmentation, as well as height. Science, this issue p. 760 The singleton density score identifies evidence of recent selection in Europeans. Detection of recent natural selection is a challenging problem in population genetics. Here we introduce the singleton density score (SDS), a method to infer very recent changes in allele frequencies from contemporary genome sequences. Applied to data from the UK10K Project, SDS reflects allele frequency changes in the ancestors of modern Britons during the past ~2000 to 3000 years. We see strong signals of selection at lactase and the major histocompatibility complex, and in favor of blond hair and blue eyes. For polygenic adaptation, we find that recent selection for increased height has driven allele frequency shifts across most of the genome. Moreover, we identify shifts associated with other complex traits, suggesting that polygenic adaptation has played a pervasive role in shaping genotypic and phenotypic variation in modern humans.

An siRNA-based functional genomics screen for the identification of regulators of ciliogenesis and ciliopathy genes

Defects in primary cilium biogenesis underlie the ciliopathies, a growing group of genetic disorders. We describe a whole-genome siRNA-based reverse genetics screen for defects in biogenesis and/or maintenance of the primary cilium, obtaining a global resource. We identify 112 candidate ciliogenesis and ciliopathy genes, including 44 components of the ubiquitin–proteasome system, 12 G-protein-coupled receptors, and 3 pre-mRNA processing factors (PRPF6, PRPF8 and PRPF31) mutated in autosomal dominant retinitis pigmentosa. The PRPFs localize to the connecting cilium, and PRPF8- and PRPF31-mutated cells have ciliary defects. Combining the screen with exome sequencing data identified recessive mutations in PIBF1, also known as CEP90, and C21orf2, also known as LRRC76, as causes of the ciliopathies Joubert and Jeune syndromes. Biochemical approaches place C21orf2 within key ciliopathy-associated protein modules, offering an explanation for the skeletal and retinal involvement observed in individuals with C21orf2 variants. Our global, unbiased approaches provide insights into ciliogenesis complexity and identify roles for unanticipated pathways in human genetic disease.

Joubert syndrome: a model for untangling recessive disorders with extreme genetic heterogeneity

Background Joubert syndrome (JS) is a recessive neurodevelopmental disorder characterised by hypotonia, ataxia, cognitive impairment, abnormal eye movements, respiratory control disturbances and a distinctive mid-hindbrain malformation. JS demonstrates substantial phenotypic variability and genetic heterogeneity. This study provides a comprehensive view of the current genetic basis, phenotypic range and gene–phenotype associations in JS. Methods We sequenced 27 JS-associated genes in 440 affected individuals (375 families) from a cohort of 532 individuals (440 families) with JS, using molecular inversion probe-based targeted capture and next-generation sequencing. Variant pathogenicity was defined using the Combined Annotation Dependent Depletion algorithm with an optimised score cut-off. Results We identified presumed causal variants in 62% of pedigrees, including the first B9D2 mutations associated with JS. 253 different mutations in 23 genes highlight the extreme genetic heterogeneity of JS. Phenotypic analysis revealed that only 34% of individuals have a ‘pure JS’ phenotype. Retinal disease is present in 30% of individuals, renal disease in 25%, coloboma in 17%, polydactyly in 15%, liver fibrosis in 14% and encephalocele in 8%. Loss of CEP290 function is associated with retinal dystrophy, while loss of TMEM67 function is associated with liver fibrosis and coloboma, but we observe no clear-cut distinction between JS subtypes. Conclusions This work illustrates how combining advanced sequencing techniques with phenotypic data addresses extreme genetic heterogeneity to provide diagnostic and carrier testing, guide medical monitoring for progressive complications, facilitate interpretation of genome-wide sequencing results in individuals with a variety of phenotypes and enable gene-specific treatments in the future.

Association of MTOR Mutations With Developmental Brain Disorders, Including Megalencephaly, Focal Cortical Dysplasia, and Pigmentary Mosaicism

IMPORTANCE Focal cortical dysplasia (FCD), hemimegalencephaly, and megalencephaly constitute a spectrum of malformations of cortical development with shared neuropathologic features. These disorders are associated with significant childhood morbidity and mortality. OBJECTIVE To identify the underlying molecular cause of FCD, hemimegalencephaly, and diffuse megalencephaly. DESIGN, SETTING, AND PARTICIPANTS Patients with FCD, hemimegalencephaly, or megalencephaly (mean age, 11.7 years; range, 2-32 years) were recruited from Pediatric Hospital A. Meyer, the University of Hong Kong, and Seattle Childrens Research Institute from June 2012 to June 2014. Whole-exome sequencing (WES) was performed on 8 children with FCD or hemimegalencephaly using standard-depth (50-60X) sequencing in peripheral samples (blood, saliva, or skin) from the affected child and their parents and deep (150-180X) sequencing in affected brain tissue. Targeted sequencing and WES were used to screen 93 children with molecularly unexplained diffuse or focal brain overgrowth. Histopathologic and functional assays of phosphatidylinositol 3-kinase-AKT (serine/threonine kinase)-mammalian target of rapamycin (mTOR) pathway activity in resected brain tissue and cultured neurons were performed to validate mutations. MAIN OUTCOMES AND MEASURES Whole-exome sequencing and targeted sequencing identified variants associated with this spectrum of developmental brain disorders. RESULTS Low-level mosaic mutations of MTOR were identified in brain tissue in 4 children with FCD type 2a with alternative allele fractions ranging from 0.012 to 0.086. Intermediate-level mosaic mutation of MTOR (p.Thr1977Ile) was also identified in 3 unrelated children with diffuse megalencephaly and pigmentary mosaicism in skin. Finally, a constitutional de novo mutation of MTOR (p.Glu1799Lys) was identified in 3 unrelated children with diffuse megalencephaly and intellectual disability. Molecular and functional analysis in 2 children with FCD2a from whom multiple affected brain tissue samples were available revealed a mutation gradient with an epicenter in the most epileptogenic area. When expressed in cultured neurons, all MTOR mutations identified here drive constitutive activation of mTOR complex 1 and enlarged neuronal size. CONCLUSIONS AND RELEVANCE In this study, mutations of MTOR were associated with a spectrum of brain overgrowth phenotypes extending from FCD type 2a to diffuse megalencephaly, distinguished by different mutations and levels of mosaicism. These mutations may be sufficient to cause cellular hypertrophy in cultured neurons and may provide a demonstration of the pattern of mosaicism in brain and substantiate the link between mosaic mutations of MTOR and pigmentary mosaicism in skin.

MIPgen: optimized modeling and design of molecular inversion probes for targeted resequencing

UNLABELLED Molecular inversion probes (MIPs) enable cost-effective multiplex targeted gene resequencing in large cohorts. However, the design of individual MIPs is a critical parameter governing the performance of this technology with respect to capture uniformity and specificity. MIPgen is a user-friendly package that simplifies the process of designing custom MIP assays to arbitrary targets. New logistic and SVM-derived models enable in silico predictions of assay success, and assay redesign exhibits improved coverage uniformity relative to previous methods, which in turn improves the utility of MIPs for cost-effective targeted sequencing for candidate gene validation and for diagnostic sequencing in a clinical setting. AVAILABILITY AND IMPLEMENTATION MIPgen is implemented in C++. Source code and accompanying Python scripts are available at http://shendurelab.github.io/MIPGEN/.

Enhancer connectome in primary human cells identifies target genes of disease-associated DNA elements

The challenge of linking intergenic mutations to target genes has limited molecular understanding of human diseases. Here we show that H3K27ac HiChIP generates high-resolution contact maps of active enhancers and target genes in rare primary human T cell subtypes and coronary artery smooth muscle cells. Differentiation of naive T cells into T helper 17 cells or regulatory T cells creates subtype-specific enhancer–promoter interactions, specifically at regions of shared DNA accessibility. These data provide a principled means of assigning molecular functions to autoimmune and cardiovascular disease risk variants, linking hundreds of noncoding variants to putative gene targets. Target genes identified with HiChIP are further supported by CRISPR interference and activation at linked enhancers, by the presence of expression quantitative trait loci, and by allele-specific enhancer loops in patient-derived primary cells. The majority of disease-associated enhancers contact genes beyond the nearest gene in the linear genome, leading to a fourfold increase in the number of potential target genes for autoimmune and cardiovascular diseases.