Deep Sequencing of Protease Inhibitor Resistant HIV Patient Isolates Reveals Patterns of Correlated Mutations in Gag and Protease

While the role of drug resistance mutations in HIV protease has been studied comprehensively, mutations in its substrate, Gag, have not been extensively cataloged. Using deep sequencing, we analyzed a unique collection of longitudinal viral samples from 93 patients who have been treated with therapies containing protease inhibitors (PIs). Due to the high sequence coverage within each sample, the frequencies of mutations at individual positions were calculated with high precision. We used this information to characterize the variability in the Gag polyprotein and its effects on PI-therapy outcomes. To examine covariation of mutations between two different sites using deep sequencing data, we developed an approach to estimate the tight bounds on the two-site bivariate probabilities in each viral sample, and the mutual information between pairs of positions based on all the bounds. Utilizing the new methodology we found that mutations in the matrix and p6 proteins contribute to continued therapy failure and have a major role in the network of strongly correlated mutations in the Gag polyprotein, as well as between Gag and protease. Although covariation is not direct evidence of structural propensities, we found the strongest correlations between residues on capsid and matrix of the same Gag protein were often due to structural proximity. This suggests that some of the strongest inter-protein Gag correlations are the result of structural proximity. Moreover, the strong covariation between residues in matrix and capsid at the N-terminus with p1 and p6 at the C-terminus is consistent with residue-residue contacts between these proteins at some point in the viral life cycle

Lost but not forgotten: MHC genotypes predict overwinter survival despite depauperate MHC diversity in a declining frog

The amphibian disease chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), has contributed to the decline of Chiricahua leopard frogs (Rana chiricahuensis), a federally threatened species native to the Southwestern United States. We characterized immunogenetic variability in R. chiricahuensis by sequencing an expressed Major Histocompatibility Complex (MHC) class IIβ gene across 13 natural populations in Arizona, USA, as well as 283 individuals that were captive reared from two egg masses. We recovered a total of five class IIβ MHC alleles compared to 84 alleles previously characterized in eight natural populations of the Arizona congener R. yavapaiensis, demonstrating reduced MHC diversity in R. chiricahuensis. One allele was fixed in five populations but none of the R. chiricahuensis alleles were closely related to R. yavapaiensis allele Q, which is significantly associated with chytridiomycosis resistance in laboratory trials. Nine of 13 R. chiricahuensis population localities were Bd positive, and bearing allele RachDRB*04 was the best genetic predictor of an individual being infected with Bd. A total of three class IIβ alleles were recovered from captive reared individuals, which were released to two natural population localities followed by recapture surveys to assess MHC-based survival over winter, the time when chytridiomycosis outbreaks are most severe. At one site, all released animals were fixed for a single allele and MHC-based survival could not be assessed. At the second site, fewer than half of the released but all of the recaptured individuals were homozygous for RachDRB*05, indicating that MHC genotype is important in determining Bd survival under natural field conditions. We conclude that the limited MHC variation in R. chiricahuensis is likely the consequence rather than the cause of natural selection favoring alleles that promote survival in the face of Bd. Our study highlights that preserving even low levels of functional genetic variation may be essential for population persistence, and that local disease adaptation may present as a reduction in genetic diversity. These finding also suggest that for populations that have declined due to a specific infectious pathogen, MHC-based genetically-informed reintroduction approaches may enhance species recovery efforts