Germline MBD4 deficiency causes a multi-tumor predisposition syndrome.

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management.

Uridine-cytidine kinase 2 potentiates the mutagenic influence of the antiviral ß-d-N4-hydroxycytidine.

Molnupiravir (EIDD-2801) is an antiviral that received approval for the treatment of severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) infection. Treatment of bacteria or cell lines with the active form of molnupiravir, ß-d-N4-hydroxycytidine (NHC, or EIDD-1931), induces mutations in DNA. Yet these results contrast in vivo genotoxicity studies conducted during registration of the drug. Using a CRISPR screen, we found that inactivating the pyrimidine salvage pathway component uridine-cytidine kinase 2 (Uck2) renders cells more tolerant of NHC. Short-term exposure to NHC increased the mutation rate in a mouse myeloid cell line, with most mutations being T:A to C:G transitions. Inactivating Uck2 impaired the mutagenic activity of NHC, whereas over-expression of Uck2 enhanced mutagenesis. UCK2 is upregulated in many cancers and cell lines. Our results suggest differences in ribonucleoside metabolism contribute to the variable mutagenicity of NHC observed in cancer cell lines and primary tissues.

Single-cell multiomics reveal the scale of multilayered adaptations enabling CLL relapse during venetoclax therapy.

Venetoclax (VEN) inhibits the prosurvival protein BCL2 to induce apoptosis and is a standard therapy for chronic lymphocytic leukemia (CLL), delivering high complete remission rates and prolonged progression-free survival in relapsed CLL but with eventual loss of efficacy. A spectrum of subclonal genetic changes associated with VEN resistance has now been described. To fully understand clinical resistance to VEN, we combined single-cell short- and long-read RNA-sequencing to reveal the previously unappreciated scale of genetic and epigenetic changes underpinning acquired VEN resistance. These appear to be multilayered. One layer comprises changes in the BCL2 family of apoptosis regulators, especially the prosurvival family members. This includes previously described mutations in BCL2 and amplification of the MCL1 gene but is heterogeneous across and within individual patient leukemias. Changes in the proapoptotic genes are notably uncommon, except for single cases with subclonal losses of BAX or NOXA. Much more prominent was universal MCL1 gene upregulation. This was driven by an overlying layer of emergent NF-κB (nuclear factor kappa B) activation, which persisted in circulating cells during VEN therapy. We discovered that MCL1 could be a direct transcriptional target of NF-κB. Both the switch to alternative prosurvival factors and NF-κB activation largely dissipate following VEN discontinuation. Our studies reveal the extent of plasticity of CLL cells in their ability to evade VEN-induced apoptosis. Importantly, these findings pinpoint new approaches to circumvent VEN resistance and provide a specific biological justification for the strategy of VEN discontinuation once a maximal response is achieved rather than maintaining long-term selective pressure with the drug.

Intact TP-53 function is essential for sustaining durable responses to BH3-mimetic drugs in leukemias.

Selective targeting of BCL-2 with the BH3-mimetic venetoclax has been a transformative treatment for patients with various leukemias. TP-53 controls apoptosis upstream of where BCL-2 and its prosurvival relatives, such as MCL-1, act. Therefore, targeting these prosurvival proteins could trigger apoptosis across diverse blood cancers, irrespective of TP53 mutation status. Indeed, targeting BCL-2 has produced clinically relevant responses in blood cancers with aberrant TP-53. However, in our study, TP53-mutated or -deficient myeloid and lymphoid leukemias outcompeted isogenic controls with intact TP-53, unless sufficient concentrations of BH3-mimetics targeting BCL-2 or MCL-1 were applied. Strikingly, tumor cells with TP-53 dysfunction escaped and thrived over time if inhibition of BCL-2 or MCL-1 was sublethal, in part because of an increased threshold for BAX/BAK activation in these cells. Our study revealed the key role of TP-53 in shaping long-term responses to BH3-mimetic drugs and reconciled the disparate pattern of initial clinical response to venetoclax, followed by subsequent treatment failure among patients with TP53-mutant chronic lymphocytic leukemia or acute myeloid leukemia. In contrast to BH3-mimetics targeting just BCL-2 or MCL-1 at doses that are individually sublethal, a combined BH3-mimetic approach targeting both prosurvival proteins enhanced lethality and durably suppressed the leukemia burden, regardless of TP53 mutation status. Our findings highlight the importance of using sufficiently lethal treatment strategies to maximize outcomes of patients with TP53-mutant disease. In addition, our findings caution against use of sublethal BH3-mimetic drug regimens that may enhance the risk of disease progression driven by emergent TP53-mutant clones.

