Mapping sites of positive selection and amino acid diversification in the HIV genome: an alternative approach to vaccine design?

Changes in selective pressures associated with human population expansion may explain metabolic and immune related pathways enriched for signatures of positive selection

Background

The study of local adaptation processes is a very important research topic in the field of population genomics. There is a particular interest in the study of human populations because they underwent a process of rapid spatial expansion and faced important environmental changes that translated into changes in selective pressures. New mutations may have been selected for in the new environment and previously existing genetic variants may have become detrimental. Immune related genes may have been released from the selective pressure exerted by pathogens in the ancestral environment and new variants may have been positively selected due to pathogens present in the newly colonized habitat. Also, variants that had a selective advantage in past environments may have become deleterious in the modern world due to external stimuli including climatic, dietary and behavioral changes, which could explain the high prevalence of some polygenic diseases such as diabetes and obesity.

Results

We performed an enrichment analysis to identify gene sets enriched for signals of positive selection in humans. We used two genome scan methods, XPCLR and iHS to detect selection using a dense coverage of SNP markers combined with two gene set enrichment approaches. We identified immune related gene sets that could be involved in the protection against pathogens especially in the African population. We also identified the glycolysis & gluconeogenesis gene set, related to metabolism, which supports the thrifty genotype hypothesis invoked to explain the current high prevalence of diseases such as diabetes and obesity. Extending our analysis to the gene level, we found signals for 23 candidate genes linked to metabolic syndrome, 13 of which are new candidates for positive selection.

Conclusions

Our study provides a list of genes and gene sets associated with immunity and metabolic syndrome that are enriched for signals of positive selection in three human populations (Europeans, Africans and Asians). Our results highlight differences in the relative importance of pathogens as drivers of local adaptation in different continents and provide new insights into the evolution and high incidence of metabolic syndrome in modern human populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2783-2) contains supplementary material, which is available to authorized users.

Darwin and Fisher meet at biotech: on the potential of computational molecular evolution in industry

Background

Today computational molecular evolution is a vibrant research field that benefits from the availability of large and complex new generation sequencing data – ranging from full genomes and proteomes to microbiomes, metabolomes and epigenomes. The grounds for this progress were established long before the discovery of the DNA structure. Specifically, Darwin’s theory of evolution by means of natural selection not only remains relevant today, but also provides a solid basis for computational research with a variety of applications. But a long-term progress in biology was ensured by the mathematical sciences, as exemplified by Sir R. Fisher in early 20th century. Now this is true more than ever: The data size and its complexity require biologists to work in close collaboration with experts in computational sciences, modeling and statistics.

Results

Natural selection drives function conservation and adaptation to emerging pathogens or new environments; selection plays key role in immune and resistance systems. Here I focus on computational methods for evaluating selection in molecular sequences, and argue that they have a high potential for applications. Pharma and biotech industries can successfully use this potential, and should take the initiative to enhance their research and development with state of the art bioinformatics approaches.

Conclusions

This review provides a quick guide to the current computational approaches that apply the evolutionary principles of natural selection to real life problems – from drug target validation, vaccine design and protein engineering to applications in agriculture, ecology and conservation.

High genetic diversity and adaptive potential of two simian hemorrhagic fever viruses in a wild primate population

Key biological properties such as high genetic diversity and high evolutionary rate enhance the potential of certain RNA viruses to adapt and emerge. Identifying viruses with these properties in their natural hosts could dramatically improve disease forecasting and surveillance. Recently, we discovered two novel members of the viral family Arteriviridae: simian hemorrhagic fever virus (SHFV)-krc1 and SHFV-krc2, infecting a single wild red colobus (Procolobus rufomitratus tephrosceles) in Kibale National Park, Uganda. Nearly nothing is known about the biological properties of SHFVs in nature, although the SHFV type strain, SHFV-LVR, has caused devastating outbreaks of viral hemorrhagic fever in captive macaques. Here we detected SHFV-krc1 and SHFV-krc2 in 40% and 47% of 60 wild red colobus tested, respectively. We found viral loads in excess of 10(6)-10(7) RNA copies per milliliter of blood plasma for each of these viruses. SHFV-krc1 and SHFV-krc2 also showed high genetic diversity at both the inter- and intra-host levels. Analyses of synonymous and non-synonymous nucleotide diversity across viral genomes revealed patterns suggestive of positive selection in SHFV open reading frames (ORF) 5 (SHFV-krc2 only) and 7 (SHFV-krc1 and SHFV-krc2). Thus, these viruses share several important properties with some of the most rapidly evolving, emergent RNA viruses.

Natural polymorphisms and unusual mutations in HIV-1 protease with potential antiretroviral resistance: a bioinformatic analysis

Background

The correlations of genotypic and phenotypic tests with treatment, clinical history and the significance of mutations in viruses of HIV-infected patients are used to establish resistance mutations to protease inhibitors (PIs). Emerging mutations in human immunodeficiency virus type 1 (HIV-1) protease confer resistance to PIs by inducing structural changes at the ligand interaction site. The aim of this study was to establish an in silico structural relationship between natural HIV-1 polymorphisms and unusual HIV-1 mutations that confer resistance to PIs.

Results

Protease sequences isolated from 151 Mexican HIV-1 patients that were naïve to, or subjected to antiretroviral therapy, were examined. We identified 41 unrelated resistance mutations with a prevalence greater than 1%. Among these mutations, nine exhibited positive selection, three were natural polymorphisms (L63S/V/H) in a codon associated with drug resistance, and six were unusual mutations (L5F, D29V, L63R/G, P79L and T91V). The D29V mutation, with a prevalence of 1.32% in the studied population, was only found in patients treated with antiretroviral drugs. Using in silico modelling, we observed that D29V formed unstable protease complexes when were docked with lopinavir, saquinavir, darunavir, tipranavir, indinavir and atazanavir.

Conclusions

The structural correlation of natural polymorphisms and unusual mutations with drug resistance is useful for the identification of HIV-1 variants with potential resistance to PIs. The D29V mutation likely confers a selection advantage in viruses; however, in silico, presence of this mutation results in unstable enzyme/PI complexes, that possibly induce resistance to PIs.

Evolutionary conservation of mammalian sperm proteins associates with overall, not tyrosine, phosphorylation in human spermatozoa

