Evolutionary analyses of viral sequences can provide insights into transmission dynamics, which in turn can optimize prevention interventions. Here, we characterized the dynamics of HIV transmission within the Mexico City metropolitan area. HIV pol sequences from persons recently diagnosed at the largest HIV clinic in Mexico City (between 2016 and 2021) were annotated with demographic/geographic metadata. A multistep phylogenetic approach was applied to identify putative transmission clades. A data set of publicly available sequences was used to assess international introductions. Clades were analyzed with a discrete phylogeographic model to evaluate the timing and intensity of HIV introductions and transmission dynamics among municipalities in the region. A total of 6,802 sequences across 96 municipalities (5,192 from Mexico City and 1,610 from the neighboring State of Mexico) were included (93.6% cisgender men, 5.0% cisgender women, and 1.3% transgender women); 3,971 of these sequences formed 1,206 clusters, involving 78 municipalities, including 89 clusters of ≥10 sequences. Discrete phylogeographic analysis revealed (i) 1,032 viral introductions into the region, over one-half of which were from the United States, and (ii) 354 migration events between municipalities with high support (adjusted Bayes factor of ≥3). The most frequent viral migrations occurred between northern municipalities within Mexico City, i.e., Cuauhtémoc to Iztapalapa (5.2% of events), Iztapalapa to Gustavo A. Madero (5.4%), and Gustavo A. Madero to Cuauhtémoc (6.5%). Our analysis illustrates the complexity of HIV transmission within the Mexico City metropolitan area but also identifies a spatially active transmission area involving a few municipalities in the north of the city, where targeted interventions could have a more pronounced effect on the entire regional epidemic. IMPORTANCE Phylogeographic investigation of the Mexico City HIV epidemic illustrates the complexity of HIV transmission in the region. An active transmission area involving a few municipalities in the north of the city, with transmission links throughout the region, is identified and could be a location where targeted interventions could have a more pronounced effect on the entire regional epidemic, compared with those dispersed in other manners.

Background

Evolutionary analyses of well-annotated human immunodeficiency virus (HIV) sequence data can provide insights into viral transmission patterns and associated factors. Here, we explored the transmission dynamics of the HIV-1 subtype B epidemic across the San Diego (US) and Tijuana (Mexico) border region to identify factors that could help guide public health policy.

Methods

HIV pol sequences were collected from people with HIV in San Diego County and Tijuana between 1996-2018. A multistep phylogenetic approach was used to characterize the dynamics of spread. The contributions of geospatial factors and HIV risk group to the local dynamics were evaluated.

Results

Phylogeographic analyses of the 2034 sequences revealed an important contribution of local transmission in sustaining the epidemic, as well as a complex viral migration network across the region. Geospatial viral dispersal between San Diego communities occurred predominantly among men who have sex with men, with central San Diego being the main source (34.9%) and recipient (39.5%) of migration events. HIV migration was more frequent from San Diego county towards Tijuana than vice versa. Migrations were best explained by the driving time between locations.

Conclusions

The US-Mexico border may not be a major barrier to the spread of HIV, which may stimulate coordinated transnational intervention approaches. Whereas a focus on central San Diego has the potential to avert most spread, the substantial viral migration independent of central San Diego shows that county-wide efforts will be more effective. Combined, this work shows that epidemiological information gleaned from pathogen genomes can uncover mechanisms that underlie sustained spread and, in turn, can be a building block of public health decision-making.

