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

Gene Rearrangement Detection in Pediatric Leukemia

The detection of gene rearrangements in pediatric leukemia is an essential component of the work-up, with implications for accurate diagnosis, proper risk stratification, and therapeutic decisions, including the use of targeted therapies. The traditional methods of karyotype and fluorescence in situ hybridization are still valuable, but many new assays are also available, with different strengths and weaknesses. These assays include next-generation sequencing-based assays that have the potential for highly multiplexed and/or unbiased detection of rearrangements.

Rhabdomyosarcoma: How Advanced Molecular Methods Are Shaping the Diagnostic and Therapeutic Paradigm

For the past 40 years, progress in rhabdomyosarcoma (RMS) has been focused on understanding its molecular basis and characterizing the mutations that drive its tumorigenesis and progression. Genetic predisposition to RMS has allowed discovery of key genetic pathways and driver mutations. Subclassification of RMS into embryonal (ERMS) and alveolar (ARMS) subtypes has shifted from histology to PAX-FOXO1 fusion status, and new driver mutations have been found in spindle cell RMS. Comprehensive molecular profiling leveraging genome-scale next-generation sequencing (NGS) indicates that the RAS/RAF/PI3K axis is mutated in the majority of ERMS and modulated by downstream effects of PAX-FOXO1 fusions in ARMS. Because of the continued poor outcome of high-risk RMS, a variety of molecular targets have been or are now being tested in current or recent therapy trials. New techniques such as single cell sequencing, spatial multi-omics, and CRISPR/Cas9 genome editing offer potential for further discovery, but a need for clinically annotated specimens persists.

Infectious triggers and novel therapeutic opportunities in childhood B cell leukaemia

B cell acute lymphoblastic leukaemia (B-ALL) is the most common form of childhood cancer. Although treatment has advanced remarkably in the past 50 years, it still fails in ~20% of patients. Recent studies revealed that more than 5% of healthy newborns carry preleukaemic clones that originate in utero, but only a small percentage of these carriers will progress to overt B-ALL. The drivers of progression are unclear, but B-ALL incidence seems to be increasing in parallel with the adoption of modern lifestyles. Emerging evidence shows that a major driver for the conversion from the preleukaemic state to the B-ALL state is exposure to immune stressors, such as infection. Here, we discuss our current understanding of the environmental triggers and genetic predispositions that may lead to B-ALL, highlighting lessons from epidemiology, the clinic and animal models, and identifying priority areas for future research.

The Landscape of Pediatric Precision Oncology: Program Design, Actionable Alterations, and Clinical Trial Development

Over the last years, various precision medicine programs have been developed for pediatric patients with high-risk, relapsed, or refractory malignancies, selecting patients for targeted treatment through comprehensive molecular profiling. In this review, we describe characteristics of these initiatives, demonstrating the feasibility and potential of molecular-driven precision medicine. Actionable events are identified in a significant subset of patients, although comparing results is complicated due to the lack of a standardized definition of actionable alterations and the different molecular profiling strategies used. The first biomarker-driven trials for childhood cancer have been initiated, but until now the effect of precision medicine on clinical outcome has only been reported for a small number of patients, demonstrating clinical benefit in some. Future perspectives include the incorporation of novel approaches such as liquid biopsies and immune monitoring as well as innovative collaborative trial design including combination strategies, and the development of agents specifically targeting aberrations in childhood malignancies.

TET2 as a tumor suppressor and therapeutic target in T-cell acute lymphoblastic leukemia

Pediatric T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive malignancy in need of novel targeted therapies to prevent relapse and lessen treatment toxicity. We reveal frequent (∼88%) transcriptional silencing or repression of the tumor suppressor TET2 in T-ALL. We show that loss of TET2 in T-ALL is correlated with hypermethylation of the TET2 promoter and that TET2 expression can be rescued by treatment with the DNA demethylating agent, 5-azacytidine (5-aza). We further reveal that the TET2 cofactor vitamin C exerts a strong synergistic effect on global transcriptional changes when added to 5-aza treatment. Importantly, 5-aza treatment results in increased cell death, specifically in T-ALL cells lacking TET2. Thus, we clearly identify 5-aza as a potentially targeted therapy for TET2-silenced T-ALL.

Pediatric T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive malignancy resulting from overproduction of immature T-cells in the thymus and is typified by widespread alterations in DNA methylation. As survival rates for relapsed T-ALL remain dismal (10 to 25%), development of targeted therapies to prevent relapse is key to improving prognosis. Whereas mutations in the DNA demethylating enzyme TET2 are frequent in adult T-cell malignancies, TET2 mutations in T-ALL are rare. Here, we analyzed RNA-sequencing data of 321 primary T-ALLs, 20 T-ALL cell lines, and 25 normal human tissues, revealing that TET2 is transcriptionally repressed or silenced in 71% and 17% of T-ALL, respectively. Furthermore, we show that TET2 silencing is often associated with hypermethylation of the TET2 promoter in primary T-ALL. Importantly, treatment with the DNA demethylating agent, 5-azacytidine (5-aza), was significantly more toxic to TET2-silenced T-ALL cells and resulted in stable re-expression of the TET2 gene. Additionally, 5-aza led to up-regulation of methylated genes and human endogenous retroviruses (HERVs), which was further enhanced by the addition of physiological levels of vitamin C, a potent enhancer of TET activity. Together, our results clearly identify 5-aza as a potential targeted therapy for TET2-silenced T-ALL.

Selective Targeting of Class I Histone Deacetylases in a Model of Human Osteosarcoma

Dysregulation of histone deacetylases (HDACs) is associated with the pathogenesis of human osteosarcoma, which may present an epigenetic vulnerability as well as a therapeutic target. Domatinostat (4SC-202) is a next-generation class I HDAC inhibitor that is currently being used in clinical research for certain cancers, but its impact on human osteosarcoma has yet to be explored. In this study, we report that 4SC-202 inhibits osteosarcoma cell growth in vitro and in vivo. By analyzing cell function in vitro, we show that the anti-tumor effect of 4SC-202 involves the combined induction of cell-cycle arrest at the G2/M phase and apoptotic program, as well as a reduction in cell invasion and migration capabilities. We also found that 4SC-202 has little capacity to promote osteogenic differentiation. Remarkably, 4SC-202 revised the global transcriptome and induced distinct signatures of gene expression in vitro. Moreover, 4SC-202 decreased tumor growth of established human tumor xenografts in immunodeficient mice in vivo. We further reveal key targets regulated by 4SC-202 that contribute to tumor cell growth and survival, and canonical signaling pathways associated with progression and metastasis of osteosarcoma. Our study suggests that 4SC-202 may be exploited as a valuable drug to promote more effective treatment of patients with osteosarcoma and provide molecular insights into the mechanism of action of class I HDAC inhibitors.

The Mutational Landscape of Myeloid Leukaemia in Down Syndrome

Children with Down syndrome (DS) are particularly prone to haematopoietic disorders. Paediatric myeloid malignancies in DS occur at an unusually high frequency and generally follow a well-defined stepwise clinical evolution. First, the acquisition of mutations in the GATA1 transcription factor gives rise to a transient myeloproliferative disorder (TMD) in DS newborns. While this condition spontaneously resolves in most cases, some clones can acquire additional mutations, which trigger myeloid leukaemia of Down syndrome (ML-DS). These secondary mutations are predominantly found in chromatin and epigenetic regulators-such as cohesin, CTCF or EZH2-and in signalling mediators of the JAK/STAT and RAS pathways. Most of them are also found in non-DS myeloid malignancies, albeit at extremely different frequencies. Intriguingly, mutations in proteins involved in the three-dimensional organization of the genome are found in nearly 50% of cases. How the resulting mutant proteins cooperate with trisomy 21 and mutant GATA1 to promote ML-DS is not fully understood. In this review, we summarize and discuss current knowledge about the sequential acquisition of genomic alterations in ML-DS.

