A benchmark of computational methods for correcting biases of established and unknown origin in CRISPR-Cas9 screening data.

BACKGROUND: CRISPR-Cas9 dropout screens are formidable tools for investigating biology with unprecedented precision and scale. However, biases in data lead to potential confounding effects on interpretation and compromise overall quality. The activity of Cas9 is influenced by structural features of the target site, including copy number amplifications (CN bias). More worryingly, proximal targeted loci tend to generate similar gene-independent responses to CRISPR-Cas9 targeting (proximity bias), possibly due to Cas9-induced whole chromosome-arm truncations or other genomic structural features and different chromatin accessibility levels. RESULTS: We benchmarked eight computational methods, rigorously evaluating their ability to reduce both CN and proximity bias in the two largest publicly available cell-line-based CRISPR-Cas9 screens to date. We also evaluated the capability of each method to preserve data quality and heterogeneity by assessing the extent to which the processed data allows accurate detection of true positive essential genes, established oncogenetic addictions, and known/novel biomarkers of cancer dependency. Our analysis sheds light on the ability of each method to correct biases under different scenarios. AC-Chronos outperforms other methods in correcting both CN and proximity biases when jointly processing multiple screens of models with available CN information, whereas CRISPRcleanR is the top performing method for individual screens or when CN information is not available. In addition, Chronos and AC-Chronos yield a final dataset better able to recapitulate known sets of essential and non-essential genes. CONCLUSIONS: Overall, our investigation provides guidance for the selection of the most appropriate bias-correction method, based on its strengths, weaknesses and experimental settings.

Next-generation sequencing of baseline genetic mutations and outcomes of eltrombopag and azacitidine therapy in patients with myelodysplastic syndromes and thrombocytopenia: Data from the SUPPORT clinical trial.

Eltrombopag has been previously shown to be effective in reversing azacitidine-mediated thrombocytopenia. This was further investigated in the SUPPORT trial, a phase III study assessing the efficacy/safety of eltrombopag plus azacitidine in patients with intermediate- to high-risk myelodysplastic syndromes and thrombocytopenia. The results did not support a clinical benefit for the addition of eltrombopag to azacitidine. We investigated if the somatic mutational profiles in the patient cohort were associated with treatment outcomes. Based on the available data, we observed no imbalance in the mutational profiles between treatment arms or a clear association between identified somatic mutations and clinical outcomes.

Initial Evidence for the Efficacy of Naporafenib in Combination With Trametinib in NRAS-Mutant Melanoma: Results From the Expansion Arm of a Phase Ib, Open-Label Study.

PURPOSE: No approved targeted therapy for the treatment of patients with neuroblastoma RAS viral (v-ras) oncogene homolog (NRAS)-mutant melanoma is currently available. PATIENTS AND METHODS: In this phase Ib escalation/expansion study (ClinicalTrials.gov identifier: NCT02974725), the safety, tolerability, and preliminary antitumor activity of naporafenib (LXH254), a BRAF/CRAF protein kinases inhibitor, were explored in combination with trametinib in patients with advanced/metastatic KRAS- or BRAF-mutant non-small-cell lung cancer (escalation arm) or NRAS-mutant melanoma (escalation and expansion arms). RESULTS: Thirty-six and 30 patients were enrolled in escalation and expansion, respectively. During escalation, six patients reported grade ≥3 dose-limiting toxicities, including dermatitis acneiform (n = 2), maculopapular rash (n = 2), increased lipase (n = 1), and Stevens-Johnson syndrome (n = 1). The recommended doses for expansion were naporafenib 200 mg twice a day plus trametinib 1 mg once daily and naporafenib 400 mg twice a day plus trametinib 0.5 mg once daily. During expansion, all 30 patients experienced a treatment-related adverse event, the most common being rash (80%, n = 24), blood creatine phosphokinase increased, diarrhea, and nausea (30%, n = 9 each). In expansion, the objective response rate, median duration of response, and median progression-free survival were 46.7% (95% CI, 21.3 to 73.4; 7 of 15 patients), 3.75 (95% CI, 1.97 to not estimable [NE]) months, and 5.52 months, respectively, in patients treated with naporafenib 200 mg twice a day plus trametinib 1 mg once daily, and 13.3% (95% CI, 1.7 to 40.5; 2 of 15 patients), 3.75 (95% CI, 2.04 to NE) months, and 4.21 months, respectively, in patients treated with naporafenib 400 mg twice a day plus trametinib 0.5 mg once daily. CONCLUSION: Naporafenib plus trametinib showed promising preliminary antitumor activity in patients with NRAS-mutant melanoma. Prophylactic strategies aimed to lower the incidence of skin-related events are under investigation.

