HIV Phylogenetics Reveals Overlapping Transmission Networks among Cities and Key Populations in Pakistan

The first case of HIV in Pakistan was documented in 1987, with multiple subtypes and circulating recombinant forms being introduced and currently circulating in the country. Since then, there has been a shift in the country from a low prevalence/high risk to a high-risk concentrated epidemic. Pakistans epidemic is concentrated among key populations at greater risk of HIV infection including people who inject drugs (PWID), Hijra sex workers (HSW), female sex workers (FSW), male sex workers (MSW), and men who have sex with men (MSM). This study focused on the geographical aspect as well as on the interactions between key populations at higher risk of contracting HIV. We aimed at understanding the behavior of these key populations at a molecular level with high granularity as well as investigating the possibility of multiple HIV-1 introductions in Pakistan. In this cross-sectional biological and behavioral survey, we collected dried blood spots (DBS) for the purposes of seroprevalence estimates and molecular epidemiology from individuals in 17 cities in Pakistan representing four key populations: PWID, HSW, MSW, and FSW. A total of 1153 envelope sequences (reference positions in HXB2: 7860-8274) of HIV were sequenced using a Sanger-based sequencing approach. To identify clusters based on the introduction of the virus in Pakistan from foreign countries we added 3623 publicly available HIV envelope sequences to our dataset. Phylogeographic inference suggests at least 15 independent introductions of the virus into Pakistan, with a total of 12 clusters ranging from 3 to 675 sequences in size containing sequences from Pakistan and neighboring countries exclusively. Our phylogenetic analysis shows a significant degree of connectivity and directionality suggesting broad and overlapping networks of HIV-1 transmission among cities and key populations in Pakistan.