Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE): a cloud-based platform for curating and classifying germline variants

Patient-derived rhabdomyosarcoma cells recapitulate the genetic and transcriptomic landscapes of primary tumors

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in childhood and adolescence. The availability of appropriate and well-characterized preclinical models for RMS is limited, posing a challenge for investigating the molecular mechanisms and evaluating new targeted compounds in preclinical settings. Here, we collected 51 RMS specimens (referred to as ZJUCH-RMS cohort) and established 9 patient-derived cells (PDCs) and validated the identity of these cells by the expression of RMS-specific markers. Whole-transcriptome analysis identified high-confidence mutations in ZJUCH-RMS cohort including RAS, TP53, ARID1A, MYOD1, and MYCN. Further studies showed that RMS PDCs retained the genetic alterations and the expression of RMS hallmark and dependency genes in matched primary tumors and acted as valuable tools to assess drug responses and pharmacogenomic interactions. Our study provides unique PDCs that are available for preclinical studies of RMS and further advances the feasibility of RMS PDCs as valuable tools for developing personalized treatments for patients.

Whole-Exome Sequencing Reveals Novel Candidate Driver Mutations and Potential Druggable Mutations in Patients with High-Risk Neuroblastoma

Neuroblastoma is the most prevalent solid tumor in early childhood, with a 5-year overall survival rate of 40-60% in high-risk cases. Therefore, the identification of novel biomarkers for the diagnosis, prognosis, and therapy of neuroblastoma is crucial for improving the clinical outcomes of these patients. In this study, we conducted the whole-exome sequencing of 48 freshly frozen tumor samples obtained from the Biobank. Somatic variants were identified and selected using a bioinformatics analysis pipeline. The mutational signatures were determined using the Mutalisk online tool. Cancer driver genes and druggable mutations were predicted using the Cancer Genome Interpreter. The most common mutational signature was single base substitution 5. MUC4, MUC16, and FLG were identified as the most frequently mutated genes. Using the Cancer Genome Interpreter, we identified five recurrent cancer driver mutations spanning MUC16, MUC4, ALK, and CTNND1, with the latter being novel and containing a missense mutation, R439C. We also identified 11 putative actionable mutations including NF1 Q1798*, Q2616*, and S636X, ALK F1174L and R1275Q, SETD2 P10L and Q1829E, BRCA1 R612S, NOTCH1 D1670V, ATR S1372L, and FGFR1 N577K. Our findings provide a comprehensive overview of the novel information relevant to the underlying molecular pathogenesis and therapeutic targets of neuroblastoma.

Mutational mechanisms in multiply relapsed pediatric acute lymphoblastic leukemia

Pediatric acute lymphoblastic leukemia (ALL) is marked by low mutational load at initial diagnosis, which increases at relapse. To determine which processes are active in (relapsed) ALL and how they behave during disease progression before and after therapy, we performed whole genome sequencing on 97 tumor samples of 29 multiply relapsed ALL patients. Mutational load increased upon relapse in 28 patients and upon every subsequent relapse in 22 patients. In addition to two clock-like mutational processes, we identified UV-like damage, APOBEC activity, reactive oxygen species, thiopurine-associated damage and an unknown therapy component as drivers of mutagenesis. Mutational processes often affected patients over longer time periods, but could also occur in isolated events, suggesting the requirement of additional triggers. Thiopurine exposure was the most prominent source of new mutations in relapse, affecting over half of the studied patients in first and/or later relapse and causing potential relapse-driving mutations in multiple patients. Our data demonstrate that multiple mutational processes frequently act in parallel as prominent secondary drivers with dynamic activity during ALL development and progression.

Novel molecular subtypes of intracranial germ cell tumors expand therapeutic opportunities

BACKGROUND

Intracranial germ cell tumors (IGCTs) are a rare group of malignancies that are clinically classified as germinomas and nongerminomatous germ cell tumors (NGGCTs). Previous studies have found that somatic mutations involving the mitogen-activated protein kinase/mTOR signaling pathway are common early events. However, a comprehensive genomic understanding of IGCTs is still lacking.

METHODS

We established a cohort including over 100 IGCTs and conducted genomic and transcriptomic sequencing.

RESULTS

