Double-Negative T Cell Number and Phenotype Alterations Before and After Effective Antiretroviral Treatment in Persons Living with HIV.

Double-negative (DN) T cells represent a small and phenotypically heterogeneous population that display regulatory functions. In HIV infection, DN T cells are decreased in peripheral blood and have been negatively associated with T cell activation. This study was aimed at describing the dynamics and phenotypic characteristics of DN T cells in peripheral blood of people living with HIV (PLHIV) before and after antiretroviral therapy (ART) initiation. We included 41 newly diagnosed, ART-naive individuals with advanced HIV infection, who were followed up for 6 months after ART initiation. The control group included 34 people without HIV (PWHIV), on preexposure prophylaxis for HIV infection. DN T cells in peripheral blood were characterized by flow cytometry. The absolute counts of DN T cells were lower in PLHIV than in PWHIV (p = 0.0223), and were particularly low in individuals with advanced HIV disease (p = 0.0311). Activation of DN T cells before ART initiation was directly associated with viral load (VL) (p = 0.0081, r = 0.4083) and inversely associated with CD4+ T cell counts (p = 0.0004, r = -0.4041). Compared with PWHIV, DN T cells of PLHIV expressed higher levels of CD57 (p = 0.0019), Ki67 (p = 0.0065), PD-1 (p = 0.0187), and CD38/HLA-DR (p < 0.0001). After 6 months on ART, expression of Ki67, PD-1, and CD38/HLA-DR on DN T cells returned to similar levels to those observed in PWHIV (p > 0.05 in all cases). However, expression of CD57 decreased only in individuals that start ART with high VL (p = 0.0127). DN T cell counts are decreased in HIV infection. Low DN T cell counts remained despite ART-induced immune reconstitution and viremia control. DN T cell phenotype is altered during chronic untreated infection with a high proportion of proliferating, activated, exhausted, and senescent cells. Most markers return to levels similar to those observed in PWHIV after ART. The impact of altered phenotype of DN T and their regulatory functions warrants further exploration.

Characterization of CD31 expression in CD4+ and CD8+T cell subpopulations in chronic untreated HIV infection.

BACKGROUND: The platelet endothelial cell adhesion molecule-1 (PECAM-1) or CD31 has been involved in regulation of T-cell tolerance, activation, survival and homing in mice cells. However, there is limited knowledge about the expression pattern and role of this molecule in human T cells, particularly in conditions of chronic immune activation. OBJECTIVES: We explored CD31 expression in T cell differentiation subsets of individuals with untreated HIV infection and in non-HIV-infected controls. We also assessed phenotypic differences between CD31+ and CD31- subsets in memory and terminally differentiated (TEMRA) CD4+ and CD8 + T cells. METHODS: Forty-one individuals with untreated HIV infection and 34 non-HIV-infected controls were included in the study. We compared the expression of CD31 in CD4+ and CD8 + T cells across stages of differentiation in the two study groups by flow cytometry. We also analyzed the expression of CD57 (a marker of senescence), Ki67 (a marker of cycling cells), PD-1 (a marker of exhaustion), and CD38/HLA-DR (a marker of immune activation) on memory and TEMRA CD31+ and CD31- T cells. RESULTS: CD31 expression was significantly higher in CD8 + T cells than in CD4 + T cells, measured as frequency, absolute numbers and median fluorescence intensity (MFI), in both study groups (p < 0.0001 in all cases). Intermediate differentiation subsets of CD4+ and CD8 + T cells expressed higher levels of CD31 in the context of HIV infection (p < 0.001 in all cases). CD31 expression frequency decreased with cellular differentiation of CD4+ and CD8 + T cells in both groups, but this decrease was steeper in individuals without HIV infection (CD4+: p < 0.001 and CD8+: p < 0.0001). As expected, memory and TEMRA CD4+ and CD8 + T cells expressed significantly higher levels of CD57, PD-1, Ki67 and CD38/HLA-DR in HIV-infected compared to non-HIV-infected individuals (p < 0.01 in all cases). CD31 expression was associated with lower activation of memory (but not TEMRA) CD4 + T cells in non-HIV-infected persons, an effect not observed in the HIV-infected group. CD31 expression on memory CD8 + T cells of HIV-infected individuals was associated higher levels of PD-1 (p = 0.0019) and CD38/HLADR (p = 0.0345), and higher PD-1 expression on CD8 + TEMRA (p = 0.0024), an effect not observed in non-HIV-infected individuals. CONCLUSION: In the context of HIV-associated chronic immune activation, specifically on memory CD8 + T cells, CD31 expression was associated with higher PD-1 and CD38/HLA-DR co-expression, suggesting that CD31 expression may result from an insufficient attempt to contain T cell exhaustion and activation. CD31-targeted therapies may contribute to modulate these cellular responses.

