Search

Detection of gene fusions is important for discovery of cancer drivers and clinical oncology testing, but existing software tools for fusion detection usually take hours to run and may fail to find lowly expressed fusions. To overcome these limitations, we developed the Fuzzion2 program, which uses pattern matching to detect known gene fusions in unmapped paired-read RNA-Seq data. Given a set of patterns representing fusion transcript breakpoints, Fuzzion2 finds every read pair matching any of the patterns. Both exact and inexact (fuzzy) matches are detected; the fuzzy matching tolerates variations caused by sequencing errors, SNVs, and indels. By employing a novel index of frequency minimizers, Fuzzion2 needs only minutes to process a sample. We have also developed pipelines to produce patterns for Fuzzion2, from fusion contig sequences, from genomic breakpoints in DNA and RNA, and from fusion protein sequences. To evaluate its applicability in clinical testing, we ran Fuzzion2 on ~2,000 RNA-seq samples profiled by the St. Jude clinical genomics program and confirmed its sensitivity in identifying lowly expressed fusions, such as KIAA1549-BRAF in low-grade glioma, which are frequently missed by commonly used fusion detection programs. Notably, Fuzzion2 detected a subclonal BCR-ABL1 fusion expressed at 1% and 6% of the wild-type BCR and ABL1 transcription level, respectively, in a B-lineage ALL sample that also has an IGH-CRLF fusion. Processing RNA-seq data from BCR-ABL1 cell lines, K562 with p210 fusion and OP1 with p190 fusion, diluted at 1:10, 1:100, and 1:1000 showed that Fuzzion2 can detect the fusion at 1:10-1:100 dilution, achieving a sensitivity 10 times greater than that of other fusion detection programs. We also evaluated the performance of Fuzzion2 for large-scale data mining in a study to compare the prevalence of gene fusions in pediatric versus adult cancers. We assembled a set of 15,474 patterns representing 5,480 fusions identified in the Pediatric Cancer Genome Project, NCI TARGET, clinical sequencing, and the COSMIC database. Fuzzion2 was deployed to the NCI Cancer Genomics Cloud and analyzed 9,464 TCGA RNA-seq samples from adult solid and brain tumors. Processing took an average of 6 minutes at a cost of only US$0.16 per sample. Among the 105 recurrent fusions identified in pediatric cancers, only 11 were also found in adult cancers. These shared fusions can be classified into two categories: 1) gene fusions present in cancers that occur in both children and young adults, e.g., synovial sarcoma, papillary thyroid cancer, and fibrolamellar hepatocellular carcinoma; and 2) kinase fusions involving ABL1, NTRK, and FGFR. Our experience with Fuzzion2 demonstrates that it is a powerful tool for time-critical clinical application and large-scale data mining. It is publicly available at https://github.com/stjude/fuzzion2. Citation Format: Stephen V. Rice, Michael N. Edmonson, Liqing Tian, Michael Rusch, David A. Wheeler, Jennifer L. Neary, Scott Newman, Lu Wang, Patrick R. Blackburn, Michael Macias, Andrew Thrasher, Jian Wang, Mark R. Wilkinson, Xin Zhou, Jinghui Zhang. Fuzzion2: Fast, sensitive detection of known gene fusions by fuzzy pattern matching for clinical testing and large-scale data mining [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 4092.