We identified novel recurrent driver genomic aberrations, including USP28 truncation mutations and high-level copy number amplification of KRAS and CRKL caused by replication of extrachromosomal DNA. Three distinct subtypes associated with unique genomic and clinical profiles were identified with transcriptome analysis: Immune-hot, MYC/E2F, and SHH. Both immune-hot and MYC/E2F were predominantly identified in germinomas and shared similar mutations involving the RAS/MAPK signaling pathway. However, the immune-hot group showed an older disease onset age and a significant immune response. MYC/E2F was characterized by a younger disease onset age and increased genomic instability, with a higher proportion of tumors showing whole-genome doubling. Additionally, the SHH subtype was mostly identified in NGGCTs.

CONCLUSIONS

Novel genomic aberrations and molecular subtypes were identified in IGCTs. These findings provide molecular basis for the potential introduction of new treatment strategies in this setting.

A new genomic framework to categorize pediatric acute myeloid leukemia

Recent studies on pediatric acute myeloid leukemia (pAML) have revealed pediatric-specific driver alterations, many of which are underrepresented in the current classification schemas. To comprehensively define the genomic landscape of pAML, we systematically categorized 887 pAML into 23 mutually distinct molecular categories, including new major entities such as UBTF or BCL11B, covering 91.4% of the cohort. These molecular categories were associated with unique expression profiles and mutational patterns. For instance, molecular categories characterized by specific HOXA or HOXB expression signatures showed distinct mutation patterns of RAS pathway genes, FLT3 or WT1, suggesting shared biological mechanisms. We show that molecular categories were strongly associated with clinical outcomes using two independent cohorts, leading to the establishment of a new prognostic framework for pAML based on these updated molecular categories and minimal residual disease. Together, this comprehensive diagnostic and prognostic framework forms the basis for future classification of pAML and treatment strategies.

Establishment and Comprehensive Characterization of a Novel Preclinical Platform of Metastatic Retinoblastoma for Therapeutic Developments

Purpose

Although there have been improvements in the management of metastatic retinoblastoma, most patients do not survive, and all patients suffer from multiple short- and long-term treatment toxicities. Reliable and informative models to assist clinicians are needed. Thus we developed and comprehensively characterized a novel preclinical platform of primary cell cultures and xenograft models of metastatic retinoblastoma to provide insights into the molecular biology underlying metastases and to perform drug screening for the identification of hit candidates with the highest potential for clinical translation.

Methods

Orbital tumor, bone marrow, cerebrospinal fluid, and lymph node tumor infiltration specimens were obtained from seven patients with metastatic retinoblastoma at diagnosis, disease progression, or relapse. Tumor specimens were engrafted in immunodeficient animals, and primary cell lines were established. Genomic, immunohistochemical/immunocytochemical, and pharmacological analysis were performed.

Results

We successfully established five primary cell lines: two derived from leptomeningeal, two from orbital, and one from lymph node tumor dissemination. After the intravitreal or intraventricular inoculation of these cells, we established cell-derived xenograft models. Both primary cell lines and xenografts accurately retained the histological and genomic features of the tumors from which they were derived and faithfully recapitulated the dissemination patterns and pharmacological sensitivity observed in the matched patients.

Conclusions

Ours is an innovative and thoroughly characterized preclinical platform of metastatic retinoblastoma developed for the understanding of tumor biology of this highly aggressive tumor and has the potential to identify drug candidates to treat patients who currently lack effective treatment options.

CR1 variants contribute to FSGS susceptibility across multiple populations

Focal segmental glomerulosclerosis (FSGS) is a common cause of nephrotic syndrome with an annual incidence in the United States in African-Americans compared to European-Americans of 24 cases and 5 cases per million, respectively. Among glomerular diseases in Europe and Latin-America, FSGS was the second most frequent diagnosis, and in Asia the fifth. We expand previous efforts in understanding genetics of FSGS by performing a case-control study involving ethnically-diverse groups FSGS cases (726) and a pool of controls (13,994), using panel sequencing of approximately 2,500 podocyte-expressed genes. Through rare variant association tests, we replicated known risk genes - KANK1, COL4A4, and APOL1. A novel significant association was observed for the gene encoding complement receptor 1 (CR1). High-risk rare variants in CR1 in the European-American cohort were commonly observed in Latin- and African-Americans. Therefore, a combined rare and common variant analysis was used to replicate the CR1 association in non-European populations. The CR1 risk variant, rs17047661, gives rise to the Sl1/Sl2 (R1601G) allele that was previously associated with protection against cerebral malaria. Pleiotropic effects of rs17047661 may explain the difference in allele frequencies across continental ancestries and suggest a possible role for genetically-driven alterations of adaptive immunity in the pathogenesis of FSGS.

