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.

Conumee 2.0: enhanced copy-number variation analysis from DNA methylation arrays for humans and mice.

MOTIVATION: Copy-number variations (CNVs) are common genetic alterations in cancer and their detection may impact tumor classification and therapeutic decisions. However, detection of clinically relevant large and focal CNVs remains challenging when sample material or resources are limited. This has motivated us to create a software tool to infer CNVs from DNA methylation arrays which are often generated as part of clinical routines and in research settings. RESULTS: We present our R package, conumee 2.0, that combines tangent normalization, an adjustable genomic binning heuristic, and weighted circular binary segmentation to utilize DNA methylation arrays for CNV analysis and mitigate technical biases and batch effects. Segmentation results were validated in a lung squamous cell carcinoma dataset from TCGA (n = 367 samples) by comparison to segmentations derived from genotyping arrays (Pearson's correlation coefficient of 0.91). We further introduce a segmented block bootstrapping approach to detect focal alternations that achieved 60.9% sensitivity and 98.6% specificity for deletions affecting CDKN2A/B (60.0% and 96.9% for RB1, respectively) in a low-grade glioma cohort from TCGA (n = 239 samples). Finally, our tool provides functionality to detect and summarize CNVs across large sample cohorts. AVAILABILITY AND IMPLEMENTATION: Conumee 2.0 is available under open-source license at: https://github.com/hovestadtlab/conumee2.

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.

Proteogenomics refines the molecular classification of chronic lymphocytic leukemia.

Cancer heterogeneity at the proteome level may explain differences in therapy response and prognosis beyond the currently established genomic and transcriptomic-based diagnostics. The relevance of proteomics for disease classifications remains to be established in clinically heterogeneous cancer entities such as chronic lymphocytic leukemia (CLL). Here, we characterize the proteome and transcriptome alongside genetic and ex-vivo drug response profiling in a clinically annotated CLL discovery cohort (n = 68). Unsupervised clustering of the proteome data reveals six subgroups. Five of these proteomic groups are associated with genetic features, while one group is only detectable at the proteome level. This new group is characterized by accelerated disease progression, high spliceosomal protein abundances associated with aberrant splicing, and low B cell receptor signaling protein abundances (ASB-CLL). Classifiers developed to identify ASB-CLL based on its characteristic proteome or splicing signature in two independent cohorts (n = 165, n = 169) confirm that ASB-CLL comprises about 20% of CLL patients. The inferior overall survival in ASB-CLL is also independent of both TP53- and IGHV mutation status. Our multi-omics analysis refines the classification of CLL and highlights the potential of proteomics to improve cancer patient stratification beyond genetic and transcriptomic profiling.

Comprehensive genomic and epigenomic analysis in cancer of unknown primary guides molecularly-informed therapies despite heterogeneity.

The benefit of molecularly-informed therapies in cancer of unknown primary (CUP) is unclear. Here, we use comprehensive molecular characterization by whole genome/exome, transcriptome and methylome analysis in 70 CUP patients to reveal substantial mutational heterogeneity with TP53, MUC16, KRAS, LRP1B and CSMD3 being the most frequently mutated known cancer-related genes. The most common fusion partner is FGFR2, the most common focal homozygous deletion affects CDKN2A. 56/70 (80%) patients receive genomics-based treatment recommendations which are applied in 20/56 (36%) cases. Transcriptome and methylome data provide evidence for the underlying entity in 62/70 (89%) cases. Germline analysis reveals five (likely) pathogenic mutations in five patients. Recommended off-label therapies translate into a mean PFS ratio of 3.6 with a median PFS1 of 2.9 months (17 patients) and a median PFS2 of 7.8 months (20 patients). Our data emphasize the clinical value of molecular analysis and underline the need for innovative, mechanism-based clinical trials.

Comparative parallel multi-omics analysis during the induction of pluripotent and trophectoderm states.

Following fertilization, it is only at the 32-64-cell stage when a clear segregation between cells of the inner cell mass and trophectoderm is observed, suggesting a 'T'-shaped model of specification. Here, we examine whether the acquisition of these two states in vitro, by nuclear reprogramming, share similar dynamics/trajectories. Using a comparative parallel multi-omics analysis (i.e., bulk RNA-seq, scRNA-seq, ATAC-seq, ChIP-seq, RRBS and CNVs) on cells undergoing reprogramming to pluripotency and TSC state we show that each reprogramming system exhibits specific trajectories from the onset of the process, suggesting 'V'-shaped model. We describe in detail the various trajectories toward the two states and illuminate reprogramming stage-specific markers, blockers, facilitators and TSC subpopulations. Finally, we show that while the acquisition of the TSC state involves the silencing of embryonic programs by DNA methylation, during the acquisition of pluripotency these regions are initially defined but retain inactive by the elimination of H3K27ac.

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.

Association of mutation signature effectuating processes with mutation hotspots in driver genes and non-coding regions.

Cancer driving mutations are difficult to identify especially in the non-coding part of the genome. Here, we present sigDriver, an algorithm dedicated to call driver mutations. Using 3813 whole-genome sequenced tumors from International Cancer Genome Consortium, The Cancer Genome Atlas Program, and a childhood pan-cancer cohort, we employ mutational signatures based on single-base substitution in the context of tri- and penta-nucleotide motifs for hotspot discovery. Knowledge-based annotations on mutational hotspots reveal enrichment in coding regions and regulatory elements for 6 mutational signatures, including APOBEC and somatic hypermutation signatures. APOBEC activity is associated with 32 hotspots of which 11 are known and 11 are putative regulatory drivers. Somatic single nucleotide variants clusters detected at hypermutation-associated hotspots are distinct from translocation or gene amplifications. Patients carrying APOBEC induced PIK3CA driver mutations show lower occurrence of signature SBS39. In summary, sigDriver uncovers mutational processes associated with known and putative tumor drivers and hotspots particularly in the non-coding regions of the genome.

Longitudinal analyses of CLL in mice identify leukemia-related clonal changes including a Myc gain predicting poor outcome in patients.

Chronic lymphocytic leukemia (CLL) is a B-cell malignancy mainly occurring at an advanced age with no single major genetic driver. Transgenic expression of TCL1 in B cells leads after a long latency to a CLL-like disease in aged Eµ-TCL1 mice suggesting that TCL1 overexpression is not sufficient for full leukemic transformation. In search for secondary genetic events and to elucidate the clonal evolution of CLL, we performed whole exome and B-cell receptor sequencing of longitudinal leukemia samples of Eµ-TCL1 mice. We observed a B-cell receptor stereotypy, as described in patients, confirming that CLL is an antigen-driven disease. Deep sequencing showed that leukemia in Eµ-TCL1 mice is mostly monoclonal. Rare oligoclonality was associated with inability of tumors to develop disease upon adoptive transfer in mice. In addition, we identified clonal changes and a sequential acquisition of mutations with known relevance in CLL, which highlights the genetic similarities and therefore, suitability of the Eµ-TCL1 mouse model for progressive CLL. Among them, a recurrent gain of chromosome 15, where Myc is located, was identified in almost all tumors in Eµ-TCL1 mice. Interestingly, amplification of 8q24, the chromosomal region containing MYC in humans, was associated with worse outcome of patients with CLL.

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.

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.

CATCH: A Prospective Precision Oncology Trial in Metastatic Breast Cancer.

PURPOSE: CATCH (Comprehensive Assessment of clinical feaTures and biomarkers to identify patients with advanced or metastatic breast Cancer for marker driven trials in Humans) is a prospective precision oncology program that uses genomics and transcriptomics to guide therapeutic decisions in the clinical management of metastatic breast cancer. Herein, we report our single-center experience and results on the basis of the first 200 enrolled patients of an ongoing trial. METHODS: From June 2017 to March 2019, 200 patients who had either primary metastatic or progressive disease, with any number of previous treatment lines and at least one metastatic site accessible to biopsy, were enrolled. DNA and RNA from tumor tissue and corresponding blood-derived nontumor DNA were profiled using whole-genome and transcriptome sequencing. Identified actionable alterations were brought into clinical context in a multidisciplinary molecular tumor board (MTB) with the aim of prioritizing personalized treatment recommendations. RESULTS: Among the first 200 enrolled patients, 128 (64%) were discussed in the MTB, of which 64 (50%) were subsequently treated according to MTB recommendation. Of 53 evaluable patients, 21 (40%) achieved either stable disease (n = 13, 25%) or partial response (n = 8, 15%). Furthermore, 16 (30%) of those patients showed improvement in progression-free survival of at least 30% while on MTB-recommended treatment compared with the progression-free survival of the previous treatment line. CONCLUSION: The initial phase of this study demonstrates that precision oncology on the basis of whole-genome and RNA sequencing is feasible when applied in the clinical management of patients with metastatic breast cancer and provides clinical benefit to a substantial proportion of patients.

Carbon ion radiotherapy eradicates medulloblastomas with chromothripsis in an orthotopic Li-Fraumeni patient-derived mouse model.

BACKGROUND: Medulloblastomas with chromothripsis developing in children with Li-Fraumeni Syndrome (germline TP53 mutations) are highly aggressive brain tumors with dismal prognosis. Conventional photon radiotherapy and DNA-damaging chemotherapy are not successful for these patients and raise the risk of secondary malignancies. We hypothesized that the pronounced homologous recombination deficiency in these tumors might offer vulnerabilities that can be therapeutically utilized in combination with high linear energy transfer carbon ion radiotherapy. METHODS: We tested high-precision particle therapy with carbon ions and protons as well as topotecan with or without PARP inhibitor in orthotopic primary and matched relapsed patient-derived xenograft models. Tumor and normal tissue underwent longitudinal morphological MRI, cellular (markers of neurogenesis and DNA damage-repair), and molecular characterization (whole-genome sequencing). RESULTS: In the primary medulloblastoma model, carbon ions led to complete response in 79% of animals irrespective of PARP inhibitor within a follow-up period of 300 days postirradiation, as detected by MRI and histology. No sign of neurologic symptoms, impairment of neurogenesis or in-field carcinogenesis was detected in repair-deficient host mice. PARP inhibitors further enhanced the effect of proton irradiation. In the postradiotherapy relapsed tumor model, median survival was significantly increased after carbon ions (96 days) versus control (43 days, P < .0001). No major change in the clonal composition was detected in the relapsed model. CONCLUSION: The high efficacy and favorable toxicity profile of carbon ions warrants further investigation in primary medulloblastomas with chromothripsis. Postradiotherapy relapsed medulloblastomas exhibit relative resistance compared to treatment-naïve tumors, calling for exploration of multimodal strategies.

The age of adult pilocytic astrocytoma cells.

Adult pilocytic astrocytomas (PAs) have been regarded as indistinguishable from pediatric PAs in terms of genome-wide expression and methylation patterns. It has been unclear whether adult PAs arise early in life and remain asymptomatic until adulthood, or whether they develop during adulthood. We sought to determine the age and origin of adult human PAs using two types of "marks" in the genomic DNA. First, we analyzed the DNA methylation patterns of adult and pediatric PAs to distinguish between PAs of different anatomic locations (n = 257 PA and control brain tissues). Second, we measured the concentration of nuclear bomb test-derived 14C in genomic DNA (n = 14 cases), which indicates the time point of the formation of human cell populations. Our data suggest that adult and pediatric PAs developing in the infratentorial brain are closely related and potentially develop from precursor cells early in life, whereas supratentorial PAs might show age and location-specific differences.

Unraveling most abundant mutational signatures in head and neck cancer.

Genomic alterations are a driving force in the multistep process of head and neck cancer (HNC) and result from the interaction of exogenous environmental exposures and endogenous cellular processes. Each of these processes leaves a characteristic pattern of mutations on the tumor genome providing the unique opportunity to decipher specific signatures of mutational processes operative during HNC pathogenesis and to address their prognostic value. Computational analysis of whole exome sequencing data of the HIPO-HNC (Heidelberg Center for Personalized Oncology-head and neck cancer) (n = 83) and TCGA-HNSC (The Cancer Genome Atlas-Head and Neck Squamous Cell Carcinoma) (n = 506) cohorts revealed five common mutational signatures (Catalogue of Somatic Mutations in Cancer [COSMIC] Signatures 1, 2, 3, 13 and 16) and demonstrated their significant association with etiological risk factors (tobacco, alcohol and HPV16). Unsupervised hierarchical clustering identified four clusters (A, B, C1 and C2) of which Subcluster C2 was enriched for cases with a higher frequency of signature 16 mutations. Tumors of Subcluster C2 had significantly lower p16INK4A expression accompanied by homozygous CDKN2A deletion in almost one half of cases. Survival analysis revealed an unfavorable prognosis for patients with tumors characterized by a higher mutation burden attributed to signature 16 as well as cases in Subcluster C2. Finally, a LASSO-Cox regression model was applied to prioritize clinically relevant signatures and to establish a prognostic risk score for head and neck squamous cell carcinoma patients. In conclusion, our study provides a proof of concept that computational analysis of somatic mutational signatures is not only a powerful tool to decipher environmental and intrinsic processes in the pathogenesis of HNC, but could also pave the way to establish reliable prognostic patterns.

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.

Chromothripsis in Human Breast Cancer.

Chromothripsis is a form of genome instability by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosomes. Widely assumed to be an early event in tumor development, this phenomenon plays a prominent role in tumor onset. In this study, an analysis of chromothripsis in 252 human breast cancers from two patient cohorts (149 metastatic breast cancers, 63 untreated primary tumors, 29 local relapses, and 11 longitudinal pairs) using whole-genome and whole-exome sequencing reveals that chromothripsis affects a substantial proportion of human breast cancers, with a prevalence over 60% in a cohort of metastatic cases and 25% in a cohort comprising predominantly luminal breast cancers. In the vast majority of cases, multiple chromosomes per tumor were affected, with most chromothriptic events on chromosomes 11 and 17 including, among other significantly altered drivers, CCND1, ERBB2, CDK12, and BRCA1. Importantly, chromothripsis generated recurrent fusions that drove tumor development. Chromothripsis-related rearrangements were linked with univocal mutational signatures, with clusters of point mutations due to kataegis in close proximity to the genomic breakpoints and with the activation of specific signaling pathways. Analyzing the temporal order of events in tumors with and without chromothripsis as well as longitudinal analysis of chromothriptic patterns in tumor pairs offered important insights into the role of chromothriptic chromosomes in tumor evolution. SIGNIFICANCE: These findings identify chromothripsis as a major driving event in human breast cancer.