Fusion Oncoproteins in Childhood Cancers: Potential Role in Targeted Therapy

Cancer remains the leading cause of death from disease in children. Historically, in contrast to their adult counterparts, the causes of pediatric malignancies have remained largely unknown, with most pediatric cancers displaying low mutational burdens. Research related to molecular genetics in pediatric cancers is advancing our understanding of potential drivers of tumorigenesis and opening new opportunities for targeted therapies. One such area is fusion oncoproteins, which are a product of chromosomal rearrangements resulting in the fusion of different genes. They have been identified as oncogenic drivers in several sarcomas and leukemias. Continued advancement in the understanding of the biology of fusion oncoproteins will contribute to the discovery and development of new therapies for childhood cancers. Here we review the current scientific knowledge on fusion oncoproteins, focusing on pediatric sarcomas and hematologic cancers, and highlight the challenges and current efforts in developing drugs to target fusion oncoproteins.

Binary pan-cancer classes with distinct vulnerabilities defined by pro- or anti-cancer YAP/TEAD activity

Cancer heterogeneity impacts therapeutic response, driving efforts to discover over-arching rules that supersede variability. Here, we define pan-cancer binary classes based on distinct expression of YAP and YAP-responsive adhesion regulators. Combining informatics with in vivo and in vitro gain- and loss-of-function studies across multiple murine and human tumor types, we show that opposite pro- or anti-cancer YAP activity functionally defines binary YAP^on or YAP^off cancer classes that express or silence YAP, respectively. YAP^off solid cancers are neural/neuroendocrine and frequently RB1^-/-, such as retinoblastoma, small cell lung cancer, and neuroendocrine prostate cancer. YAP silencing is intrinsic to the cell of origin, or acquired with lineage switching and drug resistance. The binary cancer groups exhibit distinct YAP-dependent adhesive behavior and pharmaceutical vulnerabilities, underscoring clinical relevance. Mechanistically, distinct YAP/TEAD enhancers in YAP^off or YAP^on cancers deploy anti-cancer integrin or pro-cancer proliferative programs, respectively. YAP is thus pivotal across cancer, but in opposite ways, with therapeutic implications.

svpluscnv: analysis and visualization of complex structural variation data

MOTIVATION

Despite widespread prevalence of somatic structural variations (SVs) across most tumor types, understanding of their molecular implications often remains poor. SVs are extremely heterogeneous in size and complexity, hindering the interpretation of their pathogenic role. Tools integrating large SV datasets across platforms are required to fully characterize the cancer's somatic landscape.

RESULTS

svpluscnv R package is a swiss army knife for the integration and interpretation of orthogonal datasets including copy number variant segmentation profiles and sequencing-based structural variant calls. The package implements analysis and visualization tools to evaluate chromosomal instability and ploidy, identify genes harboring recurrent SVs and detects complex rearrangements such as chromothripsis and chromoplexia. Further, it allows systematic identification of hot-spot shattered genomic regions, showing reproducibility across alternative detection methods and datasets.

AVAILABILITY AND IMPLEMENTATION

https://github.com/ccbiolab/svpluscnv.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Early-phenotype CAR-T cells for the treatment of pediatric cancers

Chimeric antigen receptor (CAR) T cell therapy is a promising approach for the treatment of childhood cancers, particularly high-risk tumors that fail to respond to standard therapies. CAR-T cells have been highly successful in treating some types of hematological malignancies. However, CAR-T cells targeting solid cancers have had limited success so far for multiple reasons, including their poor long-term persistence and proliferation. Evidence is emerging to show that maintaining CAR-T cells in an early, less differentiated state in vitro results in superior persistence, proliferation, and anti-tumor effects in vivo. Children are ideal candidates for receiving less-differentiated CAR-T cells, because their peripheral T cell pool primarily comprises naïve cells that could readily be harvested in large numbers to generate early-phenotype CAR-T cells. Although several studies have reported different approaches to successfully generate early CAR-T cells, there are only a few clinical trials testing these in adult patients. No trials are currently testing early CAR-T cells in children. Here, we summarize the different strategies used to maintain CAR-T cells in an early phenotypic stage, and present evidence suggesting that this approach may be particularly relevant to treating childhood cancers.

Mutational patterns and clonal evolution from diagnosis to relapse in pediatric acute lymphoblastic leukemia

The mechanisms driving clonal heterogeneity and evolution in relapsed pediatric acute lymphoblastic leukemia (ALL) are not fully understood. We performed whole genome sequencing of samples collected at diagnosis, relapse(s) and remission from 29 Nordic patients. Somatic point mutations and large-scale structural variants were called using individually matched remission samples as controls, and allelic expression of the mutations was assessed in ALL cells using RNA-sequencing. We observed an increased burden of somatic mutations at relapse, compared to diagnosis, and at second relapse compared to first relapse. In addition to 29 known ALL driver genes, of which nine genes carried recurrent protein-coding mutations in our sample set, we identified putative non-protein coding mutations in regulatory regions of seven additional genes that have not previously been described in ALL. Cluster analysis of hundreds of somatic mutations per sample revealed three distinct evolutionary trajectories during ALL progression from diagnosis to relapse. The evolutionary trajectories provide insight into the mutational mechanisms leading relapse in ALL and could offer biomarkers for improved risk prediction in individual patients.

TLE4 Is a Critical Mediator of Osteoblast and Runx2-Dependent Bone Development

Healthy bone homeostasis hinges upon a delicate balance and regulation of multiple processes that contribute to bone development and metabolism. While examining hematopoietic regulation by Tle4, we have uncovered a previously unappreciated role of Tle4 on bone calcification using a novel Tle4 null mouse model. Given the significance of osteoblasts in both hematopoiesis and bone development, this study investigated how loss of Tle4 affects osteoblast function. We used dynamic bone formation parameters and microCT to characterize the adverse effects of Tle4 loss on bone development. We further demonstrated loss of Tle4 impacts expression of several key osteoblastogenic genes, including Runx2, Oc, and Ap, pointing toward a potential novel mechanism for Tle4-dependent regulation of mammalian bone development in collaboration with the RUNX family members.

Canine tumor mutational burden is correlated with TP53 mutation across tumor types and breeds

Spontaneous canine cancers are valuable but relatively understudied and underutilized models. To enhance their usage, we reanalyze whole exome and genome sequencing data published for 684 cases of >7 common tumor types and >35 breeds, with rigorous quality control and breed validation. Our results indicate that canine tumor alteration landscape is tumor type-dependent, but likely breed-independent. Each tumor type harbors major pathway alterations also found in its human counterpart (e.g., PI3K in mammary tumor and p53 in osteosarcoma). Mammary tumor and glioma have lower tumor mutational burden (TMB) (median < 0.5 mutations per Mb), whereas oral melanoma, osteosarcoma and hemangiosarcoma have higher TMB (median ≥ 1 mutations per Mb). Across tumor types and breeds, TMB is associated with mutation of TP53 but not PIK3CA, the most mutated genes. Golden Retrievers harbor a TMB-associated and osteosarcoma-enriched mutation signature. Here, we provide a snapshot of canine mutations across major tumor types and breeds.