The HIV-1 accessory protein Vpu exhibits high inter- and intra- subtype genetic diversity that may influence Vpu function and possibly contribute to HIV-1 pathogenesis. However, scalable methods to evaluate genotype/phenotype relationships in natural Vpu sequences are limited, particularly those expressing the protein in CD4+ T-cells, the natural target of HIV-1 infection. A major impediment to assay scalability is the extensive genetic diversity within, and immediately upstream of, Vpu's initial 5' coding region, which has necessitated the design of oligonucleotide primers specific for each individual HIV-1 isolate (or subtype). To address this, we developed two universal forward primers, located in relatively conserved regions 38 and 90 bases upstream of Vpu, and a single universal reverse primer downstream of Vpu, which are predicted to cover the vast majority of global HIV-1 group M sequence diversity. We show that inclusion of up to 90 upstream bases of HIV-1 genomic sequence does not significantly influence in vitro Vpu expression or function when a Rev/Rev Response Element (RRE)-dependent expression system is used. We further assess the function of four diverse HIV-1 Vpu sequences, revealing reproducible and significant differences between them. Our approach represents a scalable option to measure the in vitro function of genetically diverse natural Vpu isolates in a CD4+ T-cell line.

HIV-specific CD8⁺ T cells demonstrate an exhausted phenotype associated with increased expression of inhibitory receptors, decreased functional capacity, and a skewed transcriptional profile, which are only partially restored by antiretroviral treatment (ART). Expression levels of the inhibitory receptor, T cell immunoglobulin and ITIM domain (TIGIT), the co-stimulatory receptor CD226 and their ligand PVR are altered in viral infections and cancer. However, the extent to which the TIGIT/CD226/PVR-axis is affected by HIV-infection has not been characterized. Here, we report that TIGIT expression increased over time despite early initiation of ART. HIV-specific CD8⁺ T cells were almost exclusively TIGIT⁺, had an inverse expression of the transcription factors T-bet and Eomes and co-expressed PD-1, CD160 and 2B4. HIV-specific TIGIT^hi cells were negatively correlated with polyfunctionality and displayed a diminished expression of CD226. Furthermore, expression of PVR was increased on CD4⁺ T cells, especially T follicular helper (Tfh) cells, in HIV-infected lymph nodes. These results depict a skewing of the TIGIT/CD226 axis from CD226 co-stimulation towards TIGIT-mediated inhibition of CD8⁺ T cells, despite early ART. These findings highlight the importance of the TIGIT/CD226/PVR axis as an immune checkpoint barrier that could hinder future "cure" strategies requiring potent HIV-specific CD8⁺ T cells.

The role of lymphoid tissue as a potential source of HIV-1 rebound following interruption of antiretroviral therapy (ART) is uncertain. To address this issue, we compared the latent viruses obtained from CD4+ T cells in peripheral blood and lymph nodes to viruses emerging during treatment interruption. Latent viruses were characterized by sequencing near-full-length (NFL) proviral DNA and env from viral outgrowth assays (VOAs). Five HIV-1-infected individuals on ART were studied, four of whom participated in a clinical trial of a TLR9 agonist that included an analytical treatment interruption. We found that 98% of intact or replication-competent clonal sequences overlapped between blood and lymph node. In contrast, there was no overlap between 205 latent reservoir and 125 rebound sequences in the four individuals who underwent treatment interruption. However, rebound viruses could be accounted for by recombination. The data suggest that CD4+ T cells carrying latent viruses circulate between blood and lymphoid tissues in individuals on ART and support the idea that recombination may play a role in the emergence of rebound viremia.IMPORTANCE HIV-1 persists as a latent infection in CD4+ T cells that can be found in lymphoid tissues in infected individuals during ART. However, the importance of this tissue reservoir and its contribution to viral rebound upon ART interruption are not clear. In this study, we sought to compare latent HIV-1 from blood and lymph node CD4+ T cells from five HIV-1-infected individuals. Further, we analyzed the contribution of lymph node viruses to viral rebound. We observed that the frequencies of intact proviruses were the same in blood and lymph node. Moreover, expanded clones of T cells bearing identical proviruses were found in blood and lymph node. These latent reservoir sequences did not appear to be the direct origin of rebound virus. Instead, latent proviruses were found to contribute to the rebound compartment by recombination.

Background

Pretreatment HIV drug resistance (HIVDR) to NNRTIs has consistently increased in Mexico City during the last decade.

Objectives

