Evolutionary interactions between N-linked glycosylation sites in the HIV-1 envelope

The addition of asparagine (N)-linked polysaccharide chains (i.e., glycans) to the gp120 and gp41 glycoproteins of human immunodeficiency virus type 1 (HIV-1) envelope is not only required for correct protein folding, but also may provide protection against neutralizing antibodies as a “glycan shield.” As a result, strong host-specific selection is frequently associated with codon positions where nonsynonymous substitutions can create or disrupt potential N-linked glycosylation sites (PNGSs). Moreover, empirical data suggest that the individual contribution of PNGSs to the neutralization sensitivity or infectivity of HIV-1 may be critically dependent on the presence or absence of other PNGSs in the envelope sequence. Here we evaluate how glycan–glycan interactions have shaped the evolution of HIV-1 envelope sequences by analyzing the distribution of PNGSs in a large-sequence alignment. Using a “covarion”-type phylogenetic model, we find that the rates at which individual PNGSs are gained or lost vary significantly over time, suggesting that the selective advantage of having a PNGS may depend on the presence or absence of other PNGSs in the sequence. Consequently, we identify specific interactions between PNGSs in the alignment using a new paired-character phylogenetic model of evolution, and a Bayesian graphical model. Despite the fundamental differences between these two methods, several interactions are jointly identified by both. Mapping these interactions onto a structural model of HIV-1 gp120 reveals that negative (exclusive) interactions occur significantly more often between colocalized glycans, while positive (inclusive) interactions are restricted to more distant glycans. Our results imply that the adaptive repertoire of alternative configurations in the HIV-1 glycan shield is limited by functional interactions between the N-linked glycans. This represents a potential vulnerability of rapidly evolving HIV-1 populations that may provide useful glycan-based targets for neutralizing antibodies.     

Inconsistencies in estimating the age of HIV-1 subtypes due to heterotachy

Rate heterogeneity among lineages is a common feature of molecular evolution, and it has long impeded our ability to accurately estimate the age of evolutionary divergence events. The development of relaxed molecular clocks, which model variable substitution rates among lineages, was intended to rectify this problem. Major subtypes of pandemic HIV-1 group M are thought to exemplify closely related lineages with different substitution rates. Here, we report that inferring the time of most recent common ancestor of all these subtypes in a single phylogeny under a single (relaxed) molecular clock produces significantly different dates for many of the subtypes than does analysis of each subtype on its own. We explore various methods to ameliorate this problem. We conclude that current molecular dating methods are inadequate for dealing with this type of substitution rate variation in HIV-1. Through simulation, we show that heterotachy causes root ages to be overestimated.     

Automated phylogenetic detection of recombination using a genetic algorithm

The evolution of homologous sequences affected by recombination or gene conversion cannot be adequately explained by a single phylogenetic tree. Many tree-based methods for sequence analysis, for example, those used for detecting sites evolving nonneutrally, have been shown to fail if such phylogenetic incongruity is ignored. However, it may be possible to propose several phylogenies that can correctly model the evolution of nonrecombinant fragments. We propose a model-based framework that uses a genetic algorithm to search a multiple-sequence alignment for putative recombination break points, quantifies the level of support for their locations, and identifies sequences or clades involved in putative recombination events. The software implementation can be run quickly and efficiently in a distributed computing environment, and various components of the methods can be chosen for computational expediency or statistical rigor. We evaluate the performance of the new method on simulated alignments and on an array of published benchmark data sets. Finally, we demonstrate that prescreening alignments with our method allows one to analyze recombinant sequences for positive selection.     

<Emphasis Type="Italic">α</Emphasis>-Amylase inhibitor-1 gene from <Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> expressed in <Emphasis Type="Italic">Coffea arabica</Emphasis>plants inhibits α-amylases from the coffee berry borer pest

Background  

Coffee is an important crop and is crucial to the economy of many developing countries, generating around US70 billion per year. There are 115 species in the < i > Coffea < /i > genus, but only two, < i > C. arabica < /i > and < i > C. canephora < /i > , are commercially cultivated. Coffee plants are attacked by many pathogens and insect-pests, which affect not only the production of coffee but also its grain quality, reducing the commercial value of the product. The main insect-pest, the coffee berry borer ( < i > Hypotheneumus hampei < /i > ), is responsible for worldwide annual losses of around US70 billion per year. There are 115 species in the Coffea genus, but only two, C. arabica and C. canephora, are commercially cultivated. Coffee plants are attacked by many pathogens and insect-pests, which affect not only the production of coffee but also its grain quality, reducing the commercial value of the product. The main insect-pest, the coffee berry borer (Hypotheneumus hampei), is responsible for worldwide annual losses of around US500 million. The coffee berry borer exclusively damages the coffee berries, and it is mainly controlled by organochlorine insecticides that are both toxic and carcinogenic. Unfortunately, natural resistance in the genus Coffea to H. hampei has not been documented. To overcome these problems, biotechnological strategies can be used to introduce an α-amylase inhibitor gene (α-AI1), which confers resistance against the coffee berry borer insect-pest, into C. arabica plants.     

CodonTest: Modeling Amino Acid Substitution Preferences in Coding Sequences

Codon models of evolution have facilitated the interpretation of selective forces operating on genomes. These models, however, assume a single rate of non-synonymous substitution irrespective of the nature of amino acids being exchanged. Recent developments have shown that models which allow for amino acid pairs to have independent rates of substitution offer improved fit over single rate models. However, these approaches have been limited by the necessity for large alignments in their estimation. An alternative approach is to assume that substitution rates between amino acid pairs can be subdivided into rate classes, dependent on the information content of the alignment. However, given the combinatorially large number of such models, an efficient model search strategy is needed. Here we develop a Genetic Algorithm (GA) method for the estimation of such models. A GA is used to assign amino acid substitution pairs to a series of rate classes, where is estimated from the alignment. Other parameters of the phylogenetic Markov model, including substitution rates, character frequencies and branch lengths are estimated using standard maximum likelihood optimization procedures. We apply the GA to empirical alignments and show improved model fit over existing models of codon evolution. Our results suggest that current models are poor approximations of protein evolution and thus gene and organism specific multi-rate models that incorporate amino acid substitution biases are preferred. We further anticipate that the clustering of amino acid substitution rates into classes will be biologically informative, such that genes with similar functions exhibit similar clustering, and hence this clustering will be useful for the evolutionary fingerprinting of genes.     

Scientific Authority in the Creation–Evolution Debates

At the heart of debates among creationists and evolutionists are questions about scientific integrity and rigor. Creationists often justify their rejection of biological evolution by claiming that the methodologies and interpretations of evolutionary scientists are flawed. A consideration of creationists’ critiques of the scientific data, however, reveals a deficient understanding and appreciation of the nature of the scientific process. It is essential that our schools educate students about the character of scientific inquiry. Clarifying the nature and limitation of scientific knowledge for our students will equip our students to evaluate evolutionary or creationist arguments critically. Recognizing and teaching both the strengths and limitations of the scientific process will do much to further the ongoing dialogue between science and religion.