Multi-omics profiling of longitudinal samples reveals early genomic changes in follicular lymphoma.

Follicular lymphoma (FL) is the most common indolent type of B-cell non-Hodgkin lymphoma. Advances in treatment have improved overall survival, but early relapse or transformation to aggressive disease is associated with inferior outcome. To identify early genetic events and track tumor clonal evolution, we performed multi-omics analysis of 94 longitudinal biopsies from 44 FL patients; 22 with transformation (tFL) and 22 with relapse without transformation (nFL). Deep whole-exome sequencing confirmed recurrent mutations in genes encoding epigenetic regulators (CREBBP, KMT2D, EZH2, EP300), with similar mutational landscape in nFL and tFL patients. Calculation of genomic distances between longitudinal samples revealed complex evolutionary patterns in both subgroups. CREBBP and KMT2D mutations were identified as genetic events that occur early in the disease course, and cases with CREBBP KAT domain mutations had low risk of transformation. Gains in chromosomes 12 and 18 (TCF4), and loss in 6q were identified as early and stable copy number alterations. Identification of such early and stable genetic events may provide opportunities for early disease detection and disease monitoring. Integrative analysis revealed that tumors with EZH2 mutations exhibited reduced gene expression of numerous histone genes, including histone linker genes. This might contribute to the epigenetic dysregulation in FL.

Generation of a biliary tract cancer cell line atlas reveals molecular subtypes and therapeutic targets.

Biliary tract cancers (BTCs) are a group of deadly malignancies encompassing intrahepatic and extrahepatic cholangiocarcinoma, gallbladder carcinoma, and ampullary carcinoma. Here, we present the integrative analysis of 63 BTC cell lines via multi-omics clustering and genome- scale CRISPR screens, providing a platform to illuminate BTC biology and inform therapeutic development. We identify dependencies broadly enriched in BTC compared to other cancers as well as dependencies selective to the anatomic subtypes. Notably, cholangiocarcinoma cell lines are stratified into distinct lineage subtypes based on biliary or dual biliary/hepatocyte marker signatures, associated with dependency on specific lineage survival factors. Transcriptional analysis of patient specimens demonstrates the prognostic significance of these lineage subtypes. Additionally, we delineate strategies to enhance targeted therapies or to overcome resistance in cell lines with key driver gene mutations. Furthermore, clustering based on dependencies and proteomics data elucidates unexpected functional relationships, including a BTC subgroup with partial squamous differentiation. Thus, this cell line atlas reveals potential therapeutic targets in molecularly defined BTCs, unveils biologically distinct disease subtypes, and offers a vital resource for BTC research.

Tumor cell-based liquid biopsy using high-throughput microfluidic enrichment of entire leukapheresis product.

Circulating Tumor Cells (CTCs), interrogated by sampling blood from patients with cancer, contain multiple analytes, including intact RNA, high molecular weight DNA, proteins, and metabolic markers. However, the clinical utility of tumor cell-based liquid biopsy has been limited since CTCs are very rare, and current technologies cannot process the blood volumes required to isolate a sufficient number of tumor cells for in-depth assays. We previously described a high-throughput microfluidic prototype utilizing high-flow channels and amplification of cell sorting forces through magnetic lenses. Here, we apply this technology to analyze patient-derived leukapheresis products, interrogating a mean blood volume of 5.83 liters from patients with metastatic cancer, with a median of 2,799 CTCs purified per patient. Isolation of many CTCs from individual patients enables characterization of their morphological and molecular heterogeneity, including cell and nuclear size and RNA expression. It also allows robust detection of gene copy number variation, a definitive cancer marker with potential diagnostic applications. High-volume microfluidic enrichment of CTCs constitutes a new dimension in liquid biopsies.

Mesoscale DNA features impact APOBEC3A and APOBEC3B deaminase activity and shape tumor mutational landscapes.

