DNA deaminases induce break-associated mutation showers with implication of APOBEC3B and 3A in breast cancer kataegis

Regulation, functional impact, and therapeutic targeting of APOBEC3A in cancer

Enzymes of the apolipoprotein B mRNA editing catalytic polypeptide like (APOBEC) family are cytosine deaminases that convert cytosine to uracil in DNA and RNA. Among these proteins, APOBEC3 sub-family members, APOBEC3A (A3A) and APOBEC3B (A3B), are prominent sources of mutagenesis in cancer cells. The aberrant expression of A3A and A3B in cancer cells leads to accumulation of mutations with specific single-base substitution (SBS) signatures, characterized by C→T and C→G changes, in a number of tumor types. In addition to fueling mutagenesis, A3A and A3B, particularly A3A, induce DNA replication stress, DNA damage, and chromosomal instability through their catalytic activities, triggering a range of cellular responses. Thus, A3A/B have emerged as key drivers of genome evolution during cancer development, contributing to tumorigenesis, tumor heterogeneity, and therapeutic resistance. Yet, the expression of A3A/B in cancer cells presents a cancer vulnerability that can be exploited therapeutically. In this review, we discuss the recent studies that shed light on the mechanisms regulating A3A expression and the impact of A3A in cancer. We also review recent advances in the development of A3A inhibitors and provide perspectives on the future directions of A3A research.

The cytidine deaminase APOBEC3C has unique sequence and genome feature preferences

APOBEC proteins are cytidine deaminases that restrict the replication of viruses and transposable elements. Several members of the APOBEC3 family, APOBEC3A, APOBEC3B, and APOBEC3H-I, can access the nucleus and cause what is thought to be indiscriminate deamination of the genome, resulting in mutagenesis and genome instability. Although APOBEC3C is also present in the nucleus, the full scope of its deamination target preferences is unknown. By expressing human APOBEC3C in a yeast model system, I have defined the APOBEC3C mutation signature, as well as the preferred genome features of APOBEC3C targets. The APOBEC3C mutation signature is distinct from those of the known cancer genome mutators APOBEC3A and APOBEC3B. APOBEC3C produces DNA strand-coordinated mutation clusters, and APOBEC3C mutations are enriched near the transcription start sites of active genes. Surprisingly, APOBEC3C lacks the bias for the lagging strand of DNA replication that is seen for APOBEC3A and APOBEC3B. The unique preferences of APOBEC3C constitute a mutation profile that will be useful in defining sites of APOBEC3C mutagenesis in human genomes.

The SMC5/6 complex prevents genotoxicity upon APOBEC3A-mediated replication stress

Mutational patterns caused by APOBEC3 cytidine deaminase activity are evident throughout human cancer genomes. In particular, the APOBEC3A family member is a potent genotoxin that causes substantial DNA damage in experimental systems and human tumors. However, the mechanisms that ensure genome stability in cells with active APOBEC3A are unknown. Through an unbiased genome-wide screen, we define the Structural Maintenance of Chromosomes 5/6 (SMC5/6) complex as essential for cell viability when APOBEC3A is active. We observe an absence of APOBEC3A mutagenesis in human tumors with SMC5/6 dysfunction, consistent with synthetic lethality. Cancer cells depleted of SMC5/6 incur substantial genome damage from APOBEC3A activity during DNA replication. Further, APOBEC3A activity results in replication tract lengthening which is dependent on PrimPol, consistent with re-initiation of DNA synthesis downstream of APOBEC3A-induced lesions. Loss of SMC5/6 abrogates elongated replication tracts and increases DNA breaks upon APOBEC3A activity. Our findings indicate that replication fork lengthening reflects a DNA damage response to APOBEC3A activity that promotes genome stability in an SMC5/6-dependent manner. Therefore, SMC5/6 presents a potential therapeutic vulnerability in tumors with active APOBEC3A.

Hypomorphic mutation in the large subunit of replication protein A affects mutagenesis by human APOBEC cytidine deaminases in yeast

Human APOBEC single-strand (ss) specific DNA and RNA cytidine deaminases change cytosines to uracils and function in antiviral innate immunity, RNA editing, and can cause hypermutation in chromosomes. The resulting uracils can be directly replicated, resulting in C to T mutations, or uracil-DNA glycosylase can convert the uracils to abasic (AP) sites which are then fixed as C to T or C to G mutations by translesion DNA polymerases. We noticed that in yeast and in human cancers, contributions of C to T and C to G mutations depends on the origin of ssDNA mutagenized by APOBECs. Since ssDNA in eukaryotic genomes readily binds to replication protein A (RPA) we asked if RPA could affect APOBEC-induced mutation spectrum in yeast. For that purpose, we expressed human APOBECs in the wild-type yeast and in strains carrying a hypomorph mutation rfa1-t33 in the large RPA subunit. We confirmed that the rfa1-t33 allele can facilitate mutagenesis by APOBECs. We also found that the rfa1-t33 mutation changed the ratio of APOBEC3A-induced T to C and T to G mutations in replicating yeast to resemble a ratio observed in long-persistent ssDNA in yeast and in cancers. We present the data suggesting that RPA may shield APOBEC formed uracils in ssDNA from Ung1, thereby facilitating C to T mutagenesis through the accurate copying of uracils by replicative DNA polymerases. Unexpectedly, we also found that for uracils shielded from Ung1 by wild-type RPA the mutagenic outcome is reduced in the presence of translesion DNA polymerase zeta.

Exploring the Role of Clustered Mutations in Carcinogenesis and Their Potential Clinical Implications in Cancer

Abnormal cell proliferation and growth leading to cancer primarily result from cumulative genome mutations. Single gene mutations alone do not fully explain cancer onset and progression; instead, clustered mutations-simultaneous occurrences of multiple mutations-are considered to be pivotal in cancer development and advancement. These mutations can affect different genes and pathways, resulting in cells undergoing malignant transformation with multiple functional abnormalities. Clustered mutations influence cancer growth rates, metastatic potential, and drug treatment sensitivity. This summary highlights the various types and characteristics of clustered mutations to understand their associations with carcinogenesis and discusses their potential clinical significance in cancer. As a unique mutation type, clustered mutations may involve genomic instability, DNA repair mechanism defects, and environmental exposures, potentially correlating with responsiveness to immunotherapy. Understanding the characteristics and underlying processes of clustered mutations enhances our comprehension of carcinogenesis and cancer progression, providing new diagnostic and therapeutic approaches for cancer.

APOBEC3 mutagenesis drives therapy resistance in breast cancer

Acquired genetic alterations commonly drive resistance to endocrine and targeted therapies in metastatic breast cancer 1-7 , however the underlying processes engendering these diverse alterations are largely uncharacterized. To identify the mutational processes operant in breast cancer and their impact on clinical outcomes, we utilized a well-annotated cohort of 3,880 patient samples with paired tumor-normal sequencing data. The mutational signatures associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) enzymes were highly prevalent and enriched in post-treatment compared to treatment-naïve hormone receptor-positive (HR+) cancers. APOBEC3 mutational signatures were independently associated with shorter progression-free survival on antiestrogen plus CDK4/6 inhibitor combination therapy in patients with HR+ metastatic breast cancer. Whole genome sequencing (WGS) of breast cancer models and selected paired primary-metastatic samples demonstrated that active APOBEC3 mutagenesis promoted resistance to both endocrine and targeted therapies through characteristic alterations such as RB1 loss-of-function mutations. Evidence of APOBEC3 activity in pre-treatment samples illustrated a pervasive role for this mutational process in breast cancer evolution. The study reveals APOBEC3 mutagenesis to be a frequent mediator of therapy resistance in breast cancer and highlights its potential as a biomarker and target for overcoming resistance.

