Structural and Dynamic Analyses of Pathogenic Variants in PIK3R1 Reveal a Shared Mechanism Associated among Cancer, Undergrowth, and Overgrowth Syndromes.

The PI3K enzymes modify phospholipids to regulate cell growth and differentiation. Somatic variants in PI3K are recurrent in cancer and drive a proliferative phenotype. Somatic mosaicism of PIK3R1 and PIK3CA are associated with vascular anomalies and overgrowth syndromes. Germline PIK3R1 variants are associated with varying phenotypes, including immunodeficiency or facial dysmorphism with growth delay, lipoatrophy, and insulin resistance associated with SHORT syndrome. There has been limited study of the molecular mechanism to unify our understanding of how variants in PIK3R1 drive both undergrowth and overgrowth phenotypes. Thus, we compiled genomic variants from cancer and rare vascular anomalies and sought to interpret their effects using an unbiased physics-based simulation approach for the protein complex. We applied molecular dynamics simulations to mechanistically understand how genetic variants affect PIK3R1 and its interactions with PIK3CA. Notably, iSH2 genetic variants associated with undergrowth destabilize molecular interactions with the PIK3CA receptor binding domain in simulations, which is expected to decrease activity. On the other hand, overgrowth and cancer variants lead to loss of inhibitory interactions in simulations, which is expected to increase activity. We find that all disease variants display dysfunctions on either structural characteristics or intermolecular interaction energy. Thus, this comprehensive characterization of novel mosaic somatic variants associated with two opposing phenotypes has mechanistic importance and biomedical relevance and may aid in future therapeutic developments.

A LINE-1 mediated deletion resulting in germline retinoblastoma predisposition.

Retinoblastoma is an ocular cancer associated with genomic variation in the RB1 gene. In individuals with bilateral retinoblastoma, a germline variant in RB1 is identified in virtually all cases. We describe herein an individual with bilateral retinoblastoma for whom multiple clinical lab assays performed by outside commercial laboratories failed to identify a germline RB1 variant. Paired tumor/normal exome sequencing, long-read whole genome sequencing, and long-read isoform sequencing was performed on a translational research basis ultimately identified a germline likely de novo Long Interspersed Nuclear Element (LINE)-1 mediated deletion resulting in a premature stop of translation of RB1 as the underlying genetic cause of retinoblastoma in this individual. Based on these research findings, the LINE-1 mediated deletion was confirmed via Sanger sequencing in our clinical laboratory, and results were reported in the patient's medical record to allow for appropriate genetic counseling.

Full-length isoform concatenation sequencing to resolve cancer transcriptome complexity.

BACKGROUND: Cancers exhibit complex transcriptomes with aberrant splicing that induces isoform-level differential expression compared to non-diseased tissues. Transcriptomic profiling using short-read sequencing has utility in providing a cost-effective approach for evaluating isoform expression, although short-read assembly displays limitations in the accurate inference of full-length transcripts. Long-read RNA sequencing (Iso-Seq), using the Pacific Biosciences (PacBio) platform, can overcome such limitations by providing full-length isoform sequence resolution which requires no read assembly and represents native expressed transcripts. A constraint of the Iso-Seq protocol is due to fewer reads output per instrument run, which, as an example, can consequently affect the detection of lowly expressed transcripts. To address these deficiencies, we developed a concatenation workflow, PacBio Full-Length Isoform Concatemer Sequencing (PB_FLIC-Seq), designed to increase the number of unique, sequenced PacBio long-reads thereby improving overall detection of unique isoforms. In addition, we anticipate that the increase in read depth will help improve the detection of moderate to low-level expressed isoforms. RESULTS: In sequencing a commercial reference (Spike-In RNA Variants; SIRV) with known isoform complexity we demonstrated a 3.4-fold increase in read output per run and improved SIRV recall when using the PB_FLIC-Seq method compared to the same samples processed with the Iso-Seq protocol. We applied this protocol to a translational cancer case, also demonstrating the utility of the PB_FLIC-Seq method for identifying differential full-length isoform expression in a pediatric diffuse midline glioma compared to its adjacent non-malignant tissue. Our data analysis revealed increased expression of extracellular matrix (ECM) genes within the tumor sample, including an isoform of the Secreted Protein Acidic and Cysteine Rich (SPARC) gene that was expressed 11,676-fold higher than in the adjacent non-malignant tissue. Finally, by using the PB_FLIC-Seq method, we detected several cancer-specific novel isoforms. CONCLUSION: This work describes a concatenation-based methodology for increasing the number of sequenced full-length isoform reads on the PacBio platform, yielding improved discovery of expressed isoforms. We applied this workflow to profile the transcriptome of a pediatric diffuse midline glioma and adjacent non-malignant tissue. Our findings of cancer-specific novel isoform expression further highlight the importance of long-read sequencing for characterization of complex tumor transcriptomes.