Detecting copy number alterations in RNA-Seq using SuperFreq.

MOTIVATION: Calling copy number alterations (CNAs) from RNA sequencing (RNA-Seq) is challenging, because of the marked variability in coverage across genes and paucity of single nucleotide polymorphisms (SNPs). We have adapted SuperFreq to call absolute and allele sensitive CNAs from RNA-Seq. SuperFreq uses an error-propagation framework to combine and maximize information from read counts and B-allele frequencies. RESULTS: We used datasets from The Cancer Genome Atlas (TCGA) to assess the validity of CNA calls from RNA-Seq. When ploidy estimates were consistent, we found agreement with DNA SNP-arrays for over 98% of the genome for acute myeloid leukaemia (TCGA-AML, n = 116) and 87% for colorectal cancer (TCGA-CRC, n = 377). The sensitivity of CNA calling from RNA-Seq was dependent on gene density. Using RNA-Seq, SuperFreq detected 78% of CNA calls covering 100 or more genes with a precision of 94%. Recall dropped for focal events, but this also depended on signal intensity. For example, in the CRC cohort SuperFreq identified all cases (7/7) with high-level amplification of ERBB2, where the copy number was typically >20, but identified only 6% of cases (1/17) with moderate amplification of IGF2, which occurs over a smaller interval. SuperFreq offers an integrated platform for identification of CNAs and point mutations. As evidence of how SuperFreq can be applied, we used it to reproduce the established relationship between somatic mutation load and CNA profile in CRC using RNA-Seq alone. AVAILABILITY AND IMPLEMENTATION: SuperFreq is implemented in R and the code is available through GitHub: https://github.com/ChristofferFlensburg/SuperFreq/. Data and code to reproduce the figures are available at: https://gitlab.wehi.edu.au/flensburg.c/SuperFreq_RNA_paper. Data from TCGA (phs000178) was accessed from GDC following completion of a data access request through the database of Genotypes and Phenotypes (dbGaP). Data from the Leucegene consortium was downloaded from GEO (AML samples: GSE67040; normal CD34+ cells: GSE48846). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SuperFreq: Integrated mutation detection and clonal tracking in cancer.

Analysing multiple cancer samples from an individual patient can provide insight into the way the disease evolves. Monitoring the expansion and contraction of distinct clones helps to reveal the mutations that initiate the disease and those that drive progression. Existing approaches for clonal tracking from sequencing data typically require the user to combine multiple tools that are not purpose-built for this task. Furthermore, most methods require a matched normal (non-tumour) sample, which limits the scope of application. We developed SuperFreq, a cancer exome sequencing analysis pipeline that integrates identification of somatic single nucleotide variants (SNVs) and copy number alterations (CNAs) and clonal tracking for both. SuperFreq does not require a matched normal and instead relies on unrelated controls. When analysing multiple samples from a single patient, SuperFreq cross checks variant calls to improve clonal tracking, which helps to separate somatic from germline variants, and to resolve overlapping CNA calls. To demonstrate our software we analysed 304 cancer-normal exome samples across 33 cancer types in The Cancer Genome Atlas (TCGA) and evaluated the quality of the SNV and CNA calls. We simulated clonal evolution through in silico mixing of cancer and normal samples in known proportion. We found that SuperFreq identified 93% of clones with a cellular fraction of at least 50% and mutations were assigned to the correct clone with high recall and precision. In addition, SuperFreq maintained a similar level of performance for most aspects of the analysis when run without a matched normal. SuperFreq is highly versatile and can be applied in many different experimental settings for the analysis of exomes and other capture libraries. We demonstrate an application of SuperFreq to leukaemia patients with diagnosis and relapse samples.

Finding a suitable library size to call variants in RNA-Seq.

