Human Immunodeficiency Virus-1 Diversity in the Moscow Region, Russia: Phylodynamics of the Most Common Subtypes

Features of Tat Protein in HIV-1 Sub-Subtype A6 Variants Circulating in the Moscow Region, Russia

Tat, the trans-activator of transcription, is a multifunctional HIV-1 protein that can induce chronic inflammation and the development of somatic diseases in HIV-infected patients. Natural polymorphisms in Tat can impact the propagation of the inflammatory signal. Currently, Tat is considered an object for creating new therapeutic agents. Therefore, the identification of Tat protein features in various HIV-1 variants is a relevant task. The purpose of the study was to characterize the genetic variations of Tat-A6 in virus variants circulating in the Moscow Region. The authors analyzed 252 clinical samples from people living with HIV (PLWH) with different stages of HIV infection. Nested PCR for two fragments (tat1, tat2) with subsequent sequencing, subtyping, and statistical analysis was conducted. The authors received 252 sequences for tat1 and 189 for tat2. HIV-1 sub-subtype A6 was identified in 250 samples. The received results indicated the features of Tat1-A6 in variants of viruses circulating in the Moscow Region. In PLWH with different stages of HIV infection, C31S in Tat1-A6 was detected with different occurrence rates. It was demonstrated that Tat2-A6, instead of a functional significant 78RGD80 motif, had a 78QRD80 motif. Herewith, G79R in Tat2-A6 was defined as characteristic amino acid substitution for sub-subtype A6. Tat2-A6 in variants of viruses circulating in the Moscow Region demonstrated high conservatism.

Female Sex Hormones Upregulate the Replication Activity of HIV-1 Sub-Subtype A6 and CRF02_AG but Not HIV-1 Subtype B

More than 50% of all people living with HIV worldwide are women. Globally, HIV/AIDS is the leading cause of death among women aged 15 to 44. The safe and effective methods of hormonal contraception are an essential component of preventive medical care in order to reduce maternal and infant mortality. However, there is limited knowledge regarding the effect of hormones on the rate of viral replication in HIV infection, especially non-B subtypes. The goal of the present work was to study in vitro how the female hormones β-estradiol and progesterone affect the replication of the HIV-1 subtypes A6, CRF02_AG, and B. The findings show that high doses of hormones enhanced the replication of HIV-1 sub-subtype A6 by an average of 1.75 times and the recombinant variant CRF02_AG by 1.4 times but did not affect the replication of HIV-1 subtype B. No difference was detected in the expression of CCR5 and CXCR4 co-receptors on the cell surface, either in the presence or absence of hormones. However, one of the reasons for the increased viral replication could be the modulated TLRs secretion, as it was found that high doses of estradiol and progesterone upregulated, to varying degrees, the expression of TLR2 and TLR9 genes in the PBMCs of female donors infected with HIV-1 sub-subtype A6.

[Genetic diversity of capsid protein (p24) in human immunodeficiency virus type-1 (HIV-1) variants circulating in the Russian Federation]

INTRODUCTION

The human immunodeficiency virus (HIV) protein p24 plays an important role in the life cycle of the virus, and also is a target for diagnostic tests and for new antiretroviral drugs and therapeutic vaccines. The most studied variant of HIV-1 in the world is subtype B. In Russia, the most common variant is A6, the spread of recombinant forms (CRF63_02A6, CRF03_A6B) is observed as well as circulation of G and CRF02_AG variants. However, a detailed study of the p24 protein in these variants has not yet been conducted. The aim was to study the features of the p24 protein in HIV-1 variants circulating in Russia and estimate the frequency of occurrence of pre-existing mutations associated with resistance to lenacapavir, the first antiretroviral drug in the class of capsid inhibitors.

MATERIALS AND METHODS

The objects of the study were the nucleotide sequences obtained from the Los Alamos international database and clinical samples from HIV infected patients.

RESULTS AND DISCUSSION

The features of HIV-1 variants circulating in Russia have been determined. V86A, H87Q, I91F are characteristic substitutions in A6 genome. It is shown that the presence of preexisting mutations associated with resistance to lenacapavir is unlikely.

CONCLUSION

