Phylogenetic Evidence of Local HIV-1 Transmission and Antiretroviral Drug Resistance in the Middle East and North Africa

The molecular epidemiology and antiretroviral (ARV) drug resistance of human immunodeficiency virus type 1 (HIV-1) remain incompletely outlined in the Middle East and North Africa (MENA). The aim of this retrospective molecular epidemiology study was to analyze MENA HIV-1 sequences for phylogenetic clustering and to delineate surveillance drug-resistance mutations (SDRMs) for nucleoside reverse-transcriptase inhibitors (NRTIs), non-nucleoside reverse-transcriptase inhibitors (NNRTIs), and protease inhibitors (PIs) across various periods, locations, and subtypes/circulating recombinant forms (CRFs). Viral sequences were retrieved from the Los Alamos HIV Sequence Database as of 15 April 2026. Analyses were done using multiple sub-gene regions (two env regions (n = 224 and n = 60) and PR (n = 2413) and RT (n = 2103) of the pol gene). Phylogeny construction was conducted using maximum-likelihood estimation, while ARV drug resistance analysis was conducted using the Stanford HIVdb algorithm. The HIV-1 MENA sequences showed a remarkable genetic diversity, with co-circulation of multiple subtypes/CRFs, including subtype B in the Maghreb, Levant, and Egypt sub-regions, subtypes A1, G, CRF01_AE, and CRF02_AG in the Gulf Cooperation Council (GCC) and Yemen sub-region, and subtypes C and D in the Horn of Africa and Sudan sub-region. The percentage of MENA HIV-1 sequences in clusters was 10.3% for env1, 8.3% for env2, 22.0% for PR and 37.2% for RT. Phylogenetic reconstruction hinted at a structured epidemic dominated by small transmission units, with most clusters comprising dyads (n = 260) or networks (n = 142) and a limited number of large clusters (n = 8) that were largely confined within national boundaries, with only occasional cross-border linkages (n = 8). Overall SDRM prevalence was 3.2% in the PR region and 14.9% in the RT region, with a higher percentage of NNRTI-associated mutations (10.0%) than NRTI-associated mutations (9.1%) and dual-class resistance observed in 4.1% of sequences. Phylogenetic clustering was not associated with the probability of harboring SDRMs; however, negative binomial models showed that non-clustered sequences had a greater burden of NRTI-associated mutations, whereas no such association was observed for NNRTI- or PI-associated mutations. The findings showed predominantly localized and fragmented MENA HIV-1 transmission dynamics. Heterogeneous ARV drug resistance dynamics indicated that resistance emergence might be shaped by broader epidemiologic and treatment-related factors rather than ongoing clustered transmission. There is a need for coordinated molecular surveillance and optimized ART strategies across the MENA countries.