Scientific visualization is important to cancer research as different graphs enable researchers to spot trends and detect outliers in large-scale, high-dimensional data. For example, mutation landscape maps integrate mutational and clinical data, revealing which somatic alterations are most closely associated with diagnoses, age, ethnicity, and outcome; t-SNE plots allow dissection of cancer subtypes based on gene expression or methylation profiling. These graphs are often created by researchers after extensive data curation and optimization of the layout using software tools designed for data visualization and are commonly presented as static figures in published literature. To permit seamless access to the genomic and epigenomic visualization tools created at St. Jude, we developed Visualization Community (https://viz.stjude.cloud/, VisComm), which allows researchers to explore published graphs, develop custom visualization for their private data, and also serve as a public repository for completed graphs, as part of the St. Jude Cloud data sharing ecosystem. Sixty-five (65) interactive graphs from 31 publications in 231 pediatric subtypes are currently available, including t-SNE plots of RNA-seq and methylation data, mutation landscape maps, variant visualization on protein or genomic axes, epigenome states, and chromatin contact interactions. The graphs are searchable by cancer diagnosis, molecular profiling, research interest, publication, and research organization with features designed to enable user customization. To date, 5,014 users, 25% of whom are recurrent, have accessed VisComm. On average, each public visualization has over 4,000 views while the most accessed graph has >13,000 views. We have developed the ability to work in teams with in VisComm to empower collaborative development of new visualizations and controlled release by a research team. This feature enables the professional quality visualizations serve as resources for focused scientific communities such as the Audacious Goals Initiative (AGI) - Retina), a multi-institutional cross-disciplinary research team whose goal is restoring vision to patients by regenerating retinas. Videos and textual tutorials are being prepared using data generated from programs such as the Genome for Kids (G4K), a clinical genomic profiling program, that will help guide users to apply the underlying visualization software tools and to integrate public and private data of all types. VisComm is the first platform dedicated to enabling the development and sharing of interactive high-quality data-rich visualizations, which we expect will enhance data exploration and hypothesis generation for cancer research. Citation Format: David B. Finkelstein, James Madson, Lucian Vacaroiu, Alexander M. Gout, Leonard Hirja, Stephanie Sandor, Andrew Thrasher, Colleen Reilly, Jian Wang, Xiaolong Chen, Xin Zhou, Jinghui Zhang, Clay McLeod. VisComm: A cloud-based platform for accessing and creating scientific visualization for the pediatric cancer research community [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1901.

Clinical genomic studies of pediatric cancer have primarily focused on specific tumor types or high-risk disease. In the Genomes for Kids study (NCT02530658) we used a three-platform sequencing approach, including whole genome (WGS), whole exome (WES) and RNA sequencing, to examine tumor and paired germline genomes from prospectively identified children with cancer. The goal of the study was to assess the potential of comprehensive next generation sequencing to elucidate the molecular mechanisms underlying tumor formation and investigate the potential of this information to influence clinical decision-making.The cohort, with a median age of 6 yrs, range 0 - 26 yrs, included 301 patients with newly diagnosed (85%) or relapsed/refractory (15%) cancers, unselected for tumor type or stage. Patients with hematologic malignancies accounted for 41% of cases, 31% had CNS tumors, and 28% had other non-CNS solid tumors. This cohort also included 18 patients with very rare tumor types, defined here as occurring in less than 2 cases per million person per year.Two hundred fifty three patients (84%) had sufficient tumor for three-platform sequencing and all 301 had adequate paired germline samples. Following analysis, 86% of patients harbored diagnostic (53%), prognostic (57%), therapeutically relevant (25%), and/or cancer predisposing (18%) variants. The inclusion of WGS enabled detection of oncogenic gene fusions, as well as 22 cases in which oncogenes were activated through enhancer hijacking, a particularly frequent occurrence in hematologic malignancies. In addition, WGS effectively detected clinically relevant small intragenic deletions (15% of tumors) and a variety of mutational signatures, which were not detectable through analysis of whole exome data. Evaluation of 56 pathogenic germline variants in the context of paired tumor sequence data helped establish the disease relevance of several genes that are not typically associated with the cancer type in question, providing critical insights on a case-by-case basis. Examples include a pathogenic germline variant in MUTYH in a patient with retinoblastoma whose tumor exhibited a mutation signature associated with reactive oxygen species indicative of loss of MUTYH function; and conversely, a likely pathogenic variant in PMS2 in a rare brain cancer, which did not exhibit a mutation signature associated with microsatellite instability. This study successfully demonstrated the power of this three-platform approach to interrogate and interpret the full range of genomic variants across newly diagnosed as well as relapsed/refractory pediatric cancers. As a result of these findings, we have incorporated this three-platform approach into our routine real-time clinical service at St. Jude Children's Hospital. Citation Format: David A. Wheeler, Scott Newman, Joy Nakitandwe, Chimene A. Kesserwan, Elizabeth M. Azzato, Michael C. Rusch, Sheila Shurtleff, Armita Bahrami, Brent Orr, Jeffery M. Klco, Dale J. Hedges, Kayla V. Hamilton, Scott G. Foy, Michael N. Edmonson, Andrew Thrasher, Jiali Gu, Lynn W. Harrison, Lu Wang, Roya Mostafavi, Manish Kubal, Jamie Maciaszek, Michael Clay, Annastasia Ouma, Antonina Silkov, Yanling Liu, Zhaojie Zhang, Yu Liu, Samuel W. Brady, Xin Zhou, Mark Wilkinson, Delaram Rahbarinia, Jay Knight, Jian Wang, Charles G. Mullighan, Rose B. McGee, Emily A. Quinn, Elsie L. Gerhardt, Leslie M. Taylor, Regina Nuccio, Jessica M. Valdez, Stacy J. Hines-Dowell, Alberto Pappo, Giles Robinson, Liza-Marie Johnson, Ching-Hon Pui, David W. Ellison, James R. Downing, Jinghui Zhang, Kim E. Nichols. Genomes for Kids: Comprehensive DNA and RNA sequencing defining the scope of actionable mutations in pediatric cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 642.

Next-generation sequencing-based genomic profiling is now a mainstay of pediatric oncology research and clinical testing. Correlating genomic features of patient cancer genomes with curated data extracted from large reference cohorts is critical for identifying molecular subtypes and underlying mutagenesis processes. To facilitate such investigation, we developed two user-friendly workflows on St. Jude Cloud, a data sharing ecosystem hosting genomic data for >10,000 pediatric cancer patients and survivors. These workflows leverage St. Jude Cloud comprehensive pediatric cancer genomic data, including 1,616 RNA-seq of 135 cancer subtypes and 958 whole genome sequencing (WGS) of 35 subtypes, to enable user analysis of their data in the context of St. Jude Cloud cohorts without a need to download large datasets. The RNA-Seq Expression Classification workflow enables a user to compare their patient RNA-Seq gene expression data with blood (832), brain (456), and solid tumor (319) pediatric cancer reference cohorts and PDX models (45), enabling subtype classification using t-Distributed Stochastic Neighbor Embedding (t-SNE). Reference cohorts include curated subtype-defining somatic alterations integrating genomic variant data with expression profile. Resulting interactive t-SNE plots can be explored and annotated - with options to highlight cancer subtypes or samples and display sample information (age of onset, clinical diagnosis, molecular driver). To demonstrate, we analyze PAWNXH, a Children's Oncology Group AML sample with a novel ZBTB7A-NUTM1 fusion and find it clusters with AML samples harboring KMT2A re-arrangements suggesting a potential mechanism of pathogenesis. Integrating PDX samples enables model selection for functional experiments by connecting patient subtypes with mouse models. The Mutational Signatures workflow identifies and quantifies COSMIC mutational signatures in user-uploaded somatic VCF files for comparison to reference pediatric cancer cohorts. The interactive interface enables rapid identification of signatures within the query cohort and facilitates comparison to the reference using a cohort-level summary view. Identified signatures may also be explored at the sample-level for both query and reference cohorts, enabling the user to identify samples with signatures of interest for further analysis. We show an example comparison of mutational signatures identified in pediatric and adult AML samples. These workflows enable users to leverage curated pediatric cancer data to make discoveries in their own samples. Enabling point-and-click analysis in St. Jude Cloud removes the barrier for non-computational researchers and eliminates the need to download large reference datasets for local analysis. Both workflows utilize post-processed rather than raw genomic data, reducing transfer costs for uploading user data to the cloud. Citation Format: Andrew Thrasher, Michael Macias, Alexander M. Gout, Delaram Rahbarinia, Xin Zhou, Samuel W. Brady, Clay McLeod, Michael C. Rusch, Xiaolong Chen, Soheil Meshinchi, Michael A. Dyer, Suzanne J. Baker, Martine F. Roussel, Jinghui Zhang. Empowering point-and-click genomic analysis with large pediatric genomic reference data on St. Jude Cloud [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2289.

Gene fusions are important biomarkers for cancer diagnosis, subtype classification and therapeutic decision-making. While fusion detection using RNA-seq data has become a standard practice, existing computational methods primarily focus on identifying canonical exon-to-exon fusions. However, more complex events such as multi-partner fusions, truncations, enhancer hijacking and internal tandem duplications (ITD) can also lead to abnormal function or aberrant transcription of cancer driver genes. To aid discovery of complex and diverse driver fusions, we developed CICERO (CICERO Is Clipping Extended for RNA Optimization), a local assembly-based algorithm that integrates RNA-seq reads bearing aberrant mapping signatures with extensive annotation for ranking candidate fusions. Our benchmark data set, designed to support the main application of RNA-seq fusion analysis, consists of 184 driver fusions from 170 pediatric leukemia, solid tumor and brain tumor detected by paired tumor-normal WGS and orthogonally validated by capture sequencing, RT-PCR and/or FISH. CICERO detected 95% of these fusions with an average ranking of 1.9, whereas ChimeraScan, deFuse, FusionCatcher and STAR-Fusion detected only 63%, 66%, 77% and 63% with an average ranking of 37.0, 9.0, 18.1 and 4.4, respectively. Notably, events such as ITD and rearrangements involving the highly repetitive IGH locus were detected almost exclusively by CICERO. Our re-analysis of 167 RNA-seq data from the TCGA Glioblastoma Multiforme (GBM) cohort unveiled 158 fusions of cancer genes that were not reported previously. These include kinase fusions (KLHL7-BRAF), ITD of EGFR kinase domain and a 13% prevalence of EGFR C-terminal truncation compared to the 6% reported by the TCGA Network. CICERO has greatly improved our ability to discover non-canonical fusions which are overlooked by existing fusion detection methods, and has been used to analyze >2,000 RNA-seq samples generated by the two largest pediatric cancer genomics initiatives: the St. Jude/Washington University Pediatric Cancer Genome Project (PCGP) and the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) project. We anticipate that CICERO will also improve fusion analysis for adult cancer RNA-seq data, as demonstrated through our re-analysis of TCGA-GBM and our recent discovery of MAP3K8 C-terminal truncation fusion in 2% of TCGA melanoma samples. CICERO is accessible via standard (https://github.com/stjude/Cicero) or cloud-based (https://platform.stjude.cloud/tools/rapid_rna-seq) implementation. To further improve accuracy, fusions predicted by CICERO can be curated by FusionEditor (https://proteinpaint.stjude.org/FusionEditor/), an interactive viewer allowing inspection of protein domains involved in the fusion and the gene expression status of fusion-positive samples. Citation Format: Liqing Tian, Yongjin Li, Michael N. Edmonson, Xin Zhou, Scott Newman, Clay McLeod, Yu Liu, Bo Tang, Michael C. Rusch, John Easton, Jing Ma, Austyn Trull, J. Robert Michael, Andrew Thrasher, Charles Mullighan, Suzanne J. Baker, James R. Downing, David W. Ellison, Jinghui Zhang. CICERO: An accurate method for detecting complex and diverse driver fusions using cancer transcriptome sequencing (RNA-seq) data [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5478.

Childhood cancer survivors are at increased risk of subsequent neoplasms (SN), largely considered to be therapy-related. Studies of cancer predisposition genes (CPGs) and risk of SN among long-term survivors are lacking. We characterized germline mutations in CPGs in childhood cancer survivors to determine their contribution to SN risk. Whole genome (30x) and exome (100x) sequencing was performed for 2988 5+ year survivors of childhood cancer (1629 leukemia/lymphoma, 332 CNS, 1027 other solid tumors, 53% male, median follow-up 28 [range 6-55] years). Survivors underwent a comprehensive clinical assessment, treatment exposures were abstracted from medical records, and SN were validated by pathology reports. Germline mutations in 63 CPGs were classified using the American College of Medical Genetics and Genomics guidelines as previously described (Zhang et al. NEJM 2015). Logistic regression, adjusting for age, sex and race, was used to evaluate associations between mutation status, cancer therapy and the SN risk. 1062 SNs were diagnosed in 437 survivors, of whom 98 developed ≥2 histologically distinct SNs. Median age at SN and time to first SN was 38.2 (range 3.3-67.4) and 29.2 (0.9-48.4) years, respectively. Common SNs were basal cell carcinoma (542 in 153 survivors), meningioma (201 in 100), thyroid (64 in 64), and breast cancer (58 in 50). Cumulative incidence of SN at age 45 was 25.5% (95% CI: 22.9-27.9). 169 survivors (5.7%) had a pathogenic/likely pathogenic (P/LP) mutation in a CPG, consisting of 97 single nucleotide variations, 63 insertion/deletions and 9 copy number alterations (49% of mutations not in ClinVar). Frequently mutated genes were: RB1 (n=41), NF1 (n=22), BRCA2 (n=13), BRCA1 (n=12) and TP53 (n=10). Our data confirmed known associations between CPG mutations and specific primary diagnoses including RB1 mutations in 32 of 41 (78%) of bilateral and 7 of 57 (12%) of unilateral retinoblastoma survivors, 22 NF1 (20 of 332 CNS survivors), 4 SUFU (all in medulloblastoma survivors) and 5 WT1 mutations (all in Wilms’ tumor survivors). Analyses revealed novel associations between CPG mutations and SN risk. Among 1326 survivors not exposed to radiation therapy (non-RT), 62 SNs developed in 54 survivors, of which 15 (24.2%) occurred in P/LP mutation carriers. Non-RT exposed survivors with a P/LP mutation had an increased risk of SN (OR=5.6, 95% CI=2.6-12.0, P Citation Format: Zhaoming Wang, Carmen L. Wilson, John Easton, Dale Hedges, Qi Liu, Gang Wu, Michael Rusch, Michael Edmonson, Shawn Levy, Jennifer Q. Lanctot, Eric Caron, Kyla Shelton, Kelsey Currie, Matthew Lear, Heather L. Mulder, Donald Yergeau, Celeste Rosencrance, Bhavin Vadodaria, Yadav Sapkota, Russell J. Brooke, Wonjong Moon, Evadnie Rampersaud, Xiaotu Ma, Shuoguo Wang, Ti-Cheng Chang, Stephen Rice, Andrew Thrasher, Aman Patel, Cynthia Pepper, Xin Zhou, Xiang Chen, Wenan Chen, Angela Jones, Braden Boone, Deo Kumar Srivastava, Chimene A. Kesserwan, Kim E. Nichols, James R. Downing, Melissa M. Hudson, Yutaka Yasui, Leslie L. Robison, Jinghui Zhang. Germline mutations in cancer predisposition genes and risk for subsequent neoplasms among long-term survivors of childhood cancer in the St. Jude Lifetime Cohort [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 3001. doi:10.1158/1538-7445.AM2017-3001

Next-generation sequencing (NGS) of the whole genome, whole exome, and transcriptome has enabled characterization of genetic landscapes of multiple cancers. By analyzing over 2,000 pediatric cancer patients, we have developed a comprehensive database for recurrent somatic alterations and pathogenic germline mutations as part of the St. Jude/Washington University Pediatric Cancer Genome Project. However, there is no systematic evaluation on whether NGS is able to identify germline and somatic lesions reported by existing molecular diagnostic assays and what combination of NGS platforms is best suited for clinical sequencing. Here we report the first comprehensive study that employs whole-genome sequencing at 30-45X coverage, whole-exome sequencing at 100X coverage and transcriptome sequencing using matched tumor/normal samples from cancer patients. A pilot study was carried out to perform NGS analysis on 78 children of leukemia, solid tumor or brain tumor with a total of 112 diagnostic or prognostic biomarkers previously characterized by multiple molecular diagnostic assays. We implemented an analysis pipeline that integrates the genetic lesions detected by all three NGS platforms to characterize somatic and germline single nucleotide variations (SNVs), short insertions and deletions (indels), structural variations including fusions, karyotypes, copy number alterations, loss of heterozygosity, tumor purity and tumor-in-normal contamination. The turn-around time for data analysis is 2 weeks with an overall sensitivity of 99% on detecting known biomarkers. Extensive validation of >3,000 somatic sequence mutations or structural variations from 38 cases shows that the specificity for somatic SNV, indel and structural variation is at 98%, 95% and 84% across the genome. We demonstrate that in addition to providing cross-validation, multi-platform NGS is required for detecting all genetic lesions of pathological significance including complex re-arrangements such as chromothripsis. In addition to known pathogenic or likely pathogenic mutations, our analysis has also unveiled novel pathogenic mutations (e.g. a germline deletion in TP53 in one patient with medulloblastoma) and identified multiple variants of unknown significance that may be worth further exploration (e.g. an in-frame deletion of exons 3-9 of DNMT3A in one neuroblastoma). Our study demonstrates that NGS is able to detect a wide range of genetic lesions currently characterized by multiple molecular diagnostic assays, providing critical insight into the design of clinical sequencing for ongoing studies. Citation Format: Jinghui Zhang, Michael Rusch, Joy Nakitandwe, Zhaojie Zhang, Michael N. Edmonson, Matthew Parker, Xiaotu Ma, Jared Becksfort, Andrew Thrasher, Jiali Gu, Yongjin Li, Erin Hedlund, Aman Patel, John Easton, Donald Yergeau, Bhavin Vadodaria, Xiang Chen, Tanja A. Gruber, Rose McGee, David Ellison, Sheila Shurtleff, James R. Downing. Molecular diagnosis for pediatric cancer through integrative analysis of whole-genome, whole-exome and transcriptome sequencing. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2628.