Publications

Acquired drug resistance to anticancer targeted therapies remains an unsolved clinical problem. Although many drivers of acquired drug resistance have been identified 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.

APOBEC mutagenesis, a major driver of cancer evolution, is known for targeting TpC sites in DNA. Recently, we showed that APOBEC3A (A3A) targets DNA hairpin loops. Here, we show that DNA secondary structure is in fact an orthogonal influence on A3A substrate optimality and, surprisingly, can override the TpC sequence preference. VpC (non-TpC) sites in optimal hairpins can outperform TpC sites as mutational hotspots. This expanded understanding of APOBEC mutagenesis illuminates the genomic Twin Paradox, a puzzling pattern of closely spaced mutation hotspots in cancer genomes, in which one is a canonical TpC site but the other is a VpC site, and double mutants are seen only in trans, suggesting a two-hit driver event. Our results clarify this paradox, revealing that both hotspots in these twins are optimal A3A substrates. Our findings reshape the notion of a mutation signature, highlighting the additive roles played by DNA sequence and DNA structure.

Intrinsic and acquired resistance to CDK4/6 inhibitors (CDK4/6i) is a critical challenge in hormone receptor-positive (HR+), human epidermal growth factor receptor 2-negative (HER2-) metastatic breast cancer (MBC). CDK2 inhibitors (CDK2i) show promise in this context, with several currently under early-phase clinical evaluation. Here, by analyzing tumor and liquid biopsies from patients treated with the selective CDK2i PF-07104091, we observed that PF-07104091 monotherapy achieved disease stabilization in a cohort of CDK4/6i-resistant HR+/HER2- MBC patients. Notably, responses were observed irrespective of RB expression and phosphorylation, whereas non-progressive disease appeared more frequently among patients whose tumors retained wild-type TP53. Studies in established and patient-derived cell lines further substantiated that the growth suppressive effects of CDK2i were independent of RB status in CDK4/6i-resistant models. CDK2i activity reduced DNA replication rate, increased DNA damage, and inhibited mitotic entry, with growth inhibition dependent on p53 expression. These findings indicate that CDK dependency shifts from CDK4/6 toward CDK2 as cells transition from CDK4/6i-sensitive to CDK4/6i-resistant state. The distinct mechanisms of CDK2i and CDK4/6i support enhanced synergistic activity in HR+/HER2- MBC with acquired resistance to CDK4/6i therapy.

Activation-induced deaminase (AID) and apolipoprotein B mRNA-editing enzyme catalytic subunit (APOBEC) enzymes convert cytosines to uracils, creating signature mutations that have been used to predict sites targeted by these enzymes. Mutation-based targeting maps are distorted by the error-prone or error-free repair of these uracils and by selection pressures. To directly map uracils created by AID/APOBEC enzymes, here we used uracil-DNA glycosylase and an alkoxyamine to covalently tag and sequence uracil-containing DNA fragments (UPD-Seq). We applied this technique to the genome of repair-defective, APOBEC3A-expressing bacterial cells and created a uracilation genome map, i.e. uracilome. The peak uracilated regions were in the 5′-ends of genes and operons mainly containing tRNA genes and a few protein-coding genes. We validated these findings through deep sequencing of pulldown regions and whole-genome sequencing of independent clones. The peaks were not correlated with high transcription rates or stable RNA:DNA hybrid formation. We defined the uracilation index (UI) as the frequency of occurrence of TT in UPD-Seq reads at different original TC dinucleotides. Genome-wide UI calculation confirmed that APOBEC3A modifies cytosines in the lagging-strand template during replication and in short hairpin loops. APOBEC3A's preference for tRNA genes was observed previously in yeast, and an analysis of human tumor sequences revealed that in tumors with a high percentage of APOBEC3 signature mutations, the frequency of tRNA gene mutations was much higher than in the rest of the genome. These results identify multiple causes underlying selection of cytosines by APOBEC3A for deamination, and demonstrate the utility of UPD-Seq.

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.

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.

