Search

Background and objectives Public health officials faced with a large number of transmission clusters require a rapid, scalable and unbiased way to prioritize distribution of limited resources to maximize benefits. We hypothesize that transmission cluster prioritization based on phylogenetically derived lineage-level diversification rates will perform as well as or better than commonly used growth-based prioritization measures, without need for historical data or subjective interpretation. Methodology 9822 HIV pol sequences collected during routine drug resistance genotyping were used alongside simulated sequence data to infer sets of phylogenetic transmission clusters via patristic distance threshold. Prioritized clusters inferred from empirical data were compared to those prioritized by the current public health protocols. Prioritization of simulated clusters was evaluated based on correlation of a given prioritization measure with future cluster growth, as well as the number of direct downstream transmissions from cluster members. Results Empirical data suggest diversification rate-based measures perform comparably to growth-based measures in recreating public heath prioritization choices. However, unbiased simulated data reveals phylogenetic diversification rate-based measures perform better in predicting future cluster growth relative to growth-based measures, particularly long-term growth. Diversification rate-based measures also display advantages over growth-based measures in highlighting groups with greater future transmission events compared to random groups of the same size. Furthermore, diversification rate measures were notably more robust to effects of decreased sampling proportion. Conclusions and implications Our findings indicate diversification rate-based measures frequently outperform growth-based measures in predicting future cluster growth and offer several additional advantages beneficial to optimizing the public health prioritization process.

Active surveillance of influenza A viruses of swine (IAV-S) involving 262 farms and 10 slaughterhouses in seven provinces in northern and southern Vietnam from 2010 to 2015 yielded 388 isolates from 32 farms; these viruses were classified into H1N1, H1N2, and H3N2 subtypes. Whole-genome sequencing followed by phylogenetic analysis revealed that the isolates represented 15 genotypes, according to the genetic constellation of the eight segments. All of the H1N1 viruses were entirely A(H1N1)pdm09 viruses, whereas all of the H1N2 and H3N2 viruses were reassortants among 5 distinct ancestral viruses: H1 and H3 triple-reassortant (TR) IAV-S that originated from North American pre-2009 human seasonal H1, human seasonal H3N2, and A(H1N1)pdm09 viruses. Notably, 93% of the reassortant IAV-S retained M genes that were derived from A(H1N1)pdm09, suggesting some advantage in terms of their host adaptation. Bayesian Markov chain Monte Carlo analysis revealed that multiple introductions of A(H1N1)pdm09 and TR IAV-S into the Vietnamese pig population have driven the genetic diversity of currently circulating Vietnamese IAV-S. In addition, our results indicate that a reassortant IAV-S with human-like H3 and N2 genes and an A(H1N1)pdm09 origin M gene likely caused a human case in Ho Chi Minh City in 2010. Our current findings indicate that human-to-pig transmission as well as cocirculation of different IAV-S have contributed to diversifying the gene constellations of IAV-S in Vietnam. IMPORTANCE This comprehensive genetic characterization of 388 influenza A viruses of swine (IAV-S) isolated through active surveillance of Vietnamese pig farms from 2010 through 2015 provides molecular epidemiological insight into the genetic diversification of IAV-S in Vietnam after the emergence of A(H1N1)pdm09 viruses. Multiple reassortments among A(H1N1)pdm09 viruses and enzootic IAV-S yielded 14 genotypes, 9 of which carried novel gene combinations. The reassortants that carried M genes derived from A(H1N1)pdm09 viruses became predominant, replacing those of the IAV-S that had been endemic in Vietnam since 2011. Notably, one of the novel reassortants likely caused a human case in Vietnam. Given that Vietnam is the second-largest pig-producing country in Asia, continued monitoring of IAV-S is highly important from the viewpoints of both the swine industry and human public health.