Papillary Hemangioma Harbors Somatic GNA11 and GNAQ Mutations.

Papillary hemangioma (PH) is a small, primarily dermal lesion occurring predominantly in the head and neck in both children and adults. Its signature characteristics are dilated thin-walled channels containing papillary clusters of mainly capillary-sized vessels and endothelial cytoplasmic eosinophilic inclusions. Given certain histopathologic similarities to congenital hemangioma which harbor mutations in GNAQ and GNA11 , we investigated whether similar mutations are present in PH. Seven PH specimens were studied. All presented in the first 4 years of life, with one being noted at birth. With the exception of one lesion, all were in the head and neck. Lesions were bluish and ranged in size from 0.5 to 2.8 cm. Four samples had GNA11 p.Q209L and 3 had GNAQ p.Q209L missense mutations. Mutations in GNA11 and GNAQ are associated with other types of somatic vascular lesions including capillary malformation, congenital hemangioma, anastomosing hemangioma, thrombotic anastomosing hemangioma, and hepatic small cell neoplasm. Shared mutations in GNA11 and GNAQ may account for some overlapping clinical and pathologic features in these entities, perhaps explicable by the timing of the mutation or influence of the germline phenotype.

Exome sequencing in the pediatric neuromuscular clinic leads to more frequent diagnosis of both neuromuscular and neurodevelopmental conditions.

INTRODUCTION/AIMS: Exome sequencing (ES) has proven to be a valuable diagnostic tool for neuromuscular disorders, which often pose a diagnostic challenge. The aims of this study were to investigate the clinical outcomes associated with utilization of ES in the pediatric neuromuscular clinic and to determine if specific phenotypic features or abnormal neurodiagnostic tests were predictive of a diagnostic result. METHODS: This was a retrospective medical record review of 76 pediatric neuromuscular clinic patients who underwent ES. Based upon clinical assessment prior to ES, patients were divided into two groups: affected by neuromuscular (n = 53) or non-neuromuscular (n = 23) syndromes. RESULTS: A diagnosis was made in 28/76 (36.8%), with 29 unique disorders identified. In the neuromuscular group, a neuromuscular condition was confirmed in 78% of those receiving a genetic diagnosis. Early age of symptom onset was associated with a significantly higher diagnostic yield. The most common reason neuromuscular diagnoses were not detected on prior testing was due to causative genes not being present on disease-specific panels. Changes to medical care were made in 57% of individuals receiving a diagnosis on ES. DISCUSSION: These data further support ES as a powerful diagnostic tool in the pediatric neuromuscular clinic and highlight the advantages of ES over gene panels, including the ability to identify diagnoses regardless of etiology, identify genes newly associated with disease, and identify multiple confounding diagnoses. Rapid and accurate diagnosis by ES can not only end the patient's diagnostic odyssey, but often impacts patients' medical management and genetic counseling of families.

Expanding the Clinical Utility of Targeted RNA Sequencing Panels beyond Gene Fusions to Complex, Intragenic Structural Rearrangements.

Gene fusions are a form of structural rearrangement well established as driver events in pediatric and adult cancers. The identification of such events holds clinical significance in the refinement, prognostication, and provision of treatment in cancer. Structural rearrangements also extend beyond fusions to include intragenic rearrangements, such as internal tandem duplications (ITDs) or exon-level deletions. These intragenic events have been increasingly implicated as cancer-promoting events. However, the detection of intragenic rearrangements may be challenging to resolve bioinformatically with short-read sequencing technologies and therefore may not be routinely assessed in panel-based testing. Within an academic clinical laboratory, over three years, a total of 608 disease-involved samples (522 hematologic malignancy, 86 solid tumors) underwent clinical testing using Anchored Multiplex PCR (AMP)-based RNA sequencing. Hematologic malignancies were evaluated using a custom Pan-Heme 154 gene panel, while solid tumors were assessed using a custom Pan-Solid 115 gene panel. Gene fusions, ITDs, and intragenic deletions were assessed for diagnostic, prognostic, or therapeutic significance. When considering gene fusions alone, we report an overall diagnostic yield of 36% (37% hematologic malignancy, 41% solid tumors). When including intragenic structural rearrangements, the overall diagnostic yield increased to 48% (48% hematologic malignancy, 45% solid tumor). We demonstrate the clinical utility of reporting structural rearrangements, including gene fusions and intragenic structural rearrangements, using an AMP-based RNA sequencing panel.

