A HIV-1 heterosexual transmission chain in Guangzhou, China: a molecular epidemiological study
- Zhigang Han†1,
- Tommy WC Leung†2,
- Jinkou Zhao†3,
- Ming Wang1,
- Lirui Fan1,
- Kai Li4,
- Xinli Pang4,
- Zhenbo Liang4,
- Wilina WL Lim2 and
- Huifang Xu1Email author
© Han et al; licensee BioMed Central Ltd. 2009
Received: 25 June 2009
Accepted: 25 September 2009
Published: 25 September 2009
We conducted molecular analyses to confirm four clustering HIV-1 infections (Patient A, B, C & D) in Guangzhou, China. These cases were identified by epidemiological investigation and suspected to acquire the infection through a common heterosexual transmission chain.
Env C2V3V4 region, gag p17/p24 junction and partial pol gene of HIV-1 genome from serum specimens of these infected cases were amplified by reverse transcription polymerase chain reaction (RT-PCR) and nucleotide sequenced.
Phylogenetic analyses indicated that their viral nucleotide sequences were significantly clustered together (bootstrap value is 99%, 98% and 100% in env, gag and pol tree respectively). Evolutionary distance analysis indicated that their genetic diversities of env, gag and pol genes were significantly lower than non-clustered controls, as measured by unpaired t-test (env gene comparison: p < 0.005; gag gene comparison: p < 0.005; pol gene comparison: p < 0.005).
Epidemiological results and molecular analyses consistently illustrated these four cases represented a transmission chain which dispersed in the locality through heterosexual contact involving commercial sex worker.
Epidemiology of human immunodeficiency virus type 1 (HIV-1) infection in China is changing from predominantly injecting drug use to increasingly sexual transmission , with close to 50% of new infections attributable to sexual transmission in 2005 . The unsafe sexual behavior is one of the key risk factors of HIV transmission in China. A review cited the Newsweek reported that the number of commercial sex workers (CSWs) in mainland China exceed 10 million in 2003 . Data from a nation-wide sentinel surveillance revealed that HIV prevalence among CSWs has risen from 0.02% in 1996 to 0.93% in 2004 and over 1% in some places . Recently, a comprehensive surveillance revealed that 60% of CSWs do not use condoms every time . These results suggested that CSWs serve as a bridge to transmit HIV from core risk groups into general population of the country.
What puzzles epidemiologists is the difficulties in defining the linkage amongst HIV cases using traditional epidemiological approaches based on behavioral information. Nucleotide sequence analysis presents an unique opportunity to identify possible epidemiological linkage between infected cases and to track the viral transmission from person to person through viral genome analysis, as explored in various HIV-1 studies [5, 6].
We conducted molecular analyses to confirm four HIV-1 infections from Guangzhou China which were transmitted through a heterosexual transmission chain involving CSW.
Patients, materials and methods
Serum specimens from four patients (Patients A, B, C and D) and five contacts (Contacts 1, 2, 3, 4 and 5) in this surveillance investigation were collected between January 2008 and May 2008 after obtaining their informed consent. These sera were collected prior to any antiretroviral treatment. Sera were screened for HIV-1 antibody by enzyme-linked immunosorbent assays (bioMerieux, France and Peking BGI-GBI, China) and confirmed by Western blotting (MP Biomedicals, Singapore). HIV-1 positive sera were aliquoted to avoid unnecessary freeze-thaw deterioration cycle and were stored at -70°C.
RNA extraction/PCR/nucleotide sequencing
RNA extraction, reverse transcription polymerase chain reaction (RT-PCR) amplification and nucleotide sequencing were performed in physically separated laboratories. Viral RNA was extracted by QIAamp Viral RNA extraction Mini kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. Extracted RNA was reverse transcribed by MuLV reverse transcriptase (Applied Biosystems, Inc., Foster City, CA) into cDNA using random hexamer. The cDNA was used as template for PCR amplification of two HIV-1 subgenomic regions (env C2V3V4 region and gag p17/p24 junction) by nested PCR. In addition, the cDNA from four clustered cases was further amplified at the pol gene. The primers and conditions of PCR for gag and env applied in this study were as previously described . The region of pol gene was divided into three segments (pro, rt5' and rt3'). The primers and conditions for pol were as follows: pro segment: outer sense primer (PolF20): 5'-GAG AGA CAG GCT TAT TTT TT-3', common antisense primer (PR96): 5'-CTT CCC AGA AGT CTT GAG TTC T-3', inner sense primer (PolF21): 5'-GCA GAC CAG AGC CAA CAG C-3'; rt5' segment: outer sense primer (A1): 5'-AAT TTT CCC ATT AGT CCT ATT-3', outer antisense (NE1): 5'-TAT GTC ATT GAC AGT CCA GCT-3', inner sense (NNA): 5'-AAG CCA GGA ATG GAT GGC CCA-3', inner antisense (E): 5'-CCA TTT ATC AGG ATG GAG TTC-3'; rt3' segment: outer sense (RTF31): 5'-CCA CAC CAG ACA AAA ARC ATC-3', outer antisense (PolR40): 5'-CTG TTA CTA TGT TTA CTT CT-3', inner sense (RTF32): 5'-CAT CAG AAA GAA CCY CCA TT-3', inner antisense (PolR38): 5'-TTA GCT GCC CCA TCT ACA TAG-3'. PCR products were sequenced by the ABI PRISM BigDye Terminator Cycle Sequencing Ready Reaction Kit (Version 3.1) (Applied Biosystems). The sequenced products were resolved and analyzed using an ABI PRISM 3100 or 3130 xl Sequence Detection System (Applied Biosystems) . Relevant positive and negative controls were included in each time to avoid false positive in the PCR.
Phylogenetic analyses of env, gag and pol genes
Nucleotide sequences from Patients A through D together with other non-clustered control sequences were aligned by ClustalW software  and followed by manual adjustment. Phylogenetic and molecular evolutionary analyses were conducted using MEGA version 4 . Evolutionary distances were calculated by Kimura two-parameter modeling, excluding positions with alignment gaps in any sequence. Phylogenetic dendrograms based on partial sequence of the env C2V3V4 region, gag p17/p24 junction and pol gene of the HIV-1 genome were constructed using neighbor-joining method with Kimura two-parameter modeling. The reliability of each node was evaluated by bootstrapping with 2000 replicates.
Amino acid sequence analysis
Gag p17/p24 junction nucleotide sequences from clustered cases were translated into amino acid and aligned with their closest CRF01_AE reference sequence (Accession number: AF197340) and the highest homology matched sequence detected in Yunnan China in 2002 (Accession number: EF062020). Signature analysis was conducted by VESPA tool of Los Alamos National Laboratory http://www.hiv.lanl.gov to identify unique sequence pattern present in those clustering cases as compare with the above-mentioned reference sequences .
Differences in degree of diversity (evolutionary distance) in env, gag and pol genes of clustered cases and non-clustered controls were ascertained using unpaired t-test.
HIV-1 viral RNA was extracted from the sera of Patients A through D. HIV-1 subtype determination using env, gag and pol genes of the HIV-1 genome of these patients consistently showed they belonged to CRF01_AE . Based on the genetic characterization of gag p17/p24 junction sequences, it revealed that these four CRF01_AE viral strains were equally homologous (95%) to CRF01_AE isolates detected in Thailand at 1999 (GenBank Accession: AY945731) and Yunnan at 2002 (GenBank Accession: EF062020).
Evolutionary distances within the four patients, controls and between patients and controls on env, gag and pol genes a
5.23 ± 0.85
n = 4
1.19 × 10-11
2.39 ± 0.60
n = 4
6.37 × 10-13
1.54 ± 0.22
n = 4
5.12 × 10-14
11.02 ± 0.92
n = 23
6.27 ± 0.66
n = 32
2.95 ± 0.18
n = 25
14.21 ± 1.48
8.23 ± 1.15
4.12 ± 0.37
Gag p17/p24 amino acid sequence analysis
HIV-1 is characterized by high genomic variability. Sequence diversity was observed whether among isolates within the same individual or between isolates from different infected individuals. Previous studies have shown that HIV-1 isolates from same individual may differ by up to 2% in the env gene and those from unrelated individuals may differ by 6-22% , while the degree of diversity between isolates from closely related individuals fall in between .
Phylogenetic analyses are increasingly used in clarifying the epidemiological linkage of HIV-1 transmission by comparing nucleotide sequence fragments from one or more subgenomic regions, such as gag, env and pol genes, or full-length genome sequence. Early studies have concentrated on sequence variation in the env V3 region, which contains the principal neutralization domain and determinants for biological phenotypes [13, 14]. Thus, there was argument that it could be less suitable for epidemiological study than other regions (gag and pol) of HIV-1 genome [15, 16]. However, Leitner et al. accurately reconstructed a known HIV-1 transmission history by phylogenetic analysis and demonstrated that it was at least as accurate using env V3 sequences as with gag p17 sequences . Recent data showed that the env, gag or pol genes were frequently used in epidemiological studies [18–22].
Our evolutionary distance analyses at three different genomic regions (env, gag and pol genes) of the HIV-1 viral genome demonstrated that the HIV-1 isolates from these four patients formed a cluster which was highly related with each other and differed by 5.23% in the env gene. Molecular evidence strongly supported the close epidemiological linkage among them. In addition, analyses showed that the average evolutionary distance in env, gag and pol genes of these four clustering members were significantly lower than the ones from local unrelated non-cluster CRF01_AE controls [5, 7].
By far bootstrap test is one of the most commonly used tests of the reliability of a constructed phylogenetic dendrogram. The study of Hillis and Bull showed that a bootstrap value of >70% commonly signified a probability of 95% or higher and the topology at this branch is real . Phylogenetic analyses illustrated that the sequences from the four patients clustered together in phylogenetic analyses of both env, gag and pol genes. The high bootstrap values (99% for env gene, 98% for gag gene and 100% for pol gene, respectively) strongly indicated that they are of closely related and monophyletic.
Amino acid signature was observed in some study to detect for closely related cases, such as transmission of HIV-1 from a dentist to his five patients . The study of a possible single-source sexual transmission cluster in upstate New York showed that a specific amino acid signature was present in the gag p17/p24 junction of subtype B sequence , while the study of three clusters among subtype B samples in Hong Kong also indicated that amino acid signatures were unique and cluster specific . In this study we detected a unique signature pattern of four amino acid residues in the gag p17/p24 junction of CRF01_AE sequences of these four clustering cases in comparison with selected reference sequences.
Genetic analyses of gag p17/p24 junction of these four patients revealed that these four CRF01_AE viral isolates were equally homologous (95%) to CRF01_AE isolate which predominated amongst infected heterosexuals in Thailand at 1999 and Yunnan at 2002. It is possible that CRF01_AE strain was transmitted from Thailand through Yunnan China and eventually to Guangzhou China.
In summary, we performed molecular epidemiology to track the transmission of HIV-1 infection within a group of heterosexual patients in Guangzhou. Highly related sequences from four patients indicated a transmission chain. The result complemented the epidemiological findings that the infection was sustained within the locality through heterosexual contact involving a CSW.
Nucleotide sequences of the env, gag and pol gene of Patient A - D were submitted to GenBank at the NCBI (Accession no. FJ752409 - FJ752412, FJ752413 - FJ752416 and FJ752417 - FJ752420 for env, gag and pol sequences respectively).
This study was supported by a grant provided by the Bureau of Science and Technology of Guangzhou Municipality (Grant 2006Z1-E0093). Special thanks to relevant staff at Guangzhou Municipal CDC and Huadu District CDC for immunoassays, and Patients and Contacts for their voluntary participation. Written consent for publication was obtained from the Patients and Contacts.
- State council AIDS Working committee Office, UN Theme Group on AIDS in China: A Joint Assessment of HIV-AIDS Prevention Treatment and Care in China (2007). Beijing 2008. 10.1016/S0140-6736(07)60315-8Google Scholar
- Wu ZY, Sullivan SG, Wang Y, Rotheram-Borus MJ, Detels R: Evolution of China's response to HIV/AIDS. Lancet 2007, 369: 679-690. 10.1097/01.olq.0000162360.11910.5aView ArticlePubMedGoogle Scholar
- Yang H, Li X, Stanton B, Liu H, Liu H, Wang N, Fang X, Lin D, Chen X: Heterosexual transmission of HIV in China: a systematic review of behavioral studies in the past two decades. Sex Transm Dis 2005, 32: 270-280. 10.1126/science.256.5060.1165PubMed CentralView ArticlePubMedGoogle Scholar
- Ministry of Health of China, UNAIDS, WHO: 2005 update on the HIV/AIDS epidemic and response in China. Beijing 2006.Google Scholar
- Ou CY, Ciesielski CA, Myers G, Bandea CI, Luo CC, Korber BT, Mullins JI, Schochetman G, Berkelman RL, Economou AN, Witte JJ, Furman LJ, Satten GA, Maclnnes KA, Curran JW, Jaffe HW, Laboratory Investigation Group, Epidemiologic Investigation Groupl: Molecular epidemiology of HIV transmission in a dental practice. Science 1992, 256: 1165-1171. 10.1089/aid.2007.0272View ArticlePubMedGoogle Scholar
- Brooks JT, Robbins KE, Youngpairoj AS, Rotblatt H, Kerndt PR, Taylor MM, Daar ES, Kalish ML: Molecular analysis of HIV strains from a cluster of worker infections in the adult film industry, Los Angeles 2004. AIDS 2006, 20: 923-928. 10.1093/nar/22.22.4673View ArticlePubMedGoogle Scholar
- Leung TW, Mak D, Wong KH, Wang Y, Song YH, Tsang DN, Wong C, Shao YM, Lim WL: Molecular epidemiology demonstrated three emerging clusters of human immunodeficiency virus type 1 subtype B infection in Hong Kong. AIDS Res Hum Retroviruses 2008, 24: 903-910. 10.1093/molbev/msm092View ArticlePubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994, 22: 4673-4680. 10.1089/aid.1992.8.1549PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular Biology and Evolution 2007, 24: 1596-1599. 10.1073/pnas.88.24.11236View ArticlePubMedGoogle Scholar
- Korber B, Myers G: Signature pattern analysis: a method for assessing viral sequence relatedness. AIDS Res Hum Retroviruses 1992, 9: 1549-1560. 10.1073/pnas.85.6.1932View ArticleGoogle Scholar
- Myers G, Rabson AB, Berzofsky JA, Smith TF, Wong-Staal F, (Eds): Human Retroviruses and AIDS 1990: A Compilation and Analysis of Nucleic Acid and Amino Acid Sequences. Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM;
- Burger H, Weiser B, Flaherty K, Gulla J, Nguyen PN, Gibbs RA: Evolution of human immunodeficiency virus type 1 nucleotide sequence diversity among close contacts. Proc Natl Acad Sci USA 1991, 88: 11236-11240. 10.1073/pnas.85.9.3198PubMed CentralView ArticlePubMedGoogle Scholar
- Palker TJ, Clark ME, Langlois AJ, Matthews TJ, Weinhold KJ, Randall RR, Bolognesi DP, Haynes BF: Type-specific neutralization of the human immunodeficiency virus with antibodies to env-encoded synthetic peptides. Proc Natl Acad Sci USA 1988, 85: 1932-1936. 10.1073/pnas.85.6.1932PubMed CentralView ArticlePubMedGoogle Scholar
- Rusche JR, Javaherian K, McDanal C, Petro J, Lynn DL, Grimaila R, Langlois A, Gallo RC, Arthur LO, Fischinger PJ, Bolognesi DP, Putney SD, Matthews TJ: Antibodies that inhibit fusion of human immunodeficiency virus-infected cells bind a 24-amino acid sequence of the viral envelope, gp120. Proc Natl Acad Sci USA 1988, 85: 3198-3202. 10.1038/364766b0PubMed CentralView ArticlePubMedGoogle Scholar
- Holmes EC, Zhang LQ, Robertson P, Cleland A, Harvey E, Simmonds P, Leigh Brown AJ: The molecular epidemiology of human immunodeficiency virus type 1 in Edinburgh. J Infect Dis 1995, 171: 45-53. 10.1073/pnas.93.20.10864View ArticlePubMedGoogle Scholar
- Holmes EC, Brown AJ, Simmonds P: Sequence data as evidence. Nature (London) 1993, 364: 766. 10.1038/364766b0View ArticleGoogle Scholar
- Leitner T, Escanilla D, Franzén C, Uhlén M, Albert J: Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis. Proc Natl Acad Sci USA 1996, 93: 10864-10869. 10.1073/pnas.93.20.10864PubMed CentralView ArticlePubMedGoogle Scholar
- Albert J, Wahlberg J, Leitner T, Escanilla D, Uhlén M: Analysis of a Rape Case by Direct Sequencing of the Human-immunodeficiency-virus type-1 pol and gag genes. J Virol 1994, 68: 5918-5924. 10.1073/pnas.222522599PubMed CentralPubMedGoogle Scholar
- Blanchard A, Ferris S, Chamaret S, Guétard D, Montagnier L: Molecular evidence for nosocomial transmission of human immunodeficiency virus from a surgeon to one of his patients. J Virol 1998, 72: 4537-4540. 10.1128/JVI.74.6.2525-2532.2000PubMed CentralPubMedGoogle Scholar
- Metzker ML, Mindell DP, Liu XM, Ptak RG, Gibbs RA, Hillis DM: Molecular evidence of HIV-1 transmission in a criminal case. Proc Natl Acad Sci USA 2002, 99: 14292-14297. 10.1016/j.jcv.2003.08.008PubMed CentralView ArticlePubMedGoogle Scholar
- Goujon CP, Schneider VM, Grofti J, Montigny J, Jeantils V, Astagneau P, Rozenbaum W, Lot F, Frocrain-Herchkovitch C, Delphin N, Le Gal F, Nicolas JC, Milinkovitch MC, Dény P: Phylogenetic analyses indicate an atypical nurse-to-patient transmission of human immunodeficiency virus type 1. J Virol 2000, 74: 2525-2532. 10.1128/JVI.74.6.2525-2532.2000PubMed CentralView ArticlePubMedGoogle Scholar
- Pistello M, Del Santo B, Buttò S, Bargagna M, Domenici R, Bendinelli M: Genetic and phylogenetic analyses of HIV-1 corroborate the transmission link hypothesis. J Clin Virol 2004, 30: 11-18. 10.1089/088922202320567914View ArticlePubMedGoogle Scholar
- Hillis DM, Bull JJ: An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analyses. Syst Biol 1993, 42: 182-192.View ArticleGoogle Scholar
- Robbins KE, Weidle PJ, Brown TM, Saekhou AM, Coles B, Holmberg SD, Folks TM, Kalish ML: Molecular analysis in support of an investigation of a cluster of HIV-1-infected women. AIDS Res Hum Retroviruses 2002, 18: 1157-1161. 10.1089/088922202320567914View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.