We investigated possible associations between sequence evolution of mammalian sperm proteins and their phosphorylation status in humans. As a reference, spermatozoa from three normozoospermic men were analyzed combining two-dimensional gel electrophoresis, immunoblotting, and mass spectrometry. We identified 99 sperm proteins (thereof 42 newly described) and determined the phosphorylation status for most of them. Sequence evolution was studied across six mammalian species using nonsynonymous/synonymous rate ratios (dN/dS) and amino acid distances. Site-specific purifying selection was assessed employing average ratios of evolutionary rates at phosphorylated versus nonphosphorylated amino acids (α). According to our data, mammalian sperm proteins do not show statistically significant sequence conservation difference, no matter if the human ortholog is a phosphoprotein with or without tyrosine (Y) phosphorylation. In contrast, overall phosphorylation of human sperm proteins, i.e., phosphorylation at serine (S), threonine (T), and/or Y residues, associates with above-average conservation of sequences. Complementary investigations suggest that numerous protein-protein interactants constrain sequence evolution of sperm phosphoproteins. Although our findings reject a special relevance of Y phosphorylation for sperm functioning, they still indicate that overall phosphorylation substantially contributes to proper functioning of sperm proteins. Hence, phosphorylated sperm proteins might be considered as prime candidates for diagnosis and treatment of reduced male fertility.

HIV-1 Nef sequence and functional compartmentalization in the gut is not due to differential cytotoxic T lymphocyte selective pressure

The gut is the largest lymphoid organ in the body and a site of active HIV-1 replication and immune surveillance. The gut is a reservoir of persistent infection in some individuals with fully suppressed plasma viremia on combination antiretroviral therapy (cART) although the cause of this persistence is unknown. The HIV-1 accessory protein Nef contributes to persistence through multiple functions including immune evasion and increasing infectivity. Previous studies showed that Nef’s function is shaped by cytotoxic T lymphocyte (CTL) responses and that there are distinct populations of Nef within tissue compartments. We asked whether Nef’s sequence and/or function are compartmentalized in the gut and how compartmentalization relates to local CTL immune responses. Primary nef quasispecies from paired plasma and sigmoid colon biopsies from chronically infected subjects not on therapy were sequenced and cloned into Env⁻ Vpu⁻ pseudotyped reporter viruses. CTL responses were mapped by IFN-γ ELISpot using expanded CD8+ cells from blood and gut with pools of overlapping peptides covering the entire HIV proteome. CD4 and MHC Class I Nef-mediated downregulation was measured by flow cytometry. Multiple tests indicated compartmentalization of nef sequences in 5 of 8 subjects. There was also compartmentalization of function with MHC Class I downregulation relatively well preserved, but significant loss of CD4 downregulation specifically by gut quasispecies in 5 of 7 subjects. There was no compartmentalization of CTL responses in 6 of 8 subjects, and the selective pressure on quasispecies correlated with the magnitude CTL response regardless of location. These results demonstrate that Nef adapts via diverse pathways to local selective pressures within gut mucosa, which may be predominated by factors other than CTL responses such as target cell availability. The finding of a functionally distinct population within gut mucosa offers some insight into how HIV-1 may persist in the gut despite fully suppressed plasma viremia on cART.

The evolutionary significance of certain amino acid substitutions and their consequences for HIV-1 immunogenicity toward HLA's A*0201 and B*27

In silico tools are employed to examine the evolutionary relationship to possible vaccine peptide candidates' development. This perspective sheds light on the proteomic changes affecting the creation of HLA specific T-cell stimulating peptide vaccines for HIV. Full-length sequences of the envelope protein of the HIV subtypes A, B, C and D were obtained through the NCBI Protein database were aligned using CLUSTALW. They were then analyzed using RANKPEP specific to Human Leukocyte Antigen A*02 and B*27. Geneious was used to catalogue the collected gp160 sequences and to construct a phylogenic tree. Mesquite was employed for ancestral state reconstruction to infer the order of amino acid substitutions in the epitopes examined. The results showed that consensus peptide identified SLAEKNITI had changes that indicated predicted escape mutation in strains of HIV responding to pressure exerted by CD8+ cells expressing HLA A*02. The predominating 9-mers IRIGPGQAF of gp120 are significantly less immunogenic toward HLA B*27 than to HLA A*02. The data confirms previous findings on the importance for efficacious binding, of an arginine residue at the 2(nd) position of the gag SL9 epitope, and extends this principle to other epitopes which interacts with HLA B*27. This study shows that the understanding of viral evolution relating T-cell peptide vaccine design is a development that has much relevance for the creation of personalized therapeutics for HIV treatment.

Molecular characterization of HTLV-1 gp46 glycoprotein from health carriers and HAM/TSP infected individuals

BACKGROUND

Human T-cell Leukemia Virus type 1 (HTLV-1) is the etiological agent of tropical spastic paraparesis/HTLV-associated myelopathy (HAM/TSP) that can be identified in around 0.25%-3.8% of the infected population. Disease progression can be monitored by the proviral load and may depend on genetic factors, however, it is not well understood why some HTLV-1 infected people develop the disease while others do not. The present study attempts to assess the molecular diversity of gp46 glycoprotein in HAM/TSP patients and Health Carrier (HC) individuals.

METHODS

Blood samples were collected from 10 individuals, and DNA was extracted from PBMCs to measure the HTLV-1 proviral load. The gp46 coding sequences were amplified PCR, cloned and sequenced. The molecular characterization was performed using bioinformatics tools.

RESULTS

The median HTLV-1 proviral load of HC (n = 5) and HAM/TSP (n = 5) patients was similar (average 316,227 copies/106 PBMCs). The gp46 molecular characterization of 146 clones (70 HC and 76 HAM/TSP) revealed an overall diversity, within HC and HAM/TSP clones, of 0.4% and 0.6%, respectively. Five frequent mutations were detected among groups (HAM/TSP and HC clone sequences). A single amino acid (aa) substitution (S35L) was exclusive for the HC group, and three gp46 substitutions (F14S, N42H, G72S) were exclusive for the HAM/TSP group. The remaining frequent mutation (V247I) was present in both groups (p = 0.0014). The in silico protein analysis revealed that the mutated alleles F14S and N42H represent more hydrophilic and flexible protein domains that are likely to be less antigenic. The Receptor Binding Domain is quite variable in the HAM/TSP group. Two other domains (aa 53-75 and 175-209) that contain multiple linear T-cell epitopes showed genetic diversity in both HAM/TSP and HC groups. Further analysis revealed 27 and 13 T-cell epitopes for class I HLA alleles and class II HLA alleles, when analyzing the entire gp46.

CONCLUSIONS

The most common gp46 mutations were not associated clinical status because they were found in only one individual, except for the V247I mutation, that was found at viral clones from HAM/TSP ad HC individuals. Because of this, we cannot associate any of the gp46 found mutations with the clinical profile.

Evidence for positive selection in the extracellular domain of human cytomegalovirus encoded G protein-coupled receptor US28

