HLA class I-driven evolution of human immunodeficiency virus type 1 subtype c proteome: immune escape and viral load.

Human immunodeficiency virus type 1 (HIV-1) mutations that confer escape from cytotoxic T-lymphocyte (CTL) recognition can sometimes result in lower viral fitness. These mutations can then revert upon transmission to a new host in the absence of CTL-mediated immune selection pressure restricted by the HLA alleles of the prior host. To identify these potentially critical recognition points on the virus, we assessed HLA-driven viral evolution using three phylogenetic correction methods across full HIV-1 subtype C proteomes from a cohort of 261 South Africans and identified amino acids conferring either susceptibility or resistance to CTLs. A total of 558 CTL-susceptible and -resistant HLA-amino acid associations were identified and organized into 310 immunological sets (groups of individual associations related to a single HLA/epitope combination). Mutations away from seven susceptible residues, including four in Gag, were associated with lower plasma viral-RNA loads (q < 0.2 [where q is the expected false-discovery rate]) in individuals with the corresponding HLA alleles. The ratio of susceptible to resistant residues among those without the corresponding HLA alleles varied in the order Vpr > Gag > Rev > Pol > Nef > Vif > Tat > Env > Vpu (Fisher's exact test; P < or = 0.0009 for each comparison), suggesting the same ranking of fitness costs by genes associated with CTL escape. Significantly more HLA-B (chi(2); P = 3.59 x 10(-5)) and HLA-C (chi(2); P = 4.71 x 10(-6)) alleles were associated with amino acid changes than HLA-A, highlighting their importance in driving viral evolution. In conclusion, specific HIV-1 residues (enriched in Vpr, Gag, and Rev) and HLA alleles (particularly B and C) confer susceptibility to the CTL response and are likely to be important in the development of vaccines targeted to decrease the viral load.

Phylogenetic dependency networks: inferring patterns of CTL escape and codon covariation in HIV-1 Gag.

HIV avoids elimination by cytotoxic T-lymphocytes (CTLs) through the evolution of escape mutations. Although there is mounting evidence that these escape pathways are broadly consistent among individuals with similar human leukocyte antigen (HLA) class I alleles, previous population-based studies have been limited by the inability to simultaneously account for HIV codon covariation, linkage disequilibrium among HLA alleles, and the confounding effects of HIV phylogeny when attempting to identify HLA-associated viral evolution. We have developed a statistical model of evolution, called a phylogenetic dependency network, that accounts for these three sources of confounding and identifies the primary sources of selection pressure acting on each HIV codon. Using synthetic data, we demonstrate the utility of this approach for identifying sites of HLA-mediated selection pressure and codon evolution as well as the deleterious effects of failing to account for all three sources of confounding. We then apply our approach to a large, clinically-derived dataset of Gag p17 and p24 sequences from a multicenter cohort of 1144 HIV-infected individuals from British Columbia, Canada (predominantly HIV-1 clade B) and Durban, South Africa (predominantly HIV-1 clade C). The resulting phylogenetic dependency network is dense, containing 149 associations between HLA alleles and HIV codons and 1386 associations among HIV codons. These associations include the complete reconstruction of several recently defined escape and compensatory mutation pathways and agree with emerging data on patterns of epitope targeting. The phylogenetic dependency network adds to the growing body of literature suggesting that sites of escape, order of escape, and compensatory mutations are largely consistent even across different clades, although we also identify several differences between clades. As recent case studies have demonstrated, understanding both the complexity and the consistency of immune escape has important implications for CTL-based vaccine design. Phylogenetic dependency networks represent a major step toward systematically expanding our understanding of CTL escape to diverse populations and whole viral genes.

Racial differences in the evaluation and treatment of hepatitis C among veterans: a retrospective cohort study.

OBJECTIVES: We examined the association between race and hepatitis C virus (HCV) evaluation and treatment of veterans in the Northwest Network of the Department of Veterans Affairs (VA). METHODS: In our retrospective cohort study, we used medical records to determine antiviral treatment of 4263 HCV-infected patients from 8 VA medical centers. Secondary outcomes included specialty referrals, laboratory evaluation, viral genotype testing, and liver biopsy. Multiple logistic regression was used to adjust for clinical (measured through laboratory results and International Classification of Diseases, Ninth Revision, codes) and sociodemographic factors. RESULTS: Blacks were less than half as likely as Whites to receive antiviral treatment (odds ratio [OR]=0.38; 95% confidence interval [CI]=0.23, 0.63). Both had similar odds of referral and liver biopsy. However, Blacks were significantly less likely to have complete laboratory evaluation (OR=0.67; 95% CI=0.52, 0.88) and viral genotype testing (OR=0.68; 95% CI=0.51, 0.90). CONCLUSIONS: Race is associated with receipt of medical care for various medical conditions. Further investigation is warranted to help understand whether patient preference or provider bias may explain why HCV-infected Blacks were less likely to receive medical care than Whites.

Large-scale amplification, cloning and sequencing of near full-length HIV-1 subtype C genomes.

Full-length HIV-1 genome sequencing provides important data needed to address several vaccine design, molecular epidemiologic and pathogenesis questions. A protocol is presented for obtaining near full-length genomes (NFLGs) from subjects infected with HIV-1 subtype C. This protocol was used to amplify NFLGs from 244 of 366 (67%) samples collected at two clinics in Durban, South Africa (SK and PS). Viral load was directly associated with frequency of successful NFLG amplification for both cohorts (PS; p = 0.005 and SK; p < 0.001). Seventeen of 38 initially NFLG-negative SK samples had variation within the PCR primer binding sites, however only 3 of these were successfully re-amplified using re-designed primers homologous to the target viruses. NFLGs were obtained from 7 of 24 PBMC samples processed from subjects whose plasma did not yield a NFLG. Stable plasmid clones were obtained from all 244 NFLG-positive PCR products, and both strands of each genome were sequenced, using a primary set of 46 primers. These methods thus allow the large-scale collection of HIV-1 NFLGs from populations infected primarily with subtype C. The methods are readily adaptable to other HIV-1 subtypes, and provide materials for viral functional analyses and population-based molecular epidemiology studies that include analysis of viral genome chimerization.

Amino-acid co-variation in HIV-1 Gag subtype C: HLA-mediated selection pressure and compensatory dynamics.

BACKGROUND: Despite high potential for HIV-1 genetic variation, the emergence of some mutations is constrained by fitness costs, and may be associated with compensatory amino acid (AA) co-variation. To characterize the interplay between Cytotoxic T Lymphocyte (CTL)-mediated pressure and HIV-1 evolutionary pathways, we investigated AA co-variation in Gag sequences obtained from 449 South African individuals chronically infected with HIV-1 subtype C. METHODOLOGY/PRINCIPAL FINDINGS: Individuals with CTL responses biased toward Gag presented lower viral loads than individuals with under-represented Gag-specific CTL responses. Using methods that account for founder effects and HLA linkage disequilibrium, we identified 35 AA sites under Human Leukocyte Antigen (HLA)-restricted CTL selection pressure and 534 AA-to-AA interactions. Analysis of two-dimensional distances between co-varying residues revealed local stabilization mechanisms since 40% of associations involved neighboring residues. Key features of our co-variation analysis included sites with a high number of co-varying partners, such as HLA-associated sites, which had on average 55% more connections than other co-varying sites. CONCLUSIONS/SIGNIFICANCE: Clusters of co-varying AA around HLA-associated sites (especially at typically conserved sites) suggested that cooperative interactions act to preserve the local structural stability and protein function when CTL escape mutations occur. These results expose HLA-imprinted HIV-1 polymorphisms and their interlinked mutational paths in Gag that are likely due to opposite selective pressures from host CTL-mediated responses and viral fitness constraints.

Rare HLA drive additional HIV evolution compared to more frequent alleles.

HIV-1 can evolve HLA-specific escape variants in response to HLA-mediated cellular immunity. HLA alleles that are common in the host population may increase the frequency of such escape variants at the population level. When loss of viral fitness is caused by immune escape variation, these variants may revert upon infection of a new host who does not have the corresponding HLA allele. Furthermore, additional escape variants may appear in response to the nonconcordant HLA alleles. Because individuals with rare HLA alleles are less likely to be infected by a partner with concordant HLA alleles, viral populations infecting hosts with rare HLA alleles may undergo a greater amount of evolution than those infecting hosts with common alleles due to the loss of preexisting escape variants followed by new immune escape. This hypothesis was evaluated using maximum likelihood phylogenetic trees of each gene from 272 full-length HIV-1 sequences. Recent viral evolution, as measured by the external branch length, was found to be inversely associated with HLA frequency in nef (p < 0.02), env (p < 0.03), and pol (p < or = 0.05), suggesting that rare HLA alleles provide a disproportionate force driving viral evolution compared to common alleles, likely due to the loss of preexisting escape variants during early stages postinfection.

Broad and Gag-biased HIV-1 epitope repertoires are associated with lower viral loads.

BACKGROUND: HLA class-I alleles differ in their ability to control HIV replication through cell-mediated immune responses. No consistent associations have been found between the breadth of Cytotoxic T Lymphocytes (CTL) responses and the control of HIV-1, and it is unknown whether the size or distribution of the viral proteome-wide epitope repertoire, i.e., the intrinsic ability to present fewer, more or specific viral epitopes, could affect clinical markers of disease progression. METHODOLOGY/PRINCIPAL FINDINGS: We used an epitope prediction model to identify all epitope motifs in a set of 302 HIV-1 full-length proteomes according to each individual's HLA (Human Leukocyte Antigen) genotype. The epitope repertoire, i.e., the number of predicted epitopes per HIV-1 proteome, varied considerably between HLA alleles and thus among individual proteomes. In a subgroup of 270 chronically infected individuals, we found that lower viral loads and higher CD4 counts were associated with a larger predicted epitope repertoire. Additionally, in Gag and Rev only, more epitopes were restricted by alleles associated with low viral loads than by alleles associated with higher viral loads. CONCLUSIONS/SIGNIFICANCE: This comprehensive analysis puts forth the epitope repertoire as a mechanistic component of the multi-faceted HIV-specific CTL response. The favorable impact on markers of disease status of the propensity to present more HLA binding peptides and specific proteins gives impetus to vaccine design strategies that seek to elicit responses to a broad array of HIV-1 epitopes, and suggest a particular focus on Gag.

Extensive intrasubtype recombination in South African human immunodeficiency virus type 1 subtype C infections.

Recombinant human immunodeficiency virus type 1 (HIV-1) strains containing sequences from different viral genetic subtypes (intersubtype) and different lineages from within the same subtype (intrasubtype) have been observed. A consequence of recombination can be the distortion of the phylogenetic signal. Several intersubtype recombinants have been identified; however, less is known about the frequency of intrasubtype recombination. For this study, near-full-length HIV-1 subtype C genomes from 270 individuals were evaluated for the presence of intrasubtype recombination. A sliding window schema (window, 2 kb; step, 385 bp) was used to partition the aligned sequences. The Shimodaira-Hasegawa test detected significant topological incongruence in 99.6% of the comparisons of the maximum-likelihood trees generated from each sequence partition, a result that could be explained by recombination. Using RECOMBINE, we detected significant levels of recombination using five random subsets of the sequences. With a set of 23 topologically consistent sequences used as references, bootscanning followed by the interactive informative site test defined recombination breakpoints. Using two multiple-comparison correction methods, 47% of the sequences showed significant evidence of recombination in both analyses. Estimated evolutionary rates were revised from 0.51%/year (95% confidence interval [CI], 0.39 to 0.53%) with all sequences to 0.46%/year (95% CI, 0.38 to 0.48%) with the putative recombinants removed. The timing of the subtype C epidemic origin was revised from 1961 (95% CI, 1947 to 1962) with all sequences to 1958 (95% CI, 1949 to 1960) with the putative recombinants removed. Thus, intrasubtype recombinants are common within the subtype C epidemic and these impact analyses of HIV-1 evolution.

Control of human immunodeficiency virus type 1 is associated with HLA-B*13 and targeting of multiple gag-specific CD8+ T-cell epitopes.

To better understand relationships between CD8+ T-cell specificity and the immune control of human immunodeficiency virus type 1 (HIV-1), we analyzed the role of HLA-B*13, an allele associated with low viremia, in a cohort of 578 C clade-infected individuals in Durban, South Africa. Six novel B*13-restricted cytotoxic T lymphocyte epitopes were defined from analyses of 37 B*13-positive subjects, including three Gag epitopes. These B*13-restricted epitopes contribute to a broad Gag-specific CD8+ response that is associated with the control of viremia. These data are consistent with data from studies of other HLA-class I alleles associated with HIV control that have shown that the targeting of multiple Gag epitopes is associated with relative suppression of viremia.

Differential selection pressure exerted on HIV by CTL targeting identical epitopes but restricted by distinct HLA alleles from the same HLA supertype.

HLA diversity is seen as a major challenge to CTL vaccines against HIV. One current approach focuses on "promiscuous" epitopes, presented by multiple HLA alleles from within the same HLA supertype. However, the effectiveness of such supertype vaccines depends upon the functional equivalence of CTL targeting a particular epitope, irrespective of the restricting HLA. In this study, we describe the promiscuous HIV-specific CTL epitopes presented by alleles within the B7 supertype. Substantial differences were observed in the ability of CTL to select for escape mutation when targeting the same epitope but restricted by different HLA. This observation was common to all six promiscuous B7 epitopes identified. Moreover, with one exception, there were no significant differences in the frequency, magnitude, or immunodominance of the CTL responses restricted by different HLA alleles to explain these discrepancies. This suggests that the unique peptide/MHC complexes generated by even closely related HLA induce CTL responses that are qualitatively different. This hypothesis is supported by additional differences observed between CTL targeting identical epitopes but restricted by different HLA: first, the occurrence of distinct, HLA-specific escape mutation; second, the recruitment of distinct TCR repertoires by particular peptide/MHC complexes; and, third, significant differences in the functional avidity of CTL. Taken together, these data indicate that significant functional differences exist between CTL targeting identical epitopes but restricted by different, albeit closely related HLA. These findings are of relevance to vaccine approaches that seek to exploit HLA supertypes to overcome the problem of HLA diversity.

Subtype C Is associated with increased vaginal shedding of HIV-1.

The prevalence of human immunodeficiency virus (HIV)-1-infected cells and HIV-1 RNA levels in genital secretions and breast milk and the risk of mother-to-child transmission of HIV-1 were compared among subtypes A, C, and D in a Kenyan cohort. Pregnant women infected with subtype C were significantly more likely to shed HIV-1-infected vaginal cells than were those infected with subtype A or D (odds ratio [OR], 3.6 [95% confidence interval {CI}, 1.4-8.8]; P = .006). This relationship held after adjusting for age, CD4 cell count, and plasma HIV-1 RNA load (OR, 3.1 [95% CI, 1.1-8.6]; P = .03). These observations suggest that HIV-1 subtype influences mucosal shedding of HIV-1.

Association of levels of HIV-1-infected breast milk cells and risk of mother-to-child transmission.

Understanding how the level of human immunodeficiency virus type 1 (HIV-1)-infected breast milk cells (BMCs) affects HIV transmission via breast-feeding can shed light on the mechanism of infection and aid in establishing effective interventions. The proportion of infected cells to total cells was measured in serial breast milk samples collected from 291 HIV-1-infected women in Nairobi, Kenya, by use of real-time DNA polymerase chain reaction amplification of BMCs. The number of infected BMCs per million cells was associated with levels of cell-free viral RNA in breast milk (R=.144; P=.032), levels of cell-free virus in blood plasma (R=.365; P<.001), and the detection of proviral DNA in cervical and vaginal secretions (P<.001 and P = .030, respectively). The number of infected BMCs per million cells was lower in colostrum or early milk than in mature milk (P<.001). Previous studies demonstrated that the concentration of BMCs varies throughout lactation, and we used these data to transform infected BMCs per million cells to infected BMCs per milliliter. The estimated concentration of infected BMCs per milliliter was higher in colostrum or early milk than in mature milk (P<.001). Each log10 increase in infected BMCs per milliliter was associated with a 3.19-fold-increased risk of transmission (P=.002), after adjustment for cell-free virus in plasma (hazard ratio [HR], 2.09; P=.03) and breast milk (HR, 1.01; P=1.00). This suggests that infected BMCs may play a more important role in transmission of HIV via breast-feeding than does cell-free virus.

Longitudinal analysis of human immunodeficiency virus type 1 RNA in breast milk and of its relationship to infant infection and maternal disease.

Transmission of human immunodeficiency virus type 1 (HIV-1) via breast-feeding can occur throughout lactation. Defining both fluctuation in breast-milk virus level over time and how breast-milk virus correlates with mother-to-child transmission is important for establishing effective interventions. We quantified breast-milk HIV-1 RNA levels in serial samples collected from 275 women for up to 2 years after delivery. Higher maternal plasma virus load, lower maternal CD4 T cell count, and detection of HIV-1 DNA in maternal genital secretions were significantly associated with elevated breast-milk HIV-1 RNA. Within women who breast-fed, median virus load in colostrum/early milk was significantly higher than that in mature breast milk collected 14 days after delivery (P< or =.004). Breast-feeding mothers who transmitted HIV-1 to their infants had both significantly higher breast-milk viral RNA throughout lactation and more-consistent viral shedding, compared with mothers who did not transmit HIV-1. In breast-feeding women, a 2-fold-increased risk of transmission was associated with every 10-fold increase in breast-milk virus load (95% confidence interval, 1.3-3.0; P<.001). These results indicate that the risk of infant infection from breast-feeding is influenced by breast-milk virus load, which is highest early after delivery.