Host and viral factors influence the HIV-1 infection course. Reduced Nef function has been observed in HIV-1 controllers during the chronic phase, but the kinetics and mechanisms of Nef attenuation in such individuals remain unclear. We examined plasma RNA-derived Nef clones from 10 recently infected individuals who subsequently suppressed viremia to less than 2,000 RNA copies/ml within 1 year postinfection (acute controllers) and 50 recently infected individuals who did not control viremia (acute progressors). Nef clones from acute controllers displayed a lesser ability to downregulate CD4 and HLA class I from the cell surface and a reduced ability to enhance virion infectivity compared to those from acute progressors (all P < 0.01). HLA class I downregulation activity correlated inversely with days postinfection (Spearman's R = −0.85, P = 0.004) and positively with baseline plasma viral load (Spearman's R = 0.81, P = 0.007) in acute controllers but not in acute progressors. Nef polymorphisms associated with functional changes over time were identified in follow-up samples from six controllers. For one such individual, mutational analyses indicated that four polymorphisms selected by HLA-A*31 and B*37 acted in combination to reduce Nef steady-state protein levels and HLA class I downregulation activity. Our results demonstrate that relative control of initial HIV-1 viremia is associated with Nef clones that display reduced function, which in turn may influence the course of HIV-1 infection. Transmission of impaired Nef sequences likely contributed in part to this observation; however, accumulation of HLA-associated polymorphisms in Nef that impair function also suggests that CD8 + T-cell pressures play a role in this phenomenon. IMPORTANCE Rare individuals can spontaneously control HIV-1 viremia in the absence of antiretroviral treatment. Understanding the host and viral factors that contribute to the controller phenotype may identify new strategies to design effective vaccines or therapeutics. The HIV-1 Nef protein enhances viral pathogenesis through multiple mechanisms. We examined the function of plasma HIV-1 RNA-derived Nef clones isolated from 10 recently infected individuals who subsequently controlled HIV viremia compared to the function of those from 50 individuals who failed to control viremia. Our results demonstrate that early Nef clones from HIV controllers displayed lower HLA class I and CD4 downregulation activity, as well as a reduced ability to enhance virion infectivity. The accumulation of HLA-associated polymorphisms in Nef during the first year postinfection was associated with impaired protein function in some controllers. This report highlights the potential for host immune responses to modulate HIV pathogenicity and disease outcome by targeting cytotoxic T lymphocyte (CTL) epitopes in Nef.

The extent to which HIV-1 clade B strains exhibit population-specific adaptations to host HLA alleles remains incompletely known, in part due to incomplete characterization of HLA-associated HIV-1 polymorphisms (HLA-APs) in different global populations. Moreover, it remains unknown to what extent the same HLA alleles may drive significantly different escape pathways across populations. As the Japanese population exhibits distinctive HLA class I allele distributions, comparative analysis of HLA-APs between HIV-1 clade B-infected Japanese and non-Asian cohorts could shed light on these questions. However, HLA-APs remain incompletely mapped in Japan. In a cohort of 430 treatment-naive Japanese with chronic HIV-1 clade B infection, we identified 284 HLA-APs in Gag, Pol, and Nef using phylogenetically corrected methods. The number of HLA-associated substitutions in Pol, notably those restricted by HLA-B*52:01, was weakly inversely correlated with the plasma viral load (pVL), suggesting that the transmission and persistence of B*52:01-driven Pol mutations could modulate the pVL. Differential selection of HLA-APs between HLA subtype members, including those differing only with respect to substitutions outside the peptide-binding groove, was observed, meriting further investigation as to their mechanisms of selection. Notably, two-thirds of HLA-APs identified in Japan had not been reported in previous studies of predominantly Caucasian cohorts and were attributable to HLA alleles unique to, or enriched in, Japan. We also identified 71 cases where the same HLA allele drove significantly different escape pathways in Japan versus predominantly Caucasian cohorts. Our results underscore the distinct global evolution of HIV-1 clade B as a result of host population-specific cellular immune pressures. IMPORTANCE Cytotoxic T lymphocyte (CTL) escape mutations in HIV-1 are broadly predictable based on the HLA class I alleles expressed by the host. Because HLA allele distributions differ among worldwide populations, the pattern and diversity of HLA-associated escape mutations are likely to be somewhat distinct to each race and region. HLA-associated polymorphisms (HLA-APs) in HIV-1 have previously been identified at the population level in European, North American, Australian, and African cohorts; however, large-scale analyses of HIV-1 clade B-specific HLA-APs in Asians are lacking. Differential intraclade HIV-1 adaptation to global populations can be investigated via comparative analyses of HLA-associated polymorphisms across ethnic groups, but such studies are rare. Here, we identify HLA-APs in a large Japanese HIV-1 clade B cohort using phylogenetically informed methods and observe that the majority of them had not been previously characterized in predominantly Caucasian populations. The results highlight HIV's unique adaptation to cellular immune pressures imposed by different global populations.

