APOBEC-induced mutations in human cancers are strongly enriched on the lagging DNA strand during replication

A tale of two strands: Decoding chromatin replication through strand-specific sequencing

DNA replication, a fundamental process in all living organisms, proceeds with continuous synthesis of the leading strand by DNA polymerase ε (Pol ε) and discontinuous synthesis of the lagging strand by polymerase δ (Pol δ). This inherent asymmetry at each replication fork necessitates the development of methods to distinguish between these two nascent strands in vivo. Over the past decade, strand-specific sequencing strategies, such as enrichment and sequencing of protein-associated nascent DNA (eSPAN) and Okazaki fragment sequencing (OK-seq), have become essential tools for studying chromatin replication in eukaryotic cells. In this review, we outline the foundational principles underlying these methodologies and summarize key mechanistic insights into DNA replication, parental histone transfer, epigenetic inheritance, and beyond, gained through their applications. Finally, we discuss the limitations and challenges of current techniques, highlighting the need for further technological innovations to better understand the dynamics and regulation of chromatin replication in eukaryotic cells.

Impact of APOBEC3s on the occurrence, development and prognosis of esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a severe malignant tumor of the digestive system that poses a significant threat to human health. Despite its significance, the complex molecular mechanism regulating the occurrence and development of ESCC remain elusive. The apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3 (APOBEC3) members constitute a pivotal subfamily of the APOBEC family that possess cytidine deaminase activity. In recent years, APOBEC3s (A3s) have received increasing attention due to their pivotal roles in the occurrence, development, and prognosis of ESCC. This comprehensive review systematically summarizes the latest research progress on the mechanisms of action of A3s in ESCC and discusses their impact on the development and therapeutic considerations for ESCC, with a particular focus on their potential role in immunotherapy. These insights may be of great value in continued exploration of ESCC pathogenesis and provides a theoretical foundation for the development of clinical treatment strategies for ESCC.

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

Two APOBEC DNA cytosine deaminase enzymes, APOBEC3A and APOBEC3B, generate somatic mutations in cancer, thereby 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, in squamous cell carcinoma we find that, there is expansion of GRHL3expression and activity to a subset of cells undergoing DNA replication and concomitant extension of APOBEC3A expression to proliferating cells. These findings suggest that APOBEC3A may play a functional role during keratinocyte differentiation, and offer a mechanism for acquisition of APOBEC3A mutagenic activity in tumours.

APOBEC3B expression in 293T viral producer cells drives mutations in chimeric antigen receptors and reduces CAR T cell efficacy

Chimeric antigen receptor (CAR) T cells are a clinically approved therapy for blood cancers. To produce clinical-grade CAR T cells, a retroviral or lentiviral vector is used to deliver the CAR and associated genes to patient T cells. Apolipoprotein B editing enzyme, catalytic polypeptide 3 (APOBEC3) enzymes are known to be upregulated after transfection and retroviral infection and to deaminate cytidine to uracil in nucleic acids, resulting in cytidine-to-thymine mutations in DNA. Here, we hypothesized that APOBEC3 enzymes, induced during the production of CAR T cells, impact the efficacy of the resulting CAR T cells. We demonstrated that APOBEC3 family member APOBEC3B was upregulated at the RNA and protein levels after transfection of HEK293T cells with plasmids to make lentivirus, and that APOBEC3 signature mutations were present in the CAR construct. APOBEC3B overexpression in HEK293T cells led to further mutations in the resulting CAR T cells, and significantly decreased CAR T cell killing. APOBEC3B knockout in HEK293T cells led to reduced mutations in the CAR construct and significantly increased in CAR T cell killing. These results suggest that generation of CAR-expressing viruses from producer cell lines deficient in genome-modifying proteins such as APOBEC3B could enhance the quality of CAR T cell production.

The interferon response at the intersection of genome integrity and innate immunity

In recent years, numerous reports indicated that, besides pathogen infections, DNA replication stress and defective DNA repair can trigger the innate immune response by introducing a state of viral mimicry, due to cytosolic accumulation of the self-nucleic acid species, which culminates in the activation of type I interferon (IFN) pathway. In turn, IFN upregulates a variety of factors mutually implicated in immune- and genome-related mechanisms, shedding light on the unprecedented causality between genome stability and innate immunity. Intriguingly, in addition to being induced by replication stress, IFN-regulated factors can also promote it, pinpointing IFN signaling as both a consequence and a cause of replication stress. Here, we provide an overview of the factors and molecular mechanisms implicated in the evolutionary conserved crosstalk between genome maintenance and innate immunity, highlighting the role of the IFN-stimulated gene 15 (ISG15), which appears to be at the hub of this intersection. Moreover, we discuss the potential significance and clinical implications of the immune-mediated modulation of DNA replication and repair upon pathogen infection and in human diseases such as cancer and autoinflammatory syndromes. Finally, we discuss the relevant open questions and future directions.

Characterisation of APOBEC3B-Mediated RNA editing in breast cancer cells reveals regulatory roles of NEAT1 and MALAT1 lncRNAs

