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.

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.

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.

Phylogenetic Characterization of HIV-1 Sub-Subtype A1 in Karachi, Pakistan.

(1) Background: HIV-1 sub-subtype A1 is common in parts of Africa, Russia, former Soviet Union countries, and Eastern Europe. In Pakistan, sub-subtype A1 is the predominant HIV-1 subtype. Preliminary evidence suggests that distinct strains of HIV-1 sub-subtype A1 are circulating in Pakistan; however, an in-depth molecular phylogenetic characterization of HIV-1 sub-subtype A1 strains in Pakistan have not been presented. We performed a detailed characterization of the HIV-1 sub-subtype A1 epidemic in Pakistan using state-of-the-art molecular epidemiology and phylodynamics. (2) Methods: A total of 143 HIV-1 sub-subtype A1 gag sequences, including 61 sequences generated specifically for this study from PLHIVs part of our cohort, representing all sub-subtype A1 gag sequences from Pakistan, were analyzed. Maximum-likelihood phylogenetic cluster analysis was used to determine the relationship between Pakistani sub-subtype A1 strains and pandemic sub-subtype A1 strains. Furthermore, we used signature variation, charge distribution, selection pressures, and epitope prediction analyses to characterize variations unique to Pakistani HIV-1 strains and establish the association between signature variations and Gag epitope profile. (3) Results: The HIV-1 sub-subtype A1 sequences from Pakistan formed three main clusters: two that clustered with Kenyan sequences (7 and 10 sequences, respectively) and one that formed a Pakistan-specific cluster of 123 sequences that were much less related to other sub-subtype A1 sequences available in the database. The sequences in the Pakistan-specific cluster and the Kenyan reference strains exhibited several signature variations, especially at amino acid positions 312, 319, 331, 372, 373, 383, and 402. Structural protein modeling suggested that amino acid changes in these positions result in alterations of the Gag protein structure as well as in Gag-specific T-cell epitopes. (4) Conclusions: Our results suggest that the majority of the Pakistan HIV-1 sub-subtype A1 strains were unique to Pakistan and with a specific mutation pattern in Gag.

Phylogenetic and drug- and vaccine-resistance profiles of Hepatitis B Virus among children with HIV co-infection in Pakistan.

INTRODUCTION: HIV-1 and hepatitis B virus (HBV) share common routes of transmission and therefore co-infection is common. In 2019, an HIV-1 outbreak that resulted in >1000 children being infected, predominantly through nosocomial transmission, occurred in Sindh, Pakistan. We conducted a phylogenetic and drug resistance analysis of the HBV Reverse Transcriptase (RT) gene in children with HIV-1 and HBV co-infection. METHODOLOGY: Blood samples were collected from 321 children with HIV who were recruited as part of a study to investigate the HIV-1 outbreak. All samples were tested for HBV surface antigen (HBsAg) using an ELISA assay, and positive samples were used to amplify and sequence the HBV RT gene. The phylogenetic relationship between sequences was analyzed, and drug- and vaccine- resistance mutations in the RT gene were explored. RESULTS: Of 321 samples, 23% (n = 75) were positive for HBsAg on ELISA. Phylogenetic analysis of the sequences revealed that 63.5% of HBV sequences were sub-genotype D1, while the rest were sub-genotype D2. Cluster analysis revealed grouping of sub-genotype D1 sequences exclusively with Pakistani sequences, while clustering of sub-genotypes D2 predominantly with global sequences. The 236Y mutation associated with resistance to tenofovir was observed in 2.8% of HBV sequences. Additionally, seven vaccine escape mutations were observed, the most common being 128 V. CONCLUSION: Our study suggests ongoing transmission of HBV D1 and D2 sub-genotypes in the HIV-1 co-infected population, likely nosocomially, given common routes of HVB and HIV-1 transmission. The prevalence of major HBV drug- and vaccine-resistant mutations remains low. Surveillance for further transmissions and the possible emergence of major drug- or vaccine-resistant variants is required.

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.

Viral infections in Pakistan: prevalence, factors affecting spread, and recommendations for control.

