Previous Article | Next Article ![]()
Journal of Virology, August 2008, p. 8172-8182, Vol. 82, No. 16
0022-538X/08/$08.00+0 doi:10.1128/JVI.01115-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Department of Medical Microbiology, University of Manitoba, Winnipeg, Manitoba, Canada,1 National HIV and Retrovirology Laboratories, Public Health Agency of Canada, Ottawa, Ontario, Canada,2 Department of Medical Microbiology, University of Nairobi, Nairobi, Kenya,3 Public Health Agency of Canada, Winnipeg, Manitoba, Canada4
Received 27 May 2008/ Accepted 4 June 2008
|
|
|---|
|
|
|---|
The in vivo antiviral mechanism of APOBEC3F/G is not well defined. Certainly, G-to-A hypermutation would cause detrimental mutations in the viral genome, reducing or abolishing the production of viable progeny, but other antiviral APOBEC functions have also been reported. APOBEC-mediated hypermutation may trigger degradation of the nascent viral DNA by host uracil glycosylases and apurinic/apyrimidinic endonucleases, thereby inhibiting the production of provirus (51). APOBEC has also been suggested to interfere with removal of primer tRNA, as well as to inhibit DNA strand transfer and integration (27, 33). Some antiviral effects have been reported for APOBEC in the absence of hypermutation. Guo et al. reported reduced reverse transcription priming and reduced levels of viral DNA in Vif-negative virus in the absence of APOBEC3G deamination and hypermutation (15). APOBEC3G was also found to interfere with proviral integration in the absence of hypermutation (27). However, a recent publication stated that deaminase-defective APOBEC3G had antiviral activity only when it was expressed at high levels, questioning the physiological relevance of deaminase-independent activity of APOBEC3G (35).
The roles of APOBEC3F/G in HIV disease progression and pathogenesis are similarly unclear, although a relationship between the two has been suggested (1, 10, 20, 21, 29, 36). APOBEC3G mRNA expression level has been associated with decreased viremia and long-term nonprogression (20) and with HIV-exposed seronegative individuals (1). However, these studies have not been readily replicated (10), resulting in significant controversy over the role of APOBEC in HIV disease progression in vivo. To examine the phenomenon of APOBEC-mediated hypermutation in a clinical setting, we examined proviral DNA from 240 HIV-infected subjects from Nairobi, Kenya. Significant hypermutation of the predominate provirus was identified and found to correlate with patient CD4 count, suggesting that APOBEC-mediated hypermutation influences disease progression in a clinical setting.
|
|
|---|
DNA isolation and PCR amplification. DNA was isolated from PBMCs using the QIAamp DNA mini kit (Qiagen Inc., Mississauga, Canada). Nested PCRs were performed in order to amplify the proviral vpu and the first 349 nucleotides of env. The first reaction created an amplicon that was nearly 2 kb, using primers outer-f (5'-CAAGCAGGACATAACAAGGTAG-3') and outer-r (5'-TGTTATTTCTAGATCCCCTCCTG-3') and the Expand High Fidelity PCR System (Roche Diagnostics GmbH, Mannheim, Germany) at the recommended conditions. The secondary nested reaction specifically amplified the region of interest, using primers inner-f (5'-GGCTTAGGCATCTCCTATGGCAGGAAGAAG-3') and inner-r (5'-CGAGTGGGGTTAACTTTACACATG-3'), also with the Expand High Fidelity PCR System. Nested PCR was also used to amplify the proviral vif region, as described above, using outer primers 5'-AGTTATCCCAGCAGAAACAGGAC-3' and 5'-TCGCTGTCTCCGCTTCTTCCTG-3' and inner primers 5'-CTGCAGTTAAGGCAGCCTGTTG-3' and 5'-CTTCAACTCCTGCCCAAGTATC-3'.
Sequencing and assembly. PCR amplicons were directly sequenced, using ABI Prism BigDye terminator cycle sequencing ready reaction kit v.3.1 (Applied Biosystems) and the inner PCR primers described above. The sequences were resolved on an ABI Prism 3100 genetic analyzer (Hitachi, Japan). The resulting complementary sequence electropherograms were assembled using Sequencher v.4.2 (Gene Codes Corporation, Ann Arbor, MI). Sequences were aligned using ClustalW (7).
Cloning. Samples that contained gross mutations in vpu and env were cloned into pCR4-TOPO vector (TOPO TA cloning kit for sequencing; Invitrogen Life Technologies, Carlsbad, CA). Thirteen to 20 clones, depending on cloning efficiency, were sequenced for each patient as described above using M13 forward and reverse primers. The PCR amplicon was additionally cloned for a randomly determined subset of samples, and 14 to 20 clones from each sample were sequenced and analyzed to ensure that the PCR sequence was representative of the patient proviral population.
Statistics. The statistical analyses Mann-Whitney test and test for correlation were performed using GraphPad Prism 4.
RNA isolation from plasma. Viral nucleic acids were isolated from archived heparin-plasma using NucliSens isolation reagents (bioMerieux Canada) according to the manufacturer's instructions (4). Reverse transcription-PCR (Qiagen OneStep reverse transcription-PCR kit) was performed using the specific primers outer-f and outer-r, and the viral nucleic acids were further amplified with the nested PCR strategy described above. The plasma RNA sequences were aligned with all the proviral sequences using ClustalW (7). MEGA 3 was employed for neighbor-joining phylogenetic analysis (nucleotide distance calculated by Kimura's two-parameter method) (25). MEGA 3 was employed for bootstrapping to determine whether the two independently derived sequences from the same patient significantly clustered together.
Hypermutation detection. Hypermut 2.0, available from http://www.hiv.lanl.gov/content/sequence/HYPERMUT/hypermut.html, was employed to detect APOBEC-type hypermutation (38). Sequences generated in this study were used to generate a population-specific consensus sequence for each HIV-1 clade that was represented. The sequences were then compared on a clade-by-clade basis to the consensus sequence using Hypermut to identify G-to-A hypermutation and also to identify the nucleotide context of this hypermutation.
Determination of CD4 count. Peripheral blood T-lymphocyte subset analysis was performed as previously described (13).
Determination of plasma viral load. Viral load determinations were performed on 200 µl heparin-plasma using the NucliSens HIV-1 QT assay (bioMérieux, Marcy l'Etoile, France) according to the manufacturer's instructions. The lower detection limit of the assay, when starting from this volume, is 125 copies/ml.
GenBank accession numbers. The genetic sequences for the vpu/env segment from the 240 proviral samples and 18 plasma RNA samples have been submitted to GenBank under accession numbers EU836326 to EU836565 and EU836566 to EU836593. Unique clones from the nonhypermutated sequences and the hypermutated sequences were submitted to GenBank under accession numbers EU875081 to EU875368. The genetic sequences for vif from six proviruses that were hypermutated in the vpu/env regions and eight proviruses that were not hypermutated in the vpu/env regions have been submitted to GenBank under accession numbers EU839404 to EU839417.
|
|
|---|
![]() View larger version (29K): [in a new window] |
FIG. 1. Sequence context of identified G-to-A hypermutation in the first 75 nucleotides of the env ORF from 13 dramatically mutated proviruses. The top panel displays the nucleotide sequence, while the bottom panel shows the corresponding amino acid sequence. The consensus sequence for both the DNA and protein was generated based on the 240 proviral sequences examined in this publication. The sequences are numbered from the HIV-1 env start. The proviral sequence was compared to the corresponding nonhypermutated plasma RNA sequence, where available, or else to a clade-specific consensus sequence. Darkened nucleotide letters in the sample sequences indicate changes from G to A in either a GG or GA context (i.e., proposed APOBEC3G or APOBEC3F hypermutation). Darkened amino acid letters in the sample sequences indicate where proposed hypermutation at the nucleotide level caused an amino acid change at the protein level.
|
![]() View larger version (15K): [in a new window] |
FIG. 2. Proviral adenine proportion. (A) Distribution of proviral adenine proportion in a 590-nucleotide fragment spanning vpu and the 5' end of env for 240 HIV-1 isolates. The adenine proportion is based on proviral sequences from PBMC DNA isolated from HIV-infected women from two Kenyan cohorts. The values greater than 1 standard deviation above the mean ( ± 1 S.D.) and the outliers are indicated. (B) Adenine proportion in proviral HIV-1 sequence and plasma-derived viral RNA sequence in 18 matched samples. The same 590-nucleotide HIV-1 genomic region was examined for both proviral and viral RNA sequences. The four sequences indicated with the asterisk had hypermutated proviral vpu/env sequences.
|
Directly sequenced PCR products show extensive hypermutation and are representative of proviral sequence. Direct sequencing of the PCR amplicon will reveal the major proviral sequence but may not be representative of intrapatient diversity. In order to examine proviral diversity within individual subjects, we sequenced 13 to 20 clones from 23 patients, including all 13 subjects with dramatically hypermutated proviral sequence and a random subset of 10 subjects with proviral sequences that were not dramatically hypermutated. The sequence data generated from the multiple clones for each subject confirmed that the detrimental mutations in Vpu and Env were not the result of sequencing or sampling errors but were in fact representative of the subjects' proviral population. For the majority of hypermutated samples, the clonal sequences showed levels of hypermutation similar to that observed in the sequence generated directly from the proviral PCR product (Fig. 3). Fifteen clonal populations showed hypermutation patterns identical to that of the direct PCR sequence for all the clones examined (including 10 out of the 13 extremely hypermutated sequences). Additionally, two clonal populations had the majority of cloned sequences (85% or greater) match the direct PCR sequence (the hypermutated sequences ML1578 and ML1970). In three clonal populations, half of the clones matched the direct PCR sequence, including the hypermutated sequence ML1857. In three clonal populations, the hypermutation patterns in the clones were generally similar, though not identical, to that observed in the direct PCR sequence. Overall, the vast majority of clonal and proviral sequences from the samples display almost identical patterns of hypermutation; for example, ML103 shows no hypermutation in the direct PCR sequence and only one site of hypermutation in 1 of the 14 clones examined (Fig. 3). Similarly, ML1053 shows many hypermutation sites in the direct PCR sequence, and these are shared with and between all 19 clones examined. However, the levels of hypermutation of 3 out of 16 clonal sequences from ML1578 and 9 out of 17 clonal sequences from ML1857 were lower than the level of hypermutation of the proviral sequence. Conversely, 2 out of 19 clonal sequences from ML1970 had higher levels of hypermutation than did the proviral sequence. These data suggest that the direct PCR sequence is generally representative of the predominant proviral sequence; however, as expected, subtle diversity can be observed in some patients.
![]() ![]() View larger version (186K): [in a new window] |
FIG. 3. Clones and directly sequenced PCR product are compared to a population-specific, clade-specific consensus sequence. The nucleotides are numbered across the bottom of each panel, indicating position within the 590-nucleotide region. ML103 was selected as a representative nonhypermutated sequence. Changes in the patient sequence compared to the consensus sequence are indicated by colored bars, as indicated in the legend. The representations were generated with the Highlighter tool available from http://www.hiv.lanl.gov/content/sequence/HIGHLIGHT/highlighter.html. Black circles below hypermutation sites were added manually to ease visualization. The sequences are ordered from top to bottom from the most similar to the consensus sequence to the least similar. Identical sequences are shown once, with the number of times the sequence was retrieved listed beside the clone name. Patients with an asterisk showed intraclonal diversity.
|
|
View this table: [in a new window] |
TABLE 1. Hypermutation criteria for proviral HIV-1 sequences that were significant in at least one categorya
|
To further examine the roles of these enzymes in the 13 proviral sequences with the most pronounced hypermutation, we examined the dinucleotide sequence connotation of the G-to-A hypermutation. Figure 4 shows that most of the G-to-A nucleotide changes occurred where the original sequence was a GA or GG, suggestive of APOBEC3F or APOBEC3G involvement, respectively. Comparatively fewer changes occurred at GC or GT dinucleotides. It is interesting to note that patients with a higher proportion of G-to-A nucleotide changes occurring at a GG dinucleotide context have a comparatively lower proportion of changes occurring at a GA dinucleotide context. The reverse is also true. In patient ML1970, however, the proportions of G-to-A nucleotide changes occurring at GG and GA dinucleotides are similar.
![]() View larger version (20K): [in a new window] |
FIG. 4. Dinucleotide context of G-to-A hypermutation in 13 patients. The proportion of G-to-A hypermutation in each proviral sequence occurring at the dinucleotides GG, GA, GC, and GT is indicated. Hypermutation context was determined through comparison to a clade-specific consensus sequence, using Hypermut 2.0.
|
![]() View larger version (19K): [in a new window] |
FIG. 5. Association of higher CD4 count with hypermutated sequences. CD4 counts are compared in patients with hypermutated provirus and patients without hypermutated provirus. (A) Hypermutation was defined as significant general APOBEC cytidine deaminase activity, as determined by Hypermut 2.0. (B) Hypermutation was defined as sequences with mutations in Vpu and Env, proviral adenine proportion greater than 1 standard deviation above the mean, and significant Hypermut 2.0 hypermutation. Groups are compared by a Mann-Whitney test. The height of the bar indicates the mean, and the error bar represents the standard error of the mean.
|
![]() View larger version (8K): [in a new window] |
FIG. 6. Correlation of adenine proportion with CD4 count. CD4 counts were available for 208 subjects. Adenine proportion was measured from the 590-nucleotide vpu/env HIV-1 proviral region. These two measures were significantly positively correlated (P = 0.041).
|
of 4.07 log10 copies/ml versus 4.78 log10 copies/ml, respectively), the distribution was not statistically significant (P = 0.82). However, viral load determination was conducted on archived samples, not freshly isolated plasma. Sample storage is known to affect viral load measurement, whereas the CD4 counts were measured on site from fresh blood and are likely more accurate (6). Hypermutation is not associated with Vif mutations. In order to determine the role Vif sequence polymorphisms may be playing in the observed APOBEC3F/3G-like cytosine deamination hypermutation, a subset of Vif sequences was examined from both highly hypermutated and nonhypermutated proviruses. Intact starting methionine residues, no stop substitutions, and a low or undetectable level of hypermutation was found in the vif genes of all samples examined. Amino acid polymorphisms in the Vif sequences were identified, but none in regions that had been previously identified as critical for Vif interaction with APOBEC (28, 32, 34, 39, 46, 50, 54). Additionally, only conservative substitutions were identified at single amino acid residue locations that had been suggested to have an impact on APOBEC interaction (43). Overall, Vif from highly hypermutated provirus was markedly similar to Vif from the other samples in this study and to published consensus sequences, suggesting that Vif polymorphisms are not responsible for the increased APOBEC3-type hypermutation observed within this population.
|
|
|---|
The sequence context of the observed hypermutation further suggested that the host antiviral APOBEC3F and APOBEC3G proteins were involved. The Hypermut 2.0 tool identified significant general APOBEC3-type hypermutation in these sequences, as well as specific APOBEC3F and APOBEC3G hypermutation; APOBEC3F causes GA-to-AA nucleotide changes, while APOBEC3G causes GG-to-AG transitions (26, 52). It is interesting to note that 10 samples showed significant Hypermut 2.0 hypermutation characteristic of both APOBEC3G and APOBEC3F, suggesting that these proteins may be working in concert in these patients. Examination of the dinucleotide context for the G-to-A hypermutation in 13 patients with dramatic hypermutation supported these observations. An elevated GG context of the hypermutation would be expected to be due to APOBEC3G involvement, while an elevated GA context would be caused by APOBEC3F involvement. In contrast, the GC and GT context should be representative of the background mutation rate. In most of the examined sequences, both the GG and GA context were elevated compared to GC and GT, again suggesting coinvolvement of APOBEC3F and APOBEC3G. Indeed, it has been shown that these proteins can form heteromultimers (17, 48, 49).
A subset of subjects that were determined to have dramatic levels of hypermutation (as defined by the presence of premature stop codons, elevated adenine proportion, and a significant APOBEC3F/G sequence connotation to the hypermutation) had significantly higher CD4 cell counts than the other subjects did. These subjects represent an extreme, where severely hypermutated provirus appears to be the predominant proviral species. A previous study of nine patients with long-term viral suppression due to highly active antiretroviral therapy demonstrated that hypermutation could be identified, albeit at a minority, in the proviral sequences from all the subjects (22). Entertaining the likelihood, therefore, that all HIV-positive patients have a variable amount of proviral hypermutation, adenine proportion in the dominant proviral sequence was used to estimate the level of hypermutation in all subjects and was indeed found to be positively correlated with CD4 cell count in all 240 subjects. This indicates that as the dominant proviral HIV-1 sequence is increasingly hypermutated, the subjects' CD4 counts similarly increase, both for patients with dramatic hypermutation and patients with subtle hypermutation. Although correlation does not necessarily indicate causation, this supports the hypothesis that with increasing hypermutation, fewer viable viral progeny are produced, protecting infected cells from direct viral cytopathic effects (22). This work thus also supports the antiviral role of APOBEC-mediated hypermutation, which has been recently debated (15, 27, 33, 35, 51).
An alternative explanation of the association of hypermutation with CD4 count is that a strong host cytotoxic T lymphocyte (CTL) response targets cells that are infected with replication-competent virus for destruction. As a result, cells infected with non-replication-competent HIV-1, such as virus affected by the observed hypermutation, would be comparatively enriched in patients with strong, protective immune responses and therefore higher CD4 cell levels. However, if this were true, a general enrichment of mutation would be expected, due to the error-prone nature of HIV-1 reverse transcriptase, not just G-to-A hypermutation. Yet, all the proviruses that had detrimental mutations in Vpu and Env were hypermutated, and all the mutations were due to nucleotide G-to-A hypermutation. Also, the G-to-A changes in the dinucleotide context of GG and GA, hallmarks of APOBEC3G and APOBEC3F activity, respectively, were elevated above the dinucleotide context of GC and GT, which would represent non-APOBEC-mediated mutations (26, 52). In the absence of APOBEC-specific hypermutation, the levels would presumably be similar. Furthermore, Vif, which is located at a local minimum of APOBEC hypermutation, was not significantly mutated (44). This also would not be expected if another mutating factor had a significant effect on these viruses. Nevertheless, APOBEC may work in cooperation with strong CTL responses to control HIV-1 replication, highlighting the need for continued research about the cooperation between host innate, intrinsic, and adaptive immunity.
The HIV-1 protein Vif is known to counteract APOBEC3F and APOBEC3G, but examination of Vif sequences from subjects did not reveal detrimental mutations that would explain the seemingly increased APOBEC3 hypermutation activity. Furthermore, hypermutation was not observed in the proviral vif sequences, although sequences were obtained from subjects that had significant proviral hypermutation in the vpu/env portion of the HIV-1 genome. This lack of Vif mutation is in contrast to a study by Pace et al. of 127 clade B HIV-1-infected patients from an Australian cohort, which found a correlation between hypermutated proviral sequences and viral load and attributed the relationship to stop mutations due to hypermutation in the infecting viral Vif that disabled this viral defense against cellular APOBEC (36). Priming for DNA synthesis occurs at fixed locations in the HIV-1 genome, and therefore, not all sites are single stranded for the same period of time, suggesting that as APOBEC3G acts only on single-stranded DNA, hypermutation will not occur uniformly across the genome (19, 44). In fact, Vif is located at a local minimum of predicted mutation, suggesting a mechanism that explains our findings (44).
To explain our observations, we hypothesize that subjects infected with more dramatically hypermutated HIV-1 provirus have elevated APOBEC3F/G activity. Potential causes for this elevated activity include upregulation of APOBEC3, such as by gamma interferon, which has been recently shown to be a potent inducer, allowing APOBEC to overcome Vif inhibition (3, 37, 40, 45), or the enzyme may exist in a more enzymatically active form (for example, a higher proportion of APOBEC3G may be present in its low-molecular-mass, active form) (9). Alternatively, the infecting viruses may encode a Vif protein that is inefficient or unable to target APOBEC3F and APOBEC3G for degradation or prevent their packaging into nascent viral particles, although our data do not support this hypothesis (11, 28, 32, 34, 42, 46, 50, 53, 54) In either scenario, we theorize that as a result of increased APOBEC3 or decreased Vif level/activity, the newly budded viral particles contain higher levels of APOBEC3F and APOBEC3G, ultimately resulting in hypermutation during infection of new cells. It seems likely that the resulting provirus is sufficiently hypermutated that replication-competent, hypermutated progeny cannot be generated, as illustrated by the lack of hypermutated RNA sequences. Importantly, the cells infected with hypermutated provirus would not be subjected to HIV's direct cytopathic effects, nor would they be a target of CTLs, likely explaining the observed higher CD4 counts and presumable decrease in disease progression in patients with increased proviral hypermutation (22).
Two conflicting papers highlight that the role of APOBEC3 in HIV disease progression has not yet been resolved. In their study using stimulated PBMCs, Jin et al. (20) found that APOBEC3G mRNA levels were positively correlated with patient CD4 count and inversely correlated with viral load, while Cho et al. (10), using unstimulated PBMCs, found that neither APOBEC3G nor APOBEC3F mRNA correlated with CD4 count and viral load. The hypermutation observed in the present study is likely a direct effect of APOBEC3F/G cytidine deamination, whereas mRNA expression levels, which were examined in these previous studies, may not directly correlate with enzymatic activity due to translational and posttranslational regulation. This direct assessment of sequence hypermutation, likely the result of APOBEC3 cytidine deaminase activity, found that as the level of hypermutation in the subjects' predominant proviral sequence increased, so did the subjects' CD4 counts, suggesting an in vivo role for APOBEC in disease progression.
APOBEC3G/3F proviral hypermutation likely exists in a spectrum. The 13 subjects with dramatically hypermutated provirus had detectable viral loads, intact viral RNA sequence, and significantly depressed CD4 counts due to HIV infection compared to patients not infected with HIV, suggesting that there may be a tissue or cellular compartment where APOBEC3 is not packaged into viral progeny, leading to the production of nonhypermutated virus, which in turn leads to nonhypermutated infectious viral progeny, sustaining the infection. On the other end of the spectrum, patients without hypermutation in the dominant provirus may have minority sequences that are hypermutated (22). Indeed, the examination of clonal sequences indicated some degree of variability within patients. The data presented illustrate a correlative relationship between this spectrum of hypermutation, likely the direct result of APOBEC3F/3G activity and CD4 cell count in vivo in the absence of obvious Vif polymorphism. These findings highlight the potential for enhancing host APOBEC3F and APOBEC3G for therapeutic purposes and suggest that even small increases in APOBEC activity may attenuate HIV-1 replication and disease progression.
A.M.L. is supported by a Canadian Graduate Doctoral Scholarship from CIHR and by the CIHR/ICID training program. F.A.P. is a Tier 1 CIHR Canadian Research Chair.
Published ahead of print on 11 June 2008. ![]()
|
|
|---|
A hypermutation does not result from polymerase chain reaction. AIDS Res. Hum. Retrovir. 13:985-986.[Medline]
A hypermutation in protease and reverse transcriptase regions of human immunodeficiency virus type 1 residing in resting CD4+ T cells in vivo. J. Virol. 79:1975-1980.
A hypermutation. Bioinformatics 16:400-401.
A hypermutation, p. 57-63. In G. Myers, B. T. Korber, B. T. Foley, K.-T. Jeang, J. W. Mellors, and S. Wain-Hobson (ed.), Human retroviruses and AIDS 1996. Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM.This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»