Computational analysis of sequence selection mechanisms

Mechanisms leading to gene variations are responsible for the diversity of species and are important components of the theory of evolution. One constraint on gene evolution is that of protein foldability; the three-dimensional shapes of proteins must be thermodynamically stable. We explore the impact of this constraint and calculate properties of foldable sequences using 3660 structures from the Protein Data Bank. We seek a selection function that receives sequences as input, and outputs survival probability based on sequence fitness to structure. We compute the number of sequences that match a particular protein structure with energy lower than the native sequence, the density of the number of sequences, the entropy, and the "selection" temperature. The mechanism of structure selection for sequences longer than 200 amino acids is approximately universal. For shorter sequences, it is not. We speculate on concrete evolutionary mechanisms that show this behavior.     

Sequence variation of HIV and bioinformatics

The envelope glycoprotein of human immunodeficiency virus type 1 (HIV-1) interacts with receptors on the target cell and mediates virus entry by fusing the viral and cell membranes. To maintain the viral infectivity, amino acids that interact with receptors are expected to be more conserved than the other sites on the protein surface. In contrast to the functional constraint of amino acids for the receptor binding, some amino acid changes in this protein may produce antigenic variations that enable the virus to escape from recognition of the host immune system. Therefore, both positive selection (higher fitness) and negative selection (lower fitness) against amino acid changes are taking place during evolution of surface proteins of parasites To elucidate the evolutionary mechanisms of the whole HIV-1 gp120 envelope glycoprotein at the single site level, we collected and analyzed all available sequence data for the protein. By analyzing 186 sequences of the HIV-1 gp120 (subtype B), we reevaluated amino acid variability at the single site level, and estimated the numbers of synonymous and nonsynonymous substitutions at each codon position to detect positive and negative selection. We identified 33 amino acid positions which may be under positive selection. Some of these positions may form discontinuous epitopes. We also analyzed amino acid sequences to find amino acid positions responsible for usage of the second receptor. We found that, in addition to the V3 loop, amino acid variation at residue 440 in C4 region is clearly linked with the usage of CXCR 4.     

Protein evolution with dependence among codons due to tertiary structure

Markovian models of protein evolution that relax the assumption of independent change among codons are considered. With this comparatively realistic framework, an evolutionary rate at a site can depend both on the state of the site and on the states of surrounding sites. By allowing a relatively general dependence structure among sites, models of evolution can reflect attributes of tertiary structure. To quantify the impact of protein structure on protein evolution, we analyze protein-coding DNA sequence pairs with an evolutionary model that incorporates effects of solvent accessibility and pairwise interactions among amino acid residues. By explicitly considering the relationship between nonsynonymous substitution rates and protein structure, this approach can lead to refined detection and characterization of positive selection. Analyses of simulated sequence pairs indicate that parameters in this evolutionary model can be well estimated. Analyses of lysozyme c and annexin V sequence pairs yield the biologically reasonable result that amino acid replacement rates are higher when the replacements lead to energetically favorable proteins than when they destabilize the proteins. Although the focus here is evolutionary dependence among codons that is associated with protein structure, the statistical approach is quite general and could be applied to diverse cases of evolutionary dependence where surrogates for sequence fitness can be measured or modeled.     

Demographic processes affect HIV-1 evolution in primary infection before the onset of selective processes

HIV-1 transmission and viral evolution in the first year of infection were studied in 11 individuals representing four transmitter-recipient pairs and three independent seroconverters. Nine of these individuals were enrolled during acute infection; all were men who have sex with men (MSM) infected with HIV-1 subtype B. A total of 475 nearly full-length HIV-1 genome sequences were generated, representing on average 10 genomes per specimen at 2 to 12 visits over the first year of infection. Single founding variants with nearly homogeneous viral populations were detected in eight of the nine individuals who were enrolled during acute HIV-1 infection. Restriction to a single founder variant was not due to a lack of diversity in the transmitter as homogeneous populations were found in recipients from transmitters with chronic infection. Mutational patterns indicative of rapid viral population growth dominated during the first 5 weeks of infection and included a slight contraction of viral genetic diversity over the first 20 to 40 days. Subsequently, selection dominated, most markedly in env and nef. Mutants were detected in the first week and became consensus as early as day 21 after the onset of symptoms of primary HIV infection. We found multiple indications of cytotoxic T lymphocyte (CTL) escape mutations while reversions appeared limited. Putative escape mutations were often rapidly replaced with mutually exclusive mutations nearby, indicating the existence of a maturational escape process, possibly in adaptation to viral fitness constraints or to immune responses against new variants. We showed that establishment of HIV-1 infection is likely due to a biological mechanism that restricts transmission rather than to early adaptive evolution during acute infection. Furthermore, the diversity of HIV strains coupled with complex and individual-specific patterns of CTL escape did not reveal shared sequence characteristics of acute infection that could be harnessed for vaccine design.     

