Phylogenetic and Drug-Resistance Analysis of HIV-1 Sequences From an Extensive Paediatric HIV-1 Outbreak in Larkana, Pakistan

• Published in: Frontiers in microbiology Vol. 12; p. 658186
• Main Authors: Abidi, Syed Hani, Nduva, George Makau, Siddiqui, Dilsha, Rafaqat, Wardah, Mahmood, Syed Faisal, Siddiqui, Amna Rehana, Nathwani, Apsara Ali, Hotwani, Aneeta, Shah, Sharaf Ali, Memon, Sikander, Sheikh, Saqib Ali, Khan, Palwasha, Esbjörnsson, Joakim, Ferrand, Rashida Abbas, Mir, Fatima
• Format: Journal Article
• Language: English
• Published: Switzerland Frontiers Media S.A 17.08.2021
• Subjects: Clinical Medicine; drug resistance; HIV-1; Infectious Medicine; Infektionsmedicin; Klinisk medicin; Medical and Health Sciences; Medicin och hälsovetenskap; Microbiology; outbreak investigation; paediatric [MeSH]; phylogenetic analysis; HIV-1; phylogenetic analysis; paediatric [MeSH]; drug resistance; outbreak investigation
• ISSN: 1664-302X; 1664-302X
• DOI: 10.3389/fmicb.2021.658186

In April 2019, an HIV-1 outbreak among children occurred in Larkana, Pakistan, affecting more than a thousand children. It was assumed that the outbreak originated from a single source, namely a doctor at a private health facility. In this study, we performed subtype distribution, phylogenetic and drug-resistance analysis of HIV-1 sequences from 2019 outbreak in Larkana, Pakistan. A total of 401 blood samples were collected between April-June 2019, from children infected with HIV-1 aged 0-15 years recruited into a case-control study to investigate the risk factors for HIV-1 transmission. Partial HIV-1 sequences were generated from 344 blood plasma samples to determine HIV-1 subtype and drug resistance mutations (DRM). Maximum-likelihood phylogenetics based on outbreak and reference sequences was used to identify transmission clusters and assess the relationship between outbreak and key population sequences between and within the determined clusters. Bayesian analysis was employed to identify the time to the most recent common recent ancestor (tMRCA) of the main Pakistani clusters. The HIV-1 circulating recombinant form (CRF) 02_AG and subtype A1 were most common among the outbreak sequences. Of the treatment-naïve participants, the two most common mutations were RT: E138A (8%) and RT: K219Q (8%). Four supported clusters within the outbreak were identified, and the median tMRCAs of the Larkana outbreak sequences were estimated to 2016 for both the CRF02_AG and the subtype A1 clusters. Furthermore, outbreak sequences exhibited no phylogenetic mixing with sequences from other high-risk groups of Pakistan. The presence of multiple clusters indicated a multi-source outbreak, rather than a single source outbreak from a single health practitioner as previously suggested. The multiple introductions were likely a consequence of ongoing transmission within the high-risk groups of Larkana, and it is possible that the so-called Larkana strain was introduced into the general population through poor infection prevention control practices in healthcare settings. The study highlights the need to scale up HIV-1 prevention programmes among key population groups and improving infection prevention control in Pakistan.