Systematic discovery of gene fusions in pediatric cancer by integrating RNA-seq and WGS.

BACKGROUND: Gene fusions are important cancer drivers in pediatric cancer and their accurate detection is essential for diagnosis and treatment. Clinical decision-making requires high confidence and precision of detection. Recent developments show RNA sequencing (RNA-seq) is promising for genome-wide detection of fusion products but hindered by many false positives that require extensive manual curation and impede discovery of pathogenic fusions. METHODS: We developed Fusion-sq to overcome existing disadvantages of detecting gene fusions. Fusion-sq integrates and "fuses" evidence from RNA-seq and whole genome sequencing (WGS) using intron-exon gene structure to identify tumor-specific protein coding gene fusions. Fusion-sq was then applied to the data generated from a pediatric pan-cancer cohort of 128 patients by WGS and RNA sequencing. RESULTS: In a pediatric pan-cancer cohort of 128 patients, we identified 155 high confidence tumor-specific gene fusions and their underlying structural variants (SVs). This includes all clinically relevant fusions known to be present in this cohort (30 patients). Fusion-sq distinguishes healthy-occurring from tumor-specific fusions and resolves fusions in amplified regions and copy number unstable genomes. A high gene fusion burden is associated with copy number instability. We identified 27 potentially pathogenic fusions involving oncogenes or tumor-suppressor genes characterized by underlying SVs, in some cases leading to expression changes indicative of activating or disruptive effects. CONCLUSIONS: Our results indicate how clinically relevant and potentially pathogenic gene fusions can be identified and their functional effects investigated by combining WGS and RNA-seq. Integrating RNA fusion predictions with underlying SVs advances fusion detection beyond extensive manual filtering. Taken together, we developed a method for identifying candidate gene fusions that is suitable for precision oncology applications. Our method provides multi-omics evidence for assessing the pathogenicity of tumor-specific gene fusions for future clinical decision making.

Elevated enhancer-oncogene contacts and higher oncogene expression levels by recurrent CTCF inactivating mutations in acute T cell leukemia.

Monoallelic inactivation of CCCTC-binding factor (CTCF) in human cancer drives altered methylated genomic states, altered CTCF occupancy at promoter and enhancer regions, and deregulated global gene expression. In patients with T cell acute lymphoblastic leukemia (T-ALL), we find that acquired monoallelic CTCF-inactivating events drive subtle and local genomic effects in nearly half of t(5; 14) (q35; q32.2) rearranged patients, especially when CTCF-binding sites are preserved in between the BCL11B enhancer and the TLX3 oncogene. These solitary intervening sites insulate TLX3 from the enhancer by inducing competitive looping to multiple binding sites near the TLX3 promoter. Reduced CTCF levels or deletion of the intervening CTCF site abrogates enhancer insulation by weakening competitive looping while favoring TLX3 promoter to BCL11B enhancer looping, which elevates oncogene expression levels and leukemia burden.

Molecular Characterization Reveals Subclasses of 1q Gain in Intermediate Risk Wilms Tumors.

Chromosomal alterations have recurrently been identified in Wilms tumors (WTs) and some are associated with poor prognosis. Gain of 1q (1q+) is of special interest given its high prevalence and is currently actively studied for its prognostic value. However, the underlying mutational mechanisms and functional effects remain unknown. In a national unbiased cohort of 30 primary WTs, we integrated somatic SNVs, CNs and SVs with expression data and distinguished four clusters characterized by affected biological processes: muscle differentiation, immune system, kidney development and proliferation. Combined genome-wide CN and SV profiles showed that tumors profoundly differ in both their types of 1q+ and genomic stability and can be grouped into WTs with co-occurring 1p-/1q+, multiple chromosomal gains or CN neutral tumors. We identified 1q+ in eight tumors that differ in mutational mechanisms, subsequent rearrangements and genomic contexts. Moreover, 1q+ tumors were present in all four expression clusters reflecting activation of various biological processes, and individual tumors overexpress different genes on 1q. In conclusion, by integrating CNs, SVs and gene expression, we identified subgroups of 1q+ tumors reflecting differences in the functional effect of 1q gain, indicating that expression data is likely needed for further risk stratification of 1q+ WTs.

Implementation of paediatric precision oncology into clinical practice: The Individualized Therapies for Children with cancer program 'iTHER'.

