Cited by 117RNA regulation and diseaseSARS-CoV-2 and COVID-19 ResearchEvolution and Genetic Dynamics

The DOI matches its Crossref record; 69 of the 71 references it cites were found in a registry or an index; 54 of the 79 sources its claims cite were checked for support.

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

Open access

Published version of a preprint

Abstract

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.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

This study provides experimental evidence that APOBEC3A, APOBEC1 plus A1CF, and APOBEC3G can directly edit specific SARS-CoV-2 RNA sites to cause C-to-U mutations. Rather than inhibiting SARS-CoV-2, wild-type APOBEC3A expression promoted viral RNA replication and progeny production in Caco-2 cells, suggesting the virus can use APOBEC-mediated mutations for fitness and evolution.

Key findings

  • Among 307 total C sites in seven selected viral RNA segments, significant C-to-U editing sites numbered 135 for A1 + A1CF, 67 for A3A, and 11 for A3G; editing efficiencies reached as high as 22.2% for A1 + A1CF, 4.5% for A3A, and 0.18% for A3G.
  • A3A preferred a UC motif and A1 + A1CF preferred an AC motif; among C-to-U variations analyzed, AC-to-AU accounted for 38.23% and UC-to-UU for 31.83%, versus 15.44% for GC-to-GU and 14.50% for CC-to-CU.
  • In Caco-2 cells, A1 + A1CF or A3G expression did not significantly change viral RNA levels up to 96 h post-infection, but A3A expression significantly increased viral RNA at 72 and 96 h and produced approximately 10–100-fold higher viral progeny at 72 h.
  • Database analysis found C-to-U was the predominant mutation in the 5′UTR-Orf1a region, accounting for ~55.8% of single nucleotide variants; the A3A-edited UC241 mutation became a signature of dominant strains, with 90% of circulating variants carrying it by July 2020 and almost all by March 2021.

Methods

  • Cell-based RNA editing assay: seven 200 nt SARS-CoV-2 RNA segments were inserted into reporter vectors, co-transfected with APOBEC editor vectors (A1 + A1CF, A3A, or A3G) in HEK293T cells, and analyzed by Safe-Sequencing-System (SSS) deep sequencing with UID barcodes to distinguish true editing from PCR/sequencing errors.
  • Infection experiments: Caco-2 stable cell lines expressing A1 + A1CF, A3A, or A3G under a tetracycline-controlled promoter were infected with SARS-CoV-2 (MOI = 0.05); viral RNA replication was measured by qPCR for Nsp12, S, and N regions, and progeny yield by plaque assay on Vero E6 cells.
  • Sequence analysis: publicly available SARS-CoV-2 genome data from Nextstrain/GISAID (Dec. 2019 to Jan. 22, 2022) were analyzed for mutational events and C-to-U prevalence.

Limitations

  • The paper states that SSS sequencing is several times more expensive than regular NGS, each read covers ~150–200 nt, and it was too expensive to cover the entire ~30,000 nt SARS-CoV-2 genome, so only seven viral RNA segments were selected.
  • The paper states that A3G edited only a small number of sites (n = 11) and did not show a clear motif preference, so more SARS-CoV-2 regions and/or folded RNA structures are needed to verify A3G-mediated editing; UC/AC motifs alone are insufficient for efficient editing, and the exact RNA structural features remain to be defined.
  • The paper notes that a minor deaminase-independent enhancement by inactive A3A-E72A was observed and warrants further validation, and that host RNA or DNA editing by WT A3A could also contribute to increased viral replication/progeny; it also states that multiple viral passages in cell culture or animal models are needed to better assess mutation effects on fitness.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1038/s41598-022-19067-x

Cites