The human cytomegalovirus (HCMV)-encoded chemokine receptor US28 is also a seven-transmembrane G-protein coupled receptor, whose signaling pathway is known for its involvement in host immune system evasion. HCMV infection can result in serious disease in immunocompromised individuals and is also linked to atherosclerosis and cardiovascular disease. Identifying amino acid residues that play a crucial role in successful viral adaptation in response to the host's immune defense is critical for effective drug design. In this study maximum likelihood-based codon substitution analyses were carried out to determine whether any codon of US28 has evolved adaptively. If the rate of nonsynonymous (dn) to the rate of synonymous (ds) nucleotide substitutions (ω = dn/ds) is greater than one, the codon is said to be under positive selection, indicating adaptive evolution. Although the overall ω for US28 gene was 0.154, indicating that most codon sites were subject to strong purifying selection, five codon sites are under strong positive selection. Three (E18D/L, D19A/E/G, and R267K/Q) of these positively selected sites are located in extracellular domains, the domains that play a crucial role for successful viral adaptation in response to the host's immune defense. The C-terminal (R329Q/W) and the fifth transmembrane domain (V190I), each have one positively selected site. These results suggest that relative to the extracellular domains, amino acid residues present in intracellular domains are more selectively constrained. A few amino acid residues in extracellular domains of US28 evolved more rapidly, presumably due to positive selection pressure resulting from ligand-binding and pathogen interactions of extracellular domains.

Complex positive selection pressures drive the evolution of HIV-1 with different co-receptor tropisms

HIV-1 co-receptor tropism is central for understanding the transmission and pathogenesis of HIV-1 infection. We performed a genome-wide comparison between the adaptive evolution of R5 and X4 variants from HIV-1 subtypes B and C. The results showed that R5 and X4 variants experienced differential evolutionary patterns and different HIV-1 genes encountered various positive selection pressures, suggesting that complex selection pressures are driving HIV-1 evolution. Compared with other hypervariable regions of Gp120, significantly more positively selected sites were detected in the V3 region of subtype B X4 variants, V2 region of subtype B R5 variants, and V1 and V4 regions of subtype C X4 variants, indicating an association of positive selection with co-receptor recognition/binding. Intriguingly, a significantly higher proportion (33.3% and 55.6%, P<0.05) of positively selected sites were identified in the C3 region than other conserved regions of Gp120 in all the analyzed HIV-1 variants, indicating that the C3 region might be more important to HIV-1 adaptation than previously thought. Approximately half of the positively selected sites identified in the env gene were identical between R5 and X4 variants. There were three common positively selected sites (96, 113 and 281) identified in Gp41 of all X4 and R5 variants from subtypes B and C. These sites might not only suggest a functional importance in viral survival and adaptation, but also imply a potential cross-immunogenicity between HIV-1 R5 and X4 variants, which has important implications for AIDS vaccine development.

Innovative approaches to develop prophylactic and therapeutic vaccines against HIV/AIDS

The acquired immunodeficiency syndrome (AIDS) emerged in the human population in the summer of 1981. According to the latest United Nations estimates, worldwide over 33 million people are infected with human immunodeficiency virus (HIV) and the prevalence rates continue to rise globally. To control the alarming spread of HIV, an urgent need exists for developing a safe and effective vaccine that prevents individuals from becoming infected or progressing to disease. To be effective, an HIV/AIDS vaccine should induce broad and long-lasting humoral and cellular immune responses, at both mucosal and systemic level. However, the nature of protective immune responses remains largely elusive and this represents one of the major roadblocks preventing the development of an effective vaccine. Here we summarize our present understanding of the factors responsible for resistance to infection or control of progression to disease in human and monkey that may be relevant to vaccine development and briefly review recent approaches which are currently being tested in clinical trials. Finally, the rationale and the current status of novel strategies based on nonstructural HIV-1 proteins, such as Tat, Nef and Rev, used alone or in combination with modified structural HIV-1 Env proteins are discussed.

An introduction to a novel population genetic approach for HIV characterization

The rapid evolution of the HIV genome is influenced in part by host selection pressure, which may cause parallel evolution among strains under shared selection pressures. To understand the mechanisms behind HIV-host immune escape across host populations, researchers have compared signatures of positive selection pressure on HIV codons across HIV subtypes and across phylogenetic groups of isolates within major subtypes, all relying on a criterion of phylogenetic separation. The HIV codon sites that retain diversity, evolve convergently among sets of hosts (cohorts) and diverge between cohorts may be phylogenetically undiagnostic (reveal little information about the relationship of the strains) and thus undetectable on a tree. We propose a new approach to characterizing genetic divergence among isolates using existing population genetic methods to better understand HIV response to host selection pressures. The approach combines population genetic statistical methods with codon analysis to identify putative amino acid sites evolving convergently. To illustrate the approach, we compared the C2-V3-C3 region of the envelope protein of HIV-1 clade B isolates between Haiti and USA hosts. This region showed no phylogenetic separation between host populations. Still, we identified codon sites in the C2-V3-C3 HIV-1 region that may have evolved differently between the two host populations. The sites are localized in human leukocyte antigen (HLA) class I binding epitopes, N-glycosylation motifs or both and are limited to the C2 and C3 regions. Our method provides a potential means to reveal candidate sites actively involved in HIV-1 immune escape that would otherwise be missed if a requisite for phylogenetic distinctiveness was made a priori. This strategy may prove to be a helpful way to characterize HIV genetic variation among hosts with suspected selection pressure differences, like progressors versus non-progressors.

Coalescent simulation of intracodon recombination

The coalescent with recombination is a very useful tool in molecular population genetics. Under this framework, genealogies often represent the evolution of the substitution unit, and because of this, the few coalescent algorithms implemented for the simulation of coding sequences force recombination to occur only between codons. However, it is clear that recombination is expected to occur most often within codons. Here we have developed an algorithm that can evolve coding sequences under an ancestral recombination graph that represents the genealogies at each nucleotide site, thereby allowing for intracodon recombination. The algorithm is a modification of Hudson's coalescent in which, in addition to keeping track of events occurring in the ancestral material that reaches the sample, we need to keep track of events occurring in ancestral material that does not reach the sample but that is produced by intracodon recombination. We are able to show that at typical substitution rates the number of nonsynonymous changes induced by intracodon recombination is small and that intracodon recombination does not generally result in inflated estimates of the overall nonsynonymous/synonymous substitution ratio (omega). On the other hand, recombination can bias the estimation of omega at particular codons, resulting in apparent rate variation among sites and in the spurious identification of positively selected sites. Importantly, in this case, allowing for variable synonymous rates across sites greatly reduces the false-positive rate and recovers statistical power. Finally, coalescent simulations with intracodon recombination could be used to better represent the evolution of nuclear coding genes or fast-evolving pathogens such as HIV-1.We have implemented this algorithm in a computer program called NetRecodon, freely available at http://darwin.uvigo.es.

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.

Dynamic features of the selective pressure on the human immunodeficiency virus type 1 (HIV-1) gp120 CD4-binding site in a group of long term non progressor (LTNP) subjects

