Integrative analysis of gene expression and copy number alterations using canonical correlation analysis

Supplementary Figure 1. Representation of the samples from the tuning set by their coordinates in the first two pairs of features (extracted from the tuning set) using regularized dual CCA, with regularization parameters tx = 0.9, ty = 0.3 (left panel), and PCA+CCA (right panel). We show the representations with respect to both the copy number features and the gene expression features in a superimposed way, where each sample is represented by two markers. The filled markers represent the coordinates in the features extracted from the copy number variables, and the open markers represent coordinates in the features extracted from the gene expression variables. Samples with different leukemia subtypes are shown with different colors. The first feature pair distinguishes the HD50 group from the rest, while the second feature pair represents the characteristics of the samples from the E2A/PBX1 subtype. The high canonical correlation obtained for the tuning samples with regularized dual CCA is apparent in the left panel, where the two points for each sample coincide. Nevertheless, the extracted features have a high generalization ability, as can be seen in the left panel of Figure 5, showing the representation of the validation samples. 1 Supplementary Figure 2. Representation of the samples from the tuning set by their coordinates in the first two pairs of features (extracted from the tuning set) using regularized dual CCA, with regularization parameters tx = 0, ty = 0 (left panel), and tx = 1, ty = 1 (right panel). We show the representations with respect to both the copy number features and the gene expression features in a superimposed way, where each sample is represented by tw

Proteogenomics and Hi-C reveal transcriptional dysregulation in high hyperdiploid childhood acute lymphoblastic leukemia.

Hyperdiploidy, i.e. gain of whole chromosomes, is one of the most common genetic features of childhood acute lymphoblastic leukemia (ALL), but its pathogenetic impact is poorly understood. Here, we report a proteogenomic analysis on matched datasets from genomic profiling, RNA-sequencing, and mass spectrometry-based analysis of >8,000 genes and proteins as well as Hi-C of primary patient samples from hyperdiploid and ETV6/RUNX1-positive pediatric ALL. We show that CTCF and cohesin, which are master regulators of chromatin architecture, display low expression in hyperdiploid ALL. In line with this, a general genome-wide dysregulation of gene expression in relation to topologically associating domain (TAD) borders were seen in the hyperdiploid group. Furthermore, Hi-C of a limited number of hyperdiploid childhood ALL cases revealed that 2/4 cases displayed a clear loss of TAD boundary strength and 3/4 showed reduced insulation at TAD borders, with putative leukemogenic effects

The correlation pattern of acquired copy number changes in 164 ETV6/RUNX1-positive childhood acute lymphoblastic leukemias

The ETV6/RUNX1 fusion gene, present in 25% of B-lineage childhood acute lymphoblastic leukemia (ALL), is thought to represent an initiating event, which requires additional genetic changes for leukemia development. To identify additional genetic alterations, 24 ETV6/RUNX1-positive ALLs were analyzed using 500K single nucleotide polymorphism arrays. The results were combined with previously published data sets, allowing us to ascertain genomic copy number aberrations (CNAs) in 164 cases. In total, 45 recurrent CNAs were identified with an average number of 3.5 recurrent changes per case (range 0–13). Twenty-six percent of cases displayed a set of recurrent CNAs identical to that of other cases in the data set. The majority (74%), however, displayed a unique pattern of recurrent CNAs, indicating a large heterogeneity within this ALL subtype. As previously demonstrated, alterations targeting genes involved in B-cell development were common (present in 28% of cases). However, the combined analysis also identified alterations affecting nuclear hormone response (24%) to be a characteristic feature of ETV6/RUNX1-positive ALL. Studying the correlation pattern of the CNAs allowed us to highlight significant positive and negative correlations between specific aberrations. Furthermore, oncogenetic tree models identified ETV6, CDKN2A/B, PAX5, del(6q) and +16 as possible early events in the leukemogenic process

Clonal competition within complex evolutionary hierarchies shapes AML over time

Clonal heterogeneity and evolution has major implications for disease progression and relapse in acute myeloid leukemia (AML). To model clonal dynamics in vivo, we serially transplanted 23 AML cases to immunodeficient mice and followed clonal composition for up to 15 months by whole-exome sequencing of 84 xenografts across two generations. We demonstrate vast changes in clonality that both progress and reverse over time, and define five patterns of clonal dynamics: Monoclonal, Stable, Loss, Expansion and Burst. We also show that subclonal expansion in vivo correlates with a more adverse prognosis. Furthermore, clonal expansion enabled detection of very rare clones with AML driver mutations that were undetectable by sequencing at diagnosis, demonstrating that the vast majority of AML cases harbor multiple clones already at diagnosis. Finally, the rise and fall of related clones enabled deconstruction of the complex evolutionary hierarchies of the clones that compete to shape AML over time

GRM1 is upregulated through gene fusion and promoter swapping in chondromyxoid fibroma.

Glutamate receptors are well-known actors in the central and peripheral nervous systems, and altered glutamate signaling is implicated in several neurological and psychiatric disorders. It is increasingly recognized that such receptors may also have a role in tumor growth. Here we provide direct evidence of aberrant glutamate signaling in the development of a locally aggressive bone tumor, chondromyxoid fibroma (CMF). We subjected a series of CMFs to whole-genome mate-pair sequencing and RNA sequencing and found that the glutamate receptor gene GRM1 recombines with several partner genes through promoter swapping and gene fusion events. The GRM1 coding region remains intact, and 18 of 20 CMFs (90%) showed a more than 100-fold and up to 1,400-fold increase in GRM1 expression levels compared to control tissues. Our findings unequivocally demonstrate that direct targeting of GRM1 is a necessary and highly specific driver event for CMF development