Interleukin-10 receptor signaling promotes the maintenance of a PD-1int TCF-1+ CD8+ T cell population that sustains anti-tumor immunity.

T cell exhaustion limits anti-tumor immunity and responses to immunotherapy. Here, we explored the microenvironmental signals regulating T cell exhaustion using a model of chronic lymphocytic leukemia (CLL). Single-cell analyses identified a subset of PD-1hi, functionally impaired CD8+ T cells that accumulated in secondary lymphoid organs during disease progression and a functionally competent PD-1int subset. Frequencies of PD-1int TCF-1+ CD8+ T cells decreased upon Il10rb or Stat3 deletion, leading to accumulation of PD-1hi cells and accelerated tumor progression. Mechanistically, inhibition of IL-10R signaling altered chromatin accessibility and disrupted cooperativity between the transcription factors NFAT and AP-1, promoting a distinct NFAT-associated program. Low IL10 expression or loss of IL-10R-STAT3 signaling correlated with increased frequencies of exhausted CD8+ T cells and poor survival in CLL and in breast cancer patients. Thus, balance between PD-1hi, exhausted CD8+ T cells and functional PD-1int TCF-1+ CD8+ T cells is regulated by cell-intrinsic IL-10R signaling, with implications for immunotherapy.

Hypermutation of the inactive X chromosome is a frequent event in cancer.

Mutation is a fundamental process in tumorigenesis. However, the degree to which the rate of somatic mutation varies across the human genome and the mechanistic basis underlying this variation remain to be fully elucidated. Here, we performed a cross-cancer comparison of 402 whole genomes comprising a diverse set of childhood and adult tumors, including both solid and hematopoietic malignancies. Surprisingly, we found that the inactive X chromosome of many female cancer genomes accumulates on average twice and up to four times as many somatic mutations per megabase, as compared to the individual autosomes. Whole-genome sequencing of clonally expanded hematopoietic stem/progenitor cells (HSPCs) from healthy individuals and a premalignant myelodysplastic syndrome (MDS) sample revealed no X chromosome hypermutation. Our data suggest that hypermutation of the inactive X chromosome is an early and frequent feature of tumorigenesis resulting from DNA replication stress in aberrantly proliferating cells.

Genome sequencing of pediatric medulloblastoma links catastrophic DNA rearrangements with TP53 mutations.

Genomic rearrangements are thought to occur progressively during tumor development. Recent findings, however, suggest an alternative mechanism, involving massive chromosome rearrangements in a one-step catastrophic event termed chromothripsis. We report the whole-genome sequencing-based analysis of a Sonic-Hedgehog medulloblastoma (SHH-MB) brain tumor from a patient with a germline TP53 mutation (Li-Fraumeni syndrome), uncovering massive, complex chromosome rearrangements. Integrating TP53 status with microarray and deep sequencing-based DNA rearrangement data in additional patients reveals a striking association between TP53 mutation and chromothripsis in SHH-MBs. Analysis of additional tumor entities substantiates a link between TP53 mutation and chromothripsis, and indicates a context-specific role for p53 in catastrophic DNA rearrangements. Among these, we observed a strong association between somatic TP53 mutations and chromothripsis in acute myeloid leukemia. These findings connect p53 status and chromothripsis in specific tumor types, providing a genetic basis for understanding particularly aggressive subtypes of cancer.

T-bet suppresses proliferation of malignant B cells in chronic lymphocytic leukemia.

The T-box transcription factor T-bet is known as a master regulator of T-cell response but its role in malignant B cells is not sufficiently explored. Here, we conducted single-cell resolved multi-omics analyses of malignant B cells from patients with chronic lymphocytic leukemia (CLL) and studied a CLL mouse model with genetic knockout of TBX21. We found that T-bet acts as a tumor suppressor in malignant B cells by decreasing their proliferation rate. NF-κB activity induced by inflammatory signals provided by the microenvironment, triggered T-bet expression which impacted on promoter proximal and distal chromatin co-accessibility and controlled a specific gene signature by mainly suppressing transcription. Gene set enrichment analysis identified a positive regulation of interferon signaling, and a negative control of proliferation by T-bet. In line, we showed that T-bet represses cell cycling and is associated with longer overall survival of CLL patients. Our study uncovers a novel tumor suppressive role of T-bet in malignant B cells via its regulation of inflammatory processes and cell cycling which has implications for stratification and therapy of CLL patients. Linking T-bet activity to inflammation explains the good prognostic role of genetic alterations in inflammatory signaling pathways in CLL.