Kataegis in clinical and molecular subgroups of primary breast cancer

Kataegis is a hypermutation phenomenon characterized by localized clusters of single base pair substitution (SBS) reported in multiple cancer types. Despite a high frequency in breast cancer, large-scale analyses of kataegis patterns and associations with clinicopathological and molecular variables in established breast cancer subgroups are lacking. Therefore, WGS profiled primary breast cancers (n = 791) with associated clinical and molecular data layers, like RNA-sequencing data, were analyzed for kataegis frequency, recurrence, and associations with genomic contexts and functional elements, transcriptional patterns, driver alterations, homologous recombination deficiency (HRD), and prognosis in tumor subgroups defined by ER, PR, and HER2/ERBB2 status. Kataegis frequency was highest in the HER2-positive(p) subgroups, including both ER-negative(n)/positive(p) tumors (ERnHER2p/ERpHER2p). In TNBC, kataegis was neither associated with PAM50 nor TNBC mRNA subtypes nor with distant relapse in chemotherapy-treated patients. In ERpHER2n tumors, kataegis was associated with aggressive characteristics, including PR-negativity, molecular Luminal B subtype, higher mutational burden, higher grade, and expression of proliferation-associated genes. Recurrent kataegis loci frequently targeted regions commonly amplified in ER-positive tumors, while few recurrent loci were observed in TNBC. SBSs in kataegis loci appeared enriched in regions of open chromatin. Kataegis status was not associated with HRD in any subgroup or with distinct transcriptional patterns in unsupervised or supervised analysis. In summary, kataegis is a common hypermutation phenomenon in established breast cancer subgroups, particularly in HER2p subgroups, coinciding with an aggressive tumor phenotype in ERpHER2n disease. In TNBC, the molecular implications and associations of kataegis are less clear, including its prognostic value.

The SMC5/6 complex prevents genotoxicity upon APOBEC3A-mediated replication stress

Mutational patterns caused by APOBEC3 cytidine deaminase activity are evident throughout human cancer genomes. In particular, the APOBEC3A family member is a potent genotoxin that causes substantial DNA damage in experimental systems and human tumors. However, the mechanisms that ensure genome stability in cells with active APOBEC3A are unknown. Through an unbiased genome-wide screen, we define the Structural Maintenance of Chromosomes 5/6 (SMC5/6) complex as essential for cell viability when APOBEC3A is active. We observe an absence of APOBEC3A mutagenesis in human tumors with SMC5/6 dysfunction, consistent with synthetic lethality. Cancer cells depleted of SMC5/6 incur substantial genome damage from APOBEC3A activity during DNA replication. Further, APOBEC3A activity results in replication tract lengthening which is dependent on PrimPol, consistent with re-initiation of DNA synthesis downstream of APOBEC3A-induced lesions. Loss of SMC5/6 abrogates elongated replication tracts and increases DNA breaks upon APOBEC3A activity. Our findings indicate that replication fork lengthening reflects a DNA damage response to APOBEC3A activity that promotes genome stability in an SMC5/6-dependent manner. Therefore, SMC5/6 presents a potential therapeutic vulnerability in tumors with active APOBEC3A.

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.

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.

Differentiation signals induce APOBEC3A expression via GRHL3 in squamous epithelia and squamous cell carcinoma

Two APOBEC (apolipoprotein-B mRNA editing enzyme catalytic polypeptide-like) DNA cytosine deaminase enzymes (APOBEC3A and APOBEC3B) generate somatic mutations in cancer, driving tumour development and drug resistance. Here we used single cell RNA sequencing to study APOBEC3A and APOBEC3B expression in healthy and malignant mucosal epithelia, validating key observations with immunohistochemistry, spatial transcriptomics and functional experiments. Whereas APOBEC3B is expressed in keratinocytes entering mitosis, we show that APOBEC3A expression is confined largely to terminally differentiating cells and requires Grainyhead-like transcription factor 3 (GRHL3). Thus, in normal tissue, neither deaminase appears to be expressed at high levels during DNA replication, the cell cycle stage associated with APOBEC-mediated mutagenesis. In contrast, we show that in squamous cell carcinoma tissues, there is expansion of GRHL3 expression and activity to a subset of cells undergoing DNA replication and concomitant extension of APOBEC3A expression to proliferating cells. These findings indicate a mechanism for acquisition of APOBEC3A mutagenic activity in tumours.

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.

APOBEC Mutagenesis in Cancer Development and Susceptibility

APOBEC cytosine deaminases are prominent mutators in cancer, mediating mutations in over 50% of cancers. APOBEC mutagenesis has been linked to tumor heterogeneity, persistent cell evolution, and therapy responses. While emerging evidence supports the impact of APOBEC mutagenesis on cancer progression, the understanding of its contribution to cancer susceptibility and malignant transformation is limited. We examine the existing evidence for the role of APOBEC mutagenesis in carcinogenesis on the basis of the reported associations between germline polymorphisms in genes encoding APOBEC enzymes and cancer risk, insights into APOBEC activities from sequencing efforts of both malignant and non-malignant human tissues, and in vivo studies. We discuss key knowledge gaps and highlight possible ways to gain a deeper understanding of the contribution of APOBEC mutagenesis to cancer development.

Genomic alterations driving precancerous to cancerous lesions in esophageal cancer development

Esophageal squamous cell carcinoma (ESCC) develops through a series of increasingly abnormal precancerous lesions. Previous studies have revealed the striking differences between normal esophageal epithelium and ESCC in copy number alterations (CNAs) and mutations in genes driving clonal expansion. However, due to limited data on early precancerous lesions, the timing of these transitions and which among them are prerequisites for malignant transformation remained unclear. Here, we analyze 1,275 micro-biopsies from normal esophagus, early and late precancerous lesions, and esophageal cancers to decipher the genomic alterations at each stage. We show that the frequency of TP53 biallelic inactivation increases dramatically in early precancerous lesion stage while CNAs and APOBEC mutagenesis substantially increase at late stages. TP53 biallelic loss is the prerequisite for the development of CNAs of genes in cell cycle, DNA repair, and apoptosis pathways, suggesting it might be one of the earliest steps initiating malignant transformation.

An impaired ubiquitin-proteasome system increases APOBEC3A abundance

Apolipoprotein B messenger RNA (mRNA) editing enzyme, catalytic polypeptide-like (APOBEC) cytidine deaminases cause genetic instability during cancer development. Elevated APOBEC3A (A3A) levels result in APOBEC signature mutations; however, mechanisms regulating A3A abundance in breast cancer are unknown. Here, we show that dysregulating the ubiquitin-proteasome system with proteasome inhibitors, including Food and Drug Administration-approved anticancer drugs, increased A3A abundance in breast cancer and multiple myeloma cell lines. Unexpectedly, elevated A3A occurs via an ∼100-fold increase in A3A mRNA levels, indicating that proteasome inhibition triggers a transcriptional response as opposed to or in addition to blocking A3A degradation. This transcriptional regulation is mediated in part through FBXO22, a protein that functions in SKP1-cullin-F-box ubiquitin ligase complexes and becomes dysregulated during carcinogenesis. Proteasome inhibitors increased cellular cytidine deaminase activity, decreased cellular proliferation and increased genomic DNA damage in an A3A-dependent manner. Our findings suggest that proteasome dysfunction, either acquired during cancer development or induced therapeutically, could increase A3A-induced genetic heterogeneity and thereby influence therapeutic responses in patients.