iTHER is a Dutch prospective national precision oncology program aiming to define tumour molecular profiles in children and adolescents with primary very high-risk, relapsed, or refractory paediatric tumours. Between April 2017 and April 2021, 302 samples from 253 patients were included. Comprehensive molecular profiling including low-coverage whole genome sequencing (lcWGS), whole exome sequencing (WES), RNA sequencing (RNA-seq), Affymetrix, and/or 850k methylation profiling was successfully performed for 226 samples with at least 20% tumour content. Germline pathogenic variants were identified in 16% of patients (35/219), of which 22 variants were judged causative for a cancer predisposition syndrome. At least one somatic alteration was detected in 204 (90.3%), and 185 (81.9%) were considered druggable, with clinical priority very high (6.1%), high (21.3%), moderate (26.0%), intermediate (36.1%), and borderline (10.5%) priority. iTHER led to revision or refinement of diagnosis in 8 patients (3.5%). Temporal heterogeneity was observed in paired samples of 15 patients, indicating the value of sequential analyses. Of 137 patients with follow-up beyond twelve months, 21 molecularly matched treatments were applied in 19 patients (13.9%), with clinical benefit in few. Most relevant barriers to not applying targeted therapies included poor performance status, as well as limited access to drugs within clinical trial. iTHER demonstrates the feasibility of comprehensive molecular profiling across all ages, tumour types and stages in paediatric cancers, informing of diagnostic, prognostic, and targetable alterations as well as reportable germline variants. Therefore, WES and RNA-seq is nowadays standard clinical care at the Princess Máxima Center for all children with cancer, including patients at primary diagnosis. Improved access to innovative treatments within biology-driven combination trials is required to ultimately improve survival.

Mesenchymal tumor organoid models recapitulate rhabdomyosarcoma subtypes.

Rhabdomyosarcomas (RMS) are mesenchyme-derived tumors and the most common childhood soft tissue sarcomas. Treatment is intense, with a nevertheless poor prognosis for high-risk patients. Discovery of new therapies would benefit from additional preclinical models. Here, we describe the generation of a collection of 19 pediatric RMS tumor organoid (tumoroid) models (success rate of 41%) comprising all major subtypes. For aggressive tumors, tumoroid models can often be established within 4-8 weeks, indicating the feasibility of personalized drug screening. Molecular, genetic, and histological characterization show that the models closely resemble the original tumors, with genetic stability over extended culture periods of up to 6 months. Importantly, drug screening reflects established sensitivities and the models can be modified by CRISPR/Cas9 with TP53 knockout in an embryonal RMS model resulting in replicative stress drug sensitivity. Tumors of mesenchymal origin can therefore be used to generate organoid models, relevant for a variety of preclinical and clinical research questions.

Improved Gene Fusion Detection in Childhood Cancer Diagnostics Using RNA Sequencing.

PURPOSE: Gene fusions play a significant role in cancer etiology, making their detection crucial for accurate diagnosis, prognosis, and determining therapeutic targets. Current diagnostic methods largely focus on either targeted or low-resolution genome-wide techniques, which may be unable to capture rare events or both fusion partners. We investigate if RNA sequencing can overcome current limitations with traditional diagnostic techniques to identify gene fusion events. METHODS: We first performed RNA sequencing on a validation cohort of 24 samples with a known gene fusion event, after which a prospective pan-pediatric cancer cohort (n = 244) was tested by RNA sequencing in parallel to existing diagnostic procedures. This cohort included hematologic malignancies, tumors of the CNS, solid tumors, and suspected neoplastic samples. All samples were processed in the routine diagnostic workflow and analyzed for gene fusions using standard-of-care methods and RNA sequencing. RESULTS: We identified a clinically relevant gene fusion in 83 of 244 cases in the prospective cohort. Sixty fusions were detected by both routine diagnostic techniques and RNA sequencing, and one fusion was detected only in routine diagnostics, but an additional 24 fusions were detected solely by RNA sequencing. RNA sequencing, therefore, increased the diagnostic yield by 38%-39%. In addition, RNA sequencing identified both gene partners involved in the gene fusion, in contrast to most routine techniques. For two patients, the newly identified fusion by RNA sequencing resulted in treatment with targeted agents. CONCLUSION: We show that RNA sequencing is sufficiently robust for gene fusion detection in routine diagnostics of childhood cancers and can make a difference in treatment decisions.