LifeTime and improving European healthcare through cell-based interceptive medicine.

Here we describe the LifeTime Initiative, which aims to track, understand and target human cells during the onset and progression of complex diseases, and to analyse their response to therapy at single-cell resolution. This mission will be implemented through the development, integration and application of single-cell multi-omics and imaging, artificial intelligence and patient-derived experimental disease models during the progression from health to disease. The analysis of large molecular and clinical datasets will identify molecular mechanisms, create predictive computational models of disease progression, and reveal new drug targets and therapies. The timely detection and interception of disease embedded in an ethical and patient-centred vision will be achieved through interactions across academia, hospitals, patient associations, health data management systems and industry. The application of this strategy to key medical challenges in cancer, neurological and neuropsychiatric disorders, and infectious, chronic inflammatory and cardiovascular diseases at the single-cell level will usher in cell-based interceptive medicine in Europe over the next decade.

Time-resolved, integrated analysis of clonally evolving genomes.

Clonal genome evolution is a key feature of asexually reproducing species and human cancer development. While many studies have described the landscapes of clonal genome evolution in cancer, few determine the underlying evolutionary parameters from molecular data, and even fewer integrate theory with data. We derived theoretical results linking mutation rate, time, expansion dynamics, and biological/clinical parameters. Subsequently, we inferred time-resolved estimates of evolutionary parameters from mutation accumulation, mutational signatures and selection. We then applied this framework to predict the time of speciation of the marbled crayfish, an enigmatic, globally invasive parthenogenetic freshwater crayfish. The results predict that speciation occurred between 1986 and 1990, which is consistent with biological records. We also used our framework to analyze whole-genome sequencing datasets from primary and relapsed glioblastoma, an aggressive brain tumor. The results identified evolutionary subgroups and showed that tumor cell survival could be inferred from genomic data that was generated during the resection of the primary tumor. In conclusion, our framework allowed a time-resolved, integrated analysis of key parameters in clonally evolving genomes, and provided novel insights into the evolutionary age of marbled crayfish and the progression of glioblastoma.

Longitudinal molecular trajectories of diffuse glioma in adults.

The evolutionary processes that drive universal therapeutic resistance in adult patients with diffuse glioma remain unclear1,2. Here we analysed temporally separated DNA-sequencing data and matched clinical annotation from 222 adult patients with glioma. By analysing mutations and copy numbers across the three major subtypes of diffuse glioma, we found that driver genes detected at the initial stage of disease were retained at recurrence, whereas there was little evidence of recurrence-specific gene alterations. Treatment with alkylating agents resulted in a hypermutator phenotype at different rates across the glioma subtypes, and hypermutation was not associated with differences in overall survival. Acquired aneuploidy was frequently detected in recurrent gliomas and was characterized by IDH mutation but without co-deletion of chromosome arms 1p/19q, and further converged with acquired alterations in the cell cycle and poor outcomes. The clonal architecture of each tumour remained similar over time, but the presence of subclonal selection was associated with decreased survival. Finally, there were no differences in the levels of immunoediting between initial and recurrent gliomas. Collectively, our results suggest that the strongest selective pressures occur during early glioma development and that current therapies shape this evolution in a largely stochastic manner.

Author Correction: BCAT1 restricts αKG levels in AML stem cells leading to IDHmut-like DNA hypermethylation.

In Extended Data Fig. 1a of this Letter, the flow cytometry plot depicting the surface phenotype of AML sample DD08 was a duplicate of the plot for AML sample DD06. Supplementary Data 4 has been added to the Supplementary Information of the original Letter to clarify the proteome data acquisition and presentation. The original Letter has been corrected online.

