Fast Dissemination of New HIV-1 CRF02/A1 Recombinants in Pakistan

Identification of a novel first-generation HIV-1 circulating recombinant form (CRF152_DG) among people living with HIV in Karachi, Pakistan

The objective of this study was to characterize a novel circulating recombinant form of human immunodeficiency virus type 1 (HIV-1) among people living with HIV in Karachi, Pakistan. We conducted near-full-length genome (NFLG) sequencing on eight samples exhibiting D/G recombination signals in the pol gene region. We successfully obtained NFLG sequences (790-9,614; with reference to the HXB2 genome) from four of the eight samples and then conducted phylogenetic and recombination analyses on them. The four NFLG sequences from our study and one DG unique recombinant form previously identified in the United Kingdom (GenBank accession: MF109700) formed a distinct monophyletic cluster with an Shimodaira-Hasegawa approximate likelihood ratio test node support value of 100%. Bootscan analyses of the five NFLG sequences of DG recombinants showed that all five NFLGs shared the same unique mosaic pattern of recombination breakpoints between D and G clades, with two D fragments in the pol and vif regions inserted into a G backbone. Subregion phylogenetic analyses confirmed these sequences to be a novel circulating recombinant form (CRF) composed of subtypes D and G. The DG recombinant sequences were eventually designated as CRF152_DG by the Los Alamos HIV Sequence Database staff.

IMPORTANCE

In Pakistan, the genetic diversity of human immunodeficiency virus type 1 (HIV-1) is becoming increasingly complex, compared to the early years of the epidemic that started after the detection of the first cases of HIV-1 in 1987 in Karachi. Based on the available molecular studies, two dominant HIV-1 clades, sub-subtype A1 and CRF02_AG, have been found to co-circulate with other clades, namely B, C, D, G, CRF01_AE, CRF35_A1D, and CRF56_cpx, in various urban areas of Pakistan. Several novel recombinant forms have also been detected. This first report of CRF152_DG highlights the complex nature of the HIV epidemic in Pakistan and emphasizes the importance of continual molecular surveillance (ideally based on whole-genome sequences) of HIV.

Human Immunodeficiency Virus Type-1 Genetic Diversity and Drugs Resistance Mutations among People Living with HIV in Karachi, Pakistan

The human immunodeficiency virus type-1 epidemic in Pakistan has significantly increased over the last two decades. In Karachi, Pakistan, there is a lack of updated information on the complexity of HIV-1 genetic diversity and the burden of drug resistance mutations (DRMs) that can contribute to ART failure and poor treatment outcomes. This study aimed to determine HIV-1 genetic diversity and identify drug-resistance mutations among people living with HIV in Karachi. A total of 364 HIV-positive individuals, with a median age of 36 years, were enrolled in the study. The HIV-1 partial pol gene was successfully sequenced from 268 individuals. The sequences were used to generate phylogenetic trees to determine clade diversity and also to assess the burden of DRMs. Based on the partial pol sequences, 13 distinct HIV-1 subtypes and recombinant forms were identified. Subtype A1 was the most common clade (40%), followed by CRF02_AG (33.2%). Acquired DRMs were found in 30.6% of the ART-experienced patients, of whom 70.7%, 20.7%, and 8.5% were associated with resistance to NNRTIs, NRTIs, and PIs, respectively. Transmitted DRMs were found in 5.6% of the ART-naïve patients, of whom 93% were associated with resistance against NNRTIs and 7% to PIs. The high prevalence of DRMs in ART-experienced patients poses significant challenges to the long-term benefits and sustainability of the ART program. This study emphasizes the importance of continuous HIV molecular epidemiology and drug resistance surveillance to support evidence-based HIV prevention, precise ART, and targeted AIDS care.

Characterization of HIV-1 CRF02_AG/A3/G unique recombinant forms identified among children in Larkana, Pakistan