Enhancing the didactic learning experience for Laboratory Genetics and Genomics fellows through a multi-institutional lecture series.

As a result of the pandemic, the traditional in-person didactic lecture model was adapted to a virtual learning approach. Our Laboratory Genetics and Genomics fellowship program at Nationwide Children's hospital took advantage of this opportunity to organize a multi-institutional Fellow's Conference to educate fellows from different programs on a wide range of medical genetics topics. We describe our approach of developing this lecture series utilizing subject-matter experts across institutions. In addition, we discuss the value of such an approach in reducing the amount of time individual institutions spend creating and providing didactic content for their small number of learners. Our experience could serve as a model for other educators and program directors, including genetic counseling program directors, to develop multi-institutional collaborations for didactic learning.

A novel transcriptional signature identifies T-cell infiltration in high-risk paediatric cancer.

BACKGROUND: Molecular profiling of the tumour immune microenvironment (TIME) has enabled the rational choice of immunotherapies in some adult cancers. In contrast, the TIME of paediatric cancers is relatively unexplored. We speculated that a more refined appreciation of the TIME in childhood cancers, rather than a reliance on commonly used biomarkers such as tumour mutation burden (TMB), neoantigen load and PD-L1 expression, is an essential prerequisite for improved immunotherapies in childhood solid cancers. METHODS: We combined immunohistochemistry (IHC) with RNA sequencing and whole-genome sequencing across a diverse spectrum of high-risk paediatric cancers to develop an alternative, expression-based signature associated with CD8+ T-cell infiltration of the TIME. Furthermore, we explored transcriptional features of immune archetypes and T-cell receptor sequencing diversity, assessed the relationship between CD8+ and CD4+ abundance by IHC and deconvolution predictions and assessed the common adult biomarkers such as neoantigen load and TMB. RESULTS: A novel 15-gene immune signature, Immune Paediatric Signature Score (IPASS), was identified. Using this signature, we estimate up to 31% of high-risk cancers harbour infiltrating T-cells. In addition, we showed that PD-L1 protein expression is poorly correlated with PD-L1 RNA expression and TMB and neoantigen load are not predictive of T-cell infiltration in paediatrics. Furthermore, deconvolution algorithms are only weakly correlated with IHC measurements of T-cells. CONCLUSIONS: Our data provides new insights into the variable immune-suppressive mechanisms dampening responses in paediatric solid cancers. Effective immune-based interventions in high-risk paediatric cancer will require individualised analysis of the TIME.

CHEK2 Alterations in Pediatric Malignancy: A Single-Institution Experience.

BACKGROUND: Approximately 10% of pediatric malignancies are secondary to germline alterations in cancer-predisposing genes. Checkpoint kinase 2 (CHEK2) germline loss-of-function variants have been reported in pediatric cancer patients, but clinical phenotypes and outcomes are poorly described. We present our single-institution experience of pediatric oncology patients with CHEK2 germline alterations, including clinical presentations and outcomes. METHODS: Pediatric oncology patients with CHEK2 germline alterations were identified among those assessed by clinical or translational research at the Institute for Genomic Medicine at Nationwide Children's Hospital. A chart review of disease course was conducted on identified patients. RESULTS: We identified 6 patients with germline CHEK2 variants from a cohort of 300 individuals, including 1 patient with concurrent presentation of Burkitt lymphoma and neuroblastoma, 3 patients with brain tumors, 1 patient with Ewing sarcoma, and 1 patient with myelodysplastic syndrome. Three patients had a family history of malignancies. Four patients were in remission; one was undergoing treatment; one patient had developed treatment-related meningiomas. We review prior data regarding CHEK2 variants in this population, challenges associated with variant interpretation, and genetic counseling for individuals with CHEK2 variants. CONCLUSIONS: CHEK2 germline loss-of-function alterations occur in patients with a variety of pediatric tumors. Larger multicenter studies will improve our understanding of the incidence, phenotype, and molecular biology of CHEK2 germline variants in pediatric cancers.