Recent studies indicate that murine Tregs highly express the ENTDP1, as well as the 5′-NT and thereby, suppress Teff function by extracellular adenosine production. Furthermore, CD73 seems to play a role as costimulatory molecule for T cell differentiation. In this study, we analyzed the expression of CD73 on peripheral and lymph nodal Teffs and Tregs in a cohort of 95 HIV patients at different stages of disease, including LTNP and ECs. In contrast to murine Tregs, CD73 was only expressed on a small minority (∼10%) of peripheral Tregs. In contrast, we see high expression of CD73 on peripheral CD8+ T cells. In HIV infection, CD73 is markedly reduced on all Teffs and Tregs, regardless of the memory subtype. On CD8+ T cells, a positive correlation between CD73 expression and CD4 counts (P=0.0003) was detected. CD73 expression on CD8+ T cells negatively correlated with HLA-DR (

Background The highly genetically diverse HIV-1 group M subtypes may differ in their biological properties. Nef is an important mediator of viral pathogenicity; however, to date, a comprehensive inter-subtype comparison of Nef in vitro function has not been undertaken. Here, we investigate two of Nef’s most well-characterized activities, CD4 and HLA class I downregulation, for clones obtained from 360 chronic patients infected with HIV-1 subtypes A, B, C or D. Results Single HIV-1 plasma RNA Nef clones were obtained from N=360 antiretroviral-naïve, chronically infected patients from Africa and North America: 96 (subtype A), 93 (B), 85 (C), and 86 (D). Nef clones were expressed by transfection in an immortalized CD4+ T-cell line. CD4 and HLA class I surface levels were assessed by flow cytometry. Nef expression was verified by Western blot. Subset analyses and multivariable linear regression were used to adjust for differences in age, sex and clinical parameters between cohorts. Consensus HIV-1 subtype B and C Nef sequences were synthesized and functionally assessed. Exploratory sequence analyses were performed to identify potential genotypic correlates of Nef function. Subtype B Nef clones displayed marginally greater CD4 downregulation activity (p = 0.03) and markedly greater HLA class I downregulation activity (p in vitro functionality. Inter-subtype differences in HLA class I downregulation remained statistically significant after controlling for differences in age, sex, and clinical parameters (p Conclusions Results support a functional hierarchy of subtype B > A/D > C for Nef-mediated CD4 and HLA class I downregulation. The mechanisms underlying these differences and their relevance to HIV-1 pathogenicity merit further investigation.

Human leukocyte antigen (HLA) class I-associated polymorphisms in HIV-1 that persist upon transmission to HLA-mismatched hosts may spread in the population as the epidemic progresses. Transmission of HIV-1 sequences containing such adaptations may undermine cellular immune responses to the incoming virus in future hosts. Building upon previous work, we investigated the extent of HLA-associated polymorphism accumulation in HIV-1 polymerase (Pol) through comparative analysis of linked HIV-1/HLA class I genotypes sampled during historic (1979 to 1989; n = 338) and modern (2001 to 2011; n = 278) eras from across North America (Vancouver, BC, Canada; Boston, MA; New York, NY; and San Francisco, CA). Phylogenies inferred from historic and modern HIV-1 Pol sequences were star-like in shape, with an inferred most recent common ancestor (epidemic founder virus) sequence nearly identical to the modern North American subtype B consensus sequence. Nevertheless, modern HIV-1 Pol sequences exhibited roughly 2-fold-higher patristic (tip-to-tip) genetic distances than historic sequences, with HLA pressures likely driving ongoing diversification. Moreover, the frequencies of published HLA-associated polymorphisms in individuals lacking the selecting HLA class I allele was on average ∼2.5-fold higher in the modern than in the historic era, supporting their spread in circulation, though some remained stable in frequency during this time. Notably, polymorphisms restricted by protective HLA alleles appear to be spreading to a greater relative extent than others, though these increases are generally of modest absolute magnitude. However, despite evidence of polymorphism spread, North American hosts generally remain at relatively low risk of acquiring an HIV-1 polymerase sequence substantially preadapted to their HLA profiles, even in the present era. IMPORTANCE HLA class I-restricted cytotoxic T-lymphocyte (CTL) escape mutations in HIV-1 that persist upon transmission may accumulate in circulation over time, potentially undermining host antiviral immunity to the transmitted viral strain. We studied >600 experimentally collected HIV-1 polymerase sequences linked to host HLA information dating back to 1979, along with phylogenetically reconstructed HIV-1 sequences dating back to the virus' introduction into North America. Overall, our results support the gradual spread of many—though not all—HIV-1 polymerase immune escape mutations in circulation over time. This is consistent with recent observations from other global regions, though the extent of polymorphism accumulation in North America appears to be lower than in populations with high seroprevalence, older epidemics, and/or limited HLA diversity. Importantly, the risk of acquiring an HIV-1 polymerase sequence at transmission that is substantially preadapted to one's HLA profile remains relatively low in North America, even in the present era.

An attenuation of the HIV-1 replication capacity (RC) has been observed for immune-mediated escape mutations in Gag restricted by protective HLA alleles. However, the extent to which escape mutations affect other viral proteins during natural infection is not well understood. We generated recombinant viruses encoding plasma HIV-1 RNA integrase sequences from antiretroviral-naïve individuals with early ( n = 88) and chronic ( n = 304) infections and measured the in vitro RC of each. In contrast to data from previous studies of Gag, we observed little evidence that host HLA allele expression was associated with integrase RC. A modest negative correlation was observed between the number of HLA-B-associated integrase polymorphisms and RC in chronic infection ( R = −0.2; P = 0.003); however, this effect was not driven by mutations restricted by protective HLA alleles. Notably, the integrase variants S119R, G163E, and I220L, which represent uncommon polymorphisms associated with HLA-C*05, -A*33, and -B*52, respectively, correlated with lower RC (all q < 0.2). We identified a novel C*05-restricted epitope (HTDNGSNF 114–121 ) that likely contributes to the selection of the S119R variant, the polymorphism most significantly associated with lower RC in patient sequences. An NL4-3 mutant encoding the S119R polymorphism displayed a ∼35%-reduced function that was rescued by a single compensatory mutation of A91E. Together, these data indicate that substantial HLA-driven attenuation of integrase is not a general phenomenon during HIV-1 adaptation to host immunity. However, uncommon polymorphisms selected by HLA alleles that are not conventionally regarded to be protective may be associated with impaired protein function. Vulnerable epitopes in integrase might therefore be considered for future vaccine strategies.

Rare individuals homozygous for a naturally-occurring 32 base pair deletion in the CCR5 gene (CCR5∆32/∆32) are resistant to infection by CCR5-using (“R5”) HIV-1 strains but remain susceptible to less common CXCR4-using (“X4”) strains. The evolutionary dynamics of X4 infections however, remain incompletely understood. We identified two individuals, one CCR5wt/wt and one CCR5∆32/∆32, within the Vancouver Injection Drug Users Study who were infected with a genetically similar X4 HIV-1 strain. While early-stage plasma viral loads were comparable in the two individuals (~4.5–5 log10 HIV-1 RNA copies/ml), CD4 counts in the CCR5wt/wt individual reached a nadir of 3 within 17 months but remained >250 cells/mm3 in the CCR5∆32/∆32 individual. Ancestral phylogenetic reconstructions using longitudinal envelope-V3 deep sequences suggested that both individuals were infected by a single transmitted/founder (T/F) X4 virus that differed at only one V3 site (codon 24). While substantial within-host HIV-1 V3 diversification was observed in plasma and PBMC in both individuals, the CCR5wt/wt individual’s HIV-1 population gradually reverted from 100% X4 to ~60% R5 over ~4 years whereas the CCR5∆32/∆32 individual’s remained consistently X4. Our observations illuminate early dynamics of X4 HIV-1 infections and underscore the influence of CCR5 genotype on HIV-1 V3 evolution.

Nef plays a major role in HIV-1 pathogenicity. We studied HIV-1 subtype C infected individuals in acute/early (n = 120) or chronic (n = 207) infection to investigate the relationship between Nef-mediated CD4/HLA-I down-regulation activities and disease progression, and the influence of immune-driven sequence variation on these Nef functions. A single Nef sequence per individual was cloned into an expression plasmid, followed by transfection of a T cell line and measurement of CD4 and HLA-I expression. In early infection, a trend of higher CD4 down-regulation ability correlating with higher viral load set point was observed (r = 0.19, p = 0.05), and higher HLA-I down-regulation activity was significantly associated with faster rate of CD4 decline (p = 0.02). HLA-I down-regulation function correlated inversely with the number HLA-associated polymorphisms previously associated with reversion in the absence of the selecting HLA allele (r = -0.21, p = 0.0002). These data support consideration of certain Nef regions in HIV-1 vaccine strategies designed to attenuate the infection course.

HLA-restricted immune escape mutations that persist following HIV transmission could gradually spread through the viral population, thereby compromising host antiviral immunity as the epidemic progresses. To assess the extent and phenotypic impact of this phenomenon in an immunogenetically diverse population, we genotypically and functionally compared linked HLA and HIV (Gag/Nef) sequences from 358 historic (1979–1989) and 382 modern (2000–2011) specimens from four key cities in the North American epidemic (New York, Boston, San Francisco, Vancouver). Inferred HIV phylogenies were star-like, with approximately two-fold greater mean pairwise distances in modern versus historic sequences. The reconstructed epidemic ancestral (founder) HIV sequence was essentially identical to the North American subtype B consensus. Consistent with gradual diversification of a “consensus-like” founder virus, the median “background” frequencies of individual HLA-associated polymorphisms in HIV (in individuals lacking the restricting HLA[s]) were ∼2-fold higher in modern versus historic HIV sequences, though these remained notably low overall (e.g. in Gag, medians were 3.7% in the 2000s versus 2.0% in the 1980s). HIV polymorphisms exhibiting the greatest relative spread were those restricted by protective HLAs. Despite these increases, when HIV sequences were analyzed as a whole, their total average burden of polymorphisms that were “pre-adapted” to the average host HLA profile was only ∼2% greater in modern versus historic eras. Furthermore, HLA-associated polymorphisms identified in historic HIV sequences were consistent with those detectable today, with none identified that could explain the few HIV codons where the inferred epidemic ancestor differed from the modern consensus. Results are therefore consistent with slow HIV adaptation to HLA, but at a rate unlikely to yield imminent negative implications for cellular immunity, at least in North America. Intriguingly, temporal changes in protein activity of patient-derived Nef (though not Gag) sequences were observed, suggesting functional implications of population-level HIV evolution on certain viral proteins., Author Summary Upon HIV transmission, many – though not all – immune escape mutations selected in the previous host will revert to the consensus residue. The persistence of certain escape mutations following transmission has led to concerns that these could gradually accumulate in circulating HIV sequences over time, thereby undermining host antiviral immune potential as the epidemic progresses. As certain immune-driven mutations reduce viral fitness, their spread through the population could also have consequences for the average replication capacity and/or protein function of circulating HIV sequences. Here, we characterized HIV sequences, linked to host immunogenetic information, from patients enrolled in historic (1979–1989) and modern (2000–2011) HIV cohorts from four key cities in the North American epidemic. We reconstructed the epidemic's ancestral (founder) HIV sequence and assessed the subsequent extent to which known HIV immune escape mutations have spread in the population. Our data support the gradual spread of many - though not all - immune escape mutations in HIV sequences over time, but to an extent that is unlikely to have major immediate immunologic consequences for the North American epidemic. Notably, in vitro assessments of ancestral and patient-derived HIV sequences suggested functional implications of ongoing HIV evolution for certain viral proteins.

Background The reproducible nature of HIV-1 escape from HLA-restricted CD8+ T-cell responses allows the identification of HLA-associated viral polymorphisms “at the population level” – that is, via analysis of cross-sectional, linked HLA/HIV-1 genotypes by statistical association. However, elucidating their timing of selection traditionally requires detailed longitudinal studies, which are challenging to undertake on a large scale. We investigate whether the extent and relative timecourse of immune-driven HIV adaptation can be inferred via comparative cross-sectional analysis of independent early and chronic infection cohorts. Results Similarly-powered datasets of linked HLA/HIV-1 genotypes from individuals with early (median 200/dataset), HLA class I and HIV-1 Gag/Pol/Nef diversity, were established. These datasets were first used to define a list of 162 known HLA-associated polymorphisms detectable at the population level in cohorts of the present size and host/viral genetic composition. Of these 162 known HLA-associated polymorphisms, 15% (occurring at 14 Gag, Pol and Nef codons) were already detectable via statistical association in the early infection dataset at p ≤ 0.01 (q

BACKGROUND Effective antiretroviral therapy (ART) dramatically reduces human immunodeficiency virus (HIV) transmission. However, isolated shedding of HIV type 1 (HIV-1) in semen (IHS) can occur in the absence of detectable viremia or genital infections. We hypothesized that ART intensification with medications active in semen might prevent IHS. METHODS Paired blood and semen samples were collected monthly for 6 months from HIV-infected men starting ART that was intensified (iART) with maraviroc and raltegravir in an open-label fashion. Semen parameters were compared to those of historical controls starting standard ART (sART). RESULTS Compared with 25 controls who started sART, the semen HIV-1 load in 13 subjects who started iART was more rapidly suppressed (P = .043). IHS was detected at >1 visit in 2 participants (15%) receiving iART and in 12 controls (48%) receiving sART (P = .040). Among iART recipients, IHS was associated with lower raltegravir concentrations in blood and semen, compared with complete HIV-1 suppression (P = .03). Prolonged, high-level IHS (ie, shedding of >5000 RNA copies/mL) was observed in 1 iART recipient (8%), despite rapid viremia suppression and therapeutic drug levels; for 10 months, this virus remained R5 tropic, drug susceptible, and similar in sequence to virus recovered from blood. IHS was not seen after >3 years of effective ART in a parallel, prospective cohort study. CONCLUSIONS iART transiently reduced the occurrence of IHS early after ART initiation but did not prevent high-level IHS. IHS was not seen after more prolonged sART.

HLA class I-associated polymorphisms identified at the population level mark viral sites under immune pressure by individual HLA alleles. As such, analysis of their distribution, frequency, location, statistical strength, sequence conservation, and other properties offers a unique perspective from which to identify correlates of protective cellular immunity. We analyzed HLA-associated HIV-1 subtype B polymorphisms in 1,888 treatment-naïve, chronically infected individuals using phylogenetically informed methods and identified characteristics of HLA-associated immune pressures that differentiate protective and nonprotective alleles. Over 2,100 HLA-associated HIV-1 polymorphisms were identified, approximately one-third of which occurred inside or within 3 residues of an optimally defined cytotoxic T-lymphocyte (CTL) epitope. Differential CTL escape patterns between closely related HLA alleles were common and increased with greater evolutionary distance between allele group members. Among 9-mer epitopes, mutations at HLA-specific anchor residues represented the most frequently detected escape type: these occurred nearly 2-fold more frequently than expected by chance and were computationally predicted to reduce peptide-HLA binding nearly 10-fold on average. Characteristics associated with protective HLA alleles (defined using hazard ratios for progression to AIDS from natural history cohorts) included the potential to mount broad immune selection pressures across all HIV-1 proteins except Nef, the tendency to drive multisite and/or anchor residue escape mutations within known CTL epitopes, and the ability to strongly select mutations in conserved regions within HIV's structural and functional proteins. Thus, the factors defining protective cellular immune responses may be more complex than simply targeting conserved viral regions. The results provide new information to guide vaccine design and immunogenicity studies.

Background: The extent to which HIV replication capacity (RC) has changed over the epidemic’s course, and the influence of HLA-associated immune pressure as its driving force remains unknown. We performed a comparative study of immune escape and RC in historic (1979-1989) and modern Gag subtype B sequences from North America. Methods: Using phylogenetically-informed methods, we identified HLA-associated Gag polymorphisms in a historic cohort (N=239; 1979-1989). We also generated recombinant NL4-3 viruses encoding clonal plasma RNA Gag from 80 historic and 58 modern (2002-2008) sequences. Viral RC was measured using a GFP reporter T-cell assay and results were normalized to NL4-3 controls. Results: 95% of HLA-associated polymorphisms identified in the historic cohort were consistent with published modern escape pathways. Overall, the prevalence of HLA-associated polymorphisms in the general population increased a median 1.3-fold between historic and modern sequences; however in many cases this was influenced by differences in HLA allele frequencies between HIV-infected populations examined. Of note, the prevalence of the B*27-associated R264K escape mutation increased from 0.4 to 1.3% in the general population over time despite B*27 allele frequency remaining constant at 2.5%. Modestly lower viral RC was observed for Gag recombinant viruses constructed from pre-1985 sequences (median 0.86 [IQR 0.78-0.97], N=24) compared to those from 1985-1989 (median 0.98 [IQR 0.87-1.05], N=56) and 2002-2008 (median 0.96 [IQR 0.83-0.1.10], N=58) (p=0.049). In both historic and modern cohorts, host expression of HLA-B*27 was associated with lower RC (p=0.007). Gag codons associated with lower RC, including S67A, were identified in an exploratory analysis. Conclusion: Gag-mediated viral RC may have increased modestly since the beginning of the North American epidemic, despite limited evidence for HLA-driven viral sequence evolution during this time. Although mechanisms driving RC differences remain unclear, results do not support rapid and substantial accumulation of HLA-driven escape mutations in circulating North American HIV-1 sequences.