Antiviral DNA cytosine deaminases APOBEC3A and APOBEC3B are major sources of mutations in cancer by catalyzing cytosine-to-uracil deamination. APOBEC3A preferentially targets single-stranded DNAs, with a noted affinity for DNA regions that adopt stem-loop secondary structures. However, the detailed substrate preferences of APOBEC3A and APOBEC3B have not been fully established, and the specific influence of the DNA sequence on APOBEC3A and APOBEC3B deaminase activity remains to be investigated. Here, we find that APOBEC3B also selectively targets DNA stem-loop structures, and they are distinct from those subjected to deamination by APOBEC3A. We develop Oligo-seq, an in vitro sequencing-based method to identify specific sequence contexts promoting APOBEC3A and APOBEC3B activity. Through this approach, we demonstrate that APOBEC3A and APOBEC3B deaminase activity is strongly regulated by specific sequences surrounding the targeted cytosine. Moreover, we identify the structural features of APOBEC3B and APOBEC3A responsible for their substrate preferences. Importantly, we determine that APOBEC3B-induced mutations in hairpin-forming sequences within tumor genomes differ from the DNA stem-loop sequences mutated by APOBEC3A. Together, our study provides evidence that APOBEC3A and APOBEC3B can generate distinct mutation landscapes in cancer genomes, driven by their unique substrate selectivity.

Distinguishing preferences of human APOBEC3A and APOBEC3B for cytosines in hairpin loops, and reflection of these preferences in APOBEC-signature cancer genome mutations.

The APOBEC3 enzymes convert cytosines in single-stranded DNA to uracils to protect against viruses and retrotransposons but can contribute to mutations that diversify tumors. To understand the mechanism of mutagenesis, we map the uracils resulting from expression of APOBEC3B or its catalytic carboxy-terminal domain (CTD) in Escherichia coli. Like APOBEC3A, the uracilomes of A3B and A3B-CTD show a preference to deaminate cytosines near transcription start sites and the lagging-strand replication templates and in hairpin loops. Both biochemical activities of the enzymes and genomic uracil distribution show that A3A prefers 3 nt loops the best, while A3B prefers 4 nt loops. Reanalysis of hairpin loop mutations in human tumors finds intrinsic characteristics of both the enzymes, with a much stronger contribution from A3A. We apply Hairpin Signatures 1 and 2, which define A3A and A3B preferences respectively and are orthogonal to published methods, to evaluate their contribution to human tumor mutations.

Acquired Cross-Resistance in Small Cell Lung Cancer due to Extrachromosomal DNA Amplification of MYC Paralogs.

Small cell lung cancer (SCLC) presents as a highly chemosensitive malignancy but acquires cross-resistance after relapse. This transformation is nearly inevitable in patients but has been difficult to capture in laboratory models. Here, we present a preclinical system that recapitulates acquired cross-resistance, developed from 51 patient-derived xenograft (PDX) models. Each model was tested in vivo against three clinical regimens: cisplatin plus etoposide, olaparib plus temozolomide, and topotecan. These drug-response profiles captured hallmark clinical features of SCLC, such as the emergence of treatment-refractory disease after early relapse. For one patient, serial PDX models revealed that cross-resistance was acquired through MYC amplification on extrachromosomal DNA (ecDNA). Genomic and transcriptional profiles of the full PDX panel revealed that MYC paralog amplifications on ecDNAs were recurrent in relapsed cross-resistant SCLC, and this was corroborated in tumor biopsies from relapsed patients. We conclude that ecDNAs with MYC paralogs are recurrent drivers of cross-resistance in SCLC. SIGNIFICANCE: SCLC is initially chemosensitive, but acquired cross-resistance renders this disease refractory to further treatment and ultimately fatal. The genomic drivers of this transformation are unknown. We use a population of PDX models to discover that amplifications of MYC paralogs on ecDNA are recurrent drivers of acquired cross-resistance in SCLC. This article is featured in Selected Articles from This Issue, p. 695.

