Timing and evolution of the most recent common ancestor of the Korean clade HIV subtype B based on <Emphasis Type="Italic">Nef</Emphasis> and <Emphasis Type="Italic">Vif</Emphasis> sequences

Molecular phylogenetic studies of the HIV-1 isolated from Koreans have suggested the presence of the so-called “Korean clade”, which can be defined as a cluster free of foreign isolates. The Korean clade accounts for more than 60% of Korean isolates and exerts characteristic amino acid sequences. Thus, it is merited to estimate when this Korean clade first emerged in order to understand the evolutionary pattern of the Korean clade. We analyzed and reconstructed the most recent common ancestor (MRCA) sequences from nef (n=229) and vif (n=179) Korean clade sequences. Linear regression analyses of sequence divergence estimates were plotted against sampling years to infer the year in which there was zero divergence from the MRCA sequences. MRCA sequences suggested the Korean clade was first emerged around 1984, before the first detection of HIV-1 in Korea in 1985. Further studies on synonymous and nonsynonymous substitution rates suggested positive selection event for the Korean clade, while other subtype B had undergone negative to neutral evolution.     

Genetic analysis of human immunodeficiency virus type 1 strains from patients in Cyprus: identification of a new subtype designated subtype I.

DNA sequences encoding the C2 to V3 region of envelope glycoprotein gp120 of human immunodeficiency virus type 1 (HIV-1) were amplified by PCR from uncultured peripheral blood mononuclear cells obtained from 24 of 25 HIV-1-seropositive patients from Cyprus. By using a heteroduplex mobility assay (HMA), all amplified products were studied genetically and compared with 16 previously characterized HIV-1 strains belonging to subtypes A through F. HMA results revealed that HIV-1 gp120 sequences from 15 of our patients were of subtype B of HIV-1, whereas one isolate was of subtype C. However, gp120 sequences from eight patients had no obvious similarities to the known subtypes as defined by HMA. DNA sequencing and phylogenetic analyses of molecular clones confirmed the HMA results and placed the eight undefined HIV-1 isolates into three distinct genetic clusters. On the basis of branch topology and lengths of the phylogenetic tree, we conclude that one group consisting of three clones from two patients represents a new HIV-1 env subtype, which we have termed subtype I. The remaining two sequence clusters, consisting of five sequences from four patients and two sequences from two other patients, are distally related to subtypes A and F. These data demonstrate the extensive heterogeneity of HIV-1 in Cyprus, including the presence of new subtype.     

Extensive intrasubtype recombination in South African human immunodeficiency virus type 1 subtype C infections

Recombinant human immunodeficiency virus type 1 (HIV-1) strains containing sequences from different viral genetic subtypes (intersubtype) and different lineages from within the same subtype (intrasubtype) have been observed. A consequence of recombination can be the distortion of the phylogenetic signal. Several intersubtype recombinants have been identified; however, less is known about the frequency of intrasubtype recombination. For this study, near-full-length HIV-1 subtype C genomes from 270 individuals were evaluated for the presence of intrasubtype recombination. A sliding window schema (window, 2 kb; step, 385 bp) was used to partition the aligned sequences. The Shimodaira-Hasegawa test detected significant topological incongruence in 99.6% of the comparisons of the maximum-likelihood trees generated from each sequence partition, a result that could be explained by recombination. Using RECOMBINE, we detected significant levels of recombination using five random subsets of the sequences. With a set of 23 topologically consistent sequences used as references, bootscanning followed by the interactive informative site test defined recombination breakpoints. Using two multiple-comparison correction methods, 47% of the sequences showed significant evidence of recombination in both analyses. Estimated evolutionary rates were revised from 0.51%/year (95% confidence interval [CI], 0.39 to 0.53%) with all sequences to 0.46%/year (95% CI, 0.38 to 0.48%) with the putative recombinants removed. The timing of the subtype C epidemic origin was revised from 1961 (95% CI, 1947 to 1962) with all sequences to 1958 (95% CI, 1949 to 1960) with the putative recombinants removed. Thus, intrasubtype recombinants are common within the subtype C epidemic and these impact analyses of HIV-1 evolution.     

Amino acid sequence divergence of Tat protein (exon1)of subtype B and C HIV-1 strains: Does it have implications for vaccine development?