EGFR internal tandem duplications in fusion-negative congenital and neonatal spindle cell tumors.

Next-generation sequencing (NGS) assays can sensitively detect somatic variation, and increasingly can enable the identification of complex structural rearrangements. A subset of infantile spindle cell sarcomas, particularly congenital mesoblastic nephromas with classic or mixed histology, have structural rearrangement in the form of internal tandem duplications (ITD) involving EGFR. We performed prospective analysis to identify EGFR ITD through clinical or research studies, as well as retrospective analysis to quantify the frequency of EGFR ITD in pediatric sarcomas. Within our institution, three tumors with EGFR ITD were prospectively identified, all occurring in patients less than 1 year of age at diagnosis, including two renal tumors and one mediastinal soft tissue tumor. These three cases exhibited both cellular and mixed cellular and classic histology. All patients had no evidence of disease progression off therapy, despite incomplete resection. To extend our analysis and quantify the frequency of EGFR ITD in pediatric sarcomas, we retrospectively analyzed a cohort of tumors (n = 90) that were previously negative for clinical RT-PCR-based fusion testing. We identified EGFR ITD in three analyzed cases, all in patients less than 1 year of age (n = 18; 3/18, 17%). Here we expand the spectrum of tumors with EGFR ITD to congenital soft tissue tumors and report an unusual example of an EGFR ITD in a tumor with cellular congenital mesoblastic nephroma histology. We also highlight the importance of appropriate test selection and bioinformatic analysis for identification of this genomic alteration that is unexpectedly common in congenital and infantile spindle cell tumors.

Infantile metastatic ependymoma with a novel molecular profile and favorable outcome to intensive chemotherapy without irradiation: Case-based review.

Ependymal tumors are the third most common brain tumor under 14 years old. Even though metastatic disease is a rare event, it affects mostly young children and carries an adverse prognosis. The factors associated with dissemination and the best treatment approach have not yet been established and there is limited published data on how to manage metastatic disease, especially in patients under 3 years of age. We provide a review of the literature on clinical characteristics and radiation-sparing treatments for metastatic ependymoma in children under 3 years of age treated. The majority (73%) of the identified cases were above 12 months old and had the PF as the primary site at diagnosis. Chemotherapy-based approaches, in different regimens, were used with radiation reserved for progression or relapse. The prognosis varied among the studies, with an average of 50%-58% overall survival. This study also describes the case of a 7-month-old boy with metastatic posterior fossa (PF) ependymoma, for whom we identified a novel SPECC1L-RAF1 gene fusion using a patient-centric comprehensive molecular profiling protocol. The patient was successfully treated with intensive induction chemotherapy followed by high-dose chemotherapy and autologous hematopoietic progenitor cell rescue (AuHSCR). Currently, the patient is in continuous remission 5 years after his diagnosis, without radiation therapy. The understanding of the available therapeutic approaches may assist physicians in their management of such patients. This report also opens the perspective of newly identified molecular alterations in metastatic ependymomas that might drive more chemo-sensitive tumors.

Assessments of Somatic Variant Classification Using the Association for Molecular Pathology/American Society of Clinical Oncology/College of American Pathologists Guidelines: A Report from the Association for Molecular Pathology.

To assess the clinical implementation of the 2017 Standards and Guidelines for the Interpretation and Reporting of Sequence Variants in Cancer: A Joint Consensus Recommendation of the Association for Molecular Pathology, American Society of Clinical Oncology, and College of American Pathologists, identify content that may result in classification inconsistencies, and evaluate implementation barriers, an Association for Molecular Pathology Working Group conducted variant interpretation challenges and a guideline implementation survey. A total of 134 participants participated in the variant interpretation challenges, consisting of 11 variants in four cancer cases. Results demonstrate 86% (range, 54% to 94%) of the respondents correctly classified clinically significant variants, variants of uncertain significance, and benign/likely benign variants; however, only 59% (range, 39% to 84%) of responses agreed with the working group's consensus intended responses regarding both tiers and categories of clinical significance. In the implementation survey, 71% (157/220) of respondents have implemented the 2017 guidelines for variant classification and reporting either with or without modifications. Collectively, this study demonstrates that, although they may not yet be optimized, the 2017 guideline recommendations are being adopted for standardized somatic variant classification. The working group identified significant areas for future guideline improvement, including the need for a more granular and comprehensive classification system and education resources to meet the growing needs of both laboratory professionals and medical oncologists.

