Molecular epidemiology of HIV-1 infection in the Philippines, 1985 to 1997: transmission of subtypes B and E and potential emergence of subtypes C and F

Reconstructing the phylodynamic history and geographic spread of the CRF01_AE-predominant HIV-1 epidemic in the Philippines from PR/RT sequences sampled from 2008 to 2018

The Philippines has had a rapidly growing human immunodeficiency virus (HIV) epidemic with a shift in the prevalent subtype from B to CRF01_AE. However, the phylodynamic history of CRF01_AE in the Philippines has yet to be reconstructed. We conducted a descriptive retrospective study reconstructing the history of HIV-1 CRF01_AE transmissions in the Philippines through molecular epidemiology. Partial polymerase sequences (n = 1144) collected between 2008 and 2018 from three island groups were collated from the Research Institute for Tropical Medicine drug resistance genotyping database. Estimation of the time to the most recent common ancestor (tMRCA), effective reproductive number (Re), effective viral population size (Ne), relative migration rates, and geographic spread of CRF01_AE was performed with BEAST. Re and Ne were compared between CRF01_AE and B. Most CRF01_AE sequences formed a single clade with a tMRCA of June 1996 [95 per cent highest posterior density (HPD): December 1991, October 1999]. An increasing CRF01_AE Ne was observed from the tMRCA to 2013. The CRF01_AE Re reached peaks of 2.46 [95 per cent HPD: 1.76, 3.27] in 2007 and 2.52 [95 per cent HPD: 1.83, 3.34] in 2015. A decrease of CRF01_AE Re occurred in the intervening years of 2007 to 2011, reaching as low as 1.43 [95 per cent HPD: 1.06, 1.90] in 2011, followed by a rebound. The CRF01_AE epidemic most likely started in Luzon and then spread to the other island groups of the country. Both CRF01_AE and Subtype B exhibited similar patterns of Re fluctuation over time. These results characterize the subtype-specific phylodynamic history of the largest CRF01_AE cluster in the Philippines, which contextualizes and may inform past, present, and future public health measures toward controlling the HIV epidemic in the Philippines.

The State of the HIV Epidemic in the Philippines: Progress and Challenges in 2023

In the past decade, the Philippines has gained notoriety as the country with the fastest-growing human immunodeficiency virus (HIV) epidemic in the Western Pacific region. While the overall trends of HIV incidence and acquired immunodeficiency syndrome (AIDS)-related deaths are declining globally, an increase in new cases was reported to the HIV/AIDS and ART Registry of the Philippines. From 2012 to 2023, there was a 411% increase in daily incidence. Late presentation in care remains a concern, with 29% of new confirmed HIV cases in January 2023 having clinical manifestations of advanced HIV disease at the time of diagnosis. Men having sex with men (MSM) are disproportionately affected. Various steps have been taken to address the HIV epidemic in the country. The Philippine HIV and AIDS Policy Act of 2018 (Republic Act 11166) expanded access to HIV testing and treatment. HIV testing now allows for the screening of minors 15-17 years old without parental consent. Community-based organizations have been instrumental in expanding HIV screening to include self-testing and community-based screening. The Philippines moved from centralized HIV diagnosis confirmation by Western blot to a decentralized rapid HIV diagnostic algorithm (rHIVda). Dolutegravir-based antiretroviral therapy is now the first line. Pre-exposure prophylaxis in the form of emtricitabine-tenofovir disoproxil fumarate has been rolled out. The number of treatment hubs and primary HIV care facilities continues to increase. Despite these efforts, barriers to ending the HIV epidemic remain, including continued stigma, limited harm reduction services for people who inject drugs, sociocultural factors, and political deterrents. HIV RNA quantification and drug resistance testing are not routinely performed due to associated costs. The high burden of tuberculosis and hepatitis B virus co-infection complicate HIV management. CRF_01AE is now the predominant subtype, which has been associated with poorer clinical outcomes and faster CD4 T-cell decline. The HIV epidemic in the Philippines requires a multisectoral approach and calls for sustained political commitment, community involvement, and continued collaboration among various stakeholders. In this article, we outline the current progress and challenges in curbing the HIV epidemic in the Philippines.

Prevalent subtypes and one-year outcomes of an HIV-cohort from an urban Philippine center

ABSTRACT

