Vaccinia recombination 3 - Journal of Virology

JVI Accepts, published online ahead of print on 26 February 2014
J. Virol. doi:10.1128/JVI.00022-14
Copyright © 2014, American Society for Microbiology. All Rights Reserved.
1
Title: Genome scale patterns of recombination between co-infecting vaccinia viruses
2
3
Running Title: Vaccinia recombination
4
5
Authors: Li Qin1, and David H. Evans1#
6
7
8
9
Affiliation: 1Dept. of Medical Microbiology & Immunology and
10
Li Ka Shing Institute of Virology
11
6020H Katz Group Centre
12
University of Alberta
13
Edmonton, AB
14
T6G 2H7 CANADA
15
16
17
Word Counts:
18
Abstract: 250 words
19
Importance: 150 words
20
Text (including above): 6,106 words
21
22
#Corresponding author: Dept. of Medical Microbiology & Immunology, 6020 Katz
23
Group Centre, University of Alberta, Edmonton, AB, T6G 2H7, CANADA. Phone (780)
24
492-2308; FAX (780) 492-7521; E-mail [email protected]
25
ABSTRACT
26
27
Recombination plays a critical role in virus evolution. It helps avoid genetic decline and
28
creates novel phenotypes. This promotes survival, and genome sequencing suggests that
29
recombination
30
Orthopoxviruses like variola virus. Recombination can also be used to map genes, but
31
although recombinant poxviruses are easily produced in culture, classical attempts to map
32
the vaccinia virus (VACV) genome this way, met with little success. We have sequenced
33
recombinants formed when VACV TianTan and Dryvax strains are crossed under
34
different conditions. These were a single round of growth in co-infected cells, five rounds
35
of sequential passage, or using Leporipoxvirus-mediated DNA reactivation. Our studies
36
showed that recombinants contain a patchwork of DNA, with the number of exchanges
37
increasing with passage. Further passage also selected for TianTan DNA and correlated
38
with increased plaque size. The recombinants produced through a single round of
39
co-infection contain a disproportionate number of short conversion tracks (<1 kbp) and
40
exhibited 1 exchange per 12 kbp, close to the ~1 per 8 kbp in the literature. One
41
byproduct of this study was that rare mutations were also detected, VACV replication
42
produces ~1×10-8 mutations per nucleotide copied per cycle of replication and ~1 large
43
(21 kbp) deletion per 70 rounds of passage. Viruses produced using DNA reactivation
44
appeared no different from recombinants produced using ordinary methods. An attractive
45
feature of this approach is that, when combined with selection for a particular phenotype,
46
it provides a way of mapping and dissecting more complex virus traits.
has
facilitated
the
evolution
of
human
pathogens
including
47
48
Qin et al.
Vaccinia virus recombination
Page 2 of 31
49
IMPORTANCE
50
51
When two closely related viruses co-infect the same cell, they can swap genetic
52
information through a process called recombination. Recombination produces new
53
viruses bearing different combinations of genes, and it plays an important role in virus
54
evolution. Poxviruses are a family of viruses that includes variola (or smallpox) virus and
55
although poxviruses are known to recombine, no one has previously mapped the patterns
56
of DNAs exchanged between viruses. We co-infected cells with two different vaccinia
57
poxviruses, isolated the progeny, and sequenced them.
58
recombination is a very accurate process that assembles viruses containing DNA copied
59
from both parents. In a single round of infection, DNA is swapped back and forth ~18
60
times per genome to make recombinant viruses that are a mosaic of the two parental
61
DNAs. This mixes many different genes in complex combinations and illustrates how
62
recombination can produce viruses with greatly altered disease potential.
We show that poxvirus
63
Qin et al.
Vaccinia virus recombination
Page 3 of 31
64
INTRODUCTION
65
66
Recombination plays an essential role in DNA repair and, by creating new combinations
67
of genetic traits, it averts the decline in fitness caused by the accumulation of mutations
68
(”Muller’s ratchet”) while creating the genetic diversity that is the substrate for
69
Darwinian selection. Bacteriophage were the first viruses shown to be subject to
70
recombination (1) and the phenomenon was soon also detected in co-cultures of many
71
other viruses including herpes simplex virus (2) and vaccinia virus (3). In the years
72
immediately following, research showed that in vitro recombination could also produce
73
hybrids between related poxviruses such as rabbit fibroma and myxoma viruses (4) and
74
between variola virus and cowpox and rabbitpox viruses (5, 6). The subsequent discovery
75
of another natural hybrid between rabbit fibroma and myxoma viruses, called malignant
76
rabbit virus (7), suggested that poxviruses can also recombine in co-infected animals and
77
the significance for human health is illustrated by the fact that variola minor virus may be
78
a recombinant derived from more virulent West African and Asian variola major strains
79
(8). We have recently published an analysis of some of the strain variants found in an old
80
non-clonal smallpox vaccine, Dryvax, and detected one virus bearing a small region of
81
sequence wherein the single nucleotide polymorphisms (SNPs) were more characteristic
82
of horsepox than of vaccinia virus (VACV) (9). This may represent a sequence relic that
83
has been retained in the absence of counter selection and the population bottlenecks
84
caused by periodic plaque purification.
85
Genetic crosses between viruses encoding different selectable markers were once
86
also used to try and assemble recombination-based maps of VACV (10-12), although the
87
method never proved very useful and was soon supplanted by marker rescue and DNA
88
sequencing technologies. This was due to the limited distances over which linkage is
89
retained relative to the spacing between many markers (<20 kbp over a 200 kbp genome)
Qin et al.
Vaccinia virus recombination
Page 4 of 31
90
and the difficulty of reproducibly measuring recombinant frequencies (13). Further
91
studies showed that homologous recombination can be used to genetically modify VACV
92
(14, 15), that this is an accurate process (16, 17), and that these processes also operate in
93
trans and can be detected using transfected DNAs (18). Poxviruses replicate in
94
sequestered structures called factories (19), each of which derives from a single infecting
95
particle, and it is presumed that recombinants can only form within these factories if
96
different DNAs mix in the presence of the recombination machinery. VACV uses a
97
single-strand annealing mechanism to produce recombinant molecules in a reaction
98
catalyzed by the E9 viral DNA polymerase and I3 single-strand DNA binding protein (20,
99
21) and, since both proteins primarily reside within virus factories, that is presumably
100
where recombination also occurs. We have suggested that random variations in the timing
101
and degree of mixing of virus factories within co-infected cells could explain why
102
recombinant frequencies proved difficult to measure reproducibly (13). If two or more
103
viruses infect any particular cell, but a portion of the factories don’t mix, such a process
104
would decrease the yield of recombinants in a stochastic manner relative to the number of
105
non-recombinant (i.e. fully parental) viruses produced by DNA replication.
106
Although much has been learned concerning the mechanism of poxvirus
107
recombination, questions remain regarding how these processes and physical constraints
108
might affect the patterns of DNA exchange, and thus the overall genetic composition of
109
the resulting pool of parental and recombinant viruses. How genetic linkage distances
110
relate to the actual numbers of exchanges in recombinant viruses also remains to be
111
established. In this study we have used the ~1400 SNPs that differentiate two strains of
112
VACV, a Dryvax clone and a TianTan clone, as sequence tags that can be used to track
113
the origins of the different DNA segments in recombinant progeny. Our study shows that
114
VACV recombination reactions produce genomes exhibiting a “patchwork” of exchanges,
115
some apparently derived from a succession of crossovers over the course of a single
Qin et al.
Vaccinia virus recombination
Page 5 of 31
116
infection cycle. Interestingly, viruses produced using non-genetic reactivation methods
117
(22), appear indistinguishable from recombinants produced in a more regular manner,
118
showing that such viruses are likely subjected to similar replication and developmental
119
pathways.
120
121
METHODS
122
123
Cells and viruses. Vaccinia virus strains DPP17 (GenBank JN654983), TP03 (GenBank
124
KC207810), and TP05 (Genbank KC207811) were cloned from stocks of Dryvax
125
(DPP17) and TianTan (TP03/05) viruses (9, 23). They were cultured on BSC-40 cells in
126
modified Eagle’s medium (MEM, Gibco) supplemented with 5% fetal bovine serum, 1%
127
nonessential amino acids, 1% L-glutamine, and 1% antibiotic at 37°Cin a 5%CO2
128
atmosphere. Two types of recombinant virus stocks were prepared. The DTM viruses
129
(Dryvax-TianTan mixture) were generated by co-infecting cells with DPP17 and TP05 at
130
a multiplicity of infection (MOI) of 0.02 (each 0.01), culturing the cells for two days,
131
harvesting the cell-virus mixture, and releasing the virus by freeze-thaw. A sample (10
132
μL, 0.5% of the lysate, or ~0.02 PFU/cell) was then used to infect another fresh dish of
133
cells and this was repeated for a total of five rounds of passage. Individual DTM viruses
134
were then isolated using three rounds of cloning by limiting dilution as described (9). The
135
DTH viruses (Dryvax-TianTan high MOI) were produced by co-infecting cells with
136
DPP17 and TP05 at MOI=10 (each 5) for 24 hr, followed by three rounds of cloning.
137
Plaque images were processed with ImageJ (24).
138
139
Virus DNA reactivation. A third collection of recombinant VACV were prepared using
140
DNA reactivation reactions as described previously (22). Briefly, BGMK cells were
141
grown to near confluence in 60 mm dishes and infected with Shope fibroma virus (SFV,
Qin et al.
Vaccinia virus recombination
Page 6 of 31
142
strain Kasza) at MOI=1 in PBS. After one hour at 37°C, the medium was replaced with
143
MEM containing 10% fetal bovine serum, incubated for another hour, and the MEM
144
replaced with Opti-MEM (Gibco). DPP17 and TP03 VACV DNAs were extracted from
145
sucrose gradient-purified virions using phenol chloroform, mixed in 1:1 ratio, and the
146
SFV-infected cells transfected with 5 μg of this DNA using Lipofectamine 2000 (Gibco).
147
The cells were incubated for 4 hr, the Opti-MEM media was replaced with MEM
148
containing 10% serum, and the cells cultured for another 3 days. The cells were then
149
subjected to three rounds of freeze-thaw and 0.2 mL used to infect BSC-40 cells (which
150
do not support SFV growth). The reactivated VACV were cloned using three rounds of
151
limited dilution and recombinants were identified using the PCR and the primer pairs
152
DVX-209F and DVX-226R plus DVX-004F and DVX-007R (Table 1). The 209F/226R
153
primers should produce a 1.1 kbp amplicon in reactions containing DPP17 DNA whereas
154
TP03 DNA does not serve as a substrate. The 004F/007R primer pair targets telomeric
155
repeat sequences and should produce 665 bp and 1230 bp products in reactions
156
containing TP03 and DPP17 DNAs, respectively. After cloning, the viruses were purified
157
and sequenced as described. Four DTD clones (Dryvax-TianTan DNA reactivation) were
158
sequenced.
159
160
Virus sequencing and genomic analysis. Stocks of 16 DTM clones, 15 DTH, and 4
161
DTD clones were prepared and purified over sucrose gradients. Viral DNAs were
162
extracted and sequenced as described previously (9) using a Roche 454 GS Junior system.
163
Roche GS De Novo Assembler software was used to deconvolute and assemble the raw
164
sequencing data into contigs and nearly full-length genomes were generated using CLC
165
Genomics Workbench 6. The average read redundancy was 15, which permitted the
166
assembly and mapping of all of the recombinant junctions with confidence. Multiple
167
sequence alignments were prepared using the program LAGAN (http://genome.lbl.gov)
Qin et al.
Vaccinia virus recombination
Page 7 of 31
168
(25) and Base-by-Base software (26) was used to produce a visual summary of the whole
169
genome alignments. The assembled sequence data have been assigned GenBank
170
accession numbers KJ467582 to KJ467616, inclusive (Table S1).
171
PCR and Southern blotting. Southern blotting was used to confirm the rearrangement
172
detected in clone DTM28. Virus DNA was digested with NdeI (Fermentas) and size
173
fractionated by electrophoresis on 0.7% agarose gels. The DNA was fragmented with
174
0.2M HCl, denatured with 0.4M NaOH, transferred to a nylon membrane, and fixed with
175
a UV cross-linker. Two primers (201F and 239R, Table 1) and the PCR were used to
176
prepare a probe in a reaction containing biotin-16-dUTP (Roche, 1093070), which was
177
subsequently hybridized to the prepared membrane and detected using IRDye
178
800CW-coupled streptavidin (Li-Cor; 926-32230) and a Li-Cor imager. The DTM28 and
179
DTM28Δ viruses were also differentiated using the PCR and 208F and 209R (Table 1)
180
primers. This was done in combination with 201F and 239R in a PCR reaction containing
181
all four primers.
182
183
RESULTS and DISCUSSION
184
185
Virus isolation and genome sequencing and assembly. We used three different methods
186
to produce VACV recombinants. The first method was designed to explore the effects of
187
repeated passage, at low MOI, on a seed mixture initially composed of just two different
188
genetically tagged viruses. These viruses were originally cloned from stocks of Dryvax
189
(DPP17) and TianTan (TP05) vaccines and differ in sequence by 1 SNP per ~140 bp.
190
Compared to TP05, DPP17 also encodes a 6 kbp deletion near the right terminal inverted
191
repeat boundary as well as ~150 other smaller insertions and deletions (indels) distributed
192
across the two genomes (Figure 1). The TP05 strain forms plaques that are approximately
Qin et al.
Vaccinia virus recombination
Page 8 of 31
193
twice the diameter of those formed by the DPP17 strain, which provided an opportunity
194
to explore what effect a growth bias might have on the pattern of recombinants.
195
For this first experiment, a 1:1 mixture of the two different VACV were used to
196
infect BSC-40 cells at MOI=0.02, cultured 48 hr, and a portion (10 μL or ~0.02 PFU/cell)
197
of the resulting progeny passaged again under the same conditions. This was repeated to
198
produce a total of 5 rounds of replication. Each time the infection produced overlapping
199
plaques that partly cleared the entire plate. We then plated out the diluted virus in 24-well
200
plates, and identified 36 wells each containing just a single random plaque. These 36
201
viruses were then cloned again, also by limiting dilution, and designated DTM
202
(Dryvax-TianTan mixture) strains. Using this method minimized the risk of picking
203
certain plaque types, since the only criteria we used to choose a clone was that the virus
204
had to have been diluted to the point where it was the only plaque in a well, in the first
205
round of selection. To avoid the problem of resequencing any non-recombinant parental
206
stains, the PCR and three different primer pairs were used to determine the genetic origin
207
of three different sites within each genome: within the terminal inverted repeats (primers
208
004F/007R), in the central part of the genome (primers 107F/108R), and near the junction
209
with the right terminal inverted repeat (213F/209F/226R) (Table 1). Fourteen clones were
210
selected because at least one position was recombinant with respect to either of the other
211
two sites. We also chose two additional viruses, which exhibited a parental arrangement
212
of markers at these three sites, although these were subsequently determined to also be
213
recombinants. These viruses were cloned two more times and 16 were sequenced. After
214
sequencing and assembly, these DTM recombinants exhibited a patchy pattern of SNPs
215
suggesting each virus was the product of approximately 30 exchanges over the course of
216
virus replication.
217
One expects that when viruses are passaged five times under these conditions, it
218
should provide an opportunity for repeated rounds of replication and recombination. We
Qin et al.
Vaccinia virus recombination
Page 9 of 31
219
also examined what the virus progeny would look like if they were permitted just a single
220
round of infection, although it is expected that this would still involve multiple rounds of
221
replication. To do this, we co-infected BSC-40 cells with DPP17 and TP05 viruses at
222
MOI=10 (5 pfu/cell of each virus) and cultured the viruses for just 24 hr. These viruses
223
were cloned and designated DTH strains (Dryvax-TianTan high MOI). After the first
224
round of cloning, 43 viruses were randomly selected and the PCR was used to identify
225
putative recombinants as described above. Thirteen hybrid viruses were cloned twice
226
more, along with two additional viruses (DTH13 and DTH14) that exhibited a parental
227
pattern of markers at the three positions tested by the PCR. DTH14 was subsequently
228
identified as being identical to the TP05 parent virus, while DTH13 proved to be a
229
recombinant. Ultimately, 15 DTH clones were sequenced and assembled as described
230
before. These recombinants exhibited a mean of 18 exchanges per genome.
231
We should note one caveat regarding these methods, in that single plaques isolated in
232
the first round of purification were not always pure, and this provided some limited
233
opportunity for additional rounds of replication and recombination. For example, when a
234
plaque initially identified as recombinant DTM22 was cloned a second time, and the
235
subsidiary plaques reanalyzed by PCR, it was realized that the two daughter plaques
236
(DTM22.1 and DTM22.2) were not identical. However, they are clearly “sibs”, viruses
237
sharing a common genetic origin as judged by a shared pattern of exchanges in the center
238
of the two genomes (Figure 1). We also noted one case where a single apparently
239
recombinant starting plaque resolved into two clearly unrelated recombinant clones upon
240
replating (clones DTH10 and DTH10.2, Figure 1). For simplicity, our analysis has treated
241
these particular clones as being the same as the other recombinants isolated in the study,
242
although they may have experienced some additional limited opportunities to undergo
243
recombination.
244
Qin et al.
Vaccinia virus recombination
Page 10 of 31
245
Crossovers in DTM and DTH viruses. After assembly, the sequences were aligned with
246
program LAGAN and the alignment corrected manually using Base-by-Base. Inspecting
247
these sequences we could readily identify the origin of each SNP-tagged segment of
248
DNA as belonging to either the DPP17 or TP05 parent (Figure 1). What was remarkable
249
was the very low frequency of observed mutation even though numerous SNPs and
250
smalls indels commonly differentiate clones isolated from a viral stock like Dryvax (9).
251
No mutations were detected in any of the DTH clones, compared with the two parent
252
viruses, and just two mutations in two of the DTM clones. One was a small deletion in
253
DTM29 at alignment position 900, which removed two nucleotides (Figure 2, panel A)
254
just 6 nt upstream of the ORF001 start codon. Although most small deletions are
255
associated with repeats (9, 27), this event was not. It was located immediately adjacent to
256
a SNP that differentiates DPP17 from TP05. We also discovered a point mutation in
257
DTM27 at alignment position 70,493 (Figure 2, panel B). This causes a C-to-T transition
258
mutation and an alanine-to-valine substitution in gene DVX_088 (RNA-helicase). A
259
crude estimate of the VACV replicative error rate can be calculated from the following
260
observations and assumptions. We note that there were only two independent mutations
261
detected in 16 DTM viruses over the course of 5 rounds of infection, and there are
262
~200,000 nt copied per genome per each round of infection. Each round of infection
263
typically expands the VACV titer ~10,000-fold (i.e. between 214 and 215 doublings of the
264
genome) and thus the error rate is very crudely estimated as 2÷[16×200,000×5×14.5] or
265
~1×10-8 mutations per nucleotide copied per cycle of replication. Alternatively this is
266
2÷[16×5×14.5] or ~1 mutation per 600 genomes per cycle of replication. By a similar
267
method, the absence of any mutations detected amongst the 15 DTH clones over the
268
course of a single round of infection suggests an error rate of <5×10-8. The VACV E9L
269
gene encodes a typical B-family proofreading DNA polymerase of a type encoded by a
270
variety of viruses and bacteriophage, and this error frequency resembles that reported for
Qin et al.
Vaccinia virus recombination
Page 11 of 31
271
phage (28). Although some drug-resistant E9L alleles cause altered spontaneous mutation
272
rates in vivo (29-31), for comparison purposes these cannot be converted into absolute
273
mutation rates given the uncertainties in the size or number of genetic target(s). We could
274
find no other reported absolute error rates for poxviruses in the literature.
275
Beyond these two rare mutations, the remainder of the sequences in the recombinant
276
genomes could be ascribed to having been inherited from one or the other of the two
277
parent viruses. In total, 1399 SNPs (single nucleotide polymorphisms) can be used to
278
differentiate DPP17 from TP05 and we used these SNPs to track the origins, and thus the
279
sites of crossing over, in the hybrids. The relative abundance of these variant sites (1399
280
scattered across ~200 kbp) allowed us to map the site(s) of crossing over with an average
281
resolution of ~140 nt. In general, each hybrid virus encoded variable-length blocks of
282
DNA derived from each of the two parent strains, and no uniquely conserved block (a
283
hallmark of a highly selected patch of DNA) was detected in all of the viruses. The
284
lengths of these blocks of recombined sequences varied, ranging in length from one to
285
several hundred SNPs. We detected none of the large gene duplications that have been
286
described by other authors (32, 33), but would not expect to do so given that these
287
structures are stable only in the presence of strong selection pressure.
288
To examine the pattern of crossing over in greater detail, we used “Base-by-Base”
289
software to produce a table ascribing each of the 1399 SNPs detected in each hybrid as
290
being of either DPP17 or TP05 origin (data not shown). This provided a tool for
291
calculating the number of exchanges on the assumption that each time the pattern of
292
SNPs changed from DPP17 to TP05 (or vice versa), that an exchange had occurred and
293
was counted as one crossover. It is important to note that this is still an underestimate of
294
the true physical recombination frequency, as any recombination between two parental
295
genomes (e.g. DPP17×DPP17 or TP05×TP05), or recombination occurring over distances
296
less than the distances between SNPs cannot be detected by this method. Figure 3 shows
Qin et al.
Vaccinia virus recombination
Page 12 of 31
297
the results of this analysis. We had expected that viruses given more opportunities for
298
recombination would exhibit a greater number of exchanges, and this was supported by
299
these measurements. The number of crossovers in the DTM viruses ranged from 14 to 44
300
(mean = 30±11 [SD]). In contrast, the numbers of crossovers in the DTH viruses ranged
301
from 0 to 38 (mean = 18±11 [SD]), with one virus, DTH14, identical to the TP05 parents.
302
In a previous publication (13) we conducted a meta-analysis of all of the published
303
VACV genetic recombination data and attempted to correlate these data with the known
304
physical map of the virus. From these studies we derived an estimate that half-maximal
305
recombination was detected in a single round of infection (akin to the method used to
306
produce the DTH viruses), when the distance between VACV markers was ~8 kbp.
307
Although this number is difficult to estimate with precision due to a number of
308
experimental factors (13), from our measurements of the conversion track lengths (see
309
below) we calculated that the DTH viruses exhibited a mean of about 1 physical
310
crossover per 12 kbp, a number compatible with the ~8 kbp deduced from the older
311
genetic literature. Historically, classical genetic crosses never proved a useful tool for
312
mapping VACV genomes, due to a combination of experimental noise and short linkage
313
distances relative to the distances between most VACV markers. When one considers that
314
our estimate of 1 crossover per 12 kbp is associated with a standard deviation of 19 kbp
315
(i.e. 12±19 kbp), the source of this problem is clearly apparent.
316
317
Average length of the conversion tracts. Besides measuring the numbers of exchanges
318
suffered by each of the recombinant viruses, we also examined the lengths of the DNA
319
segments exchanged between viruses (i.e. the conversion track length) in the DTH group.
320
To do this, the calculation assumed that the start and end of each exchange lay midway
321
between the SNPs flanking the two sites of exchange. The resolution of the method varies
322
depending upon the local SNP density, but with an average of 1 SNP per 140 bp we could
Qin et al.
Vaccinia virus recombination
Page 13 of 31
323
detect exchanges ranging in size from 55 to 92,000 bp. An interesting feature of VACV
324
recombination is illustrated by this analysis, which showed that there were relatively
325
more short conversion tracks than long ones (Figure 4). Thus while the mean length of a
326
conversion track was 12 kbp, the median was only 2.6 kbp. The abundance of short
327
conversion tracks would help favor intragenic recombination events, which can be
328
detected between markers spaced only 54 bp apart (34). We should note that this estimate
329
of the recombination frequency is lower than has been previously reported. For example
330
we detected a loss of linkage at distances exceeding 350 bp in one study (35). However,
331
this earlier experiment measured the yield of recombinants when DNA was transfected
332
into Shope fibroma virus-infected cells, and it is possible that the non-specific DNA
333
replication that is seen under these circumstances (36) also exposes transfected DNAs to
334
higher levels of recombination than is normally experienced by viruses.
335
Crossing over is not the only process that could produce this abundance of short
336
exchanges. Poxvirus replication and recombination reactions also produce hybrid (or
337
heteroduplex) DNA (37). Such molecules would contain mismatched bases wherever the
338
sequences differ and if a subset of mismatches were subjected to mismatch-specific and
339
directionally biased short patch DNA repair prior to further replication, it could create the
340
appearance of closely linked crossovers. However, our data provide no evidence that
341
mismatch repair is producing artifactual exchanges. If any repair bias existed, it would
342
most probably occur at sites where the hybrid DNA contained G·T mismatches, since G·T
343
mismatches are generally repaired to G·C basepairs in cells by a pathway employing a
344
thymine DNA glycolylase (38). A G·T mismatch would be formed (along with a
345
reciprocal A·C mismatch) at sites where substitution mutations differentiate the two
346
viruses. Such substitutions (i.e. GÆA, AÆG, CÆT, or TÆC) comprise 72% of the SNP
347
markers in these crosses. If one makes the simplifying assumption that A·C mismatches
348
are just replicated and not repaired, and that all G·T mismatches are converted to G·C
Qin et al.
Vaccinia virus recombination
Page 14 of 31
349
prior to replication, then these single marker exchanges should be biased 2:1 in favour of
350
forming (or retaining) a G or C. We detected 51 single exchanges at sites containing base
351
substitutions (82% of all single exchanges), in the DTM and DTH viruses. Of these, 29
352
retained a G or C, and 22 retained an A or T, which is not significantly different from a
353
1:1 split (Ȥ2=0.96, P=0.33). Although we cannot disprove the hypothesis that biased
354
mismatched repair created some of the short conversion tracts, the simplest explanation
355
for these data is that poxvirus recombination reactions produce an abundance of short
356
conversion tracks through a process formally akin to crossing over.
357
358
Biased genetic origins in progeny viruses. An interesting difference between the
359
TianTan and Dryvax clones used in this study is that TP05 forms plaques twice the size
360
of DPP17 plaques on BSC40 cells. We wondered how this phenotype might segregate
361
amongst the recombinants deriving from either the DTH or DTM crosses. We used the
362
SNPs to determine what fraction of each genome derived from TP05 or DPP17, and
363
plated all of the cloned viruses on BSC-40 cells, at the same time, to determine the
364
average plaque size. The DTH viruses, passaged just once, showed no particular
365
compositional bias, comprising about equal portions (50±27%) of each of the parental
366
viruses (Figure 5A). In contrast the DTM hybrids bear a diminished (19±11%) fraction of
367
the genome derived from DPP17 SNPs (Figure 5A). Oddly, there seems to be a simple
368
linear relationship between plaque size and the proportion of the genome derived from
369
each parent strain, with larger plaque sizes associated with a greater proportion of
370
TP05-derived DNA (Figure 5B). These data suggest that the TP05-derived DNA may
371
confer a selective growth advantage in multiple rounds of culture (i.e. DTM viruses), but
372
one round of growth (i.e. DTH viruses) provides insufficient time or selective pressure to
373
bias the composition of the recombinants.
Qin et al.
Vaccinia virus recombination
Page 15 of 31
374
What would produce this effect is not clear. Plaque size is likely determined by many
375
different genetic factors, and there are many differences in the gene composition of the
376
two parent strains. For example DPP17 contains a large deletion in the right TIR
377
compared to TP05 and this deletion bears a number of different genes (Table 2). However,
378
this deletion is not completely responsible for plaque variation since DPP25, containing
379
all the genes in this region, forms plaques only slightly larger than DPP17 (9). We
380
subsequently annotated all of the hybrid genomes using GATU (39) and evaluated the
381
differences (Table 3). There are many mutant genes segregating in complex ways
382
between the different viruses including mutant forms of I4L (40), F3L (41), E5R (42),
383
M1L, A51R, and C23L, but no obvious distribution patterns could be discerned by
384
inspection beyond the fact that the more presumably functional genes the virus encoded
385
(Table 3, grey cells), the better it grew. The non-transcribed telomeric repeats in VAC
386
telomeres cannot be assembled into contiguous sequences, due to the redundancies in the
387
repeats, but some of the fragments of junction sequences differ enough in TianTan and
388
Dryvax to deduce the origins of the telomeres. This analysis detected a trend suggesting
389
that viruses bearing TP05-like telomeres also form larger plaques. However, all that
390
could really be concluded from this analysis is that TianTan-derived sequences generally
391
contribute greater advantage in culture than does Dryvax DNA.
392
393
Large deletions formed through illegitimate recombination. Poxviruses are also
394
known to suffer deletion mutations during passage. The most extreme example is
395
probably modified vaccinia virus Ankara (MVA), which accumulated six large deletions,
396
and many smaller ones, when it was passaged >570 times in chicken embryo fibroblasts
397
(43). Over the course of these experiments we did detect one such large deletion mutation
398
when we sequenced clone DTM28. The deletion spans 21 kbp and encompasses genes
399
DVX_201 to 239 (Figure 6). In the initial assembly, we found 11 sequence reads that
Qin et al.
Vaccinia virus recombination
Page 16 of 31
400
started in gene DVX_201 and terminated in gene DVX_239 (Figure 6a), along with
401
sequence reads derived from all of the intervening genes. The deletion spans the right
402
TIR boundary, but amongst the reads were some from the unique genes DVX_202-209
403
suggesting we had sequenced a stock of virus containing the DTM28 parent contaminated
404
by a virus bearing the deletion (DTM28Δ). To confirm this interpretation of the data we
405
prepared primers 201F and 239R (which are located 21 kbp apart in genes DVX_201 and
406
DVX_239, respectively; Figure 6b) and used the PCR to detect the novel 1.2 kbp
407
amplicon that was predicted to be formed in this process (Figure 6c). We also tested
408
DNAs extracted from the virus stocks that had been archived during the process of
409
passaging these viruses 5 times, before cloning, as well as DTM27, another independent
410
clone that was purified in parallel. Only the purified DTM28 stock contained a virus
411
bearing the deletion (Figure 6c) suggesting that DTM28Δ arose during the expansion of
412
the stock. Finally we used the 1.2 kbp amplicon to probe a Southern blot of NdeI-cut
413
virus DNA, and showed that the DTM28 stock contains viruses contributing a 6.7kbp
414
band characteristic of the deletion-containing fragment as well as a 5.4 kbp band, which
415
derives from the two 5.4 kbp NdeI fragments that encode the boundaries of the deletion
416
(Figure 6d). We subsequently sub-cloned this stock and separately isolated the two
417
viruses, confirming the viability of DTM28Δ and the fact that none of the deleted genes
418
are essential (Figure 6e) in cell culture.
419
The DNA surrounding the vaccinia virus right TIR boundary is a well-established
420
hotspot for large deletion mutations (23). The mechanism is probably the same as that
421
which drives the formation of small deletion mutations, stating with the misalignment of
422
regions containing imperfect repeats (9, 27). If one aligns the reads spanning the junction
423
boundary between DVX_201 and DVX_239, one sees several small blocks of homology
424
(Figure 6a, boxed) that could have stabilized the first step in an illegitimate
425
recombination reaction. It is difficult to establish an exact rate by which such mutations
Qin et al.
Vaccinia virus recombination
Page 17 of 31
426
arise, but this stock was plaque purified three times, following bulk up and only one virus
427
in 16 DTM viruses passaged in parallel suffered a deletion of this type. This creates a rate
428
of ~0.06/4 passages or 1 deletion per 70 passages. The six large deletions introduced into
429
MVA over 570 passages are thus quite consistent with this estimate although, of course,
430
the selection pressures are very different in the two experiments.
431
432
Recombination in SFV-reactivated vaccinia viruses. A third small collection of
433
recombinant viruses was also produced using Shope fibroma virus mediated DNA
434
reactivation assays and DNAs extracted from DPP17 and TianTan strain TP03 (22). [We
435
used TP03, instead of the TP05 used for the preceding experiments, to test whether
436
viruses could also be produced encoding all three of the large TP03 and DP17 telomeric
437
deletions.] This method relies upon a replicating helper virus (SFV) to rescue or
438
“reactivate” fragments of transfected virus DNA (VACV). The SFV is subsequently
439
eliminated by passage on a cell line that supports only VACV growth. Figure 7 shows the
440
maps of the viruses that were recovered by this method. There were just four viruses
441
obtained and two (DTD03 and DTD11) are so similar that they are probably “sibs”
442
sharing a mostly common history. These viruses were too few in number, and the passage
443
history too complicated, to derive much in the way of statistics about recombination
444
patterns, but the pattern of exchanges generally resembled the lesser numbers and longer
445
conversion tracks seen in the DTH viruses. The method did also produce clones bearing
446
the three large telomeric deletions (DTD03 and DTD18, Figure 7), which left the virus
447
with TIRs just 7.3 kbp long. An important caveat is that the DTD viruses were recovered
448
from cells that had been transfected for a few hours and then incubated for three days, so
449
whether the recombinants were produced during the reactivation stage or during
450
subsequent rounds of re-infection and replication is difficult to deduce.
Qin et al.
Vaccinia virus recombination
Page 18 of 31
451
Overall, there were no strikingly unique features of these reactivated viruses that
452
would differentiate them from any other type of recombinant poxvirus. Perhaps the most
453
important conclusion that could be drawn from this brief study is that this process is very
454
accurate (no mutant viruses were recovered) and no Shope fibroma virus DNA sequences
455
were detected in any of the reactivated VACV. The two most similar genes in SFV and
456
VACV are S068R and J6R, respectively (44), which share only 73% nucleotide sequence
457
identity with no blocks of perfect alignment >17 nt. There is even less similarity between
458
VACV and fowlpox virus, another virus that has also been used to reactivate
459
Orthopoxviruses (45). This is probably insufficient sequence similarity to support
460
frequent recombination between the helper and reactivated viruses. Additionally VACV
461
hybrids may well be rare and difficult to isolate in the absence of selection, if not simply
462
inviable. Such data support the long-standing suspicion that using a heterologous helper
463
virus, like SFV or fowlpox virus, to reactivate Orthopoxviruses can be done without
464
mutation and does not produce hybrid strains.
465
466
Conclusions. Next generation DNA sequencing technologies are greatly improving our
467
understanding of the genome structures and genes encoded by large DNA viruses. Here
468
we show that these methods can also be used to characterize the structures of
469
recombinant poxviruses. These studies show that recombinant VACV are not surprisingly
470
composed of a patchwork of DNA fragments derived from the parent viruses. The
471
numbers of exchanges varies depending upon the passage history, but if one uses
472
methods like those classically used to produce VACV recombinants (a high multiplicity
473
of infection [10] and one day of co-culture), one detects about 1 physical crossover per 12
474
kbp in the DTH viruses, a number only slightly higher than the ~8 kbp we have estimated
475
from a review of the older genetic literature. However, there is a lot of noise observed in
Qin et al.
Vaccinia virus recombination
Page 19 of 31
476
this number (12±19 kbp), perhaps explaining why accurate classical recombination maps
477
were never produced for VACV.
478
Interestingly the lengths of the recombinant patches (i.e. the conversion tracts) are
479
heavily biased towards shorter sizes, something that would favour intragenic
480
recombination. What mechanism would produce such an effect is difficult to identify,
481
although we have previously used genetic methods to show that VACV replication and
482
recombination are intimately linked processes (18, 46), probably because the VACV E9
483
DNA polymerase exhibits properties characteristic of a recombinase both in vitro (21)
484
and in vivo (20). Thus recombination may just be an indirect by-product of virus
485
replication, conceivably associated with the DNA polymerase-catalyzed repair of broken
486
replication structures. Regardless of the mechanism, this process could have interesting
487
genetic consequences for virus evolution, as it would create a lot of diversity within
488
recombinant genes, not just diverse combinations of different genes. This becomes of
489
critical importance when one considers the challenge posed to viruses by rapidly evolving
490
responses to biological features like immunodominant epitopes. Short conversion tracks
491
offer a selective advantage for a virus, as they provide a mechanism for rearranging and
492
eliminating peptide epitopes while still retaining gene function.
493
These studies also show how sequencing could be used to characterize more complex
494
virus traits than those regulated by single genes. Continued passage of the DTM viruses
495
selected for viruses bearing greater proportions of the TianTan genome and this was
496
correlated in some, still unclear, manner to plaque size. By producing recombinants,
497
applying a selection strategy (perhaps in an iterative manner), and then sequencing clones
498
bearing the desired traits, it should be possible to map genes that collectively regulate the
499
phenotype of interest. This is not a novel approach of course; related methods have been
500
used for decades to map complex genetic traits in many different organisms. However,
501
the widespread availability of next generation sequencing technologies creates a tool that
Qin et al.
Vaccinia virus recombination
Page 20 of 31
502
could easily be used by many more laboratories studying gene families and gene
503
interactions in large DNA viruses.
504
Acknowledgements. We thank Dr. Wendy Magee and Mr. Rob Maranchuk for their
505
technical supports with the Roche 454 Junior sequencer. This work was supported by
506
operating grants from the Canadian Institutes for Health Research and the Alberta Cancer
507
Foundation and by an infrastructure award from the Canada Foundation for Innovation.
Qin et al.
Vaccinia virus recombination
Page 21 of 31
508
FIGURE LEGENDS
509
510
Figure 1. Patterns of DNA exchange in recombinant vaccinia viruses. The genome
511
sequences of DTM (panel A) and DTH (panel B) recombinant clones were aligned
512
against the parent genomes DPP17 and TP05 using the program “LAGAN”, and edited
513
using the program “Base-by-Base”. TP05 was used as the reference strain and any
514
differences between a given virus and TP05 are colour coded to indicate insertion,
515
substitution, and deletion mutations derived from strain DPP17 (inset at bottom). Thus
516
the blank regions represented fragments derived from TP05.
517
518
Figure 2. Rare mutations associated with replication and recombination. Panel A
519
shows a deletion mutation in DTM29, immediately adjacent to a G-to-T SNP found at
520
alignment position 70493. Panel B shows a C-to-T substitution detected only in clone
521
DTM27 at position 900.
522
523
Figure 3. Numbers of exchanges in DTM and DTH clones. Each of the hybrid viruses
524
was first aligned against strains TP05 and DPP17. Then, the program “Base-by-Base”
525
was used to determine where each cross-over was located relative to the 1399 SNPs that
526
differentiate the two strains, along with the number of such exchanges. The viruses
527
passaged five times (DTM) exhibited more exchanges per genome than the viruses
528
passed just once (DTH). That is 30±11 versus 18±11 exchanges/genome, respectively.
529
530
Figure 4. Length of the DNA segments exchanged in DTH clones. The lengths of all
531
the conversion tracks were measured in all 14 genomes using midpoints defined by the
532
four SNPs flanking the two bounding sites of exchange. The numbers of events, of a
533
given exchange length, were then determined by assignment to 200 bp bins. A
Qin et al.
Vaccinia virus recombination
Page 22 of 31
534
semi-logarithm of the bin size (i.e. exchange length) is presented because the values
535
differ so greatly in scale across the different genomes. VACV recombination appears to
536
be associated with a disproportionate number of very short exchanges. Because there is
537
approximately one SNP per 140 bp, greater resolution than ~200 bp is not achievable.
538
539
Figure 5. Biased selection for sequences associated with the TianTan parent. The
540
percent of DNA derived from each of the parental viruses was determined from the
541
fraction of SNPs derived from each parent. Panel A shows how the composition varied in
542
viruses passaged just once (DTH hybrids) or five times (DTM hybrids) prior to cloning.
543
Passage appeared to select for SNPs linked to the TP05 TianTan parent, as the percentage
544
of Dryvax DNA decreased from 50±27% to 19±11% with continued passage. Panel B
545
illustrates how the plaque size is related to the genetic origins of the hybrid. The viruses
546
forming smaller plaques, more closely resemble the DPP17 parent. To measure the
547
plaque size, each of the cloned hybrids was plated on BSC-40 cells (in parallel), cultured
548
for two days, stained with crystal violet, scanned, and the plaque area determined using
549
ImageJ(24). Twenty randomly selected plaques were measured for each virus.
550
551
Figure 6. Illegitimate recombination detected during the cloning and sequencing of
552
hybrid DTM28. During the sequencing of DTM28, 11 reads were detected linking gene
553
DVX_201 to gene DVX_239. Panel A shows an alignment of these reads to sequences
554
within the two genes, which are normally spaced 21 kbp apart. We have also identified
555
sequence identities (circles), short patches of homology (boxed) and a simple repeat
556
(underlined) common to sequences flanking the fusion site. The sequence in these reads
557
transitions cleanly from one gene to the next, with no evidence of any unrelated
558
additional sequences having been added in the process. Panel B showed a way to form
559
this deletion. Illegitimate recombination between identical parents (DTM28) excised
Qin et al.
Vaccinia virus recombination
Page 23 of 31
560
21kbp and created the virus we subsequently called DTM28Δ. Panel C shows the results
561
of a PCR analysis using primers targeting sites flanking the fusion site. These are located
562
too far apart in the parent viruses (e.g. DTM27) to amplify 21 kpb of intervening
563
sequence. The DTM28Δ virus was probably formed during the expansion of the clone
564
prior to sequencing, as it is not detected in intermediary viruses during the course of
565
passages. Panel D shows a Southern blot of NdeI-digested virus DNA showing that the
566
sequenced virus stock contained two viruses. These are the DTM28 hybrid (indicated by
567
a 5.4 kbp fragment common to both parent strains), and the DTM28Δ
568
(indicated by a 6.7 kbp fragment containing the fusion junction). Panel E shows the
569
DTM28Δ is independently viable. Six randomly selected viruses were separately plaque
570
purified from the sequenced stock and the PCR was used to detect sequences found only
571
in the deleted region in DTM28 (primers 208F+209R) or only capable of being amplified
572
if the intervening sequences are deleted (primers 201F+239R).
573
574
Figure 7. Patterns of DNA exchange in recombinant vaccinia viruses produced using
575
Leporipoxvirus-mediated reactivation reactions. The genome sequences of the DTD
576
recombinant clones were aligned against the parent genomes DPP17 and TP03 using the
577
program “LAGAN”, and edited using the program “Base-by-Base”. Because this
578
experiment used TP03 DNA, and TP05 was always used as the reference strain in all of
579
our analyses (Figure 1), the telomeric deletion mutations that differentiate TP05 from
580
TP03 show up as additional red blocks in the TP03 alignment.
581
Qin et al.
Vaccinia virus recombination
Page 24 of 31
582
TABLES
583
584
Table 1. PCR primers used in this study
585
Amplicon size (bp)
DPP17
TP05
1120
-
-
653
1230
665
225
260
712
712
Primer ID
Primer (5’-to-3’)
DVX-209F
CGAAGAAGATGATGGGGAC
DVX-226R
ATAAGAGGAAAGAGGACAC
DVX-213F
CGTTGGATGGATTCGATA
DVX-226R
ATAAGAGGAAAGAGGACAC
DVX-004F
GCAGTAGGCTAGTATCTT
DVX-007R
TACCGGCATCATAAACAC
DVX-107F
AACTGGAGTAGAGATAGC
DVX-108R
CCGAGAATATAGCTGTCC
201F
AATATGATGGTGATGAGCGA
239R
TATTGCGAGATGTGAAGG
208F
TTTCTTCTCTTCTCCCTTTC
209R
ATTCTATCCCGTACCTCT
586
Qin et al.
Vaccinia virus recombination
Page 25 of 31
587
Table 2. Gene differences in the right TIR deletion: TP05 versus DPP17.
588
Gene/virus
Gene (or feature)
Nucleotide length (bp)
Copenhagen
Dryvax
DPP17
TP05
IL-1-β-receptor
Cop-B16R
DVX_209
1002
981
Unknown
Cop-B17L
DVX_210
-
1023
Ankyrin motifs
Cop-B18R
DVX_211
-
1725
IFN-α/β-receptor
Cop-B19R
DVX_212
-
1056
Ankyrin motifs
Cop-B20R
DVX_213
-
1794
589
Qin et al.
Vaccinia virus recombination
Page 26 of 31
590
Table 3. Gene complement in parent and hybrid viruses
Virus
DPP17
DTH13
DTH22
DTH36
DTH21
DTH31
DTH08
DTH34
DTM29
DTH06
DTH27
DTH10.2
DTH41
DTM03
DTM10
DTH30
DTH28
DTM32
DTM28
DTM27
DTH10
DTM30
DTM11
DTM22.2
DTM09
DTM33
DTM22.1
DTM06
DTM04
DTM19
DTH14
DTM08
TP05
Qin et al.
Plaque
area
100%
101
103
108
111
112
115
118
124
126
131
137
137
140
141
144
146
150
156
158
160
162
168
175
176
178
179
181
186
188
199
204
206
Telo
mere
D1
D
D
D
T2
T
D
D
T
D
T
D
T
D
D
D
T
D
D
T
T
T
T
D
D
T
T
T
T
T
T
T
T
right
TIR
-3
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
F4L
C23L
M1L
F3L
E5R
A51R
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+4
+
+
+
+
+
+
+
+
+
+
+
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Vaccinia virus recombination
Page 27 of 31
591
1. D=DPP17-like telomere repeats; 2. T=TP05-like telomere repeats; 3. Gene or genes
592
are truncated or deleted; 4. Gene appears intact (gray shading).
593
REFERENCES CITED
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
Hershey AD, Rotman R. 1948. Linkage Among Genes Controlling Inhibition of
Lysis in a Bacterial Virus. Proc Natl Acad Sci U S A 34:89-96.
Wildy P. 1955. Recombination with herpes simplex virus. J Gen Microbiol
13:346-360.
Fenner F, Comben BM. 1958. [Genetic studies with mammalian poxviruses. I.
Demonstration of recombination between two strains of vaccina virus]. Virology
5:530-548.
Woodroofe GM, Fenner F. 1960. Genetic studies with mammalian poxviruses.
IV. Hybridization between several different poxviruses. Virology 12:272-282.
Bedson HS, Dumbell KR. 1964. Hybrids Derived from the Viruses of Variola
Major and Cowpox. The Journal of hygiene 62:147-158.
Bedson HS, Dumbell KR. 1964. Hybrids Derived from the Viruses of Alastrim
and Rabbit Pox. The Journal of hygiene 62:141-146.
Strayer DS, Skaletsky E, Cabirac GF, Sharp PA, Corbeil LB, Sell S, Leibowitz
JL. 1983. Malignant rabbit fibroma virus causes secondary immunosuppression in
rabbits. J Immunol 130:399-404.
Esposito JJ, Sammons SA, Frace AM, Osborne JD, Olsen-Rasmussen M,
Zhang M, Govil D, Damon IK, Kline R, Laker M, Li Y, Smith GL, Meyer H,
Leduc JW, Wohlhueter RM. 2006. Genome sequence diversity and clues to the
evolution of variola (smallpox) virus. Science 313:807-812.
Qin L, Upton C, Hazes B, Evans DH. 2011. Genomic analysis of the vaccinia
virus strain variants found in Dryvax vaccine. J Virol 85:13049-13060.
Ensinger MJ, Rovinsky M. 1983. Marker rescue of temperature-sensitive
mutations of vaccinia virus WR: correlation of genetic and physical maps. J Virol
48:419-428.
Ensinger MJ. 1982. Isolation and genetic characterization of
temperature-sensitive mutants of vaccinia virus WR. J Virol 43:778-790.
Condit RC, Motyczka A, Spizz G. 1983. Isolation, characterization, and physical
mapping of temperature-sensitive mutants of vaccinia virus. Virology
128:429-443.
Lin YC, Evans DH. 2010. Vaccinia virus particles mix inefficiently, and in a way
that would restrict viral recombination, in coinfected cells. J Virol 84:2432-2443.
Qin et al.
Vaccinia virus recombination
Page 28 of 31
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
Panicali D, Paoletti E. 1982. Construction of poxviruses as cloning vectors:
insertion of the thymidine kinase gene from herpes simplex virus into the DNA of
infectious vaccinia virus. Proc Natl Acad Sci U S A 79:4927-4931.
Mackett M, Smith GL, Moss B. 1982. Vaccinia virus: a selectable eukaryotic
cloning and expression vector. Proc Natl Acad Sci U S A 79:7415-7419.
Ball LA. 1987. High-frequency homologous recombination in vaccinia virus
DNA. J Virol 61:1788-1795.
Ball LA. 1995. Fidelity of homologous recombination in vaccinia virus DNA.
Virology 209:688-691.
Evans DH, Stuart D, McFadden G. 1988. High levels of genetic recombination
among cotransfected plasmid DNAs in poxvirus-infected mammalian cells. J
Virol 62:367-375.
Condit RC. 2007. Vaccinia, Inc.--probing the functional substructure of poxviral
replication factories. Cell Host Microbe 2:205-207.
Gammon DB, Evans DH. 2009. The 3'-to-5' exonuclease activity of vaccinia
virus DNA polymerase is essential and plays a role in promoting virus genetic
recombination. J Virol 83:4236-4250.
Willer DO, Yao XD, Mann MJ, Evans DH. 2000. In vitro concatemer formation
catalyzed by vaccinia virus DNA polymerase. Virology 278:562-569.
Yao XD, Evans DH. 2003. High-frequency genetic recombination and
reactivation of orthopoxviruses from DNA fragments transfected into
leporipoxvirus-infected cells. J Virol 77:7281-7290.
Qin L, Liang M, Evans DH. 2013. Genomic analysis of vaccinia virus strain
TianTan provides new insights into the evolution and evolutionary relationships
between Orthopoxviruses. Virology 442:59-66.
Abramoff MD, Magalhaes, P.J., Ram, S.J. 2004. Image Processing with
ImageJ. Biophotonics International 11:36-42.
Brudno M, Do CB, Cooper GM, Kim MF, Davydov E, Green ED, Sidow A,
Batzoglou S. 2003. LAGAN and Multi-LAGAN: efficient tools for large-scale
multiple alignment of genomic DNA. Genome Res 13:721-731.
Brodie R, Smith AJ, Roper RL, Tcherepanov V, Upton C. 2004.
Base-By-Base: single nucleotide-level analysis of whole viral genome alignments.
BMC Bioinformatics 5:96.
Coulson D, Upton C. 2011. Characterization of indels in poxvirus genomes.
Virus Genes 42:171-177.
Goodman MF, Fygenson KD. 1998. DNA polymerase fidelity: from genetics
toward a biochemical understanding. Genetics 148:1475-1482.
Qin et al.
Vaccinia virus recombination
Page 29 of 31
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
Taddie JA, Traktman P. 1991. Genetic characterization of the vaccinia virus
DNA polymerase: identification of point mutations conferring altered drug
sensitivities and reduced fidelity. J Virol 65:869-879.
Taddie JA, Traktman P. 1993. Genetic characterization of the vaccinia virus
DNA polymerase: cytosine arabinoside resistance requires a variable lesion
conferring phosphonoacetate resistance in conjunction with an invariant mutation
localized to the 3'-5' exonuclease domain. J Virol 67:4323-4336.
Andrei G, Gammon DB, Fiten P, De Clercq E, Opdenakker G, Snoeck R,
Evans DH. 2006. Cidofovir resistance in vaccinia virus is linked to diminished
virulence in mice. J Virol 80:9391-9401.
Slabaugh MB, Roseman NA, Mathews CK. 1989. Amplification of the
ribonucleotide reductase small subunit gene: analysis of novel joints and the
mechanism of gene duplication in vaccinia virus. Nucleic Acids Res
17:7073-7088.
Elde NC, Child SJ, Eickbush MT, Kitzman JO, Rogers KS, Shendure J,
Geballe AP, Malik HS. 2012. Poxviruses deploy genomic accordions to adapt
rapidly against host antiviral defenses. Cell 150:831-841.
Fathi Z, Dyster LM, Seto J, Condit RC, Niles EG. 1991. Intragenic and
intergenic recombination between temperature-sensitive mutants of vaccinia
virus. J Gen Virol 72 ( Pt 11):2733-2737.
Parks RJ, Evans DH. 1991. Effect of marker distance and orientation on
recombinant formation in poxvirus-infected cells. J Virol 65:1263-1272.
DeLange AM, McFadden G. 1986. Sequence-nonspecific replication of
transfected plasmid DNA in poxvirus- infected cells. Proc Natl Acad Sci U S A
83:614-618.
Fisher C, Parks RJ, Lauzon ML, Evans DH. 1991. Heteroduplex DNA
formation is associated with replication and recombination in poxvirus-infected
cells. Genetics 129:7-18.
Maiti A, Drohat AC. 2011. Dependence of substrate binding and catalysis on pH,
ionic strength, and temperature for thymine DNA glycosylase: Insights into
recognition and processing of G.T mispairs. DNA Repair (Amst) 10:545-553.
Tcherepanov V, Ehlers A, Upton C. 2006. Genome Annotation Transfer Utility
(GATU): rapid annotation of viral genomes using a closely related reference
genome. BMC Genomics 7:150.
Gammon DB, Gowrishankar B, Duraffour S, Andrei G, Upton C, Evans DH.
2010. Vaccinia virus-encoded ribonucleotide reductase subunits are differentially
required for replication and pathogenesis. PLoS Pathog 6:e1000984.
Qin et al.
Vaccinia virus recombination
Page 30 of 31
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
41.
42.
43.
44.
45.
46.
Froggatt GC, Smith GL, Beard PM. 2007. Vaccinia virus gene F3L encodes an
intracellular protein that affects the innate immune response. J Gen Virol
88:1917-1921.
Murcia-Nicolas A, Bolbach G, Blais JC, Beaud G. 1999. Identification by mass
spectroscopy of three major early proteins associated with virosomes in vaccinia
virus-infected cells. Virus Res 59:1-12.
Antoine G, Scheiflinger F, Dorner F, Falkner FG. 1998. The complete genomic
sequence of the modified vaccinia Ankara strain: comparison with other
orthopoxviruses. Virology 244:365-396.
Willer DO, McFadden G, Evans DH. 1999. The complete genome sequence of
shope (rabbit) fibroma virus. Virology 264:319-343.
Scheiflinger F, Dorner F, Falkner FG. 1992. Construction of chimeric vaccinia
viruses by molecular cloning and packaging. Proc Natl Acad Sci U S A
89:9977-9981.
Willer DO, Mann MJ, Zhang W, Evans DH. 1999. Vaccinia virus DNA
polymerase promotes DNA pairing and strand-transfer reactions. Virology
257:511-523.
Qin et al.
Vaccinia virus recombination
Page 31 of 31