A systematic analysis of genetic interactions and their underlying biology in childhood cancer.

Childhood cancer is a major cause of child death in developed countries. Genetic interactions between mutated genes play an important role in cancer development. They can be detected by searching for pairs of mutated genes that co-occur more (or less) often than expected. Co-occurrence suggests a cooperative role in cancer development, while mutual exclusivity points to synthetic lethality, a phenomenon of interest in cancer treatment research. Little is known about genetic interactions in childhood cancer. We apply a statistical pipeline to detect genetic interactions in a combined dataset comprising over 2,500 tumors from 23 cancer types. The resulting genetic interaction map of childhood cancers comprises 15 co-occurring and 27 mutually exclusive candidates. The biological explanation of most candidates points to either tumor subtype, pathway epistasis or cooperation while synthetic lethality plays a much smaller role. Thus, other explanations beyond synthetic lethality should be considered when interpreting genetic interaction test results.

Structural variant detection in cancer genomes: computational challenges and perspectives for precision oncology.

Cancer is generally characterized by acquired genomic aberrations in a broad spectrum of types and sizes, ranging from single nucleotide variants to structural variants (SVs). At least 30% of cancers have a known pathogenic SV used in diagnosis or treatment stratification. However, research into the role of SVs in cancer has been limited due to difficulties in detection. Biological and computational challenges confound SV detection in cancer samples, including intratumor heterogeneity, polyploidy, and distinguishing tumor-specific SVs from germline and somatic variants present in healthy cells. Classification of tumor-specific SVs is challenging due to inconsistencies in detected breakpoints, derived variant types and biological complexity of some rearrangements. Full-spectrum SV detection with high recall and precision requires integration of multiple algorithms and sequencing technologies to rescue variants that are difficult to resolve through individual methods. Here, we explore current strategies for integrating SV callsets and to enable the use of tumor-specific SVs in precision oncology.

An organoid biobank for childhood kidney cancers that captures disease and tissue heterogeneity.

Kidney tumours are among the most common solid tumours in children, comprising distinct subtypes differing in many aspects, including cell-of-origin, genetics, and pathology. Pre-clinical cell models capturing the disease heterogeneity are currently lacking. Here, we describe the first paediatric cancer organoid biobank. It contains tumour and matching normal kidney organoids from over 50 children with different subtypes of kidney cancer, including Wilms tumours, malignant rhabdoid tumours, renal cell carcinomas, and congenital mesoblastic nephromas. Paediatric kidney tumour organoids retain key properties of native tumours, useful for revealing patient-specific drug sensitivities. Using single cell RNA-sequencing and high resolution 3D imaging, we further demonstrate that organoid cultures derived from Wilms tumours consist of multiple different cell types, including epithelial, stromal and blastemal-like cells. Our organoid biobank captures the heterogeneity of paediatric kidney tumours, providing a representative collection of well-characterised models for basic cancer research, drug-screening and personalised medicine.

DNA Methylation Profiling Identifies Distinct Clusters in Angiosarcomas.

PURPOSE: DNA methylation profiling has previously uncovered biologically and clinically meaningful subgroups within many tumor types, but was not yet performed in angiosarcoma. Angiosarcoma is a rare sarcoma with very heterogeneous clinical presentations, which may be based on differences in biological background. In this exploratory study, DNA methylation profiling of 36 primary angiosarcoma samples from visceral, deep soft tissue, radiation-induced, and UV-induced localizations was performed. EXPERIMENTAL DESIGN: Primary angiosarcoma formalin-fixed paraffin-embedded samples from visceral, soft tissue, radiation-induced, and UV-induced origin were collected from a nationwide search for angiosarcoma in the Netherlands. DNA was extracted for methylation profiling with the Illumina Infinium MethylationEPIC array. Quality control assessment and unsupervised hierarchical clustering were performed. Copy-number profiles were generated and analyzed for chromosomal stability. Clinical data were obtained from the Netherlands Cancer Registry. RESULTS: DNA methylation profiling by unsupervised hierarchical clustering of 36 angiosarcoma samples (6 visceral, 5 soft tissue, 14 radiation-induced, 11 UV-induced) revealed two main clusters (A and B), which were divided into four subclusters. The clusters largely corresponded with clinical subtypes, showing enrichment of UV-induced cases in cluster A1 and radiation-induced cases in cluster A2. Visceral and soft tissue cases almost exclusively fell into cluster B. Cluster A showed significantly increased chromosomal instability and better overall survival (22 vs. 6 months, P = 0.046) compared with cluster B. CONCLUSIONS: In this novel methylation profiling study, we demonstrated for the first time four different angiosarcoma clusters. These clusters correlated with clinical subtype, overall survival, and chromosomal stability.

