![]()
|
|
||||||||
J. Biol. Chem., Vol. 279, Issue 35, 36625-36632, August 27, 2004
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


¶

**
From the
Unité de Régulation Enzymatique des Activités Cellulaires CNRS-URA 2185 and ||Groupe de Virologie Moléculaire et Vectorologie, Institut Pasteur, 2528 rue du Dr. Roux, 75724 Paris, France
Received for publication, May 17, 2004 , and in revised form, June 23, 2004.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Dissecting the mechanisms underlying copy choice in vivo will contribute to the understanding of how recombination drives HIV-1 genome evolution. Several efforts have been made in this sense during the last decade, and various mechanisms have been proposed, based either on the infection of cells in culture (ex vivo systems) (7, 8) or on the reconstitution of the process of reverse transcription with purified proteins and nucleic acids (in vitro systems) (912). Some hard facts have been jointly established by these two approaches, including the enhancement of template switching observed by decreasing the rate of DNA synthesis (13, 14) and the importance of a temporal coupling of RT-encoded polymerase and RNaseH activities (8, 15). However, detailed mechanistic models, mostly proposed solely on the basis of in vitro studies, still await an evaluation of their physiological relevance.
In HIV-1, template switching appears as a frequent process. Indeed, ex vivo experiments in HeLa CD4+-infected cells estimated that three events of template switching occur, on average, per replication cycle (16). Furthermore, it has recently been shown that the recombination rates are influenced by the type of cell infected, with the highest rates observed in macrophages (17). A crucial issue to address now is whether this high rate results from the presence of recombination hot spots interspersed among sequences yielding low recombination rates, or if it reflects a nearly constant frequency of strand transfer along the whole genome. Determining the existence of preferential sites for template switching is important, because they would constitute ideal sites for the study of the mechanisms of recombination, and they would pinpoint blocks of sequences with relatively independent evolutionary history along the genome. At present, indications on this issue are extremely scarce. An analysis of ten intra-subtype B recombinant clones generated in cell culture highlighted the presence of a few putative hot spots along the genome (18), but no molecular characterization of these hot spots was carried out. On the other hand, in vitro studies have highlighted the existence of local fluctuations in the frequency of strand transfer on different sequences leading to the proposal of various mechanistic sketches (1012, 19, 20). However, no indication concerning the behavior of these sequences in vivo is available.
Using a cell-free system we previously identified a sequence within the C2 region of the envelope glycoprotein gp120 where strand transfer was up to five times more frequent than in the surrounding regions (21). This behavior was correlated to the folding of this portion of RNA into a large hairpin, because destabilizing this hairpin triggered a 4-fold decrease in the rate of template switching (10). Here we have developed an experimental system to study recombination in HIV-1 after a single cycle of infection of cells in culture and used it to study recombination in the C2 region.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
R8.2 (23), coding for HIV-1 gag, pol, and accessory proteins, and pHCMV-G (24), which carries the VSV envelope protein gene.
|
Vector Particles Production and Transductions of Cells in Culture HIV-1-based vectors were produced by transient transfection of 293T cells with the genomic plasmids, an HIV-1 encapsidation plasmid (pCMV
R8.2) (23) and a VSV envelope expression plasmid (pHCMV-G) (24), using the calcium phosphate method. Cells were plated at a density of 3.5 x 106 per 100-mm-diameter dish and transfected 1620 h later. The medium was replaced 8 h after transfection, and the vector supernatants were recovered 36 h later. Non-internalized DNA was removed by treatment of the vector supernatants with DNaseI (1 µg/ml in the presence of 1 µM MgCl2) for 30 min at 37 °C. The amount of p24 present in supernatants was determined by using the HIV-1 p24 enzyme-linked immunosorbent assay kit (PerkinElmer Life Sciences), and, when necessary, vector supernatants were concentrated by using Centricon® YM-50 centrifugal filter devices (Amicon-Millipore) before transduction. MT4 cells were transduced with 200 ng of p24 antigen per 106 cells (an approximate multiplicity of infection of 20) in 35-mm dishes in a 500-µl volume. Two hours post-transduction, the cells were diluted up to a 4-ml volume with supplemented RPMI medium and maintained at 37 °C in a 5% CO2 incubator for 40 h.
Purification, Cloning, and Analysis of HIV-1 Reverse Transcription ProductsReverse transcription products (RTP) were extracted 40 h after transduction by following the method described by Hirt (25): cells were lysed by incubation for 10 min at room temperature in a buffer containing 10 mM Tris-HCl, pH 8.0, 10 mM EDTA, and 0.6% SDS. High molecular weight DNA was removed by precipitation at a high salt concentration (NaCl 1 M) during 1218 h on ice. The lysates were centrifuged at 30,000 rpm for 1 h, and the supernatants were treated with 100 µg/ml RNaseA for 1 h at 37 °C and 100 µg/ml Proteinase K for 3 h at 50 °C. After phenol/chloroform extraction, DNA was ethanolprecipitated and purified using the NucleoSpin® Extract clean-up kit (Macherey-Nagel). The purified double-stranded DNA was digested with DpnI for 2 h at 37 °C prior to PCR amplification (20 cycles) with primers BH and SH (Fig. 1). The amplified product was purified on agarose gel, digested with SacII and BamHI, ligated into an appropriate plasmid vector, and used for transformation of E. coli. Plating on IPTG/X-Gal-containing dishes allowed blue/white screening of recombinant and parental colonies, respectively (10, 21). To determine the regions where strand transfer occurred, 48 recombinant clones were analyzed by restriction mapping for each assay, as described under "Results."
Estimation of the Frequency of RecombinationFor the generation of heterozygous particles, equal amounts of pLac+ and pLac plasmids were used in the transfection step. Because the same promoter drives the expression of both types of genomic RNAs, and because the regions involved in their dimerization are strictly homologous between lac+ and lac RNAs, encapsidation of the two genomic RNA moieties is expected to be random. The presence of equivalent amounts of genomic RNAs was confirmed by slot blot analysis on the viral RNA using specific probes for lac+ and lac RNA (not shown). As currently assumed, the viral population is predicted to be constituted by 25% lac+/+ homozygous vectors, 25% lac/ homozygous vectors, and 50% heterozygous vectors (lac+/), according to the Hardy-Weinberg equation (7). Our experimental strategy allows the cloning of products issued from reverse transcription in lac/ as well as in lac+/ vectors. Assuming that only one double-stranded DNA molecule is produced per each vector particle, one-third of the bacterial colonies will be generated by reverse transcription products issued from lac/ vectors and will have a white phenotype, introducing a bias in the estimation of the frequency of recombination. Therefore, for the calculation of the frequency of recombination in heterozygous particles, the total number of colonies is multiplied by two-thirds. To accurately estimate the frequency of recombination, another factor to take into account is the background among the white colonies derived from cloning of cellular DNA co-purified with the reverse transcription products. These colonies are distinguished from those issued from the cloning of RTP by the size of the cloned insert and by their restriction pattern. Typically, 48 white colonies are analyzed in each assay and a correction factor is established, given by n/48, where n is the number of colonies resulting from cloning of RTP. In all experiments n corresponded to approximately half of the colonies analyzed. The frequency of recombination (F) is therefore given by F = b/{2/3[N(n/48) + b]}, where N and b are the total number of white and blue colonies, respectively. The recombination rates per nucleotide (f) within a given interval (i) is given by f = F(xi/X)/z, where F is as above, xi is the number of colonies analyzed where recombination was identified to have occurred within the interval considered, X is the total number of colonies on which mapping was performed, and z is the size in nucleotides of the interval. To calculate the frequency of recombination in control samples ("homozygous," see "Results"), we corrected only for the occurrence of white colonies derived from the cloning of cellular DNA.
In Vitro Recombination AssaysIn vitro recombination assays were done using the reconstituted system previously developed in our laboratory (21). RNA synthesis was performed as previously described (26). RT purification and activity tests were carried out as described by Canard and colleagues (27). Constructs used for RNA synthesis were generated following standard cloning procedures. Reverse transcription was carried out on the donor RNA (100 mM) in the presence of an equimolar amount of acceptor RNA after annealing an oligonucleotide specifically onto the donor template. For the experiments with NC (55 amino acids), the protein was added at a ratio of 1 molecule of NC for 8 nt of total RNA, and incubated for 10 min at 37 °C. Reverse transcription was started by the addition of HIV-1 RT at a final concentration of 400 nM and carried out for 60 min. Synthesis of the second DNA strand, BamHI and PstI digestion, ligation, and E. coli transformation were carried out as previously described (21).
| RESULTS |
|---|
|
|
|---|
U3, Fig. 1C), and therefore lack a functional promoter to drive transcription of new genomic RNAs, hampering the generation of viral progeny. The RTP were recovered using the Hirt technique (25), which allows isolation of low molecular weight DNA. Because our genomic RNAs were devoid of FLAP sequence, which enhances nuclear import and integration (
FLAP, Fig. 1B) (22), most RTPs would be present in the unintegrated low molecular weight fraction. The RTP were analyzed after PCR amplification using primers SH and BH, which allowed the simultaneous amplification of parental and recombinant RTP (Fig. 1C). Prior to cloning in E. coli, the amplified products were digested with BamHI, whose target sequence was present only on lac RNAs, and with SacII, whose recognition site was carried by the SH primer (Fig. 1C). This procedure allows cloning only of parental lac products and of recombinant lac+ products shown in Fig. 1C. The number of lac+ bacterial colonies over the total number of colonies leads to an estimate of the frequency of recombination, as detailed under "Experimental Procedures." For each experiment a control sample was run where homozygous lac+/+ and lac/ vectors were produced separately by transfection of 293T cells with transcomplementation plasmids and either pLac or pLac+ genomic plasmids ("homozygous sample," Table I). In this sample only parental RTP can be generated, and, after restriction with BamHI and SacII, only white bacterial colonies should be found (Fig. 1C). The frequency of blue colonies recovered provides an estimate of the background of artifactual recombinant molecules generated during the experimental procedure.
|
Recombination in the C2 Region of gp120 We first produced vector particles where the region of homology between the two genomic RNAs was constituted by the 400-nt sequence depicted in Fig. 2A. This sequence spans nt 6,6397,039 numbered according to HXB2 proviral DNA and includes the portion coding for the C2 region of the gp120 we previously identified as a recombination hot spot in a cell-free system (10). In three independent experiments the average frequency of recombination was 13.4% (standard deviation: ± 0.4) in "heterozygous" samples, compared with 0.5% (± 0.2) observed in "homozygous" samples (Table I). An extrapolation of this estimate of 0.13 recombination events to the full-length genome (nearly 10,000 nt long in contrast to the 400-nt region used here) predicts 3.25 recombination events per each infection cycle. This value is in agreement with the estimate of three recombination events per infectious cycle made by Jetzt and colleagues (16).
|
To map the positions of strand transfer along the region of homology, point mutations generating specific restriction sites were introduced in the region of homology either on the donor or on the acceptor RNA, defining sub-regions named R1 to R5 following the sense of minus DNA strand synthesis during reverse transcription (Fig. 2A). This allows determining a recombination rate per nucleotide for each individual region after restriction analysis of the cloned RTP, as described under "Experimental Procedures." Regions R1, R2, and R3 participate in the formation of the hairpin previously identified by in vitro probing (Fig. 2B). Three independent experiments yielded consistent results (Fig. 2C) with the upper part of the C2 hairpin (R2) standing out with respect to the other regions. The recombination rate per nucleotide in R2 was 13.3 x 104 (± 2.6 x 104), a value significantly higher than those found in the other regions (chi-squared test: p
0.001). In contrast, when the few blue colonies obtained in the homozygous control samples were analyzed, the distribution of breakpoints appeared random (Table II).
|
0.95), indicating that the presence of the hairpin is important for recombination also in infectious particles. To test this hypothesis we replaced the 25 nt of R3 with the 25 nt colored blue in Fig. 3B, which are expected to regenerate a structure similar to the C2 hairpin (SL RNA, for "Stem Loop," Fig. 3B). On SL RNA, the frequency of recombination in R2 was again significantly higher than in all the other regions (Fig. 3C, chi-squared test: p
0.001), strongly suggesting that the restoration of the hairpin re-established the presence of the hot spot. The predicted stability of the hairpin in SL is significantly higher than in C2 (
G = 46.1 kcal/mol for SL, and
G = 25.0 kcal/mol for C2). To address the question of the influence of the stability of this hairpin on recombination, we constructed 2b RNA (for "2 bulges") by introducing in SL the mutations shown in orange in Fig. 3B. These base substitutions are expected to generate a hairpin with a lower stability than SL (
G = 28.2 kcal/mol) and similar to the one of C2. Furthermore, the shape of the 2b hairpin better resembles the one for C2 with the presence of bulges in the lower part of the stem. With this new RNA, recombination in R2 almost doubled with respect to SL RNA, with more than 50% of the total recombinant molecules mapping in R2 (Fig. 3C, chi-squared test: p
0.001).
|
0.05). When the SL sequence was used, instead, a manifest inhibition was observed at the level of R2 (Fig. 4B, deep blue). To investigate whether the low rate observed in this region resulted from the high stability of the SL hairpin, we ran the assay in the presence of the RNA chaperone NC protein, known to favor breathing of highly stable RNA hairpins (29, 30). The global frequency of strand transfer was increased in this case, with the most pronounced enhancement observed for R2 (Fig. 4B, yellow), although this region still clearly did not constitute a hot spot (chi-squared test: p
0.95). Similarly, with 2b RNA, where the stability of the hairpin was decreased by means of the introduction of point mutations, a manifest enhancement was observed in R2 (Fig. 4B, red), although the recombination rate in R2 was not significantly higher than in the other regions (chi-squared test: p
0.95).
|
|
| DISCUSSION |
|---|
|
|
|---|
The hot region maps in the upper portion of a hairpin structure that we previously identified in vitro (10). By using a series of genomic RNAs where the stability of this hairpin was varied without altering the R2 sequence, we show that the folding of the genomic RNA is crucial for the existence of this hot spot. In fact, in all cases where a large hairpin was expected to be present (C2, SL, and 2b RNAs), the rate of recombination in R2 clearly stood out from that of the other regions (Figs. 2C and 3C). In contrast, on an RNA devoid of a stable hairpin in R2 (Delhp RNA), the rate of transfer in this region fell almost into the background of the other sequences (Fig. 3C). The wild type (C2) and 2b RNAs, which display a similar stability in the hairpin region, were the optimal substrates for recombination.
Limited base substitutions were expected either to reduce or to increase the stability of the 2b hairpin, as the seven bases changed in R1 to generate SL hairpin or the replacement of the R3 region in Delhp RNA both led to a significant decrease in the efficiency of strand transfer specifically in R2 (Fig. 3B). Altogether these results strongly suggest the existence of an optimal window of stability for R2 to constitute a hot spot, as illustrated by the plot presented in Fig. 5. Obviously other factors might also contribute to the high rate of recombination observed in this region such as the presence, location, and size of bulges in the stem portion of the hairpin, the presence of pause sites of reverse transcription, or other features of the R2 sequence itself. These factors might also explain the slight residual preference for recombination in R2 on Delhp RNA (Fig. 3C).
|
The work described here was carried out on sequences from the LAI isolate. The C2 hairpin, however, does not seem to be restricted to this isolate, nor to subtype B, because a comparison of the C2 region between LAI and more than two hundred isolates from all subtypes of the M group indicated that a comparable stem-loop structure is potentially present in all cases (not shown). An implication of the C2 region in the generation of recombinant forms has also been suggested by an analysis of recombination junctions generated after multiple cycles of infection of cells in culture between A/D and B/E subtypes (32). Interestingly, recombination breakpoints in the C2 region were also found in intersubtype recombinant forms isolated from patients, although they seemed to be lost in the endemic circulating recombinant forms (32). Understanding the bases of these differences will certainly be an attractive area of research in the future.
Several in vitro works have suggested an implication of hairpin structures in promoting copy choice by bringing donor and acceptor RNAs into close proximity (33, 34) or by favoring pausing of reverse transcription at the base of the hairpin followed by strand transfer within the hairpin region (20). Based on an in vitro study of the same C2 sequence studied here, we previously proposed that strand transfer would begin by the docking of the nascent DNA onto the acceptor RNA in the loop of the hairpin and would then proceed through a process of strand exchange similar to branch migration occurring during DNA-DNA recombination (10). For this process to be efficient it is predicted that the hairpin must be stable enough as to resist the dynamic changes that the genomic RNA is expected to undergo during reverse transcription, but not too stable as to hamper its opening during the process of branch migration. The indication of the existence of a window of optimal stability for strand transfer to occur ex vivo supports this idea. Along the same lines, the difficulty in melting the stable hairpin may also be reflected by the inhibition of in vitro strand transfer in R2 on SL RNA, because it was relieved by the RNA chaperon NC, known to favor breathing of hairpins. Furthermore, with 2b RNA, where the hairpin is significantly less stable than in SL RNA, the recombination rate in vitro in R2 was ten times higher than in SL, despite the close sequence similarity. The deleterious effect of hairpins being too stable on template switching was also previously demonstrated for terminal strand transfer in vitro on the TAR hairpin (29, 30). Of note, the possibility that the role of the C2 hairpin merely promotes strand transfer at the base of the hairpin by enhancing pausing of reverse transcription is not supported by our results, because the hot region (R2) is located at least 40 nt downstream of the base of the C2 hairpin, following the sense of reverse transcription. Furthermore, in that case, one would expect that increasing the stability of the hairpin would enhance the efficiency of strand transfer in R1, a prediction disproved by the data given in Fig. 5.
In conclusion, this work shows that recombination in a portion of the env gene of HIV-1 does not occur at constant rates but rather presents a clear hot spot region, and underscores the role of the structure of the genomic RNA as an important parameter in the process of copy choice. Determining whether hot spots are frequent along the genome and progressing in their molecular characterization will certainly improve the understanding of the dynamics leading to the generation of the recombinant forms of HIV-1 that challenge immune control in vivo.
| FOOTNOTES |
|---|
Recipient of a postdoctoral fellowship from ANRS. ![]()
¶ Present address: The Wellcome Trust, Tennis Court Rd., Cambridge CB2 1QR, United Kingdom. ![]()
** To whom correspondence should be addressed. Tel.: 33-145-688-505; Fax: 33-145-688-399; E-mail: matteo{at}pasteur.fr.
1 The abbreviations used are: RT, reverse transcriptase; RTP, reverse transcription products; HIV-1, human immunodeficiency virus type 1; VSV, vesicular stomatitis virus; nt, nucleotide(s). ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. Wu, S. L. Heilman-Miller, and J. G. Levin Effects of nucleic acid local structure and magnesium ions on minus-strand transfer mediated by the nucleic acid chaperone activity of HIV-1 nucleocapsid protein Nucleic Acids Res., June 9, 2007; 35(12): 3974 - 3987. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Purohit, B. P. Roques, B. Kim, and R. A. Bambara Mechanisms That Prevent Template Inactivation by HIV-1 Reverse Transcriptase RNase H Cleavages J. Biol. Chem., April 27, 2007; 282(17): 12598 - 12609. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Mild, J. Esbjornsson, E. M. Fenyo, and P. Medstrand Frequent Intrapatient Recombination between Human Immunodeficiency Virus Type 1 R5 and X4 Envelopes: Implications for Coreceptor Switch J. Virol., April 1, 2007; 81(7): 3369 - 3376. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. N. Minin, K. S. Dorman, F. Fang, and M. A. Suchard Phylogenetic Mapping of Recombination Hotspots in Human Immunodeficiency Virus via Spatially Smoothed Change-Point Processes Genetics, April 1, 2007; 175(4): 1773 - 1785. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Fang, J. Ding, V. N. Minin, M. A. Suchard, and K. S. Dorman cBrother: relaxing parental tree assumptions for Bayesian recombination detection Bioinformatics, February 15, 2007; 23(4): 507 - 508. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. A. Baird, R. Galetto, Y. Gao, E. Simon-Loriere, M. Abreha, J. Archer, J. Fan, D. L. Robertson, E. J. Arts, and M. Negroni Sequence determinants of breakpoint location during HIV-1 intersubtype recombination Nucleic Acids Res., October 6, 2006; (2006) gkl669v3. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Chen, D. Powell, and W.-S. Hu High Frequency of Genetic Recombination Is a Common Feature of Primate Lentivirus Replication J. Virol., October 1, 2006; 80(19): 9651 - 9658. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Bampi, A. Bibillo, M. Wendeler, G. Divita, R. J. Gorelick, S. F. J. Le Grice, and J.-L. Darlix Nucleotide Excision Repair and Template-independent Addition by HIV-1 Reverse Transcriptase in the Presence of Nucleocapsid Protein J. Biol. Chem., April 28, 2006; 281(17): 11736 - 11743. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Lanciault and J. J. Champoux Pausing during Reverse Transcription Increases the Rate of Retroviral Recombination J. Virol., March 1, 2006; 80(5): 2483 - 2494. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Galetto, V. Giacomoni, M. Veron, and M. Negroni Dissection of a Circumscribed Recombination Hot Spot in HIV-1 after a Single Infectious Cycle J. Biol. Chem., February 3, 2006; 281(5): 2711 - 2720. [Abstract] [Full Text] [PDF] |
||||