Clonal hematopoiesis detection in patients with cancer using cell-free DNA sequencing.

In the context of cancer, clonal hematopoiesis of indeterminate potential (CHIP) is associated with the development of therapy-related myeloid neoplasms and shorter overall survival. Cell-free DNA (cfDNA) sequencing is becoming widely adopted for genomic screening of patients with cancer but has not been used extensively to determine CHIP status because of a requirement for matched blood and tumor sequencing. We present an accurate classification approach to determine the CH status from cfDNA sequencing alone, applying our model to 4324 oncology clinical cfDNA samples. Using this method, we determined that 30.3% of patients in this cohort have evidence of CH, and the incidence of CH varies by tumor type. Matched RNA sequencing data show evidence of increased inflammation, especially neutrophil activation, within the tumors and tumor microenvironments of patients with CH. In addition, patients with CH had evidence of neutrophil activation systemically, pointing to a potential mechanism of action for the worse outcomes associated with CH status. Neutrophil activation may be one of many mechanisms, however, because patients with estrogen receptor-positive breast cancer harboring TET2 frameshift mutations had worse outcomes but similar neutrophil frequencies to patients without CH. Together, these data show the feasibility of detecting CH through cfDNA sequencing alone and an application of this method, demonstrating increased inflammation in patients with CH both systemically and in the tumor microenvironment.

Molecular response in newly diagnosed chronic-phase chronic myeloid leukemia: prediction modeling and pathway analysis.

Tyrosine kinase inhibitor therapy revolutionized chronic myeloid leukemia treatment and showed how targeted therapy and molecular monitoring could be used to substantially improve survival outcomes. We used chronic myeloid leukemia as a model to understand a critical question: why do some patients have an excellent response to therapy, while others have a poor response? We studied gene expression in whole blood samples from 112 patients from a large phase III randomized trial (clinicaltrials gov. Identifier: NCT00471497), dichotomizing cases into good responders (BCR::ABL1 ≤10% on the International Scale by 3 and 6 months and ≤0.1% by 12 months) and poor responders (failure to meet these criteria). Predictive models based on gene expression demonstrated the best performance (area under the curve =0.76, standard deviation =0.07). All of the top 20 pathways overexpressed in good responders involved immune regulation, a finding validated in an independent data set. This study emphasizes the importance of pretreatment adaptive immune response in treatment efficacy and suggests biological pathways that can be targeted to improve response.

Combined PD-1, BRAF and MEK inhibition in BRAFV600E colorectal cancer: a phase 2 trial.

While BRAF inhibitor combinations with EGFR and/or MEK inhibitors have improved clinical efficacy in BRAFV600E colorectal cancer (CRC), response rates remain low and lack durability. Preclinical data suggest that BRAF/MAPK pathway inhibition may augment the tumor immune response. We performed a proof-of-concept single-arm phase 2 clinical trial of combined PD-1, BRAF and MEK inhibition with sparatlizumab (PDR001), dabrafenib and trametinib in 37 patients with BRAFV600E CRC. The primary end point was overall response rate, and the secondary end points were progression-free survival, disease control rate, duration of response and overall survival. The study met its primary end point with a confirmed response rate (24.3% in all patients; 25% in microsatellite stable patients) and durability that were favorable relative to historical controls of BRAF-targeted combinations alone. Single-cell RNA sequencing of 23 paired pretreatment and day 15 on-treatment tumor biopsies revealed greater induction of tumor cell-intrinsic immune programs and more complete MAPK inhibition in patients with better clinical outcome. Immune program induction in matched patient-derived organoids correlated with the degree of MAPK inhibition. These data suggest a potential tumor cell-intrinsic mechanism of cooperativity between MAPK inhibition and immune response, warranting further clinical evaluation of optimized targeted and immune combinations in CRC. ClinicalTrials.gov registration: NCT03668431.

