Development of an unbiased statistical method for the analysis of unigenic evolution

Background  

Unigenic evolution is a powerful genetic strategy involving random mutagenesis of a single gene product to delineate functionally important domains of a protein. This method involves selection of variants of the protein which retain function, followed by statistical analysis comparing expected and observed mutation frequencies of each residue. Resultant mutability indices for each residue are averaged across a specified window of codons to identify hypomutable regions of the protein. As originally described, the effect of changes to the length of this averaging window was not fully eludicated. In addition, it was unclear when sufficient functional variants had been examined to conclude that residues conserved in all variants have important functional roles.     

Functionally important residues in the peptidyl-prolyl isomerase Pin1 revealed by unigenic evolution

Pin1 is a phosphorylation-dependent member of the parvulin family of peptidyl-prolyl isomerases exhibiting functional conservation between yeast and man. To perform an unbiased analysis of the regions of Pin1 essential for its functions, we generated libraries of randomly mutated forms of the human Pin1 cDNA and identified functional Pin1 alleles by their ability to complement the Pin1 homolog Ess1 in Saccharomyces cerevisiae. We isolated an extensive collection of functional mutant Pin1 clones harboring a total of 356 amino acid substitutions. Surprisingly, many residues previously thought to be critical in Pin1 were found to be altered in this collection of functional mutants. In fact, only 17 residues were completely conserved in these mutants and in Pin1 sequences from other eukaryotic organisms, with only two of these conserved residues located within the WW domain of Pin1. Examination of invariant residues provided new insights regarding a phosphate-binding loop that distinguishes a phosphorylation-dependent peptidyl-prolyl isomerase such as Pin1 from other parvulins. In addition, these studies led to an investigation of residues involved in catalysis including C113 that was previously implicated as the catalytic nucleophile. We demonstrate that substitution of C113 with D does not compromise Pin1 function in vivo nor does this substitution abolish catalytic activity in purified recombinant Pin1. These findings are consistent with the prospect that the function of residue 113 may not be that of a nucleophile, thus raising questions about the model of nucleophilic catalysis. Accordingly, an alternative catalytic mechanism for Pin1 is postulated.     

The evolution of senescence through decelerating selection for system reliability

Senescence is a universal phenomenon in organisms, characterized by increasing mortality and decreasing fecundity with advancing chronological age. Most proximate agents of senescence, such as reactive oxygen species and UV radiation, are thought to operate by causing a gradual build-up of bodily damage. Yet most current evolutionary theories of senescence emphasize the deleterious effects of functioning genes in late life, leaving a gap between proximate and ultimate explanations. Here, we present an evolutionary model of senescence based on reliability theory, in which beneficial genes or gene products gradually get damaged and thereby fail, rather than actively cause harm. Specifically, the model allows organisms to evolve multiple redundant copies of a gene product (or gene) that performs a vital function, assuming that organisms can avoid condition-dependent death so long as at least one copy remains undamaged. We show that organisms with low levels of extrinsic mortality, and high levels of genetic damage, tend to evolve high levels of redundancy, and that mutation-selection balance results in a stable population distribution of the number of redundant elements. In contrast to previous evolutionary models of senescence, the mortality curves that emerge from such populations match empirical senescence patterns in three key respects: they exhibit: (1) an initially low, but rapidly increasing mortality rate at young ages, (2) a plateau in mortality at advanced ages and (3) 'mortality compensation', whereby the height of the mortality plateau is independent of the environmental conditions under which different populations evolved.     

Molecular evolution of transferrin: evidence for positive selection in salmonids