Co-circulation of different human immunodeficiency virus type 1 HIV-1 subtypes among infected populations can lead to the generation of new recombinants. In Pakistan, subtype A1 and CRF02_AG are the dominant strains circulating among key populations. The high prevalence of new HIV infections among the key populations highlights the possibility of recombination between the dominant strains, which can lead to the generation of new recombinants. Here, we identified a recombinant cluster composed of CRF02_AG, sub-subtype A3, and subtype G among HIV-infected children in Larkana. For the study, 10 retrospectively collected samples, with recombination signals in the pol gene, were used to perform a near full-length genome NFLG sequencing. Of the 10 samples, NFLG was successfully sequenced from seven samples. Phylogenetic analysis of the seven NFLGs showed that all recombinants formed a distinct monophyletic cluster and were distinct from known HIV-1 circulating recombinant forms CRFs. Recombination analyses showed that all seven NFLGs shared a similar recombinant structure consisting of CRF02_AG, sub-subtype A3, and subtype G, with a sub-subtype A3 fragment inserted into pol and vif regions spanning from (HXB2: 4218-5518), and a subtype G fragment inserted into vpu, rev, tat and env regions spanning from (HXB2: 5957-8250) of the CRF02_AG backbone. The identification of unique recombinant forms may indicate the presence and transmission of several co-circulating lineages in Larkana, giving rise to newer CRFs. This study also highlights the importance of continuous molecular surveillance to fully understand HIV-1 genetic diversity in Pakistan, particularly in Larkana, which is the epicenter of HIV outbreaks.

Genetic and antiretroviral drug resistance mutations analysis of reverse transcriptase and protease gene from Pakistani people living with HIV-1

BACKGROUND

Antiretroviral therapy (ART) effectiveness is compromised by the emergence of HIV drug resistance mutations (DRM) and can lead to the failure of ART. Apart from intrinsic viral factors, non-compliance with drugs and/or the use of sub-optimum therapy can lead to the emergence of DRMs. In Pakistan HIV currently exists as a concentrated epidemic, however, ART coverage is very low, and drug adherence is poor. ART is selected assuming without baseline genotyping. Pakistan has recently seen a rise in treatment failures, but the country's actual burden of DRM is still unknown. In this study, we perform the genetic and drug resistance analysis of the pol gene from Pakistani HIV-positive ART-naïve and ART-experienced individuals.

METHODS

In this study, HIV-1 pol was sequenced from 146 HIV-1 positive individuals, divided into ART-naïve (n = 37) and ART-experienced (n = 109). The sequences were also used to determine HIV-1 subtypes, the prevalence of DRM, and pol genetic variability.

RESULTS

DRM analysis identified numerous DRMs against reverse transcriptase inhibitors in both ART-naïve and ART-experienced groups, including a few that are classified as rare. Additionally, the ART-experienced group showed mutations associated with resistance to protease inhibitors. Genetic analysis showed negative selection pressure in both groups, but a higher rate of evolution in the ART-naïve group.

CONCLUSION

High prevalence of DRMs, especially against previous first-line treatment in ART- naïve and the accumulation of DRMs in ART-experienced groups is concerning and warrants that a more extensive DRM survey be carried out to inform first-line and second-line ART regimen recommendations.

Unassigned Complex Unique Recombinant Forms Related to CRF36_cpx in Children Identified in an HIV-1 Outbreak in Pakistan

In 2019, an outbreak of HIV infection predominantly affecting children occurred in Larkana district, Pakistan. This is the largest outbreak ever reported in this age group in Pakistan. In this study, we report two HIV-1 unique recombinant forms identified during the outbreak. Blood samples were collected from HIV-positive children as part of a case-control study to investigate the outbreak. The pol gene was sequenced and used to detect HIV subtype/recombinant forms using subtype, recombination, and phylogenetic analyses. Drug resistance mutation (DRM) analysis was performed to characterize the DRMs in each sequence. We observed the emergence of two unassigned unique recombinant forms related to CRF36_cpx in 15 individuals of the 344 samples collected. Genotype analysis revealed the presence of multiple DRMs associated with resistance to reverse transcriptase inhibitors. The discovery of these unassigned unique recombinant forms in our population highlights the need for comprehensive molecular epidemiological studies to fully understand the distribution and drug resistance patterns to aid control efforts.