Cancer germline predisposing variants and late mortality from subsequent malignant neoplasms among long-term childhood cancer survivors: a report from the St Jude Lifetime Cohort and the Childhood Cancer Survivor Study

BACKGROUND

Carriers of cancer predisposing variants are at an increased risk of developing subsequent malignant neoplasms among those who have survived childhood cancer. We aimed to investigate whether cancer predisposing variants contribute to the risk of subsequent malignant neoplasm-related late mortality (5 years or more after diagnosis).

METHODS

In this analysis, data were included from two retrospective cohort studies, St Jude Lifetime Cohort (SJLIFE) and the Childhood Cancer Survivor Study (CCSS), with prospective follow-up of patients who were alive for at least 5 years after diagnosis with childhood cancer (ie, long-term childhood cancer survivors) with corresponding germline whole genome or whole exome sequencing data. Cancer predisposing variants affecting 60 genes associated with well-established autosomal-dominant cancer-predisposition syndromes were characterised. Subsequent malignant neoplasms were graded using the National Cancer Institute Common Terminology Criteria for Adverse Events (CTCAE) version 4.03 with modifications. Cause-specific late mortality was based on linkage with the US National Death Index and systematic cohort follow up. Fine-Gray subdistribution hazard models were used to estimate subsequent malignant neoplasm-related late mortality starting from the first biospecimen collection, treating non-subsequent malignant neoplasm-related deaths as a competing risk, adjusting for genetic ancestry, sex, age at diagnosis, and cancer treatment exposures. SJLIFE (NCT00760656) and CCSS (NCT01120353) are registered with ClinicalTrials.gov.

FINDINGS

12 469 (6172 male and 6297 female) participants were included, 4402 from the SJLIFE cohort (median follow-up time since collection of the first biospecimen 7·4 years [IQR 3·1-9·4]) and 8067 from the CCSS cohort (median follow-up time since collection of the first biospecimen 12·6 years [2·2-16·6]). 641 (5·1%) of 12 469 participants carried cancer predisposing variants (294 [6·7%] in the SJLIFE cohort and 347 [4·3%] in the CCSS cohort), which were significantly associated with an increased severity of subsequent malignant neoplasms (CTCAE grade ≥4 vs grade <4: odds ratio 2·15, 95% CI 1·18-4·19, p=0·0085). 263 (2·1%) subsequent malignant neoplasm-related deaths (44 [1·0%] in the SJLIFE cohort; and 219 [2·7%] in the CCSS cohort) and 426 (3·4%) other-cause deaths (103 [2·3%] in SJLIFE; and 323 [4·0%] in CCSS) occurred. Cumulative subsequent malignant neoplasm-related mortality at 10 years after the first biospecimen collection in carriers of cancer predisposing variants was 3·7% (95% CI 1·2-8·5) in SJLIFE and 6·9% (4·1-10·7) in CCSS versus 1·5% (1·0-2·1) in SJLIFE and 2·1% (1·7-2·5) in CCSS in non-carriers. Carrying a cancer predisposing variant was associated with an increased risk of subsequent malignant neoplasm-related mortality (SJLIFE: subdistribution hazard ratio 3·40 [95% CI 1·37-8·43]; p=0·0082; CCSS: 3·58 [2·27-5·63]; p<0·0001).

INTERPRETATION

Identifying participants at increased risk of subsequent malignant neoplasms via genetic counselling and clinical genetic testing for cancer predisposing variants and implementing early personalised cancer surveillance and prevention strategies might reduce the substantial subsequent malignant neoplasm-related mortality burden.

FUNDING

American Lebanese Syrian Associated Charities and US National Institutes of Health.

Case report: Somatic mutations in microtubule dynamics-associated genes in patients with WNT-medulloblastoma tumors