Hepatitis B virus infection disrupts homologous recombination in hepatocellular carcinoma by stabilizing resection inhibitor ADRM1

Many cancers harbor homologous recombination defects (HRDs). A HRD is a therapeutic target that is being successfully utilized in treatment of breast/ovarian cancer via synthetic lethality. However, canonical HRD caused by BRCAness mutations do not prevail in liver cancer. Here we report a subtype of HRD caused by the perturbation of a proteasome variant (CDW19S) in hepatitis B virus-bearing (HBV-bearing) cells. This amalgamate protein complex contained the 19S proteasome decorated with CRL4WDR70 ubiquitin ligase, and assembled at broken chromatin in a PSMD4Rpn10- and ATM-MDC1-RNF8-dependent manner. CDW19S promoted DNA end processing via segregated modules that promote nuclease activities of MRE11 and EXO1. Contrarily, a proteasomal component, ADRM1Rpn13, inhibited resection and was removed by CRL4WDR70-catalyzed ubiquitination upon commitment of extensive resection. HBx interfered with ADRM1Rpn13 degradation, leading to the imposition of ADRM1Rpn13-dependent resection barrier and consequent viral HRD subtype distinguishable from that caused by BRCA1 defect. Finally, we demonstrated that viral HRD in HBV-associated hepatocellular carcinoma can be exploited to restrict tumor progression. Our work clarifies the underlying mechanism of a virus-induced HRD subtype.

APOBEC3G Is a p53-Dependent Restriction Factor in Respiratory Syncytial Virus Infection of Human Cells Included in the p53/Immune Axis

Identifying and understanding genetic factors that influence the propagation of the human respiratory syncytial virus (RSV) can lead to health benefits and possibly augment recent vaccine approaches. We previously identified a p53/immune axis in which the tumor suppressor p53 directly regulates the expression of immune system genes, including the seven members of the APOBEC3 family of DNA cytidine deaminases (A3), which are innate immune sentinels against viral infections. Here, we examined the potential p53 and A3 influence in RSV infection, as well as the overall p53-dependent cellular and p53/immune axis responses to infection. Using a paired p53 model system of p53+ and p53- human lung tumor cells, we found that RSV infection activates p53, leading to the altered p53-dependent expression of A3D, A3F, and A3G, along with p53 site-specific binding. Focusing on A3G because of its 10-fold-greater p53 responsiveness to RSV, the overexpression of A3G can reduce RSV viral replication and syncytial formation. We also observed that RSV-infected cells undergo p53-dependent apoptosis. The study was expanded to globally address at the transcriptional level the p53/immune axis response to RSV. Nearly 100 genes can be directly targeted by the p53/immune axis during RSV infection based on our p53BAER analysis (Binding And Expression Resource). Overall, we identify A3G as a potential p53-responsive restriction factor in RSV infection. These findings have significant implications for RSV clinical and therapeutic studies and other p53-influenced viral infections, including using p53 adjuvants to boost the response of A3 genes.

A Rapid Antibody Enhancement Platform in Saccharomyces cerevisiae Using an Improved, Diversifying CRISPR Base Editor

The yeast Saccharomyces cerevisiae is commonly used to interrogate and screen protein variants and to perform directed evolution studies to develop proteins with enhanced features. While several techniques have been described that help enable the use of yeast for directed evolution, there remains a need to increase their speed and ease of use. Here we present yDBE, a yeast diversifying base editor that functions in vivo and employs a CRISPR-dCas9-directed cytidine deaminase base editor to diversify DNA in a targeted, rapid, and high-breadth manner. To develop yDBE, we enhanced the mutation rate of an initial base editor by employing improved deaminase variants and characterizing several scaffolded guide constructs. We then demonstrate the ability of the yDBE platform to improve the affinity of a displayed antibody scFv, rapidly generating diversified libraries and isolating improved binders via cell sorting. By performing high-throughput sequencing analysis of the high-activity yDBE, we show that it enables a mutation rate of 2.13 × 10-4 substitutions/bp/generation over a window of 100 bp. As yDBE functions entirely in vivo and can be easily programmed to diversify nearly any such window of DNA, we posit that it can be a powerful tool for facilitating a variety of directed evolution experiments.

Mutational impact of APOBEC3A and APOBEC3B in a human cell line and comparisons to breast cancer

A prominent source of mutation in cancer is single-stranded DNA cytosine deamination by cellular APOBEC3 enzymes, which results in signature C-to-T and C-to-G mutations in TCA and TCT motifs. Although multiple enzymes have been implicated, reports conflict and it is unclear which protein(s) are responsible. Here we report the development of a selectable system to quantify genome mutation and demonstrate its utility by comparing the mutagenic activities of three leading candidates-APOBEC3A, APOBEC3B, and APOBEC3H. The human cell line, HAP1, is engineered to express the thymidine kinase (TK) gene of HSV-1, which confers sensitivity to ganciclovir. Expression of APOBEC3A and APOBEC3B, but not catalytic mutant controls or APOBEC3H, triggers increased frequencies of TK mutation and similar TC-biased cytosine mutation profiles in the selectable TK reporter gene. Whole genome sequences from independent clones enabled an analysis of thousands of single base substitution mutations and extraction of local sequence preferences with APOBEC3A preferring YTCW motifs 70% of the time and APOBEC3B 50% of the time (Y = C/T; W = A/T). Signature comparisons with breast tumor whole genome sequences indicate that most malignancies manifest intermediate percentages of APOBEC3 signature mutations in YTCW motifs, mostly between 50 and 70%, suggesting that both enzymes contribute in a combinatorial manner to the overall mutation landscape. Although the vast majority of APOBEC3A- and APOBEC3B-induced single base substitution mutations occur outside of predicted chromosomal DNA hairpin structures, whole genome sequence analyses and supporting biochemical studies also indicate that both enzymes are capable of deaminating the single-stranded loop regions of DNA hairpins at elevated rates. These studies combine to help resolve a long-standing etiologic debate on the source of APOBEC3 signature mutations in cancer and indicate that future diagnostic and therapeutic efforts should focus on both APOBEC3A and APOBEC3B.

Non-homogeneous Poisson and renewal processes as spatial models for cancer mutation

Advances in sequencing technology assisted biologists in revealing signatures of DNA cancer mutation process and in demonstrating the mutagenesis behind. However, most of these signatures proposed by majority of work focus only on the type and frequency of mutations, without considering spatial information which is non-negligible in exploring mechanisms of mutation occurrence, e.g., Kataegis. Statistical characterization of the distance between consecutive mutations can give us relative spatial information; however, it ignores location information which is as important as distance information. In this work, we assume that DNA cancer mutations are location-dependent and that integrating the two variables, location and inter-distance, is beneficial to study DNA cancer mutation processes more accurately. Particularly, instead of following a specific distribution over the whole DNA sequence, we found out that the distribution of distance between successive mutations alternates between exponential and power-law distributions. Apart from this, the parameters of either of the two distributions vary with DNA locations. The cancers with kataegis phenomenon, a specific mutation pattern caused by abnormal activity of APOBEC protein family, are more likely to be accompanied by higher parameter values of distance distribution, implying higher occurrence rate of mutation. Therefore, we propose non-homogeneous Poisson and Renewal processes to spatially model DNA cancer mutations and to describe mutation patterns quantitatively and more accurately through a statistical perspective.

