High-resolution mapping of cancer cell networks using co-functional interactions.

Powerful new technologies for perturbing genetic elements have recently expanded the study of genetic interactions in model systems ranging from yeast to human cell lines. However, technical artifacts can confound signal across genetic screens and limit the immense potential of parallel screening approaches. To address this problem, we devised a novel PCA-based method for correcting genome-wide screening data, bolstering the sensitivity and specificity of detection for genetic interactions. Applying this strategy to a set of 436 whole genome CRISPR screens, we report more than 1.5 million pairs of correlated "co-functional" genes that provide finer-scale information about cell compartments, biological pathways, and protein complexes than traditional gene sets. Lastly, we employed a gene community detection approach to implicate core genes for cancer growth and compress signal from functionally related genes in the same community into a single score. This work establishes new algorithms for probing cancer cell networks and motivates the acquisition of further CRISPR screen data across diverse genotypes and cell types to further resolve complex cellular processes.

Deactivated CRISPR Associated Protein 9 for Minor-Allele Enrichment in Cell-Free DNA.

BACKGROUND: Cell-free DNA (cfDNA) diagnostics are emerging as a new paradigm of disease monitoring and therapy management. The clinical utility of these diagnostics is relatively limited by a low signal-to-noise ratio, such as with low allele frequency (AF) mutations in cancer. While enriching for rare alleles to increase their AF before sample analysis is one strategy that can greatly improve detection capability, current methods are limited in their generalizability, ease of use, and applicability to point mutations. METHODS: Leveraging the robust single-base-pair specificity and generalizability of the CRISPR associated protein 9 (Cas9) system, we developed a deactivated Cas9 (dCas9)-based method of minor-allele enrichment capable of efficient single-target and multiplexed enrichment. The dCas9 protein was complexed with single guide RNAs targeted to mutations of interest and incubated with cfDNA samples containing mutant strands at low abundance. Mutation-bound dCas9 complexes were isolated, dissociated, and the captured DNA purified for downstream use. RESULTS: Targeting the 3 most common epidermal growth factor receptor mutations (exon 19 deletion, T790M, L858R) found in non-small cell lung cancer (NSCLC), we achieved >20-fold increases in AF and detected mutations by use of qPCR at an AF of 0.1%. In a cohort of 18 NSCLC patient-derived cfDNA samples, our method enabled detection of 8 out of 13 mutations that were otherwise undetected by qPCR. CONCLUSIONS: The dCas9 method provides an important application of the CRISPR/Cas9 system outside the realm of genome editing and can provide a step forward for the detection capability of cfDNA diagnostics.

Caspase-8, association with Alzheimer's Disease and functional analysis of rare variants.

The accumulation of amyloid beta (Aß) peptide (Amyloid cascade hypothesis), an APP protein cleavage product, is a leading hypothesis in the etiology of Alzheimer's disease (AD). In order to identify additional AD risk genes, we performed targeted sequencing and rare variant burden association study for nine candidate genes involved in the amyloid metabolism in 1886 AD cases and 1700 controls. We identified a significant variant burden association for the gene encoding caspase-8, CASP8 (p = 8.6x10-5). For two CASP8 variants, p.K148R and p.I298V, the association remained significant in a combined sample of 10,820 cases and 8,881 controls. For both variants we performed bioinformatics structural, expression and enzymatic activity studies and obtained evidence for loss of function effects. In addition to their role in amyloid processing, caspase-8 and its downstream effector caspase-3 are involved in synaptic plasticity, learning, memory and control of microglia pro-inflammatory activation and associated neurotoxicity, indicating additional mechanisms that might contribute to AD. As caspase inhibition has been proposed as a mechanism for AD treatment, our finding that AD-associated CASP8 variants reduce caspase function calls for caution and is an impetus for further studies on the role of caspases in AD and other neurodegenerative diseases.

BRCA Testing by Single-Molecule Molecular Inversion Probes.