Circulating HIV subtypes in the Philippines have increasingly diversified, potentially affecting treatment. We monitored outcomes of a treatment-naïve cohort and their virus subtype prevalence.Retrospective/prospective study cohort.HIV-I-REACT clinic patients co-enrolled in the Virology Quality Assurance Program (RUSH-VQA) from 7/2017-6/2019 were included. Relevant demographic and laboratory information were collected. The ViroSeq HIV-1 Genotyping System v.3 and HIV-1 Integrase Genotyping Kit identified protease-reverse transcriptase and integrase drug resistance mutations (DRM). Sequence subtyping followed using the Stanford University Drug Resistance Database and the REGA HIV-1 Subtyping Tool v.3. The jpHMM HIV-1 Tool and REGA HIV-1 Subtyping Tool provided additional subtype analysis of this cohort's 5'LTR-VIF regions after Sanger sequencing. One-year outcomes included virologic suppression, mortality, and follow-up.86/88 patients were males. Median age was 30 (range 19-65) years; 61/88 were MSM. 15/85 carried baseline DRM. ViroSeq-generated sequences included subtypes CRF01_AE (66/85), B (14/85), and newer recombinants (4/85). Extensive sequencing (n = 71) of the 5'-LTR-GAG-Pol genes showed CRF01_AE (n = 50), subtype B (n = 7), and other recombinants (n = 13). Bootstrap analysis identified 7 pairs of highly related strains. Discordant DRM appeared in 2/7 pairs, where 1/2 strains displayed DRM. After 1 year, 87 individuals were alive, with 19 lost to care. Viral load (VL) was repeated for only 31/77 (40.2%). Follow-up CD4 testing for 39/77 (50.6%) showed an increase to a median of 327 cells/mm3.Our cohort currently carries subtype CRF01_AE (∼68%-70%), followed by subtype B and CRF01_AE/B recombinants. Outcomes were favorable, regardless of subtype after 1 year on cART.

The changing molecular epidemiology of HIV in the Philippines

BACKGROUND

The Philippines has one of the fastest-growing HIV epidemics in the world. Possible reasons for this include increased testing, increased local transmission, and possibly more aggressive strains of HIV. This study sought to determine whether local molecular subtypes of HIV have changed.

METHODS

Viruses from 81 newly diagnosed, treatment-naive HIV patients were genotyped using protease and reverse transcriptase genes. Demographic characteristics and CD4 count data were collected.

RESULTS

The cohort had an average age of 29 years (range 19-51 years), CD4+ count of 255 cells/mm³ (range 2-744 cells/mm³), and self-reported acquisition time of 2.42 years (range 0.17-8.17 years). All were male, including 79 men who have sex with men (MSM). The genotype distribution was 77% CRF01_AE, 22% B, and 1% C. Previous data from 1985-2000 showed that most Philippine HIV infections were caused by subtype B (71%, n=100), followed by subtype CRF01_AE (20%). Comparison with the present cohort showed a significant shift in subtype (p<0.0001). Comparison between CRF01_AE and B showed a lower CD4+ count (230 vs. 350 cells/mm³, p=0.03). Survival data showed highly significant survival associated with antiretroviral (ARV) treatment (p<0.0001), but no significant difference in mortality or CD4 count increase on ARVs between subtypes.

CONCLUSIONS

The molecular epidemiology of HIV in the Philippines has changed, with the more aggressive CRF01_AE now being the predominant subtype.

Inferring viral population structures using heteroduplex mobility and DNA sequence analyses

Heteroduplex mobility (HMA) and tracking assays (HTA) are used to assess genetic relationships between DNA molecules. While distinguishing relationships between clonal or nearly clonal molecules is relatively straightforward, inferring population structures is more complex. To address this issue, HIV-1 quasispecies with varying levels of diversity were studied using both HTA and DNA sequencing. Viral diversity estimates and the temporal features of virus evolution were found to be generally concordant between HTA and DNA sequencing. In addition, the distribution of pairwise differences and the rates of virus divergence were similar between the two methods. These findings support the use of HTA to characterize variant populations of DNA and strengthen previous inferences concerning the evolution of HIV-1 over the course of infection.

Virus fitness: concept, quantification, and application to HIV population dynamics

Viral fitness has been broadly studied during the past three decades, mainly to test evolutionary models and population theories difficult to analyze and interpret with more complex organisms. More recent studies, however, are focused in the role of fitness on viral transmission, pathogenesis, and drug resistance. Here, we used human immunodeficiency virus (HIV) as one of the most relevant models to evaluate the importance of viral quasispecies and fitness in HIV evolution, population dynamics, disease progression, and potential clinical implications.

Rapid spread of hepatitis C virus among injecting-drug users in the Philippines: Implications for HIV epidemics