The characteristics of intra-host human immunodeficiency virus type 1 (HIV-1) env evolution were evaluated in untreated HIV-1-infected subjects with different patterns of disease progression, including 2 normal progressor [NP], and 5 Long term non-progressor [LTNP] patients. High-resolution phylogenetic analysis of the C2-C5 env gene sequences of the replicating HIV-1 was performed in sequential samples collected over a 3–5 year period; overall, 301 HIV-1 genomic RNA sequences were amplified from plasma samples, cloned, sequenced and analyzed. Firstly, the evolutionary rate was calculated separately in the 3 codon positions. In all LTNPs, the 3^rd codon mutation rate was equal or even lower than that observed at the 1^st and 2^nd positions (p = 0.016), thus suggesting strong ongoing positive selection. A Bayesian approach and a maximum-likelihood (ML) method were used to estimate the rate of virus evolution within each subject and to detect positively selected sites respectively. A great number of N-linked glycosylation sites under positive selection were identified in both NP and LTNP subjects. Viral sequences from 4 of the 5 LTNPs showed extensive positive selective pressure on the CD4-binding site (CD4bs). In addition, localized pressure in the area of the IgG-b12 epitope, a broad neutralizing human monoclonal antibody targeting the CD4bs, was documented in one LTNP subject, using a graphic colour grade 3-dimensional visualization. Overall, the data shown here documenting high selective pressure on the HIV-1 CD4bs of a group of LTNP subjects offers important insights for planning novel strategies for the immune control of HIV-1 infection.

Extensive purifying selection acting on synonymous sites in HIV-1 Group M sequences

Background

Positive selection pressure acting on protein-coding sequences is usually inferred when the rate of nonsynonymous substitution is greater than the synonymous rate. However, purifying selection acting directly on the nucleotide sequence can lower the synonymous substitution rate. This could result in false inference of positive selection because when synonymous changes at some sites are under purifying selection, the average synonymous rate is an underestimate of the neutral rate of evolution. Even though HIV-1 coding sequences contain a number of regions that function at the nucleotide level, and are thus likely to be affected by purifying selection, studies of positive selection assume that synonymous substitutions can be used to estimate the neutral rate of evolution.

Results

We modelled site-to-site variation in the synonymous substitution rate across coding regions of the HIV-1 genome. Synonymous substitution rates were found to vary significantly within and between genes. Surprisingly, regions of the genome that encode proteins in more than one frame had significantly higher synonymous substitution rates than regions coding in a single frame. We found evidence of strong purifying selection pressure affecting synonymous mutations in fourteen regions with known functions. These included an exonic splicing enhancer, the rev-responsive element, the poly-purine tract and a transcription factor binding site. A further five highly conserved regions were located within known functional domains. We also found four conserved regions located in env and vpu which have not been characterized previously.

Conclusion

We provide the coordinates of genomic regions with markedly lower synonymous substitution rates, which are putatively under the influence of strong purifying selection pressure at the nucleotide level as well as regions encoding proteins in more than one frame. These regions should be excluded from studies of positive selection acting on HIV-1 coding regions.

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

Host immune responses against infectious pathogens exert strong selective pressures favouring the emergence of escape mutations that prevent immune recognition. Escape mutations within or flanking functionally conserved epitopes can occur at a significant cost to the pathogen in terms of its ability to replicate effectively. Such mutations come under selective pressure to revert to the wild type in hosts that do not mount an immune response against the epitope. Amino acid positions exhibiting this pattern of escape and reversion are of interest because they tend to coincide with immune responses that control pathogen replication effectively. We have used a probabilistic model of protein coding sequence evolution to detect sites in HIV-1 exhibiting a pattern of rapid escape and reversion. Our model is designed to detect sites that toggle between a wild type amino acid, which is susceptible to a specific immune response, and amino acids with lower replicative fitness that evade immune recognition. Through simulation, we show that this model has significantly greater power to detect selection involving immune escape and reversion than standard models of diversifying selection, which are sensitive to an overall increased rate of non-synonymous substitution. Applied to alignments of HIV-1 protein coding sequences, the model of immune escape and reversion detects a significantly greater number of adaptively evolving sites in env and nef. In all genes tested, the model provides a significantly better description of adaptively evolving sites than standard models of diversifying selection. Several of the sites detected are corroborated by association between Human Leukocyte Antigen (HLA) and viral sequence polymorphisms. Overall, there is evidence for a large number of sites in HIV-1 evolving under strong selective pressure, but exhibiting low sequence diversity. A phylogenetic model designed to detect rapid toggling between wild type and escape amino acids identifies a larger number of adaptively evolving sites in HIV-1, and can in some cases correctly identify the amino acid that is susceptible to the immune response.

Viruses, such as HIV, are able to evade host immune responses through escape mutations, yet sometimes they do so at a cost. This cost is the reduction in the ability of the virus to replicate, and thus selective pressure exists for a virus to revert to its original state in the absence of the host immune response that caused the initial escape mutation. This pattern of escape and reversion typically occurs when viruses are transmitted between individuals with different immune responses. We develop a phylogenetic model of immune escape and reversion and provide evidence that it outperforms existing models for the detection of selective pressure associated with host immune responses. Finally, we demonstrate that amino acid toggling is a pervasive process in HIV-1 evolution, such that many of the positions in the virus that evolve rapidly, under the influence of positive Darwinian selection, nonetheless display quite low sequence diversity. This highlights the limitations of HIV-1 evolution, and sites such as these are potentially good targets for HIV-1 vaccines.

Investigating selection on viruses: a statistical alignment approach

Background

Two problems complicate the study of selection in viral genomes: Firstly, the presence of genes in overlapping reading frames implies that selection in one reading frame can bias our estimates of neutral mutation rates in another reading frame. Secondly, the high mutation rates we are likely to encounter complicate the inference of a reliable alignment of genomes. To address these issues, we develop a model that explicitly models selection in overlapping reading frames. We then integrate this model into a statistical alignment framework, enabling us to estimate selection while explicitly dealing with the uncertainty of individual alignments. We show that in this way we obtain un-biased selection parameters for different genomic regions of interest, and can improve in accuracy compared to using a fixed alignment.

Results

We run a series of simulation studies to gauge how well we do in selection estimation, especially in comparison to the use of a fixed alignment. We show that the standard practice of using a ClustalW alignment can lead to considerable biases and that estimation accuracy increases substantially when explicitly integrating over the uncertainty in inferred alignments. We even manage to compete favourably for general evolutionary distances with an alignment produced by GenAl. We subsequently run our method on HIV2 and Hepatitis B sequences.

Conclusion

We propose that marginalizing over all alignments, as opposed to using a fixed one, should be considered in any parametric inference from divergent sequence data for which the alignments are not known with certainty. Moreover, we discover in HIV2 that double coding regions appear to be under less stringent selection than single coding ones. Additionally, there appears to be evidence for differential selection, where one overlapping reading frame is under positive and the other under negative selection.

