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.

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.

APOBEC2 safeguards skeletal muscle cell fate through binding chromatin and regulating transcription of non-muscle genes during myoblast differentiation.

The apolipoprotein B messenger RNA editing enzyme, catalytic polypeptide (APOBEC) family is composed of nucleic acid editors with roles ranging from antibody diversification to RNA editing. APOBEC2, a member of this family with an evolutionarily conserved nucleic acid-binding cytidine deaminase domain, has neither an established substrate nor function. Using a cellular model of muscle differentiation where APOBEC2 is inducibly expressed, we confirmed that APOBEC2 does not have the attributed molecular functions of the APOBEC family, such as RNA editing, DNA demethylation, and DNA mutation. Instead, we found that during muscle differentiation APOBEC2 occupied a specific motif within promoter regions; its removal from those regions resulted in transcriptional changes. Mechanistically, these changes reflect the direct interaction of APOBEC2 with histone deacetylase (HDAC) transcriptional corepressor complexes. We also found that APOBEC2 could bind DNA directly, in a sequence-specific fashion, suggesting that it functions as a recruiter of HDAC to specific genes whose promoters it occupies. These genes are normally suppressed during muscle cell differentiation, and their suppression may contribute to the safeguarding of muscle cell fate. Altogether, our results reveal a unique role for APOBEC2 within the APOBEC family.

Unraveling C-to-U RNA editing events from direct RNA sequencing.

In mammals, RNA editing events involve the conversion of adenosine (A) in inosine (I) by ADAR enzymes or the hydrolytic deamination of cytosine (C) in uracil (U) by the APOBEC family of enzymes, mostly APOBEC1. RNA editing has a plethora of biological functions, and its deregulation has been associated with various human disorders. While the large-scale detection of A-to-I is quite straightforward using the Illumina RNAseq technology, the identification of C-to-U events is a non-trivial task. This difficulty arises from the rarity of such events in eukaryotic genomes and the challenge of distinguishing them from background noise. Direct RNA sequencing by Oxford Nanopore Technology (ONT) permits the direct detection of Us on sequenced RNA reads. Surprisingly, using ONT reads from wild-type (WT) and APOBEC1-knock-out (KO) murine cell lines as well as in vitro synthesized RNA without any modification, we identified a systematic error affecting the accuracy of the Cs call, thereby leading to incorrect identifications of C-to-U events. To overcome this issue in direct RNA reads, here we introduce a novel machine learning strategy based on the isolation Forest (iForest) algorithm in which C-to-U editing events are considered as sequencing anomalies. Using in vitro synthesized and human ONT reads, our model optimizes the signal-to-noise ratio improving the detection of C-to-U editing sites with high accuracy, over 90% in all samples tested. Our results suggest that iForest, known for its rapid implementation and minimal memory requirements, is a promising tool to denoise ONT reads and reliably identify RNA modifications.

Regulatory variants of APOBEC3 genes potentially associate with COVID-19 severity in populations with African ancestry.

Since November 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has caused the worldwide pandemic of the coronavirus disease 2019 (COVID-19), the impact of which is huge to the lives of world populations. Many studies suggested that such situation will continue due to the endless mutations in SARS-CoV-2 genome that result in complexity of the efforts for the control of SARS-CoV-2, since the special enrichment of nucleotide substitution C>U in SARS-CoV-2 sequences were discovered mainly due to the editing by human host factors APOBEC3 genes. The observation of SARS-CoV-2 variants Beta (B.1.351) and Omicron (B.1.1.529) firstly spreading in South Africa promoted us to hypothesize that genetic variants of APOBEC3 special in African populations may be attributed to the higher mutation rate of SARS-CoV-2 variants in Africa. Current study was conducted to search for functional variants of APOBEC3 genes associate with COVID-19 hospitalization in African population. By integrating data from the 1000 Genomes Project, Genotype-Tissue Expression (GTEx), and Host Genetics Initiative (HGI) of COVID-19, we identified potential functional SNPs close to APOBEC3 genes that are associated with COVID-19 hospitalization in African but not with other populations. Our study provides new insights on the potential contribution of APOBEC3 genes on the evolution of SARS-CoV-2 mutations in African population, but further replication is needed to confirm our results.

Proportion of APOBEC3-induced defective HIV DNA after 1†year of dolutegravir + lamivudine simplification in the ANRS 167 LAMIDOL trial.