A comparative analysis of RAS variants in patients with disorders of somatic mosaicism.

PURPOSE: RAS genes (HRAS, KRAS, and NRAS) are commonly found to be mutated in cancers, and activating RAS variants are also found in disorders of somatic mosaicism (DoSM). A survey of the mutational spectrum of RAS variants in DoSM has not been performed. METHODS: A total of 938 individuals with suspected DoSM underwent high-sensitivity clinical next-generation sequencing-based testing. We investigated the mutational spectrum and genotype-phenotype associations of mosaic RAS variants. RESULTS: In this article, we present a series of individuals with DoSM with RAS variants. Classic hotspots, including Gly12, Gly13, and Gln61 constituted the majority of RAS variants observed in DoSM. Furthermore, we present 12 individuals with HRAS and KRAS in-frame duplication/insertion (dup/ins) variants in the switch II domain. Among the 18.3% individuals with RAS in-frame dup/ins variants, clinical findings were mainly associated with vascular malformations. Hotspots were associated with a broad phenotypic spectrum, including vascular tumors, vascular malformations, nevoid proliferations, segmental overgrowth, digital anomalies, and combinations of these. The median age at testing was higher and the variant allelic fraction was lower in individuals with in-frame dup/ins variants than those in individuals with mosaic RAS hotspots. CONCLUSION: Our work provides insight into the allelic and clinical heterogeneity of mosaic RAS variants in nonmalignant conditions.

Early-onset hypertension associated with extensive cutaneous capillary malformations harboring postzygotic variants in GNAQ and GNA11.

BACKGROUND AND OBJECTIVES: Cutaneous capillary malformations (CMs) describe a group of vascular birthmarks with heterogeneous presentations. CMs may present as an isolated finding or with other associations, including glaucoma and leptomeningeal angiomatosis (i.e., Sturge-Weber syndrome) or pigmentary birthmarks (i.e., phakomatosis pigmentovascularis). The use of targeted genetic sequencing has revealed that postzygotic somatic variations in GNAQ and GNA11 at codon 183 are associated with CMs. We report five patients with early-onset hypertension and discuss possible pathogenesis of hypertension. METHODS: Twenty-nine patients with CMs, confirmed GNAQ/11 postzygotic variants, and documented past medical history were identified from a multi-institutional vascular anomalies study. Early-onset hypertension was defined as hypertension before the age of 55 years. Clinical data were reviewed for evidence of hypertension, such as documentation of diagnosis or elevated blood pressure measurements. RESULTS: Five of the 29 patients identified as having GNAQ/11 postzygotic variants had documented early-onset hypertension. Three individuals harbored a GNAQ p.R183Q variant, and two individuals harbored a GNA11 p.R183C variant. All individuals had extensive cutaneous CMs involving the trunk and covering 9%-56% of their body surface area. The median age of hypertension diagnosis was 15 years (range 11-24 years), with three individuals having renal abnormalities on imaging. CONCLUSIONS: Early-onset hypertension is associated with extensive CMs harboring somatic variations in GNAQ/11. Here, we expand on the GNAQ/11 phenotype and hypothesize potential mechanisms driving hypertension. We recommend serial blood pressure measurements in patients with extensive CMs on the trunk and extremities to screen for early-onset hypertension.

Pacific Biosciences Fusion and Long Isoform Pipeline for Cancer Transcriptome-Based Resolution of Isoform Complexity.

Genomic profiling using short-read sequencing has utility in detecting disease-associated variation in both DNA and RNA. However, given the frequent occurrence of structural variation in cancer, molecular profiling using long-read sequencing improves the resolution of such events. For example, the Pacific Biosciences long-read RNA-sequencing (Iso-Seq) transcriptome protocol provides full-length isoform characterization, discernment of allelic phasing, and isoform discovery, and identifies expressed fusion partners. The Pacific Biosciences Fusion and Long Isoform Pipeline (PB_FLIP) incorporates a suite of RNA-sequencing software analysis tools and scripts to identify expressed fusion partners and isoforms. In addition, sequencing of a commercial reference (Spike-In RNA Variants) with known isoform complexity was performed and demonstrated high recall of the Iso-Seq and PB_FLIP workflow to benchmark our protocol and analysis performance. This study describes the utility of Iso-Seq and PB_FLIP analysis in improving deconvolution of complex structural variants and isoform detection within an institutional pediatric and adolescent/young adult translational cancer research cohort. The exemplar case studies demonstrate that Iso-Seq and PB_FLIP discover novel expressed fusion partners, resolve complex intragenic alterations, and discriminate between allele-specific expression profiles.