Duplicated paralogous genes subject to positive selection in the genome of Trypanosoma brucei

Background

Whole genome studies have highlighted duplicated genes as important substrates for adaptive evolution. We have investigated adaptive evolution in this class of genes in the human parasite Trypanosoma brucei, as indicated by the ratio of non-synonymous (amino-acid changing) to synonymous (amino acid retaining) nucleotide substitution rates.

Methodology/Principal Findings

We have identified duplicated genes that are most rapidly evolving in this important human parasite. This is the first attempt to investigate adaptive evolution in this species at the codon level. We identify 109 genes within 23 clusters of paralogous gene expansions to be subject to positive selection.

Conclusions/Significance

Genes identified include surface antigens in both the mammalian and insect host life cycle stage suggesting that competitive interaction is not solely with the adaptive immune system of the mammalian host. Also surface transporters related to drug resistance and genes related to developmental progression are detected. We discuss how adaptive evolution of these genes may highlight lineage specific processes essential for parasite survival. We also discuss the implications of adaptive evolution of these targets for parasite biology and control.

Positive selection on HIV accessory proteins and the analysis of molecular adaptation after interspecies transmission

Studies examining positive selection on accessory proteins of HIV are rare, although these proteins play an important role in pathogenesis in vivo. Moreover, despite the biological relevance of analyses of molecular adaptation after viral transmission between species, the issue is still poorly studied. Here we present evidence that accessory proteins are subjected to positive selective forces exclusively in HIV. This scenario suggests that accessory protein genes are under adaptive evolution in HIV clades, while in SIVcpz such a phenomenon could not be detected. As a result, we show that comparative studies are critical to carry out functional investigation of positively selected protein sites, as they might help to achieve a better comprehension of the biology of HIV pathogenesis.

Reducing the false positive rate in the non-parametric analysis of molecular coevolution

Background

The strength of selective constraints operating on amino acid sites of proteins has a multifactorial nature. In fact, amino acid sites within proteins coevolve due to their functional and/or structural relationships. Different methods have been developed that attempt to account for the evolutionary dependencies between amino acid sites. Researchers have invested a significant effort to increase the sensitivity of such methods. However, the difficulty in disentangling functional co-dependencies from historical covariation has fuelled the scepticism over their power to detect biologically meaningful results. In addition, the biological parameters connecting linear sequence evolution to structure evolution remain elusive. For these reasons, most of the evolutionary studies aimed at identifying functional dependencies among protein domains have focused on the structural properties of proteins rather than on the information extracted from linear multiple sequence alignments (MSA). Non-parametric methods to detect coevolution have been reported to be especially susceptible to produce false positive results based on the properties of MSAs. However, no formal statistical analysis has been performed to definitively test the differential effects of these properties on the sensitivity of such methods.

Results

Here we test the effect that variations on the MSA properties have over the sensitivity of non-parametric methods to detect coevolution. We test the effect that the size of the MSA (number of sequences), mean pairwise amino acid distance per site and the strength of the coevolution signal have on the ability of non-parametric methods to detect coevolution. Our results indicate that all three factors have significant effects on the accuracy of non-parametric methods. Further, introducing statistical filters improves the sensitivity and increases the statistical power of the methods to detect functional coevolution. Statistical analysis of the physico-chemical properties of amino acid sites in the context of the protein structure reveals striking dependencies among amino acid sites. Results indicate a covariation trend in the hydrophobicities and molecular weight characteristics of amino acid sites when analysing a non-redundant set of 8000 protein structures. Using this biological information as filter in coevolutionary analyses minimises the false positive rate of these methods. Application of these filters to three different proteins with known functional domains supports the importance of using biological filters to detect coevolution.

Conclusion

Coevolutionary analyses using non-parametric methods have proved difficult and highly prone to provide spurious results depending on the properties of MSAs and on the strength of coevolution between amino acid sites. The application of statistical filters to the number of pairs detected as coevolving reduces significantly the number of artifactual results. Analysis of the physico-chemical properties of amino acid sites in the protein structure context reveals their structure-dependent covariation. The application of this known biological information to the analysis of covariation greatly enhances the functional coevolutionary signal and removes historical covariation. Simultaneous use of statistical and biological data is instrumental in the detection of functional amino acid sites dependencies and compensatory changes at the protein level.

Tat protein vaccination of cynomolgus macaques influences SHIV-89.6P cy243 epitope variability

In a previous study we showed that vaccination with the native Tat protein controlled virus replication in five out of seven monkeys against challenge with the simian human immunodeficiency virus (SHIV)-89.6P cy243 and that this protection correlated with T helper (Th)-1 response and cytotoxic T lymphocyte (CTL) activity. To address the evolution of the SHIV-89.6P cy243 both in control and vaccinated infected monkeys, the sequence of the human immunodeficiency virus (HIV)-1 Tat protein and the C2-V3 Env region of the proviral-DNA-derived clones were analyzed in both control and vaccinated but unprotected animals. We also performed analysis of the T cell epitope using a predictive epitope model taking into consideration the phylogeny of the variants. Our results suggest that even though the viral evolution observed in both groups of monkeys was directed toward variations in the major histocompatibility complex (MHC)-I epitopes, in the control animals it was associated with mutational escape of such epitopes. On the contrary, it is possible that viral evolution in the vaccinated monkeys was linked to mutations that arose to keep high the viral fitness. In the vaccinated animals the reduction of epitope variability, obtained prompting the immune system by vaccination and inducing a specific immunological response against virus, was able to reduce the emergence of escape mutants. Thus the intervention of host's selective forces in driving CTL escape mutants and in modulating viral fitness appeared to be different in the two groups of monkeys. We concluded that in the vaccinated unprotected animals, vaccination with the Tat protein induced a broad antiviral response, as demonstrated by the reduced ability to develop escape mutants, which is known to help in the control of viral replication.

Protein structure and oligomerization are important for the formation of export-competent HIV-1 Rev-RRE complexes

The translation of the unspliced and partially spliced viral mRNAs that encode the late, structural proteins of HIV-1 depends on the viral-protein Rev. Oligomeric binding of Rev to the Rev response element (RRE) in these mRNAs promotes their export from the nucleus and thus controls their expression. Here, we compared the effects of hydrophobic to hydrophilic mutations within the oligomerization domain of Rev using assays for oligomeric RNA binding, protein structure, and export from the nucleus. Oligomeric RNA binding alone does not correlate well with RNA transport activity in the subset of mutants. However, protein structure as judged by CD spectroscopy does correlate well with Rev function. The oligomeric assembly of Rev-L18T is impaired but exhibits minor defects in structure and retains a basal level of activity in vivo. The prevalence of L18T in infected individuals suggests a positive selection mechanism for L18T modulation of Rev activity that may delay the onset of AIDS.

Computational comparison of T-cell epitopes of gp120 of Iranian HIV-1 with different subtypes of the virus