Transferrins are iron-binding proteins that are involved in iron storage and resistance to bacterial disease. Previous work has shown that nonsynonymous-to-synonymous-site substitution ratios (d(n)/d(s) ratios) between transferrin genes from some salmonid species were significantly greater than 1.0, providing evidence for positive selection at the transferrin gene. The purpose of the current study was to put these earlier results in a broader evolutionary context by examining variation among 25 previously published transferrin sequences from fish, amphibians, and mammals. The results of the study show that evidence for positive selection at transferrin is limited to salmonids-d(n)/d(s) ratios estimated for nonsalmonid lineages were generally less than 1.0. Within the salmonids, approximately 13% of the transferrin codon sites are estimated to be subject to positive selection, with an estimated d(n)/d(s) ratio of approximately 7. The three- dimensional locations of some of the selected sites were inferred by comparing these sites to homologous sites in the bovine lactoferrin crystallographic structure. The selected sites generally fall on the outside of the molecule, within and near areas that are bound by transferrin-binding proteins from human pathogenic bacteria. The physical locations of sites estimated to be subject to positive selection support previous speculation that competition for iron from pathogenic bacteria could be the source of positive selection.     

Estimating the reliability of evolutionary trees

Six protein sequences from the same 11 mammalian taxa were used to estimate the accuracy and reliability of phylogenetic trees using real, rather than simulated, data. A tree comparison metric was used to measure the increase in similarity of minimal trees as larger, randomly selected subsets of nucleotide positions were taken. The ratio of the observed to the expected number of incompatibilities for each nucleotide position (character) is a good predictor of the number of changes required at that position on the minimal (most-parsimonious) tree. This allows a higher weighting of nucleotide positions that have changed more slowly and should result in the minimal length tree converging to the correct tree as more sequences are obtained. An estimate was made of the smallest subset of trees that need to be considered to include the actual historical tree for a given set of data. It was concluded that it is possible to give a reasonable estimate of the reliability of the final tree, at least when several sequences are combined. With the present data, resolving the rodent- primate-lagomorph (rabbit) trichotomy is the least certain aspect of the final tree, followed then by establishing the position of dog. In our opinion, it is unreasonable to publish an evolutionary tree derived from sequence data without giving an idea of the reliability of the tree.      

The evolution of the phenotypic covariance matrix: evidence for selection and drift in Melanoplus

Phenotypic variation in trait means is a common observation for geographically separated populations. Such variation is typically retained under common garden conditions, indicating that there has been evolutionary change in the populations, as a result of selection and/or drift. Much less frequently studied is variation in the phenotypic covariance matrix (hereafter, P matrix), although this is an important component of evolutionary change. In this paper, we examine variation in the phenotypic means and P matrices in two species of grasshopper, Melanoplus sanguinipes and M. devastator. Using the P matrices estimated for 14 populations of M. sanguinipes and three populations of M. devastator we find that (1) significant differences between the sexes can be attributed to scaling effects; (2) there is no significant difference between the two species; (3) there are highly significant differences among populations that cannot be accounted for by scaling effects; (4) these differences are a consequence of statistically significant patterns of covariation with geographic and environmental factors, phenotypic variances and covariances increasing with increased temperature but decreasing with increased latitude and altitude. This covariation suggests that selection has been important in the evolution of the P matrix in these populations Finally, we find a significant positive correlation between the average difference between matrices and the genetic distance between the populations, indicating that drift has caused some of the variation in the P matrices.     

Individual selection, kin selection, and the shifting balance in the evolution of warning colours: the evidence from butterflies

It is difficult to imagine how warning colours evolve in unpalatable prey. Firstly, novel warningly coloured variants gain no protection from their colours, since predators have not previously encountered and learnt their colour patterns. This leads to a frequency-dependent disadvantage of a rare variant within a species. Secondly, novel warningly coloured variants may be more conspicuous than non-aposematic prey.
Nevertheless, it is obvious that many palatable butterflies have bright colours used in intraspecific communication and in duping predators. Other palatable butterflies are already warningly coloured. Should such butterflies evolve unpalatability, perhaps because of a host-plant shift, these bright colours would be preadapted to a warning role. Warning colours could then continue to evolve by enhancement of memorable characteristics of these patterns, or by mimicry.
Even within lineages of warningly coloured, unpalatable butterflies, colour patterns have continued to evolve rapidly. This diversity of warning colour patterns could have evolved in a number of ways, including individual and kin selection, and by the shifting balance. Evidence for these mechanisms is discussed, as are the similarities between the evolution of warning colours and more general evolutionary processes, including sexual selection and speciation.     