Molecular Heterogeneity in Pediatric Malignant Rhabdoid Tumors in Patients With Multi-Organ Involvement.

Rhabdoid tumors (RTs) of the brain (atypical teratoid/rhabdoid tumor; AT/RT) and extracranial sites (most often the kidney; RTK) are malignant tumors predominantly occurring in children, frequently those with SMARCB1 germline alterations. Here we present data from seven RTs from three pediatric patients who all had multi-organ involvement. The tumors were analyzed using a multimodal molecular approach, which included exome sequencing of tumor and germline comparator and RNA sequencing and DNA array-based methylation profiling of tumors. SMARCB1 germline alterations were identified in all patients and in all tumors. We observed a second hit in SMARCB1 via chr22 loss of heterozygosity. By methylation profiling, all tumors were classified as rhabdoid tumors with a corresponding subclassification within the MYC, TYR, or SHH AT/RT subgroups. Using RNA-seq gene expression clustering, we recapitulated the classification of known AT/RT subgroups. Synchronous brain and kidney tumors from the same patient showed different patterns of either copy number variants, single-nucleotide variants, and/or genome-wide DNA methylation, suggestive of non-clonal origin. Furthermore, we demonstrated that a lung and abdominal metastasis from two patients shared overlapping molecular features with the patient's primary kidney tumor, indicating the likely origin of the metastasis. In addition to the SMARCB1 events, we identified other whole-chromosome events and single-nucleotide variants in tumors, but none were found to be prognostic, diagnostic, or offer therapeutic potential for rhabdoid tumors. While our findings are of biological interest, there may also be clinical value in comprehensive molecular profiling in patients with multiple rhabdoid tumors, particularly given the potential prognostic and therapeutic implications for different rhabdoid tumor subgroups demonstrated in recent clinical trials and other large cohort studies.

Detection of brain somatic variation in epilepsy-associated developmental lesions.

OBJECTIVE: Epilepsy-associated developmental lesions, including malformations of cortical development and low-grade developmental tumors, represent a major cause of drug-resistant seizures requiring surgical intervention in children. Brain-restricted somatic mosaicism has been implicated in the genetic etiology of these lesions; however, many contributory genes remain unidentified. METHODS: We enrolled 50 children who were undergoing epilepsy surgery into a translational research study. Resected tissue was divided for clinical neuropathologic evaluation and genomic analysis. We performed exome and RNA sequencing to identify somatic variation and we confirmed our findings using high-depth targeted DNA sequencing. RESULTS: We uncovered candidate disease-causing somatic variation affecting 28 patients (56%), as well as candidate germline variants affecting 4 patients (8%). In agreement with previous studies, we identified somatic variation affecting solute carrier family 35 member A2 (SLC35A2) and mechanistic target of rapamycin kinase (MTOR) pathway genes in patients with focal cortical dysplasia. Somatic gains of chromosome 1q were detected in 30% (3 of 10) of patients with Type I focal cortical dysplasia (FCD)s. Somatic variation in mitogen-activated protein kinase (MAPK) pathway genes (i.e., fibroblast growth factor receptor 1 [FGFR1], FGFR2, B-raf proto-oncogene, serine/threonine kinase [BRAF], and KRAS proto-oncogene, GTPase [KRAS]) was associated with low-grade epilepsy-associated developmental tumors. RNA sequencing enabled the detection of somatic structural variation that would have otherwise been missed, and which accounted for more than one-half of epilepsy-associated tumor diagnoses. Sampling across multiple anatomic regions revealed that somatic variant allele fractions vary widely within epileptogenic tissue. Finally, we identified putative disease-causing variants in genes not yet associated with focal cortical dysplasia. SIGNIFICANCE: These results further elucidate the genetic basis of structural brain abnormalities leading to focal epilepsy in children and point to new candidate disease genes.