Features of the p24 protein in HIV-1 variants circulating in Russia allow them to be distinguished from others variants and among themselves. The prognosis for the use of lenacapavir in Russia is generally favorable. The results obtained could be taken into account in developing and using antiretroviral drugs and therapeutic vaccines.

Reconstructing the Temporal Origin and the Transmission Dynamics of the HIV Subtype B Epidemic in St. Petersburg, Russia

The HIV/AIDS epidemic in Russia is among the fastest growing in the world. HIV epidemic burden is non-uniform in different Russian regions and diverse key populations. An explosive epidemic has been documented among people who inject drugs (PWID) starting from the mid-1990s, whereas presently, the majority of new infections are linked to sexual transmission. Nationwide, HIV sub-subtype A6 (previously called A_FSU) predominates, with the increasing presence of other subtypes, namely subtype B and CRF063_02A. This study explores HIV subtype B sequences from St. Petersburg, collected from 2006 to 2020, in order to phylogenetically investigate and characterize transmission clusters, focusing on their evolutionary dynamics and potential for further growth, along with a socio-demographic analysis of the available metadata. In total, 54% (107/198) of analyzed subtype B sequences were found grouped in 17 clusters, with four transmission clusters with the number of sequences above 10. Using Bayesian MCMC inference, tMRCA of HIV-1 subtype B was estimated to be around 1986 (95% HPD 1984-1991), whereas the estimated temporal origin for the four large clusters was found to be more recent, between 2001 and 2005. The results of our study imply a complex pattern of the epidemic spread of HIV subtype B in St. Petersburg, Russia, still in the exponential growth phase, and in connection to the men who have sex with men (MSM) transmission, providing a useful insight needed for the design of public health priorities and interventions.

Romania, a Harbour of HIV-1 Subtype F1: Where Are We after 33 Years of HIV-1 Infection?

Infection with the human immunodeficiency virus (HIV) has been a major public health concern worldwide for more than 30 years, including in Romania. The F1 HIV-1 subtype was exported from Angola to Romania most probably because of the two countries' close political connections. Patients infected with HIV-1 via re-used and improperly sterilized injection equipment and through transfusions of unscreened blood, also known as the "Romanian cohort", were the most common type of HIV-1 infection in Romania in the early 1990s, when the virus's presence was recognized. Recently, subtype B started to increase in our country, mostly diagnosed in people using intravenous drugs or in men having sex with men. The evolution of the HIV-1 infection in Romania has been unique, with a dominance of the subtype F1, making it different from other countries in Europe.

Spatiotemporal dynamics of HIV-1 CRF63_02A6 sub-epidemic

HIV-1 epidemic in Russia is one of the fastest growing in the world reaching 1.14 million people living with HIV-1 (PLWH) in 2021. Since mid-1990s, the HIV-1 epidemic in Russia has started to grow substantially due to the multiple HIV-1 outbreaks among persons who inject drugs (PWID) leading to expansion of the HIV-1 sub-subtype A6 (former Soviet Union (FSU) subtype A). In 2006, a local HIV-1 sub-epidemic caused by the distribution of novel genetic lineage CRF63_02A6 was identified in Siberia. In this study, we used a comprehensive dataset of CRF63_02A6 pol gene sequences to investigate the spatiotemporal dynamic of the HIV-1 CRF63_02A6 sub-epidemic. This study includes all the available CRF63_02A6 HIV-1 pol gene sequences from Los Alamos National Laboratory (LANL) HIV Sequence Database. The HIV-1 subtypes of those sequences were conferred using phylogenetic analysis, and two automated HIV-1 subtyping tools Stanford HIVdb Program and COMET. Ancestral state reconstruction and origin date were estimated using Nextstrain. Evolutionary rate and phylodynamic analysis were estimated using BEAST v 1.10.4. CRF63_02A6 was assigned for 872 pol gene sequences using phylogenetic analysis approach. Predominant number (n = 832; 95.4%) of those sequences were from Russia; the remaining 40 (4.6%) sequences were from countries of Central Asia. Out of 872 CRF63_02A6 sequences, the corresponding genetic variant was assigned for 75.7 and 79.8% of sequences by Stanford and COMET subtyping tools, respectively. Dated phylogenetic analysis of the CRF63_02A6 sequences showed that the virus most likely originated in Novosibirsk, Russia, in 2005. Over the last two decades CRF63_02A6 has been widely distributed across Russia and has been sporadically detected in countries of Central Asia. Introduction of new genetic variant into mature sub-subtype A6 and CRF02_AGFSU epidemics could promote the increase of viral genetic diversity and emergence of new recombinant forms. Further HIV-1 studies are needed due to a continuing rapid virus distribution. Also, the implementation of HIV-1 prevention programs is required to reduce HIV-1 transmission. This study also highlights the discrepancies in HIV-1 subtyping approaches. The reference lists of HIV-1 sequences implemented in widely used HIV-1 automated subtyping tools need to be updated to provide reliable results.