Natural selection and the evolution of mtDNA-encoded peptides: evidence for intergenomic co-adaptation.

Mitochondrial DNA (mtDNA) variation is an important tool for the investigation of the population genetics of animal species. Recently, recognition of the role of mtDNA mutations in human disease has spurred increasing interest in the function and evolution of mtDNA and the 13 polypeptides it encodes. These proteins interact with a large number of peptides encoded in the nucleus to form the mitochondrial electron transport system (ETS). As the ETS is the primary energy generation system in aerobic metazoans, natural selection would be expected to favor mutations that enhance ETS function. Such mutations could occur in either the mitochondrial or nuclear genes encoding ETS proteins and would lead to positive intergenomic interactions, or co-adaptation. Direct evidence for intergenomic co-adaptation comes from functional studies of systems where nuclear-mitochondrial DNA combinations vary naturally or can be manipulated experimentally.     

Estimating the reliability of diet overlap measures

Synopsis Diet overlap measures, commonly used in studies of resource partitioning and competition in fish, are too often treated as fixed values; in fact, they are random variables. Two methods for estimating the variance of some overlap measures using stomach content data are examined here: the jackknife and the bootstrap. Simulation results indicate that the methods work well. In addition, they are valid for either continuous (biomass, percent volume) or discrete (abundance) data. Applications indicate that outliers, such as a stomach much different than others, may have large influence in determing the overlap value.     

Molecular evolution of rickettsia surface antigens: evidence of positive selection

The Rickettsia genus is a group of obligate intracellular parasitic alpha-proteobacteria that includes human pathogens responsible for the typhus disease and various types of spotted fevers. rOmpA and rOmpB are two members of the "surface cell antigen" (Sca) autotransporter (AT) protein family that may play key roles in the adhesion of the Rickettsia cells to the host tissue. These molecules are likely determinants for the pathogenicity of the Rickettsia and represent good candidates for vaccine development. We identified the 17 members of this family of outer-membrane proteins in nine fully sequenced Rickettsia genomes. The typical architecture of the Sca proteins is composed of an N-terminal signal peptide and a C-terminal AT domain that promote the export of the central passenger domain to the outside of the bacteria. A characteristic of this family is the frequent degradation of the genes, which results in different subsets of the sca genes being expressed among Rickettsia species. Here, we present a detailed analysis of their phylogenetic relationships and evolution. We provide strong evidence that rOmpA and rOmpB as well as three other members of the Sca protein family--Sca1, Sca2, and Sca4--have evolved under positive selection. The exclusive distribution of the predicted positively selected sites within the passenger domains of these proteins argues that these regions are involved in the interaction with the host and may be locked in "arms race" coevolutionary conflicts.     

Bayesian inference of character evolution

Much recent progress in evolutionary biology is based on the inference of ancestral states and past transformations in important traits on phylogenetic trees. These exercises often assume that the tree is known without error and that ancestral states and character change can be mapped onto it exactly. In reality, there is often considerable uncertainty about both the tree and the character mapping. Recently introduced Bayesian statistical methods enable the study of character evolution while simultaneously accounting for both phylogenetic and mapping uncertainty, adding much needed credibility to the reconstruction of evolutionary history.     

Estimating the rate of evolution of the rate of molecular evolution

A simple model for the evolution of the rate of molecular evolution is presented. With a Bayesian approach, this model can serve as the basis for estimating dates of important evolutionary events even in the absence of the assumption of constant rates among evolutionary lineages. The method can be used in conjunction with any of the widely used models for nucleotide substitution or amino acid replacement. It is illustrated by analyzing a data set of rbcL protein sequences.      

Estimating fluctuating selection in age-structured populations

