INTRODUCTION

Top Abstract Introduction Materials & Methods Results Discussion References

The dimeric retrovirus RNA genome consists of two plus-strand RNAs which are linked by RNA interactions primarily within the highly structured 5' untranslated region (5' 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 minus-strand 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 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 5' UTR in a region overlapping the retroviral packaging signal. Support for a kissing-loop-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 kissing-loop 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 5' 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 5' 244 nucleotides of the wild-type 476-nucleotide Akv 5' UTR. PBS-modified vectors may thus be transferred through reverse transcriptase-mediated recombination with a MLV-like endogenous virus (MLEV) involving either R-U5-mediated second-strand transfer (39) or minus-strand template shift within the 5' UTR (38). In the latter case, we registered a clustering of recombination sites within a narrow region of the 244-bp 5' 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 5' UTR. We address in this report (i) whether the overall structure of the entire 5' UTR and more specifically cis-acting elements in the downstream part of the 5' 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

Top Abstract Introduction Materials & Methods Results Discussion References

Vector construction. Vectors pPBSPro244Akv-neo and pPBSUmu244Akv-neo have been previously described (33, 38). Vector designations refer to the type of PBS present and 5' UTR length. Vectors harboring the complete 476-bp Akv 5' UTR upstream from the Tn5 transposon fragment which encompasses the neomycin phosphotransferase gene (neo) were generated as follows. A fragment containing the 5' long terminal repeat, (LTR), the proline PBS (PBS-Pro), and the 476-bp Akv 5' 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-Pro244Akv-neo; the resulting retroviral vector was designated pPBSPro476Akv-neo. PBS-Pro was replaced by the PBS-Umu modification by a two-step PCR approach previously described (33) to generate pPBSUmu476Akv-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 5' 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 -2 packaging cells. To obtain the sequence of the 3' part of the 5' 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 5' UTR. The 465-bp MLEV 5' UTR was introduced into vector context by PCR connection with a PCR fragment harboring the 5' LTR and PBS of choice and subsequent cloning of the resulting fragment. The vectors generated were designated pPBSProMLEVAkv-neo and pPBSGlnMLEVAkv-neo (see Fig. 2A).

Cells, transfections, and virus infections. Growth conditions for -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 µg of vector DNA was transfected into -2 packaging cells seeded at 5 × 10³ cells per cm² on the day before transfection. To estimate the level of marker gene expression in vectors with various 5' 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, virus-containing 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 5' LTR, the PBS, the 5' 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) (5'-TTCATAAGGCTTAGCCAGCTAACTGCAG-3'), and ON2, matching neo positions 1656 to 1683 (3) (5'-GGCGCCCCTGCGCTGACAGCCGGAACAC-3'). The resulting PCR product was sequenced by use of an upstream primer (ON3) matching Akv MLV positions 69 to 96 (64) (5'-TCCGAATCGTGGTCTCGCTGATCCTTGG-3') and, for relevant clones harboring PBS-Gln, by using a downstream primer (ON4) matching neo positions 1223 to 1244 (2) (5'-CTTCCTTTAGCAGCCCTTGCGC-3').

PCR-based screening for recombinants. A PCR-based screening for Akv-MLEV recombination within the 5' 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; 5'-GTCTTTCATTTGGAGGTCCCA-3') 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

Top Abstract Introduction Materials & Methods Results Discussion References

5' UTR sequences of Akv-MLV and MLEV. MLEV is encapsidated into virus particles released from -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 -2-derived virus particles. Resulting PCR products were sequenced with primers recognizing the PCR product termini, thereby obtaining sequences of the entire 5' UTR and part of the gag coding region. Alignment of the 5' 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 5' 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 5' UTR of an endogenous virus closely related to MLEV (13).

View larger version (40K): [in this window] [in a new window]   FIG. 1.   Alignment of Akv and MLEV 5' 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 3' 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 5' 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.

View larger version (30K): [in this window] [in a new window]   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. pPBSPro244Akv-neo and pPBSUmu244Akv-neo were previously used in studies of retroviral recombination (38). The entire Akv 5' UTR was inserted in vectors pPBSPro476Akv-neo and pPBSUmu476Akv-neo. Vectors studied in this work and previously differ in the length of the 5' UTR sequence, being 244 and 476 bp, respectively. neo expression was estimated by transfection of -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 µg 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 5' 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 5' 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 -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.

Transfer of vectors with various lengths of sequences between PBS and neo. The versatile function of the MLV 5' 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 5' UTR affect the production of a protein encoded downstream from the region, the MLEV 5' UTR was introduced into an Akv vector. The resulting vectors harboring PBS-Pro and PBS-Gln were designated pPBSProMLEVAkv-neo and pPBSGlnMLEVAkv-neo, respectively (Fig. 2A).

Recombinational repair of PBS-impaired vectors harboring 476-bp 5' UTR. Previous results with PBS-impaired vectors harboring a shortened 5' UTR indicate that the majority of the rare events of transduction of these vectors are mediated by recombination with MLEV RNA, thereby introducing PBS-Gln in the transduced provirus (38). In the present study of PBS-impaired vectors with the longer 5' UTR, we initially tested by sequence analysis the PBS composition in 38 G418-resistant 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 5' UTR minus-strand recombination with MLEV. The sequences of the proviral 5' UTR in these two clones (28 and 42) are shown in Fig. 3. The remaining proviruses were results of 5' UTR minus-strand recombination with the Moloney MLV-based 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 5' UTR, indicating that these proviruses were generated from vector RNA homodimers by aberrant reverse transcription mechanisms.

View larger version (61K): [in this window] [in a new window]   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 5' 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 5' 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.

PCR screening for recombination with MLEV. To reveal additional template switching events within the 5' 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 PBS-Gln-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.

View larger version (21K): [in this window] [in a new window]   FIG. 4.   Mapping of sites for recombinational repair of PBS-modified vectors. Template switching during minus-strand DNA synthesis within the 5' 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 5' 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).

Recombination within Akv-MLEV SIWs. Having identified the exact recombination sites in 24 distinct events involving 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 5' 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 5' 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 5' UTR but that some windows were more prone to template switching than others.

	DISCUSSION

Top Abstract Introduction Materials & Methods Results Discussion References

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 5' UTR sequences. To study template switching events within the leader region which harbors multiple cis-acting 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 3' 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 5' UTR to obtain perfect PBS-Gln complementarity during plus-strand transfer. Therefore, a transduced PBS-Gln sequence is indicative of the occurrence of recombination within the 5' 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 5' 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 5' 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 5' 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 3' end of the 5' UTR and by previous studies showing that MLEV-related leader sequences do not interfere with expression of downstream genes (13, 21).

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 5' UTR, likely within SIWXVIII+XIX.

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 SIWXVIII+IX (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).

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 5' 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.

View larger version (22K): [in this window] [in a new window]   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.

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 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 (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.