Functional genes of HIV-1 like the tat express proteins essential for viral survival and propagation. There are variations reported in levels of Tat transactivation among the different subtypes of HIV-1. This study looked at the amino acid differences in the different regions of Tat protein (exon 1) of subtype B and C strains of HIV-1 and tried to observe a molecular basis for protein function. HIV-1 sequences of subtype B (n=30) and C (n=60) strains were downloaded from HIV-1 Los Alamos data base. Among the 60 subtype C strain sequences, 30 each were from India and Africa. A HIV-1 Tat protein (exon 1) sequence, the consensus B and C sequence was obtained from the 'sequence search interface' in the Los Alamos HIV-1 sequence data. The sequences were visualized using Weblogo and the RNA binding regions of the three consensus sequences were also determined using BindN software program. Compared to subtype B, there was a high level of divergence in the auxiliary domain of tat exon 1 (amino acid positions 58- 69). The net charge of the subtype C (Indian) Tat protein (exon 1) auxiliary domain was -1.9 at pH 7 and it had an isoelectric point of 4.1. The net charge of the subtype C (African) auxiliary domain was -2.9 at pH 7 and it had an isoelectric point of 3.7 while the net charge of same region in subtype B was -0.9 at pH 7 with an isoelectric point of 4.9. The ratio of the hydrophilic residues to the total number of residues was 60% in the in both the Indian and African subtype C in the auxiliary domain while this was 50% in subtype B. The consensus subtype B sequence was found to have 36 RNA binding sites while subtype C (India) had 33 and subtype C (Africa) had 32 RNA binding sites. The HIV-1 Tat-TAR interaction is a potential target for inhibitors and being considered for its potential use in HIV-1 vaccines. Development of such inhibitor/vaccines would have to take into consideration the variation in amino acid sequence analyzed in this study as this could determine epitope presentation on MHC class I antigen for afferent immune response.     

Phenotypic and genotypic characteristics of human immunodeficiency virus type 1 from patients with AIDS in northern Thailand.

Primary human immunodeficiency virus type 1 (HIV-1) isolates were obtained from 22 patients with AIDS from northern Thailand, where HIV-1 is transmitted primarily through the heterosexual route. Viral sequences were determined for the 22 patients with AIDS, and all were subtype E HIV-1 on the basis of sequence analysis of a region from the envelope protein gp120. Syncytium-inducing (SI) viruses were detected for 16 of 22 patients with AIDS by using MT-2 cells. Characteristics of amino acid sequences in V3 which have not been reported previously for subtype B SI HIV-1 were associated with the subtype E HIV-1 SI phenotype. The SI viruses from our study population contain predominantly a GPGR or GPGH motif at the tip of the V3 loop, in contrast to the previously described subtype E HIV-1 from Thailand which contained predominantly GPGQ. All the SI viruses lost a potential N-linked glycosylation site in V3 which is highly conserved among previously described subtype E HIV-1 isolates from asymptomatic patients from Thailand. HIV-1 envelope sequences including V3 from some patients with AIDS were significantly more divergent than viruses from asymptomatic patients in Thailand characterized 2 years ago or earlier. These results suggest that emergence of subtype E SI HIV-1 variants is associated with the development of AIDS, as it is for subtype B HIV-1. The divergence of subtype E HIV-1 in patients with AIDS as the disease progresses, and the divergence of subtype E HIV-1 in the infected population as the epidemic continues in Thailand, may have important implications for vaccine development.     

广西壮族自治区HIV-1流行毒株的基因序列测定和亚型分析