69 of 69 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Figure 1 APOBEC-mediated editing test of SARS-CoV-2 RNA. (A) Diagram of the SARS-CoV-2 genomic RNA, showing the positions (box) of the seven RNA segments (1–7) selected for studying the RNA editing by APOBECs. (B) Reporter vector (top) that contain each of the seven selected viral RNA segments that are transcribed into an RNA containing an AAV intron between the eGFP and the viral RNA segment. Splicing out the AAV intron yields a mature mRNA with a new spliced junction sequence (JUNC) that differs from its coding DNA, which can be used to selectively amplify either the mature mRNA or the coding DNA. (C) Three APOBEC editor vectors (top, A1-2A-A1CF, A3A, and A3G) and the Western blot showing their expression in 293 T cells (bottom). A1-2A-A1CF is constructed as one open reading frame (ORF) with a self-cleavage peptide T2A inserted between A1 and A1CF, which produced individual A1 and A1CF proteins in a 1:1 ratio28,58. (D) Strategy of the Safe-Sequencing-System (SSS) to minimize errors from PCR amplification and sequencing. After the SARS-CoV-2 RNAs from cell extracts are reverse transcribed, the cDNAs are sequentially amplified by the UID barcode (2 cycles) and the Illumina adapter (30 cycles). This SSS approach will distinguish the C-to-U mutations caused by APOBECs from the PCR and sequencing errors (see “Methods”).
Figure 1 APOBEC-mediated editing test of SARS-CoV-2 RNA. (A) Diagram of the SARS-CoV-2 genomic RNA, showing the positions (box) of the seven RNA segments (1–7) selected for studying the RNA editing by APOBECs. (B) Reporter vector (top) that contain each of the seven selected viral RNA segments that are transcribed into an RNA containing an AAV intron between the eGFP and the viral RNA segment. Splicing out the AAV intron yields a mature mRNA with a new spliced junction sequence (JUNC) that differs from its coding DNA, which can be used to selectively amplify either the mature mRNA or the coding DNA. (C) Three APOBEC editor vectors (top, A1-2A-A1CF, A3A, and A3G) and the Western blot showing their expression in 293 T cells (bottom). A1-2A-A1CF is constructed as one open reading frame (ORF) with a self-cleavage peptide T2A inserted between A1 and A1CF, which produced individual A1 and A1CF proteins in a 1:1 ratio28,58. (D) Strategy of the Safe-Sequencing-System (SSS) to minimize errors from PCR amplification and sequencing. After the SARS-CoV-2 RNAs from cell extracts are reverse transcribed, the cDNAs are sequentially amplified by the UID barcode (2 cycles) and the Illumina adapter (30 cycles). This SSS approach will distinguish the C-to-U mutations caused by APOBECs from the PCR and sequencing errors (see “Methods”).CC BY · PubMed Central
Figure 2 Local sequence context at the APOBEC-edited C sites on SARS-CoV-2 RNA. (A) Local sequences around the significantly edited target C sites (± 5 nucleotides from target C at position 0) by A1 + A1CF, A3A, or A3G. The editing level of each C site was normalized to the Ctrl, and only sites with 3 × or higher editing levels than the normalized value were defined as significant editing sites. (B) Analysis of local sequences around the top 30% edited C sites (or hotspot editing sites), showing predominantly AC motif for A1 + A1CF, UC for A3A, and CC for A3G. (C) Comparison of the C-to-U editing rates (%) of a particular dinucleotide motif by the three APOBECs. Each dot represents the C-to-U editing level obtained from the SSS results. In panel-D, statistical significance was calculated by unpaired two-tailed student’s t-test with P-values represented as: P > 0.05 = not significant; not indicated, *P < 0.05, ***P < 0.001, ****P < 0.0001.
Figure 2 Local sequence context at the APOBEC-edited C sites on SARS-CoV-2 RNA. (A) Local sequences around the significantly edited target C sites (± 5 nucleotides from target C at position 0) by A1 + A1CF, A3A, or A3G. The editing level of each C site was normalized to the Ctrl, and only sites with 3 × or higher editing levels than the normalized value were defined as significant editing sites. (B) Analysis of local sequences around the top 30% edited C sites (or hotspot editing sites), showing predominantly AC motif for A1 + A1CF, UC for A3A, and CC for A3G. (C) Comparison of the C-to-U editing rates (%) of a particular dinucleotide motif by the three APOBECs. Each dot represents the C-to-U editing level obtained from the SSS results. In panel-D, statistical significance was calculated by unpaired two-tailed student’s t-test with P-values represented as: P > 0.05 = not significant; not indicated, *P < 0.05, ***P < 0.001, ****P < 0.0001.CC BY · PubMed Central
Figure 3 Overall features of the RNA around the most preferred APOBEC-edited sites on SARS-CoV-2. The predicted RNA secondary structures62 of the sequences near the top 3 highest editing C sites by A1 + A1CF (A), A3A (B), and A3G (C) (see related Supplementary Dataset File 2). The editing efficiency of each site is listed at the top of each panel. In the secondary structure, the target C sites are highlighted in red, and -1 positions of the target C sites are highlighted in green for A, pink for U, and blue for C, respectively. In panel (A), the proposed canonical mooring sequences for A1 + A1CF (highlighted in sky blue) contain relatively high U/A/G contents downstream of the target C.
Figure 3 Overall features of the RNA around the most preferred APOBEC-edited sites on SARS-CoV-2. The predicted RNA secondary structures62 of the sequences near the top 3 highest editing C sites by A1 + A1CF (A), A3A (B), and A3G (C) (see related Supplementary Dataset File 2). The editing efficiency of each site is listed at the top of each panel. In the secondary structure, the target C sites are highlighted in red, and -1 positions of the target C sites are highlighted in green for A, pink for U, and blue for C, respectively. In panel (A), the proposed canonical mooring sequences for A1 + A1CF (highlighted in sky blue) contain relatively high U/A/G contents downstream of the target C.CC BY · PubMed Central
Figure 4 The potential effect of APOBEC-mediated editing on SARS-CoV-2 mutations and fitness. (A) The number of mutational events (all single nucleotide variants) on SARS-CoV-2 RNA segment 5’UTR-Orf1a (segment 1 in Fig. 1A) from the SARS-CoV-2 genome sequence data (the Nextstrain datasets from Dec. 2019 to Jan. 22nd, 2022 downloaded from the GISAID database, https://www.gisaid.org/hcov19-variants/ and https://nextstrain.org/ncov/global). The C203, C222, and C241 represent many of the C-to-U mutational events (asterisks) with the A3A-editing UC motif in the SARS-CoV-2 variants. (B) The A3A-mediated C-to-U editing rate on UC motif in the same 5′UTR-Orf1a region obtained from our cell-based editing system and the SSS analysis. The Ctrl (EV) editing levels (or background error rates) of the corresponding region are presented as negative values (%). The C203, C222, and C241 (asterisks) all showed significant editing by A3A. (C) The C-to-U mutation prevalence over time at C203, C222, and C241. The sequencing frequency is represented by C in blue and U in yellow (referred to the Nextstrain datasets: https://nextstrain.org/ncov/global). This analysis showed that SARS-CoV-2 started to acquire the C-to-U mutation at C241 in January 2020. By July 2020, 90% of the circulating viral variants carry this mutation at C241. By March 2021, almost all circulating viral variants have this mutation, suggesting the C241 to U mutation in the 5’UTR is beneficial to the viral fitness (see Supplementary Fig. S14).
Figure 4 The potential effect of APOBEC-mediated editing on SARS-CoV-2 mutations and fitness. (A) The number of mutational events (all single nucleotide variants) on SARS-CoV-2 RNA segment 5’UTR-Orf1a (segment 1 in Fig. 1A) from the SARS-CoV-2 genome sequence data (the Nextstrain datasets from Dec. 2019 to Jan. 22nd, 2022 downloaded from the GISAID database, https://www.gisaid.org/hcov19-variants/ and https://nextstrain.org/ncov/global). The C203, C222, and C241 represent many of the C-to-U mutational events (asterisks) with the A3A-editing UC motif in the SARS-CoV-2 variants. (B) The A3A-mediated C-to-U editing rate on UC motif in the same 5′UTR-Orf1a region obtained from our cell-based editing system and the SSS analysis. The Ctrl (EV) editing levels (or background error rates) of the corresponding region are presented as negative values (%). The C203, C222, and C241 (asterisks) all showed significant editing by A3A. (C) The C-to-U mutation prevalence over time at C203, C222, and C241. The sequencing frequency is represented by C in blue and U in yellow (referred to the Nextstrain datasets: https://nextstrain.org/ncov/global). This analysis showed that SARS-CoV-2 started to acquire the C-to-U mutation at C241 in January 2020. By July 2020, 90% of the circulating viral variants carry this mutation at C241. By March 2021, almost all circulating viral variants have this mutation, suggesting the C241 to U mutation in the 5’UTR is beneficial to the viral fitness (see Supplementary Fig. S14).CC BY · PubMed Central
Figure 5 SARS-CoV-2 replication and virion production in cells expressing APOBECs. (A) Overview of experiments for SARS-CoV-2 replication and viral production in the presence of APOBECs. The Caco-2 stable cell lines were constructed to express A1 + A1CF, A3A, or A3G under a tetracycline-controlled promoter. The Caco-2-APOBEC stable cell lines were then infected with SARS-CoV-2 (MOI = 0.05), and the viral RNA replication and progeny production were measured at different time points. (B) Effect of each APOBEC expression on SARS-CoV-2 viral RNA replication. Measurement of relative viral RNA abundance at different time points after viral infection of the Caco-2-APOBEC stable cell lines expressing A1 + A1CF, A3A, or A3G. The viral RNA abundance was measured using real-time quantitative PCR (qPCR) to detect RNA levels by using specific primers to amplify three separate viral regions, the Nsp12, S, or N coding regions (see Methods in SI). (C) Effect of each APOBEC expression on SARS-CoV-2 progeny production. Infectious viral progeny yield harvested in the medium at 48 h and 72 h post-infection was determined by plaque assay (see “Methods”). In panel (B) and (C), statistical significance was calculated by unpaired two-tailed student’s t-test with P-values represented as: P > 0.05 = not significant, ***P < 0.001.
Figure 5 SARS-CoV-2 replication and virion production in cells expressing APOBECs. (A) Overview of experiments for SARS-CoV-2 replication and viral production in the presence of APOBECs. The Caco-2 stable cell lines were constructed to express A1 + A1CF, A3A, or A3G under a tetracycline-controlled promoter. The Caco-2-APOBEC stable cell lines were then infected with SARS-CoV-2 (MOI = 0.05), and the viral RNA replication and progeny production were measured at different time points. (B) Effect of each APOBEC expression on SARS-CoV-2 viral RNA replication. Measurement of relative viral RNA abundance at different time points after viral infection of the Caco-2-APOBEC stable cell lines expressing A1 + A1CF, A3A, or A3G. The viral RNA abundance was measured using real-time quantitative PCR (qPCR) to detect RNA levels by using specific primers to amplify three separate viral regions, the Nsp12, S, or N coding regions (see Methods in SI). (C) Effect of each APOBEC expression on SARS-CoV-2 progeny production. Infectious viral progeny yield harvested in the medium at 48 h and 72 h post-infection was determined by plaque assay (see “Methods”). In panel (B) and (C), statistical significance was calculated by unpaired two-tailed student’s t-test with P-values represented as: P > 0.05 = not significant, ***P < 0.001.CC BY · PubMed Central

Verification

Claims the paper makes

55 Own result · 28 Own interpretation · 14 Partly supported · 11 Not checked yet · 6 Own method · 5 Supported · 1 Couldn’t confirm · 1 No full text of the sources · accuracy not yet measured

Verify all 21 references — free account
  • Probably not a substantial claim. 20
  • It may state two separate claims as one. 35
  • The quote only partly states this claim. 1
51 of 51
  • The paper relies on this without support.

    There is no report on direct experimental evidence to demonstrate if APOBECs can edit the SARS-CoV-2 genome.

    “However, there is no report on the direct experimental evidence to demonstrate if APOBECs can edit the SARS-CoV-2 genome and, if so, what is the extent of editing on the viral genome by different APOBECs, and the potential effect of the APOBEC-editing on the virus.”
  • The paper relies on this without support.

    A1 + A1CF is expressed in the small intestine and liver during SARS-CoV-2 infection, enabling AC-to-AU editing in patients.

    “These results indicate that many of the AC-to-AU mutations in the SARS-CoV-2 genome from patients can be driven by A1 + A1CF mediated RNA-editing in the small intestine and liver with SARS-CoV-2 infection.”
  • The paper relies on this without support.

    The inactive-A3A mutant used lacks deaminase activity.

    “even though some increase of replication and progeny production was also observed in the inactive-A3A mutant expressing cell line (Supplementary Fig. S10).”
  • The paper relies on this without support.

    The HEK293T cell-based editing assay reflects APOBEC-mediated editing in natural SARS-CoV-2 infection.

    “Here we described the first experimental evidence demonstrating APOBEC enzymes, A1 + A1CF and A3A can target specific SARS-CoV-2 viral sequences for RNA editing, and the resulting mutations likely contribute to the viral replication and fitness.”
4 of 4

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

2 of 2

Citations

Cites

Show all 71 references

Cited by

Sorted by citations

See all 117

Metadata

Authors
Kim, K., Calabrese, P., Wang, S., Qin, C., Rao, Y., Feng, P., & Chen, X. S.
Affiliations
University of Southern California (US)
Corresponding author
Chen, Xiaojiang S.
Venue
Scientific Reports
Published
2022-09-13
Volume
12
Issue
1
Pages
14972-14972
ISSN
2045-2322 (ISSN-L)
OpenAlex
W4295346375
Cited by
115
References
71
Citations by year
2026: 13 · 2025: 24 · 2024: 41 · 2023: 27 · 2022: 10
FWCI
7.65
Type
journal article
Language
en
Open access
Locations
  • Scientific Reports (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc-nd)
  • Research Square (acceptedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • PubMed Central (submittedVersion, cc-by-nc-nd)
  • PubMed Central (submittedVersion, cc-by)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
Funders
National Institutes of Health (AI150524)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA regulation and disease
  • Health Sciences › Medicine › Infectious Diseases › SARS-CoV-2 and COVID-19 Research
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Evolution and Genetic Dynamics
Keywords
SARS-CoV-2, APOBEC3A, APOBEC1, APOBEC3G, viral RNA replication, viral fitness, viral evolution, host-virus interaction, Caco-2 cells, viral genome structure

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-11-10, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-11-10 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review