RNA editing is a crucial post-transcriptional process that influences gene expression and increases the diversity of the proteome as a result of amino acid substitution. Recently, the APOBEC3 family has emerged as a significant player in this mechanism, with APOBEC3A (A3A) having prominent roles in base editing during immune and stress responses. APOBEC3B (A3B), another family member, has gained attention for its potential role in generating genomic DNA mutations in breast cancer. In this study, we coupled an inducible expression cell model with a novel methodology for identifying differential variants in RNA (DVRs) to map A3B-mediated RNA editing sites in a breast cancer cell model. Our findings indicate that A3B engages in selective RNA editing including targeting NEAT1 and MALAT1 long non-coding RNAs that are often highly expressed in tumour cells. Notably, the binding of these RNAs sequesters A3B and suppresses global A3B activity against RNA and DNA. Release of A3B from NEAT1/MALAT1 resulted in increased A3B activity at the expense of A3A activity suggesting a regulatory feedback loop between the two family members. This research substantially advances our understanding of A3B's role in RNA editing, its mechanistic underpinnings, and its potential relevance in the pathogenesis of breast cancer.

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.

Kataegis associated mutational processes linked to adverse prostate cancer presentation in African men

Kataegis, the focal hypermutation of single base substitutions (SBS) in tumour genomes, has received little attention with respect to prostate cancer (PCa) associated molecular and clinical features. Most notably, data is lacking with regards to this tumour evolutionary phenomenon and PCa racial disparities, with African men disproportionately impacted. Here through comparison between African (n = 109) and non-African (n = 79) whole genome sequenced treatment naïve primary tumours, using a single analytical workflow we assessed for shared and unique features of kataegis. Linking kataegis to aggressive presentation, structural variant burden and copy number loss, we attributed APOBEC3 activity through higher rates of SBS2 to high-risk African tumours. While kataegis positive African patients presented with elevated prostate specific antigen levels, their tumours showed evolutionary unique trajectories marked by increased subclonal and structural variant-independent kataegis. The potential to exacerbate tumour heterogeneity emphases the significance of continued exploration of biological behaviours and environmental exposures for African patients.

DNA mismatch and damage patterns revealed by single-molecule sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other diseases1,2. Most mutations begin as nucleotide mismatches or damage in one of the two strands of the DNA before becoming double-strand mutations if unrepaired or misrepaired3,4. However, current DNA-sequencing technologies cannot accurately resolve these initial single-strand events. Here we develop a single-molecule, long-read sequencing method (Hairpin Duplex Enhanced Fidelity sequencing (HiDEF-seq)) that achieves single-molecule fidelity for base substitutions when present in either one or both DNA strands. HiDEF-seq also detects cytosine deamination-a common type of DNA damage-with single-molecule fidelity. We profiled 134 samples from diverse tissues, including from individuals with cancer predisposition syndromes, and derive from them single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumours deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples that are deficient in only polymerase proofreading. We also define a single-strand damage signature for APOBEC3A. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. As double-strand DNA mutations are only the end point of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable studies of how mutations arise in a variety of contexts, especially in cancer and ageing.

Whole genome and transcriptome integrated analyses guide clinical care of pediatric poor prognosis cancers

The role for routine whole genome and transcriptome analysis (WGTA) for poor prognosis pediatric cancers remains undetermined. Here, we characterize somatic mutations, structural rearrangements, copy number variants, gene expression, immuno-profiles and germline cancer predisposition variants in children and adolescents with relapsed, refractory or poor prognosis malignancies who underwent somatic WGTA and matched germline sequencing. Seventy-nine participants with a median age at enrollment of 8.8 y (range 6 months to 21.2 y) are included. Germline pathogenic/likely pathogenic variants are identified in 12% of participants, of which 60% were not known prior. Therapeutically actionable variants are identified by targeted gene report and whole genome in 32% and 62% of participants, respectively, and increase to 96% after integrating transcriptome analyses. Thirty-two molecularly informed therapies are pursued in 28 participants with 54% achieving a clinical benefit rate; objective response or stable disease ≥6 months. Integrated WGTA identifies therapeutically actionable variants in almost all tumors and are directly translatable to clinical care of children with poor prognosis cancers.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

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.

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.

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.

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.

Genetic inhibitors of APOBEC3B-induced mutagenesis

The cytidine deaminases APOBEC3A (A3A) and APOBEC3B (A3B) are prominent mutators of human cancer genomes. However, tumor-specific genetic modulators of APOBEC-induced mutagenesis are poorly defined. Here, we used a screen to identify 61 gene deletions that increase A3B-induced mutations in yeast. We also determined whether each deletion was epistatic with Ung1 loss, which indicated whether the encoded factors participate in the homologous recombination (HR)-dependent bypass of A3B/Ung1-dependent abasic sites or suppress A3B-catalyzed deamination by protecting against aberrant formation of single-stranded DNA (ssDNA). We found that the mutation spectra of A3B-induced mutations revealed genotype-specific patterns of strand-specific ssDNA formation and nucleotide incorporation across APOBEC-induced lesions. Combining these three metrics, we were able to establish a multifactorial signature of APOBEC-induced mutations specific to (1) failure to remove H3K56 acetylation, (2) defective CTF18-RFC complex function, and (3) defective HR-mediated bypass of APOBEC-induced lesions. We extended these results by analyzing mutation data for human tumors and found BRCA1/2-deficient breast cancers display three- to fourfold more APOBEC-induced mutations. Mirroring our results in yeast, Rev1-mediated C-to-G substitutions are mainly responsible for increased APOBEC-signature mutations in BRCA1/2-deficient tumors, and these mutations associate with lagging strand synthesis during replication. These results identify important factors that influence DNA replication dynamics and likely the abundance of APOBEC-induced mutation during tumor progression. They also highlight a novel role for BRCA1/2 during HR-dependent lesion bypass of APOBEC-induced lesions during cancer cell replication.

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.