Pakistan is endemic to a number of viral infections, owing to its humid climate, topographical variation, soaring population, and lack of education and awareness. These viruses may have several different modes of transmission, including respiratory or airborne transmission, sexual transmission, blood-borne, fecal-oral transmission, vector-borne transmission, and transmission following an organ transplant. Although several different microorganisms are responsible for causing these infections, a few viruses are found more commonly in Pakistan and are primarily responsible for causing infections. In this study, we present a review of the most recent studies on different viruses, transmitted through various transmission routes, found commonly in Pakistan, along with the prevalence of each, and recommend control measures required against these viruses.

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.

Probabilistic measures of HIV-1 transmission in different HIV-1 key population groups of Larkana, Pakistan.

OBJECTIVE: To estimate the probability of human immunodeficiency virus (HIV)-1 transmission from different key HIV population groups using probabilistic modelling. METHODS: This study was conducted in December 2020. A probabilistic model was used to estimate the probability of HIV-1 transmission from different key HIV population groups in Larkana. Our model was run on three probabilistic assumptions: 1) each replication gave two conceivable results: 'true' or 'false'; 2) the chance of giving a 'true' result is the same for each replication; and 3) the replications are independent - 'true' in one will not impact the likelihood of 'true' in another. RESULTS: The results estimated the probability of HIV transmission in key HIV population groups in Larkana to range between 0.42-0.54 per trial, where the highest probability of transmission was predicted for men who have sex with men (MSM; 0.54 per trial), followed by transgender (TG; 0.46 per trial) and people who inject drugs (PWID; 0.457 per trial). CONCLUSIONS: Our results suggest that there is a high likelihood of HIV transmission by key population groups in Larkana, such as MSM, TG, and PWID. Mathematic models, such as one proposed in our study can aid the HIV and acquired immunodeficiency syndrome (AIDS) control programmes in evaluating and optimising the strategies in controlling transmission of HIV from the key population groups.

Factors associated with HIV infection among children in Larkana District, Pakistan: a matched case-control study.