Epistasis Creates Invariant Sites and Modulates the Rate of Molecular Evolution

Invariant sites are a common feature of amino acid sequence evolution. The presence of invariant sites is frequently attributed to the need to preserve function through site-specific conservation of amino acid residues. Amino acid substitution models without a provision for invariant sites often fit the data significantly worse than those that allow for an excess of invariant sites beyond those predicted by models that only incorporate rate variation among sites (e.g., a Gamma distribution). An alternative is epistasis between sites to preserve residue interactions that can create invariant sites. Through computer-simulated sequence evolution, we evaluated the relative effects of site-specific preferences and site-site couplings in the generation of invariant sites and the modulation of the rate of molecular evolution. In an analysis of ten major families of protein domains with diverse sequence and functional properties, we find that the negative selection imposed by epistasis creates many more invariant sites than site-specific residue preferences alone. Further, epistasis plays an increasingly larger role in creating invariant sites over longer evolutionary periods. Epistasis also dictates rates of domain evolution over time by exerting significant additional purifying selection to preserve site couplings. These patterns illuminate the mechanistic role of epistasis in the processes underlying observed site invariance and evolutionary rates.     

Potent autologous and heterologous neutralizing antibody responses occur in HIV-2 infection across a broad range of infection outcomes

Few studies have explored the role of neutralizing antibody (NAb) responses in controlling HIV-2 viremia and disease progression. Using a TZM-bl neutralization assay, we assessed heterologous and autologous NAb responses from a community cohort of HIV-2-infected individuals with a broad range of disease outcomes in rural Guinea-Bissau. All subjects (n = 40) displayed exceptionally high heterologous NAb titers (50% inhibitory plasma dilution or 50% inhibitory concentration [IC(50)], 1:7,000 to 1:1,000,000) against 5 novel primary HIV-2 envelopes and HIV-2 7312A, whereas ROD A and 3 primary envelopes were relatively resistant to neutralization. Most individuals also showed high autologous NAb against contemporaneous envelopes (78% of plasma-envelope combinations in 69 envelopes from 21 subjects), with IC(50)s above 1:10,000. No association between heterologous or autologous NAb titer and greater control of HIV-2 was found. A subset of envelopes was found to be more resistant to neutralization (by plasma and HIV-2 monoclonal antibodies). These envelopes were isolated from individuals with greater intrapatient sequence diversity and were associated with changes in potential N-linked glycosylation sites but not CD4 independence or CXCR4 use. Plasma collected from up to 15 years previously was able to potently neutralize recent autologous envelopes, suggesting a lack of escape from NAb and the persistence of neutralization-sensitive variants over time, despite significant NAb pressure. We conclude that despite the presence of broad and potent NAb responses in HIV-2-infected individuals, these are not the primary forces behind the dichotomous outcomes observed but reveal a limited capacity for adaptive selection and escape from host immunity in HIV-2 infection.     

Simulating protein evolution in sequence and structure space

Naturally occurring proteins comprise a special subset of all plausible sequences and structures selected through evolution. Simulating protein evolution with simplified and all-atom models has shed light on the evolutionary dynamics of protein populations, the nature of evolved sequences and structures, and the extent to which today's proteins are shaped by selection pressures on folding, structure and function. Extensive mapping of the native structure, stability and folding rate in sequence space using lattice proteins has revealed organizational principles of the sequence/structure map important for evolutionary dynamics. Evolutionary simulations with lattice proteins have highlighted the importance of fitness landscapes, evolutionary mechanisms, population dynamics and sequence space entropy in shaping the generic properties of proteins. Finally, evolutionary-like simulations with all-atom models, in particular computational protein design, have helped identify the dominant selection pressures on naturally occurring protein sequences and structures.     

Immune-mediated positive selection drives human immunodeficiency virus type 1 molecular variation and predicts disease duration

