Objective: To identify variations in the env gene of human immunodeficiency virus type 1 (HIV-1) subtype CRF01-AE strains circulating in China and to elucidate the potential relationship between genetic variation and evolutionary pressure.
Methods: Fragments of the HIV-1 env gene were amplified by nested-polymerase chain reaction (n-PCR) from the whole blood of HIV-1 infected individuals from four provinces in Southeast China (Guangdong, Hunan, Jiangsu and Jiangxi). The PCR products were then directly sequenced by ABI 377 DNA sequencers. The sequences covering the env V3-V4 region of 34 HIV-1 subtype CRF01-AE strains were selected to analyse phylogenetic trees and amino acid mutations. The accumulation of synonymous (Ks) and antonymous (Ka) substitutions as well as Ks/Ka ratios were calculated using DIVERGE.
Results: Phylogenetic trees showed that the 34 HIV-1 subtype CRF01-AE strains from China clustered with the Chinese AE reference strain (AE.97CNGX2F), as well as with the reference strains from Thailand (AE.CM240 and AE.93TH253). The amino acid sequences of the env V4 and C3 regions in the samples were highly variable, compared with those of V3 and V3-downstream regions. The V3 loop central motif in the majority (87.5%) of the strains was GPGQ. The majority of strains did not contain positively charged amino acids at positions 306 and 320 in V3 loop. The N-linked glycosylation sites in the V3-V4 region and flanking regions in these strains were relatively conserved. Analysis of the entire region showed that the mean Ks values were significantly higher than that of the Ka values (P < 0.001), with the Ks/Ka significantly higher than 1.0 (P < 0.001). In contrast, the Ks/Ka ratio in the V4 region was significantly lower than 1.0 (P < 0.01).
Conclusions: Our study indicated that the majority of HIV-1 subtype CRF01-AE strains circulating in China were highly homogeneous. The amino acid sequences of the V4 and C3 regions were significantly more variable than those of the V3 loop. Our analysis also suggested that the phenotype of nearly all strains was likely to be non-syncytium inducing (NSI). Finally, the variation found in the V3-V4 sequence was significantly influenced by functional constraints as opposed to positive selective pressure, while the variability of the lone V4 region was strongly related to positive selective pressure.