Cancer Informatics for Cancer Centers: Scientific Drivers for Informatics, Data Science, and Care in Pediatric, Adolescent, and Young Adult Cancer

Cancer Informatics for Cancer Centers (CI4CC) is a grassroots, nonprofit 501c3 organization intended to provide a focused national forum for engagement of senior cancer informatics leaders, primarily aimed at academic cancer centers anywhere in the world but with a special emphasis on the 70 National Cancer Institute-funded cancer centers. This consortium has regularly held topic-focused biannual face-to-face symposiums. These meetings are a place to review cancer informatics and data science priorities and initiatives, providing a forum for discussion of the strategic and pragmatic issues that we faced at our respective institutions and cancer centers. Here, we provide meeting highlights from the latest CI4CC Symposium, which was delayed from its original April 2020 schedule because of the COVID-19 pandemic and held virtually over three days (September 24, October 1, and October 8) in the fall of 2020. In addition to the content presented, we found that holding this event virtually once a week for 6 hours was a great way to keep the kind of deep engagement that a face-to-face meeting engenders. This is the second such publication of CI4CC Symposium highlights, the first covering the meeting that took place in Napa, California, from October 14-16, 2019. We conclude with some thoughts about using data science to learn from every child with cancer, focusing on emerging activities of the National Cancer Institute's Childhood Cancer Data Initiative.

Gene Fusions Create Partner and Collateral Dependencies Essential to Cancer Cell Survival

Gene fusions frequently result from rearrangements in cancer genomes. In many instances, gene fusions play an important role in oncogenesis; in other instances, they are thought to be passenger events. Although regulatory element rearrangements and copy number alterations resulting from these structural variants are known to lead to transcriptional dysregulation across cancers, the extent to which these events result in functional dependencies with an impact on cancer cell survival is variable. Here we used CRISPR-Cas9 dependency screens to evaluate the fitness impact of 3,277 fusions across 645 cell lines from the Cancer Dependency Map. We found that 35% of cell lines harbored either a fusion partner dependency or a collateral dependency on a gene within the same topologically associating domain as a fusion partner. Fusion-associated dependencies revealed numerous novel oncogenic drivers and clinically translatable alterations. Broadly, fusions can result in partner and collateral dependencies that have biological and clinical relevance across cancer types. SIGNIFICANCE: This study provides insights into how fusions contribute to fitness in different cancer contexts beyond partner-gene activation events, identifying partner and collateral dependencies that may have direct implications for clinical care.

Chemotherapy and mismatch repair deficiency cooperate to fuel TP53 mutagenesis and ALL relapse

Chemotherapy is a standard treatment for pediatric acute lymphoblastic leukemia (ALL), which sometimes relapses with chemoresistant features. However, whether acquired drug-resistance mutations in relapsed ALL pre-exist or are induced by treatment remains unknown. Here we provide direct evidence of a specific mechanism by which chemotherapy induces drug-resistance-associated mutations leading to relapse. Using genomic and functional analysis of relapsed ALL we show that thiopurine treatment in mismatch repair (MMR)-deficient leukemias induces hotspot TP53 R248Q mutations through a specific mutational signature (thio-dMMR). Clonal evolution analysis reveals sequential MMR inactivation followed by TP53 mutation in some patients with ALL. Acquired TP53 R248Q mutations are associated with on-treatment relapse, poor treatment response and resistance to multiple chemotherapeutic agents, which could be reversed by pharmacological p53 reactivation. Our findings indicate that TP53 R248Q in relapsed ALL originates through synergistic mutagenesis from thiopurine treatment and MMR deficiency and suggest strategies to prevent or treat TP53-mutant relapse.

Functional and epigenetic phenotypes of humans and mice with DNMT3A Overgrowth Syndrome

Germline pathogenic variants in DNMT3A were recently described in patients with overgrowth, obesity, behavioral, and learning difficulties (DNMT3A Overgrowth Syndrome/DOS). Somatic mutations in the DNMT3A gene are also the most common cause of clonal hematopoiesis, and can initiate acute myeloid leukemia (AML). Using whole genome bisulfite sequencing, we studied DNA methylation in peripheral blood cells of 11 DOS patients and found a focal, canonical hypomethylation phenotype, which is most severe with the dominant negative DNMT3AR882H mutation. A germline mouse model expressing the homologous Dnmt3aR878H mutation phenocopies most aspects of the human DOS syndrome, including the methylation phenotype and an increased incidence of spontaneous hematopoietic malignancies, suggesting that all aspects of this syndrome are caused by this mutation.

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.

CLEC12A and CD33 coexpression as a preferential target for pediatric AML combinatorial immunotherapy

Emerging immunotherapies such as chimeric antigen receptor T cells have advanced the treatment of acute lymphoblastic leukemia. In contrast, long-term control of acute myeloid leukemia (AML) cannot be achieved by single lineage-specific targeting while sparing benign hematopoiesis. In addition, heterogeneity of AML warrants combinatorial targeting, and several suitable immunotargets (HAVCR2/CD33 and HAVCR2/CLEC12A) have been identified in adult AML. However, clinical and biologic characteristics of AML differ between children and the elderly. Here, we analyzed 36 bone marrow (BM) samples of pediatric AML patients and 13 age-matched healthy donors using whole RNA sequencing of sorted CD45dim and CD34+CD38-CD45dim BM populations and flow cytometry for surface expression of putative target antigens. Pediatric AML clusters apart from healthy myeloid BM precursors in principal-component analysis. Known immunotargets of adult AML, such as IL3RA, were not overexpressed in pediatric AML compared with healthy precursors by RNA sequencing. CD33 and CLEC12A were the most upregulated immunotargets on the RNA level and showed the highest surface expression on AML detected by flow cytometry. KMT2A-mutated infant AML clusters separately by RNA sequencing and overexpresses FLT3, and hence, CD33/FLT3 cotargeting is an additional specific option for this subgroup. CLEC12A and CD33/CLEC12Adouble-positive expression was absent in CD34+CD38-CD45RA-CD90+ hematopoietic stem cells (HSCs) and nonhematopoietic tissue, while CD33 and FLT3 are expressed on HSCs. In summary, we show that expression of immunotargets in pediatric AML differs from known expression profiles in adult AML. We identify CLEC12A and CD33 as preferential generic combinatorial immunotargets in pediatric AML and CD33 and FLT3 as immunotargets specific for KMT2A-mutated infant AML.

The Second Oncogenic Hit Determines the Cell Fate of ETV6-RUNX1 Positive Leukemia

ETV6-RUNX1 is almost exclusively associated with childhood B-cell acute lymphoblastic leukemia (B-ALL), but the consequences of ETV6-RUNX1 expression on cell lineage decisions during B-cell leukemogenesis are completely unknown. Clinically silent ETV6-RUNX1 preleukemic clones are frequently found in neonatal cord blood, but few carriers develop B-ALL as a result of secondary genetic alterations. The understanding of the mechanisms underlying the first transforming steps could greatly advance the development of non-toxic prophylactic interventions. Using genetic lineage tracing, we examined the capacity of ETV6-RUNX1 to instruct a malignant phenotype in the hematopoietic lineage by cell-specific Cre-mediated activation of ETV6-RUNX1 from the endogenous Etv6 gene locus. Here we show that, while ETV6-RUNX1 has the propensity to trigger both T- and B-lymphoid malignancies, it is the second hit that determines tumor cell identity. To instigate leukemia, both oncogenic hits must place early in the development of hematopoietic/precursor cells, not in already committed B-cells. Depending on the nature of the second hit, the resulting B-ALLs presented distinct entities that were clearly separable based on their gene expression profiles. Our findings give a novel mechanistic insight into the early steps of ETV6-RUNX1+ B-ALL development and might have major implications for the potential development of ETV6-RUNX1+ B-ALL prevention strategies.