BACKGROUND: RNA sequencing allows the study of both gene expression changes and transcribed mutations, providing a highly effective way to gain insight into cancer biology. When planning the sequencing of a large cohort of samples, library size is a fundamental factor affecting both the overall cost and the quality of the results. Here we specifically address how overall library size influences the detection of somatic mutations in RNA-seq data in two acute myeloid leukaemia datasets. RESULTS : We simulated shallower sequencing depths by downsampling 45 acute myeloid leukaemia samples (100 bp PE) that are part of the Leucegene project, which were originally sequenced at high depth. We compared the sensitivity of six methods of recovering validated mutations on the same samples. The methods compared are a combination of three popular callers (MuTect, VarScan, and VarDict) and two filtering strategies. We observed an incremental loss in sensitivity when simulating libraries of 80M, 50M, 40M, 30M and 20M fragments, with the largest loss detected with less than 30M fragments (below 90%, average loss of 7%). The sensitivity in recovering insertions and deletions varied markedly between callers, with VarDict showing the highest sensitivity (60%). Single nucleotide variant sensitivity is relatively consistent across methods, apart from MuTect, whose default filters need adjustment when using RNA-Seq. We also analysed 136 RNA-Seq samples from the TCGA-LAML cohort (50 bp PE) and assessed the change in sensitivity between the initial libraries (average 59M fragments) and after downsampling to 40M fragments. When considering single nucleotide variants in recurrently mutated myeloid genes we found a comparable performance, with a 6% average loss in sensitivity using 40M fragments. CONCLUSIONS: Between 30M and 40M 100 bp PE reads are needed to recover 90-95% of the initial variants on recurrently mutated myeloid genes. To extend this result to another cancer type, an exploration of the characteristics of its mutations and gene expression patterns is suggested.

MBD4 guards against methylation damage and germ line deficiency predisposes to clonal hematopoiesis and early-onset AML.

The tendency of 5-methylcytosine (5mC) to undergo spontaneous deamination has had a major role in shaping the human genome, and this methylation damage remains the primary source of somatic mutations that accumulate with age. How 5mC deamination contributes to cancer risk in different tissues remains unclear. Genomic profiling of 3 early-onset acute myeloid leukemias (AMLs) identified germ line loss of MBD4 as an initiator of 5mC-dependent hypermutation. MBD4-deficient AMLs display a 33-fold higher mutation burden than AML generally, with >95% being C>T in the context of a CG dinucleotide. This distinctive signature was also observed in sporadic cancers that acquired biallelic mutations in MBD4 and in Mbd4 knockout mice. Sequential sampling of germ line cases demonstrated repeated expansion of blood cell progenitors with pathogenic mutations in DNMT3A, a key driver gene for both clonal hematopoiesis and AML. Our findings reveal genetic and epigenetic factors that shape the mutagenic influence of 5mC. Within blood cells, this links methylation damage to the driver landscape of clonal hematopoiesis and reveals a conserved path to leukemia. Germ line MBD4 deficiency enhances cancer susceptibility and predisposes to AML.

Correction: The Subclonal Architecture of Metastatic Breast Cancer: Results from a Prospective Community-Based Rapid Autopsy Program "CASCADE".

[This corrects the article DOI: 10.1371/journal.pmed.1002204.].

Jarid2 regulates hematopoietic stem cell function by acting with polycomb repressive complex 2.

Polycomb repressive complex 2 (PRC2) plays a key role in hematopoietic stem and progenitor cell (HSPC) function. Analyses of mouse mutants harboring deletions of core components have implicated PRC2 in fine-tuning multiple pathways that instruct HSPC behavior, yet how PRC2 is targeted to specific genomic loci within HSPCs remains unknown. Here we use short hairpin RNA-mediated knockdown to survey the function of PRC2 accessory factors that were defined in embryonic stem cells (ESCs) by testing the competitive reconstitution capacity of transduced murine HSPCs. We find that, similar to the phenotype observed upon depletion of core subunit Suz12, depleting Jarid2 enhances the competitive transplantation capacity of both fetal and adult mouse HSPCs. Furthermore, we demonstrate that depletion of JARID2 enhances the in vitro expansion and in vivo reconstitution capacity of human HSPCs. Gene expression profiling revealed common Suz12 and Jarid2 target genes that are enriched for the H3K27me3 mark established by PRC2. These data implicate Jarid2 as an important component of PRC2 that has a central role in coordinating HSPC function.

The Subclonal Architecture of Metastatic Breast Cancer: Results from a Prospective Community-Based Rapid Autopsy Program "CASCADE".

