An APOBEC3A hypermutation signature is distinguishable from the signature of background mutagenesis by APOBEC3B in human cancers

The Impact of Sugar Conformation on the Single-Stranded DNA Selectivity of APOBEC3A and APOBEC3B Enzymes

The APOBEC3 family of polynucleotide cytidine deaminases has diverse roles as viral restriction factors and oncogenic mutators. These enzymes convert cytidine to uridine in single-stranded (ss)DNA, inducing genomic mutations that promote drug resistance and tumor heterogeneity. Of the seven human APOBEC3 members, APOBEC3A (A3A) and APOBEC3B (A3B) are most implicated in driving pro-tumorigenic mutations. How these enzymes engage and selectively deaminate ssDNA over RNA is not well understood. We previously conducted molecular dynamics (MD) simulations that support the role of sugar conformation as a key molecular determinant in nucleic acid recognition by A3B. We hypothesize that A3A and A3B selectively deaminate substrates in the 2'-endo (DNA) conformation and show reduced activity for 3'-endo (RNA) conformation substrates. Consequently, we have characterized A3A- and A3B-binding and deaminase activity with chimeric oligonucleotides containing cytidine analogues that promote either the 2'-endo or 3'-endo conformation. Using fluorescence polarization and gel-based deamination assays, we determined that sugar conformation preferentially impacts the ability of these enzymes to deaminate substrates and less so binding to substrates. Using MD simulations, we identify specific active site interactions that promote selectivity based on the 2'-endo conformation. These findings help inform the biological functions of A3A and A3B in providing antiviral innate immunity and pathogenic functions in cancer.

Prolonged persistence of mutagenic DNA lesions in somatic cells

DNA is subject to continual damage, leaving each cell with thousands of individual DNA lesions at any given moment1-3. The efficiency of DNA repair means that most known classes of lesion have a half-life of minutes to hours3,4, but the extent to which DNA damage can persist for longer durations remains unknown. Here, using high-resolution phylogenetic trees from 89 donors, we identified mutations arising from 818 DNA lesions that persisted across multiple cell cycles in normal human stem cells from blood, liver and bronchial epithelium5-12. Persistent DNA lesions occurred at increased rates, with distinctive mutational signatures, in donors exposed to tobacco or chemotherapy, suggesting that they can arise from exogenous mutagens. In haematopoietic stem cells, persistent DNA lesions, probably from endogenous sources, generated the characteristic mutational signature SBS1913; occurred steadily throughout life, including in utero; and endured for 2.2 years on average, with 15-25% of lesions lasting at least 3 years. We estimate that on average, a haematopoietic stem cell has approximately eight such lesions at any moment in time, half of which will generate a mutation with each cell cycle. Overall, 16% of mutations in blood cells are attributable to SBS19, and similar proportions of driver mutations in blood cancers exhibit this signature. These data indicate the existence of a family of DNA lesions that arise from endogenous and exogenous mutagens, are present in low numbers per genome, persist for months to years, and can generate a substantial fraction of the mutation burden of somatic cells.

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.

Genetic alteration of mRNA editing enzyme APOBEC3B in the pathogenesis of ovarian endometriosis

RESEARCH QUESTION

What are the specific genetic alterations and associated network in endometriotic cells responsible for the disease pathogenesis?

DESIGN

Case control experimental study involving 45 women with endometriosis who underwent laparoscopic surgery (case) and 45 normal samples from women undergoing total abdominal hysterectomy (control). The endometrial samples were subjected to whole exome sequencing (WES) of endometriotic tissue and copy number variation analysis. Validation of gene hits were obtained from WES using polymerase chain reaction techniques, immunological techniques, in-silico tools and transgenic cell line models.

RESULTS

Germline heterozygous deletion of mRNA editing enzyme subunit APOBEC3B was identified in about 96% of endometriosis samples. The presence of germline deletion was confirmed with blood, endometrium and normal ovary samples obtained from the same patient. APOBEC3B deletions resulted in a hybrid protein that activates A1CF. APOBEC3B deletion can be a major cause of changes in the endometriotic microenvironment, and contributes to the pathogenesis and manifestation of the disease. The effect of APOBEC3B deletion was proved by in-vitro experiments in a cell line model, which displayed endometriosis-like characteristics. APOBEC3B germline deletion plays a major role in the pathogenesis of endometriosis, which is evident by the activation of A1CF, an increase in epithelial to mesenchymal transition, cellular proliferation, inflammation markers and a decrease in apoptosis markers.

CONCLUSION

The deleterious effects caused by APOBEC3B deletion in endometriosis were identified and confirmed. These results might provide a base for identifying the complete pathogenetic mechanism of endometriosis, thereby moving a step closer to better diagnosis and treatment options.

An overview of the functions and mechanisms of APOBEC3A in tumorigenesis