Activation-induced deaminase (AID) is a DNA-cytosine deaminase that mediates maturation of antibodies through somatic hypermutation and class-switch recombination. While it causes mutations in immunoglobulin heavy and light chain genes and strand breaks in the switch regions of the immunoglobulin heavy chain gene, it largely avoids causing such damage in the rest of the genome. To help understand targeting by human AID, we expressed it in repair-deficient Escherichia coli and mapped the created uracils in the genomic DNA using uracil pull-down and sequencing, UPD-seq. We found that both AID and the human APOBEC3A preferentially target tRNA genes and transcription start sites, but do not show preference for highly transcribed genes. Unlike A3A, AID did not show a strong replicative strand bias or a preference for hairpin loops. Overlapping uracilation peaks between these enzymes contained binding sites for a protein, FIS, that helps create topological domains in the E. coli genome. To confirm whether these findings were relevant to B cells, we examined mutations from lymphoma and leukemia genomes within AID-preferred sequences. These mutations also lacked replicative strand bias or a hairpin loop preference. We propose here a model for how AID avoids causing mutations in the single-stranded DNA found within replication forks.

My PhD Dissertation on how a Novel Sequencing technique, UPD-Seq, was used to understand the genome-wide targeting of the human Aid/APOBEC3 enzymes. This work was done in the lab of Dr. Ashok S Bhagwat at Wayne State University.

Introduction of uracils in specific regions within immunoglobulin genes by the activation-induced deaminase (AID) promotes mutations and double-strand breaks (DSBs). Although uracils are repaired through multiple DNA repair pathways, previous work has used mutations or DSBs as proxies for uracils and not mapped the uracils directly. We mapped uracils in the Ig heavy chain gene, Igh, in a murine cell line, CH12F3, undergoing class-switch recombination (CSR) using the uracil pull-down and sequencing technique. These cells undergo IgM-to-IgA switch upon expression of AID but do not undergo somatic hypermutation. We mapped uracils in cells defective in uracil repair and show that AID introduces high levels of uracils only in parts of switch-mu and switch-alpha regions and not in constant regions, the Igh variable region or the light chain genes. Furthermore, the peaks of uracilation match the previously determined distribution of switch junctions, which are representative of DSBs that cause isotype switching. This work confirms that AID creates uracils in both DNA strands and shows that there is a direct correlation between uracil creation and DSBs in the relevant switch regions. We evaluate proposed mechanisms of CSR in light of these findings and show that mapping uracils provides a fresh perspective on CSR.

Current standard initial therapy for advanced, ROS proto-oncogene 1, receptor tyrosine kinase fusion (ROS1)-positive (ROS1+) non–small cell lung cancer (NSCLC) is crizotinib or entrectinib. Lorlatinib, a next-generation anaplastic lymphoma kinase/ROS1 inhibitor, recently demonstrated efficacy in ROS1+ NSCLC, including in crizotinib-pretreated patients. However, mechanisms of lorlatinib resistance in ROS1+ disease remain poorly understood. Here, we assessed mechanisms of resistance to crizotinib and lorlatinib.

Most B cell cancers overexpress the enzyme activation-induced deaminase at high levels and this enzyme converts cytosines in DNA to uracil. The constitutive expression of this enzyme in these cells greatly increases the uracil content of their genomes. We show here that these genomes also contain high levels of abasic sites presumably created during the repair of uracils through base-excision repair. We further show that three alkoxyamines with an alkyne functional group covalently link to abasic sites in DNA and kill immortalized cell lines created from B cell lymphomas, but not other cancers. They also do not kill normal B cells. Treatment of cancer cells with one of these chemicals causes strand breaks, and the sensitivity of the cells to this chemical depends on the ability of the cells to go through the S phase. However, other alkoxyamines that also link to abasic sites- but lack the alkyne functionality- do not kill cells from B cell lymphomas. This shows that the ability of alkoxyamines to covalently link to abasic sites is insufficient for their cytotoxicity and that the alkyne functionality may play a role in it. These chemicals violate the commonly accepted bioorthogonality of alkynes and are attractive prototypes for anti-B cell cancer agents.