scSNPdemux: a sensitive demultiplexing pipeline using single nucleotide polymorphisms for improved pooled single-cell RNA sequencing analysis.

BACKGROUND: Here we present scSNPdemux, a sample demultiplexing pipeline for single-cell RNA sequencing data using natural genetic variations in humans. The pipeline requires alignment files from Cell Ranger (10× Genomics), a population SNP database and genotyped single nucleotide polymorphisms (SNPs) per sample. The tool works on sparse genotyping data in VCF format for sample identification. RESULTS: The pipeline was tested on both single-cell and single-nuclei based RNA sequencing datasets and showed superior demultiplexing performance over the lipid-based CellPlex and Multi-seq sample multiplexing technique which incurs additional single cell library preparation steps. Specifically, our pipeline demonstrated superior sensitivity and specificity in cell-identity assignment over CellPlex, especially on immune cell types with low RNA content. CONCLUSIONS: We designed a streamlined pipeline for single-cell sample demultiplexing, aiming to overcome common problems in multiplexing samples using single cell libraries which might affect data quality and can be costly.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

BCAT1 restricts αKG levels in AML stem cells leading to IDHmut-like DNA hypermethylation.

The branched-chain amino acid (BCAA) pathway and high levels of BCAA transaminase 1 (BCAT1) have recently been associated with aggressiveness in several cancer entities. However, the mechanistic role of BCAT1 in this process remains largely uncertain. Here, by performing high-resolution proteomic analysis of human acute myeloid leukaemia (AML) stem-cell and non-stem-cell populations, we find the BCAA pathway enriched and BCAT1 protein and transcripts overexpressed in leukaemia stem cells. We show that BCAT1, which transfers α-amino groups from BCAAs to α-ketoglutarate (αKG), is a critical regulator of intracellular αKG homeostasis. Further to its role in the tricarboxylic acid cycle, αKG is an essential cofactor for αKG-dependent dioxygenases such as Egl-9 family hypoxia inducible factor 1 (EGLN1) and the ten-eleven translocation (TET) family of DNA demethylases. Knockdown of BCAT1 in leukaemia cells caused accumulation of αKG, leading to EGLN1-mediated HIF1α protein degradation. This resulted in a growth and survival defect and abrogated leukaemia-initiating potential. By contrast, overexpression of BCAT1 in leukaemia cells decreased intracellular αKG levels and caused DNA hypermethylation through altered TET activity. AML with high levels of BCAT1 (BCAT1high) displayed a DNA hypermethylation phenotype similar to cases carrying a mutant isocitrate dehydrogenase (IDHmut), in which TET2 is inhibited by the oncometabolite 2-hydroxyglutarate. High levels of BCAT1 strongly correlate with shorter overall survival in IDHWTTET2WT, but not IDHmut or TET2mut AML. Gene sets characteristic for IDHmut AML were enriched in samples from patients with an IDHWTTET2WTBCAT1high status. BCAT1high AML showed robust enrichment for leukaemia stem-cell signatures, and paired sample analysis showed a significant increase in BCAT1 levels upon disease relapse. In summary, by limiting intracellular αKG, BCAT1 links BCAA catabolism to HIF1α stability and regulation of the epigenomic landscape, mimicking the effects of IDH mutations. Our results suggest the BCAA-BCAT1-αKG pathway as a therapeutic target to compromise leukaemia stem-cell function in patients with IDHWTTET2WT AML.

Evidence for multi-copy Mega-NUMTs in the human genome.

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.

Modeling the B-cell receptor signaling on single cell level reveals a stable network circuit topology between nonmalignant B cells and chronic lymphocytic leukemia cells and between untreated cells and cells treated with kinase inhibitors.