The APOBEC3 (A3) family plays a pivotal role in the immune system by performing DNA/RNA single-strand deamination. Cancers mostly arise from the accumulation of chronic mutations in somatic cells, and recent research has highlighted the A3 family as a major contributor to tumor-associated mutations, with A3A being a key driver gene leading to cancer-related mutations. A3A helps to defend the host against virus-induced tumors by editing the genome of cancer-associated viruses that invade the host. However, when it is abnormally expressed, it leads to persistent, chronic mutations in the genome, thereby fueling tumorigenesis. Notably, A3A is prominently expressed in innate immune cells, particularly macrophages, thereby affecting the functional state of tumor-infiltrating immune cells and tumor growth. Furthermore, the expression of A3A in tumor cells may directly affect their proliferation and migration. A growing body of research has unveiled that A3A is closely related to various cancers, which signifies the potential significance of A3A in cancer therapy. This paper mainly classifies and summarizes the evidence of the relationship between A3A and tumorigenesis based on the potential mechanisms, aiming to provide valuable references for further research on the functions of A3A and its development in the area of cancer therapy.

APOBEC3A drives metastasis of high-grade serous ovarian cancer by altering epithelial-to-mesenchymal transition

High-grade serous ovarian cancer (HGSOC) is the most prevalent and aggressive histological subtype of ovarian cancer, and often presents with metastatic disease. The drivers of metastasis in HGSOC remain enigmatic. APOBEC3A (A3A), an enzyme that generates mutations across various cancers, has been proposed as a mediator of tumor heterogeneity and disease progression. However, the role of A3A in HGSOC has not been explored. Through analysis of genome sequencing from primary HGSOC, we observed an association between high levels of APOBEC3 mutagenesis and poor overall survival. We experimentally addressed this correlation by modeling A3A activity in HGSOC cell lines and mouse models which resulted in increased metastatic behavior of HGSOC cells in culture and distant metastatic spread in vivo . A3A activity in both primary and cultured HGSOC cells yielded consistent alterations in expression of epithelial-mesenchymal-transition (EMT) genes resulting in hybrid EMT and mesenchymal signatures, and providing a mechanism for their increased metastatic potential. Our findings define the prevalence of A3A mutagenesis in HGSOC and implicate A3A as a driver of HGSOC metastasis via EMT, underscoring its clinical relevance as a potential prognostic biomarker. Our study lays the groundwork for the development of targeted therapies aimed at mitigating the deleterious impact of A3A-driven EMT in HGSOC.

Epigenetic control of Topoisomerase 1 activity presents a cancer vulnerability