Identifying HIV-1 Transmission Clusters in Uzbekistan through Analysis of Molecular Surveillance Data

The CRF02_AG and sub-subtype A6 are currently the predominant HIV-1 variants in the Republic of Uzbekistan, but little is known about their time-spatial clustering patterns in high-risk populations. We have applied molecular evolution methods and network analyses to better understand the transmission patterns of these subtypes by analyzing 316 pol sequences obtained during the surveillance study of HIV drug resistance. Network analysis showed that about one third of the HIV infected persons were organized into clusters, including large clusters with more than 35 members. These clusters were composed mostly of injecting drug users and/or heterosexuals, with women having mainly high centrality within networks identified in both subtypes. Phylogenetic analyses of the 'Uzbek' sequences, including those publicly available, show that Russia and Ukraine played a role as the main sources of the current subtype A6 epidemic in the Republic. At the same time, Uzbekistan has been a local center of the CRF02_AG epidemic spread in the former USSR since the early 2000s. Both of these HIV-1 variants continue to spread in Uzbekistan, highlighting the importance of identifying transmission networks and transmission clusters to prevent further HIV spread, and the need for HIV prevention and education campaigns in high-risk groups.

Molecular Epidemiology of HIV-1 in Oryol Oblast, Russia

The HIV/AIDS epidemic in Russia is growing, with approximately 100,000 people infected annually. Molecular epidemiology can provide insight into the structure and dynamics of the epidemic. However, its applicability in Russia is limited by the weakness of genetic surveillance, as viral genetic data are only available for <1 per cent of cases. Here, we provide a detailed description of the HIV-1 epidemic for one geographic region of Russia, Oryol Oblast, by collecting and sequencing viral samples from about a third of its known HIV-positive population (768 out of 2,157 patients). We identify multiple introductions of HIV-1 into Oryol Oblast, resulting in eighty-two transmission lineages that together comprise 66 per cent of the samples. Most introductions are of subtype A (315/332), the predominant HIV-1 subtype in Russia, followed by CRF63 and subtype B. Bayesian analysis estimates the effective reproduction number R_e for subtype A at 2.8 [1.7–4.4], in line with a growing epidemic. The frequency of CRF63 has been growing more rapidly, with the median R_e of 11.8 [4.6–28.7], in agreement with recent reports of this variant rising in frequency in some regions of Russia. In contrast to the patterns described previously in European and North American countries, we see no overrepresentation of males in transmission lineages; meanwhile, injecting drug users are overrepresented in transmission lineages. This likely reflects the structure of the HIV-1 epidemic in Russia dominated by heterosexual and, to a smaller extent, people who inject drugs transmission. Samples attributed to men who have sex with men (MSM) transmission are associated with subtype B and are less prevalent than expected from the male-to-female ratio for this subtype, suggesting underreporting of the MSM transmission route. Together, our results provide a high-resolution description of the HIV-1 epidemic in Oryol Oblast, Russia, characterized by frequent interregional transmission, rapid growth of the epidemic, and rapid displacement of subtype A with the recombinant CRF63 variant.

Ongoing HIV transmission following a large outbreak among people who inject drugs in Athens, Greece (2014-20)

BACKGROUND AND AIMS

The human immunodeficiency virus (HIV) outbreak among people who inject drugs (PWID) in Athens, Greece in 2011-13 was the largest recent epidemic in Europe and North America. We aimed to assess trends in HIV prevalence, drug use and access to prevention among PWID in Athens to estimate HIV incidence and identify risk factors and to explore HIV-1 dispersal using molecular methods during 2014-20.