BACKGROUND: Despite advances in next generation DNA sequencing (NGS), NGS-based single gene tests for diagnostic purposes require improvements in terms of completeness, quality, speed, and cost. Single-molecule molecular inversion probes (smMIPs) are a technology with unrealized potential in the area of clinical genetic testing. In this proof-of-concept study, we selected 2 frequently requested gene tests, those for the breast cancer genes BRCA1 and BRCA2, and developed an automated work flow based on smMIPs. METHODS: The BRCA1 and BRCA2 smMIPs were validated using 166 human genomic DNA samples with known variant status. A generic automated work flow was built to perform smMIP-based enrichment and sequencing for BRCA1, BRCA2, and the checkpoint kinase 2 (CHEK2) c.1100del variant. RESULTS: Pathogenic and benign variants were analyzed in a subset of 152 previously BRCA-genotyped samples, yielding an analytical sensitivity and specificity of 100%. Following automation, blind analysis of 65 in-house samples and 267 Norwegian samples correctly identified all true-positive variants (>3000), with no false positives. Consequent to process optimization, turnaround times were reduced by 60% to currently 10-15 days. Copy number variants were detected with an analytical sensitivity of 100% and an analytical specificity of 88%. CONCLUSIONS: smMIP-based genetic testing enables automated and reliable analysis of the coding sequences of BRCA1 and BRCA2. The use of single-molecule tags, double-tiled targeted enrichment, and capturing and sequencing in duplo, in combination with automated library preparation and data analysis, results in a robust process and reduces routine turnaround times. Furthermore, smMIP-based copy number variation analysis could make independent copy number variation tools like multiplex ligation-dependent probes amplification dispensable.

PIK3CA-associated developmental disorders exhibit distinct classes of mutations with variable expression and tissue distribution.

Mosaicism is increasingly recognized as a cause of developmental disorders with the advent of next-generation sequencing (NGS). Mosaic mutations of PIK3CA have been associated with the widest spectrum of phenotypes associated with overgrowth and vascular malformations. We performed targeted NGS using 2 independent deep-coverage methods that utilize molecular inversion probes and amplicon sequencing in a cohort of 241 samples from 181 individuals with brain and/or body overgrowth. We identified PIK3CA mutations in 60 individuals. Several other individuals (n = 12) were identified separately to have mutations in PIK3CA by clinical targeted-panel testing (n = 6), whole-exome sequencing (n = 5), or Sanger sequencing (n = 1). Based on the clinical and molecular features, this cohort segregated into three distinct groups: (a) severe focal overgrowth due to low-level but highly activating (hotspot) mutations, (b) predominantly brain overgrowth and less severe somatic overgrowth due to less-activating mutations, and (c) intermediate phenotypes (capillary malformations with overgrowth) with intermediately activating mutations. Sixteen of 29 PIK3CA mutations were novel. We also identified constitutional PIK3CA mutations in 10 patients. Our molecular data, combined with review of the literature, show that PIK3CA-related overgrowth disorders comprise a discontinuous spectrum of disorders that correlate with the severity and distribution of mutations.

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 Children's 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.

Characterisation of mutations of the phosphoinositide-3-kinase regulatory subunit, PIK3R2, in perisylvian polymicrogyria: a next-generation sequencing study.

