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

Molecular profiling of 888 pediatric tumors informs future precision trials and data-sharing initiatives in pediatric cancer

To inform clinical trial design and real-world precision pediatric oncology practice, we classified diagnoses, assessed the landscape of mutations, and identified genomic variants matching trials in a large unselected institutional cohort of solid tumors patients sequenced at Dana-Farber / Boston Children's Cancer and Blood Disorders Center. Tumors were sequenced with OncoPanel, a targeted next-generation DNA sequencing panel. Diagnoses were classified according to the International Classification of Diseases for Oncology (ICD-O-3.2). Over 6.5 years, 888 pediatric cancer patients with 95 distinct diagnoses had successful tumor sequencing. Overall, 33% (n = 289/888) of patients had at least 1 variant matching a precision oncology trial protocol, and 14% (41/289) were treated with molecularly targeted therapy. This study highlights opportunities to use genomic data from hospital-based sequencing performed either for research or clinical care to inform ongoing and future precision oncology clinical trials. Furthermore, the study results emphasize the importance of data sharing to define the genomic landscape and targeted treatment opportunities for the large group of rare pediatric cancers we encounter in clinical practice.

Benefits for children with suspected cancer from routine whole-genome sequencing

Clinical whole-genome sequencing (WGS) has been shown to deliver potential benefits to children with cancer and to alter treatment in high-risk patient groups. It remains unknown whether offering WGS to every child with suspected cancer can change patient management. We collected WGS variant calls and clinical and diagnostic information from 281 children (282 tumors) across two English units (n = 152 from a hematology center, n = 130 from a solid tumor center) where WGS had become a routine test. Our key finding was that variants uniquely attributable to WGS changed the management in ~7% (20 out of 282) of cases while providing additional disease-relevant findings, beyond standard-of-care molecular tests, in 108 instances for 83 (29%) cases. Furthermore, WGS faithfully reproduced every standard-of-care molecular test (n = 738) and revealed several previously unknown genomic features of childhood tumors. We show that WGS can be delivered as part of routine clinical care to children with suspected cancer and can change clinical management by delivering unexpected genomic insights. Our experience portrays WGS as a clinically impactful assay for routine practice, providing opportunities for assay consolidation and for delivery of molecularly informed patient care.

Recent Advancements and Innovations in Pediatric Precision Oncology

Precision oncology incorporates comprehensive genomic profiling into the individualized clinical care of pediatric cancer patients. In recent years, comprehensive pan-cancer analyses have led to the successful implementation of genomics-based pediatric trials and accelerated approval of novel targeted agents. In addition, disease-specific studies have resulted in molecular subclassification of myriad cancer types with subsequent tailoring of treatment intensity based on the patient's prognostic factors. This review discusses the progress of the field and highlights developments that are leading to more personalized cancer care and improved patient outcomes. Increased understanding of the evolution of precision oncology over recent decades emphasizes the tremendous impact of improved genomic applications. New technologies and improved diagnostic modalities offer further promise for future advancements within the field.

Multi-scale variational autoencoder for imputation of missing values in untargeted metabolomics using whole-genome sequencing data

BACKGROUND

Missing data is a common challenge in mass spectrometry-based metabolomics, which can lead to biased and incomplete analyses. The integration of whole-genome sequencing (WGS) data with metabolomics data has emerged as a promising approach to enhance the accuracy of data imputation in metabolomics studies.

METHOD

In this study, we propose a novel method that leverages the information from WGS data and reference metabolites to impute unknown metabolites. Our approach utilizes a multi-scale variational autoencoder to jointly model the burden score, polygenetic risk score (PGS), and linkage disequilibrium (LD) pruned single nucleotide polymorphisms (SNPs) for feature extraction and missing metabolomics data imputation. By learning the latent representations of both omics data, our method can effectively impute missing metabolomics values based on genomic information.

RESULTS

We evaluate the performance of our method on empirical metabolomics datasets with missing values and demonstrate its superiority compared to conventional imputation techniques. Using 35 template metabolites derived burden scores, PGS and LD-pruned SNPs, the proposed methods achieved R2-scores > 0.01 for 71.55 % of metabolites.

CONCLUSION

The integration of WGS data in metabolomics imputation not only improves data completeness but also enhances downstream analyses, paving the way for more comprehensive and accurate investigations of metabolic pathways and disease associations. Our findings offer valuable insights into the potential benefits of utilizing WGS data for metabolomics data imputation and underscore the importance of leveraging multi-modal data integration in precision medicine research.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs

BACKGROUND

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models.

RESULTS