BACKGROUND: Understanding the cancer genome is seen as a key step in improving outcomes for cancer patients. Genomic assays are emerging as a possible avenue to personalised medicine in breast cancer. However, evolution of the cancer genome during the natural history of breast cancer is largely unknown, as is the profile of disease at death. We sought to study in detail these aspects of advanced breast cancers that have resulted in lethal disease. METHODS AND FINDINGS: Three patients with oestrogen-receptor (ER)-positive, human epidermal growth factor receptor 2 (HER2)-negative breast cancer and one patient with triple negative breast cancer underwent rapid autopsy as part of an institutional prospective community-based rapid autopsy program (CASCADE). Cases represented a range of management problems in breast cancer, including late relapse after early stage disease, de novo metastatic disease, discordant disease response, and disease refractory to treatment. Between 5 and 12 metastatic sites were collected at autopsy together with available primary tumours and longitudinal metastatic biopsies taken during life. Samples underwent paired tumour-normal whole exome sequencing and single nucleotide polymorphism (SNP) arrays. Subclonal architectures were inferred by jointly analysing all samples from each patient. Mutations were validated using high depth amplicon sequencing. Between cases, there were significant differences in mutational burden, driver mutations, mutational processes, and copy number variation. Within each case, we found dramatic heterogeneity in subclonal structure from primary to metastatic disease and between metastatic sites, such that no single lesion captured the breadth of disease. Metastatic cross-seeding was found in each case, and treatment drove subclonal diversification. Subclones displayed parallel evolution of treatment resistance in some cases and apparent augmentation of key oncogenic drivers as an alternative resistance mechanism. We also observed the role of mutational processes in subclonal evolution. Limitations of this study include the potential for bias introduced by joint analysis of formalin-fixed archival specimens with fresh specimens and the difficulties in resolving subclones with whole exome sequencing. Other alterations that could define subclones such as structural variants or epigenetic modifications were not assessed. CONCLUSIONS: This study highlights various mechanisms that shape the genome of metastatic breast cancer and the value of studying advanced disease in detail. Treatment drives significant genomic heterogeneity in breast cancers which has implications for disease monitoring and treatment selection in the personalised medicine paradigm.

Overexpression of Lmo2 initiates T-lymphoblastic leukemia via impaired thymocyte competition.

Cell competition has recently emerged as an important tumor suppressor mechanism in the thymus that inhibits autonomous thymic maintenance. Here, we show that the oncogenic transcription factor Lmo2 causes autonomous thymic maintenance in transgenic mice by inhibiting early T cell differentiation. This autonomous thymic maintenance results in the development of self-renewing preleukemic stem cells (pre-LSCs) and subsequent leukemogenesis, both of which are profoundly inhibited by restoration of thymic competition or expression of the antiapoptotic factor BCL2. Genomic analyses revealed the presence of Notch1 mutations in pre-LSCs before subsequent loss of tumor suppressors promotes the transition to overt leukemogenesis. These studies demonstrate a critical role for impaired cell competition in the development of pre-LSCs in a transgenic mouse model of T cell acute lymphoblastic leukemia (T-ALL), implying that this process plays a role in the ontogeny of human T-ALL.

Ganciclovir-induced mutations are present in a diverse spectrum of post-transplant malignancies.

BACKGROUND: Ganciclovir (GCV) is widely used in solid organ and haematopoietic stem cell transplant patients for prophylaxis and treatment of cytomegalovirus. It has long been considered a mutagen and carcinogen. However, the contribution of GCV to cancer incidence and other factors that influence its mutagenicity remains unknown. METHODS: This retrospective cohort study analysed genomics data for 121,771 patients who had undergone targeted sequencing compiled by the Genomics Evidence Neoplasia Information Exchange (GENIE) or Foundation Medicine (FM). A statistical approach was developed to identify patients with GCV-associated mutational signature (GCVsig) from targeted sequenced data of tumour samples. Cell line exposure models were further used to quantify mutation burden and DNA damage caused by GCV and other antiviral and immunosuppressive drugs. RESULTS: Mutational profiles from 22 of 121,771 patient samples in the GENIE and FM cohorts showed evidence of GCVsig. A diverse range of cancers was represented. All patients with detailed clinical history available had previously undergone solid organ transplantation and received GCV and mycophenolate treatment. RAS hotspot mutations associated with GCVsig were present in 9 of the 22 samples, with all samples harbouring multiple GCV-associated protein-altering mutations in cancer driver genes. In vitro testing in cell lines showed that elevated DNA damage response and GCVsig are uniquely associated with GCV but not acyclovir, a structurally similar antiviral. Combination treatment of GCV with the immunosuppressant, mycophenolate mofetil (MMF), increased the misincorporation of GCV in genomic DNA and mutations attributed to GCVsig in cell lines and organoids. CONCLUSIONS: In summary, GCV can cause a diverse range of cancers. Its mutagenicity may be potentiated by other therapies, such as mycophenolate, commonly co-prescribed with GCV for post-transplant patients. Further investigation of the optimal use of these drugs could help reduce GCV-associated mutagenesis in post-transplant patients.