In the present study, T-cell epitopes of gp120 of an Iranian isolate have been compared to different subtypes of HIV-1. At first, the amino acid sequences of gp120 were fetched from data banks. Then T-cell epitopes, disulfide bonding states, protein kinase C phosphorylation sites, cAMP-dependent protein kinase phosphorylation sites, casein kinase II phosphorylation sites, N-myristoylation sites and amidation sites were predicted using different soft wares. According to this computational analysis 6 good disulfide binding states in Iranian gp120 were predicted. From the viewpoint of cAMP-dependent protein kinase phosphorylation site (1 site) Iranian isolate was similar to clades B and F. Like subtype C 1 amidation site was predicted in the Iranian subtype. In the Iranian isolate 7 sites protein kinase C phosphorylation sites and 4 N-myristoylation sites were predicted. Since the number of individuals infected with HIV-1 in Iran, like many other countries is increasing, study of similarities and differences between the Iranian samples and different clades of HIV-1 can help us in identification of the origin and understanding the changes in the virus.

Annotation of selection strengths in viral genomes

MOTIVATION

Viral genomes tend to code in overlapping reading frames to maximize informational content. This may result in atypical codon bias and particular evolutionary constraints. Due to the fast mutation rate of viruses, there is additional strong evidence for varying selection between intra- and intergenomic regions. The presence of multiple coding regions complicates the concept of K(a)/K(s) ratio, and thus begs for an alternative approach when investigating selection strengths. Building on the paper by McCauley and Hein, we develop a method for annotating a viral genome coding in overlapping reading frames. We introduce an evolutionary model capable of accounting for varying levels of selection along the genome, and incorporate it into our prior single sequence HMM methodology, extending it now to a phylogenetic HMM. Given an alignment of several homologous viruses to a reference sequence, we may thus achieve an annotation both of coding regions as well as selection strengths, allowing us to investigate different selection patterns and hypotheses.

RESULTS

We illustrate our method by applying it to a multiple alignment of four HIV2 sequences, as well as of three Hepatitis B sequences. We obtain an annotation of the coding regions, as well as a posterior probability for each site of the strength of selection acting on it. From this we may deduce the average posterior selection acting on the different genes. Whilst we are encouraged to see in HIV2, that the known to be conserved genes gag and pol are indeed annotated as such, we also discover several sites of less stringent negative selection within the env gene. To the best of our knowledge, we are the first to subsequently provide a full selection annotation of the Hepatitis B genome by explicitly modelling the evolution within overlapping reading frames, and not relying on simple K(a)/K(s) ratios.

Towards realistic codon models: among site variability and dependency of synonymous and non-synonymous rates

Codon evolutionary models are widely used to infer the selection forces acting on a protein. The non-synonymous to synonymous rate ratio (denoted by Ka/Ks) is used to infer specific positions that are under purifying or positive selection. Current evolutionary models usually assume that only the non-synonymous rates vary among sites while the synonymous substitution rates are constant. This assumption ignores the possibility of selection forces acting at the DNA or mRNA levels. Towards a more realistic description of sequence evolution, we present a model that accounts for among-site-variation of both synonymous and non-synonymous substitution rates. Furthermore, we alleviate the widespread assumption that positions evolve independently of each other. Thus, possible sources of bias caused by random fluctuations in either the synonymous or non-synonymous rate estimations at a single site is removed. Our model is based on two hidden Markov models that operate on the spatial dimension: one describes the dependency between adjacent non-synonymous rates while the other describes the dependency between adjacent synonymous rates. The presented model is applied to study the selection pressure across the HIV-1 genome. The new model better describes the evolution of all HIV-1 genes, as compared to current codon models. Using both simulations and real data analyses, we illustrate that accounting for synonymous rate variability and dependency greatly increases the accuracy of Ka/Ks estimation and in particular of positively selected sites. Finally, we discuss the applicability of the developed model to infer the selection forces in regulatory and overlapping regions of the HIV-1 genome.

Modification of chicken avian β-defensin-8 at positively selected amino acid sites enhances specific antimicrobial activity

Antimicrobial peptides (AMPs), essential components of innate immunity, are found in a range of phylogenetically diverse species and are thought to act by disrupting the membrane integrity of microbes. In this paper, we used evolutionary signatures to identify sites that are most relevant during the functional evolution of these molecules and introduced amino acid substitutions to improve activity. We first demonstrate that the anti-microbial activity of chicken avian beta-defensin-8, previously known as gallinacin-12, can be significantly increased against Escherichia coli, Listeria monocytogenes, Salmonella typhimurium, Salmonella typhimurium phoP- mutant and Streptococcus pyogenes through targeted amino acid substitutions, which confer increased peptide charge. However, by increasing the AMP charge through amino acid substitutions at sites predicted to be subject to positive selection, antimicrobial activity against Escherichia coli was further increased. In contrast, no further increase in activity was observed against the remaining pathogens. This result suggests that charge-increasing modifications confer increased broad-spectrum activity to an AMP, whilst positive selection at particular sites is involved in directing the antimicrobial response against specific pathogens. Thus, there is potential for the rational design of novel therapeutics based on specifically targeted and modified AMPs.

Comparative annotation of viral genomes with non-conserved gene structure

MOTIVATION

Detecting genes in viral genomes is a complex task. Due to the biological necessity of them being constrained in length, RNA viruses in particular tend to code in overlapping reading frames. Since one amino acid is encoded by a triplet of nucleic acids, up to three genes may be coded for simultaneously in one direction. Conventional hidden Markov model (HMM)-based gene-finding algorithms may typically find it difficult to identify multiple coding regions, since in general their topologies do not allow for the presence of overlapping or nested genes. Comparative methods have therefore been restricted to likelihood ratio tests on potential regions as to being double or single coding, using the fact that the constrictions forced upon multiple-coding nucleotides will result in atypical sequence evolution. Exploiting these same constraints, we present an HMM based gene-finding program, which allows for coding in unidirectional nested and overlapping reading frames, to annotate two homologous aligned viral genomes. Our method does not insist on conserved gene structure between the two sequences, thus making it applicable for the pairwise comparison of more distantly related sequences.

RESULTS

We apply our method to 15 pairwise alignments of six different HIV2 genomes. Given sufficient evolutionary distance between the two sequences, we achieve sensitivity of approximately 84-89% and specificity of approximately 97-99.9%. We additionally annotate three pairwise alignments of the more distantly related HIV1 and HIV2, as well as of two different hepatitis viruses, attaining results of approximately 87% sensitivity and approximately 98.5% specificity. We subsequently incorporate prior knowledge by 'knowing' the gene structure of one sequence and annotating the other conditional on it. Boosting accuracy close to perfect we demonstrate that conservation of gene structure on top of nucleotide sequence is a valuable source of information, especially in distantly related genomes.