DNA transactions introduce torsional constraints that pose an inherent risk to genome integrity. While topoisomerase 1 (TOP1) activity is essential for removing DNA supercoiling, aberrant stabilization of TOP1:DNA cleavage complexes (TOP1ccs) can result in cytotoxic DNA lesions. What protects genomic hot spots of topological stress from aberrant TOP1 activity remains unknown. Here, we identify chromatin context as an essential means to coordinate TOP1cc resolution. Through its ability to bind poly(ADP-ribose) (PAR), a protein modification required for TOP1cc repair, the histone variant macroH2A1.1 establishes a TOP1-permissive chromatin environment, while the alternatively spliced macroH2A1.2 isoform is unable to bind PAR or protect from TOP1ccs. By visualizing transcription-induced topological stress in single cells, we find that macroH2A1.1 facilitates PAR-dependent recruitment of the TOP1cc repair effector XRCC1 to protect from ssDNA damage. Impaired macroH2A1.1 splicing, a frequent cancer feature, was predictive of increased sensitivity to TOP1 poisons in a pharmaco-genomic screen in breast cancer cells, and macroH2A1.1 inactivation mirrored this effect. Consistent with this, low macroH2A1.1 expression correlated with improved survival in cancer patients treated with TOP1 inhibitors. We propose that macroH2A1 alternative splicing serves as an epigenetic modulator of TOP1-associated genome maintenance and a potential cancer vulnerability.

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 (U's) and function in antiviral innate immunity and RNA editing and can cause hypermutation in chromosomes. The resulting U's can be directly replicated, resulting in C to T mutations, or U-DNA glycosylase can convert the U's 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 depend 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 (WT) 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 U's in ssDNA from Ung1, thereby facilitating C to T mutagenesis through the accurate copying of U's by replicative DNA polymerases. Unexpectedly, we also found that for U's shielded from Ung1 by WT RPA, the mutagenic outcome is reduced in the presence of translesion DNA polymerase zeta.

Somatic mutations in 3929 HPV positive cervical cells associated with infection outcome and HPV type

Invasive cervical cancers (ICC), caused by HPV infections, have a heterogeneous molecular landscape. We investigate the detection, timing, and HPV type specificity of somatic mutations in 3929 HPV-positive exfoliated cervical cell samples from individuals undergoing cervical screening in the U.S. using deep targeted sequencing in ICC cases, precancers, and HPV-positive controls. We discover a subset of hotspot mutations rare in controls (2.6%) but significantly more prevalent in precancers, particularly glandular precancer lesions (10.2%), and cancers (25.7%), supporting their involvement in ICC carcinogenesis. Hotspot mutations differ by HPV type, and HPV18/45-positive ICC are more likely to have multiple hotspot mutations compared to HPV16-positive ICC. The proportion of cells containing hotspot mutations is higher (i.e., higher variant allele fraction) in ICC and mutations are detectable up to 6 years prior to cancer diagnosis. Our findings demonstrate the feasibility of using exfoliated cervical cells for detection of somatic mutations as potential diagnostic biomarkers.

Next generation APOBEC3 inhibitors: Optimally designed for potency and nuclease stability

APOBEC3 (or A3) enzymes have emerged as potential therapeutic targets due to their role in introducing heterogeneity in viruses and cancer, often leading to drug resistance. Inhibiting these enzymes has remained elusive as initial phosphodiester (PO) linked DNA based inhibitors lack stability and potency. We have enhanced both potency and nuclease stability, of 2'-deoxy-zebularine (dZ), substrate-based oligonucleotide inhibitors for two critical A3's: A3A and A3G. While replacing the phosphate backbone with phosphorothioate (PS) linkages increased nuclease stability, fully PS-modified inhibitors lost potency (1.4-3.7 fold) due to the structural constraints of the active site. For both enzymes, mixed PO/PS backbones enhanced potency (2.3-9.2 fold), while also vastly improving nuclease resistance. We also strategically introduced 2'-fluoro sugar modifications, creating the first nanomolar inhibitor of A3G-CTD2. With hairpin-structured inhibitors containing optimized PS patterns and LNA sugar modifications, we characterize the first single-digit nanomolar inhibitor targeting A3A. These extremely potent A3A inhibitors, were highly resistant to nuclease degradation in serum stability assays. Overall, our optimally designed A3 oligonucleotide inhibitors show improved potency and stability, compared to previous attempts to inhibit these critical enzymes, opening the door to realize the therapeutic potential of A3 inhibition.

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.

Spatial intra-tumour heterogeneity and treatment-induced genomic evolution in oesophageal adenocarcinoma: implications for prognosis and therapy

BACKGROUND

Oesophageal adenocarcinoma (OAC) is a highly heterogeneous cancer with poor survival. Standard curative treatment is chemotherapy with or without radiotherapy followed by oesophagectomy. Genomic heterogeneity is a feature of OAC and has been linked to treatment resistance.

METHODS

Whole-genome sequencing data from 59 treatment-naïve and 18 post-treatment samples from 29 OAC patients was analysed. Twenty-seven of these were enrolled in the DOCTOR trial, sponsored by the Australasian Gastro-Intestinal Trials Group. Two biopsies from each treatment-naïve tumour were assessed to define 'shared' (between both samples) and 'private' (present in one sample) mutations.

RESULTS

Mutational signatures SBS2/13 (APOBEC) and SBS3 (BRCA) were almost exclusively detected in private mutation populations of treatment-naïve tumours. Patients presenting these signatures had significantly worse disease specific survival. Furthermore, mutational signatures associated with platinum-based chemotherapy treatment as well as high platinum enrichment scores were only detected in post-treatment samples. Additionally, clones with high putative neoantigen binding scores were detected in some treatment-naïve samples suggesting immunoediting of clones.

CONCLUSIONS

This study demonstrates the high intra-tumour heterogeneity in OAC, as well as indicators for treatment-induced changes during tumour evolution. Intra-tumour heterogeneity remains a problem for successful treatment strategies in OAC.

The role of distinct APOBEC/ADAR mRNA levels in mutational signatures linked to aging and ultraviolet radiation

The APOBEC/AID family is known for its mutator activity, and recent evidence also supports the potential impact of ADARs. Furthermore, the mutator impacts of APOBEC/ADAR mutations have not yet been investigated. Assessment of pancancer TCGA exomes identified enriched somatic variants among exomes with nonsynonymous APOBEC1, APOBEC3B, APOBEC3C, ADAR, and ADARB1 mutations, compared to exomes with synonymous ones. Principal component (PC) analysis reduced the number of potential players to eight in cancer exomes/genomes, and to five in cancer types. Multivariate regression analysis was used to assess the impact of the PCs on each COSMIC mutational signature among pancancer exomes/genomes and particular cancers, identifying several novel links, including SBS17b, SBS18, and ID7 mainly determined by APOBEC1 mRNA levels; SBS40, ID1, and ID2 by age; SBS3 and SBS16 by APOBEC3A/APOBEC3B mRNA levels; ID5 and DBS9 by DNA repair/replication (DRR) defects; and SBS7a-d, SBS38, ID4, ID8, ID13, and DBS1 by ultraviolet (UV) radiation/ADARB1 mRNA levels. APOBEC/ADAR mutations appeared to potentiate the impact of DRR defects on several mutational signatures, and some factors seemed to inversely affect certain signatures. These findings potentially implicate certain APOBEC/ADAR mutations/mRNA levels in distinct mutational signatures, particularly APOBEC1 mRNA levels in aging-related signatures and ADARB1 mRNA levels in UV radiation-related signatures.

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.

The genomic landscape of Vk*MYC myeloma highlights shared pathways of transformation between mice and humans

Multiple myeloma (MM) is a heterogeneous disease characterized by frequent MYC translocations. Sporadic MYC activation in the germinal center of genetically engineered Vk*MYC mice is sufficient to induce plasma cell tumors in which a variety of secondary mutations are spontaneously acquired and selected over time. Analysis of 119 Vk*MYC myeloma reveals recurrent copy number alterations, structural variations, chromothripsis, driver mutations, apolipoprotein B mRNA-editing enzyme, catalytic polypeptide (APOBEC) mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identify frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC. In summary, here we credential the Vk*MYC mouse as a unique resource to explore MM genomic evolution and describe a fully annotated collection of diverse and immortalized murine MM tumors.

Changes in intrahost genetic diversity according to lesion severity in longitudinal HPV16 samples

Human papillomavirus type 16 (HPV16) is the most common cause of cervical cancer, but most infections are transient with lesions not progressing to cancer. There is a lack of specific biomarkers for early cancer risk stratification. This study aimed to explore the intrahost HPV16 genomic variation in longitudinal samples from HPV16-infected women with different cervical lesion severity (normal, low-grade, and high-grade). The TaME-seq deep sequencing protocol was used to generate whole genome HPV16 sequences of 102 samples collected over time from 40 individuals. Single nucleotide variants (SNVs) and intrahost SNVs (iSNVs) were identified in the viral genomes. A majority of individuals had a unique set of SNVs and these SNVs were stable over time. Overall, the number of iSNVs and APOBEC3-induced iSNVs were significantly lower in high-grade relative to normal and low-grade samples. A significant increase in the number of APOBEC3-induced iSNVs over time was observed for normal samples when compared to high-grade. Our results indicates that the lower incidence of iSNVs and APOBEC3-induced iSNVs in high-grade lesions may have implications for novel biomarkers discoveries, potentially aiding early stratification of HPV-induced cervical precancerous lesions.

Whole-genome mapping of APOBEC mutagenesis in metastatic urothelial carcinoma identifies driver hotspot mutations and a novel mutational signature

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) enzymes mutate specific DNA sequences and hairpin-loop structures, challenging the distinction between passenger and driver hotspot mutations. Here, we characterized 115 whole genomes of metastatic urothelial carcinoma (mUC) to identify APOBEC mutagenic hotspot drivers. APOBEC-associated mutations were detected in 92% of mUCs and were equally distributed across the genome, while APOBEC hotspot mutations (ApoHMs) were enriched in open chromatin. Hairpin loops were frequent targets of didymi (twins in Greek), two hotspot mutations characterized by the APOBEC SBS2 signature, in conjunction with an uncharacterized mutational context (Ap[C>T]). Next, we developed a statistical framework that identified ApoHMs as drivers in coding and non-coding genomic regions of mUCs. Our results and statistical framework were validated in independent cohorts of 23 non-metastatic UCs and 3,744 samples of 17 metastatic cancers, identifying cancer-type-specific drivers. Our study highlights the role of APOBEC in cancer development and may contribute to developing novel targeted therapy options for APOBEC-driven cancers.