Contribution and clinical relevance of germline variation to the cancer transcriptome.

BACKGROUND: Somatic alterations in the cancer genome, some of which are associated with changes in gene expression, have been characterized in multiple studies across diverse cancer types. However, less is known about germline variants that influence tumor biology by shaping the cancer transcriptome. METHODS: We performed expression quantitative trait loci (eQTL) analyses using multi-dimensional data from The Cancer Genome Atlas to explore the role of germline variation in mediating the cancer transcriptome. After accounting for associations between somatic alterations and gene expression, we determined the contribution of inherited variants to the cancer transcriptome relative to that of somatic variants. Finally, we performed an interaction analysis using estimates of tumor cellularity to identify cell type-restricted eQTLs. RESULTS: The proportion of genes with at least one eQTL varied between cancer types, ranging between 0.8% in melanoma to 28.5% in thyroid cancer and was correlated more strongly with intratumor heterogeneity than with somatic alteration rates. Although contributions to variance in gene expression was low for most genes, some eQTLs accounted for more than 30% of expression of proximal genes. We identified cell type-restricted eQTLs in genes known to be cancer drivers including LPP and EZH2 that were associated with disease-specific mortality in TCGA but not associated with disease risk in published GWAS. Together, our results highlight the need to consider germline variation in interpreting cancer biology beyond risk prediction.

TET2-Driven Clonal Hematopoiesis and Response to Canakinumab: An Exploratory Analysis of the CANTOS Randomized Clinical Trial.

Importance: Clonal hematopoiesis of indeterminate potential (CHIP) is associated with increased risk of atherosclerotic cardiovascular disease, and mouse experiments suggest that CHIP related to Tet2 loss of function in myeloid cells accelerates atherosclerosis via augmented interleukin (IL) 1ß signaling. Objective: To assess whether individuals with CHIP have greater cardiovascular event reduction in response to IL-1ß neutralization in the Canankinumab Anti-inflammatory Thrombosis Outcomes Trial (CANTOS). Design, Setting, and Participants: This randomized clinical trial took place from April 2011 to June 2017 at more than 1000 clinical sites in 39 countries. Targeted deep sequencing of genes previously associated with CHIP in a subset of trial participants using genomic DNA prepared from baseline peripheral blood samples were analyzed. All participants had prior myocardial infarction and elevated high-sensitivity C-reactive protein level above 0.20 mg/dL. Analysis took place between June 2017 and December 2021. Interventions: Canakinumab, an anti-IL-1ß antibody, given at doses of 50, 150, and 300 mg once every 3 months. Main Outcomes and Measures: Major adverse cardiovascular events (MACE). Results: A total of 338 patients (8.6%) were identified in this subset with evidence for clonal hematopoiesis. As expected, the incidence of CHIP increased with age; the mean (SD) age of patients with CHIP was 66.3 (9.2) years and 61.5 (9.6) years in patients without CHIP. Unlike other populations that were not preselected for elevated C-reactive protein, in the CANTOS population variants in TET2 were more common than DNMT3A (119 variants in 103 patients vs 86 variants in 85 patients). Placebo-treated patients with CHIP showed a nonsignificant increase in the rate of MACE compared with patients without CHIP using a Cox proportional hazard model (hazard ratio, 1.32 [95% CI, 0.86-2.04]; P = .21). Exploratory analyses of placebo-treated patients with a somatic variant in either TET2 or DNMT3A (n = 58) showed an equivocal risk for MACE (hazard ratio, 1.65 [95% CI, 0.97-2.80]; P = .06). Patients with CHIP due to somatic variants in TET2 also had reduced risk for MACE while taking canakinumab (hazard ratio, 0.38 [95% CI, 0.15-0.96]) with equivocal difference compared with others (P for interaction = .14). Conclusions and Relevance: These results are consistent with observations of increased risk for cardiovascular events in patients with CHIP and raise the possibility that those with TET2 variants may respond better to canakinumab than those without CHIP. Future studies are required to further substantiate this hypothesis. Trial Registration: ClinicalTrials.gov Identifier: NCT01327846.