A novel role for gag as a cis-acting element regulating RNA structure, dimerization and packaging in HIV-1 lentiviral vectors

Clinical usage of lentiviral vectors is now established and increasing but remains constrained by vector titer with RNA packaging being a limiting factor. Lentiviral vector RNA is packaged through specific recognition of the packaging signal on the RNA by the viral structural protein Gag. We investigated structurally informed modifications of the 5′ leader and gag RNA sequences in which the extended packaging signal lies, to attempt to enhance the packaging process by facilitating vector RNA dimerization, a process closely linked to packaging. We used in-gel SHAPE to study the structures of these mutants in an attempt to derive structure-function correlations that could inform optimized vector RNA design. In-gel SHAPE of both dimeric and monomeric species of RNA revealed a previously unreported direct interaction between the U5 region of the HIV-1 leader and the downstream gag sequences. Our data suggest a structural equilibrium exists in the dimeric viral RNA between a metastable structure that includes a U5–gag interaction and a more stable structure with a U5–AUG duplex. Our data provide clarification for the previously unexplained requirement for the 5′ region of gag in enhancing genomic RNA packaging and provide a basis for design of optimized HIV-1 based vectors.

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

Introduction

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.

Methods

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 pol 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.

Results

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.

Conclusion

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.

Streamlined Subpopulation, Subtype, and Recombination Analysis of HIV-1 Half-Genome Sequences Generated by High-Throughput Sequencing

The highly recombinogenic nature of human immunodeficiency virus type 1 (HIV-1) leads to recombination and emergence of quasispecies. It is important to reliably identify subpopulations to understand the complexity of a viral population for drug resistance surveillance and vaccine development. High-throughput sequencing (HTS) provides improved resolution over Sanger sequencing for the analysis of heterogeneous viral subpopulations. However, current methods of analysis of HTS reads are unable to fully address accurate population reconstruction. Hence, there is a dire need for a more sensitive, accurate, user-friendly, and cost-effective method to analyze viral quasispecies. For this purpose, we have improved the HIVE-hexahedron algorithm that we previously developed with in silico short sequences to analyze raw HTS short reads. The significance of this study is that our standalone algorithm enables a streamlined analysis of quasispecies, subtype, and recombination patterns from long HIV-1 genome regions without the need of additional sequence analysis tools. Distinct viral populations and recombination patterns identified by HIVE-hexahedron are further validated by comparison with sequences obtained by single genome sequencing (SGS).

High-throughput sequencing (HTS) has been widely used to characterize HIV-1 genome sequences. There are no algorithms currently that can directly determine genotype and quasispecies population using short HTS reads generated from long genome sequences without additional software. To establish a robust subpopulation, subtype, and recombination analysis workflow, we amplified the HIV-1 3′-half genome from plasma samples of 65 HIV-1-infected individuals and sequenced the entire amplicon (∼4,500 bp) by HTS. With direct analysis of raw reads using HIVE-hexahedron, we showed that 48% of samples harbored 2 to 13 subpopulations. We identified various subtypes (17 A1s, 4 Bs, 27 Cs, 6 CRF02_AGs, and 11 unique recombinant forms) and defined recombinant breakpoints of 10 recombinants. These results were validated with viral genome sequences generated by single genome sequencing (SGS) or the analysis of consensus sequence of the HTS reads. The HIVE-hexahedron workflow is more sensitive and accurate than just evaluating the consensus sequence and also more cost-effective than SGS.