AVAILABILITY

The Java code is available from the authors.

Building collaborative networks for HIV/AIDS vaccine development: the AVIP experience

The need for an effective HIV/AIDS vaccine is imperative to halt a pandemic that involves more than 40 million individuals worldwide as of 2005 and is causing enormous socio-economic losses, especially in developing countries (DC). The overall failure of more than two decades of HIV vaccine research justifies the demands for a concerted effort for the rapid development of new and efficacious vaccines against HIV/AIDS. In this context, building international collaborative networks is a must for speeding up scientific research and optimizing the use of funding in a synergistic fashion, as resources for HIV/AIDS are limited and do not involve most of the biggest Pharmas that are more interested in drug discovery. The AIDS Vaccine Integrated Project (AVIP) consortium is an example of synergistic partnership of international European Union and DC experts with a common research goal. AVIP is a European Commission-funded (FP-6), consortium-based, 5-year program directed to the fast development of new HIV/AIDS vaccine candidates to be tested in phase I clinical trials in Europe for future advancement to phase II/III testing in DC. To ensure their rapid development, AVIP novel combined vaccines include both regulatory and structural HIV antigens, which have already been tested, as single components, in phase I clinical trials. In particular, such combination vaccines may be superior to earlier vaccine candidates, the vast majority of which are based only on either structural or regulatory HIV products. In fact, the generation of immune responses to both types of viral antigens expressed either early (regulatory products) or late (structural products) during the viral life cycle can maximize immune targeting of both primary or chronic viral infection. Further, the rational design of combined vaccines allows exploitation of immunomodulatory functions of HIV regulatory proteins, which can improve immunity against structural vaccine components. The building of the AVIP consortium and its scientific strategy will be reviewed in this paper as an example of the establishment of a consortium regulated by a specific intellectual property agreement.

Adaptation to different human populations by HIV-1 revealed by codon-based analyses

Several codon-based methods are available for detecting adaptive evolution in protein-coding sequences, but to date none specifically identify sites that are selected differentially in two populations, although such comparisons between populations have been historically useful in identifying the action of natural selection. We have developed two fixed effects maximum likelihood methods: one for identifying codon positions showing selection patterns that persist in a population and another for detecting whether selection is operating differentially on individual codons of a gene sampled from two different populations. Applying these methods to two HIV populations infecting genetically distinct human hosts, we have found that few of the positively selected amino acid sites persist in the population; the other changes are detected only at the tips of the phylogenetic tree and appear deleterious in the long term. Additionally, we have identified seven amino acid sites in protease and reverse transcriptase that are selected differentially in the two samples, demonstrating specific population-level adaptation of HIV to human populations.

Despite the efforts devoted to surveying HIV genetic diversity and the development of an effective vaccine, there is still no consensus on the extent to which the former prejudices the latter. Experimental studies show that escape from cell-mediated immunity is selected for in HIV and SIV, and sometimes very quickly. Conversely, escape mutants may be selected against at transmission, so how much does this selection within individuals influence the genotype of the circulating HIV population overall? Kosakovsky Pond, Leigh Brown, and colleagues have developed a new statistical approach to address this question. Using sequences from the globally most abundant HIV subtype (subtype C), the authors directly compared virus of the same subtype infecting genetically different human populations. They show at least half of the amino acid sites selected within individuals are not selected at a population level, and they identify six amino acid sites in the RT gene that are selected differentially between populations. We can now partition molecular adaptation between individual and population components for whatever genes may be included in candidate vaccines, in the target populations themselves. The methods are also important beyond the HIV world, where analogous issues arise in the more general question of the evolution of virulence in pathogens.

Identification of a positively evolving putative binding region with increased variability in posttranslational motifs in zonadhesin MAM domain 2

Positive selection has been shown to be pervasive in sex-related proteins of many metazoan taxa. However, we are only beginning to understand molecular evolutionary processes on the lineage to humans. To elucidate the evolution of proteins involved in human reproduction, we studied the sequence evolution of MAM domains of the sperm-ligand zonadhesin in respect to single amino acid sites, solvent accessibility, and posttranslational modification. GenBank-data were supplemented by new cDNA-sequences of a representative non-human primate panel. Solvent accessibility predictions identified a probably exposed fragment of 30 amino acids belonging to MAM domain 2 (i.e., MAM domain 3 in mouse). The fragment is characterized by significantly increased rate of positively selected amino acid sites and exhibits high variability in predicted posttranslational modification, and, thus, might represent a binding region in the mature protein. At the same time, there is a significant coincidence of positively selected amino acid sites and non-conserved posttranslational motifs. We conclude that the binding specificity of zonadhesin MAM domains, especially of the presumed epitope, is achieved by positive selection at the level of single amino acid sites and posttranslational modifications, respectively.

Bovine seminal plasma proteins and their relatives: A new expanding superfamily in mammals

BSP proteins represent three major proteins of bovine seminal plasma: BSP-A1/-A2, -A3 and -30 kDa. The BSP protein signature is characterized by two tandemly repeated fibronectin type 2 (Fn2) domains. Although classical affinity chromatography and protein sequencing have proven that the BSP protein homologs may be ubiquitous in mammals and functionally related to sperm capacitation, only the three bovine genes have been reported thus far. In this study, we report three new BSP protein-related genes from bovine, as well as other BSP protein-related DNA sequences from human, chimpanzee, mouse, rat, dog, horse and rabbit. Analysis of the relationships between all Fn2 domain-containing proteins revealed that the Fn2 domains found in BSP-related proteins have special features that distinguish them from non-BSP-related proteins. These features can be used to identify new BSP protein-related sequences. Further molecular evolutionary analysis of the BSP protein lineage revealed that all BSP proteins and their related sequences can be grouped into three subfamilies: BSPH4, BSPH5 and BSPH6, which indicates that the BSP protein family is much bigger than previously envisioned. More interestingly, the three BSP proteins in bovine within the BSPH4-subfamily were shown to evolve rapidly. The ratio of nonsynonymous to synonymous substitutions was higher than 1. The analysis also indicated that the rate of evolution was heterogeneous between the first and second Fn2 domains of the genes. These data may reflect that some amino acids in BSP proteins are under a strong positive selection after gene duplication and that each BSP protein evolves rapidly, possibly to acquire new functions.

Estimating diversifying selection and functional constraint in the presence of recombination

Models of molecular evolution that incorporate the ratio of nonsynonymous to synonymous polymorphism (dN/dS ratio) as a parameter can be used to identify sites that are under diversifying selection or functional constraint in a sample of gene sequences. However, when there has been recombination in the evolutionary history of the sequences, reconstructing a single phylogenetic tree is not appropriate, and inference based on a single tree can give misleading results. In the presence of high levels of recombination, the identification of sites experiencing diversifying selection can suffer from a false-positive rate as high as 90%. We present a model that uses a population genetics approximation to the coalescent with recombination and use reversible-jump MCMC to perform Bayesian inference on both the dN/dS ratio and the recombination rate, allowing each to vary along the sequence. We demonstrate that the method has the power to detect variation in the dN/dS ratio and the recombination rate and does not suffer from a high false-positive rate. We use the method to analyze the porB gene of Neisseria meningitidis and verify the inferences using prior sensitivity analysis and model criticism techniques.