APOBEC3B regulates R-loops and promotes transcription-associated mutagenesis in cancer

The single-stranded DNA cytosine-to-uracil deaminase APOBEC3B is an antiviral protein implicated in cancer. However, its substrates in cells are not fully delineated. Here APOBEC3B proteomics reveal interactions with a surprising number of R-loop factors. Biochemical experiments show APOBEC3B binding to R-loops in cells and in vitro. Genetic experiments demonstrate R-loop increases in cells lacking APOBEC3B and decreases in cells overexpressing APOBEC3B. Genome-wide analyses show major changes in the overall landscape of physiological and stimulus-induced R-loops with thousands of differentially altered regions, as well as binding of APOBEC3B to many of these sites. APOBEC3 mutagenesis impacts genes overexpressed in tumors and splice factor mutant tumors preferentially, and APOBEC3-attributed kataegis are enriched in RTCW motifs consistent with APOBEC3B deamination. Taken together with the fact that APOBEC3B binds single-stranded DNA and RNA and preferentially deaminates DNA, these results support a mechanism in which APOBEC3B regulates R-loops and contributes to R-loop mutagenesis in cancer.

Novel Insights in the Genomics of Anaplastic Thyroid Carcinoma: A Role for Cyclin-Dependent Kinase Inhibition?

Anaplastic thyroid carcinoma (ATC) stands as a rare but extraordinarily lethal tumor, marked by its limited treatment options [...].

Loss of SYNCRIP unleashes APOBEC-driven mutagenesis, tumor heterogeneity, and AR-targeted therapy resistance in prostate cancer

Tumor mutational burden and heterogeneity has been suggested to fuel resistance to many targeted therapies. The cytosine deaminase APOBEC proteins have been implicated in the mutational signatures of more than 70% of human cancers. However, the mechanism underlying how cancer cells hijack the APOBEC mediated mutagenesis machinery to promote tumor heterogeneity, and thereby foster therapy resistance remains unclear. We identify SYNCRIP as an endogenous molecular brake which suppresses APOBEC-driven mutagenesis in prostate cancer (PCa). Overactivated APOBEC3B, in SYNCRIP-deficient PCa cells, is a key mutator, representing the molecular source of driver mutations in some frequently mutated genes in PCa, including FOXA1, EP300. Functional screening identifies eight crucial drivers for androgen receptor (AR)-targeted therapy resistance in PCa that are mutated by APOBEC3B: BRD7, CBX8, EP300, FOXA1, HDAC5, HSF4, STAT3, and AR. These results uncover a cell-intrinsic mechanism that unleashes APOBEC-driven mutagenesis, which plays a significant role in conferring AR-targeted therapy resistance in PCa.

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.

Similar deamination activities but different phenotypic outcomes induced by APOBEC3 enzymes in breast epithelial cells

APOBEC3 (A3) enzymes deaminate cytosine to uracil in viral single-stranded DNA as a mutagenic barrier for some viruses. A3-induced deaminations can also occur in human genomes resulting in an endogenous source of somatic mutations in multiple cancers. However, the roles of each A3 are unclear since few studies have assessed these enzymes in parallel. Thus, we developed stable cell lines expressing A3A, A3B, or A3H Hap I using non-tumorigenic MCF10A and tumorigenic MCF7 breast epithelial cells to assess their mutagenic potential and cancer phenotypes in breast cells. The activity of these enzymes was characterized by γH2AX foci formation and in vitro deamination. Cell migration and soft agar colony formation assays assessed cellular transformation potential. We found that all three A3 enzymes had similar γH2AX foci formation, despite different deamination activities in vitro. Notably, in nuclear lysates, the in vitro deaminase activity of A3A, A3B, and A3H did not require digestion of cellular RNA, in contrast to that of A3B and A3H in whole-cell lysates. Their similar activities in cells, nonetheless, resulted in distinct phenotypes where A3A decreased colony formation in soft agar, A3B decreased colony formation in soft agar after hydroxyurea treatment, and A3H Hap I promoted cell migration. Overall, we show that in vitro deamination data do not always reflect cell DNA damage, all three A3s induce DNA damage, and the impact of each is different.

APOBEC3G protects the genome of human cultured cells and mice from radiation-induced damage

Cytosine deaminases AID/APOBEC proteins act as potent nucleic acid editors, playing important roles in innate and adaptive immunity. However, the mutagenic effects of some of these proteins compromise genomic integrity and may promote tumorigenesis. Here, we demonstrate that human APOBEC3G (A3G), in addition to its role in innate immunity, promotes repair of double-strand breaks (DSBs) in vitro and in vivo. Transgenic mice expressing A3G successfully survived lethal irradiation, whereas wild-type controls quickly succumbed to radiation syndrome. Mass spectrometric analyses identified the differential upregulation of a plethora of proteins involved in DSB repair pathways in A3G-expressing cells early following irradiation to facilitate repair. Importantly, we find that A3G not only accelerates DSB repair but also promotes deamination-dependent error-free rejoining. These findings have two implications: (a) strategies aimed at inhibiting A3G may improve the efficacy of genotoxic therapies used to cure malignant tumours; and (b) enhancing A3G activity may reduce acute radiation syndrome in individuals exposed to ionizing radiation.

Influence network model uncovers relations between biological processes and mutational signatures

BACKGROUND

There has been a growing appreciation recently that mutagenic processes can be studied through the lenses of mutational signatures, which represent characteristic mutation patterns attributed to individual mutagens. However, the causal links between mutagens and observed mutation patterns as well as other types of interactions between mutagenic processes and molecular pathways are not fully understood, limiting the utility of mutational signatures.

METHODS

To gain insights into these relationships, we developed a network-based method, named GENESIGNET that constructs an influence network among genes and mutational signatures. The approach leverages sparse partial correlation among other statistical techniques to uncover dominant influence relations between the activities of network nodes.

RESULTS

Applying GENESIGNET to cancer data sets, we uncovered important relations between mutational signatures and several cellular processes that can shed light on cancer-related processes. Our results are consistent with previous findings, such as the impact of homologous recombination deficiency on clustered APOBEC mutations in breast cancer. The network identified by GENESIGNET also suggest an interaction between APOBEC hypermutation and activation of regulatory T Cells (Tregs), as well as a relation between APOBEC mutations and changes in DNA conformation. GENESIGNET also exposed a possible link between the SBS8 signature of unknown etiology and the Nucleotide Excision Repair (NER) pathway.

CONCLUSIONS

GENESIGNET provides a new and powerful method to reveal the relation between mutational signatures and gene expression. The GENESIGNET method was implemented in python, and installable package, source codes and the data sets used for and generated during this study are available at the Github site https://github.com/ncbi/GeneSigNet.

In situ Raman spectroscopy and machine learning unveil biomolecular alterations in invasive breast cancer

SIGNIFICANCE

As many as 60% of patients with early stage breast cancer undergo breast-conserving surgery. Of those, 20% to 35% need a second surgery because of incomplete resection of the lesions. A technology allowing in situ detection of cancer could reduce re-excision procedure rates and improve patient survival.

AIM

Raman spectroscopy was used to measure the spectral fingerprint of normal breast and cancer tissue ex-vivo. The aim was to build a machine learning model and to identify the biomolecular bands that allow one to detect invasive breast cancer.

APPROACH

