The co-evolution of the genome and epigenome in colorectal cancer.

Colorectal malignancies are a leading cause of cancer-related death1 and have undergone extensive genomic study2,3. However, DNA mutations alone do not fully explain malignant transformation4-7. Here we investigate the co-evolution of the genome and epigenome of colorectal tumours at single-clone resolution using spatial multi-omic profiling of individual glands. We collected 1,370 samples from 30 primary cancers and 8 concomitant adenomas and generated 1,207 chromatin accessibility profiles, 527 whole genomes and 297 whole transcriptomes. We found positive selection for DNA mutations in chromatin modifier genes and recurrent somatic chromatin accessibility alterations, including in regulatory regions of cancer driver genes that were otherwise devoid of genetic mutations. Genome-wide alterations in accessibility for transcription factor binding involved CTCF, downregulation of interferon and increased accessibility for SOX and HOX transcription factor families, suggesting the involvement of developmental genes during tumourigenesis. Somatic chromatin accessibility alterations were heritable and distinguished adenomas from cancers. Mutational signature analysis showed that the epigenome in turn influences the accumulation of DNA mutations. This study provides a map of genetic and epigenetic tumour heterogeneity, with fundamental implications for understanding colorectal cancer biology.

Phenotypic plasticity and genetic control in colorectal cancer evolution.

Genetic and epigenetic variation, together with transcriptional plasticity, contribute to intratumour heterogeneity1. The interplay of these biological processes and their respective contributions to tumour evolution remain unknown. Here we show that intratumour genetic ancestry only infrequently affects gene expression traits and subclonal evolution in colorectal cancer (CRC). Using spatially resolved paired whole-genome and transcriptome sequencing, we find that the majority of intratumour variation in gene expression is not strongly heritable but rather 'plastic'. Somatic expression quantitative trait loci analysis identified a number of putative genetic controls of expression by cis-acting coding and non-coding mutations, the majority of which were clonal within a tumour, alongside frequent structural alterations. Consistently, computational inference on the spatial patterning of tumour phylogenies finds that a considerable proportion of CRCs did not show evidence of subclonal selection, with only a subset of putative genetic drivers associated with subclone expansions. Spatial intermixing of clones is common, with some tumours growing exponentially and others only at the periphery. Together, our data suggest that most genetic intratumour variation in CRC has no major phenotypic consequence and that transcriptional plasticity is, instead, widespread within a tumour.

Treatment with the apoptosis inhibitor Asunercept reduces clone sizes in patients with lower risk Myelodysplastic Neoplasms.

In low-risk Myelodysplastic Neoplasms (MDS), increased activity of apoptosis-promoting factors such as tumor necrosis factor (TNFα) and pro-apoptotic Fas ligand (CD95L) have been described as possible pathomechanisms leading to impaired erythropoiesis. Asunercept (APG101) is a novel therapeutic fusion protein blocking CD95, which has previously shown partial efficacy in reducing transfusion requirement in a clinical phase I trial for low-risk MDS patients (NCT01736436; 2012-11-26). In the current study we aimed to evaluate the effect of Asunercept therapy on the clonal bone marrow composition to identify potential biomarkers to predict response. Bone marrow samples of n = 12 low-risk MDS patients from the above referenced clinical trial were analyzed by serial deep whole exome sequencing in a total of n = 58 time points. We could distinguish a mean of 3.5 molecularly defined subclones per patient (range 2-6). We observed a molecular response defined as reductions of dominant clone sizes by a variant allele frequency (VAF) decrease of at least 10% (mean 20%, range: 10.5-39.2%) in dependency of Asunercept treatment in 9 of 12 (75%) patients. Most of this decline in clonal populations was observed after completion of 12 weeks treatment. Particularly early and pronounced reductions of clone sizes were found in subclones driven by mutations in genes involved in regulation of methylation (n = 1 DNMT3A, n = 1 IDH2, n = 1 TET2). Our results suggest that APG101 could be efficacious in reducing clone sizes of mutated hematopoietic cells in the bone marrow of Myelodysplastic Neoplasms, which warrants further investigation.

High erythroferrone expression in CD71+ erythroid progenitors predicts superior survival in myelodysplastic syndromes.

Ineffective erythropoiesis and iron overload are common in myelodysplastic syndromes (MDS). Erythroferrone (ERFE) and growth/differentiation factor 15 (GDF15) are two regulators of iron homeostasis produced by erythroid progenitors. Elevated systemic levels of ERFE and GDF15 in MDS are associated with dysregulated iron metabolism and iron overload, which is especially pronounced in MDS with SF3B1 gene mutations. However, the role of ERFE and GDF15 in MDS pathogenesis and their influence on disease progression are largely unknown. Here, we analyzed the expression of ERFE and GDF15 in CD71+ erythroid progenitors of n = 111 MDS patients and assessed their effects on patient survival. The expression of ERFE and GDF15 in MDS was highly aberrant. Unexpectedly, ERFE expression in erythroprogenitors was highly relevant for MDS prognosis and independent of International Prognostic Scoring System (IPSS) stratification. Although ERFE expression was increased in patients with SF3B1 mutations, it predicted overall survival (OS) in both the SF3B1wt and SF3B1mut subgroups. Of note, ERFE overexpression predicted superior OS in the IPSS low/Int-1 subgroup and in patients with normal karyotype. Similar observations were made for GDF15, albeit not reaching statistical significance. In summary, our results revealed a strong association between ERFE expression and MDS outcome, suggesting a possible involvement of ERFE in molecular MDS pathogenesis.