Haplotype-resolved germline and somatic alterations in renal medullary carcinomas

BACKGROUND

Renal medullary carcinomas (RMCs) are rare kidney cancers that occur in adolescents and young adults of African ancestry. Although RMC is associated with the sickle cell trait and somatic loss of the tumor suppressor, SMARCB1, the ancestral origins of RMC remain unknown. Further, characterization of structural variants (SVs) involving SMARCB1 in RMC remains limited.

METHODS

We used linked-read genome sequencing to reconstruct germline and somatic haplotypes in 15 unrelated patients with RMC registered on the Children's Oncology Group (COG) AREN03B2 study between 2006 and 2017 or from our prior study. We performed fine-mapping of the HBB locus and assessed the germline for cancer predisposition genes. Subsequently, we assessed the tumor samples for mutations outside of SMARCB1 and integrated RNA sequencing to interrogate the structural variants at the SMARCB1 locus.

RESULTS

We find that the haplotype of the sickle cell mutation in patients with RMC originated from three geographical regions in Africa. In addition, fine-mapping of the HBB locus identified the sickle cell mutation as the sole candidate variant. We further identify that the SMARCB1 structural variants are characterized by blunt or 1-bp homology events.

CONCLUSIONS

Our findings suggest that RMC does not arise from a single founder population and that the HbS allele is a strong candidate germline allele which confers risk for RMC. Furthermore, we find that the SVs that disrupt SMARCB1 function are likely repaired by non-homologous end-joining. These findings highlight how haplotype-based analyses using linked-read genome sequencing can be applied to identify potential risk variants in small and rare disease cohorts and provide nucleotide resolution to structural variants.

Future Match Making: When Pediatric Oncology Meets Organoid Technology

Unlike adult cancers that frequently result from the accumulation in time of mutational “hits” often linked to lifestyle, childhood cancers are emerging as diseases of dysregulated development through massive epigenetic alterations. The ability to reconstruct these differences in cancer models is therefore crucial for better understanding the uniqueness of pediatric cancer biology. Cancer organoids (i.e., tumoroids) represent a promising approach for creating patient-derived in vitro cancer models that closely recapitulate the overall pathophysiological features of natural tumorigenesis, including intra-tumoral heterogeneity and plasticity. Though largely applied to adult cancers, this technology is scarcely used for childhood cancers, with a notable delay in technological transfer. However, tumoroids could provide an unprecedented tool to unravel the biology of pediatric cancers and improve their therapeutic management. We herein present the current state-of-the-art of a long awaited and much needed matchmaking.

Epidemiological overview of multidimensional chromosomal and genome toxicity of cannabis exposure in congenital anomalies and cancer development

Cannabis and cannabinoids are implicated in multiple genotoxic, epigenotoxic and chromosomal-toxic mechanisms and interact with several morphogenic pathways, likely underpinning previous reports of links between cannabis and congenital anomalies and heritable tumours. However the effects of cannabinoid genotoxicity have not been assessed on whole populations and formal consideration of effects as a broadly acting genotoxin remain unexplored. Our study addressed these knowledge gaps in USA datasets. Cancer data from CDC, drug exposure data from National Survey of Drug Use and Health 2003-2017 and congenital anomaly data from National Birth Defects Prevention Network were used. We show that cannabis, THC cannabigerol and cannabichromene exposure fulfill causal criteria towards first Principal Components of both: (A) Down syndrome, Trisomies 18 and 13, Turner syndrome, Deletion 22q11.2, and (B) thyroid, liver, breast and pancreatic cancers and acute myeloid leukaemia, have mostly medium to large effect sizes, are robust to adjustment for ethnicity, other drugs and income in inverse probability-weighted models, show prominent non-linear effects, have 55/56 e-Values > 1.25, and are exacerbated by cannabis liberalization (P = 9.67 × 10^-43, 2.66 × 10^-15). The results confirm experimental studies showing that cannabinoids are an important cause of community-wide genotoxicity impacting both birth defect and cancer epidemiology at the chromosomal hundred-megabase level.

Epidemiological overview of multidimensional chromosomal and genome toxicity of cannabis exposure in congenital anomalies and cancer development

Cannabis and cannabinoids are implicated in multiple genotoxic, epigenotoxic and chromosomal-toxic mechanisms and interact with several morphogenic pathways, likely underpinning previous reports of links between cannabis and congenital anomalies and heritable tumours. However the effects of cannabinoid genotoxicity have not been assessed on whole populations and formal consideration of effects as a broadly acting genotoxin remain unexplored. Our study addressed these knowledge gaps in USA datasets. Cancer data from CDC, drug exposure data from National Survey of Drug Use and Health 2003-2017 and congenital anomaly data from National Birth Defects Prevention Network were used. We show that cannabis, THC cannabigerol and cannabichromene exposure fulfill causal criteria towards first Principal Components of both: (A) Down syndrome, Trisomies 18 and 13, Turner syndrome, Deletion 22q11.2, and (B) thyroid, liver, breast and pancreatic cancers and acute myeloid leukaemia, have mostly medium to large effect sizes, are robust to adjustment for ethnicity, other drugs and income in inverse probability-weighted models, show prominent non-linear effects, have 55/56 e-Values > 1.25, and are exacerbated by cannabis liberalization (P = 9.67 × 10^-43, 2.66 × 10^-15). The results confirm experimental studies showing that cannabinoids are an important cause of community-wide genotoxicity impacting both birth defect and cancer epidemiology at the chromosomal hundred-megabase level.

Network Analysis Reveals Synergistic Genetic Dependencies for Rational Combination Therapy in Philadelphia Chromosome-Like Acute Lymphoblastic Leukemia

PURPOSE

Systems biology approaches can identify critical targets in complex cancer signaling networks to inform new therapy combinations that may overcome conventional treatment resistance.

EXPERIMENTAL DESIGN

We performed integrated analysis of 1,046 childhood B-ALL cases and developed a data-driven network controllability-based approach to identify synergistic key regulator targets in Philadelphia chromosome-like B-acute lymphoblastic leukemia (Ph-like B-ALL), a common high-risk leukemia subtype associated with hyperactive signal transduction and chemoresistance.

RESULTS

We identified 14 dysregulated network nodes in Ph-like ALL involved in aberrant JAK/STAT, Ras/MAPK, and apoptosis pathways and other critical processes. Genetic cotargeting of the synergistic key regulator pair STAT5B and BCL2-associated athanogene 1 (BAG1) significantly reduced leukemia cell viability in vitro. Pharmacologic inhibition with dual small molecule inhibitor therapy targeting this pair of key nodes further demonstrated enhanced antileukemia efficacy of combining the BCL-2 inhibitor venetoclax with the tyrosine kinase inhibitors ruxolitinib or dasatinib in vitro in human Ph-like ALL cell lines and in vivo in multiple childhood Ph-like ALL patient-derived xenograft models. Consistent with network controllability theory, co-inhibitor treatment also shifted the transcriptomic state of Ph-like ALL cells to become less like kinase-activated BCR-ABL1-rearranged (Ph+) B-ALL and more similar to prognostically favorable childhood B-ALL subtypes.