IMPORTANCE The highly recombinogenic nature of human immunodeficiency virus type 1 (HIV-1) leads to recombination and emergence of quasispecies. It is important to reliably identify subpopulations to understand the complexity of a viral population for drug resistance surveillance and vaccine development. High-throughput sequencing (HTS) provides improved resolution over Sanger sequencing for the analysis of heterogeneous viral subpopulations. However, current methods of analysis of HTS reads are unable to fully address accurate population reconstruction. Hence, there is a dire need for a more sensitive, accurate, user-friendly, and cost-effective method to analyze viral quasispecies. For this purpose, we have improved the HIVE-hexahedron algorithm that we previously developed with in silico short sequences to analyze raw HTS short reads. The significance of this study is that our standalone algorithm enables a streamlined analysis of quasispecies, subtype, and recombination patterns from long HIV-1 genome regions without the need of additional sequence analysis tools. Distinct viral populations and recombination patterns identified by HIVE-hexahedron are further validated by comparison with sequences obtained by single genome sequencing (SGS).

Genetic Diversity of Human Immunodeficiency Virus Type 1 in Asymptomatic Blood Donors in Islamabad, Pakistan

Objective  The serological testing of human immunodeficiency virus (HIV) is mandatory under the blood safety legislation of Pakistan; hence, data exist on the prevalence of HIV in blood donors. However, little is known about the molecular epidemiology of HIV in the blood donor population. Therefore, the current study was designed to study the genetic diversity of HIV-1 infection in a population of apparently healthy treatment-naive blood donors in Islamabad, Pakistan.

Material and Methods  A total of 85,736 blood donors were tested for HIV by the chemiluminescence immunoassay. All positive donor samples were analyzed for the presence of various HIV genotypes (types and subtypes). Viral ribonucleic acid was extracted from blood samples of HIV positive donors and reverse transcribed into complementary deoxyribonucleic acid (cDNA). The cDNA of all positive donors was then analyzed for the presence of various HIV genotypes (types and subtypes) by employing subtype-specific primers in a nested polymerase chain reaction. The amplified products were run on ethidium bromide-stained 2% agarose gel and visualized using a ultraviolet transilluminator. A particular subtype was assigned to a sample if the subtype-specific reaction made a band 20% highly intense compared with the band made by the subtype-independent reaction.

Results  A total of 85,736 blood donors were screened for the presence of antibodies to HIV. Out of them, 114 were initially found reactive for HIV. The repeat testing resulted in 112 (0.13%) positive donors, 95% confidence interval 0.0014 (0.0011–0.0018). These 112 samples were analyzed for molecular typing of HIV-1. The predominant HIV-1 subtype was A ( n = 101) (90.1%) followed by subtype B ( n = 11) (9.9%).

Conclusion  These findings are key to understand the diversified HIV epidemic at the molecular level and should assist public health workers in implementing measures to lessen the further dissemination of these viruses in the country.

HIV-1 phylodynamic analysis among people who inject drugs in Pakistan correlates with trends in illicit opioid trade

Pakistan is considered by the World Health Organization to currently have a "concentrated" HIV-1 epidemic due to a rapid rise in infections among people who inject drugs (PWID). Prevalence among the country's nearly 105,000 PWID is estimated to be 37.8% but has been shown to be higher in several large urban centers. A lack of public health resources, the common use of professional injectors and unsafe injection practices are believed to have fueled the outbreak. Here we evaluate the molecular characteristics of HIV-1 sequences (n = 290) from PWID in several Pakistani cities to examine transmission dynamics and the association between rates of HIV-1 transmission with regards to regional trends in opioid trafficking. Tip-to-tip (patristic) distance based phylogenetic cluster inferences and BEAST2 Bayesian phylodynamic analyses of time-stamped data were performed on HIV-1 pol sequences generated from dried blood spots collected from 1,453 PWID as part of a cross-sectional survey conducted in Pakistan during 2014/2015. Overall, subtype A1 strains were dominant (75.2%) followed by CRF02_AG (14.1%), recombinants/unassigned (7.2%), CRF35_AD (2.1%), G (1.0%) and C (0.3%). Nearly three quarters of the PWID HIV-1 sequences belonged to one of five distinct phylogenetic clusters. Just below half (44.4%) of individuals in the largest cluster (n = 118) did seek help injecting from professional injectors which was previously identified as a strong correlate of HIV-1 infection. Spikes in estimated HIV-1 effective population sizes coincided with increases in opium poppy cultivation in Afghanistan, Pakistan's western neighbor. Structured coalescent analysis was undertaken in order to investigate the spatial relationship of HIV-1 transmission among the various cities under study. In general terms, our analysis placed the city of Larkana at the center of the PWID HIV-1 epidemic in Pakistan which is consistent with previous epidemiological data.