Genome-wide maps of UVA and UVB mutagenesis in yeast reveal distinct causative lesions and mutational strand asymmetries

Ultraviolet (UV) light primarily causes C > T substitutions in lesion-forming dipyrimidine sequences. However, many of the key driver mutations in melanoma do not fit this canonical UV signature, but are instead caused by T > A, T > C, or C > A substitutions. To what extent exposure to the UVB or UVA spectrum of sunlight can induce these noncanonical mutation classes, and the molecular mechanism involved is unclear. Here, we repeatedly exposed wild-type or repair-deficient yeast (Saccharomyces cerevisiae) to UVB or UVA light and characterized the resulting mutations by whole genome sequencing. Our data indicate that UVB induces C > T and T > C substitutions in dipyrimidines, and T > A substitutions that are often associated with thymine-adenine (TA) sequences. All of these mutation classes are induced in nucleotide excision repair-deficient cells and show transcriptional strand asymmetry, suggesting they are caused by helix-distorting UV photoproducts. In contrast, UVA exposure induces orders of magnitude fewer mutations with a distinct mutation spectrum. UVA-induced mutations are elevated in Ogg1-deficient cells, and the resulting spectrum consists almost entirely of C > A/G > T mutations, indicating they are likely derived from oxidative guanine lesions. These mutations show replication asymmetry, with elevated G > T mutations on the leading strand, suggesting there is a strand bias in the removal or bypass of guanine lesions during replication. Finally, we develop a mutation reporter to show that UVA induces a G > T reversion mutation in yeast that mimics the oncogenic NRAS Q61K mutation in melanoma. Taken together, these findings indicate that UVA and UVB exposure can induce many of the noncanonical mutation classes that cause driver mutations in melanoma.

Genetic modifiers of APOBEC-induced mutagenesis

The cytidine deaminases APOBEC3A and APOBEC3B (A3B) are prominent mutators of human cancer genomes. However, tumor-specific genetic modulators of APOBEC-induced mutagenesis are poorly defined. Here, we utilized a screen to identify 61 gene deletions that increase A3B-induced mutations in yeast. Also, we determined whether each deletion was epistatic with UNG1 loss, which indicated whether the encoded factors participate in the error-free bypass of A3B/Ung1-dependent abasic sites or suppress A3B-catalyzed deamination by protecting against aberrant formation of single stranded DNA (ssDNA). Additionally, we determined that the mutation spectra of A3B-induced mutations revealed genotype-specific patterns of strand-specific ssDNA formation and nucleotide incorporation across APOBEC-induced lesions. Combining these three metrics we were able to establish a multifactorial signature of APOBEC-induced mutations specific to (1) failure to remove H3K56 acetylation, which results in extremely high A3B-induced mutagenesis, (2) defective CTF18-RFC complex function, which results in high levels of A3B induced mutations specifically on the leading strand template that synergistically increase with loss of UNG1, and (3) defective HR-mediated bypass of APOBEC-induced lesions, which were epistatic with Ung1 loss and result from increased Rev1-mediated C-to-G substitutions. We extended these results by analyzing mutation data for human tumors and found BRCA1/2-deficient breast cancer tumors display 3- to 4-fold more APOBEC-induced mutations. Mirroring our results in yeast, for BRCA1/2 deficient tumors Rev1-mediated C-to-G substitutions are solely responsible for increased APOBEC-signature mutations and these mutations occur on the lagging strand during DNA replication. Together these results identify important factors that influence the dynamics of DNA replication and likely the abundance of APOBEC-induced mutation during tumor progression as well as a novel mechanistic role for BRCA1/2 during HR-dependent lesion bypass of APOBEC-induced lesions during cancer cell replication.

APOBEC3-mediated mutagenesis in cancer: causes, clinical significance and therapeutic potential

