Evidence for a “Founder Effect” among HIV-infected injection drug users (IDUs) in Pakistan
The first evidence of Human immunodeficiency virus (HIV) infection in Pakistan was reported in 1987 among patients receiving tainted blood or blood products . Having been relatively safe from any indigenous HIV cases for around two decades, this conservative South-east Asian nation is finally registering isolated HIV outbreaks all over the country.
In 2004, in the remote desert town of Larkana, Pakistan experienced its first HIV outbreak . The implicated populace was a community of Injection Drug Users (IDU), who documented an HIV prevalence approaching an outstanding 27%. Earlier, HIV prevalence in IDUs was reported by the National AIDS Control Program as bordering around 0.4% in December, 2003. After the Larkana episode, outbreaks were also recorded in other major cities of Pakistan . The changing trend denotes that the country is transitioning from a low prevalence to a concentrated epidemic stage.
Apart from IDUs, other high-risk populations including truck drivers , commercial sex workers, including Hijras  have also been identified in the Pakistani context. Even though very little serological data are available, various knowledge, attitude and practice studies [3,6-9] have demonstrated an increasingly vulnerable position for a devastatingly spread of HIV in these high-risk groups.
Earlier, we have reported an HIV subtype A concentrated epidemic in a community of IDUs in Karachi, Pakistan's largest city . The study presented here is an extension of that previous study on the same community of IDUs. In the previous study we showed that HIV-1 subtype A was the only subtype of HIV prevalent among this community. From the data we gathered, there were indications that this community of IDUs was infected by HIV positive overseas commercial workers (OCW) deported to Pakistan during early nineties . We now present evidence that the subtype A infection in question was a single-source infection, most likely brought into the IDU community by an OCW.
Ethical approval for this study was obtained from the Ethical Review Committee, Aga Khan University, Karachi, Pakistan. HIV-1 positive patients were selected from a community of IDUs in Karachi, the Southern Port City of Pakistan. After obtaining informed consent, 3 mL of blood was collected from 26 men (mean age, 24 years) previously diagnosed positive for HIV based on sero-testing.
DNA was extracted from the blood samples using methods previously described . PCR amplification was carried out for HIV gag gene using two sets of primers in a nested strategy. The primers used in nested PCR were GOPF (5' CTCTCGACGCAGGACTCGGCTTGC-3', nt 683-706, HXB2) and GOPR (5'- CCAATTCCCCCTATCATTTTTGG-3', nt 2382-2404) for the first round of amplification and primers GIPF (5'- GAGGCTAGAAGGAGAGAGATGGG-3', nt 772-794) and GIPR (5'-TTATTGTGACGAGGGGTCGTTGCC-3', nt 2269-2292) for the second round of PCR.
The reaction mixture of 25µl for both first and second round PCR contained 1X PCR buffer (5X Green GoTaq® Flexi Buffer, pH 8.5), 2mM MgCl2, 400µM dNTPs and 0.3U of Taq Polymerase. The first round of PCR was performed with 0.48 pmol of primers GOPF and GOPR. Thermocycle was: denaturation at 95oC for 5 min, followed by 35 cycles of denaturation at 95oC for 1 min, annealing at 58oC for 1 min and extension at 72oC for 1 min, with a final extension of at 72oC for 15 min.
1µl of the first-round PCR product along with 0.48 pmol of the primers GIPF and GIPR was used for the second-round PCR. Thermocycle was: denaturation at 95oC for 5 min, followed by 35 cycles of denaturation at 95oC for 1 min, annealing at 60oC for 1 min and extension at 72oC for 1 min, with a final extension of at 72oC for 15 min. The amplified products were electrophoresed on 1.2% agarose gel, stained by ethidium bromide and visualized under ultraviolet light.
Sequencing and Phylogenetic Analysis
Nested PCR products of gag gene were partially sequenced from Macrogen Inc, Korea, using primer GSP1 (5'- CCATCAATGAGGAAGCTGC-3', nt 1400-1418, HXB2). The nucleotide sequence spanning the p24 and p7 region of gag gene surrounded by the primers H1Gag1584 (5'- AAAGATGGATAATCCTGGG-3') and g17 (5'-TCCACATTTCCAACAGCCCTTTTT-3') , which has been considered for accurate subtyping in previous studies was analyzed. Analysis of the resulting sequences was carried out using Bioinformatic tools described previously . Briefly, the sequences (comprising 460-470 bp) were compared with the sequences from the Los Alamos HIV sequence database. This was accomplished by using the HIV BLAST Search (http://www.hiv.lanl.gov/). The samples were assigned subtypes based on the closest homology found with the subtypes in the Los Alamos database.
Using the same sequence, alignments were obtained by the Clustal X program (1.83) . After alignment, gaps were stripped and minor manual adjustments were made using MacClade . From these alignments, phylogenetic relationships were determined by constructing phylogenetic trees using neighbor-joining method with the help of PAUP* software . Pairwise genetic distances were calculated with Kimura's two parameter method .
Results and Discussion
Sequence alignment and comparison of the partial HIV-1 gag gene from all twenty six IDUs revealed highly congruent topologies. Based on the phylogenetic analysis, the HIV-1 strains from Pakistan clustered closely with HIV-1 subtype A sequences from Senegal and Uganda. The HIV-1 subtype A strains from Pakistan were similar to those circulating in Africa, as opposed to those in neighboring India. Moreover, the HIV-1 strains from Pakistan formed a monophyletic group validating our initial hypothesis of a possible ‘founder effect' among the HIV strains from IDU in Karachi.
An identity matrix constructed by BioEdit software revealed more than 98% similarity among the 26 sequences analyzed, further confirming the ‘founder effect' (Additional file 1). In an independent experiment, randomly selected 15 samples were analyzed in a similar way, this time using sequences from HIV nef gene. Once again, a phylogenetic tree was obtained showing high degree of monophylogeny among the studied samples (Additional file 2).
The HIV epidemic is still in its formative stages in Pakistan . The work presented here is the first report that analyses the molecular demographics of HIV spread in this conservative South-East nation. We have demonstrated that the strains circulating amongst this select community of IDUs in Karachi are most likely a consequence of a single-source infection. This finding may well be applicable to other IDU communities in Pakistan that are recording an explosive increase in HIV infection.
In a previous report on the same community of IDUs we indicated to the possible origin of this particular endemic being the overseas commercial workers who had lived in the Middle East over extended period during which they admitted to having had sexual contact with commercial sex workers. We now report deeper analysis of data from the same samples. Close phylogenetic clustering and strong identity among these HIV sequences indicates a “founder effect” of HIV in this community of IDUs, supporting the hypothesis that this entire group of 26 IDUs may have been infected with HIV from a single source. The HIV infection most likely got transmitted to the entire group through exchange of contaminated needles. The study subject HIV-GAG35 has previously reported to having had contact with commercial sex workers during a prolong stay in Sharja, Middle East . During the early nineties, the subject in question was deported, after being tested positive for HIV in Sharja, to Pakistan; thereafter he developed a habit for injecting drugs. The rest of the IDUs were diagnosed positive for HIV in the year 2000 and later. Based on the fact that HIV positive HIV-GAG35 sequence clustered closely with the rest of the IDU samples (Figure 1), we speculate that this patient could be a possible source of infection in this community. That the origin of this particular infection resides in the Middle East is an attractive theory that may be corroborated further through collection and analysis of HIV positive samples from the countries in that region.
During the early nineties, an influx of HIV positive Pakistani deportees from the Middle East constituted a major threat to the transmission of infection in the country . The Pakistani government did respond, albeit inadequately, to this looming disaster. Although commendable efforts were made to devise effective measures to control HIV transmission in the country, a structured program to screen and educate Pakistani migrant workers regarding sexually transmitted infections is still awaited. We have indicated before , and emphasize again, that with constant human traffic across borders from Middle East, Afghanistan, and India, and the prevalence of HIV in Pakistan on the rise, clear-cut policies on HIV prevention across borders is urgently warranted.
The HIV/AIDS epidemic is following along the same atypical lines it has followed in the rest of Asia. Beginning from infection in the high-risk populace, the virus soon traverses the barrier to the general population. Pakistan is already experiencing multiple HIV outbreaks in the IDU community in geographically diverse regions. It would not be long before a steep rise in infection occurs, not unlike its HIV-havocked neighbor, India . HIV-1 is circulating mostly within close-knit communities in Pakistan, such as the IDUs. It is therefore an early opportunity for implementing targeted HIV prevention programs in these vulnerable communities before the infection becomes an issue for the general populace.
The authors declare that they have no competing interests.
MAR: Performed PCR and sequence analysis and wrote the first draft of the paper. VRN: Developed the study and provided guidance, space and financial support for MAR's work. SKho: Performed PCR and sequence analysis for HIV gag gene. SKha: Helped in collection and extraction of samples. AR: Helped in accessing and surveying the IDU community. RY: carried out part of sequence analysis. SHA: Supervised the study and arranged for partial funding. All authors have read and approved the submitted manuscript.
1. Khanani RM, Hafeez A, Rab SM, Rasheed S. Human immunodeficiency virus-associated disorders in Pakistan. AIDS Res Hum Retroviruses 1988;4(2):149-54.
2. Shah SA, Altaf A, Mujeeb SA, Memon A. An outbreak of HIV infection among injection drug users in a small town in Pakistan: potential for national implications. Int J STD AIDS 2004;15(3):209.
3. Platt L, Vickerman P, Collumbien M, Hasan S, Lalji N, Mayhew S, Muzaffar R, Andreasen A, Hawkes S. Prevalence of HIV, HCV and sexually transmitted infections among injecting drug users in Rawalpindi and Abbottabad, Pakistan: evidence for an emerging injection-related HIV epidemic. Sex Transm Infect. 2009; 85 Suppl 2:ii17-22.
4. Agha S. Potential for HIV transmission among truck drivers in Pakistan. Aids 2000;14(15):2404-6.
5. Baqi S, Kayani N, Khan JA. Epidemiology and clinical profile of HIV/AIDS in Pakistan. Trop Doct 1999;29(3):144-8.
6. Jamal N, Khushk IA, Naeem Z. Knowledge and attitudes regarding AIDS among female commercial sex workers at hyderabad city pakistan. J Coll Physicians Surg Pak 2006;16(2):91-3.
7. Luby SP, Niaz Q, Siddiqui S, Mujeeb SA, Fisher-Hoch S. Patients' perceptions of blood transfusion risks in Karachi, Pakistan. Int J Infect Dis 2001;5(1):24-6.
8. Bokhari A, Nizamani NM, Jackson DJ, Rehan NE, Rahman M, Muzaffar R, et al. HIV risk in Karachi and Lahore, Pakistan: an emerging epidemic in injecting and commercial sex networks. Int J STD AIDS 2007;18(7):486-92.
9. Raheel H, White F, Kadir MM, Fatmi Z. Knowledge and beliefs of adolescents regarding sexually transmitted infections and HIV/AIDS in a rural district in Pakistan. J Pak Med Assoc 2007;57(1):8-11.
10. Khan S, Rai MA, Khanani MR, Khan MN, Ali SH. HIV-1 subtype A infection in a community of intravenous drug users in Pakistan. BMC Infect Dis 2006;6:164.
11. Shah SA, Khan OA, Kristensen S, Vermund SH. HIV-infected workers deported from the Gulf States: impact on Southern Pakistan. Int J STD AIDS 1999;10(12):812-4.
12. Heyndrickx L, Janssens W, Zekeng L, Musonda R, Anagonou S, Van der Auwera G, et al. Simplified strategy for detection of recombinant human immunodeficiency virus type 1 group M isolates by gag/env heteroduplex mobility assay. Study Group on Heterogeneity of HIV Epidemics in African Cities. J Virol 2000;74(1):363-70.
13. Khoja S, Ojwang P, Khan S, Okinda N, Harania R, Ali S. Genetic analysis of HIV-1 subtypes in Nairobi, Kenya. PLoS One, 2008. 3(9): p. e3191.
14. Jeanmougin F, Thompson JD, Gouy M, Higgins DG, Gibson TJ. Multiple sequence alignment with Clustal X. Trends Biochem Sci 1998;23(10):403-5.
15. Maddison DR, Maddison WP. MacClade version 4 Analysis of phylogeny and character evolution. Sunderland Massachusetts: Sinauer Associates, 2000.
16. Swofford DL. PAUP*: Phylogenetic analysis using parsimony (*and other methods) beta version 4.0b8. Sunderland, Massachusetts: Sinauer Associates, 2001.
17. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 1980;16(2):111-20.
18. Rai MA, Warraich HJ, Ali SH, Nerurkar VR. HIV/AIDS in Pakistan: the battle begins. Retrovirology 2007;4:22.
19. Chandrasekaran P, Dallabetta G, Loo V, Rao S, Gayle H, Alexander A. Containing HIV/AIDS in India: the unfinished agenda. Lancet Infect Dis 2006;6(8):508-21.