Frequent toggling between alternative amino acids is driven by selection in HIV-1

Dynamic expedition of leading mutations in SARS-CoV-2 spike glycoproteins

The continuous evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the recent pandemic, has generated countless new variants with varying fitness. Mutations of the spike glycoprotein play a particularly vital role in shaping its evolutionary trajectory, as they have the capability to alter its infectivity and antigenicity. We present a time-resolved statistical method, Dynamic Expedition of Leading Mutations (deLemus), to analyze the evolutionary dynamics of the SARS-CoV-2 spike glycoprotein. The proposed L -index of the deLemus method is effective in quantifying the mutation strength of each amino acid site and outlining evolutionarily significant sites, allowing the comprehensive characterization of the evolutionary mutation pattern of the spike glycoprotein.

Dissecting positive selection events and immunological drives during the evolution of adeno-associated virus lineages

Adeno-associated virus (AAV) serotypes from primates are being developed and clinically used as vectors for human gene therapy. However, the evolutionary mechanism of AAV variants is far from being understood, except that genetic recombination plays an important role. Furthermore, little is known about the interaction between AAV and its natural hosts, human and nonhuman primates. In this study, natural AAV capsid genes were subjected to systemic evolutionary analysis with a focus on selection drives during the diversification of AAV lineages. A number of positively selected sites were identified from these AAV lineages with functional relevance implied by their localization on the AAV structures. The selection drives of the two AAV2 capsid sites were further investigated in a series of biological experiments. These observations did not support the evolution of the site 410 of the AAV2 capsid driven by selection pressure from the human CD4+ T-cell response. However, positive selection on site 548 of the AAV2 capsid was directly related to host humoral immunity because of the profound effects of mutations at this site on the immune evasion of AAV variants from human neutralizing antibodies at both the individual and population levels. Overall, this work provides a novel interpretation of the genetic diversity and evolution of AAV lineages in their natural hosts, which may contribute to their further engineering and application in human gene therapy.

Identification of evolutionary trajectories shared across human betacoronaviruses

Comparing the evolution of distantly related viruses can provide insights into common adaptive processes related to shared ecological niches. Phylogenetic approaches, coupled with other molecular evolution tools, can help identify mutations informative on adaptation, whilst the structural contextualization of these to functional sites of proteins may help gain insight into their biological properties. Two zoonotic betacoronaviruses capable of sustained human-to-human transmission have caused pandemics in recent times (SARS-CoV-1 and SARS-CoV-2), whilst a third virus (MERS-CoV) is responsible for sporadic outbreaks linked to animal infections. Moreover, two other betacoronaviruses have circulated endemically in humans for decades (HKU1 and OC43). To search for evidence of adaptive convergence between established and emerging betacoronaviruses capable of sustained human-to-human transmission (HKU1, OC43, SARS-CoV-1 and SARS-CoV-2), we developed a methodological pipeline to classify shared non-synonymous mutations as putatively denoting homoplasy (repeated mutations that do not share direct common ancestry) or stepwise evolution (sequential mutations leading towards a novel genotype). In parallel, we look for evidence of positive selection, and draw upon protein structure data to identify potential biological implications. We find 30 mutations, with four of these [codon sites 18121 (nsp14/residue 28), 21623 (spike/21), 21635 (spike/25) and 23948 (spike/796); SARS-CoV-2 genome numbering] displaying evolution under positive selection and proximity to functional protein regions. Our findings shed light on potential mechanisms underlying betacoronavirus adaptation to the human host and pinpoint common mutational pathways that may occur during establishment of human endemicity.

Nanopore Amplicon Sequencing Reveals Molecular Convergence and Local Adaptation of Rhodopsin in Great Lakes Salmonids

Local adaptation can drive diversification of closely related species across environmental gradients and promote convergence of distantly related taxa that experience similar conditions. We examined a potential case of adaptation to novel visual environments in a species flock (Great Lakes salmonids, genus Coregonus) using a new amplicon genotyping protocol on the Oxford Nanopore Flongle and MinION. We sequenced five visual opsin genes for individuals of Coregonus artedi, Coregonus hoyi, Coregonus kiyi, and Coregonus zenithicus. Comparisons revealed species-specific differences in a key spectral tuning amino acid in rhodopsin (Tyr261Phe substitution), suggesting local adaptation of C. kiyi to the blue-shifted depths of Lake Superior. Ancestral state reconstruction demonstrates that parallel evolution and "toggling" at this amino acid residue has occurred several times across the fish tree of life, resulting in identical changes to the visual systems of distantly related taxa across replicated environmental gradients. Our results suggest that ecological differences and local adaptation to distinct visual environments are strong drivers of both evolutionary parallelism and diversification.

Evolution of respiratory syncytial virus genotype BA in Kilifi, Kenya, 15 years on

Respiratory syncytial virus (RSV) is recognised as a leading cause of severe acute respiratory disease and deaths among infants and vulnerable adults. Clinical RSV isolates can be divided into several known genotypes. RSV genotype BA, characterised by a 60-nucleotide duplication in the G glycoprotein gene, emerged in 1999 and quickly disseminated globally replacing other RSV group B genotypes. Continual molecular epidemiology is critical to understand the evolutionary processes maintaining the success of the BA viruses. We analysed 735 G gene sequences from samples collected from paediatric patients in Kilifi, Kenya, between 2003 and 2017. The virus population comprised of several genetically distinct variants (n = 56) co-circulating within and between epidemics. In addition, there was consistent seasonal fluctuations in relative genetic diversity. Amino acid changes increasingly accumulated over the surveillance period including two residues (N178S and Q180R) that mapped to monoclonal antibody 2D10 epitopes, as well as addition of putative N-glycosylation sequons. Further, switching and toggling of amino acids within and between epidemics was observed. On a global phylogeny, the BA viruses from different countries form geographically isolated clusters suggesting substantial localized variants. This study offers insights into longitudinal population dynamics of a globally endemic RSV genotype within a discrete location.

Deep Profiling Identifies Selection of Nonsynonymous Amino Acid Substitutions in HIV-1 Envelope During Early Infection

Understanding the evolutionary dynamics of the viruses within an individual at or near the moment of transmission can provide critical inputs for the design of an effective vaccine for HIV infection. In this study, high-throughput sequencing technology was employed to analyze the evolutionary rate in viruses obtained at a single time point from drug-naive recently infected infants and adults in the chronic stage of disease. Gene-wise nonsynonymous (pN) and synonymous (pS) mutation rates were estimated and compared between the two groups. Significant differences were observed in the evolutionary rates between viruses in the early and late stages of infection. Higher rates of adaptive mutations in the HIV-1 envelope gene (env) were found in the chronic viruses as compared with those in the early stages of HIV infection. Conversely, percentage of nonsynonymous substitutions in env was found to be higher in recently transmitted viruses. In addition, a positive correlation was found between mutation and the evolutionary rate, and infectivity titer in recent infection. Despite the small sample size, the study identified useful information about viral evolution on transmission-associated bottlenecks. The effect of intraindividual HIV-1 evolution at the population level was highly contemporary, and the higher percentage of nonsynonymous substitutions seen in env during recent HIV-1 infection has suggested a pattern of convergent evolution leading to a positive selection for survival fitness and disease progression.