The system was used to interrogate specimens from 20 patients undergoing lumpectomy, mastectomy, or breast reduction surgery. This resulted in 238 ex-vivo measurements spatially registered with standard histology classifying tissue as cancer, normal, or fat. A technique based on support vector machines led to the development of predictive models, and their performance was quantified using a receiver-operating-characteristic analysis.

RESULTS

Raman spectroscopy combined with machine learning detected normal breast from ductal or lobular invasive cancer with a sensitivity of 93% and a specificity of 95%. This was achieved using a model based on only two spectral bands, including the peaks associated with C-C stretching of proteins around 940    cm - 1 and the symmetric ring breathing at 1004    cm - 1 associated with phenylalanine.

CONCLUSIONS

Detection of cancer on the margins of surgically resected breast specimen is feasible with Raman spectroscopy.

The APOBEC3B cytidine deaminase is an adenovirus restriction factor

Human adenoviruses (HAdVs) are a large family of DNA viruses counting more than a hundred strains divided into seven species (A to G). HAdVs induce respiratory tract infections, gastroenteritis and conjunctivitis. APOBEC3B is a cytidine deaminase that restricts several DNA viruses. APOBEC3B is also implicated in numerous cancers where it is responsible for the introduction of clustered mutations into the cellular genome. In this study, we demonstrate that APOBEC3B is an adenovirus restriction factor acting through a deaminase-dependent mechanism. APOBEC3B introduces C-to-T clustered mutations into the adenovirus genome. APOBEC3B reduces the propagation of adenoviruses by limiting viral genome replication, progression to late phase, and production of infectious virions. APOBEC3B restriction efficiency varies between adenoviral strains, the A12 strain being more sensitive to APOBEC3B than the B3 or C2 strains. In A12-infected cells, APOBEC3B clusters in the viral replication centers. Importantly, we show that adenovirus infection leads to a reduction of the quantity and/or enzymatic activity of the APOBEC3B protein depending on the strains. The A12 strain seems less able to resist APOBEC3B than the B3 or C2 strains, a characteristic which could explain the strong depletion of the APOBEC3-targeted motifs in the A12 genome. These findings suggest that adenoviruses evolved different mechanisms to antagonize APOBEC3B. Elucidating these mechanisms could benefit the design of cancer treatments. This study also identifies adenoviruses as triggers of the APOBEC3B-mediated innate response. The involvement of certain adenoviral strains in the genesis of the APOBEC3 mutational signature observed in tumors deserves further study.

APOBEC mutagenesis is a common process in normal human small intestine

APOBEC mutational signatures SBS2 and SBS13 are common in many human cancer types. However, there is an incomplete understanding of its stimulus, when it occurs in the progression from normal to cancer cell and the APOBEC enzymes responsible. Here we whole-genome sequenced 342 microdissected normal epithelial crypts from the small intestines of 39 individuals and found that SBS2/SBS13 mutations were present in 17% of crypts, more frequent than most other normal tissues. Crypts with SBS2/SBS13 often had immediate crypt neighbors without SBS2/SBS13, suggesting that the underlying cause of SBS2/SBS13 is cell-intrinsic. APOBEC mutagenesis occurred in an episodic manner throughout the human lifespan, including in young children. APOBEC1 mRNA levels were very high in the small intestine epithelium, but low in the large intestine epithelium and other tissues. The results suggest that the high levels of SBS2/SBS13 in the small intestine are collateral damage from APOBEC1 fulfilling its physiological function of editing APOB mRNA.

Immunogenomic analysis of human brain metastases reveals diverse immune landscapes across genetically distinct tumors

Brain metastases (BrMs) are the most common form of brain tumors in adults and frequently originate from lung and breast primary cancers. BrMs are associated with high mortality, emphasizing the need for more effective therapies. Genetic profiling of primary tumors is increasingly used as part of the effort to guide targeted therapies against BrMs, and immune-based strategies for the treatment of metastatic cancer are gaining momentum. However, the tumor immune microenvironment (TIME) of BrM is extremely heterogeneous, and whether specific genetic profiles are associated with distinct immune states remains unknown. Here, we perform an extensive characterization of the immunogenomic landscape of human BrMs by combining whole-exome/whole-genome sequencing, RNA sequencing of immune cell populations, flow cytometry, immunofluorescence staining, and tissue imaging analyses. This revealed unique TIME phenotypes in genetically distinct lung- and breast-BrMs, thereby enabling the development of personalized immunotherapies tailored by the genetic makeup of the tumors.

Structure-Guided Design of a Potent and Specific Inhibitor against the Genomic Mutator APOBEC3A

Nucleic acid structure plays a critical role in governing the selectivity of DNA- and RNA-modifying enzymes. In the case of the APOBEC3 family of cytidine deaminases, these enzymes catalyze the conversion of cytosine (C) to uracil (U) in single-stranded DNA, primarily in the context of innate immunity. DNA deamination can also have pathological consequences, accelerating the evolution of viral genomes or, when the host genome is targeted by either APOBEC3A (A3A) or APOBEC3B (A3B), promoting tumor evolution leading to worse patient prognosis and chemotherapeutic resistance. For A3A, nucleic acid secondary structure has emerged as a critical determinant of substrate targeting, with a predilection for DNA that can form stem loop hairpins. Here, we report the development of a specific nanomolar-level, nucleic acid-based inhibitor of A3A. Our strategy relies on embedding the nucleobase 5-methylzebularine, a mechanism-based inhibitor, into a DNA dumbbell structure, which mimics the ideal substrate secondary structure for A3A. Structure-activity relationship studies using a panel of diverse inhibitors reveal a critical role for the stem and position of the inhibitor moiety in achieving potent inhibition. Moreover, we demonstrate that DNA dumbbell inhibitors, but not nonstructured inhibitors, show specificity against A3A relative to the closely related catalytic domain of A3B. Overall, our work demonstrates the feasibility of leveraging secondary structural preferences in inhibitor design, offering a blueprint for further development of modulators of DNA-modifying enzymes and potential therapeutics to circumvent APOBEC-driven viral and tumor evolution.

Regional mutational signature activities in cancer genomes

Cancer genomes harbor a catalog of somatic mutations. The type and genomic context of these mutations depend on their causes and allow their attribution to particular mutational signatures. Previous work has shown that mutational signature activities change over the course of tumor development, but investigations of genomic region variability in mutational signatures have been limited. Here, we expand upon this work by constructing regional profiles of mutational signature activities over 2,203 whole genomes across 25 tumor types, using data aggregated by the Pan-Cancer Analysis of Whole Genomes (PCAWG) consortium. We present GenomeTrackSig as an extension to the TrackSig R package to construct regional signature profiles using optimal segmentation and the expectation-maximization (EM) algorithm. We find that 426 genomes from 20 tumor types display at least one change in mutational signature activities (changepoint), and 306 genomes contain at least one of 54 recurrent changepoints shared by seven or more genomes of the same tumor type. Five recurrent changepoint locations are shared by multiple tumor types. Within these regions, the particular signature changes are often consistent across samples of the same type and some, but not all, are characterized by signatures associated with subclonal expansion. The changepoints we found cannot strictly be explained by gene density, mutation density, or cell-of-origin chromatin state. We hypothesize that they reflect a confluence of factors including evolutionary timing of mutational processes, regional differences in somatic mutation rate, large-scale changes in chromatin state that may be tissue type-specific, and changes in chromatin accessibility during subclonal expansion. These results provide insight into the regional effects of DNA damage and repair processes, and may help us localize genomic and epigenomic changes that occur during cancer development.