BACKGROUND: Hypermutated viruses induced by APOBEC3 (apolipoprotein B mRNA-editing, enzyme-catalytic, polypeptide-like 3) proteins comprise some of the defective viruses in the HIV reservoir. Here, we assessed the proportion of APOBEC3-induced defective proviruses in HIV-positive patients before and after receiving dolutegravir + lamivudine dual therapy. METHODS: PBMCs of virologically suppressed patients enrolled in the ANRS 167 LAMIDOL trial, evaluating a switch from triple therapy to dolutegravir + lamivudine, were collected 8 weeks before (W-8) and 48 weeks after (W48) dual-therapy initiation. The Vif and RT regions were subject to next-generation sequencing. Bioinformatic algorithms were developed to identify APOBEC3-defective sequences and APOBEC3-related drug resistance mutations (APOMuts). All hypermutated sequences and those containing at least one stop codon were considered as defective. RESULTS: One hundred and four patients were enrolled (median virological suppression duration: 4.2 years; IQR: 2.0-9.1). Proviral defective reads at W-8 and W48 were detected in Vif in 22% and 29% of patients, respectively, and in RT in 38% and 42% of patients, respectively. At least one APOMut was present in proviruses of 27% and 38% of patients at W-8 and W48, respectively. The ratio of APOMuts/number of potential APOMut sites was significantly higher at W48 (16.5%) than at W-8 (9.8%, P = 0.007). The presence of APOBEC3-defective viruses at W-8 was not associated with HIV total DNA level, nor with the third drug class received prior to switching to dolutegravir + lamivudine, nor with the duration of virological suppression. CONCLUSIONS: Whereas no significant change in the proportion of patients with APOBEC3-defective proviruses was evidenced after 1 year of dolutegravir + lamivudine maintenance, enrichment in APOMuts was observed. Further longer-term studies are needed to assess the other forms of defective viruses with dual-therapy.

APOBEC3 deaminase editing in mpox virus as evidence for sustained human transmission since at least 2016.

Historically, mpox has been characterized as an endemic zoonotic disease that transmits through contact with the reservoir rodent host in West and Central Africa. However, in May 2022, human cases of mpox were detected spreading internationally beyond countries with known endemic reservoirs. When the first cases from 2022 were sequenced, they shared 42 nucleotide differences from the closest mpox virus (MPXV) previously sampled. Nearly all these mutations are characteristic of the action of APOBEC3 deaminases, host enzymes with antiviral function. Assuming APOBEC3 editing is characteristic of human MPXV infection, we developed a dual-process phylogenetic molecular clock that-inferring a rate of ~6 APOBEC3 mutations per year-estimates that MPXV has been circulating in humans since 2016. These observations of sustained MPXV transmission present a fundamental shift to the perceived paradigm of MPXV epidemiology as a zoonosis and highlight the need for revising public health messaging around MPXV as well as outbreak management and control.

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.

Sequencing of monkeypox virus from infected patients reveals viral genomes with APOBEC3-like editing, gene inactivation, and bacterial agents of skin superinfection.

A large outbreak of Monkeypox virus (MPXV) infections has arisen in May 2022 in nonendemic countries. Here, we performed DNA metagenomics using next-generation sequencing with Illumina or Nanopore technologies for clinical samples from MPXV-infected patients diagnosed between June and July 2022. Classification of the MPXV genomes and determination of their mutational patterns were performed using Nextclade. Twenty-five samples from 25 patients were studied. A MPXV genome was obtained for 18 patients, essentially from skin lesions and rectal swabbing. All 18 genomes were classified in clade IIb, lineage B.1, and we identified four B.1 sublineages (B.1.1, B.1.10, B.1.12, B.1.14). We detected a high number of mutations (range, 64-73) relatively to a 2018 Nigerian genome (genome GenBank Accession no. NC_063383.1), which were harbored by a large part of a set of 3184 MPXV genomes of lineage B.1 recovered from GenBank and Nextstrain; and we detected 35 mutations relatively to genome ON563414.3 (a B.1 lineage reference genome). Nonsynonymous mutations occurred in genes encoding central proteins, among which transcription factors and core and envelope proteins, and included two mutations that would truncate a RNA polymerase subunit and a phospholipase d-like protein, suggesting an alternative start codon and gene inactivation, respectively. A large majority (94%) of nucleotide substitutions were G > A or C > U, suggesting the action of human APOBEC3 enzymes. Finally, >1000 reads were identified as from Staphylococcus aureus and Streptococcus pyogenes for 3 and 6 samples, respectively. These findings warrant a close genomic monitoring of MPXV to get a better picture of the genetic micro-evolution and mutational patterns of this virus, and a close clinical monitoring of skin bacterial superinfection in monkeypox patients.

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.

An APOBEC3 Mutational Signature in the Genomes of Human-Infecting Orthopoxviruses.