Apolipoprotein B mRNA-editing enzyme, catalytic polypeptides (APOBECs) are cytosine deaminases involved in innate and adaptive immunity. However, some APOBEC family members can also deaminate host genomes to generate oncogenic mutations. The resulting mutations, primarily signatures 2 and 13, occur in many tumor types and are among the most common mutational signatures in cancer. This review summarizes the current evidence implicating APOBEC3s as major mutators and outlines the exogenous and endogenous triggers of APOBEC3 expression and mutational activity. The review also discusses how APOBEC3-mediated mutagenesis impacts tumor evolution through both mutagenic and non-mutagenic pathways, including by inducing driver mutations and modulating the tumor immune microenvironment. Moving from molecular biology to clinical outcomes, the review concludes by summarizing the divergent prognostic significance of APOBEC3s across cancer types and their therapeutic potential in the current and future clinical landscapes.

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.

Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other genetic diseases1-4. Almost all of these mosaic mutations begin as nucleotide mismatches or damage in only one of the two strands of the DNA prior to becoming double-strand mutations if unrepaired or misrepaired5. However, current DNA sequencing technologies cannot resolve these initial single-strand events. Here, we developed a single-molecule, long-read sequencing method that achieves single-molecule fidelity for single-base substitutions when present in either one or both strands of the DNA. It also detects single-strand cytosine deamination events, a common type of DNA damage. We profiled 110 samples from diverse tissues, including from individuals with cancer-predisposition syndromes, and define the first single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumors deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples deficient in only polymerase proofreading. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. Since the double-strand DNA mutations interrogated by prior studies are only the endpoint of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable new studies of how mutations arise in a variety of contexts, especially in cancer and aging.

Lactobacillus for ribosome peptide editing cancer

This study reviews newly discovered insect peptide point mutations as new possible cancer research targets. To interpret newly discovered peptide point mutations in insects as new possible cancer research targets, we focused on the numerous peptide changes found in the 'CSP' family on the sex pheromone gland of the female silkworm moth Bombyx mori. We predict that the Bombyx peptide modifications will have a significant effect on cancer CUP (cancers of unknown primary) therapy and that bacterial peptide editing techniques, specifically Lactobacillus combined to CRISPR, will be used to regulate ribosomes and treat cancer in humans.

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.

Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes

DNA transfer from cytoplasmic organelles to the cell nucleus is a legacy of the endosymbiotic event-the majority of nuclear-mitochondrial segments (NUMTs) are thought to be ancient, preceding human speciation^1-3. Here we analyse whole-genome sequences from 66,083 people-including 12,509 people with cancer-and demonstrate the ongoing transfer of mitochondrial DNA into the nucleus, contributing to a complex NUMT landscape. More than 99% of individuals had at least one of 1,637 different NUMTs, with 1 in 8 individuals having an ultra-rare NUMT that is present in less than 0.1% of the population. More than 90% of the extant NUMTs that we evaluated inserted into the nuclear genome after humans diverged from apes. Once embedded, the sequences were no longer under the evolutionary constraint seen within the mitochondrion, and NUMT-specific mutations had a different mutational signature to mitochondrial DNA. De novo NUMTs were observed in the germline once in every 10⁴ births and once in every 10³ cancers. NUMTs preferentially involved non-coding mitochondrial DNA, linking transcription and replication to their origin, with nuclear insertion involving multiple mechanisms including double-strand break repair associated with PR domain zinc-finger protein 9 (PRDM9) binding. The frequency of tumour-specific NUMTs differed between cancers, including a probably causal insertion in a myxoid liposarcoma. We found evidence of selection against NUMTs on the basis of size and genomic location, shaping a highly heterogenous and dynamic human NUMT landscape.

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.

APOBEC3B: Future direction of liver cancer research

Liver cancer is one of the most common cancers in the world, and the rate of liver cancer is high due to the of its illness. The main risk factor for liver cancer is infection with the hepatitis B virus (HBV), but a considerable number of genetic and epigenetic factors are also directly or indirectly involved in the underlying pathogenesis of liver cancer. In particular, the apolipoprotein B mRNA editing enzyme, catalytic peptide-like protein (APOBEC) family (DNA or mRNA editor family), which has been the focus of virology research for more than a decade, has been found to play a significant role in the occurrence and development of various cancers, providing a new direction for the research of liver cancer. APOBEC3B is a cytosine deaminase that controls a variety of biological processes, such as protein expression, innate immunity, and embryonic development, by participating in the process of cytidine deamination to uridine in DNA and RNA. In humans, APOBEC3B has long been known as a DNA editor for limiting viral replication and transcription. APOBEC3B is widely expressed at low levels in a variety of normal tissues and organs, but it is significantly upregulated in different types of tumor tissues and tumor lines. Thus, APOBEC3B has received increasing attention in various cancers, but the role of APOBEC3B in the occurrence and development of liver cancer due to infection with HBV remains unclear. This review provides a brief introduction to the pathogenesis of hepatocellular carcinoma induced by HBV, and it further explores the latest results of APOBEC3B research in the development of HBV and liver cancer, thereby providing new directions and strategies for the treatment and prevention of liver cancer.

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.

APOBEC mutagenesis is low in most types of non-B DNA structures