Kaposiform hemangioendothelioma and tufted angioma - (epi)genetic analysis including genome-wide methylation profiling.

Kaposiform hemangioendothelioma (KHE) is a locally aggressive vascular condition of childhood and is clinicopathologically related to tufted angioma (TA), a benign skin lesion. Due to their rarity molecular data are scarce. We investigated 7 KHE and 3 TA by comprehensive mutational analysis and genome-wide methylation profiling and compared the clustering, also with vascular malformations. Lesions were from 7 females and 3 males. The age range was 2 months to 9 years with a median of 10 months. KHEs arose in the soft tissue of the thigh (n = 2), retroperitoneum (n = 1), thoracal/abdominal (n = 1), supraclavicular (n = 1) and neck (n = 1). One patient presented with multiple lesions without further information. Two patients developed a Kasabach-Merritt phenomenon. TAs originated in the skin of the shoulder (n = 2) and nose/forehead (n = 1). Of the 5 KHEs and 2 TAs investigated by DNA sequencing, one TA showed a hot spot mutation in NRAS, and one KHE a mutation in RAD50. Unsupervised hierarchical clustering analysis indicated a common methylation pattern of KHEs and TAs, which separated from the homogeneous methylation pattern of vascular malformations. In conclusion, methylation profiling provides further evidence for KHEs and TAs potentially forming a spectrum of one entity. Using next generation sequencing, heterogeneous mutations were found in a subset of cases (2/7) without the presence of GNA14 mutations, previously reported in KHE and TA.

The clinical implementation of copy number detection in the age of next-generation sequencing.

INTRODUCTION: The role of copy number variants (CNVs) in disease is now well established. In parallel NGS technologies, such as long-read technologies, there is continual development and data analysis methods continue to be refined. Clinical exome sequencing data is now a reality for many diagnostic laboratories in both congenital genetics and oncology. This provides the ability to detect and report both SNVs and structural variants, including CNVs, using a single assay for a wide range of patient cohorts. Areas covered: Currently, whole-genome sequencing is mainly restricted to research applications and clinical utility studies. Furthermore, detecting the full-size spectrum of CNVs as well as somatic events remains difficult for both exome and whole-genome sequencing. As a result, the full extent of genomic variants in an individual's genome is still largely unknown. Recently, new sequencing technologies have been introduced which maintain the long-range genomic context, aiding the detection of CNVs and structural variants. Expert commentary: The development of long-read sequencing promises to resolve many CNV and SV detection issues but is yet to become established. The current challenge for clinical CNV detection is how to fully exploit all the data which is generated by high throughput sequencing technologies.

Molecular mechanisms that distinguish TFIID housekeeping from regulatable SAGA promoters.

An important distinction is frequently made between constitutively expressed housekeeping genes versus regulated genes. Although generally characterized by different DNA elements, chromatin architecture and cofactors, it is not known to what degree promoter classes strictly follow regulatability rules and which molecular mechanisms dictate such differences. We show that SAGA-dominated/TATA-box promoters are more responsive to changes in the amount of activator, even compared to TFIID/TATA-like promoters that depend on the same activator Hsf1. Regulatability is therefore an inherent property of promoter class. Further analyses show that SAGA/TATA-box promoters are more dynamic because TATA-binding protein recruitment through SAGA is susceptible to removal by Mot1. In addition, the nucleosome configuration upon activator depletion shifts on SAGA/TATA-box promoters and seems less amenable to preinitiation complex formation. The results explain the fundamental difference between housekeeping and regulatable genes, revealing an additional facet of combinatorial control: an activator can elicit a different response dependent on core promoter class.

A comprehensive framework of E2-RING E3 interactions of the human ubiquitin-proteasome system.