Epistasis of polymorphisms related to the articular cartilage extracellular matrix in knee osteoarthritis: Analysis-based multifactor dimensionality reduction.

Osteoarthritis (OA) is a complex disease with a multifactorial etiology. The genetic component is one of the main associated factors, resulting from interactions between genes and environmental factors. The aim of this study was to identify gene-gene interactions (epistasis) of the articular cartilage extracellular matrix (ECM) in knee OA. Ninety-two knee OA patients and 147 healthy individuals were included. Participants were genotyped in order to evaluate nine variants of eight genes associated with ECM metabolism using the OpenArray technology. Epistasis was analyzed using the multifactor dimensionality reduction (MDR) method. The MDR analysis showed significant gene-gene interactions between MMP3 (rs679620) and COL3A1 (rs1800255), and between COL3A1 (rs1800255) and VEGFA (rs699947) polymorphisms, with information gain values of 3.21% and 2.34%, respectively. Furthermore, in our study we found interactions in high-risk genotypes of the HIF1AN, MMP3 and COL3A1 genes; the most representative were [AA+CC+GA], [AA+CT+GA] and [AA+CT+GG], respectively; and low-risk genotypes [AA+CC+GG], [GG+TT+GA] and [AA+TT+GA], respectively. Knowing the interactions of these polymorphisms involved in articular cartilage ECM metabolism could provide a new tool to identify individuals at high risk of developing knee OA.

Risk of Wnt/ß-catenin signalling pathway gene polymorphisms in primary Sjögren's syndrome.

OBJECTIVE: To explore genetic polymorphisms of the Wnt/ß-catenin signalling pathway in primary SS (PSS). METHODS: We included 98 patients with PSS and 165 healthy volunteers. Genomic DNA was extracted from peripheral blood samples. Through an open-array platform of low density, we genotyped 25 polymorphisms from 14 genes (WISP1, DKK1, SOST, FRZB, LRP1, LRP4, LRP5, LRP6, GSKB, ADAMTS5, GDF5, FMN2, ADIPOQ and COL11A1) involved in the Wnt/ß-catenin signalling pathway. We compared the allelic and genotypic frequencies with Fisher's exact test and logistic regression analysis adjusted by age, gender and individual admixture, as well as bootstrap-resampling analysis. We assessed the gene-gene interaction by the multifactor dimensionality reduction method. RESULTS: We found a positive significant association with four polymorphisms: LRP5 rs606989, FRZB rs409238, GSK3B rs2037547 and ADIPOQ rs2241766. All of them conferred risk for PSS, being the highest among subjects carrying three to four risk alleles (P < 0.001). According to a multifactor dimensionality reduction analysis, the best models included the LRP5 (rs606989), FRZB (rs409238) and ADIPOQ (rs2241766) polymorphisms. CONCLUSION: LRP5, FRZB and ADIPOQ genes related in the Wnt/ß-catenin signalling pathway increased the risk of PSS. Further research is needed to establish their functional role in this clinical entity.

Associations between recent thymic emigrants and CD4+ T-cell recovery after short-term antiretroviral therapy initiation.

