Santiago Avila-Rios

Geographic and Temporal Trends in the Molecular Epidemiology and Genetic Mechanisms of Transmitted HIV-1 Drug Resistance: An Individual-Patient- and Sequence-Level Meta-Analysis

Background Regional and subtype-specific mutational patterns of HIV-1 transmitted drug resistance (TDR) are essential for informing first-line antiretroviral (ARV) therapy guidelines and designing diagnostic assays for use in regions where standard genotypic resistance testing is not affordable. We sought to understand the molecular epidemiology of TDR and to identify the HIV-1 drug-resistance mutations responsible for TDR in different regions and virus subtypes. Methods and Findings We reviewed all GenBank submissions of HIV-1 reverse transcriptase sequences with or without protease and identified 287 studies published between March 1, 2000, and December 31, 2013, with more than 25 recently or chronically infected ARV-naïve individuals. These studies comprised 50,870 individuals from 111 countries. Each set of study sequences was analyzed for phylogenetic clustering and the presence of 93 surveillance drug-resistance mutations (SDRMs). The median overall TDR prevalence in sub-Saharan Africa (SSA), south/southeast Asia (SSEA), upper-income Asian countries, Latin America/Caribbean, Europe, and North America was 2.8%, 2.9%, 5.6%, 7.6%, 9.4%, and 11.5%, respectively. In SSA, there was a yearly 1.09-fold (95% CI: 1.05–1.14) increase in odds of TDR since national ARV scale-up attributable to an increase in non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance. The odds of NNRTI-associated TDR also increased in Latin America/Caribbean (odds ratio [OR] = 1.16; 95% CI: 1.06–1.25), North America (OR = 1.19; 95% CI: 1.12–1.26), Europe (OR = 1.07; 95% CI: 1.01–1.13), and upper-income Asian countries (OR = 1.33; 95% CI: 1.12–1.55). In SSEA, there was no significant change in the odds of TDR since national ARV scale-up (OR = 0.97; 95% CI: 0.92–1.02). An analysis limited to sequences with mixtures at less than 0.5% of their nucleotide positions—a proxy for recent infection—yielded trends comparable to those obtained using the complete dataset. Four NNRTI SDRMs—K101E, K103N, Y181C, and G190A—accounted for >80% of NNRTI-associated TDR in all regions and subtypes. Sixteen nucleoside reverse transcriptase inhibitor (NRTI) SDRMs accounted for >69% of NRTI-associated TDR in all regions and subtypes. In SSA and SSEA, 89% of NNRTI SDRMs were associated with high-level resistance to nevirapine or efavirenz, whereas only 27% of NRTI SDRMs were associated with high-level resistance to zidovudine, lamivudine, tenofovir, or abacavir. Of 763 viruses with TDR in SSA and SSEA, 725 (95%) were genetically dissimilar; 38 (5%) formed 19 sequence pairs. Inherent limitations of this study are that some cohorts may not represent the broader regional population and that studies were heterogeneous with respect to duration of infection prior to sampling. Conclusions Most TDR strains in SSA and SSEA arose independently, suggesting that ARV regimens with a high genetic barrier to resistance combined with improved patient adherence may mitigate TDR increases by reducing the generation of new ARV-resistant strains. A small number of NNRTI-resistance mutations were responsible for most cases of high-level resistance, suggesting that inexpensive point-mutation assays to detect these mutations may be useful for pre-therapy screening in regions with high levels of TDR. In the context of a public health approach to ARV therapy, a reliable point-of-care genotypic resistance test could identify which patients should receive standard first-line therapy and which should receive a protease-inhibitor-containing regimen.

Global epidemiology of drug resistance after failure of WHO recommended first-line regimens for adult HIV-1 infection: A multicentre retrospective cohort study