In age-structured populations, viability and fecundity selection of varying strength may occur in different age classes. On the basis of an original idea by Fisher of weighting individuals by their reproductive value, we show that the combined effect of selection on traits at different ages acts through the individual reproductive value defined as the stochastic contribution of an individual to the total reproductive value of the population the following year. The selection differential is a weighted sum of age-specific differentials that are the covariances between the phenotype and the age-specific relative fitness defined by the individual reproductive value. This enables estimation of weak selection on a multivariate quantitative character in populations with no density regulation by combinations of age-specific linear regressions of individual reproductive values on the traits. Demographic stochasticity produces random variation in fitness components in finite samples of individuals and affects the statistical inference of the temporal average directional selection as well as the magnitude of fluctuating selection. Uncertainties in parameter estimates and test power depend strongly on the demographic stochasticity. Large demographic variance results in large uncertainties in yearly estimates of selection that complicates detection of significant fluctuating selection. The method is illustrated by an analysis of age-specific selection in house sparrows on a fitness-related two-dimensional morphological trait, tarsus length and body mass of fledglings.     

The evolution of base composition and phylogenetic inference

Base composition varies at all levels of the phylogenetic hierarchy and throughout the genome, and can be caused by active selection or passive mutation pressure. This variation can make reconstructing trees difficult. However, recent tree-based analyses have shed light on the forces responsible for the evolution of base composition, forces that might be very general. More explicit tree-based work is encouraged.     

Fertility selection,genetic selection and evolution

The relationship between fertility selection as measured by the correlation in progeny number between parents and offspring, and selection at individual loci is investigated in humans. Estimates for the magnitude of fertility selection (0.1) and the rate of gene substitution (0.5 gene substitutions per generation per genome) are used in various mathematical models for selection. It is found that the observed magnitude of fertility selection cannot be explained by non‐epistatic directional selection at individual loci. A symmetric quantitative directional selection model is consistent with the observed data. But it is possible that fertility selection does not have a genetic basis.     

Kin selection and the evolution of virulence

Social interactions between conspecific parasites are partly dependent on the relatedness of interacting parasites (kin selection), which, in turn, is predicted to affect the extent of damage they cause their hosts (virulence). High relatedness is generally assumed to favour less competitive interactions, but the relationship between relatedness and virulence is crucially dependent on the social behaviour in question. Here, we discuss the rather limited body of experimental work that addresses how kin-selected social behaviours affect virulence. First, if prudent use of host resources (a form of cooperation) maximizes the transmission success of the parasite population, decreased relatedness is predicted to result in increased host exploitation and virulence. Experimental support for this well-established theoretical result is surprisingly limited. Second, if parasite within-host growth rate is a positive function of cooperation (that is, when individuals need to donate public goods, such as extracellular enzymes), virulence is predicted to increase with increasing relatedness. The limited studies testing this hypothesis are broadly consistent with this prediction. Finally, there is some empirical evidence supporting theory that suggests that spiteful behaviours are maximized at intermediate degrees of relatedness, which, in turn, leads to minimal virulence because of the reduced growth rate of the infecting population. We highlight the need for further thorough experimentation on the role of kin selection in the evolution of virulence and identify additional biological complexities to these simple frameworks.     

Positive selection in the evolution of cancer

We hypothesize that forms of antagonistic coevolution have forged strong links between positive selection at the molecular level and increased cancer risk. By this hypothesis, evolutionary conflict between males and females, mothers and foetuses, hosts and parasites, and other parties with divergent fitness interests has led to rapid evolution of genetic systems involved in control over fertilization and cellular resources. The genes involved in such systems promote cancer risk as a secondary effect of their roles in antagonistic coevolution, which generates evolutionary disequilibrium and maladaptation. Evidence from two sources: (1) studies on specific genes, including SPANX cancer/testis antigen genes, several Y-linked genes, the pem homebox gene, centromeric histone genes, the breast cancer gene BRCA1, the angiogenesis gene ANG, cadherin genes, cytochrome P450 genes, and viral oncogenes; and (2) large-scale database studies of selection on different functional categories of genes, supports our hypothesis. These results have important implications for understanding the evolutionary underpinnings of cancer and the dynamics of antagonistically-coevolving molecular systems.     

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.