使用ＰＣＲ技术对１４份广西ＨＩＶ－１阳性感染者外周血单核细胞（ＰＢＭＣｓ）样品进行扩增,获得ＨＩＶ－１膜蛋白（ｅｎｖ）基因的核酸片段,并对其Ｃ２－Ｖ３及邻区３５０－４５０个核苷酸序列进行了测定和分析。结果表明,１４份样品中９份为泰国Ｂ（Ｂ′）亚型,５份为Ｅ亚型毒株。其中Ｂ′亚型毒株的基因离散率为４．２％,与Ａ－Ｅ参考亚型及部分Ｂ亚型代表株序列相比较,与包括泰国、缅甸及云南德宏在内的Ｂ亚型毒株序列十分接近,相互之间基因离散率在３．０％－４．４％的范围内;而Ｅ亚型毒株的基因离散率为２．１％,与国际Ｅ亚型毒株的基因离散率最近,为５．６％,与其它国际参考亚型基因离散率很远,在２１．１％－２７．３％。根据以上数据及其它资料提示,广西存在Ｂ′和Ｅ两种亚型的ＨＩＶ－１的流行,且其Ｂ′亚型毒株的传入,与流行在云南德宏州的相同亚型ＨＩＶ－１毒株密切相关,而Ｅ亚型毒株则可能是由泰国经越南传入广西的     

Genetic analysis of human immunodeficiency virus type 1 and 2 (HIV-1 and HIV-2) mixed infections in India reveals a recent spread of HIV-1 and HIV-2 from a single ancestor for each of these viruses.

DNA sequences encoding the surface envelope glycoproteins of human immunodeficiency virus type 1 (HIV-1) and type 2 (HIV-2) were amplified by PCR from uncultured peripheral blood mononuclear cells obtained from patients with serologically defined HIV-1/HIV-2 mixed infections from Bombay, India. HIV-1-specific PCR products were obtained in seven of seven randomly chosen doubly reactive cases, while HIV-2-specific sequences were detected in five of seven cases (71%). DNA sequence analysis showed that the HIV-1 gp120 coding sequences were closely related to each other (nucleotide sequence divergence of between 3.1 and 6.8%). Phylogenetic tree analysis placed the Indian strains within the C subtype of HIV-1, being most similar to sequences previously found in East and South Africa. The HIV-2 sequences were also closely related to each other, with an overall sequence divergence of between 5.6 and 10.5%. The low level of nucleotide divergence among Indian HIV-1 and HIV-2 sequences suggests a fairly recent introduction of each virus into this population from a single point of entry in each case. The HIV-2 sequences reported here represent the first analysis of Asian HIV-2 strains and confirm the serological pattern previously detected in India. These data show that a substantial spread of HIV-2, together with HIV-1, has appeared outside Africa in a population hitherto unexposed to HIV. These findings imply that further spread of HIV-2 worldwide is to be expected and have important implications for future vaccine and therapy development.     

Phylogeographic analysis of HIV-1 subtype C dissemination in Southern Brazil

The HIV-1 subtype C has spread efficiently in the southern states of Brazil (Rio Grande do Sul, Santa Catarina and Paraná). Phylogeographic studies indicate that the subtype C epidemic in southern Brazil was initiated by the introduction of a single founder virus population at some time point between 1960 and 1980, but little is known about the spatial dynamics of viral spread. A total of 135 Brazilian HIV-1 subtype C pol sequences collected from 1992 to 2009 at the three southern state capitals (Porto Alegre, Florianópolis and Curitiba) were analyzed. Maximum-likelihood and Bayesian methods were used to explore the degree of phylogenetic mixing of subtype C sequences from different cities and to reconstruct the geographical pattern of viral spread in this country region. Phylogeographic analyses supported the monophyletic origin of the HIV-1 subtype C clade circulating in southern Brazil and placed the root of that clade in Curitiba (Paraná state). This analysis further suggested that Florianópolis (Santa Catarina state) is an important staging post in the subtype C dissemination displaying high viral migration rates from and to the other cities, while viral flux between Curitiba and Porto Alegre (Rio Grande do Sul state) is very low. We found a positive correlation (r(2) = 0.64) between routine travel and viral migration rates among localities. Despite the intense viral movement, phylogenetic intermixing of subtype C sequences from different Brazilian cities is lower than expected by chance. Notably, a high proportion (67%) of subtype C sequences from Porto Alegre branched within a single local monophyletic sub-cluster. These results suggest that the HIV-1 subtype C epidemic in southern Brazil has been shaped by both frequent viral migration among states and in situ dissemination of local clades.     

新疆乌鲁木齐HIV—1流行毒株膜蛋白基因C2—V3区序列测定和亚型分析