The ongoing worldwide monkeypox outbreak is caused by viral lineages (globally referred to as hMPXV1) that are related to but distinct from clade IIb MPXV viruses transmitted within Nigeria. Analysis of the genetic differences has indicated that APOBEC-mediated editing might be responsible for the unexpectedly high number of mutations observed in hMPXV1 genomes. Here, using 1,624 publicly available hMPXV1 sequences, we analyzed the mutations that accrued between 2017 and the emergence of the current predominant variant (B.1), as well as those that that have been accumulating during the 2022 outbreak. We confirmed an overwhelming prevalence of C-to-T and G-to-A mutations, with a sequence context (5'-TC-3') consistent with the preferences of several human APOBEC3 enzymes. We also found that mutations preferentially occur in highly expressed viral genes, although no transcriptional asymmetry was observed. A comparison of the mutation spectrum and context was also performed against the human-specific variola virus (VARV) and the zoonotic cowpox virus (CPXV), as well as fowlpox virus (FWPV). The results indicated that in VARV genomes, C-to-T and G-to-A changes were more common than the opposite substitutions, although the effect was less marked than for hMPXV1. Conversely, no preference toward C-to-T and G-to-A changes was observed in CPXV and FWPV. Consistently, the sequence context of C-to-T changes confirmed a preference for a T in the -1 position for VARV, but not for CPXV or FWPV. Overall, our results strongly support the view that, irrespective of the transmission route, orthopoxviruses infecting humans are edited by the host APOBEC3 enzymes. IMPORTANCE Analysis of the viral lineages responsible for the 2022 monkeypox outbreak suggested that APOBEC enzymes are driving hMPXV1 evolution. Using 1,624 public sequences, we analyzed the mutations that accumulated between 2017 and the emergence of the predominant variant and those that characterize the last outbreak. We found that the mutation spectrum of hMPXV1 has been dominated by TC-to-TT and GA-to-AA changes, consistent with the editing activity of human APOBEC3 proteins. We also found that mutations preferentially affect highly expressed viral genes, possibly because transcription exposes single-stranded DNA (ssDNA), a target of APOBEC3 editing. Notably, analysis of the human-specific variola virus (VARV) and the zoonotic cowpox virus (CPXV) indicated that in VARV genomes, TC-to-TT and GA-to-AA changes are likewise extremely frequent. Conversely, no preference toward TC-to-TT and GA-to-AA changes is observed in CPXV. These results suggest that APOBEC3 proteins have an impact on the evolution of different human-infecting orthopoxviruses.

Antiretroviral APOBEC3 cytidine deaminases alter HIV-1 provirus integration site profiles.

APOBEC3 (A3) proteins are host-encoded deoxycytidine deaminases that provide an innate immune barrier to retroviral infection, notably against HIV-1. Low levels of deamination are believed to contribute to the genetic evolution of HIV-1, while intense catalytic activity of these proteins can induce catastrophic hypermutation in proviral DNA leading to near-total HIV-1 restriction. So far, little is known about how A3 cytosine deaminases might impact HIV-1 proviral DNA integration sites in human chromosomal DNA. Using a deep sequencing approach, we analyze the influence of catalytic active and inactive APOBEC3F and APOBEC3G on HIV-1 integration site selections. Here we show that DNA editing is detected at the extremities of the long terminal repeat regions of the virus. Both catalytic active and non-catalytic A3 mutants decrease insertions into gene coding sequences and increase integration sites into SINE elements, oncogenes and transcription-silencing non-B DNA features. Our data implicates A3 as a host factor influencing HIV-1 integration site selection and also promotes what appears to be a more latent expression profile.

Clinical Implications of APOBEC3-Mediated Mutagenesis in Breast Cancer.

Over recent years, members of the APOBEC3 family of cytosine deaminases have been implicated in increased cancer genome mutagenesis, thereby contributing to intratumor and intertumor genomic heterogeneity and therapy resistance in, among others, breast cancer. Understanding the available methods for clinical detection of these enzymes, the conditions required for their (dysregulated) expression, the clinical impact they have, and the clinical implications they may offer is crucial in understanding the current impact of APOBEC3-mediated mutagenesis in breast cancer. Here, we provide a comprehensive review of recent developments in the detection of APOBEC3-mediated mutagenesis and responsible APOBEC3 enzymes, summarize the pathways that control their expression, and explore the clinical ramifications and opportunities they pose. We propose that APOBEC3-mediated mutagenesis can function as a helpful predictive biomarker in several standard-of-care breast cancer treatment plans and may be a novel target for treatment.

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.

APOBEC Mutational Signatures in Hormone Receptor-Positive Human Epidermal Growth Factor Receptor 2-Negative Breast Cancers Are Associated With Poor Outcomes on CDK4/6 Inhibitors and Endocrine Therapy.