CONCLUSIONS

Our study represents a powerful conceptual framework for combinatorial drug discovery based on systematic interrogation of synergistic vulnerability pathways with pharmacologic inhibitor validation in preclinical human leukemia models.

Degradation of CCNK/CDK12 is a druggable vulnerability of colorectal cancer

Novel treatment options for metastatic colorectal cancer (CRC) are urgently needed to improve patient outcome. Here, we screen a library of non-characterized small molecules against a heterogeneous collection of patient-derived CRC spheroids. By prioritizing compounds with inhibitory activity in a subset of-but not all-spheroid cultures, NCT02 is identified as a candidate with minimal risk of non-specific toxicity. Mechanistically, we show that NCT02 acts as molecular glue that induces ubiquitination of cyclin K (CCNK) and proteasomal degradation of CCNK and its complex partner CDK12. Knockout of CCNK or CDK12 decreases proliferation of CRC cells in vitro and tumor growth in vivo. Interestingly, sensitivity to pharmacological CCNK/CDK12 degradation is associated with TP53 deficiency and consensus molecular subtype 4 in vitro and in patient-derived xenografts. We thus demonstrate the efficacy of targeted CCNK/CDK12 degradation for a CRC subset, highlighting the potential of drug-induced proteolysis for difficult-to-treat types of cancer.

Single-cell DNA amplicon sequencing reveals clonal heterogeneity and evolution in T-cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive leukemia that is most frequent in children and is characterized by the presence of few chromosomal rearrangements and 10 to 20 somatic mutations in protein-coding regions at diagnosis. The majority of T-ALL cases harbor activating mutations in NOTCH1 together with mutations in genes implicated in kinase signaling, transcriptional regulation, or protein translation. To obtain more insight in the level of clonal heterogeneity at diagnosis and during treatment, we used single-cell targeted DNA sequencing with the Tapestri platform. We designed a custom ALL panel and obtained accurate single-nucleotide variant and small insertion-deletion mutation calling for 305 amplicons covering 110 genes in about 4400 cells per sample and time point. A total of 108 188 cells were analyzed for 25 samples of 8 T-ALL patients. We typically observed a major clone at diagnosis (>35% of the cells) accompanied by several minor clones of which some were less than 1% of the total number of cells. Four patients had >2 NOTCH1 mutations, some of which present in minor clones, indicating a strong pressure to acquire NOTCH1 mutations in developing T-ALL cells. By analyzing longitudinal samples, we detected the presence and clonal nature of residual leukemic cells and clones with a minor presence at diagnosis that evolved to clinically relevant major clones at later disease stages. Therefore, single-cell DNA amplicon sequencing is a sensitive assay to detect clonal architecture and evolution in T-ALL.

Navigating the Regulatory Landscape to Develop Pediatric Oncology Drugs: Expert Opinion Recommendations

Regulatory changes have been enacted in the United States (US) and European Union (EU) to encourage the development of new treatments for pediatric cancer. Here, we review some of the factors that have hampered the development of pediatric cancer treatments and provide a comparison of the US and EU regulations implemented to address this clinical need. We then provide some recommendations for each stage of the oncology drug development pathway to help researchers maximize their chance of successful drug development while complying with regulations. A key recommendation is the engagement of key stakeholders such as regulatory authorities, pediatric oncologists, academic researchers, patient advocacy groups, and a Pediatric Expert Group early in the drug development process. During drug target selection, sponsors are encouraged to consult the Food and Drug Administration (FDA), European Medicines Agency (EMA), and the FDA target list, in addition to relevant US and European consortia that have been established to characterize and prioritize oncology drug targets. Sponsors also need to carefully consider the resourcing requirements for preclinical testing, which include ensuring appropriate access to the most relevant databases, clinical samples, and preclinical models (cell lines and animal models). During clinical development, sponsors can account for the pharmacodynamic (PD)/pharmacokinetic (PK) considerations specific to a pediatric population by developing pediatric formulations, selecting suitable PD endpoints, and employing sparse PK sampling or modeling/simulation of drug exposures where appropriate. Additional clinical considerations include the specific design of the clinical trial, the potential inclusion of children in adult trials, and the value of cooperative group trials.

Scientific considerations for global drug development

Requiring regional or in-country confirmatory clinical trials before approval of drugs already approved elsewhere delays access to medicines in low- and middle-income countries and raises drug costs. Here, we discuss the scientific and technological advances that may reduce the need for in-country or in-region clinical trials for drugs approved in other countries and limitations of these advances that could necessitate in-region clinical studies.

Personalized prediction of overall survival in patients with AML in non-complete remission undergoing allo-HCT

Allogenic hematopoietic stem cell transplantation (allo‐HCT) is the standard treatment for acute myeloid leukemia (AML) in non‐complete remission (non‐CR); however, the prognosis is inconsistent. This study aimed to develop and validate nomograms and a web application to predict the overall survival (OS) of patients with non‐CR AML undergoing allo‐HCT (cord blood transplantation [CBT], bone marrow transplantation [BMT], and peripheral blood stem cell transplantation [PBSCT]). Data from 3052 patients were analyzed to construct and validate the prognostic models. The common significant prognostic factors among patients undergoing allo‐HCT were age, performance status, percentage of peripheral blasts, cytogenetic risk, chemotherapy response, and number of transplantations. The conditioning regimen was a significant prognostic factor only in patients undergoing CBT. Compared with cyclophosphamide/total body irradiation, a conditioning regimen of ≥3 drugs, including fludarabine, with CBT exhibited the lowest hazard ratio for mortality (0.384; 95% CI, 0.266–0.554; p < 0.0001). A conditioning regimen of ≥3 drugs with CBT also showed the best leukemia‐free survival among all conditioning regimens. Based on the results of the multivariable analysis, we developed prognostic models showing adequate calibration and discrimination (the c‐indices for CBT, BMT, and PBSCT were 0.648, 0.600, and 0.658, respectively). Our prognostic models can help in assessing individual risks and designing future clinical studies. Furthermore, our study indicates the effectiveness of multi‐drug conditioning regimens in patients undergoing CBT.

Age related gene DST represents an independent prognostic factor for MYCN non-amplified neuroblastoma

Background

MYCN amplification and age are two critical prognostic factors of pediatric neuroblastoma. Previously, we had revealed the prognosis of MYCN target genes. However, the prognostic effects of age related genes in neuroblastoma are unclear.

Methods

The prognostic significance of age and MYCN amplification was determined through multivariate cox regression and Kaplan-Meier survival analysis. Genes differentially expressed in MYCN non-amplified younger neuroblastoma patients were identified using Therapeutically Applicable Research to Generate Effective Treatments (TARGET) and Gene Expression Omnibus (GEO) datasets. The prognostic effects of age related genes ALCAM, CACNA2D3, DST, EPB41L4A and KIF1B in pediatric neuroblastoma patients were determined by Kaplan-Meier survival.

Results

In a pediatric pan-cancer analysis, age was associated with the overall survival of pediatric B-lineage acute lymphoblastic leukemia, neuroblastoma and wilms tumor in TARGET dataset. Moreover, the prognostic effects of age in neuroblastoma were validated using two independent neuroblastoma cohorts. Furthermore, age and MYCN amplification were independent prognostic factors in pediatric neuroblastoma. Compared with MYCN non-amplified older neuroblastoma patients, MYCN non-amplified younger neuroblastoma patients had better clinical outcomes. ALCAM, CACNA2D3, DST, EPB41L4A and KIF1B were highly expressed in MYCN non-amplified younger neuroblastoma patients. And the higher expression levels of ALCAM, CACNA2D3, DST, EPB41L4A or KIF1B were associated with better prognosis of MYCN non-amplified neuroblastoma patients. DST was an independent prognostic factor in MYCN non-amplified neuroblastoma patients and MYCN non-amplified neuroblastoma younger patients with higher DST expression levels had the best clinical overall survival.