Targeting APOBECs in cancer: It's about timing

APOBEC3 cytidine deaminases have emerged as key drivers of mutagenesis in a wide spectrum of tumor types and are now appreciated to play a causal role in driving tumor evolution and drug resistance. As efforts to develop APOBEC3 inhibitors progress, understanding the timing and consequences of APOBEC3-mediated mutagenesis in distinct clinical contexts will be critical for guiding the development of anti-cancer therapeutic strategies.

APOBEC shapes tumor evolution and age at onset of lung cancer in smokers

APOBEC enzymes are part of the innate immunity and are responsible for restricting viruses and retroelements by deaminating cytosine residues1,2. Most solid tumors harbor different levels of somatic mutations attributed to the off-target activities of APOBEC3A (A3A) and/or APOBEC3B (A3B)3-6. However, how APOBEC3A/B enzymes shape the tumor evolution in the presence of exogenous mutagenic processes is largely unknown. Here, by combining deep whole-genome sequencing with multi-omics profiling of 309 lung cancers from smokers with detailed tobacco smoking information, we identify two subtypes defined by low (LAS) and high (HAS) APOBEC mutagenesis. LAS are enriched for A3B-like mutagenesis and KRAS mutations, whereas HAS for A3A-like mutagenesis and TP53 mutations. Unlike APOBEC3A, APOBEC3B expression is strongly associated with an upregulation of the base excision repair pathway. Hypermutation by unrepaired A3A and tobacco smoking mutagenesis combined with TP53-induced genomic instability can trigger senescence7, apoptosis8, and cell regeneration9, as indicated by high expression of pulmonary healing signaling pathway, stemness markers and distal cell-of-origin in HAS. The expected association of tobacco smoking variables (e.g., time to first cigarette) with genomic/epigenomic changes are not observed in HAS, a plausible consequence of frequent cell senescence or apoptosis. HAS have more neoantigens, slower clonal expansion, and older age at onset compared to LAS, particularly in heavy smokers, consistent with high proportions of newly generated, unmutated cells and frequent immuno-editing. These findings show how heterogeneity in mutational burden across co-occurring mutational processes and cell types contributes to tumor development, with important clinical implications.

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.