Blood-based screening for HPV-associated cancers.

Background: HPV-associated oropharyngeal cancer (HPV+OPSCC) is the most common HPV-associated cancer in the United States yet unlike cervical cancer lacks a screening test. HPV+OPSCCs are presumed to start developing 10-15 years prior to clinical diagnosis. Circulating tumor HPV DNA (ctHPVDNA) is a sensitive and specific biomarker for HPV+OPSCC. Taken together, blood-based screening for HPV+OPSCC may be feasible years prior to diagnosis. Methods: We developed an HPV whole genome sequencing assay, HPV-DeepSeek, with 99% sensitivity and specificity at clinical diagnosis. 28 plasma samples from HPV+OPSCC patients collected 1.3-10.8 years prior to diagnosis along with 1:1 age and gender-matched controls were run on HPV-DeepSeek and an HPV serology assay. Results: 22/28 (79%) of cases and 0/28 controls screened positive for HPV+OPSCC with 100% detection within four years of diagnosis and a maximum lead time of 7.8 years. We next applied a machine learning model classifying 27/28 cases (96%) with 100% detection within 10 years. Plasma-based PIK3CA gene mutations, viral genome integration events and HPV serology were used to orthogonally validate cancer detection with 68% (19/28) of the cohort having multiple cancer signals detected. Molecular fingerprinting of HPV genomes was performed across patients demonstrating that each viral genome was unique, ruling out contamination. In patients with tumor blocks from diagnosis (15/28), molecular fingerprinting was performed within patients confirming the same viral genome across time. Conclusions: We demonstrate accurate blood-based detection of HPV-associated cancers with lead times up to 10 years before clinical cancer diagnosis and in doing so, highlight the enormous potential of ctDNA-based cancer screening.

APOBEC3A induces DNA gaps through PRIMPOL and confers gap-associated therapeutic vulnerability.

Mutation signatures associated with apolipoprotein B mRNA editing catalytic polypeptide-like 3A/B (APOBEC3A/B) cytidine deaminases are prevalent across cancers, implying their roles as mutagenic drivers during tumorigenesis and tumor evolution. APOBEC3A (A3A) expression induces DNA replication stress and increases the cellular dependency on the ataxia telangiectasia and Rad3-related (ATR) kinase for survival. Nonetheless, how A3A induces DNA replication stress remains unclear. We show that A3A induces replication stress without slowing replication forks. We find that A3A induces single-stranded DNA (ssDNA) gaps through PrimPol-mediated repriming. A3A-induced ssDNA gaps are repaired by multiple pathways involving ATR, RAD51, and translesion synthesis. Both ATR inhibition and trapping of poly(ADP-ribose) polymerase (PARP) on DNA by PARP inhibitor impair the repair of A3A-induced gaps, preferentially killing A3A-expressing cells. When used in combination, PARP and ATR inhibitors selectively kill A3A-expressing cells synergistically in a manner dependent on PrimPol-generated gaps. Thus, A3A-induced replication stress arises from PrimPol-generated ssDNA gaps, which confer a therapeutic vulnerability to gap-targeted DNA repair inhibitors.

Patterns in the tapestry of chromatin-bound RB.

The retinoblastoma protein (RB)-mediated regulation of E2F is a component of a highly conserved cell cycle machine. However, RB's tumor suppressor activity, like RB's requirement in animal development, is tissue-specific, context-specific, and sometimes appears uncoupled from cell proliferation. Detailed new information about RB's genomic distribution provides a new perspective on the complexity of RB function, suggesting that some of its functional specificity results from context-specific RB association with chromatin. Here we summarize recent evidence showing that RB targets different types of chromatin regulatory elements at different cell cycle stages. RB controls traditional RB/E2F targets prior to S-phase, but, when cells proliferate, RB redistributes to cell type-specific chromatin loci. We discuss the broad implications of the new data for RB research.