While somatic mutations are known to be enriched in genome regions with non-canonical DNA secondary structure, the impact of particular mutagens still needs to be elucidated. Here, we demonstrate that in human cancers, the APOBEC mutagenesis is not enriched in direct repeats, mirror repeats, short tandem repeats, and G-quadruplexes, and even decreased below its level in B-DNA for cancer samples with very high APOBEC activity. In contrast, we observe that the APOBEC-induced mutational density is positively associated with APOBEC activity in inverted repeats (cruciform structures), where the impact of cytosine at the 3’-end of the hairpin loop is substantial. Surprisingly, the APOBEC-signature mutation density per TC motif in the single-stranded DNA of a G-quadruplex (G4) is lower than in the four-stranded part of G4 and in B-DNA. The APOBEC mutagenesis, as well as the UV-mutagenesis in melanoma samples, are absent in Z-DNA regions, owing to the depletion of their mutational signature motifs.

•

APOBEC mutagenesis is not enriched in most non-canonical DNA structures

•

Inverted repeats (cruciform structures) show increased APOBEC mutagenesis

•

G-quadruplex’s unstructured strand has low APOBEC-induced mutation density

•

Decrease of APOBEC mutagenesis in non-B DNA possibly associated with PrimPol

HPV16 and HPV18 Genome Structure, Expression, and Post-Transcriptional Regulation

Human papillomaviruses (HPV) are a group of small non-enveloped DNA viruses whose infection causes benign tumors or cancers. HPV16 and HPV18, the two most common high-risk HPVs, are responsible for ~70% of all HPV-related cervical cancers and head and neck cancers. The expression of the HPV genome is highly dependent on cell differentiation and is strictly regulated at the transcriptional and post-transcriptional levels. Both HPV early and late transcripts differentially expressed in the infected cells are intron-containing bicistronic or polycistronic RNAs bearing more than one open reading frame (ORF), because of usage of alternative viral promoters and two alternative viral RNA polyadenylation signals. Papillomaviruses proficiently engage alternative RNA splicing to express individual ORFs from the bicistronic or polycistronic RNA transcripts. In this review, we discuss the genome structures and the updated transcription maps of HPV16 and HPV18, and the latest research advances in understanding RNA cis-elements, intron branch point sequences, and RNA-binding proteins in the regulation of viral RNA processing. Moreover, we briefly discuss the epigenetic modifications, including DNA methylation and possible APOBEC-mediated genome editing in HPV infections and carcinogenesis.

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.

OUP accepted manuscript

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.

CDKN1A/p21WAF1, RB1, ARID1A, FLG, and HRNR mutation patterns provide insights into urinary tract environmental exposure carcinogenesis and potential treatment strategies

Bladder carcinoma has a 6% 5-year survival-rate for metastatic disease, with poorly understood links between genetic and environmental drivers of disease development, progression, and treatment response. Rhode Island has among the highest annual age-adjusted incidence rate of bladder cancer at 26.0/100,000, compared to 20.0 in the US, with a paucity of known driver genes for targeted therapies or predictive biomarkers. Bladder carcinomas have the highest frequency of alterations in CDKN1A/p21^WAF1 (10%) across all cancer types analyzed in The Cancer Genome Atlas (TCGA) PanCancer Atlas Studies, displaying a predominance of truncating mutations (86%). We found that lung carcinomas lack CDKN1A truncating mutations, despite the shared role of tobacco as a risk factor for bladder cancer. Bladder carcinomas also have the highest frequency of RB1 alterations in TCGA (25%). We find that chromophobe renal cell carcinomas with CDKN1A and RB1 mutations are 100% truncating. Analysis of 1,868 bladder tumors demonstrated that truncating CDKN1A mutations co-occur with truncating RB1 mutations, suggesting an environmental exposure signature. Moreover, we found that HRNR and FLG mutations are enriched in tumors with CDKN1A alteration, suggesting potential novel roles in promoting bladder tumorigenesis. Association of HRNR with AKT activation offers possible therapeutic avenues, and FLG may provide insight into carcinogen exposure within the bladder. We suggest that because APOBEC mutations largely shape the bladder cancer mutational landscape, these events likely contribute to dysfunctional DNA repair genes, leading to frameshifts and the predominance of truncations in CDKN1A, RB1, ARID1A, or other drivers. We propose that patients with co-occurrence of CDKN1A and RB1 truncations may display enhanced responsiveness to targeted therapies combining cisplatin with ATR, ATM, CHK1, and CHK2 inhibitors, expanding therapeutic options for patients in need of improved precision treatments.

Single-stranded DNA binding proteins influence APOBEC3A substrate preference

The cytidine deaminase, APOBEC3A (A3A), is a prominent source of mutations in multiple cancer types. These APOBEC-signature mutations are non-uniformly distributed across cancer genomes, associating with single-stranded (ss) DNA formed during DNA replication and hairpin-forming sequences. The biochemical and cellular factors that influence these specificities are unclear. We measured A3A's cytidine deaminase activity in vitro on substrates that model potential sources of ssDNA in the cell and found that A3A is more active on hairpins containing 4 nt ssDNA loops compared to hairpins with larger loops, bubble structures, replication fork mimics, ssDNA gaps, or linear DNA. Despite pre-bent ssDNAs being expected to fit better in the A3A active site, we determined A3A favors a 4 nt hairpin substrate only 2- to fivefold over linear ssDNA substrates. Addition of whole cell lysates or purified RPA to cytidine deaminase assays more severely reduced A3A activity on linear ssDNA (45 nt) compared to hairpin substrates. These results indicate that the large enrichment of A3A-driven mutations in hairpin-forming sequences in tumor genomes is likely driven in part by other proteins that preferentially bind longer ssDNA regions, which limit A3A's access. Furthermore, A3A activity is reduced at ssDNA associated with a stalled T7 RNA polymerase, suggesting that potential protein occlusion by RNA polymerase also limits A3A activity. These results help explain the small transcriptional strand bias for APOBEC mutation signatures in cancer genomes and the general targeting of hairpin-forming sequences in the lagging strand template during DNA replication.