Validation of Variant Assembly Using HAPHPIPE with Next-Generation Sequence Data from Viruses

Next-generation sequencing (NGS) offers a powerful opportunity to identify low-abundance, intra-host viral sequence variants, yet the focus of many bioinformatic tools on consensus sequence construction has precluded a thorough analysis of intra-host diversity. To take full advantage of the resolution of NGS data, we developed HAplotype PHylodynamics PIPEline (HAPHPIPE), an open-source tool for the de novo and reference-based assembly of viral NGS data, with both consensus sequence assembly and a focus on the quantification of intra-host variation through haplotype reconstruction. We validate and compare the consensus sequence assembly methods of HAPHPIPE to those of two alternative software packages, HyDRA and Geneious, using simulated HIV and empirical HIV, HCV, and SARS-CoV-2 datasets. Our validation methods included read mapping, genetic distance, and genetic diversity metrics. In simulated NGS data, HAPHPIPE generated pol consensus sequences significantly closer to the true consensus sequence than those produced by HyDRA and Geneious and performed comparably to Geneious for HIV gp120 sequences. Furthermore, using empirical data from multiple viruses, we demonstrate that HAPHPIPE can analyze larger sequence datasets due to its greater computational speed. Therefore, we contend that HAPHPIPE provides a more user-friendly platform for users with and without bioinformatics experience to implement current best practices for viral NGS assembly than other currently available options.

Circulation of multiple subtypes (A, G and CRFs 02_AG) of human immunodeficiency virus type 1 (HIV-1) in selected districts of Punjab province, Pakistan

Owing to consistent genetic mutation and recombination, various escape mutants and/or drug-resistant mutants of human immunodeficiency virus (HIV-1) are now emerging worldwide. Therefore, an understanding of the genetic characteristics of prevailing strains, particularly with regard to drug-resistance-associated substitutions, is essential for devising and implementing treatments and disease control interventions in endemic settings such as Pakistan. We processed a total of 130 plasma samples originating from HIV-treatment centers in selected districts of Punjab province, Pakistan. The samples were first screened using an HIV-1 Ag/Ab Combo test followed by amplification of the pol gene (1084 bp) from samples that were positive either for the antigen or for both the antigen and antibodies simultaneously. Screening revealed that a total of 45 samples were positive (34.62%; 95% CI: 26.99-43.13) for either antigen or both antigen and antibodies (n = 18, 40%; 95% CI: 27.02-54.55) or for antibodies alone (n = 27, 60%; 95% CI: 45.45-72.98). A largest number of positive samples was from the district of Lahore (n = 19/43, 44.18%; 95% CI: 30.44-58.9) followed by Faisalabad (n= 12/36, 33.33%; 95% CI: 20.21-49.66), Gujranwala (n = 05/23, 21.7%; 95% CI: 9.66-41.9) and Sargodha (n = 09/28, 32.1%; 95% CI: 17.93-50.66). The probability of occurrence of HIV infection was significantly associated with individuals having a history of injecting drug use (68.08%; OR = 11.15; 95% CI: 53.84-79.61, p = 0.0001). Phylogenetic analysis based on the pol gene showed that the sequences from this study clustered into three distinct clades representing recombinant form 02_AG (n = 14, 77.0%; 95% CI: 54.79-91.00), and subtypes A (n = 2, 11.1%; 95% CI: 3.1-32.8) and G (n = 2, 11.1%; 95% CI: 3.1-32.8). Although we screened 18 samples for drug-resistance-associated mutations, except for an accessory mutation (M46K) in the protease (PR) region in one subject, we found a lack of drug-resistance-associated substitutions in the PR region. On the other hand, we found two subjects (2/18) carrying a resistance-associated mutation (V106I) conferring a low level of resistance against non-nucleoside reverse transcriptase inhibitors. The present study shows that multiple subtypes of HIV-1 are present in the affected population. Continuous disease surveillance coupled with evaluation of drug resistance at higher resolution should be done in future studies.

