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.

Low coverage whole genome sequencing of low-grade dysplasia strongly predicts colorectal cancer risk in ulcerative colitis.

Patients with inflammatory bowel disease (IBD) are at increased risk of colorectal cancer (CRC), and this risk increases dramatically in those who develop low-grade dysplasia (LGD). However, there is currently no accurate way to risk-stratify patients with LGD, leading to both over- and under-treatment of cancer risk. Here we show that the burden of somatic copy number alterations (CNAs) within resected LGD lesions strongly predicts future cancer development. We performed a retrospective multi-centre validated case-control study of n=122 patients (40 progressors, 82 non-progressors, 270 LGD regions). Low coverage whole genome sequencing revealed CNA burden was significantly higher in progressors than non-progressors (p=2×10-6 in discovery cohort) and was a very significant predictor of CRC risk in univariate analysis (odds ratio = 36; p=9×10-7), outperforming existing clinical risk factors such as lesion size, shape and focality. Optimal risk prediction was achieved with a multivariate model combining CNA burden with the known clinical risk factor of incomplete LGD resection. The measurement of CNAs in LGD lesions is a robust, low-cost and rapidly translatable predictor of CRC risk in IBD that can be used to direct management and so prevent CRC in high-risk individuals whilst sparing those at low-risk from unnecessary intervention.

Pre-B cell receptor-mediated activation of BCL6 induces pre-B cell quiescence through transcriptional repression of MYC.