Cell-free DNA captures tumor heterogeneity and driver alterations in rapid autopsies with pre-treated metastatic cancer.

In patients with metastatic cancer, spatial heterogeneity of somatic alterations may lead to incomplete assessment of a cancer's mutational profile when analyzing a single tumor biopsy. In this study, we perform sequencing of cell-free DNA (cfDNA) and distinct metastatic tissue samples from ten rapid autopsy cases with pre-treated metastatic cancer. We show that levels of heterogeneity in genetic biomarkers vary between patients but that gene expression signatures representative of the tumor microenvironment are more consistent. Across nine patients with plasma samples available, we are able to detect 62/62 truncal and 47/121 non-truncal point mutations in cfDNA. We observe that mutation clonality in cfDNA is correlated with the number of metastatic lesions in which the mutation is detected and use this result to derive a clonality threshold to classify truncal and non-truncal driver alterations with reasonable specificity. In contrast, mutation truncality is more often incorrectly assigned when studying single tissue samples. Our results demonstrate the utility of a single cfDNA sample relative to that of single tissue samples when treating patients with metastatic cancer.

Combined PD-1, BRAF and MEK inhibition in advanced BRAF-mutant melanoma: safety run-in and biomarker cohorts of COMBI-i.

Immune and targeted therapies achieve long-term survival in metastatic melanoma; however, new treatment strategies are needed to improve patients' outcomes1,2. We report on the efficacy, safety and biomarker analysis from the single-arm safety run-in (part 1; n = 9) and biomarker (part 2; n = 27) cohorts of the randomized, placebo-controlled, phase 3 COMBI-i trial (NCT02967692) of the anti-PD-1 antibody spartalizumab, in combination with the BRAF inhibitor dabrafenib and MEK inhibitor trametinib. Patients (n = 36) had previously untreated BRAF V600-mutant unresectable or metastatic melanoma. In part 1, the recommended phase 3 regimen was identified based on the incidence of dose-limiting toxicities (DLTs; primary endpoint): 400 mg of spartalizumab every 4 weeks plus 150 mg of dabrafenib twice daily plus 2 mg of trametinib once daily. Part 2 characterized changes in PD-L1 levels and CD8+ cells following treatment (primary endpoint), and analyzed additional biomarkers. Assessments of efficacy and safety were key secondary endpoints (median follow-up, 24.3 months). Spartalizumab plus dabrafenib and trametinib led to an objective response rate (ORR) of 78%, including 44% complete responses (CRs). Grade ≥3 treatment-related adverse events (TRAEs) were experienced by 72% of patients. All patients had temporary dose modifications, and 17% permanently discontinued all three study drugs due to TRAEs. Early progression-free survival (PFS) events were associated with low tumor mutational burden/T cell-inflamed gene expression signature (GES) or high immunosuppressive tumor microenvironment (TME) GES levels at baseline; an immunosuppressive TME may also preclude CR. Overall, the efficacy, safety and on-treatment biomarker modulations associated with spartalizumab plus dabrafenib and trametinib are promising, and biomarkers that may predict long-term benefit were identified.

