JOURNAL OF VIROLOGY, Sept. 1998, p. 6967–6978 0022-538X/98/$04.0010 Copyright © 1998, American Society for Microbiology. All Rights Reserved.
Vol. 72, No. 9
Recombination in the 59 Leader of Murine Leukemia Virus Is Accurate and Influenced by Sequence Identity with a Strong Bias toward the Kissing-Loop Dimerization Region JACOB GIEHM MIKKELSEN,1 ANDERS H. LUND,1 MOGENS DUCH,1 1,2 AND FINN SKOU PEDERSEN * Department of Molecular and Structural Biology1 and Department of Medical Microbiology and Immunology,2 University of Aarhus, DK-8000 Aarhus, Denmark Received 4 December 1997/Accepted 24 May 1998
Retroviral recombination occurs frequently during reverse transcription of the dimeric RNA genome. By a forced recombination approach based on the transduction of Akv murine leukemia virus vectors harboring a primer binding site knockout mutation and the entire 5* untranslated region, we studied recombination between two closely related naturally occurring retroviral sequences. On the basis of 24 independent template switching events within a 481-nucleotide target sequence containing multiple sequence identity windows, we found that shifting from vector RNA to an endogenous retroviral RNA template during minus-strand DNA synthesis occurred within defined areas of the genome and did not lead to misincorporations at the crossover site. The nonrandom distribution of recombination sites did not reflect a bias for specific sites due to selection at the level of marker gene expression. We address whether template switching is affected by the length of sequence identity, by palindromic sequences, and/or by putative stem-loop structures. Sixteen of 24 sites of recombination colocalized with the kissing-loop dimerization region, and we propose that RNA-RNA interactions between palindromic sequences facilitate template switching. We discuss the putative role of the dimerization domain in the overall structure of the reverse-transcribed RNA dimer and note that related mechanisms of template switching may be found in remote RNA viruses. transfer; however, the preferred transfer site was mapped within the stem structure and did not coincide with the pause site (28). From in vivo studies of recombination between marker gene cassettes of nonviral origin, it appears that the length of template sequence identity is a primary determinant in template switching of the nascent DNA strand (68). The primary dimerization domain, designated the dimer linkage site, has been mapped to the 59 UTR in a region overlapping the retroviral packaging signal. Support for a kissingloop-loop interaction at least during initiation of RNA dimerization has recently been provided through in vitro studies of murine leukemia virus (MLV) (20, 52), avian sarcoma-leukosis virus (19), and HIV (9, 14, 30, 35, 44, 48–50, 57). This kissingloop interaction model proposes base pairing between palindromic loop sequences of the dimer linkage site and subsequent isomerization of the stem-loops (20), thereby generating an interstrand RNA duplex which represents a local antiparallel linkage of the two RNA subunits. However, in vivo investigations suggest that the interaction of palindromic loop sequences is not absolutely required for retroviral replication (4, 8, 22, 31, 47); therefore, alternative yet unknown RNA interactions within or outside the 59 UTR may support dimerization (4, 31). Recent work suggests that the HIV-1 kissing-loop dimerization region may be essential also for optimal proviral DNA synthesis (47). By use of a single-cycle vector transfer protocol, we have previously developed a forced recombination system based on the strongly restricted replication of Akv MLV-derived vectors harboring a mutated primer binding site (PBS) sequence and the 59 244 nucleotides of the wild-type 476-nucleotide Akv 59 UTR. PBS-modified vectors may thus be transferred through reverse transcriptase-mediated recombination with a MLVlike endogenous virus (MLEV) involving either R-U5-medi-
The dimeric retrovirus RNA genome consists of two plusstrand RNAs which are linked by RNA interactions primarily within the highly structured 59 untranslated region (59 UTR). During reverse transcription, genetic recombination occurs frequently (23, 24, 59), thus contributing to viral variability (27) and escape from lethal mutations (60). Recombination may occur during reverse transcription of a heterodimeric genome containing viral RNA of both exogenous and endogenous origin (18). Indeed, various endogenous viral relics in the mouse genome have proved to represent a source of functional sequences that may participate in recombinational mutation repair (12, 13, 36, 38, 39, 45, 56). Reverse transcriptase-mediated recombination has been demonstrated mainly during RNA-directed minus-strand DNA synthesis, in accordance with the proposed models for forced copy choice (10) and minus-strand exchange recombination (11), but may also involve plus-strand DNA assimilation (26) or unconventional template shifting during plus-strand synthesis (39, 53). Mechanisms for template switching during minusstrand DNA synthesis have been extensively studied by in vitro approaches, and the findings show that reverse transcriptase frequently pauses during polymerization of the nascent minus strand (6, 16, 29, 66). Pausing of the enzyme is directed by sequences or secondary structures within the RNA template (17, 29) and may result in enhanced strand transfer at specific pause sites (6, 16, 66). It was therefore expected in a recent in vitro work that pausing at the human immunodeficiency virus type 1 (HIV-1) transactivation response region causes strand * Corresponding author. Mailing address: Department of Molecular and Structural Biology, University of Aarhus, C. F. Moellers Alle´, Bldg. 130, DK-8000 Aarhus, Denmark. Phone: 45 89422614. Fax: 45 86196500. E-mail:
[email protected]. 6967
6968
MIKKELSEN ET AL.
ated second-strand transfer (39) or minus-strand template shift within the 59 UTR (38). In the latter case, we registered a clustering of recombination sites within a narrow region of the 244-bp 59 UTR coinciding with the primary dimerization site, raising the possibility of a combined role of template sequence identity and RNA secondary structure in template switching between naturally occurring retroviral sequences. To assess the pattern and precision of retroviral recombination between highly structured natural viral RNAs, we have by forced recombination studied transduction of vectors carrying the full-length 476-nucleotide Akv 59 UTR. We address in this report (i) whether the overall structure of the entire 59 UTR and more specifically cis-acting elements in the downstream part of the 59 UTR influence the pattern of template switching during minus-strand DNA synthesis, (ii) if homologous recombination is a simple matter of donor and acceptor template sequence homology at the transfer site, and (iii) if reverse transcriptase template shifting is governed by a minimum length of sequence similarity between the nascent minus-strand DNA and the RNA acceptor template. We conclude from our studies that the kissing-loop dimerization domain within a large recombination target sequence consisting of multiple sequence identity windows (SIWs) is a hot spot for recombination. On this basis, we propose that close RNA-RNA interactions in the primary dimerization palindrome facilitate template switching. MATERIALS AND METHODS Vector construction. Vectors pPBSPro244CAkv-neo and pPBSUmu244CAkvneo have been previously described (33, 38). Vector designations refer to the type of PBS present and 59 UTR length. Vectors harboring the complete 476-bp Akv 59 UTR upstream from the Tn5 transposon fragment which encompasses the neomycin phosphotransferase gene (neo) were generated as follows. A fragment containing the 59 long terminal repeat, (LTR), the proline PBS (PBS-Pro), and the 476-bp Akv 59 UTR was PCR amplified from pAKR59 (32), which contains wild-type Akv cis-acting elements and coding regions. The amplified sequence was cloned by standard procedures into the appropriate position of pPBS-Pro244 CAkv-neo; the resulting retroviral vector was designated pPBSPro476CAkv-neo. PBS-Pro was replaced by the PBS-Umu modification by a two-step PCR approach previously described (33) to generate pPBSUmu476CAkv-neo. Vector constructs are illustrated in Fig. 2A. Sequence analysis and cloning of the MLEV 5* UTR. Sequencing of the upstream part of the MLEV 59 UTR was performed by various PCR-based strategies as previously described (38). Among them, a semirandom-primed PCR approach (58) was used to selectively amplify glutamine PBS (PBS-Gln)-containing sequences in viral cDNA prepared from total RNA in virus-containing medium from nontransfected C-2 packaging cells. To obtain the sequence of the 39 part of the 59 UTR, PCR products were sequenced with a primer matching the degenerate primer linker used in the semirandom PCR. The obtained MLEV sequence allowed the design of primers to specifically amplify the MLEV 59 UTR. The 465-bp MLEV 59 UTR was introduced into vector context by PCR connection with a PCR fragment harboring the 59 LTR and PBS of choice and subsequent cloning of the resulting fragment. The vectors generated were designated pPBSProMLEVCAkv-neo and pPBSGlnMLEVCAkv-neo (see Fig. 2A). Cells, transfections, and virus infections. Growth conditions for C-2 (34) and NIH 3T3 cells and transfections of packaging cells and selection for stably integrated vectors have been previously described (33, 38). Briefly, 10 mg of vector DNA was transfected into C-2 packaging cells seeded at 5 3 103 cells per cm2 on the day before transfection. To estimate the level of marker gene expression in vectors with various 59 UTR lengths, G418-resistant colonies were counted before pooling. Virus infection and determination of transductional efficiencies were carried out as previously described (38). Briefly, serially diluted virus-containing medium was transferred to NIH 3T3 cells; after 10 days of selection of recipient cells, resistant colonies were counted, individually isolated, and expanded. In some experiments, G418-resistant NIH 3T3 colonies on each plate were pooled to allow for a PCR-based screening among a large number of transduction events. In a transduction series set up for colony pooling, viruscontaining medium was diluted to obtain approximately 10 G418-resistant colonies per plate. Proviral DNA sequence analysis. Genomic DNA from G418-resistant clones and colony pools was prepared as previously described (33). Sequence analysis of individual transduced vector sequences was performed on a PCR product encompassing part of the 59 LTR, the PBS, the 59 UTR, and the upstream part of the neo gene. The PCR was performed with oligonucleotide 1 (ON1), matching Akv MLV positions 7838 to 7865 (64) (59-TTCATAAGGCTTAGCCAGCTAA CTGCAG-39), and ON2, matching neo positions 1656 to 1683 (3) (59-GGCGC
J. VIROL. CCCTGCGCTGACAGCCGGAACAC-39). The resulting PCR product was sequenced by use of an upstream primer (ON3) matching Akv MLV positions 69 to 96 (64) (59-TCCGAATCGTGGTCTCGCTGATCCTTGG-39) and, for relevant clones harboring PBS-Gln, by using a downstream primer (ON4) matching neo positions 1223 to 1244 (2) (59-CTTCCTTTAGCAGCCCTTGCGC-39). PCR-based screening for recombinants. A PCR-based screening for AkvMLEV recombination within the 59 UTR was performed on genomic DNA prepared from colony pools each obtained by pooling of all G418-resistant colonies (approximately 10 colonies per plate) obtained on a single plate. In the PCR amplification, a primer matching the MLEV PBS-Gln (ON5; 59-GTCTTT CATTTGGAGGTCCCA-39) and a neo-specific primer (ON2) were used. The resulting PCR product (if any) for each pool was sequenced with ON4 and ON5. Nucleotide sequence accession number. The MLEV sequence determined in this study has been assigned GenBank accession no. AF041383.
RESULTS 5* UTR sequences of Akv-MLV and MLEV. MLEV is encapsidated into virus particles released from C-2 packaging cells (38). The functional MLEV PBS-Gln may be used for initiation of reverse transcription as an alternative to an impaired vector PBS, and endogenous MLEV RNA may thus serve as a recombination partner in rescue of PBS-modified vectors (38, 39). To further characterize MLEV, we performed, among different PCR-based strategies used, a semirandom-primed PCR on cDNA synthesized from viral RNA from C-2-derived virus particles. Resulting PCR products were sequenced with primers recognizing the PCR product termini, thereby obtaining sequences of the entire 59 UTR and part of the gag coding region. Alignment of the 59 UTRs of Akv and MLEV demonstrated homology throughout the region (Fig. 1). A total of 101 nucleotide position differences were found dispersed between the PBS and the gag start codon; these nucleotide differences were grouped into 34 leader markers (LM1 to LM34) representing single-nucleotide differences (e.g., LM3), clusters of differences (e.g., LM24), and deletions (e.g., LM5) or insertions (e.g., LM20) in MLEV. This genetic marker-based division of the 59 UTR defines an array of SIWs (designated SIWI to SIWXX) ranging in size from 6 to 27 nucleotides (Fig. 1). MLEV harbors the glyco-gag and gag start codons at positions 375 and 639, respectively (Fig. 1). It is uncertain, however, whether functional MLEV glyco-Gag and Gag proteins are produced. Previous studies thus suggest that functional glyco-Gag protein is not encoded by the 59 UTR of an endogenous virus closely related to MLEV (13). Chemical modification studies of Moloney MLV RNA have suggested a highly ordered secondary structure of the 59 UTR (1, 42, 62). The relevant regions of Moloney MLV, Akv MLV, and MLEV were put through a computer-based analysis (performed by use of RNAdraw [37]). On this basis, we predict that similar structure models may account for MLEV and Akv (indicated schematically in Fig. 2A); the only exception is the apparent lack of a stable stem-loop 6 (SL 6) in MLEV. According to the proposed RNA secondary structure model, four potential stem-loops are included in the vectors (pPBSPro 244CAkv-neo and pPBSUmu244CAkv-neo [Fig. 2A]) which we have previously studied in replication and recombination experiments (38). These stem-loops include SL 2, involved in RNA dimerization, and SLs 3 and 4, required for RNA encapsidation (41, 43, 67). The role of SL 1 is unknown, as are the roles of putative SLs 5 to 9. SL 10 is part of an internal ribosome entry site involved in gag polyprotein translation (5, 63). SLs 5 to 9 and a putative shorter form of SL 10 have been included in vectors pPBSPro476CAkv-neo and pPBSU mu476CAkv-neo (Fig. 2A) tested in this study. We conclude that Akv and MLEV 59 UTRs are closely related and likely share a highly ordered secondary structure. The sequences differ only at scattered nucleotide genetic marker positions dispersed throughout the region. Template switch-
VOL. 72, 1998
KISSING-LOOP-MEDIATED MLV RECOMBINATION
6969
FIG. 1. Alignment of Akv and MLEV 59 UTRs. The sequence of MLEV was determined by various PCR-based methods as described previously (38) and in the text. PBS sequences are indicated in bold letters. Identical nucleotide positions are indicated by asterisks; nucleotide insertions in Akv and MLEV are indicated by introduction of colons (:) in MLEV and Akv sequences, respectively. Single-nucleotide differences and clusters of differences are underlined in the MLEV sequence, and the marker number (LM1 to LM34) is given below each genetic marker within the packaging region ranging from the 39 PBS (position 163) to the gag start codon (position 638); LM1 to LM16 correspond to markers IV to XVIII in reference 38. SIWs between markers are indicated by brackets and the designations SIWI to SIWXX. Glyco-gag (positions 375 to 377) and gag (positions 639 to 641) start codons are given in italics. Palindromes (eight nucleotides or longer) in Akv are underlined. The lengths of Akv and MLEV 59 UTR sequences are 476 and 465 bp, respectively; due to a total of five nucleotide insertions in MLEV (indicated by colons in Akv sequence), the length of the entire recombination window is 481 bp.
ing within the 59 UTR of Akv and MLEV templates during minus-strand synthesis therefore represents an in vivo model system with which to evaluate whether homologous recombination between naturally occurring retroviral sequences is affected by the presence of SIWs, palindromic sequences, and putative stem-loop structures. Transfer of vectors with various lengths of sequences between PBS and neo. The versatile function of the MLV 59 UTR suggests that differences in this region may influence both early and late events of retrovirus replication. To estimate whether differences between Akv and MLEV within the 59 UTR affect the production of a protein encoded downstream from the region, the MLEV 59 UTR was introduced into an Akv vector. The resulting vectors harboring PBS-Pro and PBS-Gln were designated pPBSProMLEVCAkv-neo and pPBSGlnMLEVC Akv-neo, respectively (Fig. 2A). neo expression was assessed for vectors harboring PBS-Pro or PBS-Gln by transfection into C-2 packaging cells and subsequent G418 selection. Thus, the number of G418-resistant colonies per transfection indicated whether neo expression was affected by the presence of the 476-bp Akv 59 UTR or the 465-bp MLEV 59 UTR compared, for example, to vectors harboring the shorter 59 UTR. We did not detect any difference for vectors harboring complete Akv and MLEV 59 UTR
sequences (Fig. 2A). A higher level, however, was seen for the vector with the shorter version of the 59 UTR. To obtain an overall estimate of the effect of altering the PBS and the length of the 59 UTR, we measured vector replication efficiencies in transductional titer assays. As expected, transduction of vectors harboring PBS-Umu was strongly diminished, with titer reductions of about 5 orders of magnitude compared to the PBS-Pro constructs (Fig. 2A). There was no significant difference in titer values obtained with pPBSUmu 244C and pPBSUmu476C constructs. For vectors harboring the wild-type PBS sequence, however, we detected about a 10fold increase in titer when the longer 59 UTR was included (Fig. 2A). Since the expression level is not influenced significantly by differences within Akv and MLEV 59 UTRs, recombination between Akv-derived vectors harboring the 476-bp Akv 59 UTR and endogenous MLEV may be studied without a bias for specific recombination sites due to selection at the level of marker gene expression. Furthermore, our transduction data may indicate that one or more cis-acting elements in the downstream part of the 59 UTR are directly involved in retroviral replication or that an overall structure of the entire region is supportive for the actions of the stem-loops in the upstream part of the 59 UTR.
FIG. 2. Vector structure and function. (A) Vectors, expression in packaging cells, and transductional titer. A modified PBS sequence, designated PBS-Umu, was introduced in Akv MLV-derived vectors harboring the neo gene embedded in the bacterial Tn5 transposon. pPBSPro244CAkv-neo and pPBSUmu244CAkv-neo were previously used in studies of retroviral recombination (38). The entire Akv 59 UTR was inserted in vectors pPBSPro476CAkv-neo and pPBSUmu476CAkv-neo. Vectors studied in this work and previously differ in the length of the 59 UTR sequence, being 244 and 476 bp, respectively. neo expression was estimated by transfection of C-2 packaging cells followed by counting of G418-resistant colonies; the resulting estimate of vector expression is given as 102 G418-resistant colonies per transfection of 10 mg of vector DNA. Transductional titers were measured by counting the G418-resistant CFU per milliliter of medium transferred from stably transfected packaging cells; titers have been normalized to 107 producer cells and represent average values for three independent experiments. ND, not determined. The schematic representation of the secondary structure of the MLV 59 UTR is based on studies by Tounekti et al. (62) and Mougel et al. (42). Ten putative stem-loops (SLs 1 to 10) are found within the region. As indicated by dotted lines, SLs 1 to 9 and a putative shortened SL 10 were included in the longer vectors harboring the complete 59 UTR, whereas only SLs 1 to 4 were included in the shorter vectors (38). By analogy with the previously presented model, stem-loops are located as follows (position numbers as given in Fig. 1): SL 1, 236 to 265; SL 2, 291 to 322; SL 3, 329 to 370; SL 4, 372 to 394; SL 5, 400 to 417; SL 6, 423 to 453; SL 7, 459 to 480; SL 8, 486 to 534; SL 9, 539 to 561; and SL 10, 605 to 647. (B) Experimental approach for studies of forced recombination. PBS-mutated vectors were tested in single-cycle transfer protocol by transfection into C-2 packaging cells and subsequent virus transfer to NIH 3T3 target cells. Recombination-based vector transduction, or forced recombination, is the result of (i) heterodimeric RNA encapsidation, (ii) initiation of minus-strand synthesis, (iii) successful plus-strand transfer, and (iv) expression of the neo gene. G418-resistant colonies were individually cloned or pooled (see Materials and Methods) in order to sequence individual transduced proviruses or to allow PCR screening for recombinants, respectively.
6970
VOL. 72, 1998
KISSING-LOOP-MEDIATED MLV RECOMBINATION
6971
TABLE 1. neo transduction by PBS-Umu vectors harboring the Akv 476-bp 59 UTR Transduction pathway Transduced PBS
No. of clones with specific PBS/38 proviruses analyzed
59 UTR minus-strand recombination with indicated recombination partnerb MLEV
PBS-Gln PBS-Pro
2 2
PBS-Umu
33
DPBS
2/2
Mo-C2
2/2
1
R-U5-mediated second-strand transfer recombinationa Recombination partner MLEV
4/33
Minus-strand DNA transferc
Mo-C2
Inter
Intra
4/33
2/2 2/2 2/4 4/4
0/2 0/2 2/4 0/4
Unknownd
25/33 1/1e
a Transduction was mediated through initial priming on MLEV RNA or packaging construct RNA (40), read-through of PBS-Umu during minus-strand synthesis, and subsequent second-strand transfer mediated by R-U5 complementarity of the plus strand and the extended minus-strand DNA acceptor template (39). b Transduced PBS-Gln originates from recombination with MLEV; transduced proviral sequences 28 and 42 are given in Fig. 3. Transduced PBS-Pro originates from recombination with the Moloney MLV-derived C-2 packaging construct (Mo-C2), which harbors a 350-bp deletion within the packaging signal (34). Sequence identities were 43 of 48 positions upstream from deletion and 28 of 52 downstream from the deletion. Recombination sites were mapped within the upstream region (40). c The character of minus-strand DNA strong-stop transfer to the vector 39 end (inter) or MLEV 39 end (intra) was determined through sequence analysis of the U3 region in the downstream LTR. The presence of MLEV U3 sequences thus indicated intrastrand transfer. d Proviruses did not contain sequences of MLEV or packaging construct origin in the regions analyzed (59 LTR, 59 UTR, and 39 LTR). e PBS potentially deleted due to aberrant second-strand transfer (53).
Recombinational repair of PBS-impaired vectors harboring 476-bp 5* UTR. Previous results with PBS-impaired vectors harboring a shortened 59 UTR indicate that the majority of the rare events of transduction of these vectors are mediated by recombination with MLEV RNA, thereby introducing PBSGln in the transduced provirus (38). In the present study of PBS-impaired vectors with the longer 59 UTR, we initially tested by sequence analysis the PBS composition in 38 G418resistant colonies, each representing an individual transductional event (Table 1). Surprisingly, we found that only two clones harbored PBS-Gln, suggesting a relatively low incidence of 59 UTR minus-strand recombination with MLEV. The sequences of the proviral 59 UTR in these two clones (28 and 42) are shown in Fig. 3. The remaining proviruses were results of 59 UTR minus-strand recombination with the Moloney MLVbased packaging construct (34), R-U5-mediated second-strand transfer recombination involving MLEV (39) or the packaging construct (40), or transfer through yet unknown transduction pathways. None of the 25 proviruses that were transduced through unknown pathways harbored sequences of MLEV or packaging construct origin in upstream or downstream LTRs or in the 59 UTR, indicating that these proviruses were generated from vector RNA homodimers by aberrant reverse transcription mechanisms. Most likely, the capacity to form Akv-MLEV heterodimers allowing for recombination is disturbed by the longer 59 UTR. This observation is consistent with the higher titers observed for vectors harboring the 476-bp 59 UTR and possibly reflects that vector RNAs with the entire UTR are more likely to generate packageable homodimers. Alternatively, pPBSUmu 476C vectors are generated at a higher level than vectors with a 244-bp 59 UTR in G418-selected packaging cells and therefore may diminish heterodimerization with MLEV, resulting in a lower incidence of vector-MLEV recombination. PCR screening for recombination with MLEV. To reveal additional template switching events within the 59 UTR, a large number of colonies were screened for Akv-MLEV recombinant proviruses harboring PBS-Gln. Twenty-five colony pools were generated by pooling of G418-resistant NIH 3T3 colonies obtained on 25 separate plates, each containing approximately 10 G418-resistant colonies. Colony pools were screened by PCR with primers matching specifically PBS-Gln
and sequences within the neo gene. PCR products were obtained in 21 of 25 amplifications, and the resulting DNA fragments were sequenced with PBS-Gln and neo primers (Fig. 3). This strategy presented the possibility that two or more proviruses could have been simultaneously amplified in the same PCR. However, in only one case (pool O), two distinct PBSGln-containing proviral sequences were evident in the product amplified from the same pool. This resulted in the generation of an A/G mixture position at LM9. Obviously, we cannot exclude the possibility that two identical proviruses were amplified from the same pool. Sequence analysis of PBS-Gln-containing proviruses revealed the specific MLEV pattern of nucleotide differences from Akv flanked downstream by Akv sequences (Fig. 3). Sequences ranged from harboring only the PBS-Gln of MLEV origin (pools S and V) to harboring markers LM1 to LM28 of MLEV origin (pools F and J). We conclude that PBS-impaired vectors with a complete 59 UTR undergo recombinational repair through initiation of reverse transcription on the functional MLEV PBS and subsequent minus-strand template switching within the 59 UTR to obtain the PBS complementarity required for efficient plus-strand transfer (Fig. 4). The genetic markers dispersed throughout the region allowed us to map specifically the site of template switching in each transduction event. Recombination sites were clustered within markers LM8 and LM10. Thus, we mapped 5 template shifts within SIWVIII and 11 within SIWIX, whereas template switching occurred within SIWI, SIWIII, and SIWXIX in two, four, and two cases, respectively (Fig. 4). We did not in any case register sequence abnormalities such as deletions, insertions, or single-base mutations in the proviral recombinants analyzed (Fig. 3). Misincorporations were thus not evident at the point of transfer of the growing DNA strand, and we propose that in vivo reverse transcriptase-mediated recombination between naturally occurring retroviral sequences is not an error-prone process. In conclusion, precise recombination sites were mapped and found clustered within specific regions of the Akv-MLEV 59 UTR. In addition, we did not in any proviral recombinant find evidence for three template shifts within the 59 UTR. Recombination within Akv-MLEV SIWs. Having identified the exact recombination sites in 24 distinct events involving
FIG. 3. Nucleotide sequences of PBS-Gln-harboring transduced proviruses and recombination partners Akv and MLEV. The sequences of individual transduced proviruses were determined by sequence analysis of PCR fragments encompassing the entire 59 UTR (here defined as the region from PBS to gag start codon). Transduced viral sequences are compared with homologous regions of Akv (top) and MLEV (bottom). Two sequences, 28 and 42, originate from analysis of individual colony clones (Table 1); the remaining sequences originate from PCR screening of colony pools obtained from separate plates and subsequent sequence analysis. Nucleotides homologous to positions in Akv are indicated by hyphens; deleted nucleotides compared to Akv are indicated by colons, whereas insertions are indicated by introduction of a colon in the Akv sequence. Genetic markers consisting of more than one-nucleotide differences are underlined. Molecular differences between Akv-neo and MLEV within the 59 UTR are designated LM1 through LM34; marker numbering is indicated below the Akv sequence. LM1 to LM16 correspond to markers IV to XVIII in reference 38. The gag start codon is indicated for convenience (position 639); however, the ATG sequence was not included in the vectors utilized. R, A/G mixed nucleotide position.
6972
VOL. 72, 1998
KISSING-LOOP-MEDIATED MLV RECOMBINATION
6973
FIG. 4. Mapping of sites for recombinational repair of PBS-modified vectors. Template switching during minus-strand DNA synthesis within the 59 UTR between the PBS and neo gene is established to obtain perfect PBS complementarity in second-strand transfer. Thin lines indicate RNA; thick lines represent DNA. The sloping line represents the 59 UTR with genetic markers LM1 to LM34. Stippled dots indicate single-nucleotide differences, whereas clusters of differences are indicated by black dots. Lengths of SIWs are given between genetic markers. Recombination sites are dispersed throughout the region as delineated by arrow width and ratios (number of proviruses with specific recombination site/total number of proviruses analyzed). In 16 of 24 analyzed proviral sequences harboring PBS-Gln, the recombination site was mapped within a region coinciding with the kissing stem-loop dimerization domain. The kissing stem-loop (SL 2 [Fig. 2A]) harbors LM8 and LM9, whereas LM10 is located within SL 3 (Fig. 2A).
MLEV, we sought to deduce whether specific sites of template switching were preferred by the reverse transcriptase. Among the 20 SIWs (defined as DNA stretches with six or more identical nucleotides) which appeared from the alignment of Akv and MLEV 59 UTRs (Fig. 1), we found recombination sites clustered within 5 (Table 2). This pattern of recombination sites did not reflect a selective bias at the level of transcription or translation of the proviral neo gene, as indicated by the similar levels of marker gene expression obtained with vectors harboring the Akv and MLEV 59 UTRs (Fig. 2A). We thus posit that downstream genetic markers of MLEV origin do not interfere with marker gene expression, and therefore selection at the level of gene expression does not result in discrimination of any sites of template switching. We conclude that the recombination sites detected were not randomly distributed throughout the 59 UTR but that some windows were more prone to template switching than others.
Recombination-prone windows ranged in size from 11 to 24 nucleotides (Table 2); thus, recombination was not detected within windows smaller than 11 nucleotides. In case of random site distribution, we would expect to observe template shifting also within the smaller SIWs. The possibility exists that a certain length of sequence homology or sequence identity is required for strand transfer. However, the fact that large SIWs (such as SIWII, SIWX, and SIWXVIII [Table 2]) did not harbor any recombination site suggested that donor-acceptor homology is not the primary determinant in template switching. As delineated in Table 2, there is no direct correlation between length of SIW and number of sites. Considering the five largest SIWs defined by allowing one single-nucleotide difference (Table 2), we likewise deduce that recombination is not a simple matter of sequence similarity for the reverse transcriptase to shift template. As a striking example, two neighboring windows (SIWVIII-IX and SIWX- XI) of similar lengths (37 and 36 nu-
6974
MIKKELSEN ET AL.
J. VIROL. TABLE 2. Recombination within Akv-MLEV sequence identity windows
Length (nta)
Location
Single I II III IV V VI VII VIII IX X XI XII XIII XIV XV XVI XVII XVIII XIX XX
11 19 21 6 16 14 12 12 24 27 8 12 11 8 7 13 8 23 23 6
PBS-LM1 LM1-LM2 LM2-LM3 LM3-LM4 LM4-LM5 LM6-LM7 LM7-LM8 LM8-LM9 LM9-LM10 LM11-LM12 LM12-LM13 LM13-LM14 LM16-LM17 LM17-LM18 LM22-LM23 LM25-LM26 LM26-LM27 LM28-LM29 LM29-LM30 LM31-LM32
Doubleb I 1 II III 1 IV VIII 1 IX X 1 XI XVIII 1 XIX
31 28 37 36 47
PBS-LM2 LM2-LM4 LM8-LM10 LM11-LM13 LM28-LM30
SIW
a b
Recombination within window (recombinants/total)
Sequential and structural characteristics
2/24 4/24
10-nt palindrome 8-nt palindrome
5/24 11/24
Upstream leg of kissing stem-loop, SL 2 (and 6-nt palindrome) Downstream leg of kissing stem-loop, SL 2 (and 6-nt palindrome) SL 3 Upstream leg of SL 4 Loop, SL 4 (and 10-nt palindrome) Upstream leg of SL 5 (and 6-nt palindrome) Downstream leg of SL 5 Downstream leg of SL 6 Downstream leg of SL 7 Upstream leg of SL 8 Upstream leg of SL 9 Downstream leg of SL 9 (and 6-nt palindrome)
2/24
2/24 4/24 16/24
10-nt palindrome Kissing loop, SL 2 (and 16-nt palindrome) Part of SL 4 SL 9 (and 6-nt palindrome)
2/24
nt, nucleotides. Windows covering two adjacent SIWs separated by one single-nucleotide difference.
cleotides, respectively) appeared to support recombination in 16 (66%) and no transductional events, respectively. However notably, our data do not exclude that sequence homology may be a prerequisite for efficient template switching. DISCUSSION Forced recombination. We have previously seen that PBS knockout Akv MLV vectors may restore function (38) or be rescued (39) through recombination with MLEV, a retroviral sequence endogenous to the murine host genome. In this study, we used impairment of the Akv PBS to study specifically by forced recombination Akv-MLEV heterodimerization and recombination between retroviral species with distinct but homologous 59 UTR sequences. To study template switching events within the leader region which harbors multiple cisacting viral elements, we exploit the dual role of the PBS in initiation of and strand transfer during proviral DNA synthesis. Hence, cDNA synthesis initiation on copackaged MLEV RNA and transfer of minus-strand strong-stop DNA to the 39 end of the vector result in subsequent copying of the MLEV PBS-Gln during plus-strand strong-stop DNA synthesis. Consequently, conventional minus-strand synthesis without strand exchange leads to a dead end due to the lack of complementarity between the mutated PBS and PBS-Gln. Second-strand transfer may in this situation be mediated by complementary sequences within R and U5 regions (39). However, another possibility is switching of the nascent minus-strand DNA from the Akv donor to the MLEV acceptor RNA template within the 59 UTR to obtain perfect PBS-Gln complementarity during plusstrand transfer. Therefore, a transduced PBS-Gln sequence is indicative of the occurrence of recombination within the 59
UTR and may, as in this study, be used for screening among a large number of transduction events. Clustering of recombination sites and template shift precision. In this work, we analyzed 24 independent events of template switching within the complete MLV 59 UTR. Nucleotide differences between the recombination partners allowed us to map the positions where the proceeding reverse transcriptase transfers from donor to acceptor template. Recombination sites were distributed nonrandomly throughout the 481-nucleotide 59 UTR recombination window. This apparent clustering of template switching sites was not caused by a selective bias at the level of marker gene expression, as demonstrated by the ability of MLEV 59 UTR-containing vectors to confer G418 resistance upon murine fibroblasts. Also, this observation is supported by the fact that in two proviruses recombination had occurred near the 39 end of the 59 UTR and by previous studies showing that MLEV-related leader sequences do not interfere with expression of downstream genes (13, 21). Recombination was precise in all recombinants analyzed and thus did not result in any nucleotide misincorporations at frequent or less frequent transfer sites. This observation is consistent with our previous observations of exact template switching within the 59 UTR (38) and with observations by Zhang and Temin demonstrating that homologous recombination in vivo is not error prone (68). Discrepancies with studies showing that junction sites harbor misincorporations in vitro (51, 66) may be due to the fact that in vitro studies do not reflect template switching in vivo or that distinct mechanisms for retroviral recombination differ in their ability to confer precise template switching. The process of nascent-strand transfer may therefore not necessarily in itself contribute to retroviral diversity.
VOL. 72, 1998
KISSING-LOOP-MEDIATED MLV RECOMBINATION
6975
Recombination governed by template homology? The nonrandom pattern of recombination potentially resides in a certain degree of homology between the two templates or in primary or RNA secondary structures within the region. Studies by Zhang and Temin (68) indicate that recombination between nonhomologous RNA templates depends on the length of an inserted sequence identity region. However, a detailed look at the entire region (Fig. 1) and its division into SIWs (Table 2) suggests that homology and template similarity in this case are not the primary determinants for retroviral recombination. If this were so, we would expect recombination sites to be more evenly distributed throughout the region due to extensive homology between the two recombination partners. Since the preferred region of recombination within SIWVIII and SIWIX coincides with the region previously found to mediate recombination in shorter vectors (38), we may also conclude that a need for a minimum length of homology between the growing DNA strand and the RNA acceptor template is not the decisive factor in template shifting. If this were the case, we would observe a substantial amount of recombination events in the downstream part of the 59 UTR, likely within SIWXVIII1XIX. Palindromic sequences involved in recombination? Palindromic sequences may in theory represent local RNA-RNA interaction sites possibly contributing to kissing-loop-mediated dimerization. This notion may be supported by the fact that dimeric RNA interacts at multiple sites along the genome (4). The presence of 16- and 10-nucleotide-long palindromic sequences within SIWXVIII1IX (hosting 16 of 24 sites of recombination) and SIWIII (hosting 4 of 24 sites of recombination), respectively, may indicate some significance of palindromic sequences in template switching and possibly dimerization. However, 8- and 10-nucleotide palindromes within SIWVI and SIWXII, respectively, did not in any case promote template switching. Therefore, we cannot from our analysis of recombination sites determine whether template shifts are facilitated by the presence of palindromes in the RNA templates. Site-specific recombination due to RNA secondary structures? Secondary structures may pause the reverse transcriptase during minus-strand synthesis (16); this pausing promotes strand transfer in vitro (6, 17, 66). In our experiments, 5 of 16 recombination events in the kissing stem-loop were mapped within the upstream leg of the stem-loop (SIWVIII [Table 2]). In addition, other secondary structures, such as the well-established double stem-loop structure (SLs 3 and 4 [Fig. 2]) (43, 67) downstream from the dimerization region, did not promote template shifting. In contrast to our previous studies of vectors containing SLs 1 to 4, we can from the present work state that the reverse transcriptase pausing and sequence homology combined do not promote template switching in the region downstream from the stable SL 4. In summary, we do not believe that pausing of the reverse transcriptase due to intramolecular RNA structures contributes significantly to the recombination process. A similar conclusion was drawn in a recent in vitro study of internal strand transfer at the HIV-1 transactivation response region situated in the R region. In this case, the structure-derived pause site did not coincide with the preferred transfer site, and it was suggested that the close interaction between donor and acceptor stem-loops is the driving force for template switching (28). Based on our results and the considerations presented above, we propose that a specific RNA structure is crucial for recombination within the MLV 59 UTR. The preferred region of template shifting thus coincides precisely with the kissing stem-loop demonstrated by in vitro studies to be crucial for MLV RNA dimerization (20, 52, 62). Considerable evidence
derived from both in vitro and in vivo experiments indicates that the corresponding stem-loop in HIV-1 is involved in, but possibly not essential for, dimerization (4, 8, 9, 14, 22, 30, 31, 35, 44, 47–49, 57). This colocalization of the hot spot recombination site and the primary dimerization domain suggests a relevant connection between the processes of dimerization and recombination. Thus, template shifting at this particular site may be favored by the proximity of RNA templates and by direct RNA-RNA interactions in this part of the genome. It is possible, however, that the palindromic character of the region, allowing perfect acceptor template match perhaps combined with pausing of the reverse transcriptase when it encounters and traverses the stable intermolecular RNA structure at this site, may support strand transfer promoted by the close proximity of the interacting RNA molecules. Also, frequent RNA breakage within certain secondary structures may promote forced copy-choice recombination (10). Model for kissing-loop-mediated retroviral recombination. In summary, we favor the idea that intimate RNA-RNA interactions due to kissing-loop-mediated dimerization are the driving force for recombination within the 481-bp 59 UTR recombination window. Dissociation of the reverse transcriptase from the donor template within this region may thus lead to a shift of template and continued DNA synthesis on the acceptor template. Recent studies suggest that the reverse transcriptase may otherwise tend to reassociate with the donor template RNA (25). Our data also indicate that heterodimerization of Akv and MLEV is mediated by kissing-loop interactions, despite the presence of a single-nucleotide difference within the six-nucleotide loop sequence. Noteworthy, perfect base pairing is restored by G (MLEV):U (Akv) base pairing, and the primary interaction may therefore not be disturbed. As illustrated in Fig. 5, we propose that the reverse transcriptase traversing the RNA dimerization region switches template due to template proximity and possibly reverse transcriptase pausing within or near the putative RNA duplex formed. It should be emphasized that we cannot in the present experimental setup determine the exact frequency of template switching within the particular region and therefore do not know whether the kissing-loop may act as a frequent mediator of recombination also in retroviral replication devoid of selective constraints. This model for kissing-loop-mediated recombination suggests that the kissing-loop interaction, and possibly RNA duplex formation, is intact also within the nucleocapsid core particle during reverse transcription of the mature RNA dimer. This may indicate that linkage in this particular area of the dimeric genome is crucial for the overall structure of the dimer and perhaps strand transfer during viral DNA synthesis, a notion supported by the fact that deletion of the dimerization stem-loop severely inhibits plus-strand DNA transfer (47). Consistent with this idea, it was recently proposed that close RNA interactions possibly within the dimer linkage site of the MLV packaging sequence may promote template switching during reverse transcription of a downstream gene (15). However, the view of a persistent kissing-loop palindrome interaction is challenged by the observation that kissing-loop mutations in HIV-1 do not influence RNA dimer stability, indicating, according to the authors, that the initial kissing-loop interaction is resolved during maturation of the dimer (4). At this stage, we cannot explain these discrepancies. Similarity to RNA structure-based template switching in other RNA viruses. The complementing action of a high mutation rate and frequent genetic recombination is a key player in maintenance of viability and variability in diverse groups of RNA viruses. Evidence suggests that related mechanisms may account for template switching during viral RNA or DNA
6976
MIKKELSEN ET AL.
J. VIROL.
FIG. 5. Model for kissing-loop-mediated recombination (not drawn to scale). The interaction of Akv and MLEV kissing loops and putative subsequent RNA duplex formation promote template switching during viral DNA synthesis, most likely due to close RNA-RNA interactions in the region. Enlargements of the interacting Akv and MLEV stem-loops and the RNA duplex supposedly generated subsequent to loop kissing are shown. It should be noted that the interaction of kissing loops represents a local antiparallel linkage; for convenience, the two interacting RNA templates are presented in antiparallel orientation.
synthesis in diverse viruses. Results obtained from studies of picornavirus (poliovirus) (54, 61), bromovirus (brome mosaic virus) (46), and now mammalian type C retrovirus (MLV) (this study) recombinants indicate that local regions of hybridization
between identical sequences are preferred sites for template switching. In addition, secondary RNA structures have been found to be crucial for recombination in various RNA viral species. Hence, in plant viruses such as brome mosaic virus
VOL. 72, 1998
KISSING-LOOP-MEDIATED MLV RECOMBINATION
(46), turnip crinkle virus (7), and tombusviruses (65) and in animal viruses such as picornaviruses (e.g., human poliovirus) (54) and coronavirus (mouse hepatitis virus) (55), stem-loop structures are primary mediators of template switching. Our results provide another example of viral recombination mediated by defined structural elements in the genome, thus adding support to the notion that recombination mechanisms crucial for viral evolution are exploited by a diverse spectrum of viral species. Common to all of these studies is that preferred recombination sites coincide with regions of predicted RNA secondary structure. Although the exact mechanism in most cases remains a matter of speculation, current knowledge indicates that stem-loop structures contribute to site-specific recombination by mediating intermolecular template interactions or by promoting polymerase pausing. The results presented in this report shed light on recombination at the primary site of viral RNA interaction and thus raise the possibility that alternative interaction sites in the RNA genome play similar roles in retrovirus replication and evolution. ACKNOWLEDGMENTS This work was supported by the Danish Biotechnology Programme, the Danish Cancer Society, the Danish Natural Sciences and Medical Research Councils, the Karen Elise Jensen Foundation, and contracts Biotech CT95-0100 and Biomed2 CT95-0675 from the European Commission. REFERENCES 1. Alford, R. L., S. Honda, C. B. Lawrence, and J. W. Belmont. 1991. RNA secondary structure analysis of the packaging signal for Moloney murine leukemia virus. Virology 183:611–619. 2. Auerswald, E. A., G. Ludwig, and H. Schaller. 1981. Structural analysis of Tn5. Cold Spring Harb. Symp. Quant. Biol. 45:107–113. 3. Beck, E., G. Ludwig, E. A. Auerswald, B. Reiss, and H. Schaller. 1982. Nucleotide sequence and exact localization of the neomycin phosphotransferase gene from transposon Tn5. Gene 19:327–336. 4. Berkhout, B., and J. L. B. van Wamel. 1996. Role of the DIS hairpin in replication of human immunodeficiency virus type 1. J. Virol. 70:6723–6732. 5. Berlioz, C., and J.-L. Darlix. An internal ribosomal entry mechanism promotes translation of murine leukemia virus gag polyprotein precursors. J. Virol. 69:2214–2222. 6. Buiser, R. G., R. A. Bambara, and P. J. Fay. 1993. Pausing by retroviral DNA polymerases promotes strand transfer from internal regions of RNA donor templates to homopolymeric acceptor templates. Biochim. Biophys. Acta 1216:20–30. 7. Cascone, P. J., T. F. Haydar, and A. E. Simon. 1993. Sequences and structures required for recombination between virus-associated RNAs. Science 260:801–805. 8. Clever, J. L., and T. G. Parslow. 1997. Mutant human immunodeficiency virus type 1 genomes with defects in RNA dimerization and encapsidation. J. Virol. 71:3407–3414. 9. Clever, J. L., M. L. Wong, and T. G. Parslow. 1996. Requirements for kissing-loop-mediated dimerization of human immunodeficiency virus RNA. J. Virol. 70:5902–5908. 10. Coffin, J. M. 1979. Structure, replication, and recombination of retrovirus genomes: some unifying hypotheses. J. Gen. Virol. 42:1–26. 11. Coffin, J. M. 1996. Retroviridae: the viruses and their replication, p. 1767– 1847. In B. N. Fields, D. M. Knipe, and P. M. Howley (ed.), Fields virology, 3rd ed. Lippincott-Raven Publishers, Philadelphia, Pa. 12. Colicelli, J., and S. P. Goff. 1987. Identification of endogenous retroviral sequences as potential donors for recombinational repair of mutant retroviruses: positions of crossover points. Virology 160:518–522. 13. Colicelli, J., and S. P. Goff. 1987. Isolation of a recombinant murine leukemia virus utilizing a new primer tRNA. J. Virol. 57:37–45. 14. Darlix, J.-L., C. Gabus, M.-T. Nugeyre, F. Clavel, and F. Barre´-Sinoussi. 1990. Cis elements and trans-acting factors involved in the RNA dimerization of the human immunodeficiency virus HIV-1. J. Mol. Biol. 216:689–699. 15. Delviks, K. A., W.-S. Hu, and V. K. Pathak. 1997. C2 vectors: murine leukemia virus-based self-inactivating and self-activating retroviral vectors. J. Virol. 71:6218–6224. 16. DeStefano, J. J., R. G. Buiser, L. M. Mallaber, P. J. Fay, and R. A. Bambara. 1992. Parameters that influence synthesis and site-specific termination by human immunodeficiency virus reverse transcriptase on RNA and DNA templates. Biochim. Biophys. Acta 1131:270–280. 17. DeStefano, J. J., L. M. Mallaber, L. Rodriguez-Rodriguez, P. J. Fay, and
18.
19.
20. 21. 22. 23. 24. 25. 26. 27. 28. 29. 30. 31.
32. 33. 34. 35.
36. 37. 38.
39.
40. 41. 42.
43.
6977
R. A. Bambara. 1992. Requirements for strand transfer between internal regions of heteropolymer templates by human immunodeficiency virus reverse transcriptase. J. Virol. 66:6370–6378. DiFronzo, N. L., and C. A. Holland. 1993. A direct demonstration of recombination between an injected virus and endogenous viral sequences, resulting in the generation of mink-cell focus-inducing viruses in AKR mice. J. Virol. 67:3763–3770. Fosse´, P., N. Motte´, A. Roumier, C. Gabus, D. Muriaux, J.-L. Darlix, and J. Paoletti. 1996. A short autocomplementary sequence plays an essential role in avian sarcoma-leukosis virus RNA dimerization. Biochemistry 35:16601– 16609. Girard, P.-M., B. Bonnet-Mathonie`re, D. Muriaux, and J. Paoletti. 1995. A short autocomplementary sequence in the 59 leader region is responsible for dimerization of MoMuLV genomic RNA. Biochemistry 34:9785–9794. Grez, M., E. Akgu ¨n, F. Hilberg, and W. Ostertag. 1990. Embryonic stem cell virus, a recombinant murine retrovirus with expression in embryonic stem cells. Proc. Natl. Acad. Sci. USA 87:9202–9206. Haddrick, M., A. L. Lear, A. J. Cann, and S. Heaphy. 1996. Evidence that a kissing loop structure facilitates genomic RNA dimerisation in HIV-1. J. Mol. Biol. 259:58–68. Hu, W.-S., and H. M. Temin. 1990. Genetic consequences of packaging two RNA genomes in one retroviral particle: pseudoploidy and high rate of genetic recombination. Proc. Natl. Acad. Sci. USA 87:1556–1560. Hu, W.-S., and H. M. Temin. 1990. Retroviral recombination and reverse transcription. Science 250:1227–1233. Hu, W.-S., E. H. Bowman, K. A. Delviks, and V. K. Pathak. 1997. Homologous recombination occurs in a distinct retroviral subpopulation and exhibits high negative interference. J. Virol. 71:6028–6036. Junghans, R. P., L. R. Boone, and A. M. Skalka. 1982. Retroviral DNA H structures: displacement-assimilation model for recombination. Cell 30:53– 62. Katz, R. A., and A. M. Skalka. 1990. Generation of diversity in retroviruses. Annu. Rev. Genet. 24:409–445. Kim, J. K., C. Palaniappan, W. Wu, P. J. Fay, and R. A. Bambara. 1997. Evidence for a unique mechanism of strand transfer from the transactivation response region of HIV-1. J. Biol. Chem. 272:16769–16777. Klarmann, G. J., C. A. Schauber, and B. D. Preston. 1993. Template-directed pausing of DNA synthesis by HIV-1 reverse transcriptase during polymerization of HIV-1 sequences in vitro. J. Biol. Chem. 268:9793–9802. Laughrea, M., and L. Jette´. 1994. A 19-nucleotide sequence upstream of the 59 major splice donor is part of the dimerization domain of human immunodeficiency virus 1 genomic RNA. Biochemistry 33:13464–13474. Laughrea, M., L. Jette´, J. Mak, L. Kleiman, C. Liang, and M. A. Wainberg. 1997. Mutations in the kissing-loop hairpin of human immunodeficiency virus type 1 reduce viral infectivity as well as genomic RNA packaging and dimerization. J. Virol. 71:3397–3406. Lenz, J., R. Crowther, A. Straceski, and W. Haseltine. 1982. Nucleotide sequence of the Akv env gene. J. Virol. 42:519–529. Lund, A. H., M. Duch, J. Lovmand, P. Jørgensen, and F. S. Pedersen. 1993. Mutated primer binding sites interacting with different tRNAs allow efficient murine leukemia virus replication. J. Virol. 67:7125–7130. Mann, R., R. C. Mulligan, and D. Baltimore. 1983. Construction of a retrovirus packaging mutant and its use to produce helper-free defective retrovirus. Cell 33:153–159. Marquet, R., J.-C. Paillart, E. Skripkin, C. Ehresmann, and B. Ehresmann. 1994. Dimerization of human immunodeficiency virus type 1 RNA involves sequences located upstream of the splice donor site. Nucleic Acids Res. 22: 145–151. Martinelli, S. C., and S. P. Goff. 1990. Rapid reversion of a deletion mutation in Moloney murine leukemia virus by recombination with a closely related endogenous provirus. Virology 174:135–144. Matzura, O., and A. Wennborg. 1996. RNAdraw: an integrated program for RNA secondary structure calculation and analysis under 32-bit Microsoft Windows. CABIOS 12:247–249. Mikkelsen, J. G., A. H. Lund, K. D. Kristensen, M. Duch, M. S. Sørensen, P. Jørgensen, and F. S. Pedersen. 1996. A preferred region for recombinational patch repair in the 59 untranslated region of primer binding site-impaired murine leukemia virus vectors. J. Virol. 70:1439–1447. Mikkelsen, J. G., A. H. Lund, K. Dybkær, M. Duch, and F. S. Pedersen. 1998. Extended minus-strand DNA as template for R-U5-mediated second-strand transfer in recombinational rescue of primer binding site-modified retroviral vectors. J. Virol. 72:2519–2525. Mikkelsen, J. G., A. H. Lund, M. Duch, and F. S. Pedersen. Unpublished results. Mougel, M., and E. Barklis. 1997. A role for two hairpin structures as a core RNA encapsidation signal in murine leukemia virus virions. J. Virol. 71: 8061–8065. Mougel, M., N. Tounekti, J.-L. Darlix, J. Paoletti, B. Ehresmann, and C. Ehresmann. 1993. Conformational analysis of the 59 leader and the gag initiation site of Mo-MuLV RNA and allosteric transitions induced by dimerization. Nucleic Acids Res. 21:4677–4684. Mougel, M., Y., Zhang, and E. Barklis. 1996. cis-active structural motifs
6978
44. 45. 46. 47.
48. 49.
50. 51. 52. 53. 54.
55.
56.
MIKKELSEN ET AL.
involved in specific encapsidation of Moloney murine leukemia virus RNA. J. Virol. 70:5043–5050. Muriaux, D., P.-M. Girard, B. Bonnet-Mathonie´re, and J. Paoletti. 1995. Dimerization of HIV-1Lai RNA at low ionic strength. J. Biol. Chem. 270: 8209–8216. Murphy, J. E., and S. P. Goff. 1994. Forced integration of Moloney murine leukemia virus DNA with a mutant integration site occurs through recombination with VL30 DNA. Virology 204:458–461. Nagy, P. D., and J. J. Bujarski. 1993. Targeting the site of RNA-RNA recombination in brome mosaic virus with antisense sequences. Proc. Natl. Acad. Sci. USA 90:6390–6394. Paillart, J.-C., L. Berthoux, M. Ottman, J.-L. Darlix, R. Marquet, B. Ehresmann, and C. Ehresmann. 1996. A dual role of the putative dimerization initiation site of human immunodeficiency virus type 1 in genomic RNA packaging and proviral DNA synthesis. J. Virol. 70:8348–8354. Paillart, J.-C., E. Skripkin, B. Ehresmann, C. Ehresmann, and R. Marquet. 1996. A loop-loop “kissing” complex is the essential part of the dimer linkage of genomic HIV-1 RNA. Proc. Natl. Acad. Sci. USA 93:5572–5577. Paillart, J.-C., R. Marquet, E. Skripkin, B. Ehresmann, and C. Ehresmann. 1994. Mutational analysis of the bipartite dimer linkage structure of human immunodeficiency virus type 1 genomic RNA. J. Biol. Chem. 269:27486– 27493. Paillart, J.-C., E. Westhof, C. Ehresmann, B. Ehresmann, and R. Marquet. 1997. Non-canonical interactions in a kissing-loop complex: the dimerization initiation site of HIV-1 genomic RNA. J. Mol. Biol. 270:36–49. Peliska, J. A., and S. J. Benkovic. 1992. Mechanism of DNA strand transfer reactions catalyzed by HIV-1 reverse transcriptase. Science 258:1112–1118. Prats, A.-C., C. Roy, P. Wang, M. Erard, V. Housset, C. Gabus, C. Paoletti, and J.-L. Darlix. 1990. cis elements and trans-acting factors involved in dimer formation of murine leukemia virus RNA. J. Virol. 64:774–783. Pulsinelli, G. A., and H. M. Temin. 1991. Characterization of large deletions occurring during a single round of retrovirus vector replication: novel deletion mechanism involving errors in strand transfer. J. Virol. 65:4786–4797. Romanova, L. I., V. M. Blinov, E. A. Tolskaya, E. G. Viktorova, M. S. Kolesnikova, E. A. Guseva, and V. I. Agol. 1986. The primary structure of crossover regions of intertypic poliovirus recombinants: a model of recombination between RNA genomes. Virology 155:202–213. Rowe, C. L., J. O. Fleming, M. J. Nathan, J.-Y. Sgro, A. C. Palmenberg, and S. C. Baker. 1997. Generation of coronavirus spike deletion variants by high-frequency recombination at regions of predicted RNA secondary structure. J. Virol. 71:6183–6190. Schwartzberg, P., J. Colicelli, and S. P. Goff. 1985. Recombination between
J. VIROL.
57. 58. 59. 60. 61.
62.
63.
64.
65. 66.
67. 68.
a defective retrovirus and homologous sequences in host DNA: reversion by patch repair. J. Virol. 53:719–726. Skripkin, E., J.-C. Paillart, R. Marquet, B. Ehresmann, and C. Ehresmann. 1994. Identification of the primary site of the human immunodeficiency virus type 1 RNA dimerization in vitro. Proc. Natl. Acad. Sci. USA 91:4945–4949. Sørensen, A. B., M. Duch, P. Jørgensen, and F. S. Pedersen. 1993. Amplification and sequence analysis of DNA flanking integrated proviruses by a simple two-step polymerase chain reaction method. J. Virol. 67:7118–7124. Stuhlmann, H., and P. Berg. 1992. Homologous recombination of copackaged retrovirus RNAs during reverse transcription. J. Virol. 66:2378–2388. Temin, H. M. 1993. Retrovirus variation and reverse transcription: abnormal strand transfers result in retrovirus genetic variation. Proc. Natl. Acad. Sci. USA 90:6900–6903. Tolskaya, E. A., L. I. Romanova, V. M. Blinov, E. G. Viktorova, A. N. Sinyakov, M. S. Kolesnikova, and V. I. Agol. 1987. Studies on the recombination between RNA genomes of poliovirus: the primary structure and nonrandom distribution of crossover regions in the genomes of intertypic poliovirus recombinants. Virology 161:54–61. Tounekti, N., M. Mougel, C. Roy, R. Marquet, J.-L. Darlix, J. Paoletti, B. Ehresmann, and C. Ehresmann. 1992. Effect of dimerization on the conformation of the encapsidation psi domain of Moloney murine leukemia virus RNA. J. Mol. Biol. 223:205–220. Vagner, S., A. Waysbort, M. Marenda, M.-C. Gensac, F. Amalric, and A.-C. Prats. 1995. Alternative translation initiation of the Moloney murine leukemia virus mRNA controlled by internal ribosome entry involving the p57/ PTB splicing factor. J. Biol. Chem. 270:20376–20383. Van Beveren, C., J. Coffin, and S. Hughes. 1985. Nucleotide sequences complemented with functional and structural analysis, p. 790–805. In R. Weiss, N. Teich, H. Varmus, and J. Coffin (ed.), RNA tumor viruses, vol. 2. Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y. White, K. A., and T. J. Morris. 1995. RNA determinants of junction site selection in RNA virus recombinants and defective interfering RNAs. RNA 1:1029–1040. Wu, W., B. M. Blumberg, P. J. Fay, and R. A. Bambara. 1995. Strand transfer mediated by human immunodeficiency virus reverse transcriptase in vitro is promoted by pausing and results in misincorporation. J. Biol. Chem. 270: 325–332. Yang, S., and H. M. Temin. 1994. A double hairpin structure is necessary for the efficient encapsidation of spleen necrosis virus retroviral RNA. EMBO J. 13:713–726. Zhang, J., and H. M. Temin. 1994. Retrovirus recombination depends on the length of sequence identity and is not error prone. J. Virol. 68:2409–2414.