Competition for DNA binding between the genome protector replication protein A and the genome modifying APOBEC3 single-stranded DNA deaminases

The human APOBEC family of eleven cytosine deaminases use RNA and single-stranded DNA (ssDNA) as substrates to deaminate cytosine to uracil. This deamination event has roles in lipid metabolism by altering mRNA coding, adaptive immunity by causing evolution of antibody genes, and innate immunity through inactivation of viral genomes. These benefits come at a cost where some family members, primarily from the APOBEC3 subfamily (APOBEC3A-H, excluding E), can cause off-target deaminations of cytosine to form uracil on transiently single-stranded genomic DNA, which induces mutations that are associated with cancer evolution. Since uracil is only promutagenic, the mutations observed in cancer genomes originate only when uracil is not removed by uracil DNA glycosylase (UNG) or when the UNG-induced abasic site is erroneously repaired. However, when ssDNA is present, replication protein A (RPA) binds and protects the DNA from nucleases or recruits DNA repair proteins, such as UNG. Thus, APOBEC enzymes must compete with RPA to access their substrate. Certain APOBEC enzymes can displace RPA, bind and scan ssDNA efficiently to search for cytosines, and can become highly overexpressed in tumor cells. Depending on the DNA replication conditions and DNA structure, RPA can either be in excess or deficient. Here we discuss the interplay between these factors and how despite RPA, multiple cancer genomes have a mutation bias at cytosines indicative of APOBEC activity.

Addressing the benefits of inhibiting APOBEC3-dependent mutagenesis in cancer

Mutational signatures associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC)3 cytosine deaminase activity have been found in over half of cancer types, including some therapy-resistant and metastatic tumors. Driver mutations can occur in APOBEC3-favored sequence contexts, suggesting that mutagenesis by APOBEC3 enzymes may drive cancer evolution. The APOBEC3-mediated signatures are often detected in subclonal branches of tumor phylogenies and are acquired in cancer cell lines over long periods of time, indicating that APOBEC3 mutagenesis can be ongoing in cancer. Collectively, these and other observations have led to the proposal that APOBEC3 mutagenesis represents a disease-modifying process that could be inhibited to limit tumor heterogeneity, metastasis and drug resistance. However, critical aspects of APOBEC3 biology in cancer and in healthy tissues have not been clearly defined, limiting well-grounded predictions regarding the benefits of inhibiting APOBEC3 mutagenesis in different settings in cancer. We discuss the relevant mechanistic gaps and strategies to address them to investigate whether inhibiting APOBEC3 mutagenesis may confer clinical benefits in cancer.

Exploring ABOBEC3A and APOBEC3B substrate specificity and their role in HPV positive head and neck cancer

APOBEC3 family members are cytidine deaminases catalyzing conversion of cytidine to uracil. Many studies have established a link between APOBEC3 expression and cancer development and progression, especially APOBEC3A (A3A) and APOBEC3B (A3B). Preclinical studies with human papillomavirus positive (HPV+) head and neck squamous cell carcinoma (HNSCC) and clinical trial specimens revealed induction of A3B, but not A3A expression after demethylation. We examined the kinetic features of the cytidine deaminase activity for full length A3B and found that longer substrates and a purine at −2 position favored by A3B, whereas A3A prefers shorter substrates and an adenine or thymine at −2 position. The importance and biological significance of A3B catalytic activity rather than A3A and a preference for purine at the −2 position was also established in HPV+ HNSCCs. Our study explored factors influencing formation of A3A and A3B-related cancer mutations that are essential for understanding APOBEC3-related carcinogenesis and facilitating drug discovery.

•

A3B is upregulated after 5-AzaC treatment and related to 5-AzaC sensitivity in HPV+ HNSCC

•

Full-length A3B prefers longer substrates and a purine at −2 site biochemically

•

A3B also prefers a purine at −2 site in both HPV+ and HPV− HNSCC cells

•

A3B signature at -2 site linked to poor patient survival in HPV+ HNSCC low smokers

APOBEC3: Friend or Foe in Human Papillomavirus Infection and Oncogenesis?

Human papillomavirus (HPV) infection is a causative agent of multiple human cancers, including cervical and head and neck cancers. In these HPV-positive tumors, somatic mutations are caused by aberrant activation of DNA mutators such as members of the apolipoprotein B messenger RNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) family of cytidine deaminases. APOBEC3 proteins are most notable for their restriction of various viruses, including anti-HPV activity. However, the potential role of APOBEC3 proteins in HPV-induced cancer progression has recently garnered significant attention. Ongoing research stems from the observations that elevated APOBEC3 expression is driven by HPV oncogene expression and that APOBEC3 activity is likely a significant contributor to somatic mutagenesis in HPV-positive cancers. This review focuses on recent advances in the study of APOBEC3 proteins and their roles in HPV infection and HPV-driven oncogenesis. Further, we discuss critical gaps and unanswered questions in our understanding of APOBEC3 in virus-associated cancers.

Identification of a signature of evolutionarily conserved stress-induced mutagenesis in cancer

The clustering of mutations observed in cancer cells is reminiscent of the stress-induced mutagenesis (SIM) response in bacteria. Bacteria deploy SIM when faced with DNA double-strand breaks in the presence of conditions that elicit an SOS response. SIM employs DinB, the evolutionary precursor to human trans-lesion synthesis (TLS) error-prone polymerases, and results in mutations concentrated around DNA double-strand breaks with an abundance that decays with distance. We performed a quantitative study on single nucleotide variant calls for whole-genome sequencing data from 1950 tumors, non-inherited mutations from 129 normal samples, and acquired mutations in 3 cell line models of stress-induced adaptive mutation. We introduce statistical methods to identify mutational clusters, quantify their shapes and tease out the potential mechanism that produced them. Our results show that mutations in both normal and cancer samples are indeed clustered and have shapes indicative of SIM. Clusters in normal samples occur more often in the same genomic location across samples than in cancer suggesting loss of regulation over the mutational process during carcinogenesis. Additionally, the signatures of TLS contribute the most to mutational cluster formation in both patient samples as well as experimental models of SIM. Furthermore, a measure of cluster shape heterogeneity was associated with cancer patient survival with a hazard ratio of 5.744 (Cox Proportional Hazard Regression, 95% CI: 1.824-18.09). Our results support the conclusion that the ancient and evolutionary-conserved adaptive mutation response found in bacteria is a source of genomic instability in cancer. Biological adaptation through SIM might explain the ability of tumors to evolve in the face of strong selective pressures such as treatment and suggests that the conventional 'hit it hard' approaches to therapy could prove themselves counterproductive.

Mechanisms of APOBEC3 mutagenesis in human cancer cells

The APOBEC3 family of cytosine deaminases has been implicated in some of the most prevalent mutational signatures in cancer^1-3. However, a causal link between endogenous APOBEC3 enzymes and mutational signatures in human cancer genomes has not been established, leaving the mechanisms of APOBEC3 mutagenesis poorly understood. Here, to investigate the mechanisms of APOBEC3 mutagenesis, we deleted implicated genes from human cancer cell lines that naturally generate APOBEC3-associated mutational signatures over time⁴. Analysis of non-clustered and clustered signatures across whole-genome sequences from 251 breast, bladder and lymphoma cancer cell line clones revealed that APOBEC3A deletion diminished APOBEC3-associated mutational signatures. Deletion of both APOBEC3A and APOBEC3B further decreased APOBEC3 mutation burdens, without eliminating them. Deletion of APOBEC3B increased APOBEC3A protein levels, activity and APOBEC3A-mediated mutagenesis in some cell lines. The uracil glycosylase UNG was required for APOBEC3-mediated transversions, whereas the loss of the translesion polymerase REV1 decreased overall mutation burdens. Together, these data represent direct evidence that endogenous APOBEC3 deaminases generate prevalent mutational signatures in human cancer cells. Our results identify APOBEC3A as the main driver of these mutations, indicate that APOBEC3B can restrain APOBEC3A-dependent mutagenesis while contributing its own smaller mutation burdens and dissect mechanisms that translate APOBEC3 activities into distinct mutational signatures.