A de novo approach to inferring within-host fitness effects during untreated HIV-1 infection

In the absence of effective antiviral therapy, HIV-1 evolves in response to the within-host environment, of which the immune system is an important aspect. During the earliest stages of infection, this process of evolution is very rapid, driven by a small number of CTL escape mutations. As the infection progresses, immune escape variants evolve under reduced magnitudes of selection, while competition between an increasing number of polymorphic alleles (i.e., clonal interference) makes it difficult to quantify the magnitude of selection acting upon specific variant alleles. To tackle this complex problem, we developed a novel multi-locus inference method to evaluate the role of selection during the chronic stage of within-host infection. We applied this method to targeted sequence data from the p24 and gp41 regions of HIV-1 collected from 34 patients with long-term untreated HIV-1 infection. We identify a broad distribution of beneficial fitness effects during infection, with a small number of variants evolving under strong selection and very many variants evolving under weaker selection. The uniquely large number of infections analysed granted a previously unparalleled statistical power to identify loci at which selection could be inferred to act with statistical confidence. Our model makes no prior assumptions about the nature of alleles under selection, such that any synonymous or non-synonymous variant may be inferred to evolve under selection. However, the majority of variants inferred with confidence to be under selection were non-synonymous in nature, and in most cases were have previously been associated with either CTL escape in p24 or neutralising antibody escape in gp41. We also identified a putative new CTL escape site (residue 286 in gag), and a region of gp41 (including residues 644, 648, 655 in env) likely to be associated with immune escape. Sites inferred to be under selection in multiple hosts have high within-host and between-host diversity although not all sites with high between-host diversity were inferred to be under selection at the within-host level. Our identification of selection at sites associated with resistance to broadly neutralising antibodies (bNAbs) highlights the need to fully understand the role of selection in untreated individuals when designing bNAb based therapies.

Bat adeno-associated viruses as gene therapy vectors with the potential to evade human neutralizing antibodies

The prevalence of adeno-associated virus (AAV) has been investigated in bat populations, but little is known about the biological properties of this virus. In this study, four full-length bat AAV capsid genes were isolated in China, with their amino acid sequences sharing 61% identity with those of AAV2 on average. These capsid genes could package AAV particles in combination with AAV2 rep and ITRs, albeit at a lower efficiency. Bat AAVs could only slightly infect mouse liver but could transduce mouse muscle to some extent after systemic administration with a higher muscle/liver ratio than that of primate AAVs. Bat AAV 10HB showed moderate muscle transduction, similar to that of AAV2, during direct intramuscular injection and, compared with other AAV serotypes, was also relatively efficient in resisting human antibody neutralization after intramuscular injection. Evolutionary analysis revealed a number of codons in bat AAV capsid genes subject to positive selection, with sites corresponding to V259 and N691 in 10HB capsids being localized on the surface of the AAV2 capsid. Mutagenesis studies indicated that the positive selection in bat AAV capsids is driven by their tropism evolution in host species. Taken together, the results of this study indicate that bat AAV 10HB vector has the possible applications for muscular gene therapy, especially in the presence of human AAV neutralizing antibodies.

Genotypic and Mechanistic Characterization of Subtype-Specific HIV Adaptation to Host Cellular Immunity

The extent to which viral genetic context influences HIV adaptation to human leukocyte antigen (HLA) class I-restricted immune pressures remains incompletely understood. The Ugandan HIV epidemic, where major pandemic group M subtypes A1 and D cocirculate in a single host population, provides an opportunity to investigate this question. We characterized plasma HIV RNA gag, pol, and nef sequences, along with host HLA genotypes, in 464 antiretroviral-naive individuals chronically infected with HIV subtype A1 or D. Using phylogenetically informed statistical approaches, we identified HLA-associated polymorphisms and formally compared their strengths of selection between viral subtypes. A substantial number (32%) of HLA-associated polymorphisms identified in subtype A1 and/or D had previously been reported in subtype B, C, and/or circulating recombinant form 01_AE (CRF01_AE), confirming the shared nature of many HLA-driven escape pathways regardless of viral genetic context. Nevertheless, 34% of the identified HLA-associated polymorphisms were significantly differentially selected between subtypes A1 and D. Experimental investigation of select examples of subtype-specific escape revealed distinct underlying mechanisms with important implications for vaccine design: whereas some were attributable to subtype-specific sequence variation that influenced epitope-HLA binding, others were attributable to differential mutational barriers to immune escape. Overall, our results confirm that HIV genetic context is a key modulator of viral adaptation to host cellular immunity and highlight the power of combined bioinformatic and mechanistic studies, paired with knowledge of epitope immunogenicity, to identify appropriate viral regions for inclusion in subtype-specific and universal HIV vaccine strategies.IMPORTANCE The identification of HIV polymorphisms reproducibly selected under pressure by specific HLA alleles and the elucidation of their impact on viral function can help identify immunogenic viral regions where immune escape incurs a fitness cost. However, our knowledge of HLA-driven escape pathways and their functional costs is largely limited to HIV subtype B and, to a lesser extent, subtype C. Our study represents the first characterization of HLA-driven adaptation pathways in HIV subtypes A1 and D, which dominate in East Africa, and the first statistically rigorous characterization of differential HLA-driven escape across viral subtypes. The results support a considerable impact of viral genetic context on HIV adaptation to host HLA, where HIV subtype-specific sequence variation influences both epitope-HLA binding and the fitness costs of escape. Integrated bioinformatic and mechanistic characterization of these and other instances of differential escape could aid rational cytotoxic T-lymphocyte-based vaccine immunogen selection for both subtype-specific and universal HIV vaccines.

Relative Evolutionary Rates in Proteins Are Largely Insensitive to the Substitution Model