Conclusions

Age related gene DST was an independent prognostic factor in MYCN non-amplified neuroblastoma. MYCN non-amplified younger neuroblastoma patients with higher DST expression levels had the best clinical overall survival.

Using primary teeth and archived dried spots for exposomic studies in children: Exploring new paths in the environmental epidemiology of pediatric cancer

It is estimated that 300,000 children 0-14 years of age are diagnosed with cancer worldwide each year. While the absolute risk of cancer in children is low, it is the leading cause of death due to disease in children in high-income countries. In spite of this, the etiologies of pediatric cancer are largely unknown. Environmental exposures have long been thought to play an etiologic role. However, to date, there are few well-established environmental risk factors for pediatric malignancies, likely due to technical barriers in collecting biological samples prospectively in pediatric populations for direct measurements. In this review, we propose the use of novel or underutilized biospecimens (dried blood spots and teeth) and molecular approaches for exposure assessment (epigenetics, metabolomics, and somatic mutational profiles). Future epidemiologic studies of pediatric cancer should incorporate novel exposure assessment methodologies, data on molecular features of tumors, and a more complete assessment of gene-environment interactions.

A review of the biological and clinical implications of RAS-MAPK pathway alterations in neuroblastoma

Neuroblastoma is the most common extra-cranial solid tumor in children, representing approximately 8% of all malignant childhood tumors and 15% of pediatric cancer-related deaths. Recent sequencing and transcriptomics studies have demonstrated the RAS-MAPK pathway’s contribution to the development and progression of neuroblastoma. This review compiles up-to-date evidence of this pathway’s involvement in neuroblastoma. We discuss the RAS-MAPK pathway’s general functioning, the clinical implications of its deregulation in neuroblastoma, and current promising therapeutics targeting proteins involved in signaling.

TSCCA: A tensor sparse CCA method for detecting microRNA-gene patterns from multiple cancers

Existing studies have demonstrated that dysregulation of microRNAs (miRNAs or miRs) is involved in the initiation and progression of cancer. Many efforts have been devoted to identify microRNAs as potential biomarkers for cancer diagnosis, prognosis and therapeutic targets. With the rapid development of miRNA sequencing technology, a vast amount of miRNA expression data for multiple cancers has been collected. These invaluable data repositories provide new paradigms to explore the relationship between miRNAs and cancer. Thus, there is an urgent need to explore the complex cancer-related miRNA-gene patterns by integrating multi-omics data in a pan-cancer paradigm. In this study, we present a tensor sparse canonical correlation analysis (TSCCA) method for identifying cancer-related miRNA-gene modules across multiple cancers. TSCCA is able to overcome the drawbacks of existing solutions and capture both the cancer-shared and specific miRNA-gene co-expressed modules with better biological interpretations. We comprehensively evaluate the performance of TSCCA using a set of simulated data and matched miRNA/gene expression data across 33 cancer types from the TCGA database. We uncover several dysfunctional miRNA-gene modules with important biological functions and statistical significance. These modules can advance our understanding of miRNA regulatory mechanisms of cancer and provide insights into miRNA-based treatments for cancer.

MicroRNAs (miRNAs) are a class of small non-coding RNAs. Previous studies have revealed that miRNA-gene regulatory modules play key roles in the occurrence and development of cancer. However, little has been done to discover miRNA-gene regulatory modules from a pan-cancer view. Thus, it is urgently needed to develop new methods to explore the complex cancer-related miRNA-gene patterns by integrating multi-omics data of multi-cancers. To build the connections between miRNA-gene regulatory modules across different cancer types, we propose a tensor sparse canonical correlation analysis (TSCCA) method. Our specific contributions are two-fold: (1) We propose a sparse statistical learning model TSCCA and an efficient block-coordinate descent algorithm to solve it. (2) We apply TSCCA to a multi-omics data set of 33 cancer types from TCGA and identify some cancer-related miRNA-gene modules with important biological functions and statistical significance.

Genomic characterization of relapsed acute myeloid leukemia reveals novel putative therapeutic targets

Relapse is the leading cause of death of adult and pediatric patients with acute myeloid leukemia (AML). Numerous studies have helped to elucidate the complex mutational landscape at diagnosis of AML, leading to improved risk stratification and new therapeutic options. However, multi-whole-genome studies of adult and pediatric AML at relapse are necessary for further advances. To this end, we performed whole-genome and whole-exome sequencing analyses of longitudinal diagnosis, relapse, and/or primary resistant specimens from 48 adult and 25 pediatric patients with AML. We identified mutations recurrently gained at relapse in ARID1A and CSF1R, both of which represent potentially actionable therapeutic alternatives. Further, we report specific differences in the mutational spectrum between adult vs pediatric relapsed AML, with MGA and H3F3A p.Lys28Met mutations recurrently found at relapse in adults, whereas internal tandem duplications in UBTF were identified solely in children. Finally, our study revealed recurrent mutations in IKZF1, KANSL1, and NIPBL at relapse. All of the mentioned genes have either never been reported at diagnosis in de novo AML or have been reported at low frequency, suggesting important roles for these alterations predominantly in disease progression and/or resistance to therapy. Our findings shed further light on the complexity of relapsed AML and identified previously unappreciated alterations that may lead to improved outcomes through personalized medicine.

Multimodal analysis of cell-free DNA whole-genome sequencing for pediatric cancers with low mutational burden

Sequencing of cell-free DNA in the blood of cancer patients (liquid biopsy) provides attractive opportunities for early diagnosis, assessment of treatment response, and minimally invasive disease monitoring. To unlock liquid biopsy analysis for pediatric tumors with few genetic aberrations, we introduce an integrated genetic/epigenetic analysis method and demonstrate its utility on 241 deep whole-genome sequencing profiles of 95 patients with Ewing sarcoma and 31 patients with other pediatric sarcomas. Our method achieves sensitive detection and classification of circulating tumor DNA in peripheral blood independent of any genetic alterations. Moreover, we benchmark different metrics for cell-free DNA fragmentation analysis, and we introduce the LIQUORICE algorithm for detecting circulating tumor DNA based on cancer-specific chromatin signatures. Finally, we combine several fragmentation-based metrics into an integrated machine learning classifier for liquid biopsy analysis that exploits widespread epigenetic deregulation and is tailored to cancers with low mutation rates. Clinical associations highlight the potential value of cfDNA fragmentation patterns as prognostic biomarkers in Ewing sarcoma. In summary, our study provides a comprehensive analysis of circulating tumor DNA beyond recurrent genetic aberrations, and it renders the benefits of liquid biopsy more readily accessible for childhood 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.

