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.

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.

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.

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.

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.

Diagnostic Utility of Next-Generation Sequencing for Disorders of Somatic Mosaicism: A Five-Year Cumulative Cohort.

Disorders of somatic mosaicism (DoSM) are a diverse group of syndromic and non-syndromic conditions caused by mosaic variants in genes that regulate cell survival and proliferation. Despite overlap in gene space and technical requirements, few clinical labs specialize in DoSM compared to oncology. We adapted a high-sensitivity next-generation sequencing cancer assay for DoSM in 2014. Some 343 individuals have been tested over the past 5 years, 58% of which had pathogenic and likely pathogenic (P/LP) findings, for a total of 206 P/LP variants in 22 genes. Parameters associated with the high diagnostic yield were: (1) deep sequencing (∼2,000× coverage), (2) a broad gene set, and (3) testing affected tissues. Fresh and formalin-fixed paraffin embedded tissues performed equivalently for identification of P/LP variants (62% and 71% of individuals, respectively). Comparing cultured fibroblasts to skin biopsies suggested that culturing might boost the allelic fraction of variants that confer a growth advantage, specifically gain-of-function variants in PIK3CA. Buccal swabs showed high diagnostic sensitivity in case subjects where disease phenotypes manifested in the head or brain. Peripheral blood was useful as an unaffected comparator tissue to determine somatic versus constitutional origin but had poor diagnostic sensitivity. Descriptions of all tested individuals, specimens, and P/LP variants included in this cohort are available to further the study of the DoSM population.

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.

KRIT1-positive hyperkeratotic cutaneous capillary venous malformation.

Cerebral cavernous malformations (CCM) may present in sporadic or familial forms, with different cutaneous manifestations including deep blue nodules, capillary malformations, and hyperkeratotic cutaneous capillary venous malformations (HCCVM). We report the case of an infant with a KRIT1-positive HCCVM associated with familial CCM. Moreover, histopathology showed positive immunohistochemical stain with GLUT1, further expanding the differential diagnosis of GLUT1-positive vascular anomalies.

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.

Novel morphologic findings in PLAG1-rearranged soft tissue tumors.

Oncogenesis in PLAG1-rearranged tumors often results from PLAG1 transcription factor overexpression driven by promoter-swapping between constitutively expressed fusion partners. PLAG1-rearranged tumors demonstrate diverse morphologies. This study adds to this morphologic heterogeneity by introducing two tumors with PLAG1 rearrangements that display distinct histologic features. The first arose in the inguinal region of a 3-year-old, appeared well-circumscribed with a multinodular pattern, and harbored two fusions: ZFHX4-PLAG1 and CHCHD7-PLAG1. The second arose in the pelvic cavity of a 15-year-old girl, was extensively infiltrative and vascularized with an adipocytic component, and demonstrated a COL3A1-PLAG1 fusion. Both showed low-grade cytomorphology, scarce mitoses, no necrosis, and expression of CD34 and desmin. The ZFHX4-/CHCHD7-PLAG1-rearranged tumor showed no evidence of recurrence after 5 months. By contrast, the COL3A1-PLAG1-rearranged tumor quickly recurred following primary excision with positive margins; subsequent re-excision with adjuvant chemotherapy resulted in no evidence of recurrence after 2 years. While both tumors show overlap with benign and malignant fibroblastic and fibrovascular neoplasms, they also display divergent features. These cases highlight the importance of appropriate characterization in soft tissue tumors with unusual clinical and histologic characteristics.

Gastroblastoma with a novel EWSR1-CTBP1 fusion presenting in adolescence.

Gastroblastomas are rare tumors with a biphasic epithelioid/spindle cell morphology that typically present in early adulthood and have recurrent MALAT1-GLI1 fusions. We describe an adolescent patient with Wiskott-Aldrich syndrome who presented with a large submucosal gastric tumor with biphasic morphology. Despite histologic features consistent with gastroblastoma, a MALAT1-GLI1 fusion was not found in this patient's tumor; instead, comprehensive molecular profiling identified a novel EWSR1-CTBP1 fusion and no other significant genetic alterations. The tumor also overexpressed NOTCH and FGFR by RNA profiling. The novel fusion and expression profile suggest a role for epithelial-mesenchymal transition in this tumor, with potential implications for the pathogenesis of biphasic gastric tumors such as gastroblastoma.

Discovery of clinically relevant fusions in pediatric cancer.