Using likelihood-based evolutionary methods, we demonstrate that the broad genetic diversity of human immunodeficiency virus type 1 (HIV-1) in an infected individual is a consequence of site-specific positive selection for diversity, a likely consequence of immune recognition. In particular, the extent of positive selection appears to be a good predictor of disease duration. Positively selected sites along HIV-1 partial env sequences are numerous but not distributed uniformly. In a sample of eight patients studied longitudinally, the proportion of sites per sample under positive selection was a statistically significant predictor of disease duration. Among long-term progressors, positive selection persisted at sites over time and appears to be associated with helper T-cell epitopes. In contrast, sites under positive selection shifted from one longitudinal sample to the next in short-term progressors. Our study is consistent with the hypothesis that a broad and persistent immunologic response is associated with a slower rate of disease progression. In contrast, narrow, shifting immune responses characterize short-term progressors.     

Selective escape from CD8+ T-cell responses represents a major driving force of human immunodeficiency virus type 1 (HIV-1) sequence diversity and reveals constraints on HIV-1 evolution

The sequence diversity of human immunodeficiency virus type 1 (HIV-1) represents a major obstacle to the development of an effective vaccine, yet the forces impacting the evolution of this pathogen remain unclear. To address this issue we assessed the relationship between genome-wide viral evolution and adaptive CD8+ T-cell responses in four clade B virus-infected patients studied longitudinally for as long as 5 years after acute infection. Of the 98 amino acid mutations identified in nonenvelope antigens, 53% were associated with detectable CD8+ T-cell responses, indicative of positive selective immune pressures. An additional 18% of amino acid mutations represented substitutions toward common clade B consensus sequence residues, nine of which were strongly associated with HLA class I alleles not expressed by the subjects and thus indicative of reversions of transmitted CD8 escape mutations. Thus, nearly two-thirds of all mutations were attributable to CD8+ T-cell selective pressures. A closer examination of CD8 escape mutations in additional persons with chronic disease indicated that not only did immune pressures frequently result in selection of identical amino acid substitutions in mutating epitopes, but mutating residues also correlated with highly polymorphic sites in both clade B and C viruses. These data indicate a dominant role for cellular immune selective pressures in driving both individual and global HIV-1 evolution. The stereotypic nature of acquired mutations provides support for biochemical constraints limiting HIV-1 evolution and for the impact of CD8 escape mutations on viral fitness.     

Selection on the human immunodeficiency virus type 1 proteome following primary infection

Typically during human immunodeficiency virus type 1 (HIV-1) infection, a nearly homogeneous viral population first emerges and then diversifies over time due to selective forces that are poorly understood. To identify these forces, we conducted an intensive longitudinal study of viral genetic changes and T-cell immunity in one subject at < or =17 time points during his first 3 years of infection, and in his infecting partner near the time of transmission. Autologous peptides covering amino acid sites inferred to be under positive selection were powerful for identifying HIV-1-specific cytotoxic-T-lymphocyte (CTL) epitopes. Positive selection and mutations resulting in escape from CTLs occurred across the viral proteome. We detected 25 CTL epitopes, including 14 previously unreported. Seven new epitopes mapped to the viral Env protein, emphasizing Env as a major target of CTLs. One-third of the selected sites were associated with epitopic mutational escapes from CTLs. Most of these resulted from replacement with amino acids found at low database frequency. Another one-third represented acquisition of amino acids found at high database frequency, suggesting potential reversions of CTL epitopic sites recognized by the immune system of the transmitting partner and mutation toward improved viral fitness in the absence of immune targeting within the recipient. A majority of the remaining selected sites occurred in the envelope protein and may have been subjected to humoral immune selection. Hence, a majority of the amino acids undergoing selection in this subject appeared to result from fitness-balanced CTL selection, confirming CTLs as a dominant selective force in HIV-1 infection.     

Dynamic Correlation between Intrahost HIV-1 Quasispecies Evolution and Disease Progression