The relative evolutionary rates at individual sites in proteins are informative measures of conservation or adaptation. Often used as evolutionarily aware conservation scores, relative rates reveal key functional or strongly selected residues. Estimating rates in a phylogenetic context requires specifying a protein substitution model, which is typically a phenomenological model trained on a large empirical data set. A strong emphasis has traditionally been placed on selecting the "best-fit" model, with the implicit understanding that suboptimal or otherwise ill-fitting models might bias inferences. However, the pervasiveness and degree of such bias has not been systematically examined. We investigated how model choice impacts site-wise relative rates in a large set of empirical protein alignments. We compared models designed for use on any general protein, models designed for specific domains of life, and the simple equal-rates Jukes Cantor-style model (JC). As expected, information theoretic measures showed overwhelming evidence that some models fit the data decidedly better than others. By contrast, estimates of site-specific evolutionary rates were impressively insensitive to the substitution model used, revealing an unexpected degree of robustness to potential model misspecification. A deeper examination of the fewer than 5% of sites for which model inferences differed in a meaningful way showed that the JC model could uniquely identify rapidly evolving sites that models with empirically derived exchangeabilities failed to detect. We conclude that relative protein rates appear robust to the applied substitution model, and any sensible model of protein evolution, regardless of its fit to the data, should produce broadly consistent evolutionary rates.

Evolution of HIV-1 within untreated individuals and at the population scale in Uganda

HIV-1 undergoes multiple rounds of error-prone replication between transmission events, resulting in diverse viral populations within and among individuals. In addition, the virus experiences different selective pressures at multiple levels: during the course of infection, at transmission, and among individuals. Disentangling how these evolutionary forces shape the evolution of the virus at the population scale is important for understanding pathogenesis, how drug- and immune-escape variants are likely to spread in populations, and the development of preventive vaccines. To address this, we deep-sequenced two regions of the HIV-1 genome (p24 and gp41) from 34 longitudinally-sampled untreated individuals from Rakai District in Uganda, infected with subtypes A, D, and inter-subtype recombinants. This dataset substantially increases the availability of HIV-1 sequence data that spans multiple years of untreated infection, in particular for different geographical regions and viral subtypes. In line with previous studies, we estimated an approximately five-fold faster rate of evolution at the within-host compared to the population scale for both synonymous and nonsynonymous substitutions, and for all subtypes. We determined the extent to which this mismatch in evolutionary rates can be explained by the evolution of the virus towards population-level consensus, or the transmission of viruses similar to those that establish infection within individuals. Our findings indicate that both processes are likely to be important.

Phylodynamics of classical swine fever virus with emphasis on Ecuadorian strains

Classic swine fever virus (CSFV) is a Pestivirus from the Flaviviridae family that affects pigs worldwide and is endemic in several Latin American countries. However, there are still some countries in the region, including Ecuador, for which CSFV molecular information is lacking. To better understand the epidemiology of CSFV in the Americas, sequences from CSFVs from Ecuador were generated and a phylodynamic analysis of the virus was performed. Sequences for the full-length glycoprotein E2 gene of twenty field isolates were obtained and, along with sequences from strains previously described in the Americas and from the most representative strains worldwide, were used to analyse the phylodynamics of the virus. Bayesian methods were used to test several molecular clock and demographic models. A calibrated ultrametric tree and a Bayesian skyline were constructed, and codons associated with positive selection involving immune scape were detected. The best model according to Bayes factors was the strict molecular clock and Bayesian skyline model, which shows that CSFV has an evolution rate of 3.2 × 10^-4 substitutions per site per year. The model estimates the origin of CSFV in the mid-1500s. There is a strong spatial structure for CSFV in the Americas, indicating that the virus is moving mainly through neighbouring countries. The genetic diversity of CSFV has increased constantly since its appearance, with a slight decrease in mid-twentieth century, which coincides, with eradication campaigns in North America. Even though there is no evidence of strong directional evolution of the E2 gene in CSFV, codons 713, 761, 762 and 975 appear to be selected positively and could be related to virulence or pathogenesis. These results reveal how CSFV has spread and evolved since it first appeared in the Americas and provide important information for attaining the goal of eradication of this virus in Latin America.

Small-scale intraspecific patterns of adaptive immunogenetic polymorphisms and neutral variation in Lake Superior lake trout

Many fishes express high levels of intraspecific variability, often linked to resource partitioning. Several studies show that a species' evolutionary trajectory of adaptive divergence can undergo reversals caused by changes in its environment. Such a reversal in neutral genetic and morphological variation among lake trout Salvelinus namaycush ecomorphs appears to be underway in Lake Superior. However, a water depth gradient in neutral genetic divergence was found to be associated with intraspecific diversity in the lake. To investigate patterns of adaptive immunogenetic variation among lake trout ecomorphs, we used Illumina high-throughput sequencing. The population's genetic structure of the major histocompatibility complex (MHC Class IIβ exon 2) and 18 microsatellite loci were compared to disentangle neutral and selective processes at a small geographic scale. Both MHC and microsatellite variation were partitioned more by water depth stratum than by ecomorph. Several metrics showed strong clustering by water depth in MHC alleles, but not microsatellites. We report a 75% increase in the number of MHC alleles shared between the predominant shallow and deep water ecomorphs since a previous lake trout MHC study at the same locale (c. 1990s data). This result is consistent with the reverse speciation hypothesis, although adaptive MHC polymorphisms persist along an ecological gradient. Finally, results suggested that the lake trout have multiple copies of the MHC II locus consistent with a historic genomic duplication event. Our findings indicated that conservation approaches for this species could focus on managing various ecological habitats by depth, in addition to regulating the fisheries specific to ecomorphs.

A Comparison of One-Rate and Two-Rate Inference Frameworks for Site-Specific dN/dS Estimation

Two broad paradigms exist for inferring [Formula: see text] the ratio of nonsynonymous to synonymous substitution rates, from coding sequences: (i) a one-rate approach, where [Formula: see text] is represented with a single parameter, or (ii) a two-rate approach, where [Formula: see text] and [Formula: see text] are estimated separately. The performances of these two approaches have been well studied in the specific context of proper model specification, i.e., when the inference model matches the simulation model. By contrast, the relative performances of one-rate vs. two-rate parameterizations when applied to data generated according to a different mechanism remain unclear. Here, we compare the relative merits of one-rate and two-rate approaches in the specific context of model misspecification by simulating alignments with mutation-selection models rather than with [Formula: see text]-based models. We find that one-rate frameworks generally infer more accurate [Formula: see text] point estimates, even when [Formula: see text] varies among sites. In other words, modeling [Formula: see text] variation may substantially reduce accuracy of [Formula: see text] point estimates. These results appear to depend on the selective constraint operating at a given site. For sites under strong purifying selection ([Formula: see text]), one-rate and two-rate models show comparable performances. However, one-rate models significantly outperform two-rate models for sites under moderate-to-weak purifying selection. We attribute this distinction to the fact that, for these more quickly evolving sites, a given substitution is more likely to be nonsynonymous than synonymous. The data will therefore be relatively enriched for nonsynonymous changes, and modeling [Formula: see text] contributes excessive noise to [Formula: see text] estimates. We additionally find that high levels of divergence among sequences, rather than the number of sequences in the alignment, are more critical for obtaining precise point estimates.