Increased predominance of HIV-1 CRF01_AE and its recombinants in the Philippines

The growth rate of new HIV infections in the Philippines was the fastest of any countries in the Asia-Pacific region between 2010 and 2016. To date, HIV-1 subtyping results in the Philippines have been determined by characterizing only partial viral genome sequences. It is not known whether recombination occurs in the majority of unsequenced genome regions. Near-full-length genome (NFLG) sequences were obtained by amplifying two overlapping half genomes from plasma samples collected between 2015 and 2017 from 23 newly diagnosed infected individuals in the Philippines. Phylogenetic analysis showed that the newly characterized sequences were CRF01_AE (14), subtype B (3), CRF01/B recombinants (5) and a CRF01/CRF07/B recombinant (1). All 14 CRF01_AE formed a tight cluster, suggesting that they were derived from a single introduction. The time to the most recent common ancestor (tMRCA) for CRF01_AE in the Philippines was 1995 (1992-1998), about 10-15 years later than that of CRF01_AE in China and Thailand. All five CRF01/B recombinants showed distinct recombination patterns, suggesting ongoing recombination between the two predominant circulating viruses. The identification of partial CRF07_BC sequences in one CRF01/CRF07/B recombinant, not reported previously in the Philippines, indicated that CRF07_BC may have been recently introduced into that country from China, where CRF07_BC is prevalent. Our results show that the major epidemic strains may have shifted to an increased predominance of CRF01_AE and its recombinants, and that other genotypes such as CRF07_BC may have been introduced into the Philippines.

Late diagnosis of human immunodeficiency virus infections in high-risk groups in Karachi, Pakistan

Human immunodeficiency virus (HIV) infection prevalence in Pakistan has been increasing in high-risk groups, including people who inject drugs (PWID) and transgender hijra sex workers (TG-HSWs) nationwide. Effective control of HIV requires early diagnosis of the infection. We investigated recency of HIV infections in newly-diagnosed cases in PWID and TG-HSWs. This was an observational study with convenience sampling. Overall, 210 HIV-positive subjects comprising an equal number of PWID and TG-HSWs were included. Antibody avidity was tested using the Maxim HIV-1 Limiting Antigen Avidity (LAg) EIA (Maxim Biomedical, Inc. Rockville, Maryland, USA). The mean age of study subjects was 29.5 years: PWID, 28.5 years and TG-HSWs, 30.4 years. Study subjects were married, 27%, or unmarried. Eighteen percent of individuals had recently-acquired HIV infections: 19% of PWID and 17% of TG-HSWs. Eighty-two percent of individuals had long-term HIV infections: 81% of PWID and 83% of TG-HSWs. This is the first study identification of recent HIV-1 infections in Pakistan. We show that most newly-diagnosed HIV patients in the high-risk groups studied had long-term infections. There is an urgent need for intervention in these groups to facilitate early diagnosis and treatment of HIV infection to reduce transmission in Pakistan.