The impact of rare germline variants on human somatic mutation processes

Somatic mutations are an inevitable component of ageing and the most important cause of cancer. The rates and types of somatic mutation vary across individuals, but relatively few inherited influences on mutation processes are known. We perform a gene-based rare variant association study with diverse mutational processes, using human cancer genomes from over 11,000 individuals of European ancestry. By combining burden and variance tests, we identify 207 associations involving 15 somatic mutational phenotypes and 42 genes that replicated in an independent data set at a false discovery rate of 1%. We associate rare inherited deleterious variants in genes such as MSH3, EXO1, SETD2, and MTOR with two phenotypically different forms of DNA mismatch repair deficiency, and variants in genes such as EXO1, PAXIP1, RIF1, and WRN with deficiency in homologous recombination repair. In addition, we identify associations with other mutational processes, such as APEX1 with APOBEC-signature mutagenesis. Many of the genes interact with each other and with known mutator genes within cellular sub-networks. Considered collectively, damaging variants in the identified genes are prevalent in the population. We suggest that rare germline variation in diverse genes commonly impacts mutational processes in somatic cells.

Migrating bubble synthesis promotes mutagenesis through lesions in its template

Break-induced replication (BIR) proceeds via a migrating D-loop for hundreds of kilobases and is highly mutagenic. Previous studies identified long single-stranded (ss) nascent DNA that accumulates during leading strand synthesis to be a target for DNA damage and a primary source of BIR-induced mutagenesis. Here, we describe a new important source of mutagenic ssDNA formed during BIR: the ssDNA template for leading strand BIR synthesis formed during D-loop migration. Specifically, we demonstrate that this D-loop bottom template strand (D-BTS) is susceptible to APOBEC3A (A3A)-induced DNA lesions leading to mutations associated with BIR. Also, we demonstrate that BIR-associated ssDNA promotes an additional type of genetic instability: replication slippage between microhomologies stimulated by inverted DNA repeats. Based on our results we propose that these events are stimulated by both known sources of ssDNA formed during BIR, nascent DNA formed by leading strand synthesis, and the D-BTS that we describe here. Together we report a new source of mutagenesis during BIR that may also be shared by other homologous recombination pathways driven by D-loop repair synthesis.

APOBEC3 mutational signatures are associated with extensive and diverse genomic instability across multiple tumour types

BACKGROUND

The APOBEC3 (apolipoprotein B mRNA editing enzyme catalytic polypeptide 3) family of cytidine deaminases is responsible for two mutational signatures (SBS2 and SBS13) found in cancer genomes. APOBEC3 enzymes are activated in response to viral infection, and have been associated with increased mutation burden and TP53 mutation. In addition to this, it has been suggested that APOBEC3 activity may be responsible for mutations that do not fall into the classical APOBEC3 signatures (SBS2 and SBS13), through generation of double strand breaks.Previous work has mainly focused on the effects of APOBEC3 within individual tumour types using exome sequencing data. Here, we use whole genome sequencing data from 2451 primary tumours from 39 different tumour types in the Pan-Cancer Analysis of Whole Genomes (PCAWG) data set to investigate the relationship between APOBEC3 and genomic instability (GI).

RESULTS AND CONCLUSIONS

We found that the number of classical APOBEC3 signature mutations correlates with increased mutation burden across different tumour types. In addition, the number of APOBEC3 mutations is a significant predictor for six different measures of GI. Two GI measures (INDELs attributed to INDEL signatures ID6 and ID8) strongly suggest the occurrence and error prone repair of double strand breaks, and the relationship between APOBEC3 mutations and GI remains when SNVs attributed to kataegis are excluded.We provide evidence that supports a model of cancer genome evolution in which APOBEC3 acts as a causative factor in the development of diverse and widespread genomic instability through the generation of double strand breaks. This has important implications for treatment approaches for cancers that carry APOBEC3 mutations, and challenges the view that APOBECs only act opportunistically at sites of single stranded DNA.

Examining clustered somatic mutations with SigProfilerClusters

MOTIVATION

Clustered mutations are found in the human germline as well as in the genomes of cancer and normal somatic cells. Clustered events can be imprinted by a multitude of mutational processes, and they have been implicated in both cancer evolution and development disorders. Existing tools for identifying clustered mutations have been optimized for a particular subtype of clustered event and, in most cases, relied on a predefined inter-mutational distance (IMD) cutoff combined with a piecewise linear regression analysis.

RESULTS

Here, we present SigProfilerClusters, an automated tool for detecting all types of clustered mutations by calculating a sample-dependent IMD threshold using a simulated background model that takes into account extended sequence context, transcriptional strand asymmetries and regional mutation densities. SigProfilerClusters disentangles all types of clustered events from non-clustered mutations and annotates each clustered event into an established subclass, including the widely used classes of doublet-base substitutions, multi-base substitutions, omikli and kataegis. SigProfilerClusters outputs non-clustered mutations and clustered events using standard data formats as well as provides multiple visualizations for exploring the distributions and patterns of clustered mutations across the genome.

AVAILABILITY AND IMPLEMENTATION

SigProfilerClusters is supported across most operating systems and made freely available at https://github.com/AlexandrovLab/SigProfilerClusters with an extensive documentation located at https://osf.io/qpmzw/wiki/home/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Prospectively defined patterns of APOBEC3A mutagenesis are prevalent in human cancers

Mutational signatures defined by single base substitution (SBS) patterns in cancer have elucidated potential mutagenic processes that contribute to malignancy. Two prevalent mutational patterns in human cancers are attributed to the APOBEC3 cytidine deaminase enzymes. Among the seven human APOBEC3 proteins, APOBEC3A is a potent deaminase and proposed driver of cancer mutagenesis. In this study, we prospectively examine genome-wide aberrations by expressing human APOBEC3A in avian DT40 cells. From whole-genome sequencing, we detect hundreds to thousands of base substitutions per genome. The APOBEC3A signature includes widespread cytidine mutations and a unique insertion-deletion (indel) signature consisting largely of cytidine deletions. This multi-dimensional APOBEC3A signature is prevalent in human cancer genomes. Our data further reveal replication-associated mutations, the rate of stem-loop and clustered mutations, and deamination of methylated cytidines. This comprehensive signature of APOBEC3A mutagenesis is a tool for future studies and a potential biomarker for APOBEC3 activity in cancer.

Spectrum of DNA mismatch repair failures viewed through the lens of cancer genomics and implications for therapy