应用ＰＣＲ方法对７份１９９６年３－５月采集于新疆乌鲁木齐ＨＩＶ－１阳性静脉吸毒者的外周血单核细胞样品进行扩增,获得了ＨＩＶ－１膜蛋白基因的核酸片段,并对其Ｃ２－Ｖ３区及邻区３５０－４５０个核苷酸序列进行了测定和分析。     

Origin and dynamics of HIV-1 subtype C infection in India

Objective

To investigate the geographical origin and evolution dynamics of HIV-1 subtype C infection in India.

Design

Ninety HIV-1 subtype C env gp120 subtype C sequences from India were compared with 312 env gp120 reference subtype C sequences from 27 different countries obtained from Los Alamos HIV database. All the HIV-1 subtype C env gp120 sequences from India were used for the geographical origin analysis and 61 subtype C env gp120 sequences with known sampling year (from 1991 to 2008) were employed to determine the origin of HIV infection in India.

Methods

Phylogenetic analysis of HIV-1 env sequences was used to investigate the geographical origin and tMRCA of Indian HIV-1 subtype C. Evolutionary parameters including origin date and demographic growth patterns of Indian subtype C were estimated using a Bayesian coalescent-based approach under relaxed molecular clock models.

Findings

The majority of the analyzed Indian and South African HIV-1 subtype C sequences formed a single monophyletic cluster. The most recent common ancestor date was calculated to be 1975.56 (95% HPD, 1968.78–1981.52). Reconstruction of the effective population size revealed three phases of epidemic growth: an initial slow growth, followed by exponential growth, and then a plateau phase approaching present time. Stabilization of the epidemic growth phase correlated with the foundation of National AIDS Control Organization in India.

Interpretation

Indian subtype C originated from a single South African lineage in the middle of 1970s. The current study emphasizes not only the utility of HIV-1 sequence data for epidemiological studies but more notably highlights the effectiveness of community or government intervention strategies in controlling the trend of the epidemic.     

Recent Evolutionary History of Human Immunodeficiency Virus Type 1 Subtype B: Reconstruction of Epidemic Onset Based on Sequence Distances to the Common Ancestor

We obtained and studied HIV-1 sequences with a known sampling year from three outbreaks of the HIV-1 epidemic: 141 env V3 (270 nt) sampled between 1984 and 1992 and 117 pol prot/RT (804 nt) sequences sampled between 1986 and 1999 from Dutch homosexual men and injecting drug users (IDUs), as well as 77 env V3 sequences sampled between 1983 and 1994 in the United States. Since retrospective serological and/or epidemiological data on these populations are available, providing estimates of the dates of the onset of the HIV-1 epidemics, we had the opportunity to test different phylogenetic models for their accuracy in deriving the recent evolutionary history of HIV-1 subtype B and the onset date of the HIV-1 epidemic. We observed that, in any given year, individual sequences vary widely in their distances to the common ancestor, and sequences close to the ancestors were found decades after the onset of the epidemic. Nevertheless, the mean evolutionary distances of virus strains to ancestors were increasing significantly during the course of the studied epidemics, which indicates that the molecular clock is operational in the recent evolution of HIV-1. When the relationship between the sampling years of sequences and their nucleotide distances to the common ancestor was extrapolated to the past, analysis of pol sequences provided accurate estimates of the onset years of the epidemics, whereas analysis of V3 sequences by the maximum-likelihood or neighbor-joining methods led to an overestimation of the age of the epidemics. Separate analysis of nonsynonymous and synonymous distances revealed that this overestimation results from nonsynonymous substitutions, whose numbers were not increasing significantly in all three virus populations over the observation period. In contrast, analysis of synonymous env V3 distances provided accurate estimates of the onset years for the outbreaks we studied. Received: 26 October 2001 / Accepted: 8 November 2001     

Dating the time of viral subtype divergence

Precise dating of viral subtype divergence enables researchers to correlate divergence with geographic and demographic occurrences. When historical data are absent (that is, the overwhelming majority), viral sequence sampling on a time scale commensurate with the rate of substitution permits the inference of the times of subtype divergence. Currently, researchers use two strategies to approach this task, both requiring strong conditions on the molecular clock assumption of substitution rate. As the underlying structure of the substitution rate process at the time of subtype divergence is not understood and likely highly variable, we present a simple method that estimates rates of substitution, and from there, times of divergence, without use of an assumed molecular clock. We accomplish this by blending estimates of the substitution rate for triplets of dated sequences where each sequence draws from a distinct viral subtype, providing a zeroth-order approximation for the rate between subtypes. As an example, we calculate the time of divergence for three genes among influenza subtypes A-H3N2 and B using subtype C as an outgroup. We show a time of divergence approximately 100 years ago, substantially more recent than previous estimates which range from 250 to 3800 years ago.     