Covalent attachment of ubiquitin to substrates is crucial to protein degradation, transcription regulation and cell signalling. Highly specific interactions between ubiquitin-conjugating enzymes (E2) and ubiquitin protein E3 ligases fulfil essential roles in this process. We performed a global yeast-two hybrid screen to study the specificity of interactions between catalytic domains of the 35 human E2s with 250 RING-type E3s. Our analysis showed over 300 high-quality interactions, uncovering a large fraction of new E2-E3 pairs. Both within the E2 and the E3 cohorts, several members were identified that are more versatile in their interaction behaviour than others. We also found that the physical interactions of our screen compare well with reported functional E2-E3 pairs in in vitro ubiquitination experiments. For validation we confirmed the interaction of several versatile E2s with E3s in in vitro protein interaction assays and we used mutagenesis to alter the E3 interactions of the E2 specific for K63 linkages, UBE2N(Ubc13), towards the K48-specific UBE2D2(UbcH5B). Our data provide a detailed, genome-wide overview of binary E2-E3 interactions of the human ubiquitination system.

Cotranslational assembly of the yeast SET1C histone methyltransferase complex.

While probing the role of RNA for the function of SET1C/COMPASS histone methyltransferase, we identified SET1RC (SET1 mRNA-associated complex), a complex that contains SET1 mRNA and Set1, Swd1, Spp1 and Shg1, four of the eight polypeptides that constitute SET1C. Characterization of SET1RC showed that SET1 mRNA binding did not require associated Swd1, Spp1 and Shg1 proteins or RNA recognition motifs present in Set1. RNA binding was not observed when Set1 protein and SET1 mRNA were derived from independent genes or when SET1 transcripts were restricted to the nucleus. Importantly, the protein-RNA interaction was sensitive to EDTA, to the translation elongation inhibitor puromycin and to the inhibition of translation initiation in prt1-1 mutants. Taken together, our results support the idea that SET1 mRNA binding was dependent on translation and that SET1RC assembled on nascent Set1 in a cotranslational manner. Moreover, we show that cellular accumulation of Set1 is limited by the availability of certain SET1C components, such as Swd1 and Swd3, and suggest that cotranslational protein interactions may exert an effect in the protection of nascent Set1 from degradation.

CCR4/NOT complex associates with the proteasome and regulates histone methylation.

The proteasome regulates histone lysine methylation and gene transcription, but how it does so is poorly understood. To better understand this process, we used the epistatic miniarray profile (E-MAP) approach to identify factors that genetically interact with proteasomal subunits. In addition to members of the Set1 complex that mediate histone H3 lysine 4 methylation (H3K4me), we found that deleting members of the CCR4/NOT mRNA processing complex exhibit synthetic phenotypes when combined with proteasome mutants. Further biochemical analyses revealed physical associations between CCR4/NOT and the proteasome in vivo. Consistent with the genetic and biochemical interactions linking CCR4/NOT with proteasome and Set1-mediated methylation, we find that loss of Not4 decreases global and gene-specific H3K4 trimethylation (H3K4me3) and decreases 19S proteasome recruitment to the PMA1 gene. Similar to proteasome regulation of histone methylation, loss of CCR4/NOT members does not affect ubiquitinated H2B. Mapping of Not4 identified the RING finger domain as essential for H3K4me3, suggesting a role for ubiquitin in this process. Consistent with this idea, loss of the Not4-interacting protein Ubc4, a known ubiquitin-conjugating enzyme, decreases H3K4me3. These studies implicate CCR4/NOT in the regulation of H3K4me3 through a ubiquitin-dependent pathway that likely involves the proteasome.

Transcriptome-based identification of pro- and antioxidative gene expression in kidney cortex of nitric oxide-depleted rats.