B-cell receptor (BCR) signaling is central for the pathomechanism of chronic lymphocytic leukemia (CLL), and inhibitors of BCR signaling have substantially improved treatment options. To model malignant and nonmalignant BCR signaling, we quantified five components of BCR signaling (ZAP70/SYK, BTK, PLCγ2, AKT, ERK1/2) in single cells from primary human leukemic cells and from nonmalignant tissue. We measured signaling activity in a time-resolved manner after stimulation with BCR crosslinking by anti-IgM and/or anti-CD19 and with or without inhibition of phosphatases with H2 O2 . The phosphorylation of BCR signaling components was increased in malignant cells compared to nonmalignant cells and in IGHV unmutated CLL cells compared to IGHV mutated CLL cells. Intriguingly, inhibition of phosphatases with H2 O2 led to higher phosphorylation levels of BCR components in CLL cells with mutated IGHV compared to unmutated IGHV. We modeled the connectivity of the cascade components by correlating signal intensities across single cells. The network topology remained stable between malignant and nonmalignant cells. To additionally test for the impact of therapeutic compounds on the network topology, we challenged the BCR signaling cascade with inhibitors for BTK (ibrutinib), PI3K (idelalisib), LYN (dasatinib) and SYK (entospletinib). Idelalisib treatment resulted in similar effects in malignant and nonmalignant cells, whereas ibrutinib was mostly active on CLL cells. Idelalisib and ibrutinib had complementary effects on the BCR signaling cascade whose activity was further reduced upon dasatinib and entospletinib treatment. The characterization of the molecular circuitry of leukemic BCR signaling will allow a more refined targeting of this Achilles heel.

A synergistic interaction between HDAC- and PARP inhibitors in childhood tumors with chromothripsis.

Chromothripsis is a form of genomic instability characterized by the occurrence of tens to hundreds of clustered DNA double-strand breaks in a one-off catastrophic event. Rearrangements associated with chromothripsis are detectable in numerous tumor entities and linked with poor prognosis in some of these, such as Sonic Hedgehog medulloblastoma, neuroblastoma and osteosarcoma. Hence, there is a need for therapeutic strategies eliminating tumor cells with chromothripsis. Defects in DNA double-strand break repair, and in particular homologous recombination repair, have been linked with chromothripsis. Targeting DNA repair deficiencies by synthetic lethality approaches, we performed a synergy screen using drug libraries (n = 375 compounds, 15 models) combined with either a PARP inhibitor or cisplatin. This revealed a synergistic interaction between the HDAC inhibitor romidepsin and PARP inhibition. Functional assays, transcriptome analyses and in vivo validation in patient-derived xenograft mouse models confirmed the efficacy of the combinatorial treatment.

Clonal evolution in chronic lymphocytic leukemia is scant in relapsed but accelerated in refractory cases after chemo(immune) therapy.

Clonal evolution is involved in the progression of chronic lymphocytic leukemia (CLL). In order to link evolutionary patterns to different disease courses, we performed a long-term longitudinal mutation profiling study of CLL patients. Tracking somatic mutations and their changes in allele frequency over time and assessing the underlying cancer cell fraction revealed highly distinct evolutionary patterns. Surprisingly, in long-term stable disease and in relapse after long-lasting clinical response to treatment, clonal shifts are minor. In contrast, in refractory disease major clonal shifts occur although there is little impact on leukemia cell counts. As this striking pattern in refractory cases is not linked to a strong contribution of known CLL driver genes, the evolution is mostly driven by treatment-induced selection of sub-clones, underlining the need for novel, non-genotoxic treatment regimens.

Active medulloblastoma enhancers reveal subgroup-specific cellular origins.

Medulloblastoma is a highly malignant paediatric brain tumour, often inflicting devastating consequences on the developing child. Genomic studies have revealed four distinct molecular subgroups with divergent biology and clinical behaviour. An understanding of the regulatory circuitry governing the transcriptional landscapes of medulloblastoma subgroups, and how this relates to their respective developmental origins, is lacking. Here, using H3K27ac and BRD4 chromatin immunoprecipitation followed by sequencing (ChIP-seq) coupled with tissue-matched DNA methylation and transcriptome data, we describe the active cis-regulatory landscape across 28 primary medulloblastoma specimens. Analysis of differentially regulated enhancers and super-enhancers reinforced inter-subgroup heterogeneity and revealed novel, clinically relevant insights into medulloblastoma biology. Computational reconstruction of core regulatory circuitry identified a master set of transcription factors, validated by ChIP-seq, that is responsible for subgroup divergence, and implicates candidate cells of origin for Group 4. Our integrated analysis of enhancer elements in a large series of primary tumour samples reveals insights into cis-regulatory architecture, unrecognized dependencies, and cellular origins.