BACKGROUND: Bilateral perisylvian polymicrogyria (BPP), the most common form of regional polymicrogyria, causes the congenital bilateral perisylvian syndrome, featuring oromotor dysfunction, cognitive impairment, and epilepsy. The causes of BPP are heterogeneous, but only a few genetic causes have been reported. The aim of this study was to identify additional genetic causes of BPP and characterise their frequency in this population. METHODS: Children (aged ≤18 years) with polymicrogyria were enrolled into our research programme from July, 1980, to October, 2015, at two centres (Florence, Italy, and Seattle, WA, USA). We obtained samples (blood and saliva) throughout this period at both centres and did whole-exome sequencing on DNA from eight trios (two parents and one affected child) with BPP in 2014. After the identification of mosaic PIK3R2 mutations in two of these eight children, we performed targeted screening of PIK3R2 by two methods in a cohort of 118 children with BPP. First, we performed targeted sequencing of the entire PIK3R2 gene by single molecule molecular inversion probes (smMIPs) on 38 patients with BPP with normal to large head size. Second, we did amplicon sequencing of the recurrent PIK3R2 mutation (Gly373Arg) in 80 children with various types of polymicrogyria including BPP. One additional patient had clinical whole-exome sequencing done independently, and was included in this study because of the phenotypic similarity to our cohort. FINDINGS: We identified a mosaic mutation (Gly373Arg) in a regulatory subunit of the PI3K-AKT-mTOR pathway, PIK3R2, in two children with BPP. Of the 38 patients with BPP and normal to large head size who underwent targeted next-generation sequencing by smMIPs, we identified constitutional and mosaic PIK3R2 mutations in 17 additional children. In parallel, one patient had the recurrent PIK3R2 mutation identified by clinical whole-exome sequencing. Seven of these 20 patients had BPP alone, and 13 had BPP in association with features of the megalencephaly-polymicrogyria-polydactyly-hydrocephalus (MPPH) syndrome. 19 patients had the same mutation (Gly373Arg), and one had a nearby missense mutation (Lys376Glu). Mutations were constitutional in 12 patients and mosaic in eight patients. In patients with mosaic mutations, we noted substantial variation in alternate (mutant) allele levels, ranging from ten (3%) of 377 reads to 39 (37%) of 106 reads, equivalent to 5-73% of cells analysed. Levels of mosaicism varied from undetectable to 37 (17%) of 216 reads in blood-derived DNA compared with 2030 (29%) of 6889 reads to 275 (43%) of 634 reads in saliva-derived DNA. INTERPRETATION: Constitutional and mosaic mutations in the PIK3R2 gene are associated with developmental brain disorders ranging from BPP with a normal head size to the MPPH syndrome. The phenotypic variability and low-level mosaicism, which challenge conventional molecular methods, have important implications for genetic testing and counselling. FUNDING: US National Institutes of Health.

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.

KIAA0586 is Mutated in Joubert Syndrome.

Joubert syndrome (JS) is a recessive neurodevelopmental disorder characterized by a distinctive mid-hindbrain malformation. JS is part of a group of disorders called ciliopathies based on their overlapping phenotypes and common underlying pathophysiology linked to primary cilium dysfunction. Biallelic mutations in one of 28 genes, all encoding proteins localizing to the primary cilium or basal body, can cause JS. Despite this large number of genes, the genetic cause can currently be determined in about 62% of individuals with JS. To identify novel JS genes, we performed whole exome sequencing on 35 individuals with JS and found biallelic rare deleterious variants (RDVs) in KIAA0586, encoding a centrosomal protein required for ciliogenesis, in one individual. Targeted next-generation sequencing in a large JS cohort identified biallelic RDVs in eight additional families for an estimated prevalence of 2.5% (9/366 JS families). All affected individuals displayed JS phenotypes toward the mild end of the spectrum.

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.

De novo mutations in NALCN cause a syndrome characterized by congenital contractures of the limbs and face, hypotonia, and developmental delay.