OBJECTIVE: Around 20-30% of HIV-infected individuals (HIV+) on successful antiretroviral therapy (ART) fail to normalize their CD4 T-cell counts. Various factors could contribute to the lack of immune reconstitution, one of them being thymic insufficiency. We aimed to explore associations between recent thymic emigrants (RTEs) and CD4 T-cell recovery. DESIGN: ART-naive HIV+ individuals who started ART with advanced AIDS were selected. Good versus poor immune reconstitution was defined by CD4 gains above or below 100 CD4 T cells/µl. The follow-up period was 6 months. METHODS: Peripheral blood mononuclear cells were isolated and flow cytometry was used to characterize RTEs as the fraction of naive CD4 T cells expressing CD31, the platelet endothelial cell adhesion molecule. Markers of cellular activation, senescence, exhaustion and cycling were also assessed. RESULTS: After 6 months on ART, HIV+ individuals with good immune reconstitution had higher absolute numbers of RTEs, compared with those with poor immune reconstitution, and these strongly correlated with CD4 gains in those individuals with good immune reconstitution but not with poor immune reconstitution. We also found that CD8 T-cell immune activation decreased as early as 2 months post-ART initiation in individuals with good immune reconstitution, but only at month 6 post-ART in individuals with poor immune reconstitution. Levels of immune activation were inversely correlated with the absolute numbers of RTEs in both groups, but more strongly so in individuals with poor immune reconstitution. CONCLUSION: We show that RTEs are linked to CD4 T-cell recovery and that the degree of immune reconstitution is not directly linked to persistent immune activation.

Impact of the gene-gene interactions related to the HIF-1α signaling pathway with the knee osteoarthritis development.

INTRODUCTION/OBJECTIVES: Articular cartilage is the target tissue of osteoarthritis (OA), and because it lacks capillary networks, the microenvironment is hypoxic. Hypoxia inducible factor-1 alpha (HIF-1α) regulates the homeostasis of this tissue. The aim of this study was to investigate whether genetic polymorphisms of the HIF-1α signaling pathway are involved in the development of knee OA. METHOD: We performed a case-control association study and genotyped 134 knee OA patients and 267 healthy controls. All participants were genotyped in order to evaluate 42 SNPs from 22 genes involved in the HIF-1α signaling pathway using the OpenArray technology. Gene-gene interactions (epistasis) were analyzed using the multifactor dimensionality reduction (MDR) method. RESULTS: The MDR analysis showed epistasis between AKT2 (rs8100018) and IGF1 (rs2288377), AKT2 (rs8100018) and IGF1 (rs35767), IGF1 (rs35767) and COL2A1 (rs1793953), and between GSK3B (rs6438552) and IGF1 (rs35767) polymorphisms, with information gain values of 21.24%, 8.37%, 9.93%, and 5.73%, respectively. Additionally, our model allowed us to identify high- and low-risk genotypes among COL2A1 rs1793953, GSK3B rs6438552, AKT2 rs8100018, and IGF1 rs35767 polymorphisms. CONCLUSIONS: Knowing the interactions of these polymorphisms involved in HIF-1α signaling pathway could provide a new diagnostic support tool to identify individuals at high risk of developing knee OA.

Common gene variants interactions related to uric acid transport are associated with knee osteoarthritis susceptibility.

BACKGROUND: The presence of genetic variants in uric acid (UA) transporters can be associated with hyperuricemia, and therefore with an increased risk of monosodium urate (MSU) crystal precipitation. The inflammatory process triggered by these crystals leads to cartilage damage, which, in turn, could promote knee osteoarthritis (KOA). OBJECTIVE: To determine whether genetic polymorphisms of the UA transporters and their interactions are associated with KOA. MATERIALS AND METHODS: Two hundred forty-three unrelated Mexican-mestizo individuals were recruited for this case-control study. Ninety-three of them were KOA patients but without gout, and one hundred and fifty healthy individuals with no symptoms or signs of KOA were recruited as controls. Forty-one single-nucleotide polymorphisms (SNPs) involved in the UA transporters were genotyped with OpenArray technology in a QuantStudio 12K flex-System with both cases and controls. RESULTS: After adjusting by age, gender, BMI, and ancestry, significant associations were found for eight SNPs: rs1260326 (GCKR), rs780093 (GCKR), rs17050272 (INHBB), rs1471633 (PDZK1), rs12129861 (PDZK1), rs7193778 (IGF1R), rs17786744 (STC1), and rs1106766 (R3HDM2). With respect to gene-gene interactions, the pairwise interactions of rs112129861 (PDZK1) and rs7193778 (IGF1R); rs17050272 (INHBB) and rs1106766 (R3HDM2); rs1106766 (R3HDM2) and rs780093 (GCKR); rs1260326 (GCKR) and rs17786744 (STC1); and rs17786744 (STC1) and rs1106766 (R3HDM2) make it possible to visualize the synergistic or antagonistic effect of their genotypes or alleles on KOA development. CONCLUSIONS: Our preliminary results show that the common gene variants related to UA transport are associated with KOA in the Mexican population. Further studies must be carried out to corroborate it.

