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.

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.

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.