The origin of human mutation in light of genomic data

Despite years of active research into the role of DNA repair and replication in mutagenesis, surprisingly little is known about the origin of spontaneous human mutation in the germ line. With the advent of high-throughput sequencing, genome-scale data have revealed statistical properties of mutagenesis in humans. These properties include variation of the mutation rate and spectrum along the genome at different scales in relation to epigenomic features and dependency on parental age. Moreover, mutations originated in mothers are less frequent than mutations originated in fathers and have a distinct genomic distribution. Statistical analyses that interpret these patterns in the context of known biochemistry can provide mechanistic models of mutagenesis in humans.

Immune inactivation by APOBEC3B enrichment predicts response to chemotherapy and survival in gastric cancer

Apolipoprotein B mRNA editing enzyme catalytic polypeptide 3B (APOBEC3B) plays an important role in tumor mutagenesis. However, its clinical significance in gastric cancer (GC) remains largely unknown. We enrolled a total of 482 GC patients from Zhongshan Hospital, Fudan University for immunohistochemistry (IHC) staining to evaluate the prognostic and predictive values of APOBEC3B. Genomic and phenotypic datasets from the Cancer Genome Atlas (TCGA) and Asian Cancer Research Group (ACRG) cohort were downloaded for external validation and complementary bioinformatic analysis. Fresh specimens of additional 60 patients from Zhongshan Hospital, Fudan University were collected to detect CD8⁺ T cell phenotype with flow cytometry (FCM). The high expression of APOBEC3B indicated inferior overall survival (OS, P < .001 and P = .003) and disease-free survival (DFS, P < .001 and P < .001), yet superior therapeutic responsiveness to fluorouracil-based adjuvant chemotherapy (ACT) in TNM stage II patients. The tumor microenvironment (TME) of APOBEC3B-enriched tumors was characterized by reduced infiltration of tumor reactive CD8⁺ T cells expressing both effector molecules and immune checkpoints. APOBEC3B ^high CD8⁺ T cell ^high GC patients were most likely to benefit from ACT and PD-1 blockade. Our study demonstrates that APOBEC3B was an independent prognostic and predictive factor in GC. The potential interplay between APOBEC3B and CD8⁺ T cells merited further investigations.

APOBEC3C, a nucleolar protein induced by genotoxins, is excluded from DNA damage sites

The human genome contains 11 APOBEC (apolipoprotein B mRNA editing catalytic polypeptide‐like) cytidine deaminases classified into four families. These proteins function mainly in innate antiviral immunity and can also restrict endogenous retrotransposable element multiplication. The present study focuses on APOBEC3C (A3C), a member of the APOBEC3 subfamily. Some APOBEC3 proteins use their enzymatic activity on genomic DNA, inducing mutations and DNA damage, while other members facilitate DNA repair. Our results show that A3C is highly expressed in cells treated with DNA‐damaging agents. Its expression is regulated by p53. Depletion of A3C slightly decreases proliferation and does not affect DNA repair via homologous recombination or nonhomologous end joining. The A3C interactomes obtained from control cells and cells exposed to the genotoxin etoposide indicated that A3C is a nucleolar protein. This was confirmed by the detection of either endogenous or ectopic A3C in nucleoli. Interestingly, we show that A3C is excluded from areas of DNA breaks in live cells. Our data also indicate that the C‐terminal part of A3C is responsible for its nucleolar localization and exclusion from DNA damage sites.

APOBEC Mutagenesis Is Concordant between Tumor and Viral Genomes in HPV-Positive Head and Neck Squamous Cell Carcinoma

APOBEC is a mutagenic source in human papillomavirus (HPV)-mediated malignancies, including HPV+ oropharyngeal squamous cell carcinoma (HPV + OPSCC), and in HPV genomes. It is unknown why APOBEC mutations predominate in HPV + OPSCC, or if the APOBEC-induced mutations observed in both human cancers and HPV genomes are directly linked. We performed sequencing of host somatic exomes, transcriptomes, and HPV16 genomes from 79 HPV + OPSCC samples, quantifying APOBEC mutational burden and activity in both host and virus. APOBEC was the dominant mutational signature in somatic exomes. In viral genomes, there was a mean of five (range 0-29) mutations per genome. The mean of APOBEC mutations in viral genomes was one (range 0-5). Viral APOBEC mutations, compared to non-APOBEC mutations, were more likely to be low-variant allele fraction mutations, suggesting that APOBEC mutagenesis actively occurrs in viral genomes during infection. HPV16 APOBEC-induced mutation patterns in OPSCC were similar to those previously observed in cervical samples. Paired host and viral analyses revealed that APOBEC-enriched tumor samples had higher viral APOBEC mutation rates (p = 0.028), and APOBEC-associated RNA editing (p = 0.008), supporting the concept that APOBEC mutagenesis in host and viral genomes is directly linked and occurrs during infection. Using paired sequencing of host somatic exomes, transcriptomes, and viral genomes, we demonstrated for the first-time definitive evidence of concordance between tumor and viral APOBEC mutagenesis. This finding provides a missing link connecting APOBEC mutagenesis in host and virus and supports a common mechanism driving APOBEC dysregulation.