Mutation accumulation in cancer genes relates to nonoptimal outcome in chronic myeloid leukemia

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm accounting for ∼15% of all leukemia. Progress of the disease from an indolent chronic phase to the more aggressive accelerated phase or blast phase (BP) occurs in a minority of cases and is associated with an accumulation of somatic mutations. We performed genetic profiling of 85 samples and transcriptome profiling of 12 samples from 59 CML patients. We identified recurrent somatic mutations in ABL1 (37%), ASXL1 (26%), RUNX1 (16%), and BCOR (16%) in the BP and observed that mutation signatures in the BP resembled those of acute myeloid leukemia (AML). We found that mutation load differed between the indolent and aggressive phases and that nonoptimal responders had more nonsilent mutations than did optimal responders at the time of diagnosis, as well as in follow-up. Using RNA sequencing, we identified other than BCR-ABL1 cancer-associated hybrid genes in 6 of the 7 BP samples. Uncovered expression alterations were in turn associated with mechanisms and pathways that could be targeted in CML management and by which somatic alterations may emerge in CML. Last, we showed the value of genetic data in CML management in a personalized medicine setting.

Immune subgroup analysis for non-small cell lung cancer may be a good choice for evaluating therapeutic efficacy and prognosis

Due to its effectiveness, cancer immunotherapy has attracted widespread attention from clinicians and scientific researchers. Numerous studies have proven that effective stratification of cancer patients would promote the personalized application of immunotherapy. Therefore, we used the transcriptome data of nearly 1,000 patients with non-small cell lung cancer (NSCLC) to construct a new immune subgroup. We found that the new immune subgroup, named cluster 2, was a mixture of lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC), and showed poor overall survival, which was further verified in the independent validation set. Immune infiltration correlation analysis showed that the Mast cell type and its status subdivisions had a predictive effect on the prognosis of NSCLC, especially in LUAD. Phenotypic analysis suggested that epithelial-mesenchymal transition (EMT) was positively correlated with immunosuppression, supporting the correlation between tumor phenotype and immune background. Although immune subtypes failed to significantly distinguish the progression-free survival (PFS) of immunotherapy patients, they showed the expected trend; the sample size needs to be further expanded for verification. In addition, some results indicated that the two cancer types, LUAD and LUSC, might require independent analyses.

Molecular evolution and the decline of purifying selection with age

Life history theory predicts that the intensity of selection declines with age, and this trend should impact how genes expressed at different ages evolve. Here we find consistent relationships between a gene’s age of expression and patterns of molecular evolution in two mammals (the human Homo sapiens and the mouse Mus musculus) and two insects (the malaria mosquito Anopheles gambiae and the fruit fly Drosophila melanogaster). When expressed later in life, genes fix nonsynonymous mutations more frequently, are more polymorphic for nonsynonymous mutations, and have shorter evolutionary lifespans, relative to those expressed early. The latter pattern is explained by a simple evolutionary model. Further, early-expressed genes tend to be enriched in similar gene ontology terms across species, while late-expressed genes show no such consistency. In humans, late-expressed genes are more likely to be linked to cancer and to segregate for dominant disease-causing mutations. Last, the effective strength of selection (N_es) decreases and the fraction of beneficial mutations increases with a gene’s age of expression. These results are consistent with the diminishing efficacy of purifying selection with age, as proposed by Medawar’s classic hypothesis for the evolution of senescence, and provide links between life history theory and molecular evolution.

A fundamental principle of evolutionary theory is that the force of natural selection is weaker on traits expressed late in life relative to traits expressed early. Here, the authors find strong and consistent patterns of molecular evolution reflecting this principle in four species of animals, including humans.

Detection of new drivers of frequent B-cell lymphoid neoplasms using an integrated analysis of whole genomes

B-cell lymphoproliferative disorders exhibit a diverse spectrum of diagnostic entities with heterogeneous behaviour. Multiple efforts have focused on the determination of the genomic drivers of B-cell lymphoma subtypes. In the meantime, the aggregation of diverse tumors in pan-cancer genomic studies has become a useful tool to detect new driver genes, while enabling the comparison of mutational patterns across tumors. Here we present an integrated analysis of 354 B-cell lymphoid disorders. 112 recurrently mutated genes were discovered, of which KMT2D, CREBBP, IGLL5 and BCL2 were the most frequent, and 31 genes were putative new drivers. Mutations in CREBBP, TNFRSF14 and KMT2D predominated in follicular lymphoma, whereas those in BTG2, HTA-A and PIM1 were more frequent in diffuse large B-cell lymphoma. Additionally, we discovered 31 significantly mutated protein networks, reinforcing the role of genes such as CREBBP, EEF1A1, STAT6, GNA13 and TP53, but also pointing towards a myriad of infrequent players in lymphomagenesis. Finally, we report aberrant expression of oncogenes and tumor suppressors associated with novel noncoding mutations (DTX1 and S1PR2), and new recurrent copy number aberrations affecting immune check-point regulators (CD83, PVR) and B-cell specific genes (TNFRSF13C). Our analysis expands the number of mutational drivers of B-cell lymphoid neoplasms, and identifies several differential somatic events between disease subtypes.

What Animal Cancers teach us about Human Biology

Cancers in animals present a large, underutilized reservoir of biomedical information with critical implication for human oncology and medicine in general. Discussing two distinct areas of tumour biology in non-human hosts, we highlight the importance of these findings for our current understanding of cancer, before proposing a coordinated strategy to harvest biomedical information from non-human resources and translate it into a clinical setting. First, infectious cancers that can be transmitted as allografts between individual hosts, have been identified in four distinct, unrelated groups, dogs, Tasmanian devils, Syrian hamsters and, surprisingly, marine bivalves. These malignancies might hold the key to improving our understanding of the interaction between tumour cell and immune system and, thus, allow us to devise novel treatment strategies that enhance anti-cancer immunosurveillance, as well as suggesting more effective organ and stem cell transplantation strategies. The existence of these malignancies also highlights the need for increased scrutiny when considering the existence of infectious cancers in humans. Second, it has long been understood that no linear relationship exists between the number of cells within an organism and the cancer incidence rate. To resolve what is known as Peto's Paradox, additional anticancer strategies within different species have to be postulated. These naturally occurring idiosyncrasies to avoid carcinogenesis represent novel potential therapeutic strategies.

DNA methylation as a key epigenetic player for hepatoblastoma characterization

BACKGROUND

Hepatoblastoma (HB) is a rare embryonal liver tumor of children. Although intrinsic biological differences between tumors can affect prognosis, few groups have studied these differences. Given the recent increased attention to epigenetic mechanisms in the genesis and progression of these tumors, we aimed to classify HB samples according to the stages of liver development and DNA methylation machinery.

BASIC PROCEDURES

We evaluated the expression of 24 genes associated with DNA methylation and stages of hepatocyte differentiation and global DNA methylation. Using bioinformatics tools and expression data, we propose a stratification model for HB.

MAIN FINDINGS

Tumors clustered into three groups that presented specific gene expression profiles of the panel of DNA methylation enzymes and hepatocyte differentiation markers. In addition to reinforcing these embryonal tumors' molecular heterogeneity, we propose that a panel of 13 genes can stratify HBs (TET1, TET2, TET3, DNMT1, DNMT3A, UHRF1, ALB, CYP3A4, TDO2, UGT1A1, AFP, HNF4A, and FOXA2). DNA methylation machinery participates in the characterization of HBs, directly reflected in diverse DNA methylation content. The data suggested that a subset of HBs were similar to differentiated livers, with upregulation of mature hepatocyte markers, decreased expression of DNA methylation enzymes, and higher global methylation levels; these findings might predict worse outcomes.

CONCLUSIONS

HBs are heterogeneous tumors. Despite using a small cohort of 21 HB samples, our findings reinforce that DNA methylation is a robust biomarker for this tumor type.

Integrative genomic analysis of pediatric T-cell lymphoblastic lymphoma reveals candidates of clinical significance