From the trends of human immunodeficiency virus (HIV) epidemics in South and Southeast Asia, it was postulated that an HIV epidemic would start as a blood-borne infection among injecting-drug users in the Philippines. In 2002, 560 individuals were recruited in Metro Cebu, Philippines and tested for HIV, hepatitis C virus (HCV), and hepatitis B virus (HBV) infections. The seroprevalence of anti-HCV among injecting-drug users (70.1%, 61/87) was significantly higher than those among inhalation drug users (16.3%, 7/43; P = 0.00; OR = 12), sex workers (0%, 0/130; P = 0.00; OR = infinity), antenatal clinic attendees (0%, 0/100; P = 0.00; OR = infinity), and students/health care workers (2%, 4/200; P = 0.00; OR = 115). The seroprevalence of HBsAg among injecting-drug users (10.3%, 9/87) was significantly higher than those among sex workers (2.3%, 3/130; P = 0.01; OR = 4.9), and antenatal clinic attendees (3%, 3/100; P = 0.04; OR = 3.7), but was not statistically different from those among inhalation drug users (9.3%, 4/43; P = 0.9) and students/health care workers (4.5%, 9/200; P = 0.06). None of the study population was reactive to anti-HIV antibody. The HCV strains obtained from the injecting-drug users belonged to either genotype 1a or 2b and the strains in each genotype clustered closely to each other. There was no dual infection with genotype 1a and 2b. These results suggest that the HCV infection in injecting-drug users may be emanating rapidly from limited number individuals in Metro Cebu, Philippines.

HIV/AIDS in Asia

HIV (ie, HIV-1) epidemics in Asia show great diversity, both in severity and timing. But epidemics in Asia are far from over and several countries including China, Indonesia, and Vietnam have growing epidemics. Several factors affect the rate and magnitude of growth of HIV prevalence, but two of the most important are the size of the sex worker population and the frequency with which commercial sex occurs. In view of the present state of knowledge, even countries with low prevalence of infection might still have epidemics affecting a small percentage of the population. Once HIV infection has become established, growing needs for care and treatment are unavoidable and even the so-called prevention-successful countries of Thailand and Cambodia are seeing burgeoning care needs. The manifestations of HIV disease in the region are discussed with the aim of identifying key issues in medical management and care of HIV/AIDS. In particular, issues relevant to developing appropriate highly active antiretroviral treatment programmes in the region are discussed. Although access to antiretroviral therapy is increasing globally, making it work effectively while simultaneously expanding prevention programmes to stem the flow of new infections remains a real challenge in Asia. Genuine political interest and commitment are essential foundations for success, demanding advocacy at all levels to drive policy, mobilise sufficient resources, and take effective action.

Comparing the ex vivo fitness of CCR5-tropic human immunodeficiency virus type 1 isolates of subtypes B and C

Continual human immunodeficiency virus type 1 (HIV-1) evolution and expansion within the human population have led to unequal distribution of HIV-1 group M subtypes. In particular, recent outgrowth of subtype C in southern Africa, India, and China has fueled speculation that subtype C isolates may be more fit in vivo. In this study, nine subtype B and six subtype C HIV-1 isolates were added to peripheral blood mononuclear cell cultures for a complete pairwise competition experiment. All subtype C HIV-1 isolates were less fit than subtype B isolates (P < 0.0001), but intrasubtype variations in HIV-1 fitness were not significant. Increased fitness of subtype B over subtype C was also observed in primary CD4(+) T cells and macrophages from different human donors but not in skin-derived human Langerhans cells. Detailed analysis of the retroviral life cycle during several B and C virus competitions indicated that the efficiency of host cell entry may have a significant impact on relative fitness. Furthermore, phyletic analyses of fitness differences suggested that, for a recombined subtype B/C HIV-1 isolate, higher fitness mapped to the subtype B env gene rather than the subtype C gag and pol genes. These results suggest that subtype B and C HIV-1 may be transmitted with equal efficiency (Langerhans cell data) but that subtype C isolates may be less fit following initial infection (T-cell and macrophage data) and may lead to slower disease progression.

Introduction of human immunodeficiency virus type 2 infection in the Philippines

The aim of this study was to describe and document the first case of human immunodeficiency virus type 2 (HIV-2) in the Philippines by using serological and molecular techniques and to compare the diversity of this strain to that of strains from other countries. With the introduction of HIV-2 into the country and the presence of diversified strains of HIV-1, the use of highly sensitive assays to detect all these strains is recommended.

Rapid detection of human immunodeficiency virus type 1 subtype e infection by PCR