ATR inhibition induces synthetic lethality in mismatch repair-deficient cells and augments immunotherapy.

The mismatch repair (MMR) deficiency of cancer cells drives mutagenesis and offers a useful biomarker for immunotherapy. However, many MMR-deficient (MMR-d) tumors do not respond to immunotherapy, highlighting the need for alternative approaches to target MMR-d cancer cells. Here, we show that inhibition of the ATR kinase preferentially kills MMR-d cancer cells. Mechanistically, ATR inhibitor (ATRi) imposes synthetic lethality on MMR-d cells by inducing DNA damage in a replication- and MUS81 nuclease-dependent manner. The DNA damage induced by ATRi is colocalized with both MSH2 and PCNA, suggesting that it arises from DNA structures recognized by MMR proteins during replication. In syngeneic mouse models, ATRi effectively reduces the growth of MMR-d tumors. Interestingly, the antitumor effects of ATRi are partially due to CD8+ T cells. In MMR-d cells, ATRi stimulates the accumulation of nascent DNA fragments in the cytoplasm, activating the cGAS-mediated interferon response. The combination of ATRi and anti-PD-1 antibody reduces the growth of MMR-d tumors more efficiently than ATRi or anti-PD-1 alone, showing the ability of ATRi to augment the immunotherapy of MMR-d tumors. Thus, ATRi selectively targets MMR-d tumor cells by inducing synthetic lethality and enhancing antitumor immunity, providing a promising strategy to complement and augment MMR deficiency-guided immunotherapy.

The RNA m5C modification in R-loops as an off switch of Alt-NHEJ.

The roles of R-loops and RNA modifications in homologous recombination (HR) and other DNA double-stranded break (DSB) repair pathways remain poorly understood. Here, we find that DNA damage-induced RNA methyl-5-cytosine (m5C) modification in R-loops plays a crucial role to regulate PARP1-mediated poly ADP-ribosylation (PARylation) and the choice of DSB repair pathways at sites of R-loops. Through bisulfite sequencing, we discover that the methyltransferase TRDMT1 preferentially generates m5C after DNA damage in R-loops across the genome. In the absence of m5C, R-loops activate PARP1-mediated PARylation both in vitro and in cells. Concurrently, m5C promotes transcription-coupled HR (TC-HR) while suppressing PARP1-dependent alternative non-homologous end joining (Alt-NHEJ), favoring TC-HR over Alt-NHEJ in transcribed regions as the preferred repair pathway. Importantly, simultaneous disruption of both TC-HR and Alt-NHEJ with TRDMT1 and PARP or Polymerase θ inhibitors prevents alternative DSB repair and exhibits synergistic cytotoxic effects on cancer cells, suggesting an effective strategy to exploit genomic instability in cancer therapy.

Mesoscale DNA Features Impact APOBEC3A and APOBEC3B Deaminase Activity and Shape Tumor Mutational Landscapes.

Antiviral DNA cytosine deaminases APOBEC3A and APOBEC3B are major sources of mutations in cancer by catalyzing cytosine-to-uracil deamination. APOBEC3A preferentially targets singlestranded DNAs, with a noted affinity for DNA regions that adopt stem-loop secondary structures. However, the detailed substrate preferences of APOBEC3A and APOBEC3B have been fully established, and the specific influence of the DNA sequence on APOBEC3A APOBEC3B deaminase activity remains to be investigated. Here, we find that APOBEC3B selectively targets DNA stem-loop structures, and they are distinct from those subjected deamination by APOBEC3A. We develop Oligo-seq, a novel in vitro sequencing-based to identify specific sequence contexts promoting APOBEC3A and APOBEC3B activity. Through this approach, we demonstrate that APOBEC3A an APOBEC3B deaminase activity is strongly regulated by specific sequences surrounding the targeted cytosine. Moreover, we identify structural features of APOBEC3B and APOBEC3A responsible for their substrate preferences. Importantly, we determine that APOBEC3B-induced mutations in hairpin-forming sequences within tumor genomes differ from the DNA stem-loop sequences mutated by APOBEC3A. Together, our study provides evidence that APOBEC3A and APOBEC3B can generate mutation landscapes in cancer genomes, driven by their unique substrate selectivity.