Interaction with the CCT chaperonin complex limits APOBEC3A cytidine deaminase cytotoxicity

The APOBEC3 cytidine deaminases are implicated as the cause of a prevalent somatic mutation pattern found in cancer genomes. The APOBEC3 enzymes act as viral restriction factors by mutating viral genomes. Mutation of the cellular genome is presumed to be an off-target activity of the enzymes, although the regulatory measures for APOBEC3 expression and activity remain undefined. It is therefore difficult to predict circumstances that enable APOBEC3 interaction with cellular DNA that leads to mutagenesis. The APOBEC3A (A3A) enzyme is the most potent deaminase of the family. Using proteomics, we evaluate protein interactors of A3A to identify potential regulators. We find that A3A interacts with the chaperonin-containing TCP-1 (CCT) complex, a cellular machine that assists in protein folding and function. Importantly, depletion of CCT results in A3A-induced DNA damage and cytotoxicity. Evaluation of cancer genomes demonstrates an enrichment of A3A mutational signatures in cancers with silencing mutations in CCT subunit genes. Together, these data suggest that the CCT complex interacts with A3A, and that disruption of CCT function results in increased A3A mutational activity.

Replication catastrophe induced by cyclic hypoxia leads to increased APOBEC3B activity

Tumor heterogeneity includes variable and fluctuating oxygen concentrations, which result in the accumulation of hypoxic regions in most solid tumors. Tumor hypoxia leads to increased therapy resistance and has been linked to genomic instability. Here, we tested the hypothesis that exposure to levels of hypoxia that cause replication stress could increase APOBEC activity and the accumulation of APOBEC-mediated mutations. APOBEC-dependent mutational signatures have been well-characterized, although the physiological conditions which underpin them have not been described. We demonstrate that fluctuating/cyclic hypoxic conditions which lead to replication catastrophe induce the expression and activity of APOBEC3B. In contrast, stable/chronic hypoxic conditions which induce replication stress in the absence of DNA damage are not sufficient to induce APOBEC3B. Most importantly, the number of APOBEC-mediated mutations in patient tumors correlated with a hypoxia signature. Together, our data support the conclusion that hypoxia-induced replication catastrophe drives genomic instability in tumors, specifically through increasing the activity of APOBEC3B.

Action-at-a-distance mutations at 5'-GpA-3' sites induced by oxidized guanine in WRN-knockdown cells

G:C sites distant from 8-oxo-7,8-dihydroguanine (G O, 8-hydroxyguanine) are frequently mutated when the lesion-bearing plasmid DNA is replicated in human cells with reduced Werner syndrome (WRN) protein. To detect the untargeted mutations preferentially, the oxidized guanine base was placed downstream of the reporter supF gene and the plasmid DNA was introduced into WRN-knockdown cells. The total mutant frequency seemed higher in the WRN-knockdown cells as compared to the control cells. Mutation analyses revealed that substitution mutations occurred at the G:C pairs of 5'-GpA-3'/5'-TpC-3' sites, the preferred sequence for the apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3 (APOBEC3)-family cytosine deaminases, in the supF gene in both control and knockdown cells. These mutations were observed more frequently at G sites than C sites on the DNA strand where the G O base was originally located. This tendency was promoted by the knockdown of the WRN protein. The present results imply the possible involvement of APOBEC3-family cytosine deaminases in the action-at-a-distance (untargeted) mutations at G:C (or G) sites induced by G O and in cancer initiation by oxidative stress.

Analysis of gene expression and mutation data points on contribution of transcription to the mutagenesis by APOBEC enzymes

Since the discovery of the role of the APOBEC enzymes in human cancers, the mechanisms of this type of mutagenesis remain little understood. Theoretically, targeting of single-stranded DNA by the APOBEC enzymes could occur during cellular processes leading to the unwinding of DNA double-stranded structure. Some evidence points to the importance of replication in the APOBEC mutagenesis, while the role of transcription is still underexplored. Here, we analyzed gene expression and whole genome sequencing data from five types of human cancers with substantial APOBEC activity to estimate the involvement of transcription in the APOBEC mutagenesis and compare its impact with that of replication. Using the TCN motif as the mutation signature of the APOBEC enzymes, we observed a correlation of active APOBEC mutagenesis with gene expression, confirmed the increase of APOBEC-induced mutations in early-replicating regions and estimated the relative impact of transcription and replication on the APOBEC mutagenesis. We also found that the known effect of higher density of APOBEC-induced mutations on the lagging strand was highest in middle-replicating regions and observed higher APOBEC mutation density on the sense strand, the latter bias positively correlated with the gene expression level.