Action-at-a-distance mutations induced by 8-oxo-7,8-dihydroguanine are dependent on APOBEC3

DNA oxidation is a serious threat to genome integrity and is involved in mutations and cancer initiation. The G base is most frequently damaged, and 8-oxo-7,8-dihydroguanine (GO, 8-hydroxyguanine) is one of the predominant damaged bases. In human cells, GO causes a G:C→T:A transversion mutation at the modified site, and also induces untargeted substitution mutations at the G bases of 5'-GpA-3' dinucleotides (action-at-a-distance mutations). The 5'-GpA-3' sequences are complementary to the 5'-TpC-3' sequences, the preferred substrates for apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3 (APOBEC3) cytosine deaminases, and thus their contribution to mutagenesis has been considered. In this study, APOBEC3B, the most abundant APOBEC3 protein in human U2OS cells, was knocked down in human U2OS cells, and a GO-shuttle plasmid was then transfected into the cells. The action-at-a-distance mutations were reduced to ~25% by the knockdown, indicating that GO-induced action-at-a-distance mutations are highly dependent on APOBEC3B in this cell line.

DEEPOMICS FFPE, a deep neural network model, identifies DNA sequencing artifacts from formalin fixed paraffin embedded tissue with high accuracy

Formalin-fixed, paraffin-embedded (FFPE) tissue specimens are routinely used in pathological diagnosis, but their large number of artifactual mutations complicate the evaluation of companion diagnostics and analysis of next-generation sequencing data. Identification of variants with low allele frequencies is challenging because existing FFPE filtering tools label all low-frequency variants as artifacts. To address this problem, we aimed to develop DEEPOMICS FFPE, an AI model that can classify a true variant from an artifact. Paired whole exome sequencing data from fresh frozen and FFPE samples from 24 tumors were obtained from public sources and used as training and validation sets at a ratio of 7:3. A deep neural network model with three hidden layers was trained with input features using outputs of the MuTect2 caller. Contributing features were identified using the SHapley Additive exPlanations algorithm and optimized based on training results. The performance of the final model (DEEPOMICS FFPE) was compared with those of existing models (MuTect filter, FFPolish, and SOBDetector) by using well-defined test datasets. We found 41 discriminating properties for FFPE artifacts. Optimization of property quantification improved the model performance. DEEPOMICS FFPE removed 99.6% of artifacts while maintaining 87.1% of true variants, with an F1-score of 88.3 in the entire dataset not used for training, which is significantly higher than those of existing tools. Its performance was maintained even for low-allele-fraction variants with a specificity of 0.995, suggesting that it can be used to identify subclonal variants. Different from existing methods, DEEPOMICS FFPE identified most of the sequencing artifacts in the FFPE samples while retaining more of true variants, including those of low allele frequencies. The newly developed tool DEEPOMICS FFPE may be useful in designing capture panels for personalized circulating tumor DNA assay and identifying candidate neoepitopes for personalized vaccine design. DEEPOMICS FFPE is freely available on the web ( http://deepomics.co.kr/ffpe ) for research.

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.

Inhibition of Uracil DNA Glycosylase Alters Frequency and Spectrum of Action-at-a-Distance Mutations Induced by 8-Oxo-7,8-dihydroguanine

The generation of DNA damage causes mutations and consequently cancer. Reactive oxygen species are important sources of DNA damage and some mutation signatures found in human cancers. 8-Oxo-7,8-dihydroguanine (GO, 8-hydroxyguanine) is one of the most abundant oxidized bases and induces a G→T transversion mutation at the modified site. The damaged G base also causes untargeted base substitution mutations at the G bases of 5'-GpA-3' dinucleotides (action-at-a-distance mutations) in human cells, and the cytosine deaminase apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3 (APOBEC3) is involved in the mutation process. The deaminated cytosine, i.e., uracil, bases are expected to be removed by uracil DNA glycosylase. Most of the substitution mutations at the G bases of 5'-GpA-3' might be caused by abasic sites formed by the glycosylase. In this study, we expressed the uracil DNA glycosylase inhibitor from Bacillus subtilis bacteriophage PBS2 in human U2OS cells and examined the effects on the GO-induced action-at-a-distance mutations. The inhibition of uracil DNA glycosylase increased the mutation frequency, and in particular, the frequency of G→A transitions. These results indicated that uracil DNA glycosylase, in addition to APOBEC3, is involved in the untargeted mutation process induced by GO.

The role of APOBEC3B in lung tumor evolution and targeted cancer therapy resistance

In this study, the impact of the apolipoprotein B mRNA-editing catalytic subunit-like (APOBEC) enzyme APOBEC3B (A3B) on epidermal growth factor receptor (EGFR)-driven lung cancer was assessed. A3B expression in EGFR mutant (EGFRmut) non-small-cell lung cancer (NSCLC) mouse models constrained tumorigenesis, while A3B expression in tumors treated with EGFR-targeted cancer therapy was associated with treatment resistance. Analyses of human NSCLC models treated with EGFR-targeted therapy showed upregulation of A3B and revealed therapy-induced activation of nuclear factor kappa B (NF-κB) as an inducer of A3B expression. Significantly reduced viability was observed with A3B deficiency, and A3B was required for the enrichment of APOBEC mutation signatures, in targeted therapy-treated human NSCLC preclinical models. Upregulation of A3B was confirmed in patients with NSCLC treated with EGFR-targeted therapy. This study uncovers the multifaceted roles of A3B in NSCLC and identifies A3B as a potential target for more durable responses to targeted cancer therapy.

Molecular landscapes of oral cancers of unknown etiology

The incidence of the mobile tongue cancer in young patients has been rising. This oral cancer (OC) type has no identified risk factors (NIRF), no established molecular markers and is not yet recognized as a distinct clinical entity. To understand this emerging malignancy, we innovatively analyzed the public head and neck cancer multi-omics data. We identified mutational signatures that successfully stratified 307 OC and 109 laryngeal cancer cases according to their clinico-pathological characteristics. The NIRF OCs exhibited significantly increased activities of endogenous clock-like and APOBEC-associated mutagenesis, alongside specific cancer driver gene mutations, distinct methylome patterns and prominent antimicrobial transcriptomic responses. Furthermore, we show that mutational signature SBS16 in OCs reflects the combined effects of alcohol drinking and tobacco smoking. Our study characterizes the unique disease histories and molecular programs of the NIRF OCs revealing that this emerging cancer subtype is likely driven by increased endogenous mutagenesis correlated with responses to microbial insults.

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.

Myelodysplasia after clonal hematopoiesis with APOBEC3-mediated CYBB inactivation in retroviral gene therapy for X-CGD

Stem cell gene therapy using the MFGS-gp91phox retroviral vector was performed on a 27-year-old patient with X-linked chronic granulomatous disease (X-CGD) in 2014. The patient's refractory infections were resolved, whereas the oxidase-positive neutrophils disappeared within 6 months. Thirty-two months after gene therapy, the patient developed myelodysplastic syndrome (MDS), and vector integration into the MECOM locus was identified in blast cells. The vector integration into MECOM was detectable in most myeloid cells at 12 months after gene therapy. However, the patient exhibited normal hematopoiesis until the onset of MDS, suggesting that MECOM transactivation contributed to clonal hematopoiesis, and the blast transformation likely arose after the acquisition of additional genetic lesions. In whole-genome sequencing, the biallelic loss of the WT1 tumor suppressor gene, which occurred immediately before tumorigenesis, was identified as a potential candidate genetic alteration. The provirus CYBB cDNA in the blasts contained 108 G-to-A mutations exclusively in the coding strand, suggesting the occurrence of APOBEC3-mediated hypermutations during the transduction of CD34-positive cells. A hypermutation-mediated loss of oxidase activity may have facilitated the survival and proliferation of the clone with MECOM transactivation. Our data provide valuable insights into the complex mechanisms underlying the development of leukemia in X-CGD gene therapy.

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.

Mutational processes of tobacco smoking and APOBEC activity generate protein-truncating mutations in cancer genomes

Mutational signatures represent a genomic footprint of endogenous and exogenous mutational processes through tumor evolution. However, their functional impact on the proteome remains incompletely understood. We analyzed the protein-coding impact of single-base substitution (SBS) signatures in 12,341 cancer genomes from 18 cancer types. Stop-gain mutations (SGMs) (i.e., nonsense mutations) were strongly enriched in SBS signatures of tobacco smoking, APOBEC cytidine deaminases, and reactive oxygen species. These mutational processes alter specific trinucleotide contexts and thereby substitute serines and glutamic acids with stop codons. SGMs frequently affect cancer hallmark pathways and tumor suppressors such as TP53, FAT1, and APC. Tobacco-driven SGMs in lung cancer correlate with smoking history and highlight a preventable determinant of these harmful mutations. APOBEC-driven SGMs are enriched in YTCA motifs and associate with APOBEC3A expression. Our study exposes SGM expansion as a genetic mechanism by which endogenous and carcinogenic mutational processes directly contribute to protein loss of function, oncogenesis, and tumor heterogeneity.

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.

Structure-guided inhibition of the cancer DNA-mutating enzyme APOBEC3A

The normally antiviral enzyme APOBEC3A is an endogenous mutagen in human cancer. Its single-stranded DNA C-to-U editing activity results in multiple mutagenic outcomes including signature single-base substitution mutations (isolated and clustered), DNA breakage, and larger-scale chromosomal aberrations. APOBEC3A inhibitors may therefore comprise a unique class of anti-cancer agents that work by blocking mutagenesis, slowing tumor evolvability, and preventing detrimental outcomes such as drug resistance and metastasis. Here we reveal the structural basis of competitive inhibition of wildtype APOBEC3A by hairpin DNA bearing 2'-deoxy-5-fluorozebularine in place of the cytidine in the TC substrate motif that is part of a 3-nucleotide loop. In addition, the structural basis of APOBEC3A's preference for YTCD motifs (Y = T, C; D = A, G, T) is explained. The nuclease-resistant phosphorothioated derivatives of these inhibitors have nanomolar potency in vitro and block APOBEC3A activity in human cells. These inhibitors may be useful probes for studying APOBEC3A activity in cellular systems and leading toward, potentially as conjuvants, next-generation, combinatorial anti-mutator and anti-cancer therapies.

Benchmarking variational AutoEncoders on cancer transcriptomics data

Deep generative models, such as variational autoencoders (VAE), have gained increasing attention in computational biology due to their ability to capture complex data manifolds which subsequently can be used to achieve better performance in downstream tasks, such as cancer type prediction or subtyping of cancer. However, these models are difficult to train due to the large number of hyperparameters that need to be tuned. To get a better understanding of the importance of the different hyperparameters, we examined six different VAE models when trained on TCGA transcriptomics data and evaluated on the downstream tasks of cluster agreement with cancer subtypes and survival analysis. We studied the effect of the latent space dimensionality, learning rate, optimizer, initialization and activation function on the quality of subsequent downstream tasks on the TCGA samples. We found β-TCVAE and DIP-VAE to have a good performance, on average, despite being more sensitive to hyperparameters selection. Based on these experiments, we derived recommendations for selecting the different hyperparameters settings. To ensure generalization, we tested all hyperparameter configurations on the GTEx dataset. We found a significant correlation (ρ = 0.7) between the hyperparameter effects on clustering performance in the TCGA and GTEx datasets. This highlights the robustness and generalizability of our recommendations. In addition, we examined whether the learned latent spaces capture biologically relevant information. Hereto, we measured the correlation and mutual information of the different representations with various data characteristics such as gender, age, days to metastasis, immune infiltration, and mutation signatures. We found that for all models the latent factors, in general, do not uniquely correlate with one of the data characteristics nor capture separable information in the latent factors even for models specifically designed for disentanglement.

Mutational signatures and their association with survival and gene expression in urological carcinomas

Different sources of mutagenesis cause consistently identifiable patterns of mutations and mutational signatures that mirror the various carcinogenetic processes. We used publicly available data from the Cancer Genome Atlas to evaluate the associations between the activity of the mutational signatures and various survival endpoints in six types of urological cancers after adjusting for established prognostic factors. The predictive power of the signatures was evaluated with dynamic area under curve models. In addition, links between mutational signature activities and differences in gene expression patterns were analysed. APOBEC-related signature SBS2 was associated with improved overall survival (OS) and disease-specific survival (DSS) in bladder carcinomas in the multivariate analysis, while clock-like signature SBS1 predicted shortened DSS and progression-free interval (PFI) in clear cell renal cell carcinomas (ccRCC). In papillary renal cell carcinomas (pRCC), SBS45 was a predictor of improved outcomes, and APOBEC-related SBS13 was a predictor of worse outcomes. Gene expression analyses revealed various enriched pathways between the low- and high-signature groups. Interestingly, in both the ccRCC and pRCC cohorts, the genes of several members of the melanoma antigen (MAGE) family were highly upregulated in the signatures, which predicted poor outcomes, and downregulated in signatures, which were associated with improved survival. To summarize, SBS signatures provide substantial prognostic value compared with just the traditional prognostic factors in certain cancer types. APOBEC-related SBS2 and SBS13 seem to provide robust prognostic information for particular urological cancers, maybe driven by the expression of specific groups of genes, including the MAGE gene family.

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.

Mutational signature assignment heterogeneity is widespread and can be addressed by ensemble approaches

Single-base substitution (SBS) mutational signatures have become standard practice in cancer genomics. In lieu of de novo signature extraction, reference signature assignment allows users to estimate the activities of pre-established SBS signatures within individual malignancies. Several tools have been developed for this purpose, each with differing methodologies. However, due to a lack of standardization, there may be inter-tool variability in signature assignment. We deeply characterized three assignment strategies and five SBS signature assignment tools. We observed that assignment strategy choice can significantly influence results and interpretations. Despite varying recommendations by tools, Refit performed best by reducing overfitting and maximizing reconstruction of the original mutational spectra. Even after uniform application of Refit, tools varied remarkably in signature assignments both qualitatively (Jaccard index = 0.38-0.83) and quantitatively (Kendall tau-b = 0.18-0.76). This phenomenon was exacerbated for 'flat' signatures such as the homologous recombination deficiency signature SBS3. An ensemble approach (EnsembleFit), which leverages output from all five tools, increased SBS3 assignment accuracy in BRCA1/2-deficient breast carcinomas. After generating synthetic mutational profiles for thousands of pan-cancer tumors, EnsembleFit reduced signature activity assignment error 15.9-24.7% on average using Catalogue of Somatic Mutations In Cancer and non-standard reference signature sets. We have also released the EnsembleFit web portal (https://www.ensemblefit.pittlabgenomics.com) for users to generate or download ensemble-based SBS signature assignments using any strategy and combination of tools. Overall, we show that signature assignment heterogeneity across tools and strategies is non-negligible and propose a viable, ensemble solution.

Glycidamide-induced hypermutation in yeast single-stranded DNA reveals a ubiquitous clock-like mutational motif in humans

Mutagens often prefer specific nucleotides or oligonucleotide motifs that can be revealed by studying the hypermutation spectra in single-stranded (ss) DNA. We utilized a yeast model to explore mutagenesis by glycidamide, a simple epoxide formed endogenously in humans from the environmental toxicant acrylamide. Glycidamide caused ssDNA hypermutation in yeast predominantly in cytosines and adenines. The most frequent mutations in adenines occurred in the nAt→nGt trinucleotide motif. Base substitutions A→G in this motif relied on Rev1 translesion polymerase activity. Inactivating Rev1 did not alter the nAt trinucleotide preference, suggesting it may be an intrinsic specificity of the chemical reaction between glycidamide and adenine in the ssDNA. We found this mutational motif enriched in published sequencing data from glycidamide-treated mouse cells and ubiquitous in human cancers. In cancers, this motif was positively correlated with the single base substitution (SBS) smoking-associated SBS4 signature, with the clock-like signatures SBS1, SBS5, and was strongly correlated with smoking history and with age of tumor donors. Clock-like feature of the motif was also revealed in cells of human skin and brain. Given its pervasiveness, we propose that this mutational motif reflects mutagenic lesions to adenines in ssDNA from a potentially broad range of endogenous and exogenous agents.

Unraveling the Enzyme-Substrate Properties for APOBEC3A-Mediated RNA Editing

The APOBEC3 family of human cytidine deaminases is involved in various cellular processes, including the innate and acquired immune system, mostly through inducing C-to-U in single-stranded DNA and/or RNA mutations. Although recent studies have examined RNA editing by APOBEC3A (A3A), its intracellular target specificity are not fully characterized. To address this gap, we performed in-depth analysis of cellular RNA editing using our recently developed sensitive cell-based fluorescence assay. Our findings demonstrate that A3A and an A3A-loop1-containing APOBEC3B (A3B) chimera are capable of RNA editing. We observed that A3A prefers to edit specific RNA substrates which are not efficiently deaminated by other APOBEC members. The editing efficiency of A3A is influenced by the RNA sequence contexts and distinct stem-loop secondary structures. Based on the identified RNA specificity features, we predicted potential A3A-editing targets in the encoding region of cellular mRNAs and discovered novel RNA transcripts that are extensively edited by A3A. Furthermore, we found a trend of increased synonymous mutations at the sites for more efficient A3A-editing, indicating evolutionary adaptation to the higher editing rate by A3A. Our results shed light on the intracellular RNA editing properties of A3A and provide insights into new RNA targets and potential impact of A3A-mediated RNA editing.

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.

Distinct mutational processes shape selection of MHC class I and class II mutations across primary and metastatic tumors

Disruption of antigen presentation via loss of major histocompatibility complex (MHC) expression is a strategy whereby cancer cells escape immune surveillance and develop resistance to immunotherapy. Here, we develop the personalized genomics algorithm Hapster and accurately call somatic mutations within the MHC genes of 10,001 primary and 2,199 metastatic tumors, creating a catalog of 1,663 non-synonymous mutations that provide key insights into MHC mutagenesis. We find that MHC class I genes are among the most frequently mutated genes in both primary and metastatic tumors, while MHC class II mutations are more restricted. Recurrent deleterious mutations are found within haplotype- and cancer-type-specific hotspots associated with distinct mutational processes. Functional classification of MHC residues reveals significant positive selection for mutations disruptive to the B2M, peptide, and T cell binding interfaces, as well as to MHC chaperones.

Sequence dependencies and mutation rates of localized mutational processes in cancer

BACKGROUND

Cancer mutations accumulate through replication errors and DNA damage coupled with incomplete repair. Individual mutational processes often show nucleotide sequence and functional region preferences. As a result, some sequence contexts mutate at much higher rates than others, with additional variation found between functional regions. Mutational hotspots, with recurrent mutations across cancer samples, represent genomic positions with elevated mutation rates, often caused by highly localized mutational processes.

METHODS

We count the 11-mer genomic sequences across the genome, and using the PCAWG set of 2583 pan-cancer whole genomes, we associate 11-mers with mutational signatures, hotspots of single nucleotide variants, and specific genomic regions. We evaluate the mutation rates of individual and combined sets of 11-mers and derive mutational sequence motifs.

RESULTS

We show that hotspots generally identify highly mutable sequence contexts. Using these, we show that some mutational signatures are enriched in hotspot sequence contexts, corresponding to well-defined sequence preferences for the underlying localized mutational processes. This includes signature 17b (of unknown etiology) and signatures 62 (POLE deficiency), 7a (UV), and 72 (linked to lymphomas). In some cases, the mutation rate and sequence preference increase further when focusing on certain genomic regions, such as signature 62 in transcribed regions, where the mutation rate is increased up to 9-folds over cancer type and mutational signature average.

CONCLUSIONS

We summarize our findings in a catalog of localized mutational processes, their sequence preferences, and their estimated mutation rates.

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.

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.