Initial cell surface expression of the pre-B cell receptor induces proliferation. After 2 to 5 divisions, however, large pre-BII (Fraction C') cells exit cell cycle to become resting, small pre-BII cells (Fraction D). The mechanism by which pre-BII cells exit cell cycle, however, is currently unclear. The checkpoint at the Fraction C'-D transition is critical for immunoglobulin light chain gene recombination and to prevent malignant transformation into acute lymphoblastic leukemia. Here we demonstrate that inducible activation of pre-B cell receptor signaling induces cell-cycle exit through up-regulation of the transcriptional repressor BCL6. Inducible activation of BCL6 downstream of the pre-B cell receptor results in transcriptional repression of MYC and CCND2. Hence, pre-B cell receptor-mediated activation of BCL6 limits pre-B cell proliferation and induces cellular quiescence at the small pre-BII (Fraction D) stage.

SNP array analysis of acute promyelocytic leukemia may be of prognostic relevance and identifies a potential high risk group with recurrent deletions on chromosomal subband 1q31.3.

To search for new copy number alterations (CNAs) in acute promyelocytic leukemia (APL), we analyzed DNA from leukemic blasts of 93 acute promyelocytic leukemia (APL) patients with Genome-Wide SNP 6.0 arrays (SNP-A). We identified 259 CNAs consisting of 170 heterozygous deletions, 82 amplifications, and 7 regions of copy number neutral loss of heterozygosity. One of the most common CNAs was a deletion on chromosomal subband 1q31.3 in 13 of 93 (14%) patients encompassing the coding regions for the microRNAs mir181a1/b1. In multivariable analysis with the covariates age, white blood cell count, platelet count, and FLT3-ITD/FLT3 D835 mutations we found that after adjustment for patients' age (P<0.0001), patients with 2 or more CNAs detected by SNP-A had a higher risk of death (hazard ratio=5.942, P=0.0015) than patients with 0 or 1 CNA. Deletions of 1q31.3 were associated with a higher number of CNAs (median 2 vs. 8, P<0.0001) and were a strong independent prognostic factor for an increased risk of relapse (hazard ratio=28.9, P=0.0031). This study presents a comprehensive assessment of new CNAs as pathomechanistically relevant targets and possible prognostic factors which could refine risk stratification of APL.

Adaptive therapy achieves long-term control of chemotherapy resistance in high grade ovarian cancer.

Drug resistance results in poor outcomes for most patients with metastatic cancer. Adaptive Therapy (AT) proposes to address this by exploiting presumed fitness costs incurred by drug-resistant cells when drug is absent, and prescribing dose reductions to allow fitter, sensitive cells to re-grow and re-sensitise the tumour. However, empirical evidence for treatment-induced fitness change is lacking. We show that fitness costs in chemotherapy-resistant ovarian cancer cause selective decline and apoptosis of resistant populations in low-resource conditions. Moreover, carboplatin AT caused fluctuations in sensitive/resistant tumour population size in vitro and significantly extended survival of tumour-bearing mice. In sequential blood-derived cell-free DNA and tumour samples obtained longitudinally from ovarian cancer patients during treatment, we inferred resistant cancer cell population size through therapy and observed it correlated strongly with disease burden. These data have enabled us to launch a multicentre, phase 2 randomised controlled trial (ACTOv) to evaluate AT in ovarian cancer.

Long-term control of HIV by CCR5 Delta32/Delta32 stem-cell transplantation.

Infection with the human immunodeficiency virus type 1 (HIV-1) requires the presence of a CD4 receptor and a chemokine receptor, principally chemokine receptor 5 (CCR5). Homozygosity for a 32-bp deletion in the CCR5 allele provides resistance against HIV-1 acquisition. We transplanted stem cells from a donor who was homozygous for CCR5 delta32 in a patient with acute myeloid leukemia and HIV-1 infection. The patient remained without viral rebound 20 months after transplantation and discontinuation of antiretroviral therapy. This outcome demonstrates the critical role CCR5 plays in maintaining HIV-1 infection.

Expression of IGFBP7 in acute leukemia is regulated by DNA methylation.

The important role of insulin-like growth factor binding protein 7 (IGFBP7) as a tumor suppressor in solid tumors has been revealed in several studies. Interestingly, in a recent study IGFBP7 was also shown to be aberrantly expressed in acute leukemia. Moreover, in acute T-lymphoblastic leukemia (T-ALL), high IGFBP7 expression predicts primary therapy resistance. In order to elucidate the mechanisms underlying aberrant IGFBP7 expression, we used pyrosequencing technology to investigate the DNA methylation of IGFBP7 in 109 T-ALL patient samples. Aberrant methylation was shown and hypomethylation was associated with an early immunophenotype and co-expression of the stem cell markers CD117 (P < 0.001) and CD34 (P < 0.001). In concordance, gene expression profiles of 86 T-ALL patients revealed upregulation of stem cell markers (CD34 and CD133) as well as genes associated with poor outcome and pathogenesis of leukemia (MN1, BAALC, FLT3) in the high IGFBP7 expression group. In conclusion, aberrant IGFBP7 expression is regulated by DNA methylation in acute leukemia. Hypomethylation of the gene is likely to characterize an immature and a more malignant subtype of the disease.

The role of single-cell sequencing in studying tumour evolution.

Tumour evolution is a complex interplay between the acquisition of somatic (epi)genomic changes in tumour cells and the phenotypic consequences they cause, all in the context of a changing microenvironment. Single-cell sequencing offers a window into this dynamic process at the ultimate resolution of individual cells. In this review, we discuss the transformative insight offered by single-cell sequencing technologies for understanding tumour evolution.

LiquidCNA: Tracking subclonal evolution from longitudinal liquid biopsies using somatic copy number alterations.

Cell-free DNA (cfDNA) measured via liquid biopsies provides a way for minimally invasive monitoring of tumor evolutionary dynamics during therapy. Here we present liquidCNA, a method to track subclonal evolution from longitudinally collected cfDNA samples sequenced through cost-effective low-pass whole-genome sequencing. LiquidCNA utilizes somatic copy number alteration (SCNA) to simultaneously genotype and quantify the size of the dominant subclone without requiring B-allele frequency information, matched-normal samples, or prior knowledge on the genetic identity of the emerging clone. We demonstrate the accuracy of liquidCNA in synthetically generated sample sets and in vitro mixtures of cancer cell lines. In vivo application in patients with metastatic lung cancer reveals the progressive emergence of a novel tumor subpopulation. LiquidCNA is straightforward to use, is computationally inexpensive, and enables continuous monitoring of subclonal evolution to understand and control-therapy-induced resistance.

Bone marrow derived stromal cells from myelodysplastic syndromes are altered but not clonally mutated in vivo.

The bone marrow (BM) stroma in myeloid neoplasms is altered and it is hypothesized that this cell compartment may also harbor clonal somatically acquired mutations. By exome sequencing of in vitro expanded mesenchymal stromal cells (MSCs) from n = 98 patients with myelodysplastic syndrome (MDS) and n = 28 healthy controls we show that these cells accumulate recurrent mutations in genes such as ZFX (n = 8/98), RANK (n = 5/98), and others. MDS derived MSCs display higher mutational burdens, increased replicative stress, senescence, inflammatory gene expression, and distinct mutational signatures as compared to healthy MSCs. However, validation experiments in serial culture passages, chronological BM aspirations and backtracking of high confidence mutations by re-sequencing primary sorted MDS MSCs indicate that the discovered mutations are secondary to in vitro expansion but not present in primary BM. Thus, we here report that there is no evidence for clonal mutations in the BM stroma of MDS patients.

Reconstructing single-cell karyotype alterations in colorectal cancer identifies punctuated and gradual diversification patterns.

Central to tumor evolution is the generation of genetic diversity. However, the extent and patterns by which de novo karyotype alterations emerge and propagate within human tumors are not well understood, especially at single-cell resolution. Here, we present 3D Live-Seq-a protocol that integrates live-cell imaging of tumor organoid outgrowth and whole-genome sequencing of each imaged cell to reconstruct evolving tumor cell karyotypes across consecutive cell generations. Using patient-derived colorectal cancer organoids and fresh tumor biopsies, we demonstrate that karyotype alterations of varying complexity are prevalent and can arise within a few cell generations. Sub-chromosomal acentric fragments were prone to replication and collective missegregation across consecutive cell divisions. In contrast, gross genome-wide karyotype alterations were generated in a single erroneous cell division, providing support that aneuploid tumor genomes can evolve via punctuated evolution. Mapping the temporal dynamics and patterns of karyotype diversification in cancer enables reconstructions of evolutionary paths to malignant fitness.