Abstract
Characterization of HIV-1 strains is important for surveillance of the HIV-1 epidemic. In Vietnam HIV-1-infected pregnant women often fail to receive the care they are entitled to. Here, we analyzed phylogenetically HIV-1 env sequences from 37 HIV-1-infected pregnant women from Ha Noi (n=22) and Hai Phong (n=15), where they delivered in 2005–2007. All carried CRF01_AE in the gp120 V3 region. In 21 women CRF01_AE was also found in the reverse transcriptase gene. We compared their env gp120 V3 sequences phylogenetically in a maximum likelihood tree to those of 198 other CRF01_AE sequences in Vietnam and 229 from neighboring countries, predominantly Thailand, from the HIV-1 database. Altogether 464 sequences were analyzed. All but one of the maternal sequences colocalized with sequences from northern Vietnam. The maternal sequences had evolved the least when compared to sequences collected in Ha Noi in 2002, as shown by analysis of synonymous and nonsynonymous changes, than to other Vietnamese sequences collected earlier and/or elsewhere. Since the HIV-1 epidemic in women in Vietnam may still be underestimated, characterization of HIV-1 in pregnant women is important to observe how HIV-1 has evolved and follow its molecular epidemiology.
M
In Vietnam HIV-1 infection is now endemic in some of the largest cities, such as Ho Chi Minh City, with a prevalence of 1.2%. 1,2 The most profound prevalence of the HIV-1 epidemic in northern Vietnam is in the coastal province of Quang Ninh, which borders southern China, in the port city Hai Phong (1.1%), and in the capital Ha Noi (0.9%). 2 In pregnant women the prevalence in these cities was 0.6% and 0.1%, respectively, in 2006. 3
The genetic epidemiology of HIV-1 has become increasingly important in understanding how, wherefrom, and whereto the HIV-1 epidemic is changing. The env gene is the most variable and is commonly used for the definition and characterization of genotypes. 4 The major (M) group of HIV-1 genotypes is divided into subtypes A–D, F–H, J–K, and 49 circulating recombinant forms (CRF). In Vietnam, CRF01_AE has been reported to be the dominant HIV-1 genotype, even more common than in the neighboring countries. 5 –8
Due to the widespread discrimination against HIV-1-infected people in Vietnam, 9 –11 HIV-1-infected pregnant women fear being revealed. 12 –15 The number of pregnant women having an HIV test in antenatal care decreased from 22% in 2004 to about 16% in 2006, 16 but has since increased. 17 Still, only 25% received the test in 2009. 17 HIV-1-infected pregnant women choose elective abortion in close to 70%, 18 which may be a reflection of their fear.
Thus HIV-1 in pregnant women in Vietnam has not been studied, most likely due to the difficulty of obtaining samples from them. It is not known to what extent HIV-1 in pregnant women is a reflection of HIV-1 in different risk groups, such as intravenous drug users or commercial sex workers, or has spread beyond this group. Therefore we wanted to explore the phylogenetic relation of HIV-1 in pregnant women to other strains of HIV-1 from northern Vietnam and elsewhere.
The 37 HIV-infected mothers studied here were part of a prospective study in 2004–2007 concerning follow-up of children to 135 HIV-1-infected mothers, the first in Vietnam. 19 They were followed prospectively after a gynecological check up and delivery in 2005–2007 at the Obstetrics and Gynecology hospitals in Ha Noi (n=22) and Hai Phong (n=15), respectively. 19 All pregnant women were offered an HIV test and counseling, if found to be HIV infected. All consenting women (182/243) were offered participation for a thorough follow-up of their children up to a year of age; 135 mothers participated. Most women did not provide a sample from themselves. We have had access to 37 samples for the present study. The small number may restrict the conclusions that could be drawn, but still represents the first report from this important group. They included most (7/11) mothers of infected children and a selection of mothers of uninfected children, collected during the same time and from the same cities. Their age ranged from 17 to 36 with a median of 25, which was similar to other groups of pregnant women. The CD4 counts were less than 200×106/liter in 5/7 mothers of infected children and over 200×106/liter in the others tested. Six women were treated with combination therapy within 3 weeks before delivery. None of those women had a child who contracted HIV-1. The others were treated with nevirapine at delivery (n=25), or not at all (n=6). None of the children was breastfed.
The National Institute of Hygiene and Epidemiology and the Obstetric and Gynecology Hospitals in Ha Noi and Hai Phong had given permission for the study. All subjects were treated respectfully and signed consents were obtained, according to the Helsinki declaration of human rights from 2000. Also the regional board of ethical vetting in Stockholm had approved of the study on February 26, 2008: 2207/1496-31/3. The study was registered at the Clinical trials.gov with ID: NCT 00669604 at (
The HIV status of the women was discovered when they went to the hospital for an antenatal check up late in gestation or in most cases at delivery. Three milliliters of defibrinated venous blood was taken 1 or 2 days after delivery from the mothers and processed within 8 h. Plasma was separated from the cellular fraction. Peripheral blood mononuclear cells were isolated by Lymphoprep (Axis-Shield PoC AS, Oslo, Norway). Two million cells were lysed in 0.2 ml buffer containing protease K and detergent as a source of HIV-1 DNA. RNA was extracted from 140 μl plasma (nine samples) by QIAamp (QIAGEN Gmbh, Hilden Germany). Viral RNA was used as a template to run cDNA synthesis by Omniscript kit (QIAGEN, Hilden, Germany). To avoid contamination of samples tubes were not opened at the same time. A positive and a negative control of HIV-1-infected cells and water, respectively, were included in each run. The sequences gave no indication that contamination had occurred.
The V3 loop region of the HIV-1 env gene was amplified in a nested PCR, using the outer primers JA167 and JA170 and inner primers JA168 and JA169. 4 In samples with remaining material the reverse transcriptase of the pol gene was also amplified for different HIV-1 subtypes in a nested PCR, using the outer primers GAG2, PR1, RT137, and RT3303. The inner PCR primers RT1 and RT4 spanned the reverse transcriptase gene, amino acids 30 to 227 of the pol gene. 20 Sequencing was performed using Big Dye Terminator 3.1 (AB Applied Biosystems, Foster City, CA). Sequencing reactions were purified by the Dye Ex 2.0 spin kit (QIAGEN). The product was run and sequenced on the ABI PRISM 3700 DNA sequencer (Applied Biosystems, Foster City, CA) or submitted to MWG Operon Eurofins Company, Erbesberg, Germany. The length of the V3 sequences was about 350 bp for the env gene and 646 bp for the reverse transcriptase gene. The chromatograms were evaluated using the program Sequencher 4.1 (Gene Codes, Ann Arbor, MI).
To determine the phylogenetic relation between the new and present CRF01_AEenvV3 sequences from the pregnant women from northern Vietnam to sequences characterized in the past in Vietnam and neighboring countries, 14 sequences, obtained from intravenous drug users in Ha Noi from 2002, 21 were added, as well as all other (n=413) available sequences of CRF01_AE from the Los Alamos HIV database, retrieved in June 2008. The criteria of selection were that they should include the C2V3C3 region of the env protein. One sequence per individual was selected at random to avoid bias, reducing the number of sequences used. The number of sequences from each region was over or around 100 and could be expected to represent the epidemic in that area. The numbers of maternal sequences and of intravenous drug users were less, but may still be a reflection of their risk group. Sequences from pregnant women are the first reported from this group in Vietnam. The papers, from which the Vietnamese and Chinese sequences were obtained, were reviewed. 5 –7,21 –25
A maximum likelihood phylogenetic tree was inferred using PhyML v3.0 (PHYlogenetic inferences using maximum likelihood), using a GTR+G+I model with BioNJ program initial tree and an NNI search. A tree diagram was plotted with FigTree v1.1.2 by Andrew Rambaut, available at [
A total of 464 sequences were analyzed; 159 came from Thailand, 43 from China, and 198 from Vietnam. The remaining came from other countries such as Japan (n=8), Indonesia (n=7), Singapore (n=2), and Korea (n=9). The sequences were presented in a phylogenetic tree (Fig. 1), which displayed the nodes in ascending order. Each sequence can be identified in the supplementary file (Supplementary Fig. S1; Supplementary Data are available online at

Phylogenetic tree of the CRF01_AE envV3 molecular epidemic in Vietnam and related countries in Southeast Asia. The phylogenetic tree shows the sequences from Thailand in green, from China in orange, from the south of Vietnam in blue, and from the north of Vietnam in purple. Sequences from other countries are shown in black. The indicated scale bar is in units of substitutions/site according to a GTR+I+G model. The maternal sequences are indicated by filled black circles.
The closest relationships between the sequences from northern Vietnam and a foreign source was found with the HIV-1 Pingxiang strains of HIV-1-infected individuals from the Guangxi province 5,6 in southern China, which borders with the Lang Son province in northern Vietnam.
There was an interesting separation of the Chinese sequences between the upper and lower part of the tree, as all 22 of the Chinese sequences in the upper part came from areas at the border between China and Vietnam and colocalized with sequences from northern Vietnam, while 22/23 (95.7%) of the Chinese sequences in the lower part of the tree, dominated by sequences from Thailand, came from Liaonging and Fuijan 25 in the east of China. These latter ones were, similar to the samples from the pregnant women, of relatively recent date (from 2000 to 2006), indicating that there may be two (ongoing) paths of entry of CRF01_AE into China.
The visual impression of the maternal sequences (black circles) was that they were scattered throughout the other sequences, collected from northern Vietnam and southern China, apart from one sequence, which appeared together with sequences from southern Vietnam.
Another way to study the relation to other groups of HIV-infected people would be to describe the synonymous (ds) and nonsynonymous (dn) amino acid changes, where the synonymous changes would represent a background of mutations and nonsynonymous would increase in relation to different pressures. Both would accumulate over time. Studying the individual nucleotide changes according to whether they were ds or dn showed that the maternal sequences had significantly larger numbers of synonymous changes than those of the intravenous drug users in Ha Noi in 2002 and others from northern Vietnam, respectively (p<0.001) See Table 1. The nonsynonymous changes were not significantly different from the intravenous drug users, p>0.001, but were significantly higher than the older ones from northern Vietnam (p<0.001). Altogether, the comparisons of the respective dn and ds means may reflect the time-related evolution of the maternal sequences versus the older ones from Vietnam.
Extending the analysis to geographically more distant groups, the maternal sequences had less synonymous changes than those from southern Vietnam, but more with regard to nonsynonymous changes. With regard to sequences from Thailand, these were significantly more evolved than the maternal sequences and all others with regard to both synonymous and nonsynonymous changes (p<0.001). This may be a sign of the longstanding epidemic in Thailand, broader range of sources, and a wider range in time of collection.
In Lang Son in northern Vietnam the HIV-1 prevalence rates among intravenous drug users increased dramatically between 1995 and 1998 from 0% to 12% 29 and in Pingxiang, Guangxi, from 8% to 42% in a year. 30 This was interpreted to be related to the major heroin transshipment routes from the Golden Triangle [Myanmar (Burma), Thailand, and Laos] to Hong Kong, by separate routes for different subtypes. 23,30,31 For CRF01_AE this route goes from Laos into Ha Noi and then turns north. 30,31 Since there were no sequences from Myanmar and Laos, a relation of HIV-1 to those countries may be unrevealed.
Recently the epidemic of HIV CRF01A_AE in Vietnam was analyzed phylogenetically, involving near full length sequences, collected in 1997–1998. 28 The epidemic first seemed to come from Thailand by heterosexual spread into the south of Vietnam and spread to intravenous drug users in the late 1980s and subsequently to intravenous drug users in the north during the mid-1990s. It is of interest that the sequences in our tree displayed a similar pattern. However, also their study was hampered by a lack of sequences from Myanmar and Laos.
It would be important to know if the epidemic now may be moving into the general population, here represented by pregnant women. This would be expected to be associated with a greater degree of evolution than the rapid spread associated with intravenous drug use. 32,33 Thus it has been shown that HIV-1 evolves less among intravenous drug users than among heterosexuals. 32 This has also been observed in the context of CRF01_AE in northern Vietnam and Southern China since the 1990s to more recently. 5,6,21,30
The analysis of synonymous and nonsynonymous nucleotide changes indicated that the maternal sequences had evolved when comparing them to the older sequences from northern Vietnam. It was not possible to discern whether HIV-1 in pregnant women has evolved as a result of spreading into the general population, or is still closely related to the epidemic among intravenous drug users, a group that dominates as a source of HIV sequences in this part of the country. In a recent study on the pol gene phylogeny of sequences from Northern Vietnam, there were signs of HIV spread in different groups, such as among intravenous drug users and from sexual spread. 34
The samples here come from an important and sensitive group, from which it is difficult to receive samples and information. Since the HIV epidemic in women in Vietnam may still be underestimated, 9 characterization of HIV-1 in pregnant women is important to assist in the surveillance of mother-to-child transmission of HIV-1.
Accession Numbers
The following accession numbers of previously reported sequences from intravenous drug users in Ha Noi were used: AY283413–AY283427. The following accession numbers were given to the new HIV-1 env sequences created for this publication: FJ009330–FJ009366.
Footnotes
Acknowledgments
Economic support was given by SIDA (Swedish International Development Cooperation Agency) for costs in northern Vietnam and to A.E., SWE 2006-467. Tran Ha was given a grant from the Karolinska Institutet Research Training Program (KIRT) and the Swedish Nutrition Foundation. Simani Gaseitsiwe was awarded a Ph.D. training grant and Tran Ha and Rozina Caridha received a 6 month training grant each, all from EU Marie Curie Early Training program, and support from the Karolinska Institutet. P.V.H. received travel grants from the Maud and Birger Gustavsson's Foundation. Tran Ha held an NIH-DOE interagency agreement grant: Y1-AI-1500-07. We would like to thank Thomas Leitner for scientific support in generating the phylogenetic tree, Ina Maljkovic Berry for similar support and for critical evaluations and discussions, and her and Joakim Coster for reviewing the manuscript. R. Caridha and T.T.T. Ha contributed equally to this work and shared first authorship.
Author Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