To infer the HIV genetic transmission network in Mexico City to describe the dynamics of the local HIV epidemic and spread of HIVDR.

Patients and methods

HIV pol sequences were obtained by next-generation sequencing from 2447 individuals before initiation of ART at the largest HIV clinic in Mexico City (April 2016 to June 2018). Pretreatment HIVDR was estimated using the Stanford algorithm at a Sanger-like threshold (≥20%). Genetic networks were inferred with HIV-TRACE, establishing putative transmission links with genetic distances <1.5%. We examined demographic associations among linked individuals with shared drug resistance mutations (DRMs) using a ≥ 2% threshold to include low-frequency variants.

Results

Pretreatment HIVDR reached 14.8% (95% CI 13.4%-16.2%) in the cohort overall and 9.6% (8.5%-10.8%) to NNRTIs. Putative links with at least one other sequence were found for 963/2447 (39%) sequences, forming 326 clusters (2-20 individuals). The inferred network was assortative by age and municipality (P < 0.001). Clustering individuals were younger [adjusted OR (aOR) per year = 0.96, 95% CI 0.95-0.97, P < 0.001] and less likely to include women (aOR = 0.46, 95% CI 0.28-0.75, P = 0.002). Among clustering individuals, 175/963 (18%) shared DRMs (involving 66 clusters), of which 66/175 (38%) shared K103N/S (24 clusters). Eight municipalities (out of 75) harboured 65% of persons sharing DRMs. Among all persons sharing DRMs, those sharing K103N were younger (aOR = 0.93, 95% CI 0.88-0.98, P = 0.003).

Conclusions

Our analyses suggest age- and geographically associated transmission of DRMs within the HIV genetic network in Mexico City, warranting continuous monitoring and focused interventions.

On September 10, 2021, a special tribunal established by the French government launched an inquiry into the activities of former health minister Dr. Agnes Buzyn who was charged with "endangering the lives of others". It is surprising to learn of this accusation and inquiry into the actions of a public health official whose response to the epidemic was, to all appearances, exemplary.

Background

Migration and travel are major drivers of the spread of infectious diseases. Geographic proximity and a common language facilitate travel and migration in Mesoamerica, which in turn could affect the spread of HIV in the region.

Methods

6092 HIV-1 subtype B partial pol sequences sampled from unique antiretroviral treatment-naïve individuals from Mexico (40.7%), Guatemala (24.4%), Honduras (19%), Panama (8.2%), Nicaragua (5.5%), Belize (1.4%), and El Salvador (0.7%) between 2011 and 2016 were included. Phylogenetic and genetic network analyses were performed to infer putative relationships between HIV sequences. The demographic and geographic associations with clustering were analyzed and viral migration patterns were inferred using the Slatkin-Maddison approach on 100 iterations of random subsets of equal number of sequences per location.

Results

A total of 1685/6088 (27.7%) of sequences linked with at least one other sequence, forming 603 putative transmission clusters (range: 2-89 individuals). Clustering individuals were significantly more likely to be younger (median age 29 vs 33years, p<0.01) and men-who-have-sex-with-men (40.4% vs 30.3%, p<0.01). Of the 603 clusters, 30 (5%) included sequences from multiple countries with commonly observed linkages between Mexican and Honduran sequences. Eight of the 603 clusters included >10 individuals, including two comprised exclusively of Guatemalans (52 and 89 individuals). Phylogenetic and migration analyses suggested that the Central and Southern regions of Mexico along with Belize were major sources of HIV throughout the region (p<0.01) with genetic flow southward from Mexico to the other nations of Mesoamerica. We also found evidence of significant viral migration within Mexico.

Conclusion

International clusters were infrequent, suggesting moderate migration between HIV epidemics of the different Mesoamerican countries. Nevertheless, we observed important sources of transnational HIV spread in the region, including Southern and Central Mexico and Belize.

Background

Pretreatment HIV drug resistance (HIVDR) to NNRTIs has consistently increased in Mexico City during the last decade.

Objectives