BACKGROUND: Pediatric cancers typically have a distinct genomic landscape when compared to adult cancers and frequently carry somatic gene fusion events that alter gene expression and drive tumorigenesis. Sensitive and specific detection of gene fusions through the analysis of next-generation-based RNA sequencing (RNA-Seq) data is computationally challenging and may be confounded by low tumor cellularity or underlying genomic complexity. Furthermore, numerous computational tools are available to identify fusions from supporting RNA-Seq reads, yet each algorithm demonstrates unique variability in sensitivity and precision, and no clearly superior approach currently exists. To overcome these challenges, we have developed an ensemble fusion calling approach to increase the accuracy of identifying fusions. RESULTS: Our Ensemble Fusion (EnFusion) approach utilizes seven fusion calling algorithms: Arriba, CICERO, FusionMap, FusionCatcher, JAFFA, MapSplice, and STAR-Fusion, which are packaged as a fully automated pipeline using Docker and Amazon Web Services (AWS) serverless technology. This method uses paired end RNA-Seq sequence reads as input, and the output from each algorithm is examined to identify fusions detected by a consensus of at least three algorithms. These consensus fusion results are filtered by comparison to an internal database to remove likely artifactual fusions occurring at high frequencies in our internal cohort, while a "known fusion list" prevents failure to report known pathogenic events. We have employed the EnFusion pipeline on RNA-Seq data from 229 patients with pediatric cancer or blood disorders studied under an IRB-approved protocol. The samples consist of 138 central nervous system tumors, 73 solid tumors, and 18 hematologic malignancies or disorders. The combination of an ensemble fusion-calling pipeline and a knowledge-based filtering strategy identified 67 clinically relevant fusions among our cohort (diagnostic yield of 29.3%), including RBPMS-MET, BCAN-NTRK1, and TRIM22-BRAF fusions. Following clinical confirmation and reporting in the patient's medical record, both known and novel fusions provided medically meaningful information. CONCLUSIONS: The EnFusion pipeline offers a streamlined approach to discover fusions in cancer, at higher levels of sensitivity and accuracy than single algorithm methods. Furthermore, this method accurately identifies driver fusions in pediatric cancer, providing clinical impact by contributing evidence to diagnosis and, when appropriate, indicating targeted therapies.

Somatic PIK3R1 variation as a cause of vascular malformations and overgrowth.

PURPOSE: Somatic activating variants in the PI3K-AKT pathway cause vascular malformations with and without overgrowth. We previously reported an individual with capillary and lymphatic malformation harboring a pathogenic somatic variant in PIK3R1, which encodes three PI3K complex regulatory subunits. Here, we investigate PIK3R1 in a large cohort with vascular anomalies and identify an additional 16 individuals with somatic mosaic variants in PIK3R1. METHODS: Affected tissue from individuals with vascular lesions and overgrowth recruited from a multisite collaborative network was studied. Next-generation sequencing targeting coding regions of cell-signaling and cancer-associated genes was performed followed by assessment of variant pathogenicity. RESULTS: The phenotypic and variant spectrum associated with somatic variation in PIK3R1 is reported herein. Variants occurred in the inter-SH2 or N-terminal SH2 domains of all three PIK3R1 protein products. Phenotypic features overlapped those of the PIK3CA-related overgrowth spectrum (PROS). These overlapping features included mixed vascular malformations, sandal toe gap deformity with macrodactyly, lymphatic malformations, venous ectasias, and overgrowth of soft tissue or bone. CONCLUSION: Somatic PIK3R1 variants sharing attributes with cancer-associated variants cause complex vascular malformations and overgrowth. The PIK3R1-associated phenotypic spectrum overlaps with PROS. These data extend understanding of the diverse phenotypic spectrum attributable to genetic variation in the PI3K-AKT pathway.

Targeted Therapy in a Young Adult With a Novel Epithelioid Tumor Driven by a PRRC2B-ALK Fusion.

This case report describes an 18-year-old woman with an unusual epithelioid tumor of the omentum with a novel PRRC2B-ALK fusion. Although the atypical pathologic features raised significant diagnostic challenges, expression of CD30 on tumor cells and detection of an ALK rearrangement provided critical information for selecting targeted therapy in a patient not suitable for surgical resection. Despite an initially promising therapeutic response, the patient died. The efficacy of treatment was confirmed by the lack of viable tumor cells at autopsy. This case highlights the role of timely targeted therapy in patients with rare tumors and novel actionable molecular targets.

Genomic heterogeneity of ALK fusion breakpoints in non-small-cell lung cancer.