Medulloblastoma (MB) is the most common pediatric brain tumor which accounts for about 20% of all pediatric brain tumors and 63% of intracranial embryonal tumors. MB is considered to arise from precursor cell populations present during an early brain development. Most cases (~70%) of MB occur at the age of 1-4 and 5-9, but are also infrequently found in adults. Total annual frequency of pediatric tumors is about 5 cases per 1 million children. WNT-subtype of MB is characterized by a high probability of remission, with a long-term survival rate of about 90%. However, in some rare cases there may be increased metastatic activity, which dramatically reduces the likelihood of a favorable outcome. Here we report two cases of MB with a histological pattern consistent with desmoplastic/nodular (DP) and classic MB, and genetically classified as WNT-MB. Both cases showed putative causal somatic protein truncating mutations identified in microtubule-associated genes: ARID2, TUBB4A, and ANK3.

The genomic landscape of pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Here, using whole-genome, exome and transcriptome sequencing of 2,754 childhood patients with ALL, we find that, despite a generally low mutation burden, ALL cases harbor a median of four putative somatic driver alterations per sample, with 376 putative driver genes identified varying in prevalence across ALL subtypes. Most samples harbor at least one rare gene alteration, including 70 putative cancer driver genes associated with ubiquitination, SUMOylation, noncoding transcripts and other functions. In hyperdiploid B-ALL, chromosomal gains are acquired early and synchronously before ultraviolet-induced mutation. By contrast, ultraviolet-induced mutations precede chromosomal gains in B-ALL cases with intrachromosomal amplification of chromosome 21. We also demonstrate the prognostic significance of genetic alterations within subtypes. Intriguingly, DUX4- and KMT2A-rearranged subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of the ALL genomic landscape and associated outcomes.

Hiplot: a comprehensive and easy-to-use web service for boosting publication-ready biomedical data visualization