BACKGROUND: In April, 2019, an HIV outbreak predominantly affecting children occurred in Larkana District, Sindh, Pakistan. By December, 2019, 881 (4·0%) of 21 962 children screened for HIV had tested positive. We aimed to assess factors associated with HIV infection in this outbreak. METHODS: In this individually matched case-control study, we sampled 406 cases (individuals aged <16 years who had registered for paediatric HIV care at the HIV Treatment Centre at Shaikh Zayed Children's Hospital in Larkana City, Pakistan) and 406 controls (individuals without HIV matched by age, sex, and neighbourhood residence, recruited through doorknocking at houses adjacent to case participants). An interviewer-administered questionnaire was used to collect data on possible risk factors for HIV acquisition and a blood sample was collected from all participants for hepatitis B and hepatitis C serology. Mothers of all participants underwent HIV testing. Odds ratios were estimated using conditional logistic regression to assess factors associated with HIV infection. FINDINGS: 406 case-control pairs were recruited between July 3 and Dec 26, 2019. Five pairs were excluded (three pairs had an age mismatch and two pairs were duplicate cases) and 401 were analysed. The prevalence of hepatitis B surface antigen was 18·2% (95% CI 14·5-22·3) among cases and 5·2% (3·3-7·9) among controls, and the prevalence of hepatitis C antibodies was 6·5% (95% CI 4·3-9·4) among cases and 1·0% (0·3-2·5) among controls. 28 (7%) of 397 mothers of cases for whom we had data, and no mothers of 394 controls, were HIV positive. In the 6 months before recruitment, 226 (56%) of 401 cases and 32 (8%) of 401 controls reported having more than ten injections, and 291 (73%) cases and 78 (19%) controls had received an intravenous infusion. At least one blood transfusion was reported in 56 (14%) cases and three (1%) controls in the past 2 years. HIV infection was associated with a history of more injections and infusions (adjusted odds ratio 1·63; 95% CI 1·30-2·04, p<0·0001), blood transfusion (336·75; 23·69-4787·01, p<0·0001), surgery (399·75, 13·99-11 419·39, p=0·0005), the child's mother being HIV positive or having died (3·13, 1·20-8·20, p=0·020), and increased frequency of private clinic (p<0·0001) and government hospital visits (p<0·0001), adjusting for confounders. INTERPRETATION: The predominant mode of HIV transmission in this outbreak was parenteral, probably due to unsafe injection practices and poor blood safety practices. General practitioners across Pakistan need training and systems support in reducing injection use, and in providing safe injections and transfusions only when necessary. FUNDING: Department of Pediatrics and Child Health, the Aga Khan University, Karachi, Pakistan.

Types and risk factors of violence experienced by people living with HIV, Pakistan: a cross-sectional study.

BACKGROUND: In Pakistan, data are lacking on the violence experienced by people living with HIV. AIMS: This study determined the prevalence and risk factors of violence (physical, psychological and sexual) in people living with HIV in Karachi, Pakistan. METHODS: This was a cross-sectional study in 2016 of people living with HIV attending clinics of Bridge Consultants Foundation, a community-based care provider. Date were collected using an interview-based questionnaire. Multivariate logistic regression analysis was done to assess the risk factors for violence with adjusted odds ratios (aOR) and 95% confidence intervals (CI) presented. RESULTS: The sample included 250 people living with HIV; 183 were men, 60 were women and 7 were transgender. The mean (standard deviation) age of the participants was 30 (6.5) years. The prevalence rates of psychological, sexual and physical violence were 79.6%, 74.8% and 64.4%, respectively. More women experienced physical violence than men (76.2% versus 60.7%). Psychological violence was associated with injecting drug use (aOR = 2.64, 95% CI: 1.27-5.50) and being married (aOR = 0.46, 95% CI: 0.24-0.90). Marriage (aOR = 2.30, 95% CI: 1.27-4.16) and having an HIV-positive partner (aOR = 2.07, 95% CI: 1.09-3.92) were risk factors for sexual violence. Physical violence was associated with young age (aOR = 0.95; 95% CI: 0.91-0.99) and having an HIV-positive partner (aOR = 2.17, 95% CI: 1.10-4.26). CONCLUSION: Violence is an important public health problem affecting people living with HIV in Pakistan. This issue needs to be addressed by the government and nongovernmental organizations.

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

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

Spectroscopic and photostability study of water-soluble hypericin encapsulated with polyvinylpyrrolidone.

Hypericin has gained great attention as a powerful photosensitizing and fluorescent agent for photodynamic therapy (PDT) and fluorescence diagnosis (FD) of cancer. However, native hypericin is hydrophobic and nearly insoluble in aqueous media which hinders its photobiological activity. Herein, we demonstrate the encapsulation of hypericin and polyvinylpyrrolidone (hypericin@PVP) as an attractive class of water-soluble formula of hypericin with improved absorption and emission characteristics in water. The absorption and fluorescence properties of the water-soluble hypericin@PVP were studied at room temperature. Also, the photostability of the prepared hypericin@PVP was studied under visible light irradiation. The absorbance and emission measurements confirm the association and binding of hypericin and PVP with a binding constant (Kb) of 1.2 × 105 M-1. The interaction between hypericin and PVP in water could lead to the dissociation of aggregated hypericin into their monomeric state which is crucial for effective photobiological implementation in PDT and FD. Upon encapsulation with PVP, hypericin showed a significant increase in the fluorescence properties with an enhanced emission intensity of 300% at a PVP concentration of 1 × 10-4 M. Moreover, water-soluble hypericin@PVP demonstrated high photostability under visible light irradiation with an irradiance of 15 mW/cm2 and exposure time up to 150 min. This enhancement in the absorption, emission, and photostability of hypericin in water is related to the effects of encapsulation with PVP and the unique spectroscopic properties of the formulated hypericin@PVP.

Correlation of Subclinical Hypothyroidism With Polycystic Ovary Syndrome (PCOS).

Aim The correlation of subclinical hypothyroidism (SCH) and polycystic ovary syndrome (PCOS) is a still insufficiently explored entity. The aim of this study was to determine the correlation between SCH and PCOS along with the impact of SCH on metabolic and hormonal parameters in women with PCOS. Methodology This cross-sectional study was conducted at the Gynecology Outpatient Department of Ziauddin Hospital Kemari, Karachi, Pakistan, from June 2019 to December 2019. A total of 90 diagnosed cases of PCOS were enrolled in the study. A non-probability consecutive sampling technique was used. After taking informed consent, participants were evaluated through clinical interviews, a questionnaire, and anthropometric measurements. The participants underwent the following assessments, i.e., transabdominal ultrasonography, hormonal profile (free testosterone, follicle-stimulating hormone, luteinizing hormone), and fasting blood sugar. Participants were divided into two groups based on thyroid-stimulating hormone (TSH) into the euthyroid group and subclinical hypothyroid (SCH) group. The Mann-Whitney test was used for comparing the two groups. Results Our results showed a significant difference in weight, body mass index (BMI), insulin, homeostatic model assessment of insulin resistance (HOMA-IR), and TSH were found in the SCH group as compared to the euthyroid group. A significant correlation of TSH with waist-hip ratio (WHR), weight, body mass index (BMI), insulin, and the homeostatic model assessment of insulin resistance (HOMA-IR) in PCOS patients. Conclusion This study showed a significant correlation of subclinical hypothyroidism with polycystic ovary syndrome. We found subclinical hypothyroidism may aggravate the insulin resistance; therefore, PCOS patients must be screened with a thyroid profile.

Correlates of Knowledge of HIV Transmission Among Incident Cases of HIV in a Cohort of Injection Drug Users Receiving Harm Reduction Services at Karachi, Pakistan.

BACKGROUND: Knowledge of risk factors for HIV transmission in high-risk population plays a critical role in averting the risk of HIV transmission. In Pakistan, injection drug users (IDUs) constitute the core risk group of HIV prevalence, where the epidemic has transitioned to a "concentrated level." Still nothing is known about the role of knowledge in HIV transmission and HIV sero-conversion among IDUs in Pakistan. METHODS: From 2009 to 2011, a nested case-control study was conducted in a cohort of 636 IDUs receiving harm reduction services in the mega city of Karachi. RESULTS: In multivariable regression analysis, 3 factors, namely HIV does not spread through unprotected sex (adjusted odds ratio [AOR]: 3.1, 95% confidence interval [CI] 1.39-6.90, P value .01), HIV does not transmit by sharing syringes (AOR: 3.5, 95% CI 1.97-6.40, P value <.00), and the risk of HIV cannot be minimized by using new syringe every time (AOR: 2.0, 95% CI 1.16-3.60, P value .01), were significantly associated with the incident cases of HIV. CONCLUSION: The study findings suggest the association between knowledge of HIV transmission and HIV sero-incident cases.

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

A number of HIV-1 subtypes are identified in Pakistan by characterization of partial viral gene sequences. Little is known whether new recombinants are generated and how they disseminate since whole genome sequences for these viruses have not been characterized. Near full-length genome (NFLG) sequences were obtained by amplifying two overlapping half genomes or next generation sequencing from 34 HIV-1-infected individuals in Pakistan. Phylogenetic tree analysis showed that the newly characterized sequences were 16 subtype As, one subtype C, and 17 A/G recombinants. Further analysis showed that all 16 subtype A1 sequences (47%), together with the vast majority of sequences from Pakistan from other studies, formed a tight subcluster (A1a) within the subtype A1 clade, suggesting that they were derived from a single introduction. More in-depth analysis of 17 A/G NFLG sequences showed that five shared similar recombination breakpoints as in CRF02 (15%) but were phylogenetically distinct from the prototype CRF02 by forming a tight subcluster (CRF02a) while 12 (38%) were new recombinants between CRF02a and A1a or a divergent A1b viruses. Unique recombination patterns among the majority of the newly characterized recombinants indicated ongoing recombination. Interestingly, recombination breakpoints in these CRF02/A1 recombinants were similar to those in prototype CRF02 viruses, indicating that recombination at these sites more likely generate variable recombinant viruses. The dominance and fast dissemination of new CRF02a/A1 recombinants over prototype CRF02 suggest that these recombinant have more adapted and may become major epidemic strains in Pakistan.