A versatile system to introduce clusters of genomic double-strand breaks in large cell populations.

In vitro assays for clustered DNA lesions will facilitate the analysis of the mechanisms underlying complex genome rearrangements such as chromothripsis, including the recruitment of repair factors to sites of DNA double-strand breaks (DSBs). We present a novel method generating localized DNA DSBs using UV irradiation with photomasks. The size of the damage foci and the spacing between lesions are fully adjustable, making the assay suitable for different cell types and targeted areas. We validated this setup with genomically stable epithelial cells, normal fibroblasts, pluripotent stem cells, and patient-derived primary cultures. Our method does not require a specialized device such as a laser, making it accessible to a broad range of users. Sensitization by 5-bromo-2-deoxyuridine incorporation is not required, which enables analyzing the DNA damage response in post-mitotic cells. Irradiated cells can be cultivated further, followed by time-lapse imaging or used for downstream biochemical analyses, thanks to the high throughput of the system. Importantly, we showed genome rearrangements in the irradiated cells, providing a proof of principle for the induction of structural variants by localized DNA lesions.

Combining ibrutinib and checkpoint blockade improves CD8+ T-cell function and control of chronic lymphocytic leukemia in Em-TCL1 mice.

Ibrutinib is a bruton's tyrosine kinase (BTK) inhibitor approved for the treatment of multiple B-cell malignancies, including chronic lymphocytic leukemia (CLL). In addition to blocking B-cell receptor signaling and chemokine receptor-mediated pathways in CLL cells, that are known drivers of disease, ibrutinib also affects the microenvironment in CLL via targeting BTK in myeloid cells and IL-2-inducible T-cell kinase (ITK) in T-cells. These non-BTK effects were suggested to contribute to the success of ibrutinib in CLL. By using the Eµ-TCL1 adoptive transfer mouse model of CLL, we observed that ibrutinib effectively controls leukemia development, but also results in significantly lower numbers of CD8+ effector T-cells, with lower expression of activation markers, as well as impaired proliferation and effector function. Using CD8+ T-cells from a T-cell receptor (TCR) reporter mouse, we verified that this is due to a direct effect of ibrutinib on TCR activity, and demonstrate that co-stimulation via CD28 overcomes these effects. Most interestingly, combination of ibrutinib with blocking antibodies targeting PD-1/PD-L1 axis in vivo improved CD8+ T-cell effector function and control of CLL. In sum, these data emphasize the strong immunomodulatory effects of ibrutinib and the therapeutic potential of its combination with immune checkpoint blockade in CLL.

TBET-expressing Th1 CD4+ T cells accumulate in chronic lymphocytic leukaemia without affecting disease progression in Eµ-TCL1 mice.

Chronic lymphocytic leukaemia (CLL) is associated with alterations in T cell number, subset distribution and function. Among these changes, an increase in CD4+ T cells was reported. CD4+ T cells are a heterogeneous population and distinct subsets have been described to exert pro- and anti-tumour functions. In CLL, controversial reports describing the dominance of IFNγ-expressing Th1 T cells or of IL-4-producing Th2 T cells exist. Our study shows that blood of CLL patients is enriched in Th1 T cells producing high amounts of IFNγ. Moreover, we observed that their frequency remains relatively stable in CLL patients over a time course of five years. Furthermore, we provide evidence for an accumulation of Th1 T cells in the Eµ-TCL1 mouse model of CLL. As TBET (encoded by Tbx21) is a crucial transcription factor for Th1 polarization, we generated Tbx21-/- bone marrow chimaeric mice which showed a lower number of IFNγ-producing Th1 T cells, and used them for adoptive transfer of Eµ-TCL1 leukaemia. Disease development in these mice was, however, comparable to that in wild-type controls, excluding a major role for TBET-expressing Th1 cells in Eµ-TCL1 leukaemia. Collectively, our data highlight that Th1 T cells accumulate in CLL but reducing their number has no impact on disease development.