BRAF-Mutant Transcriptional Subtypes Predict Outcome of Combined BRAF, MEK, and EGFR Blockade with Dabrafenib, Trametinib, and Panitumumab in Patients with Colorectal Cancer.

PURPOSE: The influence of the transcriptional and immunologic context of mutations on therapeutic outcomes with targeted therapy in cancer has not been well defined. BRAF V600E-mutant (BM) colorectal cancer comprises two main transcriptional subtypes, BM1 and BM2. We sought to determine the impact of BM subtype, as well as distinct biological features of those subtypes, on response to BRAF/MEK/EGFR inhibition in patients with colorectal cancer. PATIENTS AND METHODS: Paired fresh tumor biopsies were acquired at baseline and on day 15 of treatment from all consenting patients with BM colorectal cancer enrolled in a phase II clinical trial of dabrafenib, trametinib, and panitumumab. For each sample, BM subtype, cell cycle, and immune gene signature expression were determined using RNA-sequencing (RNA-seq), and a Cox proportional hazards model was applied to determine association with progression-free survival (PFS). RESULTS: Confirmed response rates, median PFS, and median overall survival (OS) were higher in BM1 subtype patients compared with BM2 subtype patients. Evaluation of immune contexture identified greater immune reactivity in BM1, whereas cell-cycle signatures were more highly expressed in BM2. A multivariate model of PFS incorporating BM subtype plus immune and cell-cycle signatures revealed that BM subtype encompasses the majority of the effect. CONCLUSIONS: BM subtype is significantly associated with the outcome of combination dabrafenib, trametinib, and panitumumab therapy and may serve as a standalone predictive biomarker beyond mutational status. Our findings support a more nuanced approach to targeted therapeutic decisions that incorporates assessment of transcriptional context.

Adjuvant dabrafenib plus trametinib versus placebo in patients with resected, BRAFV600-mutant, stage III melanoma (COMBI-AD): exploratory biomarker analyses from a randomised, phase 3 trial.

BACKGROUND: Adjuvant dabrafenib plus trametinib reduced the risk of relapse versus placebo in patients with resected, BRAFV600-mutant, stage III melanoma in the phase 3 COMBI-AD trial. This prespecified exploratory biomarker analysis aimed to evaluate potential prognostic or predictive factors and mechanisms of resistance to adjuvant targeted therapy. METHODS: COMBI-AD is a randomised, double-blind, placebo-controlled, phase 3 trial comparing dabrafenib 150 mg orally twice daily plus trametinib 2 mg orally once daily versus two matched placebos. Study participants were at least 18 years of age and underwent complete resection of stage IIIA (lymph node metastases >1 mm), IIIB, or IIIC cutaneous melanoma, per American Joint Committee on Cancer 7th edition criteria, with a BRAFV600E or BRAFV600K mutation. Patients were randomly assigned (1:1) to the two treatment groups by an interactive voice response system, stratified by mutation type and disease stage. Patients, physicians, and the investigators who analysed data were masked to treatment allocation. The primary outcome was relapse-free survival, defined as the time from randomisation to disease recurrence or death from any cause. Biomarker assessment was a prespecified exploratory outcome of the trial. We assessed intrinsic tumour genomic features by use of next-generation DNA sequencing and characteristics of the tumour microenvironment by use of a NanoString RNA assay, which might provide prognostic and predictive information. This trial is registered with ClinicalTrials.gov, number NCT01682083, and is ongoing but no longer recruiting participants. FINDINGS: Between Jan 31, 2013, and Dec 11, 2014, 870 patients were enrolled in the trial. Median follow-up at data cutoff (April 30, 2018) was 44 months (IQR 38-49) in the dabrafenib plus trametinib group and 42 months (21-49) in the placebo group. Intrinsic tumour genomic features were assessed in 368 patients (DNA sequencing set) and tumour microenvironment characteristics were assessed in 507 patients (NanoString biomarker set). MAPK pathway genomic alterations at baseline did not affect treatment benefit or clinical outcome. An IFNγ gene expression signature higher than the median was prognostic for prolonged relapse-free survival in both treatment groups. Tumour mutational burden was independently prognostic for relapse-free survival in the placebo group (high TMB, top third; hazard ratio [HR] 0·56, 95% CI 0·37-0·85, p=0·0056), but not in the dabrafenib plus trametinib group (0·83, 95% CI 0·53-1·32, p=0·44). Patients with tumour mutational burden in the lower two terciles seem to derive a substantial long-term relapse-free survival benefit from targeted therapy (HR [versus placebo] 0·49, 95% CI 0·35-0·68, p<0·0001). However, patients with high tumour mutational burden seem to have a less pronounced benefit with targeted therapy (HR [versus placebo] 0·75, 95% CI 0·44-1·26, p=0·27), especially if they had an IFNγ signature lower than the median (HR 0·88 [95% CI 0·40-1·93], p=0·74). INTERPRETATION: Tumour mutational burden alone or in combination with IFNγ gene expression signature or other markers for an adaptive immune response might be of relevance for identifying patients with stage III melanoma who might derive clinical benefit from targeted therapy. Further validation in prospective clinical trials is warranted. FUNDING: Novartis Pharmaceuticals.