PURPOSE: APOBEC mutagenesis underlies somatic evolution and accounts for tumor heterogeneity in several cancers, including breast cancer (BC). In this study, we evaluated the characteristics of a real-world cohort for time-to-treatment discontinuation (TTD) and overall survival on CDK4/6 inhibitors (CDK4/6i) plus endocrine therapy (ET) and immune checkpoint inhibitors. METHODS: Comprehensive genomic profiling results from 29,833 BC samples were analyzed for tumor mutational burden and APOBEC signatures. For clinical outcomes, a deidentified nationwide (United States-based) BC Clinico-Genomic Database (CGDB) was evaluated with log-rank and Cox models. Patients with hormone receptor-positive (HR+) human epidermal growth factor receptor 2-negative (HER2-) BC who received first-line ET and CDK4/6i were included. Eligible patients from Mayo Clinic and Duke University were HR+ HER2- BC with sequencing data between September 2013 and July 2020. RESULTS: Of 29,833 samples sequenced, 7.9% were APOBEC+ with a high rate in invasive lobular carcinoma (16.7%) and in metastatic tumors (9.7%) relative to locally biopsied BC (4.3%; P < .001). In CGDB, 857 patients with HR+ HER2- BC received ET plus CDK4/6i in the first line. APOBEC+ patients had significantly shorter TTD on ET plus CDK4/6i than APOBEC- patients, 7.8 (95% CI, 4.3 to 14.6) versus 12.4 months (95% CI, 11.2 to 14.1; hazard ratio, 1.6; 95% CI, 1.03 to 2.39; P = .0036). Clinical benefit to immune checkpoint inhibitors was observed in HR+ HER2-, APOBEC+, tumor mutational burden-high patients, with four of nine CGDB patients (TTD 0.3-11.3 months) and four of six patients in Duke/Mayo cohorts (TTD 0.9-40.5 months) with a TTD of ≥ 3 months. CONCLUSION: APOBEC+ HR+ HER2- patients had shorter TTD on first-line ET plus CDK4/6i relative to APOBEC- patients. Further research is needed to optimize the treatment of APOBEC+ HR+ HER2- BC and to investigate the efficacy of immunotherapeutic strategies in this population.

Differences in integration frequencies and APOBEC3 profiles of five high-risk HPV types adheres to phylogeny.

Persistent infection with Human Papillomavirus (HPV) is responsible for almost all cases of cervical cancers, and HPV16 and HPV18 associated with the majority of these. These types differ in the proportion of viral minor nucleotide variants (MNVs) caused by APOBEC3 mutagenesis as well as integration frequencies. Whether these traits extend to other types remains uncertain. This study aimed to investigate and compare genomic variability and chromosomal integration in the two phylogenetically distinct Alpha-7 and Alpha-9 clades of carcinogenic HPV types. The TaME-seq protocol was employed to sequence cervical cell samples positive for HPV31, HPV33 or HPV45 and combine these with data from a previous study on HPV16 and HPV18. APOBEC3 mutation signatures were found in Alpha-9 (HPV16/31/33) but not in Alpha-7 (HPV18/45). HPV45 had significantly more MNVs compared to the other types. Alpha-7 had higher integration frequency compared to Alpha-9. An increase in integration frequency with increased diagnostic severity was found for Alpha-7. The results highlight important differences and broaden our understanding of the molecular mechanisms behind cervical cancer induced by high-risk HPV types from the Alpha-7 and Alpha-9 clades.

The roles of APOBEC-mediated RNA editing in SARS-CoV-2 mutations, replication and fitness.

During COVID-19 pandemic, mutations of SARS-CoV-2 produce new strains that can be more infectious or evade vaccines. Viral RNA mutations can arise from misincorporation by RNA-polymerases and modification by host factors. Analysis of SARS-CoV-2 sequence from patients showed a strong bias toward C-to-U mutation, suggesting a potential mutational role by host APOBEC cytosine deaminases that possess broad anti-viral activity. We report the first experimental evidence demonstrating that APOBEC3A, APOBEC1, and APOBEC3G can edit on specific sites of SARS-CoV-2 RNA to produce C-to-U mutations. However, SARS-CoV-2 replication and viral progeny production in Caco-2 cells are not inhibited by the expression of these APOBECs. Instead, expression of wild-type APOBEC3 greatly promotes viral replication/propagation, suggesting that SARS-CoV-2 utilizes the APOBEC-mediated mutations for fitness and evolution. Unlike the random mutations, this study suggests the predictability of all possible viral genome mutations by these APOBECs based on the UC/AC motifs and the viral genomic RNA structure.