Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours.

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.

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.

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

Variant interpretation in the era of massively parallel sequencing is challenging. Although many resources and guidelines are available to assist with this task, few integrated end-to-end tools exist. Here, we present the Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE), a web- and cloud-based platform for annotation, identification, and classification of variations in known or putative disease genes. Starting from a set of variants in variant call format (VCF), variants are annotated, ranked by putative pathogenicity, and presented for formal classification using a decision-support interface based on published guidelines from the American College of Medical Genetics and Genomics (ACMG). The system can accept files containing millions of variants and handle single-nucleotide variants (SNVs), simple insertions/deletions (indels), multiple-nucleotide variants (MNVs), and complex substitutions. PeCanPIE has been applied to classify variant pathogenicity in cancer predisposition genes in two large-scale investigations involving >4000 pediatric cancer patients and serves as a repository for the expert-reviewed results. PeCanPIE was originally developed for pediatric cancer but can be easily extended for use for nonpediatric cancers and noncancer genetic diseases. Although PeCanPIE's web-based interface was designed to be accessible to non-bioinformaticians, its back-end pipelines may also be run independently on the cloud, facilitating direct integration and broader adoption. PeCanPIE is publicly available and free for research use.

Inhibition of SF3B1 by molecules targeting the spliceosome results in massive aberrant exon skipping.

The recent identification of compounds that interact with the spliceosome (sudemycins, spliceostatin A, and meayamycin) indicates that these molecules modulate aberrant splicing via SF3B1 inhibition. Through whole transcriptome sequencing, we have demonstrated that treatment of Rh18 cells with sudemycin leads to exon skipping as the predominant aberrant splicing event. This was also observed following reanalysis of published RNA-seq data sets derived from HeLa cells after spliceostatin A exposure. These results are in contrast to previous reports that indicate that intron retention was the major consequence of SF3B1 inhibition. Analysis of the exon junctions up-regulated by these small molecules indicated that these sequences were absent in annotated human genes, suggesting that aberrant splicing events yielded novel RNA transcripts. Interestingly, the length of preferred downstream exons was significantly longer than the skipped exons, although there was no difference between the lengths of introns flanking skipped exons. The reading frame of the aberrantly skipped exons maintained a ratio of 2:1:1, close to that of the cassette exons (3:1:1) present in naturally occurring isoforms, suggesting negative selection by the nonsense-mediated decay (NMD) machinery for out-of-frame transcripts. Accordingly, genes involved in NMD and RNAs encoding proteins involved in the splicing process were enriched in both data sets. Our findings, therefore, further elucidate the mechanisms by which SF3B1 inhibition modulates pre-mRNA splicing.

H3.3 K27M depletion increases differentiation and extends latency of diffuse intrinsic pontine glioma growth in vivo.

Histone H3 K27M mutation is the defining molecular feature of the devastating pediatric brain tumor, diffuse intrinsic pontine glioma (DIPG). The prevalence of histone H3 K27M mutations indicates a critical role in DIPGs, but the contribution of the mutation to disease pathogenesis remains unclear. We show that knockdown of this mutation in DIPG xenografts restores K27M-dependent loss of H3K27me3 and delays tumor growth. Comparisons of matched DIPG xenografts with and without K27M knockdown allowed identification of mutation-specific effects on the transcriptome and epigenome. The resulting transcriptional changes recapitulate expression signatures from K27M primary DIPG tumors and are strongly enriched for genes associated with nervous system development. Integrated analysis of ChIP-seq and expression data showed that genes upregulated by the mutation are overrepresented in apparently bivalent promoters. Many of these targets are associated with more immature differentiation states. Expression profiles indicate K27M knockdown decreases proliferation and increases differentiation within lineages represented in DIPG. These data suggest that K27M-mediated loss of H3K27me3 directly regulates a subset of genes by releasing poised promoters, and contributes to tumor phenotype and growth by limiting differentiation. The delayed tumor growth associated with knockdown of H3 K27M provides evidence that this highly recurrent mutation is a relevant therapeutic target.

Correction to: H3.3 K27M depletion increases differentiation and extends latency of diffuse intrinsic pontine glioma growth in vivo.

The original article can be found online.

CREST maps somatic structural variation in cancer genomes with base-pair resolution.

We developed 'clipping reveals structure' (CREST), an algorithm that uses next-generation sequencing reads with partial alignments to a reference genome to directly map structural variations at the nucleotide level of resolution. Application of CREST to whole-genome sequencing data from five pediatric T-lineage acute lymphoblastic leukemias (T-ALLs) and a human melanoma cell line, COLO-829, identified 160 somatic structural variations. Experimental validation exceeded 80%, demonstrating that CREST had a high predictive accuracy.