T-cell lymphoblastic lymphoma (T-LBL) is a heterogeneous malignancy of lymphoblasts committed to T-cell lineage. The dismal outcomes (15%-30%) after T-LBL relapse warrant establishing risk-based treatment. To our knowledge, this study presents the first comprehensive, systematic, integrated, genome-wide analysis including relapsed cases that identifies molecular markers of prognostic relevance for T-LBL. NOTCH1 was identified as the putative driver for T-LBL. An activated NOTCH/PI3K-AKT signaling axis and alterations in cell cycle regulators constitute the core oncogenic program for T-LBL. Mutated KMT2D was identified as a prognostic marker. The cumulative incidence of relapse was 47% ± 17% in patients with KMT2D mutations, compared with 14% ± 3% in wild-type KMT2D. Structural analysis of the mutated domains of KMT2D revealed a plausible impact on structure and functional consequences. These findings provide new insights into the pathogenesis of T-LBL, including high translational potential. The ongoing LBL 2018 trial (www.clinicaltrials.gov #NCT04043494) allows for prospective validation and subsequent fine tuning of the stratification criteria for T-LBL risk groups to improve survival of pediatric patients.

Integrated genomic analysis identifies driver genes and cisplatin-resistant progenitor phenotype in pediatric liver cancer

Pediatric liver cancers (PLCs) comprise diverse diseases affecting infants, children and adolescents. Despite overall good prognosis, PLCs display heterogeneous response to chemotherapy. Integrated genomic analysis of 126 pediatric liver tumors showed a continuum of driver mechanisms associated with patient age, including new targetable oncogenes. In 10% of hepatoblastoma patients, all before 3 years old, we identified a mosaic premalignant clonal expansion of cells altered at the 11p15.5 locus. Analysis of spatial and longitudinal heterogeneity revealed an important plasticity between 'Hepatocytic', 'Liver Progenitor' and 'Mesenchymal' molecular subgroups of hepatoblastoma. We showed that during chemotherapy, 'Liver Progenitor' cells accumulated massive loads of cisplatin-induced mutations with a specific mutational signature, leading to the development of heavily mutated relapses and metastases. Drug screening in PLC cell lines identified promising targets for cisplatin-resistant progenitor cells, validated in mouse xenograft experiments. These data provide new insights into cisplatin resistance mechanisms in PLC and suggest alternative therapies.

Unraveling Gene Fusions for Drug Repositioning in High-Risk Neuroblastoma

High-risk neuroblastoma (NB) remains a significant therapeutic challenge facing current pediatric oncology patients. Structural variants such as gene fusions have shown an initial promise in enhancing mechanistic understanding of NB and improving survival rates. In this study, we performed a comprehensive in silico investigation on the translational ability of gene fusions for patient stratification and treatment development for high-risk NB patients. Specifically, three state-of-the-art gene fusion detection algorithms, including ChimeraScan, SOAPfuse, and TopHat-Fusion, were employed to identify the fusion transcripts in a RNA-seq data set of 498 neuroblastoma patients. Then, the 176 high-risk patients were further stratified into four different subgroups based on gene fusion profiles. Furthermore, Kaplan-Meier survival analysis was performed, and differentially expressed genes (DEGs) for the redefined high-risk group were extracted and functionally analyzed. Finally, repositioning candidates were enriched in each patient subgroup with drug transcriptomic profiles from the LINCS L1000 Connectivity Map. We found the number of identified gene fusions was increased from clinical the low-risk stage to the high-risk stage. Although the technical concordance of fusion detection algorithms was suboptimal, they have a similar biological relevance concerning perturbed pathways and regulated DEGs. The gene fusion profiles could be utilized to redefine high-risk patient subgroups with significant onset age of NB, which yielded the improved survival curves (Log-rank p value ≤ 0.05). Out of 48 enriched repositioning candidates, 45 (93.8%) have antitumor potency, and 24 (50%) were confirmed with either on-going clinical trials or literature reports. The gene fusion profiles have a discrimination power for redefining patient subgroups in high-risk NB and facilitate precision medicine-based drug repositioning implementation.

Pan-Cancer Analysis of FURIN as a Potential Prognostic and Immunological Biomarker

Background

Furin is a calcium-dependent protease that processes various precursor proteins through diverse secretory pathways. The deregulation of FURIN correlated with the prognosis of patients in numerous diseases. However, the role of FURIN in human pan-cancer is still largely unknown.

Methods

Multiple bioinformatic methods were employed to comprehensively analyze the correlation of FURIN expression with prognosis, mismatch repair (MMR), microsatellite instability (MSI), tumor mutational burden (TMB), DNA methylation, tumor immune infiltration, and common immune checkpoint inhibitors (ICIs) from the public database, and aim to evaluate the potential prognostic value of FURIN across cancers.

Results

FURIN was aberrantly expressed and was strongly correlated with MMR, MSI, TMB, and DNA methylation in multiple types of cancer. Moreover, survival analysis across cancers revealed that FURIN expression was correlated with overall survival (OS) in four cancers, disease-specific survival (DSS) in five cancers, progression-free interval (PFI) in seven cancers, and disease-free interval (DFI) in two cancers. Also, FURIN expression was related to immune cell infiltration in 6 cancers and ImmuneScore/StromalScore in 10 cancers, respectively. In addition, FURIN expression also showed strong association between expression levels and immune checkpoint markers in three cancers.

Conclusion

FURIN can serve as a significant prognostic biomarker and correlate with tumor immunity in human pan-cancer.

Reconciling the Biological and Transcriptional Variability of Hepatoblastoma with Its Mutational Uniformity

Simple Summary

Hepatoblastoma (HB), the most common form of childhood liver cancer, is associated with dual mutation and/or dysregulation of the Wnt/β-catenin and Hippo pathways in ~50% of cases. However, this mutational simplicity cannot explain HB’s biological and histologic diversity. This discussion focuses upon recent work showing that specific β-catenin mutants are key determinants of this HB variability as well as their metabolic and transcriptional signatures. Dysregulation of the anti-oxidant NFE2L2 pathway also contributes to tumorigenesis by being directly transforming in association with either of the other two factors. The transcriptional overlap of tumors generated by pairs of factors identifies crucial targets that likely mediate HB tumorigenesis, behavior and appearance.

Abstract

Hepatoblastoma (HB), the most common childhood liver cancer, is associated with seven distinct histologic subtypes and variable degrees of clinical aggressiveness and presentation. Yet it is among the least genomically altered tumors known, with about half of HBs showing mutation and/or dysregulation of the Wnt/β-catenin and Hippo pathways. This raises the question of how this mutational simplicity can generate such biological and histologic complexity. Recent work shows that the identity of the underlying β-catenin mutation is a major contributor. Mutation or over-expression of the NFE2L2/NRF2 transcription factor, previously thought only to promote anti-oxidant responses, has also recently been shown to accelerate the growth of HBs generated by mutations in the Wnt/β-catenin and Hippo pathways while imparting novel features such as the tumor-associated cysts and necrosis. Moreover, patient-associated NFE2L2 mutations are overtly transforming when co-expressed with either mutant β-catenin or a Hippo pathway effector. The finding that tumorigenesis can be driven by any two arms of the β-catenin/Hippo/NFE2L2 axis has permitted the identification of a small subset of coordinately regulated tumor-specific transcripts, some of whose levels correlate with inferior long-term outcomes in HB and other cancers. Collectively, these findings begin to provide for more refined and molecularly based classification, survival algorithms and design of chemotherapeutic regimens.