Evolutionary history and spatiotemporal dynamics of the HIV-1 subtype B epidemic in Guatemala.

Different explanations exist on how HIV-1 subtype B spread in Central America, but the role of Guatemala, the Central American country with the highest number of people living with the virus, in this scenario is unknown. We investigated the evolutionary history and spatiotemporal dynamics of HIV-1 subtype B in Guatemala. A total of 1,047 HIV-1 subtype B pol sequences, from newly diagnosed ART-naïve, HIV-infected Guatemalan subjects enrolled between 2011 and 2013 were combined with published subtype B sequences from other Central American countries (n = 2,101) and with reference sequences representative of the BPANDEMIC and BCAR lineages from the United States (n = 465), France (n = 344) and the Caribbean (n = 238). Estimates of evolutionary, demographic, and phylogeographic parameters were obtained from sequence data using maximum likelihood and Bayesian coalescent-based methods. The majority of Guatemalan sequences (98.9%) belonged to the BPANDEMIC clade, and 75.2% of these sequences branched within 10 monophyletic clades: four also included sequences from other Central American countries (BCAM-I to BCAM-IV) and six were mostly (>99%) composed by Guatemalan sequences (BGU clades). Most clades mainly comprised sequences from heterosexual individuals. Bayesian coalescent-based analyses suggested that BGU clades originated during the 1990s and 2000s, whereas BCAM clades originated between the late 1970s and mid 1980s. The major hub of dissemination of all BGU, and of BCAM-II, and BCAM-IV clades was traced to the Department of Guatemala, while the root location of BCAM-I and BCAM-III was traced to Honduras. Most Guatemalan clades experienced initial phases of exponential growth (0.23 and 3.6 year-1), followed by recent growth declines. Our observations suggest that the Guatemalan HIV-1 subtype B epidemic is driven by dissemination of multiple BPANDEMIC founder viral strains, some restricted to Guatemala and others widely disseminated in the Central American region, with Guatemala City identified as a major hub of viral dissemination. Our results also suggest the existence of different sub-epidemics within Guatemala for which different targeted prevention efforts might be needed.

Gene-gene interactions of the Wnt/ß-catenin signaling pathway in knee osteoarthritis.

This study was designed to investigate whether genetic polymorphisms of the Wnt/ß-catenin signaling pathway and its interactions are involved in the development of knee osteoarthritis (KOA). Patients with KOA (n = 131) and healthy individuals (n = 190) with different ancestry from two Mexican populations (Mexico City and Guadalajara City) were analyzed. Twenty-five SNPs from thirteen genes (WISP1, DKK1, SOST, FRZB, LRP1, LRP4, LRP5, LRP6, GSKB, ADAMTS5, GDF5, FMN2 and COL11A1) involved in the Wnt/ß-catenin signaling pathway were genotyped. Genetic and allelic frequencies and gene-gene interactions were performed for this study. After adjusting for age, sex, BMI and admixture, significant associations were found for five SNPs in Mexico City: LRP6 rs12314259 (G/G genotype OR 0.22, P = 0.029; and G allele OR 0.48, P = 0.022), SOST rs851054 (C/T genotype OR 0.42, P = 0.027; and T allele OR 0.62, P = 0.026), FMN2 rs986690 (G/A genotype OR 0.42, P = 0.034; and A allele OR 0.50, P = 0.015), FRZB rs409238 (A/G genotype, OR 2.41, P = 0.022), and COL11A1 rs2615977 (A/C genotype OR 2.39, P = 0.024); no associations for Guadalajara City were found. With respect to gene-gene interactions, the pairwise interactions of WISP1-COL11A1, COL11A1-FRZB, FRZB-SOST and WISP1-FMN2 make it possible to visualize the synergistic or antagonistic effect of their genotypes or alleles in both populations. These results suggest that gene-gene interactions in the Wnt/ß-catenin signaling pathway play a role in the etiology of KOA.