Exploration of Coding and Non-coding Variants in Cancer Using GenomePaint.

GenomePaint (https://genomepaint.stjude.cloud/) is an interactive visualization platform for whole-genome, whole-exome, transcriptome, and epigenomic data of tumor samples. Its design captures the inter-relatedness between DNA variations and RNA expression, supporting in-depth exploration of both individual cancer genomes and full cohorts. Regulatory non-coding variants can be inspected and analyzed along with coding variants, and their functional impact further explored by examining 3D genome data from cancer cell lines. Further, GenomePaint correlates mutation and expression patterns with patient outcomes, and supports custom data upload. We used GenomePaint to unveil aberrant splicing that disrupts the RING domain of CREBBP, discover cis activation of the MYC oncogene by duplication of the NOTCH1-MYC enhancer in B-lineage acute lymphoblastic leukemia, and explore the inter- and intra-tumor heterogeneity at EGFR in adult glioblastomas. These examples demonstrate that deep multi-omics exploration of individual cancer genomes enabled by GenomePaint can lead to biological insights for follow-up validation.

Shortened Leukocyte Telomere Length Associates with an Increased Prevalence of Chronic Health Conditions among Survivors of Childhood Cancer: A Report from the St. Jude Lifetime Cohort.

PURPOSE: We aimed to analyze and compare leukocyte telomere length (LTL) and age-dependent LTL attrition between childhood cancer survivors and noncancer controls, and to evaluate the associations of LTL with treatment exposures, chronic health conditions (CHC), and health behaviors among survivors. EXPERIMENTAL DESIGN: We included 2,427 survivors and 293 noncancer controls of European ancestry, drawn from the participants in St. Jude Lifetime Cohort Study (SJLIFE), a retrospective hospital-based study with prospective follow-up (2007-2016). Common nonneoplastic CHCs (59 types) and subsequent malignant neoplasms (5 types) were clinically assessed. LTL was measured with whole-genome sequencing data. RESULTS: After adjusting for age at DNA sampling, gender, genetic risk score based on 9 SNPs known to be associated with telomere length, and eigenvectors, LTL among survivors was significantly shorter both overall [adjusted mean (AM) = 6.20 kb; SE = 0.03 kb] and across diagnoses than controls (AM = 6.69 kb; SE = 0.07 kb). Among survivors, specific treatment exposures associated with shorter LTL included chest or abdominal irradiation, glucocorticoid, and vincristine chemotherapies. Significant negative associations of LTL with 14 different CHCs, and a positive association with subsequent thyroid cancer occurring out of irradiation field were identified. Health behaviors were significantly associated with LTL among survivors aged 18 to 35 years (P trend = 0.03). CONCLUSIONS: LTL is significantly shorter among childhood cancer survivors than noncancer controls, and is associated with CHCs and health behaviors, suggesting LTL as an aging biomarker may be a potential mechanistic target for future intervention studies designed to prevent or delay onset of CHCs in childhood cancer survivors.See related commentary by Walsh, p. 2281.

Molecular Mechanism of Telomere Length Dynamics and Its Prognostic Value in Pediatric Cancers.

BACKGROUND: We aimed to systematically evaluate telomere dynamics across a spectrum of pediatric cancers, search for underlying molecular mechanisms, and assess potential prognostic value. METHODS: The fraction of telomeric reads was determined from whole-genome sequencing data for paired tumor and normal samples from 653 patients with 23 cancer types from the Pediatric Cancer Genome Project. Telomere dynamics were characterized as the ratio of telomere fractions between tumor and normal samples. Somatic mutations were gathered, RNA sequencing data for 330 patients were analyzed for gene expression, and Cox regression was used to assess the telomere dynamics on patient survival. RESULTS: Telomere lengthening was observed in 28.7% of solid tumors, 10.5% of brain tumors, and 4.3% of hematological cancers. Among 81 samples with telomere lengthening, 26 had somatic mutations in alpha thalassemia/mental retardation syndrome X-linked gene, corroborated by a low level of the gene expression in the subset of tumors with RNA sequencing. Telomerase reverse transcriptase gene amplification and/or activation was observed in 10 tumors with telomere lengthening, including two leukemias of the E2A-PBX1 subtype. Among hematological cancers, pathway analysis for genes with expressions most negatively correlated with telomere fractions suggests the implication of a gene ontology process of antigen presentation by Major histocompatibility complex class II. A higher ratio of telomere fractions was statistically significantly associated with poorer survival for patients with brain tumors (hazard ratio = 2.18, 95% confidence interval = 1.37 to 3.46). CONCLUSION: Because telomerase inhibitors are currently being explored as potential agents to treat pediatric cancer, these data are valuable because they identify a subpopulation of patients with reactivation of telomerase who are most likely to benefit from this novel therapeutic option.

Discovery of regulatory noncoding variants in individual cancer genomes by using cis-X.

We developed cis-X, a computational method for discovering regulatory noncoding variants in cancer by integrating whole-genome and transcriptome sequencing data from a single cancer sample. cis-X first finds aberrantly cis-activated genes that exhibit allele-specific expression accompanied by an elevated outlier expression. It then searches for causal noncoding variants that may introduce aberrant transcription factor binding motifs or enhancer hijacking by structural variations. Analysis of 13 T-lineage acute lymphoblastic leukemias identified a recurrent intronic variant predicted to cis-activate the TAL1 oncogene, a finding validated in vivo by chromatin immunoprecipitation sequencing of a patient-derived xenograft. Candidate oncogenes include the prolactin receptor PRLR activated by a focal deletion that removes a CTCF-insulated neighborhood boundary. cis-X may be applied to pediatric and adult solid tumors that are aneuploid and heterogeneous. In contrast to existing approaches, which require large sample cohorts, cis-X enables the discovery of regulatory noncoding variants in individual cancer genomes.

CICERO: a versatile method for detecting complex and diverse driver fusions using cancer RNA sequencing data.

To discover driver fusions beyond canonical exon-to-exon chimeric transcripts, we develop CICERO, a local assembly-based algorithm that integrates RNA-seq read support with extensive annotation for candidate ranking. CICERO outperforms commonly used methods, achieving a 95% detection rate for 184 independently validated driver fusions including internal tandem duplications and other non-canonical events in 170 pediatric cancer transcriptomes. Re-analysis of TCGA glioblastoma RNA-seq unveils previously unreported kinase fusions (KLHL7-BRAF) and a 13% prevalence of EGFR C-terminal truncation. Accessible via standard or cloud-based implementation, CICERO enhances driver fusion detection for research and precision oncology. The CICERO source code is available at https://github.com/stjude/Cicero.

Pathogenic Germline Mutations in DNA Repair Genes in Combination With Cancer Treatment Exposures and Risk of Subsequent Neoplasms Among Long-Term Survivors of Childhood Cancer.

PURPOSE: To investigate cancer treatment plus pathogenic germline mutations (PGMs) in DNA repair genes (DRGs) for identification of childhood cancer survivors at increased risk of subsequent neoplasms (SNs). METHODS: Whole-genome sequencing was performed on blood-derived DNA from survivors in the St Jude Lifetime Cohort. PGMs were evaluated in 127 genes from 6 major DNA repair pathways. Cumulative doses of chemotherapy and body region-specific radiotherapy (RT) were abstracted from medical records. Relative rates (RRs) and 95% CIs of SNs by mutation status were estimated using multivariable piecewise exponential models. RESULTS: Of 4,402 survivors, 495 (11.2%) developed 1,269 SNs. We identified 538 PGMs in 98 DRGs (POLG, MUTYH, ERCC2, and BRCA2, among others) in 508 (11.5%) survivors. Mutations in homologous recombination (HR) genes were significantly associated with an increased rate of subsequent female breast cancer (RR, 3.7; 95% CI, 1.8 to 7.7), especially among survivors with chest RT ≥ 20 Gy (RR, 4.4; 95% CI, 1.6 to 12.4), or with a cumulative dose of anthracyclines in the second or third tertile (RR, 4.4; 95% CI, 1.7 to 11.4). Mutations in HR genes were also associated with an increased rate of subsequent sarcoma among those who received alkylating agent doses in the third tertile (RR, 14.9; 95% CI, 4.0 to 38.0). Mutations in nucleotide excision repair genes were associated with subsequent thyroid cancer for those treated with neck RT ≥ 30 Gy (RR, 12.9; 95% CI, 1.6 to 46.6) with marginal statistical significance. CONCLUSION: Our study provides novel insights regarding the contribution of genetics, in combination with known treatment-related risks, for the development of SNs. These findings have the potential to facilitate identification of high-risk survivors who may benefit from genetic counseling and/or testing of DRGs, which may further inform personalized cancer surveillance and prevention strategies.

Polygenic Determinants for Subsequent Breast Cancer Risk in Survivors of Childhood Cancer: The St Jude Lifetime Cohort Study (SJLIFE).

PURPOSE: The risk of subsequent breast cancer among female childhood cancer survivors is markedly elevated. We aimed to determine genetic contributions to this risk, focusing on polygenic determinants implicated in breast cancer susceptibility in the general population. EXPERIMENTAL DESIGN: Whole-genome sequencing (30×) was performed on survivors in the St Jude Lifetime Cohort, and germline mutations in breast cancer predisposition genes were classified for pathogenicity. A polygenic risk score (PRS) was constructed for each survivor using 170 established common risk variants. Relative rate (RR) and 95% confidence interval (95% CI) of subsequent breast cancer incidence were estimated using multivariable piecewise exponential regression. RESULTS: The analysis included 1,133 female survivors of European ancestry (median age at last follow-up = 35.4 years; range, 8.4-67.4), of whom 47 were diagnosed with one or more subsequent breast cancers (median age at subsequent breast cancer = 40.3 years; range, 24.5-53.0). Adjusting for attained age, age at primary diagnosis, chest irradiation, doses of alkylating agents and anthracyclines, and genotype eigenvectors, RRs for survivors with PRS in the highest versus lowest quintiles were 2.7 (95% CI, 1.0-7.3), 3.0 (95% CI, 1.1-8.1), and 2.4 (95% CI, 0.1-81.1) for all survivors and survivors with and without chest irradiation, respectively. Similar associations were observed after excluding carriers of pathogenic/likely pathogenic mutations in breast cancer predisposition genes. Notably, the PRS was associated with the subsequent breast cancer rate under the age of 45 years (RR = 3.2; 95% CI, 1.2-8.3). CONCLUSIONS: Genetic profiles comprised of small-effect common variants and large-effect predisposing mutations can inform personalized breast cancer risk and surveillance/intervention in female childhood cancer survivors.

Genetic Risk for Subsequent Neoplasms Among Long-Term Survivors of Childhood Cancer.

Purpose Childhood cancer survivors are at increased risk of subsequent neoplasms (SNs), but the germline genetic contribution is largely unknown. We assessed the contribution of pathogenic/likely pathogenic (P/LP) mutations in cancer predisposition genes to their SN risk. Patients and Methods Whole-genome sequencing (30-fold) was performed on samples from childhood cancer survivors who were ≥ 5 years since initial cancer diagnosis and participants in the St Jude Lifetime Cohort Study, a retrospective hospital-based study with prospective clinical follow-up. Germline mutations in 60 genes known to be associated with autosomal dominant cancer predisposition syndromes with moderate to high penetrance were classified by their pathogenicity according to the American College of Medical Genetics and Genomics guidelines. Relative rates (RRs) and 95% CIs of SN occurrence by mutation status were estimated using multivariable piecewise exponential regression stratified by radiation exposure. Results Participants were 3,006 survivors (53% male; median age, 35.8 years [range, 7.1 to 69.8 years]; 56% received radiotherapy), 1,120 SNs were diagnosed among 439 survivors (14.6%), and 175 P/LP mutations were identified in 5.8% (95% CI, 5.0% to 6.7%) of survivors. Mutations were associated with significantly increased rates of breast cancer (RR, 13.9; 95% CI, 6.0 to 32.2) and sarcoma (RR, 10.6; 95% CI, 4.3 to 26.3) among irradiated survivors and with increased rates of developing any SN (RR, 4.7; 95% CI, 2.4 to 9.3), breast cancer (RR, 7.7; 95% CI, 2.4 to 24.4), nonmelanoma skin cancer (RR, 11.0; 95% CI, 2.9 to 41.4), and two or more histologically distinct SNs (RR, 18.6; 95% CI, 3.5 to 99.3) among nonirradiated survivors. Conclusion The findings support referral of all survivors for genetic counseling for potential clinical genetic testing, which should be prioritized for nonirradiated survivors with any SN and for those with breast cancer or sarcoma in the field of prior irradiation.

The landscape of somatic mutations in epigenetic regulators across 1,000 paediatric cancer genomes.

Studies of paediatric cancers have shown a high frequency of mutation across epigenetic regulators. Here we sequence 633 genes, encoding the majority of known epigenetic regulatory proteins, in over 1,000 paediatric tumours to define the landscape of somatic mutations in epigenetic regulators in paediatric cancer. Our results demonstrate a marked variation in the frequency of gene mutations across 21 different paediatric cancer subtypes, with the highest frequency of mutations detected in high-grade gliomas, T-lineage acute lymphoblastic leukaemia and medulloblastoma, and a paucity of mutations in low-grade glioma and retinoblastoma. The most frequently mutated genes are H3F3A, PHF6, ATRX, KDM6A, SMARCA4, ASXL2, CREBBP, EZH2, MLL2, USP7, ASXL1, NSD2, SETD2, SMC1A and ZMYM3. We identify novel loss-of-function mutations in the ubiquitin-specific processing protease 7 (USP7) in paediatric leukaemia, which result in decreased deubiquitination activity. Collectively, our results help to define the landscape of mutations in epigenetic regulatory genes in paediatric cancer and yield a valuable new database for investigating the role of epigenetic dysregulations in cancer.