Here, we generate single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We develop an unsupervised machine learning approach ("automatic consensus nonnegative matrix factorization" (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirm a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly, however, this weak-mesenchymal-like program is maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 h, suggesting an uncharacterized therapy-escape mechanism.

CONCLUSIONS

Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

Survival trend and outcome prediction for pediatric Hodgkin and non-Hodgkin lymphomas based on machine learning

Pediatric Hodgkin and non-Hodgkin lymphomas differ from adult cases in biology and management, yet there is a lack of survival analysis tailored to pediatric lymphoma. We analyzed lymphoma data from 1975 to 2018, comparing survival trends between 7,871 pediatric and 226,211 adult patients, identified key risk factors for pediatric lymphoma survival, developed a predictive nomogram, and utilized machine learning to predict long-term lymphoma-specific mortality risk. Between 1975 and 2018, we observed substantial increases in 1-year (19.3%), 5-year (41.9%), and 10-year (48.8%) overall survival rates in pediatric patients with lymphoma. Prognostic factors such as age, sex, race, Ann Arbor stage, lymphoma subtypes, and radiotherapy were incorporated into the nomogram. The nomogram exhibited excellent predictive performance with area under the curve (AUC) values of 0.766, 0.724, and 0.703 for one-year, five-year, and ten-year survival, respectively, in the training cohort, and AUC values of 0.776, 0.712, and 0.696 in the validation cohort. Importantly, the nomogram outperformed the Ann Arbor staging system in survival prediction. Machine learning models achieved AUC values of approximately 0.75, surpassing the conventional method (AUC =  ~ 0.70) in predicting the risk of lymphoma-specific death. We also observed that pediatric lymphoma survivors had a substantially reduced risk of lymphoma after ten years b,ut faced an increasing risk of non-lymphoma diseases. The study highlights substantial improvements in pediatric lymphoma survival, offers reliable predictive tools, and underscores the importance of long-term monitoring for non-lymphoma health issues in pediatric patients.

Immune Microenvironment in Childhood Cancers: Characteristics and Therapeutic Challenges

The tumor immune microenvironment is pivotal in cancer initiation, advancement, and regulation. Its molecular and cellular composition is critical throughout the disease, as it can influence the balance between suppressive and cytotoxic immune responses within the tumor's vicinity. Studies on the tumor immune microenvironment have enriched our understanding of the intricate interplay between tumors and their immunological surroundings in various human cancers. These studies illuminate the role of significant components of the immune microenvironment, which have not been extensively explored in pediatric tumors before and may influence the responsiveness or resistance to therapeutic agents. Our deepening understanding of the pediatric tumor immune microenvironment is helping to overcome challenges related to the effectiveness of existing therapeutic strategies, including immunotherapies. Although in the early stages, targeted therapies that modulate the tumor immune microenvironment of pediatric solid tumors hold promise for improved outcomes. Focusing on various aspects of tumor immune biology in pediatric patients presents a therapeutic opportunity that could improve treatment outcomes. This review offers a comprehensive examination of recent literature concerning profiling the immune microenvironment in various pediatric tumors. It seeks to condense research findings on characterizing the immune microenvironment in pediatric tumors and its impact on tumor development, metastasis, and response to therapeutic modalities. It covers the immune microenvironment's role in tumor development, interactions with tumor cells, and its impact on the tumor's response to immunotherapy. The review also discusses challenges targeting the immune microenvironment for pediatric cancer therapies.

Comprehensive single cell transcriptomics analysis of murine osteosarcoma uncovers Skp2 function in metastasis, genomic instability and immune activation and reveals additional target pathways

Osteosarcoma (OS) is the most common primary pediatric bone malignancy. One promising new therapeutic target is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase responsible for ubiquitination and proteasome degradation of substrate p27, thus driving cellular proliferation. We have shown previously that knockout of Skp2 in an immunocompetent transgenic mouse model of OS improved survival, drove apoptosis, and induced tumor inflammation. Here, we applied single-cell RNA-sequencing (scRNA-seq) to study primary OS tumors derived from Osx-Cre driven conditional knockout of Rb1 and Trp53. We showed that murine OS models recapitulate the tumor heterogeneity and microenvironment complexity observed in patient tumors. We further compared this model with OS models with functional disruption of Skp2: one with Skp2 knockout and the other with the Skp2-p27 interaction disrupted (resulting in p27 overexpression). We found reduction of T cell exhaustion and upregulation of interferon activation, along with evidence of replicative and endoplasmic reticulum-related stress in the Skp2 disruption models, and showed that interferon induction was correlated with improved survival in OS patients. Additionally, our scRNA-seq analysis uncovered decreased activities of metastasis-related gene signatures in the Skp2-disrupted OS, which we validated by observation of a strong reduction in lung metastasis in the Skp2 knockout mice. Finally, we report several potential mechanisms of escape from targeting Skp2 in OS, including upregulation of Myc targets, DNA copy number amplification and overexpression of alternative E3 ligase genes, and potential alternative lineage activation. These mechanistic insights into OS tumor biology and Skp2 function suggest novel targets for new, synergistic therapies, while the data and our comprehensive analysis may serve as a public resource for further big data-driven OS research.

Clonal hematopoiesis of indeterminate potential is rare in pediatric patients undergoing autologous stem cell transplantation

Clonal hematopoiesis of indeterminate potential (CHIP) describes recurrent somatic gene mutations in the blood of healthy individuals, associated with higher risk for hematological malignancies and higher all-cause mortality by cardiovascular disease. CHIP increases with age and is more common in adult patients after chemotherapy or radiation for cancer. Furthermore, in some adult patients undergoing autologous stem cell transplantation (ASCT) or thereafter, CHIP has been identified. In children and adolescents, it remains unclear how cellular stressors such as cytotoxic therapy influence the incidence and expansion of CHIP. We conducted a retrospective study on 33 pediatric patients mostly with solid tumors undergoing ASCT for presence of CHIP. We analyzed CD34+ selected peripheral blood stem cell grafts after several cycles of chemotherapy, prior to cell infusion, by next-generation sequencing including 18 "CHIP-genes". Apart from a somatic variant in TP53 in one patient no other variants indicative of CHIP were identified. As a CHIP-unrelated finding, germline variants in CHEK2 and in ATM were identified in two and four patients, respectively. In conclusion, we could not detect "typical" CHIP variants in our cohort of pediatric cancer patients undergoing ASCT. However, more studies with larger patient numbers are necessary to assess if chemotherapy in the pediatric setting contributes to an increased CHIP incidence and at what time point.

Promise and Perils of Precision Oncology for Patients With Pediatric and Young Adult Sarcomas

The completion of multiple national pediatric precision oncology platform trials and the incorporation of standardized molecular profiling into the diagnostic care of pediatric and young adult patients with sarcomas have proven the feasibility and potential of the approach. In this work, we explore the current state of the art of precision oncology for pediatric and young adults with sarcoma. We highlight important lessons learned and the challenges that should be addressed in the next generation of trials. The chapter outlines current efforts to improve standardization of molecular assays, harmonization of data collection, and novel molecular tools such as cell-free DNA analyses. Finally, we discuss the impacts and psychosocial outcomes experienced by patients and communication strategies for providers.

Developmental origins shape the paediatric cancer genome

In the past two decades, technological advances have brought unprecedented insights into the paediatric cancer genome revealing characteristics distinct from those of adult cancer. Originating from developing tissues, paediatric cancers generally have low mutation burden and are driven by variants that disrupt the transcriptional activity, chromatin state, non-coding cis-regulatory regions and other biological functions. Within each tumour, there are multiple populations of cells with varying states, and the lineages of some can be tracked to their fetal origins. Genome-wide genetic screening has identified vulnerabilities associated with both the cell of origin and transcription deregulation in paediatric cancer, which have become a valuable resource for designing new therapeutic approaches including those for small molecules, immunotherapy and targeted protein degradation. In this Review, we present recent findings on these facets of paediatric cancer from a pan-cancer perspective and provide an outlook on future investigations.

Proteomic insights into paediatric cancer: Unravelling molecular signatures and therapeutic opportunities

Survival rates in some paediatric cancers have improved greatly over recent decades, in part due to the identification of diagnostic, prognostic and predictive molecular signatures, and the development of risk-directed therapies. However, other paediatric cancers have proved difficult to treat, and there is an urgent need to identify novel biomarkers that reveal therapeutic opportunities. The proteome is the total set of expressed proteins present in a cell or tissue at a point in time, and is vastly more dynamic than the genome. Proteomics holds significant promise for cancer research, as proteins are ultimately responsible for cellular phenotype and are the target of most anticancer drugs. Here, we review the discoveries, opportunities and challenges of proteomic analyses in paediatric cancer, with a focus on mass spectrometry (MS)-based approaches. Accelerating incorporation of proteomics into paediatric precision medicine has the potential to improve survival and quality of life for children with cancer.

Genomic Newborn Screening for Pediatric Cancer Predisposition Syndromes: A Holistic Approach

As next-generation sequencing (NGS) has become more widely used, germline and rare genetic variations responsible for inherited illnesses, including cancer predisposition syndromes (CPSs) that account for up to 10% of childhood malignancies, have been found. The CPSs are a group of germline genetic disorders that have been identified as risk factors for pediatric cancer development. Excluding a few "classic" CPSs, there is no agreement regarding when and how to conduct germline genetic diagnostic studies in children with cancer due to the constant evolution of knowledge in NGS technologies. Various clinical screening tools have been suggested to aid in the identification of individuals who are at greater risk, using diverse strategies and with varied outcomes. We present here an overview of the primary clinical and molecular characteristics of various CPSs and summarize the existing clinical genomics data on the prevalence of CPSs in pediatric cancer patients. Additionally, we discuss several ethical issues, challenges, limitations, cost-effectiveness, and integration of genomic newborn screening for CPSs into a healthcare system. Furthermore, we assess the effectiveness of commonly utilized decision-support tools in identifying patients who may benefit from genetic counseling and/or direct genetic testing. This investigation highlights a tailored and systematic approach utilizing medical newborn screening tools such as the genome sequencing of high-risk newborns for CPSs, which could be a practical and cost-effective strategy in pediatric cancer care.

Comparative analysis of integrative classification methods for multi-omics data

Recent advances in sequencing, mass spectrometry, and cytometry technologies have enabled researchers to collect multiple 'omics data types from a single sample. These large datasets have led to a growing consensus that a holistic approach is needed to identify new candidate biomarkers and unveil mechanisms underlying disease etiology, a key to precision medicine. While many reviews and benchmarks have been conducted on unsupervised approaches, their supervised counterparts have received less attention in the literature and no gold standard has emerged yet. In this work, we present a thorough comparison of a selection of six methods, representative of the main families of intermediate integrative approaches (matrix factorization, multiple kernel methods, ensemble learning, and graph-based methods). As non-integrative control, random forest was performed on concatenated and separated data types. Methods were evaluated for classification performance on both simulated and real-world datasets, the latter being carefully selected to cover different medical applications (infectious diseases, oncology, and vaccines) and data modalities. A total of 15 simulation scenarios were designed from the real-world datasets to explore a large and realistic parameter space (e.g. sample size, dimensionality, class imbalance, effect size). On real data, the method comparison showed that integrative approaches performed better or equally well than their non-integrative counterpart. By contrast, DIABLO and the four random forest alternatives outperform the others across the majority of simulation scenarios. The strengths and limitations of these methods are discussed in detail as well as guidelines for future applications.

The Neuroblastoma Microenvironment, Heterogeneity and Immunotherapeutic Approaches

Neuroblastoma is a peripheral nervous system tumor that almost exclusively occurs in young children. Although intensified treatment modalities have led to increased patient survival, the prognosis for patients with high-risk disease is still around 50%, signifying neuroblastoma as a leading cause of cancer-related deaths in children. Neuroblastoma is an embryonal tumor and is shaped by its origin from cells within the neural crest. Hence, neuroblastoma usually presents with a low mutational burden and is, in the majority of cases, driven by epigenetically deregulated transcription networks. The recent development of Omic techniques has given us detailed knowledge of neuroblastoma evolution, heterogeneity, and plasticity, as well as intra- and intercellular molecular communication networks within the neuroblastoma microenvironment. Here, we discuss the potential of these recent discoveries with emphasis on new treatment modalities, including immunotherapies which hold promise for better future treatment regimens.

Evaluation of Hi-C sequencing for the detection of gene fusions in hematologic and solid pediatric cancer samples

HiC sequencing is a DNA-based next-generation sequencing method that preserves the 3D conformation of the genome and has shown promise in detecting genomic rearrangements in translational research studies. To evaluate HiC as a potential clinical diagnostic platform, analytical concordance with routine laboratory testing was assessed using primary pediatric leukemia and sarcoma specimens previously positive for clinically significant genomic rearrangements. Archived specimen types tested included viable and nonviable frozen leukemic cells, as well as formalin-fixed paraffin-embedded (FFPE) tumor tissues. Initially, pediatric acute myeloid leukemia (AML) and alveolar rhabdomyosarcoma (A-RMS) specimens with known genomic rearrangements were subjected to HiC analysis to assess analytical concordance. Subsequently, a discovery cohort consisting of AML and acute lymphoblastic leukemia (ALL) cases with no known genomic rearrangements based on prior clinical diagnostic testing were evaluated to determine whether HiC could detect rearrangements. Using a standard sequencing depth of 50 million raw read-pairs per sample, or approximately 5X raw genomic coverage, 100% concordance was observed between HiC and previous clinical cytogenetic and molecular testing. In the discovery cohort, a clinically relevant gene fusion was detected in 45% of leukemia cases (5/11). This study demonstrates the value of HiC sequencing to medical diagnostic testing as it identified several clinically significant rearrangements, including those that might have been missed by current clinical testing workflows.

Key points

HiC sequencing is a DNA-based next-generation sequencing method that preserves the 3D conformation of the genome, facilitating detection of genomic rearrangements.HiC was 100% concordant with clinical diagnostic testing workflows for detecting clinically significant genomic rearrangements in pediatric leukemia and rhabdomyosarcoma specimens.HiC detected clinically significant genomic rearrangements not previously detected by prior clinical cytogenetic and molecular testing.HiC performed well with archived non-viable and viable frozen leukemic cell samples, as well as archived formalin-fixed paraffin-embedded tumor tissue specimens.

Single-Cell multiomics reveals ENL mutation perturbs kidney developmental trajectory by rewiring gene regulatory landscape

Cell differentiation during organogenesis relies on precise epigenetic and transcriptional control. Disruptions to this regulation can result in developmental abnormalities and malignancies, yet the underlying mechanisms are not well understood. Wilms tumors, a type of embryonal tumor closely linked to disrupted organogenesis, harbor mutations in epigenetic regulators in 30-50% of cases. However, the role of these regulators in kidney development and pathogenesis remains unexplored. By integrating mouse modeling, histological characterizations, and single-cell transcriptomics and chromatin accessibility profiling, we show that a Wilms tumor-associated mutation in the chromatin reader protein ENL disrupts kidney development trajectory by rewiring the gene regulatory landscape. Specifically, the mutant ENL promotes the commitment of nephron progenitors while simultaneously restricting their differentiation by dysregulating key transcription factor regulons, particularly the HOX clusters. It also induces the emergence of abnormal progenitor cells that lose their chromatin identity associated with kidney specification. Furthermore, the mutant ENL might modulate stroma-nephron interactions via paracrine Wnt signaling. These multifaceted effects caused by the mutation result in severe developmental defects in the kidney and early postnatal mortality in mice. Notably, transient inhibition of the histone acetylation binding activity of mutant ENL with a small molecule displaces transcriptional condensates formed by mutant ENL from target genes, abolishes its gene activation function, and restores developmental defects in mice. This work provides new insights into how mutations in epigenetic regulators can alter the gene regulatory landscape to disrupt kidney developmental programs at single-cell resolution in vivo . It also offers a proof-of-concept for the use of epigenetics-targeted agents to rectify developmental defects.

Gene Amplification of Mediator Subunit 30 Redirects the MYC Transcriptional Program and Oncogenesis

Understanding the molecular mechanisms underlying tumorigenesis is crucial for developing effective cancer therapies. Here, we investigate the co-amplification of MED30 and MYC across diverse cancer types and its impact on oncogenic transcriptional programs. Transcriptional profiling of MYC and MED30 single or both overexpression/amplification revealed the over amount of MED30 lead MYC to a new transcriptional program that associate with poor prognosis. Mechanistically, MED30 overexpression/amplification recruits other Mediator components and binding of MYC to a small subset of novel genomic regulatory sites, changing the epigenetic marks and inducing the formation of new enhancers, which drive the expression of target genes crucial for cancer progression. In vivo studies in pancreatic ductal adenocarcinoma (PDAC) further validate the oncogenic potential of MED30, as its overexpression promotes tumor growth and can be attenuated by knockdown of MYC. Using another cancer type as an example, MED30 knockdown reduces tumor growth particularly in MYC high-expressed glioblastoma (GBM) cell lines. Overall, our study elucidates the critical role of MED30 overexpression in orchestrating oncogenic transcriptional programs and highlights its potential as a therapeutic target for MYC-amplified cancer.

A systematic pan-cancer analysis identifies LDHA as a novel predictor for immunological, prognostic, and immunotherapy resistance

Lactate dehydrogenase A (LDHA), a critical enzyme involved in glycolysis, is broadly involved multiple biological functions in human cancers. It is reported that LDHA can impact tumor immune surveillance and induce the transformation of tumor-associated macrophages, highlighting its unnoticed function of LDHA in immune system. However, in human cancers, the role of LDHA in prognosis and immunotherapy hasn't been investigated. In this study, we analyzed the expression pattern and prognostic value of LDHA in pan-cancer and explored its association between tumor microenvironment (TME), immune infiltration subtype, stemness scores, tumor mutation burden (TMB), and immunotherapy resistance. We found that LDHA expression is tumor heterogeneous and that its high expression is associated with poor prognosis in multiple human cancers. In addition, LDHA expression was positively correlated with the presence of mononuclear/macrophage cells, and also promoted the infiltration of a range of immune cells. Genomic alteration of LDHA was common in different types of cancer, while with prognostic value in pan-cancers. Pan-cancer analysis revealed that the significant correlations existed between LDHA expression and tumor microenvironment (including stromal cells and immune cells) as well as stemness scores (DNAss and RNAss) across cancer types. Drug sensitivity analysis also revealed that LDHA was able to predict response to chemotherapy and immunotherapy. Furthermore, it was confirmed that knockdown of LDHA reduced proliferation and migration ability of lung cancer cells. Taken together, LDHA could serve as a prognostic biomarker and a potential immunotherapy marker.

A human neural crest model reveals the developmental impact of neuroblastoma-associated chromosomal aberrations

Early childhood tumours arise from transformed embryonic cells, which often carry large copy number alterations (CNA). However, it remains unclear how CNAs contribute to embryonic tumourigenesis due to a lack of suitable models. Here we employ female human embryonic stem cell (hESC) differentiation and single-cell transcriptome and epigenome analysis to assess the effects of chromosome 17q/1q gains, which are prevalent in the embryonal tumour neuroblastoma (NB). We show that CNAs impair the specification of trunk neural crest (NC) cells and their sympathoadrenal derivatives, the putative cells-of-origin of NB. This effect is exacerbated upon overexpression of MYCN, whose amplification co-occurs with CNAs in NB. Moreover, CNAs potentiate the pro-tumourigenic effects of MYCN and mutant NC cells resemble NB cells in tumours. These changes correlate with a stepwise aberration of developmental transcription factor networks. Together, our results sketch a mechanistic framework for the CNA-driven initiation of embryonal tumours.

The landscape of alterations affecting epigenetic regulators in T-cell acute lymphoblastic leukemia: Roles in leukemogenesis and therapeutic opportunities

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive malignancy accounting for 10%-15% of pediatric and 20%-25% of adult ALL cases. Epigenetic irregularities in T-ALL include alterations in both DNA methylation and the post-translational modifications on histones which together play a critical role in the initiation and development of T-ALL. Characterizing the oncogenic mutations that result in these epigenetic changes combined with the reversibility of epigenetic modifications represents an opportunity for the development of epigenetic therapies. Oncogenic mutations and deregulated expression of DNA methyltransferases (DNMTs), Ten-Eleven Translocation dioxygenases (TETs), Histone acetyltransferases (HATs) and members of Polycomb Repressor Complex 2 (PRC2) have all been identified in T-ALL. This review focuses on the current understanding of how these mutations lead to epigenetic changes in T-ALL, their association with disease pathogenesis and the current efforts to exploit these clinically through the development of epigenetic therapies in T-ALL treatment.

An imbalance between proliferation and differentiation underlies the development of microRNA-defective pineoblastoma

Mutations in the microRNA processing genes DICER1 and DROSHA drive several cancers that resemble embryonic progenitors. To understand how microRNAs regulate tumorigenesis, we ablated Drosha or Dicer1 in the developing pineal gland to emulate the pathogenesis of pineoblastoma, a brain tumor that resembles undifferentiated precursors of the pineal gland. Accordingly, these mice develop pineal tumors marked by loss of microRNAs, including the let-7/miR-98-5p family, and de-repression of microRNA target genes. Pineal tumors driven by loss of Drosha or Dicer1 mimic tumors driven by Rb1 loss, as they exhibit upregulation of S-phase genes and homeobox transcription factors that regulate pineal development. Blocking proliferation of these tumors facilitates expression of pinealocyte maturation markers, with a concomitant reduction in embryonic markers. Select embryonic markers remain elevated, however, as the microRNAs that normally repress these target genes remain absent. One such microRNA target gene is the oncofetal transcription factor Plagl2, which regulates expression of pro-growth genes, and inhibiting their signaling impairs tumor growth. Thus, we demonstrate that tumors driven by loss of microRNA processing may be therapeutically targeted by inhibiting downstream drivers of proliferation.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. Here, we generated single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We developed an unsupervised machine learning approach ('automatic consensus nonnegative matrix factorization' (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirmed a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly however, this weak-mesenchymal-like program was maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 hours, suggesting an uncharacterized therapy-escape mechanism. Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

Discovery of Novel Potential Prognostic Markers and Targeted Therapy to Overcome Chemotherapy Resistance in an Advanced-Stage Wilms Tumor

Wilms tumor (WT), the most prevalent type of renal cancer in children, exhibits overall survival rates exceeding 90%. However, chemotherapy resistance, which occurs in approximately 10% of WT cases, is a major challenge for the treatment of WT, particularly for advanced-stage patients. In this study, we aimed to discover potential mutation markers and drug targets associated with chemotherapy resistance in advanced-stage WT. We performed exome sequencing to detect somatic mutations and molecular targets in 43 WT samples, comprising 26 advanced-stage WTs, of which 7 cases were chemotherapy-resistant. Our analysis revealed four genes (ALPK2, C16orf96, PRKDC, and SVIL) that correlated with chemotherapy resistance and reduced disease-free survival in advanced-stage WT. Additionally, we identified driver mutations in 55 genes within the chemotherapy-resistant group, including 14 druggable cancer driver genes. Based on the mutation profiles of the resistant WT samples, we propose potential therapeutic strategies involving platinum-based agents, PARP inhibitors, and antibiotic/antineoplastic agents. Our findings provide insights into the genetic landscape of WT and offer potential avenues for targeted treatment, particularly for patients with chemotherapy resistance.

Clinical, pathologic and molecular findings in 2 Rottweiler littermates with appendicular osteosarcoma

Appendicular osteosarcoma was diagnosed and treated in a pair of littermate Rottweiler dogs, resulting in distinctly different clinical outcomes despite similar therapy within the context of a prospective, randomized clinical trial (NCI-COTC021/022). Histopathology, immunohistochemistry, mRNA sequencing, and targeted DNA hotspot sequencing techniques were applied to both dogs' tumors to define factors that could underpin their differential response to treatment. We describe the comparison of their clinical, histologic and molecular features, as well as those from a companion cohort of Rottweiler dogs, providing new insight into potential prognostic biomarkers for canine osteosarcoma.

Elimusertib has Antitumor Activity in Preclinical Patient-Derived Pediatric Solid Tumor Models

The small-molecule inhibitor of ataxia telangiectasia and Rad3-related protein (ATR), elimusertib, is currently being tested clinically in various cancer entities in adults and children. Its preclinical antitumor activity in pediatric malignancies, however, is largely unknown. We here assessed the preclinical activity of elimusertib in 38 cell lines and 32 patient-derived xenograft (PDX) models derived from common pediatric solid tumor entities. Detailed in vitro and in vivo molecular characterization of the treated models enabled the evaluation of response biomarkers. Pronounced objective response rates were observed for elimusertib monotherapy in PDX, when treated with a regimen currently used in clinical trials. Strikingly, elimusertib showed stronger antitumor effects than some standard-of-care chemotherapies, particularly in alveolar rhabdomysarcoma PDX. Thus, elimusertib has strong preclinical antitumor activity in pediatric solid tumor models, which may translate to clinically meaningful responses in patients.

Precision Oncology in Pediatric Cancer Surgery

Pediatric precision oncology has provided a greater understanding of the wide range of molecular alterations in difficult-to-treat or rare tumors with the aims of increasing survival as well as decreasing toxicity and morbidity from current cytotoxic therapies. In this article, the authors discuss the current state of pediatric precision oncology which has increased access to novel targeted therapies while also providing a framework for clinical implementation in this unique population. The authors evaluate the targetable mutations currently under investigation-with a focus on pediatric solid tumors-and discuss the key surgical implications associated with novel targeted therapies.

Pediatric Cancer Drug Development: Leveraging Insights in Cancer Biology and the Evolving Regulatory Landscape to Address Challenges and Guide Further Progress

The discovery and development of anticancer drugs for pediatric patients have historically languished when compared to both past and recent activity in drug development for adult patients, notably the dramatic spike of targeted and immune-oncology therapies. The reasons for this difference are multifactorial. Recent changes in the regulatory landscape surrounding pediatric cancer drug development and the understanding that some pediatric cancers are driven by genetic perturbations that also drive disparate adult cancers afford new opportunities. The unique cancer-initiating events and dependencies of many pediatric cancers, however, require additional pediatric-specific strategies. Research efforts to unravel the underlying biology of pediatric cancers, innovative clinical trial designs, model-informed drug development, extrapolation from adult data, addressing the unique considerations in pediatric patients, and use of pediatric appropriate formulations, should all be considered for efficient development and dosage optimization of anticancer drugs for pediatric patients.

Comprehensive analysis of ADGRE5 gene in human tumors: Clinical relevance, prognostic implications, and potential for personalized immunotherapy

Purpose

The Adhesion G protein receptor E5 (ADGRE5) gene is involved in a wide range of biological functions in human tumors; however, its specific molecular mechanism and significance in the analysis of human tumors have not yet been determined. Here, we provide a comprehensive genomic architecture of ADGRE5 in the tumor immune microenvironment and its clinical relevance across a broad range of solid tumors.

Methods

In this study, we used publicly available bioinformatics databases, with a primary focus on The Cancer Genome Atlas (TCGA) database and GTEx data, to conduct a comprehensive analysis of the impact on patient prognosis associated with ADGRE5.

Results

Statistics of more than 30 solid tumors from TCGA and Cancer Cell Line Encyclopedia (CCLE) were examined. ADGRE5 was differentially expressed in several cancers and was significantly associated with survival outcomes. Higher ADGRE5 levels were associated with worse prognosis in adrenocortical carcinoma, low grade glioma of the brain (LGG), lung squamous cell carcinoma, liver hepatocellular carcinoma, and uveal melanoma (UVM). Additionally, ADGRE5 was found to be an independent risk factor for LGG and UVM. The clinical relevance of ADGRE5 in tumor immunogenicity was further investigated. The expression level of ADGRE5 was not only strongly associated with tumor infiltration, such as tumor-infiltrating immune cells and immune subtypes, but also with tumor mutation burden, pyroptosis, and epithelial-mesenchymal transition in various types of cancer (P < 0.05). Furthermore, we noted that ADGRE5 exhibited a positive association with targeted drug sensitivity and conversely, a negative association with traditional chemotherapeutic drug sensitivity. Thus, ADGRE5 is expected to be a guiding marker gene for clinical prognosis and personalized tumor immunotherapy.

Rapid multiple protein sequence search by parallel and heterogeneous computation

MOTIVATION

Protein sequence database search and multiple sequence alignment generation is a fundamental task in many bioinformatics analyses. As the data volume of sequences continues to grow rapidly, there is an increasing need for efficient and scalable multiple sequence query algorithms for super-large databases without expensive time and computational costs.

RESULTS

We introduce Chorus, a novel protein sequence query system that leverages parallel model and heterogeneous computation architecture to enable users to query thousands of protein sequences concurrently against large protein databases on a desktop workstation. Chorus achieves over 100× speedup over BLASTP without sacrificing sensitivity. We demonstrate the utility of Chorus through a case study of analyzing a ∼1.5-TB large-scale metagenomic datasets for novel CRISPR-Cas protein discovery within 30 min.

AVAILABILITY AND IMPLEMENTATION

Chorus is open-source and its code repository is available at https://github.com/Bio-Acc/Chorus.

Integrative analysis revealed a correlation of PIAS family genes expression with prognosis, immunomodulation and chemotherapy

BACKGROUND

Protein inhibitor of activated STATs (PIAS) has pleiotropic biological effects, such as protein post-translational modification, transcriptional coregulation and gene editing. It is reported that PIAS family genes are also correlated with immune cells infiltration in cancers that highlights their unnoticed biological role in tumor progression. However, the relationship of their expression with prognosis, immune cell infiltration, tumor microenvironment, and immunotherapy in pan-cancer has been rarely reported.

METHODS

The multi-omics data were used to investigate the expression level of PIAS family members in pan-cancer, and the prognostic value of their expression in different tumors was analyzed by univariate Cox regression and Kaplan-Meier. Correlation analysis was used to investigate the relationship of PIAS gene expression with tumor microenvironment, immune infiltrating subtypes, stemness score and drug sensitivity. In addition, we also used wound healing and transwell assays to verify the biological effects of PIAS family gene expression on invasion and metastasis of HCC cells.

RESULTS

We found that PIAS family genes expression is significantly heterogeneous in tumors by multi-genomic analysis, and associated with poor prognosis in patients with multiple types of cancer. Furthermore, we also found that genetic alterations of PIAS family genes were not only common in different types of human tumors, but were also significantly associated with disease-free survival (DFS) across pan-cancer. Single-cell analysis revealed that PIAS family genes were mainly distributed in monocytes/macrophages. Additionally, we also found that their expression was associated with tumor microenvironment (including stromal cells and immune cells) and stemness score (DNAss and RNAss). Drug sensitivity analysis showed that PIAS family genes were able to predict the response to chemotherapy and immunotherapy. PIAS family genes expression is closely related to tumor metastasis, especially PIAS3. High PIAS3 expression significantly promotes the migration and invasion of liver cancer cell lines (HCC-LM3 and MHCC97-H).

CONCLUSIONS

Taking together, these findings contribute to determine whether the PIAS family genes are a potential oncogenic target gene, which have important contribution for the development of cancer immunotherapy.

Targeting pediatric cancers via T-cell recognition of the monomorphic MHC class I-related protein MR1

Human leukocyte antigen (HLA) restriction of conventional T-cell targeting introduces complexity in generating T-cell therapy strategies for patients with cancer with diverse HLA-backgrounds. A subpopulation of atypical, major histocompatibility complex-I related protein 1 (MR1)-restricted T-cells, distinctive from mucosal-associated invariant T-cells (MAITs), was recently identified recognizing currently unidentified MR1-presented cancer-specific metabolites. It is hypothesized that the MC.7.G5 MR1T-clone has potential as a pan-cancer, pan-population T-cell immunotherapy approach. These cells are irresponsive to healthy tissue while conferring T-cell receptor(TCR) dependent, HLA-independent cytotoxicity to a wide range of adult cancers. Studies so far are limited to adult malignancies. Here, we investigated the potential of MR1-targeting cellular therapy strategies in pediatric cancer. Bulk RNA sequencing data of primary pediatric tumors were analyzed to assess MR1 expression. In vitro pediatric tumor models were subsequently screened to evaluate their susceptibility to engineered MC.7.G5 TCR-expressing T-cells. Targeting capacity was correlated with qPCR-based MR1 mRNA and protein overexpression. RNA expression of MR1 in primary pediatric tumors varied widely within and between tumor entities. Notably, embryonal tumors exhibited significantly lower MR1 expression than other pediatric tumors. In line with this, most screened embryonal tumors displayed resistance to MR1T-targeting in vitro MR1T susceptibility was observed particularly in pediatric leukemia and diffuse midline glioma models. This study demonstrates potential of MC.7.G5 MR1T-cell immunotherapy in pediatric leukemias and diffuse midline glioma, while activity against embryonal tumors was limited. The dismal prognosis associated with relapsed/refractory leukemias and high-grade brain tumors highlights the promise to improve survival rates of children with these cancers.

Conditional c-MYC activation in catecholaminergic cells drives distinct neuroendocrine tumors: neuroblastoma vs somatostatinoma

The MYC proto-oncogenes (c-MYC, MYCN , MYCL ) are among the most deregulated oncogenic drivers in human malignancies including high-risk neuroblastoma, 50% of which are MYCN -amplified. Genetically engineered mouse models (GEMMs) based on the MYCN transgene have greatly expanded the understanding of neuroblastoma biology and are powerful tools for testing new therapies. However, a lack of c-MYC-driven GEMMs has hampered the ability to better understand mechanisms of neuroblastoma oncogenesis and therapy development given that c-MYC is also an important driver of many high-risk neuroblastomas. In this study, we report two transgenic murine neuroendocrine models driven by conditional c-MYC induction in tyrosine hydroxylase (Th) and dopamine β-hydroxylase (Dbh)-expressing cells. c-MYC induction in Th-expressing cells leads to a preponderance of Pdx1 + somatostatinomas, a type of pancreatic neuroendocrine tumor (PNET), resembling human somatostatinoma with highly expressed gene signatures of δ cells and potassium channels. In contrast, c-MYC induction in Dbh-expressing cells leads to onset of neuroblastomas, showing a better transforming capacity than MYCN in a comparable C57BL/6 genetic background. The c-MYC murine neuroblastoma tumors recapitulate the pathologic and genetic features of human neuroblastoma, express GD2, and respond to anti-GD2 immunotherapy. This model also responds to DFMO, an FDA-approved inhibitor targeting ODC1, which is a known MYC transcriptional target. Thus, establishing c-MYC-overexpressing GEMMs resulted in different but related tumor types depending on the targeted cell and provide useful tools for testing immunotherapies and targeted therapies for these diseases.

Targeted inhibition of SCFSKP2 confers anti-tumor activities resulting in a survival benefit in osteosarcoma

Osteosarcoma(OS) is a highly aggressive bone cancer for which treatment has remained essentially unchanged for decades. Although OS is characterized by extensive genomic heterogeneity and instability, RB1 and TP53 have been shown to be the most commonly inactivated tumor suppressors in OS. We previously generated a mouse model with a double knockout (DKO) of Rb1 and Trp53 within cells of the osteoblastic lineage, which largely recapitulates human OS with nearly complete penetrance. SKP2 is a repression target of pRb and serves as a substrate recruiting subunit of the SCFSKP2 complex. In addition, SKP2 plays a central role in regulating the cell cycle by ubiquitinating and promoting the degradation of p27. We previously reported the DKOAA transgenic model, which harbored a knock-in mutation in p27 that impaired its binding to SKP2. Here, we generated a novel p53-Rb1-SKP2 triple-knockout model (TKO) to examine SKP2 function and its potential as a therapeutic target in OS. First, we observed that OS tumorigenesis was significantly delayed in TKO mice and their overall survival was markedly improved. In addition, the loss of SKP2 also promoted an apoptotic microenvironment and reduced the stemness of DKO tumors. Furthermore, we found that small-molecule inhibitors of SKP2 exhibited anti-tumor activities in vivo and in OS organoids as well as synergistic effects when combined with a standard chemotherapeutic agent. Taken together, our results suggest that SKP2 inhibitors may reduce the stemness plasticity of OS and should be leveraged as next-generation adjuvants in this cancer.

BRG1 establishes the neuroectodermal chromatin landscape to restrict dorsal cell fates

Cell fate decisions are achieved with gene expression changes driven by lineage-specific transcription factors (TFs). These TFs depend on chromatin remodelers including the Brahma-related gene 1 (BRG1)-associated factor (BAF) complex to activate target genes. BAF complex subunits are essential for development and frequently mutated in cancer. Thus, interrogating how BAF complexes contribute to cell fate decisions is critical for human health. We examined the requirement for the catalytic BAF subunit BRG1 in neural progenitor cell (NPC) specification from human embryonic stem cells. During the earliest stages of differentiation, BRG1 was required to establish chromatin accessibility at neuroectoderm-specific enhancers. Depletion of BRG1 dorsalized NPCs and promoted precocious neural crest specification and enhanced neuronal differentiation. These findings demonstrate that BRG1 mediates NPC specification by ensuring proper expression of lineage-specific TFs and appropriate activation of their transcriptional programs.

The evolutionary impact of childhood cancer on the human gene pool

Germline pathogenic variants associated with increased childhood mortality must be subject to natural selection. Here, we analyze publicly available germline genetic metadata from 4,574 children with cancer [11 studies; 1,083 whole exome sequences (WES), 1,950 whole genome sequences (WGS), and 1,541 gene panel] and 141,456 adults [125,748 WES and 15,708 WGS]. We find that pediatric cancer predisposition syndrome (pCPS) genes [n = 85] are highly constrained, harboring only a quarter of the loss-of-function variants that would be expected. This strong indication of selective pressure on pCPS genes is found across multiple lines of germline genomics data from both pediatric and adult cohorts. For six genes [ELP1, GPR161, VHL and SDHA/B/C], a clear lack of mutational constraint calls the pediatric penetrance and/or severity of associated cancers into question. Conversely, out of 23 known pCPS genes associated with biallelic risk, two [9%, DIS3L2 and MSH2] show significant constraint, indicating that they may monoallelically increase childhood cancer risk. In summary, we show that population genetic data provide empirical evidence that heritable childhood cancer leads to natural selection powerful enough to have significantly impacted the present-day gene pool.

The dark side of stemness - the role of hematopoietic stem cells in development of blood malignancies

Hematopoietic stem cells (HSCs) produce all blood cells throughout the life of the organism. However, the high self-renewal and longevity of HSCs predispose them to accumulate mutations. The acquired mutations drive preleukemic clonal hematopoiesis, which is frequent among elderly people. The preleukemic state, although often asymptomatic, increases the risk of blood cancers. Nevertheless, the direct role of preleukemic HSCs is well-evidenced in adult myeloid leukemia (AML), while their contribution to other hematopoietic malignancies remains less understood. Here, we review the evidence supporting the role of preleukemic HSCs in different types of blood cancers, as well as present the alternative models of malignant evolution. Finally, we discuss the clinical importance of preleukemic HSCs in choosing the therapeutic strategies and provide the perspective on further studies on biology of preleukemic HSCs.

Pathophysiology of Acute Myeloid Leukemia

BACKGROUND

Acute myeloid leukemia (AML) is a biologically heterogenous disease arising in clonally proliferating hematopoietic stem cells. Sequential acquisition of mutations leads to expanded proliferation of clonal myeloid progenitors and failure of differentiation, leading to fulminant AML.

SUMMARY

Here, we review the pathophysiology of AML with a focus on factors predisposing to AML development, including prior chemo- and radiation therapy, environmental factors, and germline predisposition.

KEY MESSAGE

Increasing genomic characterization of AML and insight into mechanisms of its development will be critical to improvement in AML prognostication and therapy.

SJPedPanel: A pan-cancer gene panel for childhood malignancies

Background

Large scale genomics projects have identified driver alterations for most childhood cancers that provide reliable biomarkers for clinical diagnosis and disease monitoring using targeted sequencing. However, there is lack of a comprehensive panel that matches the list of known driver genes. Here we fill this gap by developing SJPedPanel for childhood cancers.

Results

SJPedPanel covers 5,275 coding exons of 357 driver genes, 297 introns frequently involved in rearrangements that generate fusion oncoproteins, commonly amplified/deleted regions (e.g., MYCN for neuroblastoma, CDKN2A and PAX5 for B-/T-ALL, and SMARCB1 for AT/RT), and 7,590 polymorphism sites for interrogating tumors with aneuploidy, such as hyperdiploid and hypodiploid B-ALL or 17q gain neuroblastoma. We used driver alterations reported from an established real-time clinical genomics cohort (n=253) to validate this gene panel. Among the 485 pathogenic variants reported, our panel covered 417 variants (86%). For 90 rearrangements responsible for oncogenic fusions, our panel covered 74 events (82%). We re-sequenced 113 previously characterized clinical specimens at an average depth of 2,500X using SJPedPanel and recovered 354 (91%) of the 389 reported pathogenic variants. We then investigated the power of this panel in detecting mutations from specimens with low tumor purity (as low as 0.1%) using cell line-based dilution experiments and discovered that this gene panel enabled us to detect ∼80% variants with allele fraction of 0.2%, while the detection rate decreases to ∼50% when the allele fraction is 0.1%. We finally demonstrate its utility in disease monitoring on clinical specimens collected from AML patients in morphologic remission.

Conclusions

SJPedPanel enables the detection of clinically relevant genetic alterations including rearrangements responsible for subtype-defining fusions for childhood cancers by targeted sequencing of ∼0.15% of human genome. It will enhance the analysis of specimens with low tumor burdens for cancer monitoring and early detection.

The role of abnormal ubiquitination in hepatocellular carcinoma pathology

Primary liver cancer is known for its high incidence and fatality rate. Over the years, therapeutic strategies for primary liver cancer have advanced significantly. Nonetheless, a substantial number of patients have not benefited from these methods, underscoring the pressing need for new and effective treatments for primary liver cancer. Ubiquitination is a critical post-translational modification that enables proteins to fulfill their normal biological functions and maintain their expression stability within cells. Importantly, increasing evidence suggests that the progression of liver cancer cells is often accompanied by disruptions in protein ubiquitination and deubiquitination processes. In this comprehensive review, we have compiled pertinent research about dysregulated ubiquitination in hepatocellular carcinoma (HCC) to broaden our understanding in this field. We elucidate the connections between the ubiquitination proteasome system, deubiquitination, and HCC. Furthermore, we shed light on the role of ubiquitination in cells situated within the tumor microenvironment of HCC including its involvement in mediating the activation of oncogenic pathways, reprogramming metabolic processes, and perturbing normal cellular functions. In conclusion, targeting the dysregulation of ubiquitination in HCC holds promise as a prospective and complementary therapeutic approach to existing treatments.

SKP2 Knockout in Rb1/p53-Deficient Mouse Models of Osteosarcoma Induces Immune Infiltration and Drives a Transcriptional Program with a Favorable Prognosis

Osteosarcoma is an aggressive bone malignancy with a poor prognosis. One putative proto-oncogene in osteosarcoma is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase. We previously demonstrated that Skp2 knockout in murine osteosarcoma improved survival and delayed tumorigenesis. Here, we performed RNA sequencing (RNA-seq) on tumors from a transgenic osteosarcoma mouse model with conditional Trp53 and Rb1 knockouts in the osteoblast lineage ("DKO": Osx1-Cre;Rb1lox/lox;p53lox/lox) and a triple-knockout model with additional Skp2 germline knockout ("TKO": Osx1-Cre;Rb1lox/lox;p53lox/lox;Skp2-/-), followed by qPCR and immunohistochemistry validation. To investigate the clinical implications of our results, we analyzed a human osteosarcoma patient cohort ("NCI-TARGET OS") with RNA-seq and clinical data. We found large differences in gene expression after SKP2 knockout. Surprisingly, we observed increased expression of genes related to immune microenvironment infiltration in TKO tumors, especially the signature genes for macrophages and to a lesser extent, T cells, B cells, and vascular cells. We also uncovered a set of relevant transcription factors that may mediate these changes. In osteosarcoma patient cohorts, high expression of genes upregulated in TKO was correlated with favorable overall survival, which was largely explained by the macrophage gene signatures. This relationship was further supported by our finding that SKP2 expression was negatively correlated with macrophage infiltration in the NCI-TARGET osteosarcoma and the TCGA Sarcoma cohorts. Overall, our findings indicate that SKP2 may mediate immune exclusion from the osteosarcoma tumor microenvironment, suggesting that SKP2 modulation in osteosarcoma may induce antitumor immune activation.

An Exome Capture-Based RNA-Sequencing Assay for Genome-Wide Identification and Prioritization of Clinically Important Fusions in Pediatric Tumors

This study reports the development of an exome capture-based RNA-sequencing assay to detect recurring and novel fusions in hematologic, solid, and central nervous system tumors. The assay used Twist Comprehensive Exome capture with either fresh or formalin-fixed samples and a bioinformatic platform that provides fusion detection, prioritization, and downstream curation. A minimum of 50 million uniquely mapped reads, a consensus read alignment/fusion calling approach using four callers (Arriba, FusionCatcher, STAR-Fusion, and Dragen), and custom software were used to integrate, annotate, and rank the candidate fusion calls. In an evaluation of 50 samples, the number of calls varied substantially by caller, from a mean of 24.8 with STAR-Fusion to 259.6 with FusionCatcher; only 1.1% of calls were made by all four callers. Therefore a filtering and ranking algorithm was developed based on multiple criteria, including number of supporting reads, calling consensus, genes involved, and cross-reference against databases of known cancer-associated or likely false-positive fusions. This approach was highly effective in pinpointing known clinically relevant fusions, ranking them first in 47 of 50 samples (94%). Detection of pathogenic gene fusions in three diagnostically challenging cases highlights the importance of a genome-wide and nontargeted method for fusion detection in pediatric cancer.

Divergent HLA variations and heterogeneous expression but recurrent HLA loss-of- heterozygosity and common HLA-B and TAP transcriptional silencing across advanced pediatric solid cancers

The human leukocyte antigen (HLA) system is a major factor controlling cancer immunosurveillance and response to immunotherapy, yet its status in pediatric cancers remains fragmentary. We determined high-confidence HLA genotypes in 576 children, adolescents and young adults with recurrent/refractory solid tumors from the MOSCATO-01 and MAPPYACTS trials, using normal and tumor whole exome and RNA sequencing data and benchmarked algorithms. There was no evidence for narrowed HLA allelic diversity but discordant homozygosity and allele frequencies across tumor types and subtypes, such as in embryonal and alveolar rhabdomyosarcoma, neuroblastoma MYCN and 11q subtypes, and high-grade glioma, and several alleles may represent protective or susceptibility factors to specific pediatric solid cancers. There was a paucity of somatic mutations in HLA and antigen processing and presentation (APP) genes in most tumors, except in cases with mismatch repair deficiency or genetic instability. The prevalence of loss-of-heterozygosity (LOH) ranged from 5.9 to 7.7% in HLA class I and 8.0 to 16.7% in HLA class II genes, but was widely increased in osteosarcoma and glioblastoma (~15-25%), and for DRB1-DQA1-DQB1 in Ewing sarcoma (~23-28%) and low-grade glioma (~33-50%). HLA class I and HLA-DR antigen expression was assessed in 194 tumors and 44 patient-derived xenografts (PDXs) by immunochemistry, and class I and APP transcript levels quantified in PDXs by RT-qPCR. We confirmed that HLA class I antigen expression is heterogeneous in advanced pediatric solid tumors, with class I loss commonly associated with the transcriptional downregulation of HLA-B and transporter associated with antigen processing (TAP) genes, whereas class II antigen expression is scarce on tumor cells and occurs on immune infiltrating cells. Patients with tumors expressing sufficient HLA class I and TAP levels such as some glioma, osteosarcoma, Ewing sarcoma and non-rhabdomyosarcoma soft-tissue sarcoma cases may more likely benefit from T cell-based approaches, whereas strategies to upregulate HLA expression, to expand the immunopeptidome, and to target TAP-independent epitopes or possibly LOH might provide novel therapeutic opportunities in others. The consequences of HLA class II expression by immune cells remain to be established. Immunogenetic profiling should be implemented in routine to inform immunotherapy trials for precision medicine of pediatric cancers.

Mobilizing China and the Global Community to Confront the Treatment Desert for Pediatric Solid Tumors

SUMMARY

Pediatric solid tumors are distinct clinical entities that impose heavy socioeconomic burden and while their incidence has increased in recent years, treatment options are often limited, with only 27 drugs approved for pediatric solid tumors in the United States, and fewer still, 13, in China. The scale of the unmet medical need is immense and new efforts are urgently needed to develop efficient therapeutics and improve these children's lives.

BARD1 germline variants induce haploinsufficiency and DNA repair defects in neuroblastoma

BACKGROUND

High-risk neuroblastoma is a complex genetic disease that is lethal in more than 50% of patients despite intense multimodal therapy. Through genome-wide association studies (GWAS) and next-generation sequencing, we have identified common single nucleotide polymorphisms and rare, pathogenic or likely pathogenic germline loss-of-function variants in BARD1 enriched in neuroblastoma patients. The functional implications of these findings remain poorly understood.

METHODS

We correlated BARD1 genotype with expression in normal tissues and neuroblastomas, along with the burden of DNA damage in tumors. To validate the functional consequences of germline pathogenic or likely pathogenic BARD1 variants, we used CRISPR-Cas9 to generate isogenic neuroblastoma (IMR-5) and control (RPE1) cellular models harboring heterozygous BARD1 loss-of-function variants (R112*, R150*, E287fs, and Q564*) and quantified genomic instability in these cells via next-generation sequencing and with functional assays measuring the efficiency of DNA repair.

RESULTS

Both common and rare neuroblastoma-associated BARD1 germline variants were associated with lower levels of BARD1 mRNA and an increased burden of DNA damage. Using isogenic heterozygous BARD1 loss-of-function variant cellular models, we functionally validated this association with inefficient DNA repair. BARD1 loss-of-function variant isogenic cells exhibited reduced efficiency in repairing Cas9-induced DNA damage, ineffective RAD51 focus formation at DNA double-strand break sites, and enhanced sensitivity to cisplatin and poly (ADP-ribose) polymerase (PARP) inhibition both in vitro and in vivo.

CONCLUSIONS

Taken together, we demonstrate that germline BARD1 variants disrupt DNA repair fidelity. This is a fundamental molecular mechanism contributing to neuroblastoma initiation that may have important therapeutic implications.

DriverDBv4: a multi-omics integration database for cancer driver gene research

Advancements in high-throughput technology offer researchers an extensive range of multi-omics data that provide deep insights into the complex landscape of cancer biology. However, traditional statistical models and databases are inadequate to interpret these high-dimensional data within a multi-omics framework. To address this limitation, we introduce DriverDBv4, an updated iteration of the DriverDB cancer driver gene database (http://driverdb.bioinfomics.org/). This updated version offers several significant enhancements: (i) an increase in the number of cohorts from 33 to 70, encompassing approximately 24 000 samples; (ii) inclusion of proteomics data, augmenting the existing types of omics data and thus expanding the analytical scope; (iii) implementation of multiple multi-omics algorithms for identification of cancer drivers; (iv) new visualization features designed to succinctly summarize high-context data and redesigned existing sections to accommodate the increased volume of datasets and (v) two new functions in Customized Analysis, specifically designed for multi-omics driver identification and subgroup expression analysis. DriverDBv4 facilitates comprehensive interpretation of multi-omics data across diverse cancer types, thereby enriching the understanding of cancer heterogeneity and aiding in the development of personalized clinical approaches. The database is designed to foster a more nuanced understanding of the multi-faceted nature of cancer.

MAJIQlopedia: an encyclopedia of RNA splicing variations in human tissues and cancer

Quantification of RNA splicing variations based on RNA-Sequencing can reveal tissue- and disease-specific splicing patterns. To study such splicing variations, we introduce MAJIQlopedia, an encyclopedia of splicing variations that encompasses 86 human tissues and 41 cancer datasets. MAJIQlopedia reports annotated and unannotated splicing events for a total of 486 175 alternative splice junctions in normal tissues and 338 317 alternative splice junctions in cancer. This database, available at https://majiq.biociphers.org/majiqlopedia/, includes a user-friendly interface that provides graphical representations of junction usage quantification for each junction across all tissue or cancer types. To demonstrate case usage of MAJIQlopedia, we review splicing variations in genes WT1, MAPT and BIN1, which all have known tissue or cancer-specific splicing variations. We also use MAJIQlopedia to highlight novel splicing variations in FDX1 and MEGF9 in normal tissues, and we uncover a novel exon inclusion event in RPS6KA6 that only occurs in two cancer types. Users can download the database, request the addition of data to the webtool, or install a MAJIQlopedia server to integrate proprietary data. MAJIQlopedia can serve as a reference database for researchers seeking to understand what splicing variations exist in genes of interest, and those looking to understand tissue- or cancer-specific splice isoform usage.

Prognostic value and immunological role of PD-L1 gene in pan-cancer

OBJECTIVE

PD-L1, a target of immune checkpoint blockade, has been proven to take the role of an oncogene in most human tumors. However, the role of PD-L1 in human pan-cancers has not yet been fully investigated.

MATERIALS AND METHODS

Pan-cancer analysis was conducted to analyze expression, genetic alterations, prognosis analysis, and immunological characteristics of PD-L1. Estimating the correlation between PD-L1 expression and survival involved using pooled odds ratios and hazard ratios with 95% CI. The Kaplan-Meier (K-M) technique, COX analysis, and receiver operating characteristic (ROC) curves were applied to the survival analysis. Additionally, we investigated the relationships between PD-L1 and microsatellite instability (MSI), tumor mutational burden (TMB), DNA methyltransferases (DNMTs), the associated genes of mismatch repair (MMR), and immune checkpoint biomarkers using Spearman's correlation analysis. Also, immunohistochemical analysis and qRT-PCR were employed in evaluating PD-L1's protein and mRNA expression in pan-caner.

RESULTS

PD-L1 showed abnormal mRNA and protein expression in a variety of cancers and predicted prognosis in cancer patients. Furthermore, across a variety of cancer types, the aberrant PD-L1 expression was connected to the MSI, MMR, TMB, drug sensitivity, and tumor immune microenvironment (TIME). Moreover, PD-L1 was significantly correlated with infiltrating levels of immune cells (T cell CD8 + , neutrophil, and so on).

CONCLUSION

Our study provides a better theoretical basis and guidance for the clinical treatment of PD-L1.

Non-genetic mechanisms of drug resistance in acute leukemias

Acute leukemia is characterized by clonal heterogeneity that contributes to poor drug responses in patients. Despite treatment advances, the occurrence of relapse remains a major barrier to achieving cures as current therapeutic approaches are inadequate to effectively prevent or overcome resistance. Given that only a few genetic mutations are associated with relapse in acute leukemia patients, there is a growing focus on 'non-genetic' mechanisms that affect the hallmarks of cancer to allow leukemic cells to survive post therapy. In this review, we provide an overview of the therapeutic landscape in acute leukemias. Importantly, we discuss non-genetic mechanisms exploited by leukemic cells to promote their survival after treatment. Last, we present current strategies to prevent or overcome drug resistance in this disease.