A3G-induced mutations show a low prevalence and exhibit plus-strand regional distribution in hepatitis B virus DNA from patients with non-hepatocellular carcinoma (HCC) and HCC

APOBEC3G (A3G) cytidine deaminase is an innate immune restriction factor that can edit and inhibit hepatitis B virus (HBV) replication. The preferred target of A3G is deamination of the third cytosine of 5'CCC to form a mutant marker 5'CC C → K. However, the distribution of A3G-induced mutations on HBV DNA during infection is not well characterized. To provide clarity, we obtained the HBV DNA sequences from HBV infected individuals with and without hepatocellular carcinoma (HCC and non-HCC, respectively), from the NCBI database, and calculated the r values of A3G-induced 5'CC C → K mutation prevalence in HBV DNA. A3G-induced mutations were weakly prevalent and mainly distributed in the plus strand of HBV DNA (r = 1.407). The mutations on the minus strand were weaker (r = .8189). There were A3G-induced mutation regions in the 1200 to 2000 nt region of the plus strand and the 1600 to 1500 nt region of the minus strand. There was no significant difference in the r values of A3G-induced mutations in HBV DNA between the HCC and non-HCC groups. However, the rvalue of the plus strand 2400 to 2800 nt regions of HCC derived HBV DNA (r = 4.2) was significantly higher than that of the same regions of non-HCC derived HBV DNA (r = 1.21). These findings clarify the weak prevalence and preferred plus-strand distribution of A3G-induced mutations on HBV DNA from HCC and non-HCC. These findings may provide valuable clues regarding the interaction mechanism between A3G and HBV DNA and inform HCC screening.

DNA Methylation, Deamination, and Translesion Synthesis Combine to Generate Footprint Mutations in Cancer Driver Genes in B-Cell Derived Lymphomas and Other Cancers

Cancer genomes harbor numerous genomic alterations and many cancers accumulate thousands of nucleotide sequence variations. A prominent fraction of these mutations arises as a consequence of the off-target activity of DNA/RNA editing cytosine deaminases followed by the replication/repair of edited sites by DNA polymerases (pol), as deduced from the analysis of the DNA sequence context of mutations in different tumor tissues. We have used the weight matrix (sequence profile) approach to analyze mutagenesis due to Activation Induced Deaminase (AID) and two error-prone DNA polymerases. Control experiments using shuffled weight matrices and somatic mutations in immunoglobulin genes confirmed the power of this method. Analysis of somatic mutations in various cancers suggested that AID and DNA polymerases η and θ contribute to mutagenesis in contexts that almost universally correlate with the context of mutations in A:T and G:C sites during the affinity maturation of immunoglobulin genes. Previously, we demonstrated that AID contributes to mutagenesis in (de)methylated genomic DNA in various cancers. Our current analysis of methylation data from malignant lymphomas suggests that driver genes are subject to different (de)methylation processes than non-driver genes and, in addition to AID, the activity of pols η and θ contributes to the establishment of methylation-dependent mutation profiles. This may reflect the functional importance of interplay between mutagenesis in cancer and (de)methylation processes in different groups of genes. The resulting changes in CpG methylation levels and chromatin modifications are likely to cause changes in the expression levels of driver genes that may affect cancer initiation and/or progression.

HMCES Maintains Replication Fork Progression and Prevents Double-Strand Breaks in Response to APOBEC Deamination and Abasic Site Formation

5-Hydroxymethylcytosine (5hmC) binding, ES-cell-specific (HMCES) crosslinks to apurinic or apyrimidinic (AP, abasic) sites in single-strand DNA (ssDNA). To determine whether HMCES responds to the ssDNA abasic site in cells, we exploited the activity of apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3A (APOBEC3A). APOBEC3A preferentially deaminates cytosines to uracils in ssDNA, which are then converted to abasic sites by uracil DNA glycosylase. We find that HMCES-deficient cells are hypersensitive to nuclear APOBEC3A localization. HMCES relocalizes to chromatin in response to nuclear APOBEC3A and protects abasic sites from processing into double-strand breaks (DSBs). Abasic sites induced by APOBEC3A slow both leading and lagging strand synthesis, and HMCES prevents further slowing of the replication fork by translesion synthesis (TLS) polymerases zeta (Polζ) and kappa (Polκ). Thus, our study provides direct evidence that HMCES responds to ssDNA abasic sites in cells to prevent DNA cleavage and balance the engagement of TLS polymerases.

Mehta et al. use APOBEC3A to demonstrate that HMCES responds to ssDNA abasic sites in cells and prevents replication fork collapse. APOBEC3A-induced abasic sites slow both leading and lagging strand polymerization, and HMCES engagement prevents further fork slowing because of the action of TLS polymerases zeta (Polζ) and kappa (Polκ).

An extended APOBEC3A mutation signature in cancer

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.