Distinguishing preferences of human APOBEC3A and APOBEC3B for cytosines in hairpin loops, and reflection of these preferences in APOBEC-signature cancer genome mutations.

The APOBEC3 family of enzymes convert cytosines in single-stranded DNA to uracils thereby causing mutations. These enzymes protect human cells against viruses and retrotransposons, but in many cancers they contribute to mutations that diversify the tumors and help them escape anticancer drug treatments. To understand the mechanism of mutagenesis by APOBEC3B, we expressed the complete enzyme or its catalytic carboxy-terminal domain (CTD) in repair-deficient Eschericia coli and mapped the resulting uracils using uracil pull-down and sequencing technology. The uracilomes of A3B-full and A3B-CTD showed peaks in many of the same regions where APOBEC3A also created uracilation peaks. Like A3A, the two A3B enzymes also preferred to deaminate cytosines near transcription start sites and in the lagging-strand template at replication forks. In contrast to an earlier report that A3B does not favor hairpin loops over linear DNA, we found that both A3B-full and A3B-CTD showed a strong preference for cytosines in hairpin loops. The major difference between A3A and A3B was that while the former enzyme prefers 3 nt loops the best, A3B prefers loops of 4 nt over those of other lengths. Furthermore, within 5 nt loops, A3A prefers cytosine to be in the penultimate position, while A3B prefers it to be at the 3' end of the loop. We confirmed these loop size and sequence preferences experimentally using purified A3A and A3B-CTD proteins. Reanalysis of hairpin loop mutations in human tumors using the size, sequence and position preferences of the two enzymes found that the tumors displayed mutations with intrinsic characteristics of both the enzymes with a stronger contribution from A3A.

Therapy-induced APOBEC3A drives evolution of persistent cancer cells.

Acquired drug resistance to anticancer targeted therapies remains an unsolved clinical problem. Although many drivers of acquired drug resistance have been identified1-4, the underlying molecular mechanisms shaping tumour evolution during treatment are incompletely understood. Genomic profiling of patient tumours has implicated apolipoprotein B messenger RNA editing catalytic polypeptide-like (APOBEC) cytidine deaminases in tumour evolution; however, their role during therapy and the development of acquired drug resistance is undefined. Here we report that lung cancer targeted therapies commonly used in the clinic can induce cytidine deaminase APOBEC3A (A3A), leading to sustained mutagenesis in drug-tolerant cancer cells persisting during therapy. Therapy-induced A3A promotes the formation of double-strand DNA breaks, increasing genomic instability in drug-tolerant persisters. Deletion of A3A reduces APOBEC mutations and structural variations in persister cells and delays the development of drug resistance. APOBEC mutational signatures are enriched in tumours from patients with lung cancer who progressed after extended responses to targeted therapies. This study shows that induction of A3A in response to targeted therapies drives evolution of drug-tolerant persister cells, suggesting that suppression of A3A expression or activity may represent a potential therapeutic strategy in the prevention or delay of acquired resistance to lung cancer targeted therapy.

Acquired Cross-resistance in Small Cell Lung Cancer due to Extrachromosomal DNA Amplification of MYC paralogs.

