CoChat
Sign Up
Cited by 192interferon and immune responsesEvolution and Genetic DynamicsRNA and protein synthesis mechanisms

The DOI matches its Crossref record; 57 of the 57 references it cites were found in a registry or an index.

The influence of CpG and UpA dinucleotide frequencies on RNA virus replication and characterization of the innate cellular pathways underlying virus attenuation and enhanced replication

Open access
Atkinson, Nicky J., Witteveldt, Jeroen, Evans, David J., Simmonds, Peter·2014·Nucleic Acids Research
DOI 10.1093/nar/gku075PMC PMC3985648
Open this paper in CoChat — free accountOpen PDF

Abstract

Most RNA viruses infecting mammals and other vertebrates show profound suppression of CpG and UpA dinucleotide frequencies. To investigate this functionally, mutants of the picornavirus, echovirus 7 (E7), were constructed with altered CpG and UpA compositions in two 1.1–1.3 Kbase regions. Those with increased frequencies of CpG and UpA showed impaired replication kinetics and higher RNA/infectivity ratios compared with wild-type virus. Remarkably, mutants with CpGs and UpAs removed showed enhanced replication, larger plaques and rapidly outcompeted wild-type virus on co-infections. Luciferase-expressing E7 sub-genomic replicons with CpGs and UpAs removed from the reporter gene showed 100-fold greater luminescence. E7 and mutants were equivalently sensitive to exogenously added interferon-β, showed no evidence for differential recognition by ADAR1 or pattern recognition receptors RIG-I, MDA5 or PKR. However, kinase inhibitors roscovitine and C16 partially or entirely reversed the attenuated phenotype of high CpG and UpA mutants, potentially through inhibition of currently uncharacterized pattern recognition receptors that respond to RNA composition. Generating viruses with enhanced replication kinetics has applications in vaccine production and reporter gene construction. More fundamentally, the findings introduce a new evolutionary paradigm where dinucleotide composition of viral genomes is subjected to selection pressures independently of coding capacity and profoundly influences host–pathogen interactions.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1093/nar/gku075

Cites