Freeman-Sheldon syndrome, or distal arthrogryposis type 2A (DA2A), is an autosomal-dominant condition caused by mutations in MYH3 and characterized by multiple congenital contractures of the face and limbs and normal cognitive development. We identified a subset of five individuals who had been putatively diagnosed with "DA2A with severe neurological abnormalities" and for whom congenital contractures of the limbs and face, hypotonia, and global developmental delay had resulted in early death in three cases; this is a unique condition that we now refer to as CLIFAHDD syndrome. Exome sequencing identified missense mutations in the sodium leak channel, non-selective (NALCN) in four families affected by CLIFAHDD syndrome. We used molecular-inversion probes to screen for NALCN in a cohort of 202 distal arthrogryposis (DA)-affected individuals as well as concurrent exome sequencing of six other DA-affected individuals, thus revealing NALCN mutations in ten additional families with "atypical" forms of DA. All 14 mutations were missense variants predicted to alter amino acid residues in or near the S5 and S6 pore-forming segments of NALCN, highlighting the functional importance of these segments. In vitro functional studies demonstrated that NALCN alterations nearly abolished the expression of wild-type NALCN, suggesting that alterations that cause CLIFAHDD syndrome have a dominant-negative effect. In contrast, homozygosity for mutations in other regions of NALCN has been reported in three families affected by an autosomal-recessive condition characterized mainly by hypotonia and severe intellectual disability. Accordingly, mutations in NALCN can cause either a recessive or dominant condition characterized by varied though overlapping phenotypic features, perhaps based on the type of mutation and affected protein domain(s).

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.

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.

Deep sequencing of multiple regions of glial tumors reveals spatial heterogeneity for mutations in clinically relevant genes.

BACKGROUND: The extent of intratumoral mutational heterogeneity remains unclear in gliomas, the most common primary brain tumors, especially with respect to point mutation. To address this, we applied single molecule molecular inversion probes targeting 33 cancer genes to assay both point mutations and gene amplifications within spatially distinct regions of 14 glial tumors. RESULTS: We find evidence of regional mutational heterogeneity in multiple tumors, including mutations in TP53 and RB1 in an anaplastic oligodendroglioma and amplifications in PDGFRA and KIT in two glioblastomas (GBMs). Immunohistochemistry confirms heterogeneity of TP53 mutation and PDGFRA amplification. In all, 3 out of 14 glial tumors surveyed have evidence for heterogeneity for clinically relevant mutations. CONCLUSIONS: Our results underscore the need to sample multiple regions in GBM and other glial tumors when devising personalized treatments based on genomic information, and furthermore demonstrate the importance of measuring both point mutation and copy number alteration while investigating genetic heterogeneity within cancer samples.

Germline missense variants in the BTNL2 gene are associated with prostate cancer susceptibility.

BACKGROUND: Rare, inherited mutations account for 5% to 10% of all prostate cancer cases. However, to date, few causative mutations have been identified. METHODS: To identify rare mutations for prostate cancer, we conducted whole-exome sequencing (WES) in multiple kindreds (n = 91) from 19 hereditary prostate cancer (HPC) families characterized by aggressive or early-onset phenotypes. Candidate variants (n = 130) identified through family- and bioinformatics-based filtering of WES data were then genotyped in an independent set of 270 HPC families (n = 819 prostate cancer cases; n = 496 unaffected relatives) for replication. Two variants with supportive evidence were subsequently genotyped in a population-based case-control study (n = 1,155 incident prostate cancer cases; n = 1,060 age-matched controls) for further confirmation. All participants were men of European ancestry. RESULTS: The strongest evidence was for two germline missense variants in the butyrophilin-like 2 (BTNL2) gene (rs41441651, p.Asp336Asn and rs28362675, p.Gly454Cys) that segregated with affection status in two of the WES families. In the independent set of 270 HPC families, 1.5% (rs41441651; P = 0.0032) and 1.2% (rs28362675; P = 0.0070) of affected men, but no unaffected men, carried a variant. Both variants were associated with elevated prostate cancer risk in the population-based study (rs41441651: OR, 2.7; 95% CI, 1.27-5.87; P = 0.010; rs28362675: OR, 2.5; 95% CI, 1.16-5.46; P = 0.019). CONCLUSIONS: Results indicate that rare BTNL2 variants play a role in susceptibility to both familial and sporadic prostate cancer. IMPACT: Results implicate BTNL2 as a novel prostate cancer susceptibility gene.