Quantifying the dynamics of intrahost HIV-1 sequence evolution is one means of uncovering information about the interaction between HIV-1 and the host immune system. In the chronic phase of infection, common dynamics of sequence divergence and diversity have been reported. We developed an HIV-1 sequence evolution model that simulated the effects of mutation and fitness of sequence variants. The amount of evolution was described by the distance from the founder strain, and fitness was described by the number of offspring a parent sequence produces. Analysis of the model suggested that the previously observed saturation of divergence and decrease of diversity in later stages of infection can be explained by a decrease in the proportion of offspring that are mutants as the distance from the founder strain increases rather than due to an increase of viral fitness. The prediction of the model was examined by performing phylogenetic analysis to estimate the change in the rate of evolution during infection. In agreement with our modeling, in 13 out of 15 patients (followed for 3–12 years) we found that the rate of intrahost HIV-1 evolution was not constant but rather slowed down at a rate correlated with the rate of CD4+ T-cell decline. The correlation between the dynamics of the evolutionary rate and the rate of CD4+ T-cell decline, coupled with our HIV-1 sequence evolution model, explains previously conflicting observations of the relationships between the rate of HIV-1 quasispecies evolution and disease progression.     

Recombination Enhances HIV-1 Envelope Diversity by Facilitating the Survival of Latent Genomic Fragments in the Plasma Virus Population

HIV-1 is subject to immune pressure exerted by the host, giving variants that escape the immune response an advantage. Virus released from activated latent cells competes against variants that have continually evolved and adapted to host immune pressure. Nevertheless, there is increasing evidence that virus displaying a signal of latency survives in patient plasma despite having reduced fitness due to long-term immune memory. We investigated the survival of virus with latent envelope genomic fragments by simulating within-host HIV-1 sequence evolution and the cycling of viral lineages in and out of the latent reservoir. Our model incorporates a detailed mutation process including nucleotide substitution, recombination, latent reservoir dynamics, diversifying selection pressure driven by the immune response, and purifying selection pressure asserted by deleterious mutations. We evaluated the ability of our model to capture sequence evolution in vivo by comparing our simulated sequences to HIV-1 envelope sequence data from 16 HIV-infected untreated patients. Empirical sequence divergence and diversity measures were qualitatively and quantitatively similar to those of our simulated HIV-1 populations, suggesting that our model invokes realistic trends of HIV-1 genetic evolution. Moreover, reconstructed phylogenies of simulated and patient HIV-1 populations showed similar topological structures. Our simulation results suggest that recombination is a key mechanism facilitating the persistence of virus with latent envelope genomic fragments in the productively infected cell population. Recombination increased the survival probability of latent virus forms approximately 13-fold. Prevalence of virus with latent fragments in productively infected cells was observed in only 2% of simulations when we ignored recombination, while the proportion increased to 27% of simulations when we allowed recombination. We also found that the selection pressures exerted by different fitness landscapes influenced the shape of phylogenies, diversity trends, and survival of virus with latent genomic fragments. Our model predicts that the persistence of latent genomic fragments from multiple different ancestral origins increases sequence diversity in plasma for reasonable fitness landscapes.     

Comparison of SIV and HIV-1 Genomic RNA Structures Reveals Impact of Sequence Evolution on Conserved and Non-Conserved Structural Motifs

RNA secondary structure plays a central role in the replication and metabolism of all RNA viruses, including retroviruses like HIV-1. However, structures with known function represent only a fraction of the secondary structure reported for HIV-1_NL4-3. One tool to assess the importance of RNA structures is to examine their conservation over evolutionary time. To this end, we used SHAPE to model the secondary structure of a second primate lentiviral genome, SIVmac239, which shares only 50% sequence identity at the nucleotide level with HIV-1_NL4-3. Only about half of the paired nucleotides are paired in both genomic RNAs and, across the genome, just 71 base pairs form with the same pairing partner in both genomes. On average the RNA secondary structure is thus evolving at a much faster rate than the sequence. Structure at the Gag-Pro-Pol frameshift site is maintained but in a significantly altered form, while the impact of selection for maintaining a protein binding interaction can be seen in the conservation of pairing partners in the small RRE stems where Rev binds. Structures that are conserved between SIVmac239 and HIV-1_NL4-3 also occur at the 5′ polyadenylation sequence, in the plus strand primer sites, PPT and cPPT, and in the stem-loop structure that includes the first splice acceptor site. The two genomes are adenosine-rich and cytidine-poor. The structured regions are enriched in guanosines, while unpaired regions are enriched in adenosines, and functionaly important structures have stronger base pairing than nonconserved structures. We conclude that much of the secondary structure is the result of fortuitous pairing in a metastable state that reforms during sequence evolution. However, secondary structure elements with important function are stabilized by higher guanosine content that allows regions of structure to persist as sequence evolution proceeds, and, within the confines of selective pressure, allows structures to evolve.     