Novel HLA class I associations with HIV-1 control in a unique genetically admixed population.

Associations between HLA class I alleles and HIV progression in populations exhibiting Amerindian and Caucasian genetic admixture remain understudied. Using univariable and multivariable analyses we evaluated HLA associations with five HIV clinical parameters in 3,213 HIV clade B-infected, ART-naïve individuals from Mexico and Central America (MEX/CAM cohort). A Canadian cohort (HOMER, n = 1622) was used for comparison. As expected, HLA allele frequencies in MEX/CAM and HOMER differed markedly. In MEX/CAM, 13 HLA-A, 24 HLA-B, and 14 HLA-C alleles were significantly associated with at least one clinical parameter. These included previously described protective (e.g. B*27:05, B*57:01/02/03 and B*58:01) and risk (e.g. B*35:02) alleles, as well as novel ones (e.g. A*03:01, B*15:39 and B*39:02 identified as protective, and A*68:03/05, B*15:30, B*35:12/14, B*39:01/06, B*39:05~C*07:02, and B*40:01~C*03:04 identified as risk). Interestingly, both protective (e.g. B*39:02) and risk (e.g. B*39:01/05/06) subtypes were identified within the common and genetically diverse HLA-B*39 allele group, characteristic to Amerindian populations. While HLA-HIV associations identified in MEX and CAM separately were similar overall (Spearman's rho = 0.33, p = 0.03), region-specific associations were also noted. The identification of both canonical and novel HLA/HIV associations provides a first step towards improved understanding of HIV immune control among unique and understudied Mestizo populations.

Weaker HLA Footprints on HIV in the Unique and Highly Genetically Admixed Host Population of Mexico.

HIV circumvents HLA class I-restricted CD8+ T-cell responses through selection of escape mutations that leave characteristic mutational "footprints," also known as HLA-associated polymorphisms (HAPs), on HIV sequences at the population level. While many HLA footprints are universal across HIV subtypes and human populations, others can be region specific as a result of the unique immunogenetic background of each host population. Using a published probabilistic phylogenetically informed model, we compared HAPs in HIV Gag and Pol (PR-RT) in 1,612 subtype B-infected, antiretroviral treatment-naive individuals from Mexico and 1,641 individuals from Canada/United States. A total of 252 HLA class I allele subtypes were represented, including 140 observed in both cohorts, 67 unique to Mexico, and 45 unique to Canada/United States. At the predefined statistical threshold of a q value of <0.2, 358 HAPs (201 in Gag, 157 in PR-RT) were identified in Mexico, while 905 (534 in Gag and 371 in PR-RT) were identified in Canada/United States. HAPs identified in Mexico included both canonical HLA-associated escape pathways and novel associations, in particular with HLA alleles enriched in Amerindian and mestizo populations. Remarkably, HLA footprints on HIV in Mexico were not only fewer but also, on average, significantly weaker than those in Canada/United States, although some exceptions were noted. Moreover, exploratory analyses suggested that the weaker HLA footprint on HIV in Mexico may be due, at least in part, to weaker and/or less reproducible HLA-mediated immune pressures on HIV in this population. The implications of these differences for natural and vaccine-induced anti-HIV immunity merit further investigation.IMPORTANCE HLA footprints on HIV identify viral regions under intense and consistent pressure by HLA-restricted immune responses and the common mutational pathways that HIV uses to evade them. In particular, HLA footprints can identify novel immunogenic regions and/or epitopes targeted by understudied HLA alleles; moreover, comparative analyses across immunogenetically distinct populations can illuminate the extent to which HIV immunogenic regions and escape pathways are shared versus population-specific pathways, information which can in turn inform the design of universal or geographically tailored HIV vaccines. We compared HLA-associated footprints on HIV in two immunogenetically distinct North American populations, those of Mexico and Canada/United States. We identify both shared and population-specific pathways of HIV adaptation but also make the surprising observation that HLA footprints on HIV in Mexico overall are fewer and weaker than those in Canada/United States, raising the possibility that HLA-restricted antiviral immune responses in Mexico are weaker, and/or escape pathways somewhat less consistent, than those in other populations.

Killer cell immunoglobulin-like receptor and human leukocyte antigen gene profiles in a cohort of HIV-infected Mexican Mestizos.

Killer cell immunoglobulin-like receptors (KIRs) represent the most polymorphic genes responsible for natural killer cell function, while human leukocyte antigen (HLA) class I molecules define and restrict cytotoxic T lymphocyte responses. Specific KIR, HLA, or KIR-HLA combinations have been implicated in the outcome of human immunodeficiency virus (HIV) disease. The remarkable polymorphism of KIR and HLA genes warrants descriptive gene frequency studies in different populations, as well as their impact on HIV disease progression in different immunogenetic contexts. We report KIR and HLA class I gene profiles of 511 unrelated HIV-infected Mexican Mestizo individuals from 18 states for whom genetic ancestry proportions were assessed. KIR and HLA gene profiles were compared between individuals from the north and central-south regions of the country and between individuals with higher European (EUR) or Amerindian (AMI) genetic ancestry component. A total of 65 KIR genotypes were observed, 11 harboring novel KIR gene combinations. A total of 164 HLA alleles were observed: 43 HLA-A, 87 HLA-B, and 34 HLA-C. Differences in the distribution of 12 HLA alleles were observed between individuals with higher AMI or EUR ancestry components (p < 0.05, q < 0.2). After correcting for genetic ancestry, only individual HLA alleles were associated with HIV disease progression, including a novel association with A*02:06, an Amerindian HLA allele associated with lower CD4+ T cell counts. No KIR effects were significant. Our results highlight the advantages of considering a detailed genetic stratification within populations when studying genetic profiles that could be implicated in disease-association studies.

HIV-1 drug resistance surveillance in antiretroviral treatment-naive individuals from a reference hospital in Guatemala, 2010-2013.

The recent expansion of antiretroviral treatment (ART) coverage in middle/low-income countries has been associated with increasing prevalence of HIV pre-ART drug resistance (PDR). We assessed PDR prevalence, patterns, and trends in Guatemala. Blood samples from 1,084 ART-naive individuals, enrolled from October 2010 to December 2013 at the Roosevelt Hospital in Guatemala City, were obtained. PDR was evaluated using the WHO mutation list for transmitted drug resistance (TDR) surveillance. An overall PDR prevalence of 7.3% (95% CI 5.8-9.0%) was observed for the whole study period. TDR to nonnucleoside reverse transcriptase inhibitors (NNRTI) was the highest (4.9%, p<0.001), followed by nucleoside RT inhibitors (1.8%) and protease inhibitors (1.0%). No significant trends in PDR prevalence were observed during the study period. However, higher NNRTI PDR levels were found in individuals with >500 and 350-500 CD4(+) T cells/µl (7.4% and 8.7%, respectively) compared to individuals with <350 CD4(+) T cells/µl (3.7%; p=0.039 and p=0.007, respectively), as well as a tendency of higher levels of NNRTI transmitted drug resistance (DR) in individuals with recent infection determined by HIV incidence tests (9.7%), suggesting increasing trends in time. Clusters of viruses with NNRTI PDR suggesting complex transmission networks were observed. No associations between PDR and demographic variables were found. PDR in Guatemala remains at an intermediate level. Nevertheless, we have shown evidence suggesting increasing trends in NNRTI PDR, which need to be taken into account in national HIV management policies.

Impact of HLA-B*35 subtype differences on HIV disease outcome in Mexico.

HLA-B35 has consistently been associated with rapid HIV disease progression, particularly alleles of the Px group. As B35 is the most prevalent HLA-B in Mexico, we investigated HIV disease outcome in relation to HLA expression in a large cohort (n=976) of Mexicans. Contrary to the previous studies, no impact on viral load or CD4 cell count was observed in association with the B35 PY/Px groups. However, we observed differences in HIV disease outcome associated with specific HLA-B35 alleles.

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.

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

BACKGROUND: Mounting 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. METHODS: The 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. RESULTS: Forty 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. CONCLUSION: Our 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.