Mutational hierarchies in myelodysplastic syndromes dynamically adapt and evolve upon therapy response and failure.

Clonal evolution is believed to be a main driver for progression of various types of cancer and implicated in facilitating resistance to drugs. However, the hierarchical organization of malignant clones in the hematopoiesis of myelodysplastic syndromes (MDS) and its impact on response to drug therapy remain poorly understood. Using high-throughput sequencing of patient and xenografted cells, we evaluated the intratumoral heterogeneity (n= 54) and reconstructed mutational trajectories (n = 39) in patients suffering from MDS (n = 52) and chronic myelomonocytic leukemia-1 (n = 2). We identified linear and also branching evolution paths and confirmed on a patient-specific level that somatic mutations in epigenetic regulators and RNA splicing genes frequently constitute isolated disease-initiating events. Using high-throughput exome- and/or deep-sequencing, we analyzed 103 chronologically acquired samples from 22 patients covering a cumulative observation time of 75 years MDS disease progression. Our data revealed highly dynamic shaping of complex oligoclonal architectures, specifically upon treatment with lenalidomide and other drugs. Despite initial clinical response to treatment, patients' marrow persistently remained clonal with rapid outgrowth of founder-, sub-, or even fully independent clones, indicating an increased dynamic rate of clonal turnover. The emergence and disappearance of specific clones frequently correlated with changes of clinical parameters, highlighting their distinct and far-reaching functional properties. Intriguingly, increasingly complex mutational trajectories are frequently accompanied by clinical progression during the course of disease. These data substantiate a need for regular broad molecular monitoring to guide clinical treatment decisions in MDS.

Accurate quantification of chromosomal lesions via short tandem repeat analysis using minimal amounts of DNA.

BACKGROUND: Cytogenetic aberrations such as deletion of chromosome 5q (del(5q)) represent key elements in routine clinical diagnostics of haematological malignancies. Currently established methods such as metaphase cytogenetics, FISH or array-based approaches have limitations due to their dependency on viable cells, high costs or semi-quantitative nature. Importantly, they cannot be used on low abundance DNA. We therefore aimed to establish a robust and quantitative technique that overcomes these shortcomings. METHODS: For precise determination of del(5q) cell fractions, we developed an inexpensive multiplex-PCR assay requiring only nanograms of DNA that simultaneously measures allelic imbalances of 12 independent short tandem repeat markers. RESULTS: Application of this method to n=1142 samples from n=260 individuals revealed strong intermarker concordance (R²=0.77-0.97) and reproducibility (mean SD: 1.7%). Notably, the assay showed accurate quantification via standard curve assessment (R²>0.99) and high concordance with paired FISH measurements (R²=0.92) even with subnanogram amounts of DNA. Moreover, cytogenetic response was reliably confirmed in del(5q) patients with myelodysplastic syndromes treated with lenalidomide. While the assay demonstrated good diagnostic accuracy in receiver operating characteristic analysis (area under the curve: 0.97), we further observed robust correlation between bone marrow and peripheral blood samples (R²=0.79), suggesting its potential suitability for less-invasive clonal monitoring. CONCLUSIONS: In conclusion, we present an adaptable tool for quantification of chromosomal aberrations, particularly in problematic samples, which should be easily applicable to further tumour entities.

Concomitant MDS with isolated 5q deletion and MGUS: case report and review of molecular aspects.

Patients with monoclonal gammopathy of undetermined significance (MGUS) have a higher risk for the development of concomitant primary cancers such as multiple myeloma (MM) and myelodysplastic syndrome (MDS). We report the case of patient initially suffering from MGUS of the IgG lambda subtype for more than 10 yr, which evolved to MM and MDS with deletion (5q) with severe pancytopenia. Due to pancytopenia, he received dose-reduced treatment with lenalidomide and dexamethasone. He achieved an ongoing transfusion independency after about 1 month of treatment. Bone marrow taken 14 months after start of treatment showed a complete cytogenetic response of the del(5q) clone and a plasma cell infiltration below 5%. In contrast to the development of MM in MGUS patients, the subsequent occurrence of MDS after diagnosis of MGUS is infrequent. Moreover, the biological association of MDS with MGUS is not sufficiently understood, but the non-treatment-related occurrence supports the pathogenetic role of pre-existing alterations of stem cells. Here, we summarize data on concomitant MDS and MGUS/MM with particular emphasis on molecular aspects.