Are residues in a protein folding nucleus evolutionarily conserved?

Protein is the working molecule of the cell, and evolution is the hallmark of life. It is important to understand how protein folding and evolution influence each other. Several studies correlating experimental measurement of residue participation in folding nucleus and sequence conservation have reached different conclusions. These studies are based on assessment of sequence conservation at folding nucleus sites using entropy or relative entropy measurement derived from multiple sequence alignment. Here we report analysis of conservation of folding nucleus using an evolutionary model alternative to entropy-based approaches. We employ a continuous time Markov model of codon substitution to distinguish mutation fixed by evolution and mutation fixed by chance. This model takes into account bias in codon frequency, bias-favoring transition over transversion, as well as explicit phylogenetic information. We measure selection pressure using the ratio omega of synonymous versus non-synonymous substitution at individual residue site. The omega-values are estimated using the PAML method, a maximum-likelihood estimator. Our results show that there is little correlation between the extent of kinetic participation in protein folding nucleus as measured by experimental phi-value and selection pressure as measured by omega-value. In addition, two randomization tests failed to show that folding nucleus residues are significantly more conserved than the whole protein, or the median omega value of all residues in the protein. These results suggest that at the level of codon substitution, there is no indication that folding nucleus residues are significantly more conserved than other residues. We further reconstruct candidate ancestral residues of the folding nucleus and suggest possible test tube mutation studies for testing folding behavior of ancient folding nucleus.     

Impaired RNA incorporation and dimerization in live attenuated leader-variants of SIVmac239

Background

The ability of emerging pathogens to infect new species is likely related to the diversity of pathogen variants present in existing reservoirs and their degree of genomic plasticity, which determines their ability to adapt to new environments. Certain simian immunodeficiency viruses (SIVcpz, SIVsm) have demonstrated tremendous success in infecting new species, including humans, resulting in the HIV-1 and HIV-2 epidemics. Although SIV diversification has been studied on a population level, the essential substrates for cross-species transmission, namely SIV sequence diversity and the types and extent of viral diversification present in individual reservoir animals have not been elucidated. To characterize this intra-host SIV diversity, we performed sequence analyses of clonal viral envelope (env) V1V2 and gag p27 variants present in individual SIVsm-infected sooty mangabeys over time.

Results

SIVsm demonstrated extensive intra-animal V1V2 length variation and amino acid diversity (le 38%), and continual variation in V1V2 N-linked glycosylation consensus sequence frequency and location. Positive selection was the predominant evolutionary force. Temporal sequence shifts suggested continual selection, likely due to evolving antibody responses. In contrast, gag p27 was predominantly under purifying selection. SIVsm V1V2 sequence diversification is at least as great as that in HIV-1 infected humans, indicating that extensive viral diversification in and of itself does not inevitably lead to AIDS.

Conclusion

Positive diversifying selection in this natural reservoir host is the engine that has driven the evolution of the uniquely adaptable SIV/HIV envelope protein. These studies emphasize the importance of retroviral diversification within individual host reservoir animals as a critical substrate in facilitating cross-species transmission.     

Predicted functional RNAs within coding regions constrain evolutionary rates of yeast proteins

Functional RNAs (fRNAs) are being recognized as an important regulatory component in biological processes. Interestingly, recent computational studies suggest that the number and biological significance of functional RNAs within coding regions (coding fRNAs) may have been underestimated. We hypothesized that such coding fRNAs will impose additional constraint on sequence evolution because the DNA primary sequence has to simultaneously code for functional RNA secondary structures on the messenger RNA in addition to the amino acid codons for the protein sequence. To test this prediction, we first utilized computational methods to predict conserved fRNA secondary structures within multiple species alignments of Saccharomyces sensu strico genomes. We predict that as much as 5% of the genes in the yeast genome contain at least one functional RNA secondary structure within their protein-coding region. We then analyzed the impact of coding fRNAs on the evolutionary rate of protein-coding genes because a decrease in evolutionary rate implies constraint due to biological functionality. We found that our predicted coding fRNAs have a significant influence on evolutionary rates (especially at synonymous sites), independent of other functional measures. Thus, coding fRNA may play a role on sequence evolution. Given that coding regions of humans and flies contain many more predicted coding fRNAs than yeast, the impact of coding fRNAs on sequence evolution may be substantial in genomes of higher eukaryotes.