METHODS

Two community-based HIV/hepatitis C programmes on PWID were implemented in 2012-13 (n = 3320) and 2018-20 (n = 1635) through consecutive respondent-driven sampling (RDS) rounds. PWID were uniquely identified among rounds/programmes. We obtained RDS-weighted HIV prevalence estimates per round for 2018-20 and compared them to 2012-13. We assessed changes in HIV status, behaviours and access to prevention in PWID participating in both periods. We estimated HIV incidence in a cohort of seronegative PWID as the number of HIV seroconversions/100 person-years during 2014-20 and used Cox regression to identify associated risk factors. Molecular sequencing and phylogenetic analysis were performed in HIV seroconverters.

RESULTS

HIV prevalence per round ranged between 12.0 and 16.2% in 2012-13 and 10.7 and 11.3% in 2018-20 with overlapping 95% confidence intervals (95% CI). Among PWID participating in both programmes, HIV prevalence (95% CI) increased from 14.2% (11.7-17.1%) in 2012-13 to 22.0% (19.0-25.3%) in 2018-20 (P < 0.001). There was a deterioration in socio-economic characteristics such as homelessness [from 16.2% (95% CI = 13.5-19.2%) to 25.6% (22.3-29.0%)], a shift in cocaine use [16.6% (13.9-19.6%) versus 28.1% (24.7-31.7%], reduced access to free syringes [51.8% (48.0-55.7%) versus 44.5% (40.7-48.3%)] and a decrease in daily injecting [36.2% (32.6-39.9%) versus 28.5% (25.2-32.1%)]. HIV incidence (95% CI) in 2014-20 was 1.94 (1.50-2.52) new cases/100 person-years and younger age, lower educational level, larger injection network and daily injecting were risk factors. Almost 9% of HIV seroconversions occurred within a newly expanding phylogenetic cluster.

CONCLUSIONS

In Athens, Greece, compared with the period 2012-13, in the period 2018-20 there was a deterioration in socio-economic conditions among people who inject drugs, an increase in the use of cocaine, reduced access to needle and syringe programmes and stable low levels of human immunodeficiency virus testing. Ongoing human immunodeficiency virus transmission was documented during 2014-20 in existing as well as new transmission clusters.

Construction and Characterization of HIV-1 env-Pseudoviruses of the Recombinant Form CRF63_02A and Subtype A6

HIV-1 env-pseudoviruses are a useful tool in the search for antiviral drugs (entry inhibitors) and evaluation of the efficacy of HIV-1 vaccines. Given the high genetic variability of HIV-1, it is necessary to regularly update the panels of pseudoviruses in accordance with the emergence of new strains. Based on genetic variants of HIV-1 circulating in the regions of the Siberian Federal District, 13 HIV-1 env-pseudoviruses of recombinant form CRF63_02A and subtype A6 were obtained. Most pseudoviruses have been shown to be sensitive to neutralization by bnAbs VRC01, PGT126, and 10E8, moderately sensitive to bnAbs PG9 and 4E10, and resistant to bnAbs 2G12, PG16, and 2F5. All obtained variants of pseudoviruses are CCR5-tropic.

[Analysis of Tat protein characteristics in human immunodeficiency virus type 1 sub-subtype A6 (Retroviridae: Orthoretrovirinae: Lentivirus: Human immunodeficiency virus-1)]

INTRODUCTION

Tat protein is a major factor of HIV (human immunodeficiency virus) transcription regulation and has other activities. Tat is characterized by high variability, with some amino acid substitutions, including subtypespecific ones, being able to influence on its functionality. HIV type 1 (HIV-1) sub-subtype A6 is the most widespread in Russia. Previous studies of the polymorphisms in structural regions of the A6 variant have shown numerous characteristic features; however, Tat polymorphism in A6 has not been studied.Goals and tasks. The main goal of the work was to analyze the characteristics of Tat protein in HIV-1 A6 variant, that is, to identify substitutions characteristic for A6 and A1 variants, as well as to compare the frequency of mutations in functionally significant domains in sub-subtype A6 and subtype B.

MATERIAL AND METHODS

The nucleotide sequences of HIV-1 sub-subtypes A6, A1, A2, A3, A4, subtype B and the reference nucleotide sequence were obtained from the Los Alamos international database.

RESULTS AND DISCUSSION

Q54H and Q60H were identified as characteristic substitutions. Essential differences in natural polymorphisms between sub-subtypes A6 and A1 have been demonstrated. In the CPP-region, there were detected mutations (R53K, Q54H, Q54P, R57G) which were more common in sub-subtype A6 than in subtype B.

CONCLUSION

Tat protein of sub-subtype A6 have some characteristics that make it possible to reliably distinguish it from other HIV-1 variants. Mutations identified in the CPP region could potentially alter the activity of Tat. The data obtained could form the basis for the drugs and vaccines development.

Genetic Diversity and Drug Resistance Mutations in Reverse Transcriptase and Protease Genes of HIV-1 Isolates from Southwestern Siberia

The analysis of a pol gene fragment encoding protease and part of reverse transcriptase was carried out for 55 sera collected in 2016 and 2018 from HIV-1-infected patients diagnosed in 2014-2018 living in the south of Western Siberia, Russia: Altai Territory (n = 11), Republic of Altai (n = 15), Kemerovo region (n = 18), and Novosibirsk region (n = 11). CRF63_02A was the dominant genetic form (>70%) in the Altai Territory and Kemerovo and Novosibirsk regions, with subsubtype A6 comprising <30% of samples. In the Altai Republic, subsubtype A6 was predominant (53%), with 33% of viruses belonging to CRF63_02A. Four CRF63_02A/A6 unique recombinant forms were identified in the Altai Territory, Kemerovo Region, and the Altai Republic. A majority (11 of 15) of CRF63_02A viruses from Kemerovo were grouped in a cluster. Antiretroviral (ARV) drug resistance mutations were found in 6 (14%) of 43 drug-naive patients. This study provides new insights in HIV-1 molecular epidemiology and prevalence of transmitted ARV drug resistance mutations in Southwestern Siberia.

High-throughput sequencing reveals a high prevalence of pretreatment HIV-1 drug resistance in Sweden

OBJECTIVES

HIV-1 pretreatment drug resistance (PDR) is a global concern. Our aim was to evaluate high-throughput sequencing (HTS) for HIV-1 resistance testing and describe PDR in Sweden, where 75% of diagnosed individuals are foreign-born.

DESIGN

Cross-sectional study.

METHODS

Individuals entering HIV-1 care in Sweden 2017 to March 2019 (n = 400) were included if a viremic sample was available (n = 220). HTS was performed using an in-house assay. Drug resistance mutations (DRMs) (based on Stanford HIV DB vs. 8.7) at levels 1-5%, 5-19% and at least 20% of the viral population were described. Results from HTS and routine Sanger sequencing were compared.

RESULTS

HTS was successful in 88% of patients, 92% when viral load was at least 1000 copies/ml. DRMs at any level in protease and/or reverse transcriptase were detected in 95 individuals (49%), whereas DRMs at least 20% in 35 (18%) individuals. DRMs at least 20% correlated well to findings in routine Sanger sequencing. Protease/reverse transcriptase (PR/RT) DRMs at least 20% were predicted by treatment exposure; adjusted OR 9.28 (95% CI 2.24-38.43; P = 0.002) and origin in Asia; adjusted OR 20.65 (95% CI 1.66-256.24; P = 0.02). Nonnucleoside reverse transcriptase inhibitor (NNRTI) DRMs at least 20% were common (16%) and over-represented in individuals originating from sub-Saharan Africa or Asia. Low-level integrase strand transfer inhibitor (INSTI) DRMs less than 20% were detected in 15 individuals (8%) with no association with INSTI exposure.

CONCLUSION

Our HTS can efficiently detect PDR and findings of DRMs at least 20% compare well to routine Sanger sequencing. The high prevalence of PDR was because of NNRTI DRMs and associated with migration from areas with emerging PDR.

Genetic Features of HIV-1 Integrase Sub-Subtype A6 Predominant in Russia and Predicted Susceptibility to INSTIs

The increasing use of the integrase strand transfer inhibitor (INSTI) class for the treatment of HIV-infection has pointed to the importance of analyzing the features of HIV-1 subtypes for an improved understanding of viral genetic variability in the occurrence of drug resistance (DR). In this study, we have described the prevalence of INSTI DR in a Russian cohort and the genetic features of HIV-1 integrase sub-subtype A6. We included 408 HIV infected patients who were not exposed to INSTI. Drug resistance mutations (DRMs) were detected among 1.3% of ART-naïve patients and among 2.7% of INSTI-naïve patients. The prevalence of 12 polymorphic mutations was significantly different between sub-subtypes A6 and A1. Analysis of the genetic barriers determined two positions in which subtype A (A1 and A6) showed a higher genetic barrier (G140C and V151I) compared with subtype B, and one position in which subtypes A1 and B displayed a higher genetic barrier (L74M and L74I) than sub-subtype A6. Additionally, we confirmed that the L74I mutation was selected at the early stage of the epidemic and subsequently spread as a founder effect in Russia. Our data have added to the overall understanding of the genetic features of sub-subtype A6 in the context of drug resistance.

Molecular Surveillance of HIV-1 Infection in Krasnoyarsk Region, Russia: Epidemiology, Phylodynamics and Phylogeography

BACKGROUND

The information about the dynamics of the viral population and migration events that affect the epidemic in different parts of the Russia is insufficient. Possibly, the huge size of the country and limited transport accessibility to certain territories may determine unique traits of the HIV-1 evolutionary history in different regions.

OBJECTIVE

The aim of this study was to explore the genetic diversity of HIV-1 in the Krasnoyarsk region and reconstruct spatial-temporal dynamics of the infection in the region.

METHODS

The demographic and virologic data from 281 HIV-infected individuals in Krasnoyarsk region collected during 2011-2016 were analyzed. The time to the most recent common ancestor, evolutionary rates, population growth, and ancestral geographic movements was estimated using Bayesian coalescent-based methods.

RESULTS

The study revealed moderate diversity of the HIV-1 subtypes found in the region, which included A6 (92.3%), CRF063_02A (4.3%), B (1.1%), and unique recombinants (2.5%). Phylogenetic reconstruction revealed that the A6 subtype was introduced into Krasnoyarsk region by one viral lineage, which arose around 1996.9 (1994.5-1999.5). The phylogeography analysis pointed to Krasnoyarsk city as the geographical center of the epidemic, which further spread to central neighboring districts of the region. At least two epidemic growth phases of subtype A6 were identified which included exponential growth in early-2000s followed by the decline in the mid/late 2010s.

CONCLUSION

This study demonstrates a change in the genetic diversity of HIV-1 in the Krasnoyarsk region. At the beginning of the epidemic, subtype A6 prevailed, subtypes B and CRF063_02A appeared in the region later.

HIV-1 Unique Recombinant Forms Identified in Slovenia and Their Characterization by Near Full-Length Genome Sequencing

Surveillance of HIV circulating recombinant forms (CRFs) is important because HIV diversity can affect various aspects of HIV infection from prevention to diagnosis and patient management. A comprehensive collection of pol sequences obtained from individuals diagnosed with HIV-1 from 2000 to 2016 in Slovenia was subtyped to identify possible unique recombinant forms (URFs). Selected samples were subjected to near full-length genome (NFLG) sequencing and detailed recombination analyses. Discordant subtyping results were observed for 68/387 (17.6%) sequences and 20 sequences were identified as the most probable URFs and selected for NFLG characterization. Further, 11 NFLGs and two sequences of >7000 base pairs were obtained. Seven sequences were identified as “pure” subtypes or already characterized CRFs: subtype B (n = 5), sub-subtype A6 (n = 1), and CRF01_AE (n = 1). The remaining six sequences were determined to be URFs; four displayed a single recombination event and two exhibited a complex recombination pattern involving several subtypes or CRFs. Finally, three HIV strains were recognized as having epidemic potential and could be further characterized as new CRFs. Our study shows that the identification of new CRFs is possible, even in countries where HIV diversity is considered limited, emphasizing the importance of the surveillance of HIV recombinant forms.

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.