Neutralizing antibody responses drive the evolution of human immunodeficiency virus type 1 envelope during recent HIV infection

HIV type 1 (HIV-1) can rapidly escape from neutralizing antibody responses. The genetic basis of this escape in vivo is poorly understood. We compared the pattern of evolution of the HIV-1 env gene between individuals with recent HIV infection whose virus exhibited either a low or a high rate of escape from neutralizing antibody responses. We demonstrate that the rate of viral escape at a phenotypic level is highly variable among individuals, and is strongly correlated with the rate of amino acid substitutions. We show that dramatic escape from neutralizing antibodies can occur in the relative absence of changes in glycosylation or insertions and deletions ("indels") in the envelope; conversely, changes in glycosylation and indels occur even in the absence of neutralizing antibody responses. Comparison of our data with the predictions of a mathematical model support a mechanism in which escape from neutralizing antibodies occurs via many amino acid substitutions, with low cross-neutralization between closely related viral strains. Our results suggest that autologous neutralizing antibody responses may play a pivotal role in the diversification of HIV-1 envelope during the early stages of infection.

Sequence evolution, processing, and posttranslational modification of zonadhesin D domains in primates, as inferred from cDNA data

Zonadhesin is a mammalian transmembrane sperm ligand. Precursor zonadhesin essentially consists of MAM (meprin/A5 antigen/mu receptor tyrosine phosphatase) domains, a mucin-like repeat, and D domains (homologous to von Willebrand D). Recent immunovisualization and binding assays indicate that zonadhesin D domains 1-3 bind postacrosomally to the zona pellucida. This feature has attracted considerable interest in the evolution of zonadhesin and its possible biological and biomedical implications. Previous molecular evolutionary analyses, however, were confined to cDNA sequences of only few distantly related species. Moreover, except for rabbit and pig, little is known about zonadhesin's processing. To delineate the situation in primates including humans, we analyze here the evolution of zonadhesin on the basis of D domain encoding cDNAs of about 4900 base pairs (bp) length from a representative primate sampling (1 Strepsirhini, 3 Cercopithecidae, 3 Platyrrhini, and human; 7 new sequences) plus GenBank data from mouse, rabbit, and pig. Site-specific (CODEML and HyPhy) analysis indicates positive evolution of zonadhesin. Moreover, moving window analysis (CRANN) points to a positive correlation of sequence evolution and sperm-competition. Significant accumulations of positively selected sites across interspecifically variable motifs (identified by PROSITE) suggest that positive selection promotes differences between species by amino acid exchanges and changes in posttranslational modification. In the case of zonadhesin D domains, positive selection might thus contribute to the species-specific binding of zonadhesin and zona pellucida. A high conservation of processing and dimerization motifs of primate zonadhesin in analogy to pig, on the other hand, illustrates that zonadhesin's backbone needs to meet basic requirements in order to retain function.

Evidence of positive Darwinian selection in putative meningococcal vaccine antigens

Meningococcal meningitidis is a life-threatening disease. In Europe and the United States the majority of cases are caused by virulent meningococcal strains belonging to serogroup B. Presently there is no effective vaccine against serogroup B strains, as traditional vaccine antigens such as polysaccharide capsules are unusable as they lead to autoimmunity. The year 2000 saw the publication of the complete genome of Neisseria meningitidis MC58, a virulent serogroup B bacterium. Working in conjunction with the sequencing project, researchers endeavored to locate highly conserved membrane-associated proteins that elicit an immune response. It is hoped that these proteins will provide a basis for novel vaccines against serogroup B strains. A number of potential vaccine antigens have been located and are presently in phase I clinical trials. Recently many reports pertaining to the evidence of positive Darwinian selection in membrane proteins of pathogens have been reported. This study utilized in silico methods to test for evidence of historical positive Darwinian selection in seven such vaccine candidates. We found that two of these proteins show signatures of adaptive evolution, while the remaining proteins show evidence of strong purifying selection. This has significant implications for the design of a vaccine against serogroup B strains, as it has been shown that vaccines that target epitopes that are under strong purifying selection are better than those that target variable epitopes.

BioAfrica's HIV-1 proteomics resource: combining protein data with bioinformatics tools

Most Internet online resources for investigating HIV biology contain either bioinformatics tools, protein information or sequence data. The objective of this study was to develop a comprehensive online proteomics resource that integrates bioinformatics with the latest information on HIV-1 protein structure, gene expression, post-transcriptional/post-translational modification, functional activity, and protein-macromolecule interactions. The BioAfrica HIV-1 Proteomics Resource http://bioafrica.mrc.ac.za/proteomics/index.html is a website that contains detailed information about the HIV-1 proteome and protease cleavage sites, as well as data-mining tools that can be used to manipulate and query protein sequence data, a BLAST tool for initiating structural analyses of HIV-1 proteins, and a proteomics tools directory. The Proteome section contains extensive data on each of 19 HIV-1 proteins, including their functional properties, a sample analysis of HIV-1HXB2, structural models and links to other online resources. The HIV-1 Protease Cleavage Sites section provides information on the position, subtype variation and genetic evolution of Gag, Gag-Pol and Nef cleavage sites. The HIV-1 Protein Data-mining Tool includes a set of 27 group M (subtypes A through K) reference sequences that can be used to assess the influence of genetic variation on immunological and functional domains of the protein. The BLAST Structure Tool identifies proteins with similar, experimentally determined topologies, and the Tools Directory provides a categorized list of websites and relevant software programs. This combined database and software repository is designed to facilitate the capture, retrieval and analysis of HIV-1 protein data, and to convert it into clinically useful information relating to the pathogenesis, transmission and therapeutic response of different HIV-1 variants. The HIV-1 Proteomics Resource is readily accessible through the BioAfrica website at: http://bioafrica.mrc.ac.za/proteomics/index.html.

Nonstructural HIV proteins as targets for prophylactic or therapeutic vaccines

By the end of 2004, more than 20 HIV-1 vaccine candidates will have entered clinical testing in at least 30 trials worldwide. Almost half of these vaccines include nonstructural HIV-1 gene products. This represents an important innovation in the HIV vaccine field, because until 9 years ago not even preclinical testing in small animal models had been carried out with such immunogens. This review briefly discusses the experimental evidence that provides the rationale for the use of nonstructural HIV-1 gene products as vaccine antigens, and summarizes the current status and the future development of these novel vaccines.