Summary Background Antiretroviral therapy (ART) is crucial for controlling HIV-1 infection through wide-scale treatment as prevention and pre-exposure prophylaxis (PrEP). Potent tenofovir disoproxil fumarate-containing regimens are increasingly used to treat and prevent HIV, although few data exist for frequency and risk factors of acquired drug resistance in regions hardest hit by the HIV pandemic. We aimed to do a global assessment of drug resistance after virological failure with first-line tenofovir-containing ART. Methods The TenoRes collaboration comprises adult HIV treatment cohorts and clinical trials of HIV drug resistance testing in Europe, Latin and North America, sub-Saharan Africa, and Asia. We extracted and harmonised data for patients undergoing genotypic resistance testing after virological failure with a first-line regimen containing tenofovir plus a cytosine analogue (lamivudine or emtricitabine) plus a non-nucleotide reverse-transcriptase inhibitor (NNRTI; efavirenz or nevirapine). We used an individual participant-level meta-analysis and multiple logistic regression to identify covariates associated with drug resistance. Our primary outcome was tenofovir resistance, defined as presence of K65R/N or K70E/G/Q mutations in the reverse transcriptase (RT) gene. Findings We included 1926 patients from 36 countries with treatment failure between 1998 and 2015. Prevalence of tenofovir resistance was highest in sub-Saharan Africa (370/654 [57%]). Pre-ART CD4 cell count was the covariate most strongly associated with the development of tenofovir resistance (odds ratio [OR] 1·50, 95% CI 1·27–1·77 for CD4 cell count <100 cells per μL). Use of lamivudine versus emtricitabine increased the risk of tenofovir resistance across regions (OR 1·48, 95% CI 1·20–1·82). Of 700 individuals with tenofovir resistance, 578 (83%) had cytosine analogue resistance (M184V/I mutation), 543 (78%) had major NNRTI resistance, and 457 (65%) had both. The mean plasma viral load at virological failure was similar in individuals with and without tenofovir resistance (145 700 copies per mL [SE 12 480] versus 133 900 copies per mL [SE 16 650; p=0·626]). Interpretation We recorded drug resistance in a high proportion of patients after virological failure on a tenofovir-containing first-line regimen across low-income and middle-income regions. Effective surveillance for transmission of drug resistance is crucial. Funding The Wellcome Trust.Summary Background Antiretroviral therapy (ART) is crucial for controlling HIV-1 infection through wide-scale treatment as prevention and pre-exposure prophylaxis (PrEP). Potent tenofovir disoproxil fumarate-containing regimens are increasingly used to treat and prevent HIV, although few data exist for frequency and risk factors of acquired drug resistance in regions hardest hit by the HIV pandemic. We aimed to do a global assessment of drug resistance after virological failure with first-line tenofovir-containing ART. Methods The TenoRes collaboration comprises adult HIV treatment cohorts and clinical trials of HIV drug resistance testing in Europe, Latin and North America, sub-Saharan Africa, and Asia. We extracted and harmonised data for patients undergoing genotypic resistance testing after virological failure with a first-line regimen containing tenofovir plus a cytosine analogue (lamivudine or emtricitabine) plus a non-nucleotide reverse-transcriptase inhibitor (NNRTI; efavirenz or nevirapine). We used an individual participant-level meta-analysis and multiple logistic regression to identify covariates associated with drug resistance. Our primary outcome was tenofovir resistance, defined as presence of K65R/N or K70E/G/Q mutations in the reverse transcriptase ( RT ) gene. Findings We included 1926 patients from 36 countries with treatment failure between 1998 and 2015. Prevalence of tenofovir resistance was highest in sub-Saharan Africa (370/654 [57%]). Pre-ART CD4 cell count was the covariate most strongly associated with the development of tenofovir resistance (odds ratio [OR] 1·50, 95% CI 1·27–1·77 for CD4 cell count Interpretation We recorded drug resistance in a high proportion of patients after virological failure on a tenofovir-containing first-line regimen across low-income and middle-income regions. Effective surveillance for transmission of drug resistance is crucial. Funding The Wellcome Trust.

Highly Pathogenic Avian Influenza A(H7N3) Virus in Poultry Workers, Mexico, 2012

We identified 2 poultry workers with conjunctivitis caused by highly pathogenic avian influenza A(H7N3) viruses in Jalisco, Mexico. Genomic and antigenic analyses of 1 isolate indicated relatedness to poultry and wild bird subtype H7N3 viruses from North America. This isolate had a multibasic cleavage site that might have been derived from recombination with host rRNA.

Severe 2009 pandemic influenza A (H1N1) infection and increased mortality in patients with late and advanced HIV disease.

Objective:To describe the clinical course of infection by 2009 (H1N1) influenza virus in different stages of HIV disease. Design:Prospective, observational study. Methods:During the pandemic period, HIV-infected patients presenting respiratory symptoms at a third level referral hospital in Mexico City were tested for 2009 influenza A (H1N1) viral RNA. Clinical files were prospectively analyzed. Results:Infection by H1N1 was confirmed in 30 (23.8%) of the total 126 HIV-infected patients studied. In the group of patients with 2009 H1N1 virus infection, 16 (53.3%) were hospitalized, 12 (40%) had active opportunistic infections and six (20%) died. In the group of 96 patients not infected with 2009 H1N1 virus, 54 (56.25%) were hospitalized with opportunistic infections and 12 (12.5%) died. For all hospitalized patients, being on HAART and having undetectable HIV viral loads at hospitalization was associated with higher survival (P = 0.019). Patients with 2009 H1N1 virus infection had a higher mortality rate, even after adjusting for HAART (P = 0.043). Coinfection by HIV and H1N1 2009 virus was more severe in patients with opportunistic infections, as shown by longer hospital stays (P = 0.0013), higher rates of hospitalization (P < 0.0001), use of mechanical ventilation (P = 0.0086) and death (P = 0.026). Delayed administration of oseltamivir in hospitalized patients was significantly associated with mortality (P = 0.0022). Conclusion:Our data suggest that infection by 2009 H1N1 is more severe in HIV-infected patients with late and advanced HIV disease than in well controlled patients under HAART.

Unique features of HLA-mediated HIV evolution in a Mexican cohort: a comparative study

BackgroundMounting evidence indicates that HLA-mediated HIV evolution follows highly stereotypic pathways that result in HLA-associated footprints in HIV at the population level. However, it is not known whether characteristic HLA frequency distributions in different populations have resulted in additional unique footprints.MethodsThe phylogenetic dependency network model was applied to assess HLA-mediated evolution in datasets of HIV pol sequences from free plasma viruses and peripheral blood mononuclear cell (PBMC)-integrated proviruses in an immunogenetically unique cohort of Mexican individuals. Our data were compared with data from the IHAC cohort, a large multi-center cohort of individuals from Canada, Australia and the USA.ResultsForty three different HLA-HIV codon associations representing 30 HLA-HIV codon pairs were observed in the Mexican cohort (q < 0.2). Strikingly, 23 (53%) of these associations differed from those observed in the well-powered IHAC cohort, strongly suggesting the existence of unique characteristics in HLA-mediated HIV evolution in the Mexican cohort. Furthermore, 17 of the 23 novel associations involved HLA alleles whose frequencies were not significantly different from those in IHAC, suggesting that their detection was not due to increased statistical power but to differences in patterns of epitope targeting. Interestingly, the consensus differed in four positions between the two cohorts and three of these positions could be explained by HLA-associated selection. Additionally, different HLA-HIV codon associations were seen when comparing HLA-mediated selection in plasma viruses and PBMC archived proviruses at the population level, with a significantly lower number of associations in the proviral dataset.ConclusionOur data support universal HLA-mediated HIV evolution at the population level, resulting in detectable HLA-associated footprints in the circulating virus. However, it also strongly suggests that unique genetic backgrounds in different HIV-infected populations may influence HIV evolution in a particular direction as particular HLA-HIV codon associations are determined by specific HLA frequency distributions. Our analysis also suggests a dynamic HLA-associated evolution in HIV with fewer HLA-HIV codon associations observed in the proviral compartment, which is likely enriched in early archived HIV sequences, compared to the plasma virus compartment. These results highlight the importance of comparative HIV evolutionary studies in immunologically different populations worldwide.

National Prevalence and Trends of HIV Transmitted Drug Resistance in Mexico

Background Transmitted drug resistance (TDR) remains an important concern for the management of HIV infection, especially in countries that have recently scaled-up antiretroviral treatment (ART) access. Methodology/Principal Findings We designed a study to assess HIV diversity and transmitted drug resistance (TDR) prevalence and trends in Mexico. 1655 ART-naïve patients from 12 Mexican states were enrolled from 2005 to 2010. TDR was assessed from plasma HIV pol sequences using Stanford scores and the WHO TDR surveillance mutation list. TDR prevalence fluctuations over back-projected dates of infection were tested. HIV subtype B was highly prevalent in Mexico (99.9%). TDR prevalence (Stanford score>15) in the country for the study period was 7.4% (95% CI, 6.2∶8.8) and 6.8% (95% CI, 5.7∶8.2) based on the WHO TDR surveillance mutation list. NRTI TDR was the highest (4.2%), followed by NNRTI (2.5%) and PI (1.7%) TDR. Increasing trends for NNRTI (p = 0.0456) and PI (p = 0.0061) major TDR mutations were observed at the national level. Clustering of viruses containing minor TDR mutations was observed with some apparent transmission pairs and geographical effects. Conclusions TDR prevalence in Mexico remains at the intermediate level and is slightly lower than that observed in industrialized countries. Whether regional variations in TDR trends are associated with differences in antiretroviral drug usage/ART efficacy or with local features of viral evolution remains to be further addressed.

Molecular Characterization of the Predominant Influenza A(H1N1)pdm09 Virus in Mexico, December 2011–February 2012

When the A(H1N1)pdm09 pandemic influenza virus moved into the post-pandemic period, there was a worldwide predominance of the seasonal influenza A(H3N2) and B viruses. However, A(H1N1)pdm09 became the prevailing subtype in the 2011–2012 influenza season in Mexico and most of Central America. During this season, we collected nasopharyngeal swabs of individuals presenting with influenza-like illness at our institution in Mexico City. Samples were tested for seasonal A(H3N2) and B influenza viruses, as well as A(H1N1)pdm09 by real-time reverse transcription–polymerase chain reaction. Of 205 samples tested, 46% were positive to influenza, all of them A(H1N1)pdm09. The clinical characteristics of patients showed a similar pattern to the 2009 pandemic cases. Using next generation sequencing, we obtained whole genome sequences of viruses from 4 different patients, and in 8 additional viruses we performed partial Sanger sequencing of the HA segment. Non-synonymous changes found in the Mexican isolates with respect to the prototype isolate H1N1 (A/California/04/2009) included HA S69T, K163R and N260D unique to 2012 Mexican and North American isolates and located within or adjacent to HA antigenic sites; HA S143G, S185T, A197T and S203T previously reported in viruses from the 2010–2011 season, located within or adjacent to HA antigenic sites; and HA E374K located in a relevant site for membrane fusion. All Mexican isolates had an oseltamivir-sensitive genotype. Phylogenetic analysis with all 8 influenza gene segments showed that 2012 Mexican sequences formed a robust, distinct cluster. In all cases, 2012 Mexican sequences tended to group with 2010–2011 Asian and European sequences, but not with 2009 Mexican sequences, suggesting a possible recent common ancestor between these latter regions and the 2012 Mexican viruses. It remains to be defined if these viral changes represent an important antigenic drift that would enable viral immune evasion and/or affect influenza vaccine effectiveness.

Pretreatment HIV-drug resistance in Mexico and its impact on the effectiveness of first-line antiretroviral therapy: a nationally representative 2015 WHO survey

BACKGROUND WHO has developed a global HIV-drug resistance surveillance strategy, including assessment of pretreatment HIV-drug resistance. We aimed to do a nationally representative survey of pretreatment HIV-drug resistance in Mexico using WHO-recommended methods. METHODS Among 161 Ministry of Health antiretroviral therapy (ART) clinics in Mexico, the largest, including 90% of ART initiators within the Ministry of Health (66 in total), were eligible for the survey. We used a probability-proportional-to-size design method to sample 25 clinics throughout the country. Consecutive ART-naive patients with HIV about to initiate treatment were invited to participate in the survey; individuals with previous exposure to ART were excluded. We assessed pretreatment HIV-drug resistance by Sanger sequencing and next-generation sequencing of viruses from plasma specimens from eligible participants with Stanford University HIV Drug Resistance Database methods. We obtained follow-up data for a median of 9·4 months (range 6-12) after enrolment. We investigated possible relations between demographic variables and pretreatment drug resistance with univariate and multivariate logistic regression. FINDINGS Between Feb 3 and July 30, 2015, we screened 288 patients in 25 clinics, from whom 264 provided successfully sequenced viruses with no evidence of current exposure to antiretroviral drugs. With the Sanger method, of these 264 participants, 41 (15·5%, 95% CI 11·4-20·5) had pretreatment resistance to any antiretroviral drug and 28 (10·6%, 7·2-15·0) had pretreatment resistance to non-nucleoside reverse transcriptase inhibitors (NNRTIs). At least low-level pretreatment resistance (Stanford penalty score ≥15) was noted in 13 (4 · 9%) of participants to efavirenz and in 23 (8·7%) to the combination tenofovir plus emtricitabine plus efavirenz. With next-generation sequencing, of 264 participants, 38 (14·4%, 95% CI 10·4-19·2) had pretreatment resistance to any antiretroviral drug and 26 (9·8%, 6·5-14·1) had pretreatment resistance to NNRTIs. After median follow-up of 8 months (IQR 6·5-9·4, range 5-11) after ART initiation, 97 (72%) of 135 NNRTI initiators achieved viral suppression (<50 copies per mL) compared with ten (40%) of 25 individuals who started with protease inhibitor-based regimens (p=0·0045). After multivariate regression considering pretreatment resistance and initial ART regimen as composite variables, people starting NNRTIs with pretreatment drug resistance achieved significantly lower viral suppression (odds ratio 0·24, 95% CI 0·07-0·74; p=0·014) than patients without NNRTI resistance. INTERPRETATION High levels of pretreatment drug resistance were noted in Mexico, and NNRTI pretreatment drug resistance significantly reduced the effectiveness of first-line ART regimens based on these drugs. Baseline HIV-drug resistance testing for initial ART follow-up and decision making should be considered. FUNDING The Mexican Government and Consejo Nacional de Ciencia y Tecnología.

Activation of NK cells is associated with HIV-1 disease progression

The main predictor of HIV‐1 disease progression is CD8+ T cell activation, characterized by elevated expression of CD38 and HLA‐DR. NK cells are also activated in viremic HIV‐1‐infected individuals. However, the relationship between NK cell activation and HIV‐1 disease progression remains undefined. We characterized NK cell activation and its association with disease progression in treatment of naive HIV‐1‐infected individuals, who naturally maintained low/undetectable viremia (elite and viremic controllers), compared with progressors and AIDS subjects, and treated individuals. Our results show that CD38 expression on NK cells, predominantly in the cytotoxic CD56dimCD16+ subset, is associated with HIV‐1 disease progression (CD4+ T cell count and pVL), T cell activation (percentage of CD38+HLA‐DR+ T cells), sCD14, inflammation, and innate immune activation. Moreover, NK cell activation is increased in HIV‐1‐infected subjects progressing to AIDS but not in elite and viremic controllers. ART partially reduces the proportion of activated NK cells. Furthermore, our results show that individuals, who naturally control viremia, maintain low levels of innate immune activation similar to those of uninfected controls.

Human Immunodeficiency Virus Type 1 (HIV-1) Subtype B Epidemic in Panama Is Mainly Driven by Dissemination of Country-Specific Clades

The Human immunodeficiency virus type-1 (HIV-1) subtype B is the most predominant clade in Central America; but information about the evolutionary history of this virus in this geographic region is scarce. In this study, we reconstructed the spatiotemporal and population dynamics of the HIV-1 subtype B epidemic in Panama. A total of 761 HIV-1 subtype B pol sequences obtained in Panama between 2004 and 2013 were combined with subtype B pol sequences from the Americas and Europe. Maximum Likelihood phylogenetic analyses revealed that HIV-1 subtype B infections in Panama derived from the dissemination of multiple founder viruses. Most Panamanian subtype B viruses (94.5%) belong to the pandemic viral strain proposed as originated in the US, whereas others (5.5%) were intermixed among non-pandemic Caribbean strains. The bulk (76.6%) of subtype B sequences from Panama grouped within 12 country-specific clades that were not detected in other Central American countries. Bayesian coalescent-based analyses suggest that most Panamanian clades probably originated between the early 1970s and the early 1980s. The root location of major Panamanian clades was traced to the most densely populated districts of Panama province. Major Panamanian clades appear to have experienced one or two periods of exponential growth of variable duration between the 1970s and the 2000s, with median growth rates from 0.2 to 0.4 year− 1. Thus, the HIV-1 subtype B epidemic in Panama is driven by the expansion of local viral strains that were introduced from the Caribbean and other American countries at an early stage of the AIDS pandemic.