Tracking the Evolution of Resistance to ALK Tyrosine Kinase Inhibitors through Longitudinal Analysis of Circulating Tumor DNA.

PURPOSE: ALK rearrangements predict for sensitivity to ALK tyrosine kinase inhibitors (TKIs). However, responses to ALK TKIs are generally short-lived. Serial molecular analysis is an informative strategy for identifying genetic mediators of resistance. Although multiple studies support the clinical benefits of repeat tissue sampling, the clinical utility of longitudinal circulating tumor DNA analysis has not been established in ALK-positive lung cancer. METHODS: Using a 566-gene hybrid-capture next-generation sequencing (NGS) assay, we performed longitudinal analysis of plasma specimens from 22 ALK-positive patients with acquired resistance to ALK TKIs to track the evolution of resistance during treatment. To determine tissue-plasma concordance, we compared plasma findings to results of repeat biopsies. RESULTS: At progression, we detected an ALK fusion in plasma from 19 (86%) of 22 patients, and identified ALK resistance mutations in plasma specimens from 11 (50%) patients. There was 100% agreement between tissue- and plasma-detected ALK fusions. Among 16 cases where contemporaneous plasma and tissue specimens were available, we observed 100% concordance between ALK mutation calls. ALK mutations emerged and disappeared during treatment with sequential ALK TKIs, suggesting that plasma mutation profiles were dependent on the specific TKI administered. ALK G1202R, the most frequent plasma mutation detected after progression on a second-generation TKI, was consistently suppressed during treatment with lorlatinib. CONCLUSIONS: Plasma genotyping by NGS is an effective method for detecting ALK fusions and ALK mutations in patients progressing on ALK TKIs. The correlation between plasma ALK mutations and response to distinct ALK TKIs highlights the potential for plasma analysis to guide selection of ALK-directed therapies.

PureCN: copy number calling and SNV classification using targeted short read sequencing.