Subtypes of human immunodeficiency virus type 1 and disease stage among women in Nairobi, Kenya

In sub-Saharan Africa, where the effects of human immunodeficiency virus type 1 (HIV-1) have been most devastating, there are multiple subtypes of this virus. The distribution of different subtypes within African populations is generally not linked to particular risk behaviors. Thus, Africa is an ideal setting in which to examine the diversity and mixing of viruses from different subtypes on a population basis. In this setting, it is also possible to address whether infection with a particular subtype is associated with differences in disease stage. To address these questions, we analyzed the HIV-1 subtype, plasma viral loads, and CD4 lymphocyte levels in 320 women from Nairobi, Kenya. Subtype was determined by a combination of heteroduplex mobility assays and sequence analyses of envelope genes, using geographically diverse subtype reference sequences as well as envelope sequences of known subtype from Kenya. The distribution of subtypes in this population was as follows: subtype A, 225 (70.3%); subtype D, 65 (20.5%); subtype C, 22 (6.9%); and subtype G, 1 (0.3%). Intersubtype recombinant envelope genes were detected in 2.2% of the sequences analyzed. Given that the sequences analyzed represented only a small fraction of the proviral genome, this suggests that intersubtype recombinant viral genomes may be very common in Kenya and in other parts of Africa where there are multiple subtypes. The plasma viral RNA levels were highest in women infected with subtype C virus, and women infected with subtype C virus had significantly lower CD4 lymphocyte levels than women infected with the other subtypes. Together, these data suggest that women in Kenya who are infected with subtype C viruses are at more advanced stages of immunosuppression than women infected with subtype A or D. There are at least two models to explain the data from this cross-sectional study; one is that infection with subtype C is associated with a more rapid disease progression, and the second is that subtype C represents an older epidemic in Kenya. Discriminating between these possibilities in a longitudinal study will be important for increasing our understanding of the role of specific subtypes in the transmission and pathogenesis of HIV-1.     

Genetic heterogeneity of HIV-1 in Greece

The aim of this study was to detect and determine the genetic variation of HIV-1 in Greece and to analyze the phylogenetic relationships and transmission dynamics of identified variants. Eighty-six blood samples from HIV-1 seroconverted patients of different risk groups were collected from the AIDS clinic, AHEPA Hospital, Thessaloniki, Greece. Retroviral DNA was extracted from uncultured peripheral blood mononuclear cells. HIV-1 DNA sequences encoding a 500-bp fragment of the gp120 C2-C3 region were amplified from each study subject, and they were genetically subtyped by heteroduplex mobility assay and DNA sequencing. Genetic distances and phylogenetic relationships of DNA sequences were estimated using PHYLIP software. Our results revealed that 82 out of 86 (95.3%) subjects carried subtype B sequences, while four (4.7%) carried subtype A sequences. Subtype A in Greek individuals not having traveled abroad was documented. An average of intrasubtype B genetic divergence of 15% was noted. Our findings demonstrate the presence of at least two genetic subtypes of HIV-1 in northern Greece--subtype B and subtype A. The predominant subtype is subtype B, which was transmitted into Greece by multiple sources. Our observations lend support to the argument that the distribution of HIV-1 subtypes is determined by founder effects or other processes rather than any tropism for particular cell types or mode of transmission.     

Evolutionary modeling of rate shifts reveals specificity determinants in HIV-1 subtypes