Functional evolution of the OAS1 viral sensor: Insights from old world primates

Infections with viral pathogens impose considerable selective pressure on host defensive genes. Those genes at the forefront, responsible for identifying and binding exogenous molecular viral components, will carry the hallmarks of this struggle. Oligoadenylate synthetase (OAS) enzymes play a major role in the innate defense against a large number of viruses by acting as sensors of viral infections. Following their up-regulation by the interferon pathway, OASs bind viral dsRNA and then signal ribonuclease L (RNase L) to degrade RNA, shutting down viral and host protein synthesis. We have investigated the evolution of OAS1 in twenty-two Old World monkey species. We identified a total of 35 codons with the earmarks of positive selection and we performed a comprehensive analysis of their functional significance using in silico modeling of the OAS1 protein. Subdividing OAS1 into functional domains revealed intense purifying selection in the active domain but significant positive directional selection in the RNA-binding domain (RBD), the region where OAS1 binds viral dsRNA. The modeling analysis revealed a concentration of rapidly evolving residues in one region of the RBD suggestive of the sub-functionalization of different regions of the RBD. This analysis also identified several positively selected residues circumscribing the entry to the active site suggesting adaptive evasion of viral antagonism and/or selection for production of oligoadenylate of different length.

Modelling the role of immunity in reversion of viral antigenic sites

Antigenic sites in viral pathogens exhibit distinctive evolutionary dynamics due to their role in evading recognition by host immunity. Antigenic selection is known to drive higher rates of non-synonymous substitution; less well understood is why differences are observed between viruses in their propensity to mutate to a novel or previously encountered amino acid. Here, we present a model to explain patterns of antigenic reversion and forward substitution in terms of the epidemiological and molecular processes of the viral population. We develop an analytical three-strain model and extend the analysis to a multi-site model to predict characteristics of observed sequence samples. Our model provides insight into how the balance between selection to escape immunity and to maintain viability is affected by the rate of mutational input. We also show that while low probabilities of reversion may be due to either a low cost of immune escape or slowly decaying host immunity, these two scenarios can be differentiated by the frequency patterns at antigenic sites. Comparison between frequency patterns of human influenza A (H3N2) and human RSV-A suggests that the increased rates of antigenic reversion in RSV-A is due to faster decaying immunity and not higher costs of escape.

Interspecific adaptation by binary choice at de novo polyomavirus T antigen site through accelerated codon-constrained Val-Ala toggling within an intrinsically disordered region

It is common knowledge that conserved residues evolve slowly. We challenge generality of this central tenet of molecular biology by describing the fast evolution of a conserved nucleotide position that is located in the overlap of two open reading frames (ORFs) of polyomaviruses. The de novo ORF is expressed through either the ALTO protein or the Middle T antigen (MT/ALTO), while the ancestral ORF encodes the N-terminal domain of helicase-containing Large T (LT) antigen. In the latter domain the conserved Cys codon of the LXCXE pRB-binding motif constrains codon evolution in the overlapping MT/ALTO ORF to a binary choice between Val and Ala codons, termed here as codon-constrained Val-Ala (COCO-VA) toggling. We found the rate of COCO-VA toggling to approach the speciation rate and to be significantly accelerated compared to the baseline rate of chance substitution in a large monophyletic lineage including all viruses encoding MT/ALTO and three others. Importantly, the COCO-VA site is located in a short linear motif (SLiM) of an intrinsically disordered region, a typical characteristic of adaptive responders. These findings provide evidence that the COCO-VA toggling is under positive selection in many polyomaviruses, implying its critical role in interspecific adaptation, which is unprecedented for conserved residues.

A genotype network reveals homoplastic cycles of convergent evolution in influenza A (H3N2) haemagglutinin

Networks of evolving genotypes can be constructed from the worldwide time-resolved genotyping of pathogens like influenza viruses. Such genotype networks are graphs where neighbouring vertices (viral strains) differ in a single nucleotide or amino acid. A rich trove of network analysis methods can help understand the evolutionary dynamics reflected in the structure of these networks. Here, I analyse a genotype network comprising hundreds of influenza A (H3N2) haemagglutinin genes. The network is rife with cycles that reflect non-random parallel or convergent (homoplastic) evolution. These cycles also show patterns of sequence change characteristic for strong and local evolutionary constraints, positive selection and mutation-limited evolution. Such cycles would not be visible on a phylogenetic tree, illustrating that genotype network analysis can complement phylogenetic analyses. The network also shows a distinct modular or community structure that reflects temporal more than spatial proximity of viral strains, where lowly connected bridge strains connect different modules. These and other organizational patterns illustrate that genotype networks can help us study evolution in action at an unprecedented level of resolution.

Viral CD8 T cell epitope nucleotide composition shows evidence of short- and long-term evolutionary strategies

Viral epitopes have a distinct codon usage that reflects their dual role in infection and immunity. On the one hand, epitopes are part of proteins important to viral function; on the other hand, they are targets of the immune response. Studies of selection are most commonly based on changes of amino acid and seen through the accumulation of non-synonymous mutations. An independent measure of selection is the codon usage and underlying changeability of the nucleotide sequences. We here use multiple tools and a large-scale analysis of viral genomes to demonstrate that viral epitopes have a distinct codon usage and that this codon usage reflects distinct short- and long-term types of selection during viral evolution. We show that CD8(+) T cell epitopes are encoded by codons more distant from stop codons and more changeable than codons outside epitopes. This biased codon usage reflects the viral population toggling back and forth from a wild-type sequence to an escape mode, which enable them to avoid immune detection when needed, and go back to the functionally favorable form when the threat is removed (i.e., in a new host).

Genotypic and functional impact of HIV-1 adaptation to its host population during the North American epidemic

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.

Genetic imprint of vaccination on simian/human immunodeficiency virus type 1 transmitted viral genomes in rhesus macaques

Understanding the genetic, antigenic and structural changes that occur during HIV-1 infection in response to pre-existing immunity will facilitate current efforts to develop an HIV-1 vaccine. Much is known about HIV-1 variation at the population level but little with regard to specific changes occurring in the envelope glycoprotein within a host in response to immune pressure elicited by antibodies. The aim of this study was to track and map specific early genetic changes occurring in the viral envelope gene following vaccination using a highly controlled viral challenge setting in the SHIV macaque model. We generated 449 full-length env sequences from vaccinees, and 63 from the virus inoculum. Analysis revealed a different pattern in the distribution and frequency of mutations in the regions of the envelope gene targeted by the vaccine as well as different patterns of diversification between animals in the naïve control group and vaccinees. Given the high stringency of the model it is remarkable that we were able to identify genetic changes associated with the vaccination. This work provides insight into the characterization of breakthrough viral populations in less than fully efficacious vaccines and illustrates the value of HIV-1 Env SHIV challenge model in macaques to unravel the mechanisms driving HIV-1 envelope genetic diversity in the presence of vaccine induced-responses.