BACKGROUND: Matched sequencing of both tumor and normal tissue is routinely used to classify variants of uncertain significance (VUS) into somatic vs. germline. However, assays used in molecular diagnostics focus on known somatic alterations in cancer genes and often only sequence tumors. Therefore, an algorithm that reliably classifies variants would be helpful for retrospective exploratory analyses. Contamination of tumor samples with normal cells results in differences in expected allelic fractions of germline and somatic variants, which can be exploited to accurately infer genotypes after adjusting for local copy number. However, existing algorithms for determining tumor purity, ploidy and copy number are not designed for unmatched short read sequencing data. RESULTS: We describe a methodology and corresponding open source software for estimating tumor purity, copy number, loss of heterozygosity (LOH), and contamination, and for classification of single nucleotide variants (SNVs) by somatic status and clonality. This R package, PureCN, is optimized for targeted short read sequencing data, integrates well with standard somatic variant detection pipelines, and has support for matched and unmatched tumor samples. Accuracy is demonstrated on simulated data and on real whole exome sequencing data. CONCLUSIONS: Our algorithm provides accurate estimates of tumor purity and ploidy, even if matched normal samples are not available. This in turn allows accurate classification of SNVs. The software is provided as open source (Artistic License 2.0) R/Bioconductor package PureCN (http://bioconductor.org/packages/PureCN/).

Interchromosomal core duplicons drive both evolutionary instability and disease susceptibility of the Chromosome 8p23.1 region.

Recurrent rearrangements of Chromosome 8p23.1 are associated with congenital heart defects and developmental delay. The complexity of this region has led to inconsistencies in the current reference assembly, confounding studies of genetic variation. Using comparative sequence-based approaches, we generated a high-quality 6.3-Mbp alternate reference assembly of an inverted Chromosome 8p23.1 haplotype. Comparison with nonhuman primates reveals a 746-kbp duplicative transposition and two separate inversion events that arose in the last million years of human evolution. The breakpoints associated with these rearrangements map to an ape-specific interchromosomal core duplicon that clusters at sites of evolutionary inversion (P = 7.8 × 10-5). Refinement of microdeletion breakpoints identifies a subgroup of patients that map to the same interchromosomal core involved in the evolutionary formation of the duplication blocks. Our results define a higher-order genomic instability element that has shaped the structure of specific chromosomes during primate evolution contributing to rearrangements associated with inversion and disease.

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.

Functional significance of single nucleotide polymorphisms in the lactase gene in diverse US patients and evidence for a novel lactase persistence allele at -13909 in those of European ancestry.

OBJECTIVES: Recent data from mainly homogeneous European and African populations implicate a 140-bp region 5' to the transcriptional start site of LCT (the lactase gene) as a regulatory site for lactase persistence and nonpersistence. Because there are no studies of US nonhomogeneous populations, we performed genotype/phenotype analysis of the -13910 and -22018 LCT single nucleotide polymorphisms (SNPs) in New England children, mostly of European ancestry. METHODS: Duodenal biopsies were processed for disaccharidase activities, RNA quantification by reverse transcription polymerase chain reaction (RT-PCR), allelic expression ratios by PCR, and genotyping and SNP analysis. Results were compared with clinical information. RESULTS: Lactase activity and mRNA levels, and sucrase-to-lactase ratios of enzyme activity and mRNA, showed robust correlations with genotype. None of the other LCT SNPs showed as strong a correlation with enzyme or mRNA levels as did -13910. Data were consistent, with the -13910 being the causal sequence variant instead of -22018. Four individuals heterozygous for -13910T/C had allelic expression patterns similar to individuals with -13910C/C genotypes; of these, 2 showed equal LCT expression from the 2 alleles and a novel variant (-13909C>A) associated with lactase persistence. CONCLUSIONS: The identification of -13910C/C genotype is likely to predict lactase nonpersistence, consistent with prior published studies. A -13910T/T genotype will frequently, but not perfectly, predict lactase persistence in this mixed European-ancestry population; a -13910T/C genotype will not predict the phenotype. A long, rare haplotype in 2 individuals with -13910T/C genotype but equal allele-specific expression contains a novel lactase persistence allele present at -13909.

Genome-wide association study and admixture mapping reveal new loci associated with total IgE levels in Latinos.

BACKGROUND: IgE is a key mediator of allergic inflammation, and its levels are frequently increased in patients with allergic disorders. OBJECTIVE: We sought to identify genetic variants associated with IgE levels in Latinos. METHODS: We performed a genome-wide association study and admixture mapping of total IgE levels in 3334 Latinos from the Genes-environments & Admixture in Latino Americans (GALA II) study. Replication was evaluated in 454 Latinos, 1564 European Americans, and 3187 African Americans from independent studies. RESULTS: We confirmed associations of 6 genes identified by means of previous genome-wide association studies and identified a novel genome-wide significant association of a polymorphism in the zinc finger protein 365 gene (ZNF365) with total IgE levels (rs200076616, P = 2.3 × 10(-8)). We next identified 4 admixture mapping peaks (6p21.32-p22.1, 13p22-31, 14q23.2, and 22q13.1) at which local African, European, and/or Native American ancestry was significantly associated with IgE levels. The most significant peak was 6p21.32-p22.1, where Native American ancestry was associated with lower IgE levels (P = 4.95 × 10(-8)). All but 22q13.1 were replicated in an independent sample of Latinos, and 2 of the peaks were replicated in African Americans (6p21.32-p22.1 and 14q23.2). Fine mapping of 6p21.32-p22.1 identified 6 genome-wide significant single nucleotide polymorphisms in Latinos, 2 of which replicated in European Americans. Another single nucleotide polymorphism was peak-wide significant within 14q23.2 in African Americans (rs1741099, P = 3.7 × 10(-6)) and replicated in non-African American samples (P = .011). CONCLUSION: We confirmed genetic associations at 6 genes and identified novel associations within ZNF365, HLA-DQA1, and 14q23.2. Our results highlight the importance of studying diverse multiethnic populations to uncover novel loci associated with total IgE levels.

Whole-genome sequencing of individuals from a founder population identifies candidate genes for asthma.

Asthma is a complex genetic disease caused by a combination of genetic and environmental risk factors. We sought to test classes of genetic variants largely missed by genome-wide association studies (GWAS), including copy number variants (CNVs) and low-frequency variants, by performing whole-genome sequencing (WGS) on 16 individuals from asthma-enriched and asthma-depleted families. The samples were obtained from an extended 13-generation Hutterite pedigree with reduced genetic heterogeneity due to a small founding gene pool and reduced environmental heterogeneity as a result of a communal lifestyle. We sequenced each individual to an average depth of 13-fold, generated a comprehensive catalog of genetic variants, and tested the most severe mutations for association with asthma. We identified and validated 1960 CNVs, 19 nonsense or splice-site single nucleotide variants (SNVs), and 18 insertions or deletions that were out of frame. As follow-up, we performed targeted sequencing of 16 genes in 837 cases and 540 controls of Puerto Rican ancestry and found that controls carry a significantly higher burden of mutations in IL27RA (2.0% of controls; 0.23% of cases; nominal p = 0.004; Bonferroni p = 0.21). We also genotyped 593 CNVs in 1199 Hutterite individuals. We identified a nominally significant association (p = 0.03; Odds ratio (OR) = 3.13) between a 6 kbp deletion in an intron of NEDD4L and increased risk of asthma. We genotyped this deletion in an additional 4787 non-Hutterite individuals (nominal p = 0.056; OR = 1.69). NEDD4L is expressed in bronchial epithelial cells, and conditional knockout of this gene in the lung in mice leads to severe inflammation and mucus accumulation. Our study represents one of the early instances of applying WGS to complex disease with a large environmental component and demonstrates how WGS can identify risk variants, including CNVs and low-frequency variants, largely untested in GWAS.

Properties and rates of germline mutations in humans.

All genetic variation arises via new mutations; therefore, determining the rate and biases for different classes of mutation is essential for understanding the genetics of human disease and evolution. Decades of mutation rate analyses have focused on a relatively small number of loci because of technical limitations. However, advances in sequencing technology have allowed for empirical assessments of genome-wide rates of mutation. Recent studies have shown that 76% of new mutations originate in the paternal lineage and provide unequivocal evidence for an increase in mutation with paternal age. Although most analyses have focused on single nucleotide variants (SNVs), studies have begun to provide insight into the mutation rate for other classes of variation, including copy number variants (CNVs), microsatellites, and mobile element insertions (MEIs). Here, we review the genome-wide analyses for the mutation rate of several types of variants and suggest areas for future research.