A hallmark of the human immunodeficiency virus 1 (HIV-1) is its rapid rate of evolution within and among its various subtypes. Two complementary hypotheses are suggested to explain the sequence variability among HIV-1 subtypes. The first suggests that the functional constraints at each site remain the same across all subtypes, and the differences among subtypes are a direct reflection of random substitutions, which have occurred during the time elapsed since their divergence. The alternative hypothesis suggests that the functional constraints themselves have evolved, and thus sequence differences among subtypes in some sites reflect shifts in function. To determine the contribution of each of these two alternatives to HIV-1 subtype evolution, we have developed a novel Bayesian method for testing and detecting site-specific rate shifts. The RAte Shift EstimatoR (RASER) method determines whether or not site-specific functional shifts characterize the evolution of a protein and, if so, points to the specific sites and lineages in which these shifts have most likely occurred. Applying RASER to a dataset composed of large samples of HIV-1 sequences from different group M subtypes, we reveal rampant evolutionary shifts throughout the HIV-1 proteome. Most of these rate shifts have occurred during the divergence of the major subtypes, establishing that subtype divergence occurred together with functional diversification. We report further evidence for the emergence of a new sub-subtype, characterized by abundant rate-shifting sites. When focusing on the rate-shifting sites detected, we find that many are associated with known function relating to viral life cycle and drug resistance. Finally, we discuss mechanisms of covariation of rate-shifting sites.     

Substitution Model of Sequence Evolution for the Human Immunodeficiency Virus Type 1 Subtype B gp120 Gene over the C2-V5 Region

Phylogenetic analyses frequently rely on models of sequence evolution that detail nucleotide substitution rates, nucleotide frequencies, and site-to-site rate heterogeneity. These models can influence hypothesis testing and can affect the accuracy of phylogenetic inferences. Maximum likelihood methods of simultaneously constructing phylogenetic tree topologies and estimating model parameters are computationally intensive, and are not feasible for sample sizes of 25 or greater using personal computers. Techniques that initially construct a tree topology and then use this non-maximized topology to estimate ML substitution rates, however, can quickly arrive at a model of sequence evolution. The accuracy of this two-step estimation technique was tested using simulated data sets with known model parameters. The results showed that for a star-like topology, as is often seen in human immunodeficiency virus type 1 (HIV-1) subtype B sequences, a random starting topology could produce nucleotide substitution rates that were not statistically different than the true rates. Samples were isolated from 100 HIV-1 subtype B infected individuals from the United States and a 620 nt region of the env gene was sequenced for each sample. The sequence data were used to obtain a substitution model of sequence evolution specific for HIV-1 subtype B env by estimating nucleotide substitution rates and the site-to-site heterogeneity in 100 individuals from the United States. The method of estimating the model should provide users of large data sets with a way to quickly compute a model of sequence evolution, while the nucleotide substitution model we identified should prove useful in the phylogenetic analysis of HIV-1 subtype B env sequences. Received: 4 October 2000 / Accepted: 1 March 2001     

Evolution of the human immunodeficiency virus type 1 subtype-specific V3 domain is confined to a sequence space with a fixed distance to the subtype consensus.

Human immunodeficiency virus type 1 (HIV-1) strains can be separated into genetic subtypes based on phylogenetic analysis of the envelope gene. Once it had been shown that population-wide intrasubtype genetic variation of HIV-1 strains increases in the course of the AIDS epidemic, it remained uncertain whether HIV-1 subtypes are phenotypic entities spreading as distinct virus populations. To examine this, we applied Eigen's concepts of sequence geometry and fitness topography to the analysis of intrasubtype evolution of the gp120 V3 domain of HIV-1 subtypes A, B, C, and D in the course of the global AIDS epidemic. We observed that despite the high evolution rate of HIV-1, the nonsynonymous distances to the subtype consensus of sequences obtained early in the epidemic are similar to those obtained more than 10 years later, in contrast to the synonymous distances, which increased steadily over time. For HIV-1 subtype B, we observed that the evolution rate of the individual sequences is independent of their distance from the subtype B consensus, but for the individual sequences most distant from the consensus evolution away from the consensus is constrained. As a result, individual HIV-1 genomes fluctuate within a sequence space with fixed distance to the subtype consensus. Our findings suggest that the evolution of the V3 domain of HIV-1 subtypes A, B, C, and D is confined to an area in sequence space within a fixed distance to the consensus of a respective subtype. This in turn indicates that each HIV-1 subtype is a distinct viral quasispecies that is well adapted to the present environment, able to maintain its identity in the V3 region over time, and unlikely to merge during progression of the AIDS epidemic.