The CRF01_AE (subtype E) strain of human immunodeficiency virus type 1 (HIV-1), originally reported in Thailand, spread rapidly to and showed prevalence in several countries in Southeast Asia, including Taiwan. This strain was also found in other regions of the world. Based on sequence analysis of the vpu gene, a nested PCR assay including an outer primer pair and a subtype E-specific inner primer pair was developed in this study for rapid detection of subtype E viruses. It was tested with 397 HIV-1-positive samples of known subtypes. For these samples, the sensitivity of detection of subtype E viruses was 100% (127 of 127), and the specificity was 97.8% (264 of 270). Although six samples of either subtype A or G showed a positive PCR, most of the cross-reactivity could be reduced by raising the annealing temperature from 54 degrees C to 63 degrees C. When tested with 195 HIV-positive samples of unknown subtypes, the assay had a sensitivity of 98.0% and a specificity of 98.6%. This is a simple, convenient, and sensitive method for rapid detection of subtype E viruses, especially in regions in which viruses of subtypes B and E are predominant.

Importation of multiple HIV type 1 strains into West Papua, Indonesia (Irian Jaya)

HIV-1 from 16 sexually transmitted disease clinic patients in Timika, West Papua, Indonesia was amplified by RT-PCR and subtyped by a combination of envelope and gag region heteroduplex mobility analysis (HMA) and direct PCR DNA sequencing. HMA showed the presence of 14 subtype E (CRF01_AE) and 2 subtype B HIV-1. Phylogenetic analysis of a 540-bp V3-V4 region of gp120 showed that 9 of 10 CRF01_AE variants clustered tightly with a median distance of 1.3% (range, 0.5 to 2.2%) whereas 1 CRF01_AE variant diverged significantly from the others (median distance, 10.7%; range, 10.1 to 11.8%). One subtype B virus envelope was typical of United States/European strains whereas the other appeared to be related to Thai subtype B' variants. These results reflect the independent introduction of multiple HIV-1 strains into West Papua, with the rapid spread in the majority of infected patients tested of a single strain of HIV-1E (CRF01_AE).

Independent introduction of transmissible F/D recombinant HIV-1 from Africa into Belgium and The Netherlands

Most HIV-1 subtype F viruses described so far have been isolated from individuals originating in South America, Romania, or Central Africa. Previous studies have shown that subtype F viruses from these three areas can be distinguished by phylogenetic tree analysis of various parts of the HIV genome. Subtype F strains circulating in Central Africa and classified as subgroup F2 and F3 have relatively large nucleotide distances from strains of subgroup F1, which includes some African strains, along with strains from Romania and South America. Subtype F strains have now appeared in Europe. In this study, we analyzed the complete gag gene and a large fragment of the pol gene of seven strains of African origin that represent the three F subgroups. At least five of the seven strains appear to be intersubtype recombinants. Of four strains circulating in Belgium and the Netherlands, three were F/D mosaics and the fourth harboured a G(gag)/GH(pol)/F3(env) recombinant structure. Two of the three F/D mosaics showed identical breakpoints and were independently introduced in Belgium and the Netherlands. At least two of the mosaics were further transmitted. The remaining three strains of the seven we studied were isolated from individuals in Cameroon. Two included large or smaller F1 fragments in gag and pol. The third strain was subtype D along the entire gag and pol fragment. A parental African subtype F that showed no evidence for recombination was not found.

HIV-1 subtyping in Salvador, Bahia, Brazil: a city with African sociodemographic characteristics

To investigate the prevalence of the HIV-1 subtypes in different populations from Salvador, Bahia, Brazil, blood samples from 72 HIV-1-seropositive injecting drug users (IDUs) and 62 individuals infected sexually were analyzed using the heteroduplex mobility assay (HMA). In the IDU group, 89.5% were classified as subtype B, 3% as subtype F, and 7.5% showed a B/F HMA profile. In the sexual transmission (ST) group, 95% were identified as B subtype, 3.4% showed a B/F profile, and 1.6% a B/C/E HMA profile. All Brazilian samples that showed multiple reactivities in the HMA analysis clustered on sequencing with B North American/ European HIV-1 isolates in the phylogenetic analysis, whereas the F subtypes clustered with F Brazilian HIV-I isolates. Serologic reactivities of IDU's sera were examined using a panel of synthetic V3 loop peptides representative of the different HIV-1 subtypes. No difference in serologic reactivity between F and B subtype plasma could be observed. Predominance of HIV-I subtype B was identified in both study groups, whereas subtype F was detected only among IDUs in a frequency lower than described for other Brazilian regions.