Viral proteome size and CD8+ T cell epitope density are correlated: the effect of complexity on selection

•

We analyze the relation between viral complexity and their adaptation to the host immune system.

•

Viruses with few proteins and low number of nucleotides remove more CD8+ T cell epitopes.

•

Within a virus, short proteins (with fewer amino acids) adapt better than long ones.

•

The relation between total size and adaptation is host specific.

•

Complexity limits genetic adaptation in the high-mutation rate strong selection regime.

The relation between the complexity of organisms and proteins and their evolution rates has been discussed in the context of multiple generic models. The main robust claim from most such models is the negative relation between complexity and the accumulation rate of mutations.

Viruses accumulate escape mutations in their epitopes to avoid detection and destruction of their host cell by CD8+ T cells. The extreme regime of immune escape, namely, strong selection and high mutation rate, provide an opportunity to extend and validate the existing models of relation between complexity and evolution rate as proposed by Fisher and Kimura.

Using epitope prediction algorithms to compute the epitopes presented on the most frequent human HLA alleles in over 100 fully sequenced human viruses, and over 900 non-human viruses, we here study the correlation between viruses/proteins complexity (as measured by the number of proteins in the virus and the length of each protein, respectively) and the rate of accumulation of escape mutation. The latter is evaluated by measuring the normalized epitope density of viral proteins.

If the virus/protein complexity prevents the accumulation of escape mutations, the epitope density is expected to be positively correlated with both the number of proteins in the virus and the length of proteins. We show that such correlations are indeed observed for most human viruses. For non-human viruses the correlations were much less significant, indicating that the correlation is indeed induced by human HLA molecules.

Contrasted evolutionary constraints on secreted and non-secreted proteomes of selected Actinobacteria

BACKGROUND

Actinobacteria have adapted to contrasted ecological niches such as the soil, and among others to plants or animals as pathogens or symbionts. Mycobacterium genus contains mostly pathogens that cause a variety of mammalian diseases, among which the well-known leprosy and tuberculosis, it also has saprophytic relatives. Streptomyces genus is mostly a soil microbe known for its secondary metabolites, it contains also plant pathogens, animal pathogens and symbionts. Frankia, a nitrogen-fixing actinobacterium establishes a root symbiosis with dicotyledonous pionneer plants. Pathogens and symbionts live inside eukaryotic cells and tissues and interact with their cellular environment through secreted proteins and effectors transported through transmembrane systems; nevertheless they also need to avoid triggering host defense reactions. A comparative genome analysis of the secretomes of symbionts and pathogens allows a thorough investigation of selective pressures shaping their evolution. In the present study, the rates of silent mutations to non-silent mutations in secretory proteins were assessed in different strains of Frankia, Streptomyces and Mycobacterium, of which several genomes have recently become publicly available.

RESULTS

It was found that secreted proteins as a whole have a stronger purifying evolutionary rate (non-synonymous to synonymous substitutions or Ka/Ks ratio) than the non-secretory proteins in most of the studied genomes. This difference becomes statistically significant in cases involving obligate symbionts and pathogens. Amongst the Frankia, secretomes of symbiotic strains were found to have undergone evolutionary trends different from those of the mainly saprophytic strains. Even within the secretory proteins, the signal peptide part has a higher Ka/Ks ratio than the mature part. Two contrasting trends were noticed amongst the Frankia genomes regarding the relation between selection strength (i.e. Ka/Ks ratio) and the codon adaptation index (CAI), a predictor of the expression rate, in all the genes belonging to the core genome as well as the core secretory protein genes. The genomes of pathogenic Mycobacterium and Streptomyces also had reduced secretomes relative to saprophytes, as well as in general significant pairwise Ka/Ks ratios in their secretomes.

CONCLUSION

In marginally free-living facultative symbionts or pathogenic organisms under consideration, secretory protein genes as a whole evolve at a faster rate than the rest and this process may be an adaptive life-strategy to counter the host selection pressure. The higher evolutionary rate of signal peptide part compared to mature protein provides an indication that signal peptide parts may be under relaxed purifying selection, indicative of the signal peptides not being secreted into host cells. Codon usage analysis suggests that in actinobacterial strains under host selection pressure such as symbiotic Frankia, ACN, FD and the pathogenic Mycobacterium, codon usage bias was negatively correlated to the selective pressure exerted on the secretory protein genes.

Comparison of viral Env proteins from acute and chronic infections with subtype C human immunodeficiency virus type 1 identifies differences in glycosylation and CCR5 utilization and suggests a new strategy for immunogen design

Understanding human immunodeficiency virus type 1 (HIV-1) transmission is central to developing effective prevention strategies, including a vaccine. We compared phenotypic and genetic variation in HIV-1 env genes from subjects in acute/early infection and subjects with chronic infections in the context of subtype C heterosexual transmission. We found that the transmitted viruses all used CCR5 and required high levels of CD4 to infect target cells, suggesting selection for replication in T cells and not macrophages after transmission. In addition, the transmitted viruses were more likely to use a maraviroc-sensitive conformation of CCR5, perhaps identifying a feature of the target T cell. We confirmed an earlier observation that the transmitted viruses were, on average, modestly underglycosylated relative to the viruses from chronically infected subjects. This difference was most pronounced in comparing the viruses in acutely infected men to those in chronically infected women. These features of the transmitted virus point to selective pressures during the transmission event. We did not observe a consistent difference either in heterologous neutralization sensitivity or in sensitivity to soluble CD4 between the two groups, suggesting similar conformations between viruses from acute and chronic infection. However, the presence or absence of glycosylation sites had differential effects on neutralization sensitivity for different antibodies. We suggest that the occasional absence of glycosylation sites encoded in the conserved regions of env, further reduced in transmitted viruses, could expose specific surface structures on the protein as antibody targets.

On the statistical interpretation of site-specific variables in phylogeny-based substitution models

Phylogeny-based modeling of heterogeneity across the positions of multiple-sequence alignments has generally been approached from two main perspectives. The first treats site specificities as random variables drawn from a statistical law, and the likelihood function takes the form of an integral over this law. The second assigns distinct variables to each position, and, in a maximum-likelihood context, adjusts these variables, along with global parameters, to optimize a joint likelihood function. Here, it is emphasized that while the first approach directly enjoys the statistical guaranties of traditional likelihood theory, the latter does not, and should be approached with particular caution when the site-specific variables are high dimensional. Using a phylogeny-based mutation-selection framework, it is shown that the difference in interpretation of site-specific variables explains the incongruities in recent studies regarding distributions of selection coefficients.