Skewed X-inactivation patterns in ageing healthy and myelodysplastic haematopoiesis determined by a pyrosequencing based transcriptional clonality assay.

BACKGROUND: Investigation of X-chromosome inactivation patterns (XCIP) by determination of differential CpG-methylation has been widely applied for investigation of female cell clonality. Using this approach the clonal origin of various tumours has been corroborated. Controversially, strong age-related increase of peripheral blood (PB) cell clonality in haematologically healthy female subjects was reported. Recently, transcriptional XCIP ratio analysis challenged these results and questioned the suitability of methylation based clonality assays. METHODS: To reinvestigate XCIP-skewing in CD34, low-density mononuclear bone marrow (BM) as well as PB cells from healthy female subjects and patients with myelodysplastic syndromes (MDS), we established a transcriptional assay using pyrosequencing technique for quantification of single nucleotide polymorphism allele frequencies, representative for XCIP ratios. RESULTS: Our assay provides high sensitivity for XCIP ratio assessment as determined by standard curves, reproducibility, inter-marker correlation as well as correlation with the DNA-methylation based human androgen receptor (HUMARA) assay. Notably, in agreement with most studies investigating this issue, significant age-related increase of XCIP skewing in PB cells from healthy elderly female subjects was confirmed. Moreover, XCIP ratio analysis suggests even stronger clonal manifestation in BM and CD34 cells. In MDS, XCIP skewing levels were distinctively elevated as compared with controls of similar age and higher degrees were associated with poor clinical outcome. CONCLUSIONS: Transcriptional clonal profiling via pyrosequencing allows accurate assessment of XCIP ratios, confirms the validity of the DNA-methylation based HUMARA assay and reveals important insights into ageing healthy and myelodysplastic haematopoiesis.

Defining an Optimized Workflow for Enriching and Analyzing Residual Tumor Populations Using Intracellular Markers.

Tumor relapse is well recognized to arise from treatment-resistant residual populations. Strategies enriching such populations for in-depth downstream analyses focus on tumor-specific surface markers; however, enrichment using intracellular biomarkers remains challenging. Using B-cell lymphoma as an exemplar, we demonstrate feasibility to enrich B-cell lymphoma 2 (BCL2)high populations, a surrogate marker for t(14;18)+ lymphomas, for use in downstream applications. Different fixation protocols were assessed for impact on antibody expression and RNA integrity; glyoxal fixation demonstrated superior results regarding minimal effects on surface and intracellular expression, and RNA quality, compared with alternative fixatives evaluated. Furthermore, t(14;18)+ B cells were effectively detected using intracellular BCL2 overexpression to facilitate tumor cell enrichment. Tumor cell populations were enriched using the cellenONE F1.4 single-cell sorting platform, which detected and dispensed BCL2high-expressing cells directly into library preparation reagents for transcriptome analyses. Sorted glyoxal-fixed cells generated good quality sequencing libraries, with high concordance between live and fixed single-cell transcriptomic profiles, discriminating cell populations predominantly on B-cell biology. Overall, we successfully developed a proof-of-concept workflow employing a robust cell preparation protocol for intracellular markers combined with cell enrichment using the cellenONE platform, providing an alternative to droplet-based technologies when cellular input is low or requires prior enrichment to detect rare populations. This workflow has wider prognostic and therapeutic potential to study residual cells in a pan-cancer setting.

Epigenome and early selection determine the tumour-immune evolutionary trajectory of colorectal cancer.

Immune system control is a major hurdle that cancer evolution must circumvent. The relative timing and evolutionary dynamics of subclones that have escaped immune control remain incompletely characterized, and how immune-mediated selection shapes the epigenome has received little attention. Here, we infer the genome- and epigenome-driven evolutionary dynamics of tumour-immune coevolution within primary colorectal cancers (CRCs). We utilise our existing CRC multi-region multi-omic dataset that we supplement with high-resolution spatially-resolved neoantigen sequencing data and highly multiplexed imaging of the tumour microenvironment (TME). Analysis of somatic chromatin accessibility alterations (SCAAs) reveals frequent somatic loss of accessibility at antigen presenting genes, and that SCAAs contribute to silencing of neoantigens. We observe that strong immune escape and exclusion occur at the outset of CRC formation, and that within tumours, including at the microscopic level of individual tumour glands, additional immune escape alterations have negligible consequences for the immunophenotype of cancer cells. Further minor immuno-editing occurs during local invasion and is associated with TME reorganisation, but that evolutionary bottleneck is relatively weak. Collectively, we show that immune evasion in CRC follows a "Big Bang" evolutionary pattern, whereby genetic, epigenetic and TME-driven immune evasion acquired by the time of transformation defines subsequent cancer-immune evolution.