Complex biomedical data generated during clinical, omics and mechanism-based experiments have increasingly been exploited through cloud- and visualization-based data mining techniques. However, the scientific community still lacks an easy-to-use web service for the comprehensive visualization of biomedical data, particularly high-quality and publication-ready graphics that allow easy scaling and updatability according to user demands. Therefore, we propose a community-driven modern web service, Hiplot (https://hiplot.org), with concise and top-quality data visualization applications for the life sciences and biomedical fields. This web service permits users to conveniently and interactively complete a few specialized visualization tasks that previously could only be conducted by senior bioinformatics or biostatistics researchers. It covers most of the daily demands of biomedical researchers with its equipped 240+ biomedical data visualization functions, involving basic statistics, multi-omics, regression, clustering, dimensional reduction, meta-analysis, survival analysis, risk modelling, etc. Moreover, to improve the efficiency in use and development of plugins, we introduced some core advantages on the client-/server-side of the website, such as spreadsheet-based data importing, cross-platform command-line controller (Hctl), multi-user plumber workers, JavaScript Object Notation-based plugin system, easy data/parameters, results and errors reproduction and real-time updates mode. Meanwhile, using demo/real data sets and benchmark tests, we explored statistical parameters, cancer genomic landscapes, disease risk factors and the performance of website based on selected native plugins. The statistics of visits and user numbers could further reflect the potential impact of this web service on relevant fields. Thus, researchers devoted to life and data sciences would benefit from this emerging and free web service.

Transmembrane Transporter Sema3D Serves as a Tumor Suppressor in Localized Clear Cell Renal Cell Carcinoma

Background

The transmembrane transporter Sema3D is a vital molecule involved in axon guidance and carcinogenesis of variant malignancies. However, the relationship between Sema3D and clear cell renal cell carcinoma (ccRCC) is barely reported and remains unclear.

Methods

Sema3D expression and the connection of clinical and histological characteristics were first analyzed with transcriptome data in the TCGA repository. We then located and examined the Sema3D expression in ccRCC patients by using immunofluorescence staining in the tissue microarray. The prognostic value of Sema3D in localized ccRCC was evaluated by Cox proportional hazard analysis. Functional and gene set enrichment analysis (GSEA), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed to describe the potential mechanisms of Sema3D in ccRCC. Correlation analysis between Sema3D and tumor-infiltrating lymphocytes was calculated by ssGSEA.

Results

In 86 ccRCC patients, Sema3D mRNA and protein expression were downregulated in tumor tissues than the para-tumor tissues, and Sema3D was dominantly expressed in the extracellular space. Low expression of Sema3D was associated with advanced tumor stage, advanced histological grade, and poor prognosis in ccRCC. In the subgroup analysis of 81 localized ccRCC patients, Sema3D expression level was an independent protective prognostic factor for overall survival (OS) (HR = 0.125, p=0.043). Coagulation, complement, estrogen response, and KRAS signaling hallmark gene sets were identified as Sema3D-related signaling pathways. The expression level of Sema3D was significantly correlated with a high abundance of several immune cells (neutrophils, eosinophils, and T helper cells).

Conclusions

Transmembrane transporter Sema3D is an efficient prognostic biomarker for localized ccRCC patients, by playing the role of tumor suppressor in ccRCC. Sema3D can be a novel therapeutic target for ccRCC.

Integrated Genomic Analysis Identifies UBTF Tandem Duplications as a Recurrent Lesion in Pediatric Acute Myeloid Leukemia

The genetics of relapsed pediatric acute myeloid leukemia (AML) has yet to be comprehensively defined. Here, we present the spectrum of genomic alterations in 136 relapsed pediatric AMLs. We identified recurrent exon 13 tandem duplications (TD) in upstream binding transcription factor (UBTF) in 9% of relapsed AML cases. UBTF-TD AMLs commonly have normal karyotype or trisomy 8 with cooccurring WT1 mutations or FLT3-ITD but not other known oncogenic fusions. These UBTF-TD events are stable during disease progression and are present in the founding clone. In addition, we observed that UBTF-TD AMLs account for approximately 4% of all de novo pediatric AMLs, are less common in adults, and are associated with poor outcomes and MRD positivity. Expression of UBTF-TD in primary hematopoietic cells is sufficient to enhance serial clonogenic activity and to drive a similar transcriptional program to UBTF-TD AMLs. Collectively, these clinical, genomic, and functional data establish UBTF-TD as a new recurrent mutation in AML.

SIGNIFICANCE

We defined the spectrum of mutations in relapsed pediatric AML and identified UBTF-TDs as a new recurrent genetic alteration. These duplications are more common in children and define a group of AMLs with intermediate-risk cytogenetic abnormalities, FLT3-ITD and WT1 alterations, and are associated with poor outcomes. See related commentary by Hasserjian and Nardi, p. 173. This article is highlighted in the In This Issue feature, p. 171.

Distinct genomic landscape of Chinese pediatric acute myeloid leukemia impacts clinical risk classification

Studies have revealed key genomic aberrations in pediatric acute myeloid leukemia (AML) based on Western populations. It is unknown to what extent the current genomic findings represent populations with different ethnic backgrounds. Here we present the genomic landscape of driver alterations of Chinese pediatric AML and discover previously undescribed genomic aberrations, including the XPO1-TNRC18 fusion. Comprehensively comparing between the Chinese and Western AML cohorts reveal a substantially distinct genomic alteration profile. For example, Chinese AML patients more commonly exhibit mutations in KIT and CSF3R, and less frequently mutated of genes in the RAS signaling pathway. These differences in mutation frequencies lead to the detection of previously uncharacterized co-occurring mutation pairs. Importantly, the distinct driver profile is clinical relevant. We propose a refined prognosis risk classification model which better reflected the adverse event risk for Chinese AML patients. These results emphasize the importance of genetic background in precision medicine.

The genomic landscape of pediatric acute myeloid leukemia (AML) has mostly been characterised for Western populations. Here, the authors identify potential driver alterations in Chinese pediatric AML, which differ from Western populations, and propose a prognostic risk classification model.

Germline Variants in Cancer Genes from Young Breast Cancer Mexican Patients

Breast cancer (BC) is one of the most frequent cancer types in women worldwide. About 7% is diagnosed in young women (YBC) less than 40 years old. In Mexico, however, YBC reaches 15% suggesting a higher genetic susceptibility. There have been some reports of germline variants in YBC across the world. However, there is only one report from a Mexican population, which is not restricted by age and limited to a panel of 143 genes resulting in 15% of patients carrying putatively pathogenic variants. Nevertheless, expanding the analysis to whole exome involves using more complex tools to determine which genes and variants could be pathogenic. We used germline whole exome sequencing combined with the PeCanPie tool to analyze exome variants in 115 YBC patients. Our results showed that we were able to identify 49 high likely pathogenic variants involving 40 genes on 34% of patients. We noted many genes already reported in BC and YBC worldwide, such as BRCA1, BRCA2, ATM, CHEK2, PALB2, and POLQ, but also others not commonly reported in YBC in Latin America, such as CLTCL1, DDX3X, ERCC6, FANCE, and NFKBIE. We show further supporting and controversial evidence for some of these genes. We conclude that exome sequencing combined with robust annotation tools and further analysis, can identify more genes and more patients affected by germline mutations in cancer.

Immunotherapy for SMARCB1-Deficient Sarcomas: Current Evidence and Future Developments

Mutations in subunits of the SWItch Sucrose Non-Fermentable (SWI/SNF) complex occur in 20% of all human tumors. Among these, the core subunit SMARCB1 is the most frequently mutated, and SMARCB1 loss represents a founder driver event in several malignancies, such as malignant rhabdoid tumors (MRT), epithelioid sarcoma, poorly differentiated chordoma, and renal medullary carcinoma (RMC). Intriguingly, SMARCB1-deficient pediatric MRT and RMC have recently been reported to be immunogenic, despite their very simple genome and low tumor mutational burden. Responses to immune checkpoint inhibitors have further been reported in some SMARCB1-deficient diseases. Here, we will review the preclinical data and clinical data that suggest that immunotherapy, including immune checkpoint inhibitors, may represent a promising therapeutic strategy for SMARCB1-defective tumors. We notably discuss the heterogeneity that exists among the spectrum of malignancies driven by SMARCB1-loss, and highlight challenges that are at stake for developing a personalized immunotherapy for these tumors, notably using molecular profiling of the tumor and of its microenvironment.

Integrated genomic analyses identify high-risk factors and actionable targets in T-cell acute lymphoblastic leukemia

T cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematologic malignancy often associated with poor outcomes. To identify high-risk factors and potential actionable targets for T-ALL, we perform integrated genomic and transcriptomic analyses on samples from 165 Chinese pediatric and adult T-ALL patients, of whom 85% have outcome information. The genomic mutation landscape of this Chinese cohort is very similar to the Western cohort published previously, except that the rate of NOTCH1 mutations is significant lower in the Chinese T-ALL patients. Among 47 recurrently mutated genes in 7 functional categories, we identify RAS pathway and PTEN mutations as poor survival factors for non-TAL and TAL subtypes, respectively. Mutations in the PI3K pathway are mutually exclusive with mutations in the RAS and NOTCH1 pathways as well as transcription factors. Further analysis demonstrates that approximately 43% of the high-risk patients harbor at least one potential actionable alteration identified in this study, and T-ALLs with RAS pathway mutations are hypersensitive to MEKi in vitro and in vivo. Thus, our integrated genomic analyses not only systematically identify high-risk factors but suggest that these high-risk factors are promising targets for T-ALL therapies.

Diffuse leptomeningeal glioneuronal tumour (DLGNT) in children: the emerging role of genomic analysis

Diffuse leptomeningeal glioneuronal tumours (DLGNT) represent rare enigmatic CNS tumours of childhood. Most patients with this disease share common radiological and histopathological features but the clinical course of this disease is variable. A radiological hallmark of this disease is widespread leptomeningeal enhancement that may involve the entire neuroaxis with predilection for the posterior fossa and spine. The classic pathologic features include low- to moderate-density cellular lesions with OLIG2 expression and evidence of ‘oligodendroglioma-like’ appearance. The MAPK/ERK signaling pathway has recently been reported as a potential driver of tumourigenesis in up to 80% of DLGNT with KIAA1549:BRAF fusions being the most common event seen. Until now, limited analysis of the biological drivers of tumourigenesis has been undertaken via targeted profiling, chromosomal analysis and immunohistochemistry. Our study represents the first examples of comprehensive genomic sequencing in DLGNT and shows that it is not only feasible but crucial to our understanding of this rare disease. Moreover, we demonstrate that DLGNT may be more genomically complex than single-event MAPK/ERK signaling pathway tumours.

Genomes for Kids: The scope of pathogenic mutations in pediatric cancer revealed by comprehensive DNA and RNA sequencing

Genomic studies of pediatric cancer have primarily focused on specific tumor types or high-risk disease. Here, we used a three-platform sequencing approach, including whole genome (WGS), exome, and RNA sequencing, to examine tumor and germline genomes from 309 prospectively identified children with newly diagnosed (85%) or relapsed/refractory (15%) cancers, unselected for tumor type. Eighty-six percent of patients harbored diagnostic (53%), prognostic (57%), therapeutically-relevant (25%), and/or cancer predisposing (18%) variants. Inclusion of WGS enabled detection of activating gene fusions and enhancer hijacks (36% and 8% of tumors, respectively), small intragenic deletions (15% of tumors) and mutational signatures revealing of pathogenic variant effects. Evaluation of paired tumor-normal data revealed relevance to tumor development for 55% of pathogenic germline variants. This study demonstrates the power of a three-platform approach that incorporates WGS to interrogate and interpret the full range of genomic variants across newly diagnosed as well as relapsed/refractory pediatric cancers.

Data-driven approaches to advance research and clinical care for pediatric cancer

Pediatric cancer is a rare disease with a distinct etiology and mutational landscape compared with adult cancer. Multi-omics profiling of retrospective and prospective cohorts coupled with innovative computational analysis have been instrumental in uncovering mechanisms of tumorigenesis and drug resistance that are now informing pediatric cancer clinical therapy. In this review we present the major data resources of pediatric cancer and actionable insights into pediatric cancer etiology stemming from the identification of oncogenic gene fusions, mutational signature analysis, systems biology, cancer predisposition and survivorship studies - that have led to improved clinical diagnosis, discovery of new drug-targets, pharmacological therapy, and screening for genetic predisposition. Ultimately, integration of large-scale omics datasets generated through international collaboration is required to maximize the power of data-driven approaches to advance pediatric cancer research informing clinical therapy.

St. Jude Cloud: A Pediatric Cancer Genomic Data-Sharing Ecosystem

Effective data sharing is key to accelerating research to improve diagnostic precision, treatment efficacy, and long-term survival in pediatric cancer and other childhood catastrophic diseases. We present St. Jude Cloud (https://www.stjude.cloud), a cloud-based data-sharing ecosystem for accessing, analyzing, and visualizing genomic data from >10,000 pediatric patients with cancer and long-term survivors, and >800 pediatric sickle cell patients. Harmonized genomic data totaling 1.25 petabytes are freely available, including 12,104 whole genomes, 7,697 whole exomes, and 2,202 transcriptomes. The resource is expanding rapidly, with regular data uploads from St. Jude's prospective clinical genomics programs. Three interconnected apps within the ecosystem-Genomics Platform, Pediatric Cancer Knowledgebase, and Visualization Community-enable simultaneously performing advanced data analysis in the cloud and enhancing the Pediatric Cancer knowledgebase. We demonstrate the value of the ecosystem through use cases that classify 135 pediatric cancer subtypes by gene expression profiling and map mutational signatures across 35 pediatric cancer subtypes. SIGNIFICANCE: To advance research and treatment of pediatric cancer, we developed St. Jude Cloud, a data-sharing ecosystem for accessing >1.2 petabytes of raw genomic data from >10,000 pediatric patients and survivors, innovative analysis workflows, integrative multiomics visualizations, and a knowledgebase of published data contributed by the global pediatric cancer community.This article is highlighted in the In This Issue feature, p. 995.

RNAIndel: discovering somatic coding indels from tumor RNA-Seq data

MOTIVATION

Reliable identification of expressed somatic insertions/deletions (indels) is an unmet need due to artifacts generated in PCR-based RNA-Seq library preparation and the lack of normal RNA-Seq data, presenting analytical challenges for discovery of somatic indels in tumor transcriptome.

RESULTS

We present RNAIndel, a tool for predicting somatic, germline and artifact indels from tumor RNA-Seq data. RNAIndel leverages features derived from indel sequence context and biological effect in a machine-learning framework. Except for tumor samples with microsatellite instability, RNAIndel robustly predicts 88-100% of somatic indels in five diverse test datasets of pediatric and adult cancers, even recovering subclonal (VAF range 0.01-0.15) driver indels missed by targeted deep-sequencing, outperforming the current best-practice for RNA-Seq variant calling which had 57% sensitivity but with 14 times more false positives.

AVAILABILITY AND IMPLEMENTATION

RNAIndel is freely available at https://github.com/stjude/RNAIndel.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pre-Existing Tumoral B Cell Infiltration and Impaired Genome Maintenance Correlate with Response to Chemoradiotherapy in Locally Advanced Rectal Cancer

Locally advanced rectal cancer (LARC) remains a medical challenge. Reliable biomarkers to predict which patients will significantly respond to neoadjuvant chemoradiotherapy (nCRT) have not been identified. We evaluated baseline genomic and transcriptomic features to detect differences that may help predict response to nCRT. Eligible LARC patients received nCRT (3D-LCRT 50.4 Gy plus capecitabine 825 mg/m2/bid), preceded by three cycles of CAPOX in high systemic-relapse risk tumors, and subsequent surgery. Frozen tumor biopsies at diagnosis were sequenced using a colorectal cancer panel. Transcriptomic data was used for pathway and cell deconvolution inferential algorithms, coupled with immunohistochemical validation. Clinical and molecular data were analyzed according to nCRT outcome. Pathways related to DNA repair and proliferation (p < 0.005), and co-occurrence of RAS and TP53 mutations (p = 0.001) were associated with poor response. Enrichment of expression signatures related to enhanced immune response, particularly B cells and interferon signaling (p < 0.005), was detected in good responders. Immunohistochemical analysis of CD20+ cells validated the association of good response with B cell infiltration (p = 0.047). Findings indicate that the presence of B cells is associated with successful tumor regression following nCRT in LARC. The prevalence of simultaneous RAS and TP53 mutations along with a proficient DNA repair system that may counteract chemoradio-induced DNA damage was associated with poor response.

Metabolomics in acute myeloid leukemia

Acute myeloid leukemia (AML) is a complex, heterogenous hematological malignancy caused by mutations in myeloid differentiation and proliferation. Response to therapy and long-term outcomes vary widely based on chromosomal and molecular aberrations. Many platforms have been used to characterize and stratify AML. Metabolomics, the global profiling of small molecules in a biological sample, has emerged in the last decade as an important tool for studying the metabolic dependency of cancer cells. Metabolic reprogramming is not only an important manifestation of AML but clinically relevant for diagnosis, risk stratification and targeted drug development. In this review, we discuss notable metabolic studies of the last decade and their application to novel therapies.

Frequency of pathogenic/likely pathogenic germline variants in cancer-related genes among children with acute leukemia in Saudi Arabia

BACKGROUND

The frequency of pathogenic/likely pathogenic (P/LP) germline mutations in cancer-related genes among children with cancer in highly consanguineous populations is not well studied.

METHODS

Whole-exome sequencing of germline DNA was performed in 60 children with acute leukemia. We used the St. Jude Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE) data portal for the classification of germline variants by the St. Jude Medal Ceremony pipeline.

RESULTS

Fifty-seven patients had acute lymphoblastic leukemia (ALL) and three patients had acute myeloid leukemia. Parental consanguinity was present in 27 (45%) patients. All patients were of Arab ancestry. Three patients (5%) had a history of cancer in their siblings. Five patients (8.3%) had P/LP germline mutations in cancer-related genes. Three patients with B-ALL had heterozygous pathogenic mutations in TP53, BRCA1, and BRCA2; one patient with B-ALL had homozygous pathogenic mutation in PMS2; and one patient with T-ALL had LP homozygous mutation in AK2 that was associated with reticular dysgenesis. Among patients who had history of parental consanguinity, three (11%) had P/LP germline mutations compared with two (8%) in the absence of parental consanguinity. Fourteen (23%) patients had gold medal variants in cancer-related genes, 13 were heterozygous, and one was homozygous. Silver medal variants were present in 35 (58%) patients; all were heterozygous except one homozygous.

CONCLUSIONS

Children with acute leukemia in Saudi Arabia had low frequency of P/LP mutations in cancer-related genes despite the high rate of consanguinity. Larger studies using whole-genome sequencing are needed to further explore the heritability of childhood leukemia.