In lung adenocarcinoma, canonical EML4-ALK inversion results in a fusion protein with a constitutively active ALK kinase domain. Evidence of ALK rearrangement occurs in a minority (2-7%) of lung adenocarcinoma, and only ~60% of these patients will respond to targeted ALK inhibition by drugs such as crizotinib and ceritinib. Clinically, targeted anti-ALK therapy is often initiated based on evidence of an ALK genomic rearrangement detected by fluorescence in situ hybridization (FISH) of interphase cells in formalin-fixed, paraffin-embedded tissue sections. At the genomic level, however, ALK rearrangements are heterogeneous, with multiple potential breakpoints in EML4, and alternate fusion partners. Using next-generation sequencing of DNA and RNA together with ALK immunohistochemistry, we comprehensively characterized genomic breakpoints in 33 FISH-positive lung adenocarcinomas. Of these 33 cases, 29 (88%) had detectable DNA level ALK rearrangements involving EML4, KIF5B, or non-canonical partners including ASXL2, ATP6V1B1, PRKAR1A, and SPDYA. A subset of 12 cases had material available for RNA-Seq. Of these, eight of eight (100%) cases with DNA rearrangements showed ALK fusion transcripts from RNA-Seq; three of four cases (75%) without detectable DNA rearrangements were similarly negative by RNA-Seq, and one case was positive by RNA-Seq but negative by DNA next-generation sequencing. By immunohistochemistry, 17 of 19 (89%) tested cases were clearly positive for ALK protein expression; the remaining cases had no detectable DNA level rearrangement or had a non-canonical rearrangement not predicted to form a fusion protein. Survival analysis of patients treated with targeted ALK inhibitors demonstrates a significant difference in mean survival between patients with next-generation sequencing confirmed EML4-ALK rearrangements, and those without (20.6 months vs 5.4 months, P<0.01). Together, these data demonstrate abundant genomic heterogeneity among ALK-rearranged lung adenocarcinoma, which may account for differences in treatment response with targeted ALK inhibitors.

Routine use of clinical exome-based next-generation sequencing for evaluation of patients with thrombotic microangiopathies.

Next-generation sequencing is increasingly used for clinical evaluation of patients presenting with thrombotic microangiopathies because it allows for simultaneous interrogation of multiple complement and coagulation pathway genes known to be associated with disease. However, the diagnostic yield is undefined in routine clinical practice. Historic studies relied on case-control cohorts, did not apply current guidelines for variant pathogenicity assessment, and used targeted gene enrichment combined with next-generation sequencing. A clinically enhanced exome, targeting ~54 Mb, was sequenced for 73 patients. Variant analysis and interpretation were performed on genes with biological relevance in thrombotic microangiopathy (C3,CD46, CFB, CFH, CFI, DGKE, and THBD). CFHR3-CFHR1 deletion status was also assessed using multiplex ligation-dependent probe amplification. Variants were classified using American College of Medical Genetics and Genomics guidelines. We identified 5 unique novel and 14 unique rare variants in 25% (18/73) of patients, including a total of 5 pathogenic, 4 likely pathogenic, and 15 variants of uncertain clinical significance. Nine patients had homozygous deletions in CFHR3-CFHR1. The diagnostic yield, defined as the presence of a pathogenic variant, likely pathogenic variant or homozygous deletion of CFHR3-CFHR1, was 25% for all patients tested. Variants of uncertain clinical significance were identified in 21% (15/73) of patients.These results illustrate the expected diagnositic yield in the setting of thrombotic microangiopathies through the application of standardized variant interpretation, and highlight the utility of such an approach. Sequencing a clinically enhanced exome to enable targeted, disease-specific variant analysis is a viable approach. The moderate rate of variants of uncertain clinical significance highlights the paucity of data surrounding the variants in our cohort and illustrates the need for expanded variant curation resources to aid in thrombotic microangiopathy-related disease variant classification.

Spectrum of mutations in leiomyosarcomas identified by clinical targeted next-generation sequencing.

Recurrent genomic mutations in uterine and non-uterine leiomyosarcomas have not been well established. Using a next generation sequencing (NGS) panel of common cancer-associated genes, 25 leiomyosarcomas arising from multiple sites were examined to explore genetic alterations, including single nucleotide variants (SNV), small insertions/deletions (indels), and copy number alterations (CNA). Sequencing showed 86 non-synonymous, coding region somatic variants within 151 gene targets in 21 cases, with a mean of 4.1 variants per case; 4 cases had no putative mutations in the panel of genes assayed. The most frequently altered genes were TP53 (36%), ATM and ATRX (16%), and EGFR and RB1 (12%). CNA were identified in 85% of cases, with the most frequent copy number losses observed in chromosomes 10 and 13 including PTEN and RB1; the most frequent gains were seen in chromosomes 7 and 17. Our data show that deletions in canonical cancer-related genes are common in leiomyosarcomas. Further, the spectrum of gene mutations observed shows that defects in DNA repair and chromosomal maintenance are central to the biology of leiomyosarcomas, and that activating mutations observed in other common cancer types are rare in leiomyosarcomas.