Genome sequencing can be used to detect DNA repair failures in tumors and learn about underlying mechanisms. Here, we synthesize findings from genomic studies that examined deficiencies of the DNA mismatch repair (MMR) pathway. The impairment of MMR results in genome-wide hypermutation and in the 'microsatellite instability' (MSI) phenotype-occurrence of indel mutations at short tandem repeat (microsatellite) loci. The MSI status of tumors was traditionally assessed by molecular testing of a selected set of MS loci or by measuring MMR protein expression levels. Today, genomic data can provide a more complete picture of the consequences on genomic instability. Multiple computational studies examined somatic mutation distributions that result from failed DNA repair pathways in tumors. These include analyzing the commonly studied trinucleotide mutational spectra of single-nucleotide variants (SNVs), as well as of other features such as indels, structural variants, mutation clusters and regional mutation rate redistribution. The identified mutation patterns can be used to rigorously measure prevalence of MMR failures across cancer types, and potentially to subcategorize the MMR deficiencies. Diverse data sources, genomic and pre-genomic, from human and from experimental models, suggest there are different ways in which MMR can fail, and/or that the cell-type or genetic background may result in different types of MMR mutational patterns. The spectrum of MMR failures may direct cancer evolution, generating particular sets of driver mutations. Moreover, MMR affects outcomes of therapy by DNA damaging drugs, antimetabolites, nonsense-mediated mRNA decay (NMD) inhibitors, and immunotherapy by promoting either resistance or sensitivity, depending on the type of therapy.

The Base-Editing Enzyme APOBEC3A Catalyzes Cytosine Deamination in RNA with Low Proficiency and High Selectivity

Human APOBEC3A (A3A) is a nucleic acid-modifying enzyme that belongs to the cytidine deaminase family. Canonically, A3A catalyzes the deamination of cytosine into uracil in single-stranded DNA, an activity that makes A3A both a critical antiviral defense factor and a useful tool for targeted genome editing. However, mutagenesis by A3A has also been readily detected in both cellular DNA and RNA, activities that have been implicated in cancer. Given the importance of substrate discrimination for the physiological, pathological, and biotechnological activities of A3A, here we explore the mechanistic basis for its preferential targeting of DNA over RNA. Using a chimeric substrate containing a target ribocytidine within an otherwise DNA backbone, we demonstrate that a single hydroxyl at the sugar of the target base acts as a major selectivity determinant for deamination. To assess the contribution of bases neighboring the target cytosine, we show that overall RNA deamination is greatly reduced relative to that of DNA but can be observed when ideal features are present, such as preferred sequence context and secondary structure. A strong dependence on idealized substrate features can also be observed with a mutant of A3A (eA3A, N57G), which has been employed for genome editing due to altered selectivity for DNA over RNA. Altogether, our work reveals a relationship between the overall decreased reactivity of A3A and increased substrate selectivity, and our results hold implications both for characterizing off-target mutagenesis and for engineering optimized DNA deaminases for base-editing technologies.

Intratumor genetic heterogeneity and clonal evolution to decode endometrial cancer progression

Analyzing different tumor regions by next generation sequencing allows the assessment of intratumor genetic heterogeneity (ITGH), a phenomenon that has been studied widely in some tumor types but has been less well explored in endometrial carcinoma (EC). In this study, we sought to characterize the spatial and temporal heterogeneity of 9 different ECs using whole-exome sequencing, and by performing targeted sequencing validation of the 42 primary tumor regions and 30 metastatic samples analyzed. In addition, copy number alterations of serous carcinomas were assessed by comparative genomic hybridization arrays. From the somatic mutations, identified by whole-exome sequencing, 532 were validated by targeted sequencing. Based on these data, the phylogenetic tree reconstructed for each case allowed us to establish the tumors’ evolution and correlate this to tumor progression, prognosis, and the presence of recurrent disease. Moreover, we studied the genetic landscape of an ambiguous EC and the molecular profile obtained was used to guide the selection of a potential personalized therapy for this patient, which was subsequently validated by preclinical testing in patient-derived xenograft models. Overall, our study reveals the impact of analyzing different tumor regions to decipher the ITGH in ECs, which could help make the best treatment decision.

Mapping clustered mutations in cancer reveals APOBEC3 mutagenesis of ecDNA

Clustered somatic mutations are common in cancer genomes and previous analyses reveal several types of clustered single-base substitutions, which include doublet- and multi-base substitutions1-5, diffuse hypermutation termed omikli6, and longer strand-coordinated events termed kataegis3,7-9. Here we provide a comprehensive characterization of clustered substitutions and clustered small insertions and deletions (indels) across 2,583 whole-genome-sequenced cancers from 30 types of cancer10. Clustered mutations were highly enriched in driver genes and associated with differential gene expression and changes in overall survival. Several distinct mutational processes gave rise to clustered indels, including signatures that were enriched in tobacco smokers and homologous-recombination-deficient cancers. Doublet-base substitutions were caused by at least 12 mutational processes, whereas most multi-base substitutions were generated by either tobacco smoking or exposure to ultraviolet light. Omikli events, which have previously been attributed to APOBEC3 activity6, accounted for a large proportion of clustered substitutions; however, only 16.2% of omikli matched APOBEC3 patterns. Kataegis was generated by multiple mutational processes, and 76.1% of all kataegic events exhibited mutational patterns that are associated with the activation-induced deaminase (AID) and APOBEC3 family of deaminases. Co-occurrence of APOBEC3 kataegis and extrachromosomal DNA (ecDNA), termed kyklonas (Greek for cyclone), was found in 31% of samples with ecDNA. Multiple distinct kyklonic events were observed on most mutated ecDNA. ecDNA containing known cancer genes exhibited both positive selection and kyklonic hypermutation. Our results reveal the diversity of clustered mutational processes in human cancer and the role of APOBEC3 in recurrently mutating and fuelling the evolution of ecDNA.

Qualitative and Quantitative Analysis of DNA Cytidine Deaminase Activity

The human genome encodes eleven DNA cytidine deaminases in the AID/APOBEC family, which encompass endogenous roles ranging from genetic diversification of the immunoglobulin locus to virus restriction. All AID/APOBEC functions are enabled by their catalyzation of cytidine deamination resulting in mutations and DNA damage. When acting aberrantly, deaminases can cause off-target mutations in the cellular genome resulting in somatic mutations, DNA damage, and genome instability. An association between cytidine deaminase-induced mutations and human cancers has been recognized over the last decade, necessitating assays for investigation of intracellular deaminase activity. Here we present two assays for deamination activity which enable in vitro evaluation of in vivo events. We define both a qualitative assay to confirm deaminase activity within cells as well as a quantitative assay for granular evaluation and comparisons of deamination activity across different cell populations or experimental conditions. The two procedures are customizable assays which can easily be adapted to individual labs and experiments.

Genome Instability in Multiple Myeloma: Facts and Factors

Multiple myeloma (MM) is a malignant neoplasm of terminally differentiated immunoglobulin-producing B lymphocytes called plasma cells. MM is the second most common hematologic malignancy, and it poses a heavy economic and social burden because it remains incurable and confers a profound disability to patients. Despite current progress in MM treatment, the disease invariably recurs, even after the transplantation of autologous hematopoietic stem cells (ASCT). Biological processes leading to a pathological myeloma clone and the mechanisms of further evolution of the disease are far from complete understanding. Genetically, MM is a complex disease that demonstrates a high level of heterogeneity. Myeloma genomes carry numerous genetic changes, including structural genome variations and chromosomal gains and losses, and these changes occur in combinations with point mutations affecting various cellular pathways, including genome maintenance. MM genome instability in its extreme is manifested in mutation kataegis and complex genomic rearrangements: chromothripsis, templated insertions, and chromoplexy. Chemotherapeutic agents used to treat MM add another level of complexity because many of them exacerbate genome instability. Genome abnormalities are driver events and deciphering their mechanisms will help understand the causes of MM and play a pivotal role in developing new therapies.