To infer the HIV genetic transmission network in Mexico City to describe the dynamics of the local HIV epidemic and spread of HIVDR.

Patients and methods

HIV pol sequences were obtained by next-generation sequencing from 2447 individuals before initiation of ART at the largest HIV clinic in Mexico City (April 2016 to June 2018). Pretreatment HIVDR was estimated using the Stanford algorithm at a Sanger-like threshold (≥20%). Genetic networks were inferred with HIV-TRACE, establishing putative transmission links with genetic distances <1.5%. We examined demographic associations among linked individuals with shared drug resistance mutations (DRMs) using a ≥ 2% threshold to include low-frequency variants.

Results

Pretreatment HIVDR reached 14.8% (95% CI 13.4%-16.2%) in the cohort overall and 9.6% (8.5%-10.8%) to NNRTIs. Putative links with at least one other sequence were found for 963/2447 (39%) sequences, forming 326 clusters (2-20 individuals). The inferred network was assortative by age and municipality (P < 0.001). Clustering individuals were younger [adjusted OR (aOR) per year = 0.96, 95% CI 0.95-0.97, P < 0.001] and less likely to include women (aOR = 0.46, 95% CI 0.28-0.75, P = 0.002). Among clustering individuals, 175/963 (18%) shared DRMs (involving 66 clusters), of which 66/175 (38%) shared K103N/S (24 clusters). Eight municipalities (out of 75) harboured 65% of persons sharing DRMs. Among all persons sharing DRMs, those sharing K103N were younger (aOR = 0.93, 95% CI 0.88-0.98, P = 0.003).

Conclusions

Our analyses suggest age- and geographically associated transmission of DRMs within the HIV genetic network in Mexico City, warranting continuous monitoring and focused interventions.

Introduction

Molecular surveillance systems could provide public health benefits to focus strategies to improve the HIV care continuum. Here, we infer the HIV genetic network of Mexico City in 2020, and identify actively growing clusters that could represent relevant targets for intervention.

Methods

All new diagnoses, referrals from other institutions, as well as persons returning to care, enrolling at the largest HIV clinic in Mexico City were invited to participate in the study. The network was inferred from HIV pol sequences, using pairwise genetic distance methods, with a locally hosted, secure version of the HIV-TRACE tool: Seguro HIV-TRACE. Socio-demographic, clinical and behavioural metadata were overlaid across the network to design focused prevention interventions.

Results

A total of 3168 HIV sequences from unique individuals were included. One thousand and one-hundred and fifty (36%) sequences formed 1361 links within 386 transmission clusters in the network. Cluster size varied from 2 to 14 (63% were dyads). After adjustment for covariates, lower age (adjusted odds ratio [aOR]: 0.37, p<0.001; >34 vs. <24 years), being a man who has sex with men (MSM) (aOR: 2.47, p = 0.004; MSM vs. cisgender women), having higher viral load (aOR: 1.28, p<0.001) and higher CD4+ T cell count (aOR: 1.80, p<0.001; ≥500 vs. <200 cells/mm³ ) remained associated with higher odds of clustering. Compared to MSM, cisgender women and heterosexual men had significantly lower education (none or any elementary: 59.1% and 54.2% vs. 16.6%, p<0.001) and socio-economic status (low income: 36.4% and 29.0% vs. 18.6%, p = 0.03) than MSM. We identified 10 (2.6%) clusters with constant growth, for prioritized intervention, that included intersecting sexual risk groups, highly connected nodes and bridge nodes between possible sub-clusters with high growth potential.

Conclusions

HIV transmission in Mexico City is strongly driven by young MSM with higher education level and recent infection. Nevertheless, leveraging network inference, we identified actively growing clusters that could be prioritized for focused intervention with demographic and risk characteristics that do not necessarily reflect the ones observed in the overall clustering population. Further studies evaluating different models to predict growing clusters are warranted. Focused interventions will have to consider structural and risk disparities between the MSM and the heterosexual populations.