Temporal association of HLA-B*81:01- and HLA-B*39:10-mediated HIV-1 p24 sequence evolution with disease progression

HLA-B*81:01 and HLA-B*39:10 alleles have been associated with viremic control in HIV-1 subtype C infection. Both alleles restrict the TL9 epitope in p24 Gag, and cytotoxic-T-lymphocyte (CTL)-mediated escape mutations in this epitope have been associated with an in vitro fitness cost to the virus. We investigated the timing and impact of mutations in the TL9 epitope on disease progression in five B*81:01- and two B*39:10-positive subtype C-infected individuals. Whereas both B*39:10 participants sampled at 2 months postinfection had viruses with mutations in the TL9 epitope, in three of the five (3/5) B*81:01 participants, TL9 escape mutations were only detected 10 months after infection, taking an additional 10 to 15 months to reach fixation. In the two remaining B*81:01 individuals, one carried a TL9 escape variant at 2 weeks postinfection, whereas no escape mutations were detected in the virus from the other participant for up to 33 months postinfection, despite CTL targeting of the epitope. In all participants, escape mutations in TL9 were linked to coevolving residues in the region of Gag known to be associated with host tropism. Late escape in TL9, together with coevolution of putative compensatory mutations, coincided with a spontaneous increase in viral loads in two individuals who were otherwise controlling the infection. These results provide in vivo evidence of the detrimental impact of B*81:01-mediated viral evolution, in a single Gag p24 epitope, on the control of viremia.

Detecting individual sites subject to episodic diversifying selection

The imprint of natural selection on protein coding genes is often difficult to identify because selection is frequently transient or episodic, i.e. it affects only a subset of lineages. Existing computational techniques, which are designed to identify sites subject to pervasive selection, may fail to recognize sites where selection is episodic: a large proportion of positively selected sites. We present a mixed effects model of evolution (MEME) that is capable of identifying instances of both episodic and pervasive positive selection at the level of an individual site. Using empirical and simulated data, we demonstrate the superior performance of MEME over older models under a broad range of scenarios. We find that episodic selection is widespread and conclude that the number of sites experiencing positive selection may have been vastly underestimated.

Identifying regions of protein coding genes that have undergone adaptive evolution is important to answering many questions in evolutionary biology and genetics. In order to tease out genetic evidence for natural selection, genes from a diverse array of taxa must be analyzed, only a subset of which may have undergone adaptive evolution; the same gene region may be under stabilizing or relaxed selection in lineages leading to other taxa. Most current computational methods designed to detect the imprint of natural selection at a site in a protein coding gene assume the strength and direction of natural selection is constant across all lineages. Here, we present a method to detect adaptive evolution, even when the selective forces are not constant across taxa. Using a variety of well-characterized genes, we find evidence suggesting that natural selection is generally episodic and that modeling it as such reveals that many more sites are subject to episodic positive selection than previously appreciated.

Viral evolution in deep time: lentiviruses and mammals

Lentiviruses are a distinctive genus of retroviruses that cause chronic, persistent infections in mammals, including humans. The emergence of pandemic HIV type-1 (HIV-1) infection during the late 20th century shaped a view of lentiviruses as 'modern' viruses. However, recent research has revealed an entirely different perspective, elucidating aspects of an evolutionary relationship with mammals that extends across many millions of years. Such deep evolutionary history is likely to be typical of many host-virus systems, fundamentally underpinning their interactions in the present day. For this reason, establishing the deep history of virus and host interaction is key to developing a fully informed approach to tackling viral diseases. Here, I use the example of lentiviruses to illustrate how paleovirological, geographic and genetic calibrations allow observations of virus and host interaction across a wide range of temporal and spatial scales to be integrated into a coherent ecological and evolutionary framework.

Is there a twelfth protein-coding gene in the genome of influenza A? A selection-based approach to the detection of overlapping genes in closely related sequences

Protein-coding genes often contain long overlapping open-reading frames (ORFs), which may or may not be functional. Current methods that utilize the signature of purifying selection to detect functional overlapping genes are limited to the analysis of sequences from divergent species, thus rendering them inapplicable to genes found only in closely related sequences. Here, we present a method for the detection of selection signatures on overlapping reading frames by using closely related sequences, and apply the method to several known overlapping genes, and to an overlapping ORF on the negative strand of segment 8 of influenza A virus (NEG8), for which the suggestion has been made that it is functional. We find no evidence that NEG8 is under selection, suggesting that the intact reading frame might be non-functional, although we cannot fully exclude the possibility that the method is not sensitive enough to detect the signature of selection acting on this gene. We present the limitations of the method using known overlapping genes and suggest several approaches to improve it in future studies. Finally, we examine alternative explanations for the sequence conservation of NEG8 in the absence of selection. We show that overlap type and genomic context affect the conservation of intact overlapping ORFs and should therefore be considered in any attempt of estimating the signature of selection in overlapping genes.

Analysis of a stochastic predator-prey model with applications to intrahost HIV genetic diversity

During an infection, HIV experiences strong selection by immune system T cells. Recent experimental work has shown that MHC escape mutations form an important pathway for HIV to avoid such selection. In this paper, we study a model of MHC escape mutation. The model is a predator-prey model with two prey, composed of two HIV variants, and one predator, the immune system CD8 cells. We assume that one HIV variant is visible to CD8 cells and one is not. The model takes the form of a system of stochastic differential equations. Motivated by well-known results concerning the short life-cycle of HIV intrahost, we assume that HIV population dynamics occur on a faster time scale then CD8 population dynamics. This separation of time scales allows us to analyze our model using an asymptotic approach. Using this model we study the impact of an MHC escape mutation on the population dynamics and genetic evolution of the intrahost HIV population. From the perspective of population dynamics, we show that the competition between the visible and invisible HIV variants can reach steady states in which either a single variant exists or in which coexistence occurs depending on the parameter regime. We show that in some parameter regimes the end state of the system is stochastic. From a genetics perspective, we study the impact of the population dynamics on the lineages of an HIV sample taken after an escape mutation occurs. We show that the lineages go through severe bottlenecks and that in certain parameter regimes the lineage distribution can be characterized by a Kingman coalescent. Our results depend on methods from diffusion theory and coalescent theory.

Purifying selection can obscure the ancient age of viral lineages