Comprehensive characterization of single-cell full-length isoforms in human and mouse with long-read sequencing.

A modified Chromium 10x droplet-based protocol that subsamples cells for both short-read and long-read (nanopore) sequencing together with a new computational pipeline (FLAMES) is developed to enable isoform discovery, splicing analysis, and mutation detection in single cells. We identify thousands of unannotated isoforms and find conserved functional modules that are enriched for alternative transcript usage in different cell types and species, including ribosome biogenesis and mRNA splicing. Analysis at the transcript level allows data integration with scATAC-seq on individual promoters, improved correlation with protein expression data, and linked mutations known to confer drug resistance to transcriptome heterogeneity.

Acquisition of the Recurrent Gly101Val Mutation in BCL2 Confers Resistance to Venetoclax in Patients with Progressive Chronic Lymphocytic Leukemia.

The BCL2 inhibitor venetoclax induces high rates of durable remission in patients with previously treated chronic lymphocytic leukemia (CLL). However, despite continuous daily treatment, leukemia recurs in most patients. To investigate the mechanisms of secondary resistance, we analyzed paired pre-venetoclax and progression samples from 15 patients with CLL progression enrolled on venetoclax clinical trials. The novel Gly101Val mutation in BCL2 was identified at progression in 7 patients, but not at study entry. It was first detectable after 19 to 42 months of therapy, and its emergence anticipated clinical disease progression by many months. Gly101Val reduces the affinity of BCL2 for venetoclax by ∼180-fold in surface plasmon resonance assays, thereby preventing the drug from displacing proapoptotic mediators from BCL2 in cells and conferring acquired resistance in cell lines and primary patient cells. This mutation provides new insights into the pathobiology of venetoclax resistance and provides a potential biomarker of impending clinical relapse. SIGNIFICANCE: Why CLL recurs in patients who achieve remission with the BCL2 inhibitor venetoclax has been unknown. We provide the first description of an acquired point mutation in BCL2 arising recurrently and exclusively in venetoclax-treated patients. The mutation reduces venetoclax binding and is sufficient to confer resistance.See related commentary by Thangavadivel and Byrd, p. 320.This article is highlighted in the In This Issue feature, p. 305.

A comparison of control samples for ChIP-seq of histone modifications.

The advent of high-throughput sequencing has allowed genome wide profiling of histone modifications by Chromatin ImmunoPrecipitation (ChIP) followed by sequencing (ChIP-seq). In this assay the histone mark of interest is enriched through a chromatin pull-down assay using an antibody for the mark. Due to imperfect antibodies and other factors, many of the sequenced fragments do not originate from the histone mark of interest, and are referred to as background reads. Background reads are not uniformly distributed and therefore control samples are usually used to estimate the background distribution at any given genomic position. The Encyclopedia of DNA Elements (ENCODE) Consortium guidelines suggest sequencing a whole cell extract (WCE, or "input") sample, or a mock ChIP reaction such as an IgG control, as a background sample. However, for a histone modification ChIP-seq investigation it is also possible to use a Histone H3 (H3) pull-down to map the underlying distribution of histones. In this paper we generated data from a hematopoietic stem and progenitor cell population isolated from mouse fetal liver to compare WCE and H3 ChIP-seq as control samples. The quality of the control samples is estimated by a comparison to pull-downs of histone modifications and to expression data. We find minor differences between WCE and H3 ChIP-seq, such as coverage in mitochondria and behavior close to transcription start sites. Where the two controls differ, the H3 pull-down is generally more similar to the ChIP-seq of histone modifications. However, the differences between H3 and WCE have a negligible impact on the quality of a standard analysis.