57 of 57 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1. Genome organization of E7 and positions of mutated insert regions. Insert positions are compared with genome diagram and a plot of sequence variability within species B enteroviruses at synonymous sites (blue line) and folding energies indicative of RNA secondary structure (red and pink lines). Variability at synonymous sites (left y-axis) was computed at each codon position in alignments, plotted with a window size of 41 codons. MFED values (right y-axis) for sense and antisense RNA sequences were calculated for 200 base fragments, incrementing by 48 bases; values plotted represent mean values of five consecutive fragments. Nucleotide positions were calculated relative to the pT7:E7 cDNA sequence.
Figure 1. Genome organization of E7 and positions of mutated insert regions. Insert positions are compared with genome diagram and a plot of sequence variability within species B enteroviruses at synonymous sites (blue line) and folding energies indicative of RNA secondary structure (red and pink lines). Variability at synonymous sites (left y-axis) was computed at each codon position in alignments, plotted with a window size of 41 codons. MFED values (right y-axis) for sense and antisense RNA sequences were calculated for 200 base fragments, incrementing by 48 bases; values plotted represent mean values of five consecutive fragments. Nucleotide positions were calculated relative to the pT7:E7 cDNA sequence.CC BY · PubMed Central
Figure 2. RNA to infectivity ratios of WT and viruses with modified CpG/UpA frequencies. WT and mutant viruses were recovered from RD cells and TCID50 titres determined (note log scale). The number of viral genome copies was determined through qRT-PCR and compared with the infectivity titre. Results are the mean and standard error from three separate extractions. RNA/infectivity ratios were additionally calculated for C|C and U|U mutants in the presence of C16.
Figure 2. RNA to infectivity ratios of WT and viruses with modified CpG/UpA frequencies. WT and mutant viruses were recovered from RD cells and TCID50 titres determined (note log scale). The number of viral genome copies was determined through qRT-PCR and compared with the infectivity titre. Results are the mean and standard error from three separate extractions. RNA/infectivity ratios were additionally calculated for C|C and U|U mutants in the presence of C16.CC BY · PubMed Central
Figure 3. Replication kinetics and titres at 24 h of WT and modified viruses infected at a low MOI. RD cells were infected with E7 WT, permuted (P|P), CpG- and UpA-high mutants (C|C and U|U; A and C) or CpG- and UpA-low mutants (c|c, u|u, cu|cu; B and D) viruses at an MOI of 0.01. Infectious titre of supernatants was quantified at indicated time points by TCID50 (A and B) and mean titres and SEMs at 24 h. Results are the mean of three biological replicates.
Figure 3. Replication kinetics and titres at 24 h of WT and modified viruses infected at a low MOI. RD cells were infected with E7 WT, permuted (P|P), CpG- and UpA-high mutants (C|C and U|U; A and C) or CpG- and UpA-low mutants (c|c, u|u, cu|cu; B and D) viruses at an MOI of 0.01. Infectious titre of supernatants was quantified at indicated time points by TCID50 (A and B) and mean titres and SEMs at 24 h. Results are the mean of three biological replicates.CC BY · PubMed Central
Figure 4. Plaque morphology of E7 WT and double region mutant viruses. RD cell monolayers in 10-cm plates were infected with a similar infectious titre of virus and incubated for 96 h at 37°C. (A) Plaque appearance. (B) Plaques sizes of WT and mutant E7 viruses calculated from 25 plaques for each virus using ImageJ software (mean values and SEMs relative to WT control shown as bar heights and error bars).
Figure 4. Plaque morphology of E7 WT and double region mutant viruses. RD cell monolayers in 10-cm plates were infected with a similar infectious titre of virus and incubated for 96 h at 37°C. (A) Plaque appearance. (B) Plaques sizes of WT and mutant E7 viruses calculated from 25 plaques for each virus using ImageJ software (mean values and SEMs relative to WT control shown as bar heights and error bars).CC BY · PubMed Central
Figure 5. Synchronized infection with equal viral genome copies. RD cells were synchronously infected with 1000 genome copies of WT, R1/R2 CpG-high or R1/R2 UpA-high virus, as calculated using qRT-PCR. Cells were trypsinized and washed 1 or 4 h post infection and the intracellular viral load determined by qRT-PCR. Results are the mean and standard error of three biological replicates.
Figure 5. Synchronized infection with equal viral genome copies. RD cells were synchronously infected with 1000 genome copies of WT, R1/R2 CpG-high or R1/R2 UpA-high virus, as calculated using qRT-PCR. Cells were trypsinized and washed 1 or 4 h post infection and the intracellular viral load determined by qRT-PCR. Results are the mean and standard error of three biological replicates.CC BY · PubMed Central
Figure 6. Analysis of luciferase expression driven by E7 replicons with reduced CpG/UpA frequencies. Replicons were generated with reduced CpG/UpA frequencies, based on the backbone pRiboE7luc replicon, in which the structural genes of E7 are replaced by an insect luciferase gene. In the cu|-|W replicon the luciferase gene itself was modified to minimize both CpG and UpA frequency; in the cu|-|c and cu|-|u replicons Region 2 was additionally modified to further reduce either CpG or UpA frequency. RNA was generated from replicons and transfected into RD cells. Luminescence was measured relative to the mock-transfected control. Results are the mean and standard error of three biological replicates.
Figure 6. Analysis of luciferase expression driven by E7 replicons with reduced CpG/UpA frequencies. Replicons were generated with reduced CpG/UpA frequencies, based on the backbone pRiboE7luc replicon, in which the structural genes of E7 are replaced by an insect luciferase gene. In the cu|-|W replicon the luciferase gene itself was modified to minimize both CpG and UpA frequency; in the cu|-|c and cu|-|u replicons Region 2 was additionally modified to further reduce either CpG or UpA frequency. RNA was generated from replicons and transfected into RD cells. Luminescence was measured relative to the mock-transfected control. Results are the mean and standard error of three biological replicates.CC BY · PubMed Central
Figure 7. Fitness determination by competition assays between WT and modified viruses. RD cells were infected with an equal MOI of WT and modified virus, and the supernatant serially passaged through cells. RNA was isolated and the composition of each virus determined through selective restriction digests. Images show the virus composition in the starting inoculum and in three biological replicates following passage of CpG- and UpA-low mutants and of U|W.
Figure 7. Fitness determination by competition assays between WT and modified viruses. RD cells were infected with an equal MOI of WT and modified virus, and the supernatant serially passaged through cells. RNA was isolated and the composition of each virus determined through selective restriction digests. Images show the virus composition in the starting inoculum and in three biological replicates following passage of CpG- and UpA-low mutants and of U|W.CC BY · PubMed Central
Figure 8. Pairwise fitness comparison between E7 WT and mutant with varying degrees of CpG and/or UpA under-representation. RD cells were infected with an equal MOI of two viruses represented in columns and rows of the matrix and the supernatant serially passaged. The composition of each virus was determined through restriction endonuclease cleavage (see Figure 7) and outcome displayed by colour shading. The key refers to population representations of viruses listed in columns (for example, all but one variants were out-competed by the cu|cu mutant). The fitness ranking deduced from these results is shown underneath the figure.
Figure 8. Pairwise fitness comparison between E7 WT and mutant with varying degrees of CpG and/or UpA under-representation. RD cells were infected with an equal MOI of two viruses represented in columns and rows of the matrix and the supernatant serially passaged. The composition of each virus was determined through restriction endonuclease cleavage (see Figure 7) and outcome displayed by colour shading. The key refers to population representations of viruses listed in columns (for example, all but one variants were out-competed by the cu|cu mutant). The fitness ranking deduced from these results is shown underneath the figure.CC BY · PubMed Central
Figure 9. Effect of exogenous IFN-β on viral replication. RD cells were treated with 5, 50 and 500 U/ml human IFN-β, or mock treated, then infected with WT E7, C|C, U|U and cu|cu mutants at an MOI of 1. RNA was isolated after 8 h and viral load determined by qRT-PCR. Results are the mean and standard error of two biological replicates.
Figure 9. Effect of exogenous IFN-β on viral replication. RD cells were treated with 5, 50 and 500 U/ml human IFN-β, or mock treated, then infected with WT E7, C|C, U|U and cu|cu mutants at an MOI of 1. RNA was isolated after 8 h and viral load determined by qRT-PCR. Results are the mean and standard error of two biological replicates.CC BY · PubMed Central
Figure 10. Analysis of gene expression in cells infected with WT E7 and mutants. RD cells were infected with E7 at an MOI of 10 and changes in expression at 8 h of a range of innate and adaptive immune response genes determined using a PCR Array. Responses were compared with those induced in cells infected by high CpG- and UpA mutants at equal copy number to WT or equal MOI. Parallel infections with Sendai virus (SV) and transfection of poly-I:C were performed. Non-specific induction of genes was controlled for by parallel assays of cells exposed to u/v inactivated supernatant corresponding in amount to that added in the equal copy number (C1) and equal MOI (C2) assays.
Figure 10. Analysis of gene expression in cells infected with WT E7 and mutants. RD cells were infected with E7 at an MOI of 10 and changes in expression at 8 h of a range of innate and adaptive immune response genes determined using a PCR Array. Responses were compared with those induced in cells infected by high CpG- and UpA mutants at equal copy number to WT or equal MOI. Parallel infections with Sendai virus (SV) and transfection of poly-I:C were performed. Non-specific induction of genes was controlled for by parallel assays of cells exposed to u/v inactivated supernatant corresponding in amount to that added in the equal copy number (C1) and equal MOI (C2) assays.CC BY · PubMed Central
Figure 11. Replication of WT and dinucleotide-modified viruses in an IRF3-blocked cell line. A549 cells expressing the viral protein Npro and non-expressing control cells were infected with WT, C|C (R1/R2 CpG-high) and U|U (R1/R2 UpA-high) at an MOI of 0.01. RNA was extracted from cell culture supernatant at 24 h and viral load determined by qRT-PCR. Results are the mean and standard error of three biological replicates.
Figure 11. Replication of WT and dinucleotide-modified viruses in an IRF3-blocked cell line. A549 cells expressing the viral protein Npro and non-expressing control cells were infected with WT, C|C (R1/R2 CpG-high) and U|U (R1/R2 UpA-high) at an MOI of 0.01. RNA was extracted from cell culture supernatant at 24 h and viral load determined by qRT-PCR. Results are the mean and standard error of three biological replicates.CC BY · PubMed Central
Figure 12. Effect of inhibiting PKR, RIG-I and MDA-5 expression on E7 replication. (A) Comparison of virus replication in control A549 cells at 8 h compared with that in cells expressing shRNAs targeting PKR, RIG-I and MDA-5 (reduction in mRNA levels achieved is shown in Supplementary Figure S3). (B) RD cells were pre-treated for 24 h with an siRNA directed against PKR and infected with WT virus or CpG-high, UpA-high or CPg/UpA-low mutants at an MOI of 0.03. RNA was isolated after 22 h and intracellular viral loads determined by qRT-PCR. Results are shown relative to the replication rate in cells treated with a control (validated non-target) siRNA.
Figure 12. Effect of inhibiting PKR, RIG-I and MDA-5 expression on E7 replication. (A) Comparison of virus replication in control A549 cells at 8 h compared with that in cells expressing shRNAs targeting PKR, RIG-I and MDA-5 (reduction in mRNA levels achieved is shown in Supplementary Figure S3). (B) RD cells were pre-treated for 24 h with an siRNA directed against PKR and infected with WT virus or CpG-high, UpA-high or CPg/UpA-low mutants at an MOI of 0.03. RNA was isolated after 22 h and intracellular viral loads determined by qRT-PCR. Results are shown relative to the replication rate in cells treated with a control (validated non-target) siRNA.CC BY · PubMed Central
Figure 13. Effect of the kinase inhibitors, C16, 2-AP and Roscovitine on E7 replication. Cells were pre-treated for 45 min with 2 μM C16, 5 mM 2-AP or 40 μM Roscovitine and infected with WT virus or CpG-high, UpA-high or CpG/UpA-low mutants at an MOI of 0.1. RNA was isolated at 24 h, and intracellular viral loads were determined by qRT-PCR. Results are shown relative to the replication rate in cells untreated cells (note log scale).
Figure 13. Effect of the kinase inhibitors, C16, 2-AP and Roscovitine on E7 replication. Cells were pre-treated for 45 min with 2 μM C16, 5 mM 2-AP or 40 μM Roscovitine and infected with WT virus or CpG-high, UpA-high or CpG/UpA-low mutants at an MOI of 0.1. RNA was isolated at 24 h, and intracellular viral loads were determined by qRT-PCR. Results are shown relative to the replication rate in cells untreated cells (note log scale).CC BY · PubMed Central
Figure 14. Influence of the PKR inhibitor E7 cytopathology and infectivity. (A) Cytopathology of WT, C|C and U|U mutants in RD cells infected at MOIs of 0.1 and 0.01 in 2 μM C16 or untreated control. Monolayers of untreated cells remain intact on infection with C|C and U|U mutants while those treated with C16 display a complete CPE even at low MOI. (B) Infectivity determination of pre-titred virus stock supernatant (3000 PFU/ml) of WT virus or CpG-high, UpA-high or CPg/UpA-low mutants in the presence of C16. The enhanced infectivity of C|C and U|U mutants restored the RNA/infectivity ratios to those of WT virus (Figure 2).
Figure 14. Influence of the PKR inhibitor E7 cytopathology and infectivity. (A) Cytopathology of WT, C|C and U|U mutants in RD cells infected at MOIs of 0.1 and 0.01 in 2 μM C16 or untreated control. Monolayers of untreated cells remain intact on infection with C|C and U|U mutants while those treated with C16 display a complete CPE even at low MOI. (B) Infectivity determination of pre-titred virus stock supernatant (3000 PFU/ml) of WT virus or CpG-high, UpA-high or CPg/UpA-low mutants in the presence of C16. The enhanced infectivity of C|C and U|U mutants restored the RNA/infectivity ratios to those of WT virus (Figure 2).CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

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