Statistical methods for molecular dating of viral origins have been used extensively to infer the time of most common recent ancestor for many rapidly evolving pathogens. However, there are a number of cases, in which epidemiological, historical, or genomic evidence suggests much older viral origins than those obtained via molecular dating. We demonstrate how pervasive purifying selection can mask the ancient origins of recently sampled pathogens, in part due to the inability of nucleotide-based substitution models to properly account for complex patterns of spatial and temporal variability in selective pressures. We use codon-based substitution models to infer the length of branches in viral phylogenies; these models produce estimates that are often considerably longer than those obtained with traditional nucleotide-based substitution models. Correcting the apparent underestimation of branch lengths suggests substantially older origins for measles, Ebola, and avian influenza viruses. This work helps to reconcile some of the inconsistencies between molecular dating and other types of evidence concerning the age of viral lineages.

Universal peptide vaccines - optimal peptide vaccine design based on viral sequence conservation

Rapidly mutating viruses such as the hepatitis C virus (HCV), the human immunodeficiency virus (HIV), or influenza viruses (Flu) call for highly effective universal peptide vaccines, i.e. vaccines that do not only yield broad population coverage but also broad coverage of various viral strains. The efficacy of such vaccines is determined by multiple properties of the epitopes they comprise. Beyond the specific properties of each epitope, properties of the corresponding source antigens are of great importance. If a response is mounted against viral proteins with a low copy number within the cell or against proteins expressed very late, this response may fail to induce lysis of the infected cells before budding can take place. We here propose a novel methodology to optimize the epitope composition and assembly in order to induce maximum protection. In order for a peptide vaccine to yield the best possible universal protection, several conditions should be met: (a) an optimal choice of target antigens, (b) an optimal choice of highly conserved epitopes, (c) maximum coverage of the target population, and (d) the proper ordering of the epitopes in the final vaccine to ensure favorable cleavage. We propose a mathematical formalism for epitope selection and ordering that balances the constraints imposed by these different conditions. Focusing on HCV, HIV, and Flu, we show that not all of the conditions can be satisfied for all viruses. Depending on the virus, different constraints are harder to fulfill: for Flu, the conservation constraint is violated first, while for HIV, it is difficult to focus the response at the optimal target antigens. The proposed methodology can be applied to any virus to assess the feasibility of optimally combining the above-mentioned constraints.

Statistics and truth in phylogenomics

Phylogenomics refers to the inference of historical relationships among species using genome-scale sequence data and to the use of phylogenetic analysis to infer protein function in multigene families. With rapidly decreasing sequencing costs, phylogenomics is becoming synonymous with evolutionary analysis of genome-scale and taxonomically densely sampled data sets. In phylogenetic inference applications, this translates into very large data sets that yield evolutionary and functional inferences with extremely small variances and high statistical confidence (P value). However, reports of highly significant P values are increasing even for contrasting phylogenetic hypotheses depending on the evolutionary model and inference method used, making it difficult to establish true relationships. We argue that the assessment of the robustness of results to biological factors, that may systematically mislead (bias) the outcomes of statistical estimation, will be a key to avoiding incorrect phylogenomic inferences. In fact, there is a need for increased emphasis on the magnitude of differences (effect sizes) in addition to the P values of the statistical test of the null hypothesis. On the other hand, the amount of sequence data available will likely always remain inadequate for some phylogenomic applications, for example, those involving episodic positive selection at individual codon positions and in specific lineages. Again, a focus on effect size and biological relevance, rather than the P value, may be warranted. Here, we present a theoretical overview and discuss practical aspects of the interplay between effect sizes, bias, and P values as it relates to the statistical inference of evolutionary truth in phylogenomics.

A random effects branch-site model for detecting episodic diversifying selection

Adaptive evolution frequently occurs in episodic bursts, localized to a few sites in a gene, and to a small number of lineages in a phylogenetic tree. A popular class of "branch-site" evolutionary models provides a statistical framework to search for evidence of such episodic selection. For computational tractability, current branch-site models unrealistically assume that all branches in the tree can be partitioned a priori into two rigid classes--"foreground" branches that are allowed to undergo diversifying selective bursts and "background" branches that are negatively selected or neutral. We demonstrate that this assumption leads to unacceptably high rates of false positives or false negatives when the evolutionary process along background branches strongly deviates from modeling assumptions. To address this problem, we extend Felsenstein's pruning algorithm to allow efficient likelihood computations for models in which variation over branches (and not just sites) is described in the random effects likelihood framework. This enables us to model the process at every branch-site combination as a mixture of three Markov substitution models--our model treats the selective class of every branch at a particular site as an unobserved state that is chosen independently of that at any other branch. When benchmarked on a previously published set of simulated sequences, our method consistently matched or outperformed existing branch-site tests in terms of power and error rates. Using three empirical data sets, previously analyzed for episodic selection, we discuss how modeling assumptions can influence inference in practical situations.

Two-stepping through time: mammals and viruses

Recent studies have identified ancient virus genomes preserved as fossils within diverse animal genomes. These fossils have led to the revelation that a broad range of mammalian virus families are older and more ubiquitous than previously appreciated. Long-term interactions between viruses and their hosts often develop into genetic arms races where both parties continually jockey for evolutionary dominance. It is difficult to imagine how mammalian hosts have kept pace in the evolutionary race against rapidly evolving viruses over large expanses of time, given their much slower evolutionary rates. However, recent data has begun to reveal the evolutionary strategy of slowly-evolving hosts. We review these data and suggest a modified arms race model where the evolutionary possibilities of viruses are relatively constrained. Such a model could allow more accurate forecasting of virus evolution.

Fluidity of HIV-1-specific T-cell responses during acute and early subtype C HIV-1 infection and associations with early disease progression

Deciphering immune events during early stages of human immunodeficiency virus type 1 (HIV-1) infection is critical for understanding the course of disease. We characterized the hierarchy of HIV-1-specific T-cell gamma interferon (IFN-γ) enzyme-linked immunospot (ELISPOT) assay responses during acute subtype C infection in 53 individuals and associated temporal patterns of responses with disease progression in the first 12 months. There was a diverse pattern of T-cell recognition across the proteome, with the recognition of Nef being immunodominant as early as 3 weeks postinfection. Over the first 6 months, we found that there was a 23% chance of an increased response to Nef for every week postinfection (P = 0.0024), followed by a nonsignificant increase to Pol (4.6%) and Gag (3.2%). Responses to Env and regulatory proteins appeared to remain stable. Three temporal patterns of HIV-specific T-cell responses could be distinguished: persistent, lost, or new. The proportion of persistent T-cell responses was significantly lower (P = 0.0037) in individuals defined as rapid progressors than in those progressing slowly and who controlled viremia. Almost 90% of lost T-cell responses were coincidental with autologous viral epitope escape. Regression analysis between the time to fixed viral escape and lost T-cell responses (r = 0.61; P = 0.019) showed a mean delay of 14 weeks after viral escape. Collectively, T-cell epitope recognition is not a static event, and temporal patterns of IFN-γ-based responses exist. This is due partly to viral sequence variation but also to the recognition of invariant viral epitopes that leads to waves of persistent T-cell immunity, which appears to associate with slower disease progression in the first year of infection.

Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology

Datamonkey is a popular web-based suite of phylogenetic analysis tools for use in evolutionary biology. Since the original release in 2005, we have expanded the analysis options to include recently developed algorithmic methods for recombination detection, evolutionary fingerprinting of genes, codon model selection, co-evolution between sites, identification of sites, which rapidly escape host-immune pressure and HIV-1 subtype assignment. The traditional selection tools have also been augmented to include recent developments in the field. Here, we summarize the analyses options currently available on Datamonkey, and provide guidelines for their use in evolutionary biology. Availability and documentation: http://www.datamonkey.org.

Evolutionary trends of A(H1N1) influenza virus hemagglutinin since 1918

The Pandemic (H1N1) 2009 is spreading to numerous countries and causing many human deaths. Although the symptoms in humans are mild at present, fears are that further mutations in the virus could lead to a potentially more dangerous outbreak in subsequent months. As the primary immunity-eliciting antigen, hemagglutinin (HA) is the major agent for host-driven antigenic drift in A(H3N2) virus. However, whether and how the evolution of HA is influenced by existing immunity is poorly understood for A(H1N1). Here, by analyzing hundreds of A(H1N1) HA sequences since 1918, we show the first evidence that host selections are indeed present in A(H1N1) HAs. Among a subgroup of human A(H1N1) HAs between 1918∼2008, we found strong diversifying (positive) selection at HA₁ 156 and 190. We also analyzed the evolutionary trends at HA₁ 190 and 225 that are critical determinants for receptor-binding specificity of A(H1N1) HA. Different A(H1N1) viruses appeared to favor one of these two sites in host-driven antigenic drift: epidemic A(H1N1) HAs favor HA₁ 190 while the 1918 pandemic and swine HAs favor HA₁ 225. Thus, our results highlight the urgency to understand the interplay between antigenic drift and receptor binding in HA evolution, and provide molecular signatures for monitoring future antigenically drifted 2009 pandemic and seasonal A(H1N1) influenza viruses.

Adaptive changes in HIV-1 subtype C proteins during early infection are driven by changes in HLA-associated immune pressure

It is unresolved whether recently transmitted human immunodeficiency viruses (HIV) have genetic features that specifically favour their transmissibility. To identify potential "transmission signatures", we compared 20 full-length HIV-1 subtype C genomes from primary infections, with 66 sampled from ethnically and geographically matched individuals with chronic infections. Controlling for recombination and phylogenetic relatedness, we identified 39 sites at which amino acid frequency spectra differed significantly between groups. These sites were predominantly located within Env, Pol and Gag (14/39, 9/39 and 6/39 respectively) and were significantly clustered (33/39) within known immunoreactive peptides. Within 6 months of infection, we detected reversion-to-consensus mutations at 14 sites and potential CTL escape mutations at seven. Here we provide evidence that frequent reversion mutations probably allows the virus to recover replicative fitness which, together with immune escape driven by the HLA alleles of the new hosts, differentiate sequences from chronic infections from those sampled shortly after transmission.

Evidence of HIV-1 adaptation to host HLA alleles following chimp-to-human transmission

BACKGROUND

The cytotoxic T-lymphocyte immune response is important in controlling HIV-1 replication in infected humans. In this immune pathway, viral peptides within infected cells are presented to T-lymphocytes by the polymorphic human leukocyte antigens (HLA). HLA alleles exert selective pressure on the peptide regions and immune escape mutations that occur at some of the targeted sites can enable the virus to adapt to the infected host. The pattern of ongoing immune escape and reversion associated with several human HLA alleles has been studied extensively. Such mutations revert upon transmission to a host without the HLA allele because the escape mutation incurs a fitness cost. However, to-date there has been little attempt to study permanent loss of CTL epitopes due to escape mutations without an effect on fitness.

RESULTS

Here, we set out to determine the extent of adaptation of HIV-1 to three well-characterized HLA alleles during the initial exposure of the virus to the human cytotoxic immune responses following transmission from chimpanzee. We generated a chimpanzee consensus sequence to approximate the virus sequence that was initially transmitted to the human host and used a method based on peptide binding affinity to HLA crystal structures to predict peptides that were potentially targeted by the HLA alleles on this sequence. Next, we used codon-based phylogenetic models to quantify the average selective pressure that acted on these regions during the period immediately following the zoonosis event, corresponding to the branch of the phylogenetic tree leading to the common ancestor of all of the HIV-1 sequences. Evidence for adaptive evolution during this period was observed at regions recognised by HLA A*6801 and A*0201, both of which are common in African populations. No evidence of adaptive evolution was observed at sites targeted by HLA-B*2705, which is a rare allele in African populations.

CONCLUSION

Our results suggest that the ancestral HIV-1 virus experienced a period of positive selective pressure due to immune responses associated with HLA alleles that were common in the infected human population. We propose that this resulted in permanent escape from immune responses targeting unconstrained regions of the virus.

HLA footprints on human immunodeficiency virus type 1 are associated with interclade polymorphisms and intraclade phylogenetic clustering

The selection of escape mutations has a major impact on immune control of infections with viruses such as human immunodeficiency virus (HIV). Viral evasion of CD8(+) T-cell responses leaves predictable combinations of escape mutations, termed HLA "footprints." The most clearly defined footprints are those associated with HLA alleles that are linked with successful control of HIV, such as HLA-B*57. Here we investigated the extent to which HLA footprint sites in HIV type 1 (HIV-1) are associated with viral evolution among and within clades. First, we examined the extent to which amino acid differences between HIV-1 clades share identity with sites of HLA-mediated selection pressure and observed a strong association, in particular with respect to sites of HLA-B selection (P < 10(-6)). Similarly, the sites of amino acid variability within a clade were found to overlap with sites of HLA-selected mutation. Second, we studied the impact of HLA selection on interclade phylogeny. Removing the sites of amino acid variability did not significantly affect clade-specific clustering, reflecting the central role of founder effects in establishing distinct clades. However, HLA footprints may underpin founder strains, and we show that amino acid substitutions between clades alter phylogeny, underlining a potentially substantial role for HLA in driving ongoing viral evolution. Finally, we investigated the impact of HLA selection on within-clade phylogeny and demonstrate that even a single HLA allele footprint can result in significant phylogenetic clustering of sequences. In conclusion, these data highlight the fact that HLA can be a strong selection force for both intra- and interclade HIV evolution at a population level.