Small cell lung cancer (SCLC) presents as a highly chemosensitive malignancy but acquires cross-resistance after relapse. This transformation is nearly inevitable in patients but has been difficult to capture in laboratory models. Here we present a pre-clinical system that recapitulates acquired cross-resistance in SCLC, developed from 51 patient-derived xenografts (PDXs). Each model was tested for in vivo sensitivity to three clinical regimens: cisplatin plus etoposide, olaparib plus temozolomide, and topotecan. These functional profiles captured hallmark clinical features, such as the emergence of treatment-refractory disease after early relapse. Serially derived PDX models from the same patient revealed that cross-resistance was acquired through a MYC amplification on extrachromosomal DNA (ecDNA). Genomic and transcriptional profiles of the full PDX panel revealed that this was not unique to one patient, as MYC paralog amplifications on ecDNAs were recurrent among cross-resistant models derived from patients after relapse. We conclude that ecDNAs with MYC paralogs are recurrent drivers of cross-resistance in SCLC. SIGNIFICANCE: SCLC is initially chemosensitive, but acquired cross-resistance renders this disease refractory to further treatment and ultimately fatal. The genomic drivers of this transformation are unknown. We use a population of PDX models to discover that amplifications of MYC paralogs on ecDNA are recurrent drivers of acquired cross-resistance in SCLC.

DNA hypomethylation silences anti-tumor immune genes in early prostate cancer and CTCs.

Cancer is characterized by hypomethylation-associated silencing of large chromatin domains, whose contribution to tumorigenesis is uncertain. Through high-resolution genome-wide single-cell DNA methylation sequencing, we identify 40 core domains that are uniformly hypomethylated from the earliest detectable stages of prostate malignancy through metastatic circulating tumor cells (CTCs). Nested among these repressive domains are smaller loci with preserved methylation that escape silencing and are enriched for cell proliferation genes. Transcriptionally silenced genes within the core hypomethylated domains are enriched for immune-related genes; prominent among these is a single gene cluster harboring all five CD1 genes that present lipid antigens to NKT cells and four IFI16-related interferon-inducible genes implicated in innate immunity. The re-expression of CD1 or IFI16 murine orthologs in immuno-competent mice abrogates tumorigenesis, accompanied by the activation of anti-tumor immunity. Thus, early epigenetic changes may shape tumorigenesis, targeting co-located genes within defined chromosomal loci. Hypomethylation domains are detectable in blood specimens enriched for CTCs.

TGF-ß in the microenvironment induces a physiologically occurring immune-suppressive senescent state.

TGF-ß induces senescence in embryonic tissues. Whether TGF-ß in the hypoxic tumor microenvironment (TME) induces senescence in cancer and how the ensuing senescence-associated secretory phenotype (SASP) remodels the cellular TME to influence immune checkpoint inhibitor (ICI) responses are unknown. We show that TGF-ß induces a deeper senescent state under hypoxia than under normoxia; deep senescence correlates with the degree of E2F suppression and is marked by multinucleation, reduced reentry into proliferation, and a distinct 14-gene SASP. Suppressing TGF-ß signaling in tumors in an immunocompetent mouse lung cancer model abrogates endogenous senescent cells and suppresses the 14-gene SASP and immune infiltration. Untreated human lung cancers with a high 14-gene SASP display immunosuppressive immune infiltration. In a lung cancer clinical trial of ICIs, elevated 14-gene SASP is associated with increased senescence, TGF-ß and hypoxia signaling, and poor progression-free survival. Thus, TME-induced senescence may represent a naturally occurring state in cancer, contributing to an immune-suppressive phenotype associated with immune therapy resistance.

Chromatin-bound protein colocalization analysis using bedGraph2Cluster and PanChIP.

Computational pipelines for chromatin immunoprecipitation sequencing analysis can neglect colocalization events that occur in a mere subset of the genome. Here, we detail a streamlined approach for assessing colocalization of chromatin-bound proteins using the bedGraph2Cluster and PanChIP algorithms. Using histone modifications as an example, bedGraph2Cluster performs clustering analysis on chromatin binding patterns of target proteins. PanChIP then compares these clusters with a reference library of chromatin binding patterns and measures the overlap in peaks, capturing the heterogeneity in chromatin binding and colocalization patterns. For complete details on the use and execution of this protocol, please refer to Sanidas et al. (2022).1.