Nitric oxide (NO) depletion in rats induces severe endothelial dysfunction within 4 days. Subsequently, hypertension and renal injury develop, which are ameliorated by alpha-tocopherol (VitE) cotreatment. The hypothesis of the present study was that NO synthase (NOS) inhibition induces a renal cortical antioxidative transcriptional response and invokes pro-oxidative and proinflammatory gene expression due to elimination of dampening effects of NO and enhanced oxidative stress. Male Sprague-Dawley rats received NOS inhibitor N(omega)-nitro-l-arginine (l-NNA, 500 mg/l water) for 4 (4d-LNNA), 21 (21d-LNNA), or 21 days with VitE in chow (0.7 g/kg body wt/day). Renal cortical RNA was applied to oligonucleotide rat arrays. In 4d-LNNA, 21d-LNNA, and 21d-LNNA+VitE, 120, 320, and 184 genes were differentially expressed, respectively. Genes related to glutathione and bilirubin synthesis were suppressed during 4d and 21d-LNNA and not corrected by VitE. Proteinuria, tubulointerstitial macrophages, and heme-oxygenase-1 (HO-1) expression were strongly correlated. Remarkably, pro-oxidative genes were not induced. Inflammation- and injury-related genes, including kidney injury molecule-1 and osteopontin, were unchanged at day 4, induced at 21d, and partly corrected by VitE. Superimposing HO-1 inhibition on NOS inhibition had no impact on the development of hypertension. To summarize, renal expression of genes involved in synthesis of the antioxidants glutathione and bilirubin seemed directly NO dependent, but there were no direct effects of NO depletion on pro-oxidant systems. This indicates that renal transcriptional regulation of two defense systems, glutathione and bilirubin syntheses, seems to depend upon adequate NO synthesis. Interaction between NO synthesis and heme degradation pathways for blood pressure regulation was not found.

Multiple robust signatures for detecting lymph node metastasis in head and neck cancer.

Genome-wide mRNA expression measurements can identify molecular signatures of cancer and are anticipated to improve patient management. Such expression profiles are currently being critically evaluated based on an apparent instability in gene composition and the limited overlap between signatures from different studies. We have recently identified a primary tumor signature for detection of lymph node metastasis in head and neck squamous cell carcinomas. Before starting a large multicenter prospective validation, we have thoroughly evaluated the composition of this signature. A multiple training approach was used for validating the original set of predictive genes. Based on different combinations of training samples, multiple signatures were assessed for predictive accuracy and gene composition. The initial set of predictive genes is a subset of a larger group of 825 genes with predictive power. Many of the predictive genes are interchangeable because of a similar expression pattern across the tumor samples. The head and neck metastasis signature has a more stable gene composition than previous predictors. Exclusion of the strongest predictive genes could be compensated by raising the number of genes included in the signature. Multiple accurate predictive signatures can be designed using various subsets of predictive genes. The absence of genes with strong predictive power can be compensated by including more genes with lower predictive power. Lack of overlap between predictive signatures from different studies with the same goal may be explained by the fact that there are more predictive genes than required to design an accurate predictor.

An expression profile for diagnosis of lymph node metastases from primary head and neck squamous cell carcinomas.

Metastasis is the process by which cancers spread to distinct sites in the body. It is the principal cause of death in individuals suffering from cancer. For some types of cancer, early detection of metastasis at lymph nodes close to the site of the primary tumor is pivotal for appropriate treatment. Because it can be difficult to detect lymph node metastases reliably, many individuals currently receive inappropriate treatment. We show here that DNA microarray gene-expression profiling can detect lymph node metastases for primary head and neck squamous cell carcinomas that arise in the oral cavity and oropharynx. The predictor, established with an 82-tumor training set, outperforms current clinical diagnosis when independently validated. The 102 predictor genes offer unique insights into the processes underlying metastasis. The results show that the metastatic state can be deciphered from the primary tumor gene-expression pattern and that treatment can be substantially improved.

Monitoring global messenger RNA changes in externally controlled microarray experiments.

Expression profiling is a universal tool, with a range of applications that benefit from the accurate determination of differential gene expression. To allow normalization using endogenous transcript levels, current microarray analyses assume that relatively few transcripts vary, or that any changes that occur are balanced. When normalization using endogenous genes is carried out, changes in expression levels are calculated relative to the behaviour of most of the transcripts. This does not reflect absolute changes if global shifts in messenger RNA populations occur. Using external RNA controls, we have set up microarray experiments to monitor global changes. The levels of most mRNAs were found to change during yeast stationary phase and human heat shock when external controls were included. Even small global changes had a significant effect on the number of genes reported as being differentially expressed. This suggests that global mRNA changes occur more frequently than is assumed at present, and shows that monitoring such effects may be important for the accurate determination of changes in gene expression.