The self-check needs the paper’s full text.

Citations

Cites

Show all 57 references

Cited by

Sorted by citations

See all 192

Related

Metadata

Authors
Atkinson, N. J., Witteveldt, J., Evans, D. J., & Simmonds, P.
Affiliations
  • Roslin Institute (GB)
  • 1 Infection and Immunity Division, Roslin Institute, University of Edinburgh, Easter Bush, Edinburgh EH25 9RG, UK and 2School of Life Sciences, University of Warwick, Coventry CV4 7AL, UK
ORCID
Venue
Nucleic Acids Research
Publisher
Oxford University Press (OUP)
Published
2014-01-25
Volume
42
Issue
7
Pages
4527-4545
ISSN
  • 0305-1048 (ISSN-L)
  • 1362-4962
DOI
10.1093/nar/gku075
PMC
PMC3985648
OpenAlex
W2057088297
Cited by
192
References
57
Citations by year
2026: 12 · 2025: 15 · 2024: 9 · 2023: 8 · 2022: 11 · 2021: 36 · 2020: 26 · 2019: 17 · 2018: 11 · 2017: 20 · 2016: 10 · 2015: 13 · 2014: 3
FWCI
7.01
Type
journal article
Language
en
Open access
Licence
Locations
  • Nucleic Acids Research (publishedVersion, cc-by)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Warwick Research Archive Portal (University of Warwick) (acceptedVersion)
  • Edinburgh Research Explorer (University of Edinburgh) (submittedVersion, cc-by)
  • St Andrews Research Repository (St Andrews Research Repository) (submittedVersion, cc-by)
Funders
  • Wellcome Trust (095831, WT087628MA)
  • University of St Andrews
  • Medical Research Council
  • Biotechnology and Biological Sciences Research Council (BBS/E/D/20241864)
Topics
  • Life Sciences › Immunology and Microbiology › Immunology › interferon and immune responses
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Evolution and Genetic Dynamics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA and protein synthesis mechanisms
Keywords
CpG, RNA virus replication, echovirus 7, picornavirus, viral genome evolution, RNA composition, innate immune recognition, pattern recognition receptors, interferon-beta, roscovitine, RNA replicon, vaccine production, host-pathogen interaction
MeSH terms
Base Composition, Cell Line, Mutation, RNA, Viral, Virus Replication, Dinucleoside Phosphates, GC Rich Sequence, Enterovirus B, Human

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.

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