Adaptive protein evolution and regulatory divergence in Drosophila

Two recent studies demonstrated a positive correlation between divergence in gene expression and protein sequence in Drosophila. This correlation could be driven by positive selection or variation in functional constraint. To distinguish between these alternatives, we compared patterns of molecular evolution for 1,862 genes with two previously reported estimates of expression divergence in Drosophila. We found a slight negative trend (nonsignificant) between positive selection on protein sequence and divergence in expression levels between Drosophila melanogaster and Drosophila simulans. Conversely, shifts in expression patterns during Drosophila development showed a positive association with adaptive protein evolution, though as before the relationship was weak and not significant. Overall, we found no strong evidence for an increase in the incidence of positive selection on protein-coding regions in genes with divergent expression in Drosophila, suggesting that the previously reported positive association between protein and regulatory divergence primarily reflects variation in functional constraint.     

Relation between protein stability, evolution and structure, as probed by carboxylic acid mutations

Native proteins are marginally stable. Low thermodynamic stability may actually be advantageous, although the accumulation of neutral, destabilizing mutations may have also contributed to it. In any case, once marginal stability has been reached, it appears plausible that mutations at non-constrained positions become fixed in the course of evolution (due to random drift) with frequencies that roughly reflect the mutation effects on stability ("pseudo-equilibrium hypothesis"). We have found that all glutamate-->aspartate mutations in wild-type Escherichia coli thioredoxin are destabilizing, as well as most of the aspartate-->glutamate mutations. Furthermore, the effect of these mutations on thioredoxin thermodynamic stability shows a robust correlation with the frequencies of occurrence of the involved residues in several-hundred sequence alignments derived from a BLAST search. These results provide direct and quantitative experimental evidence for the pseudo-equilibrium hypothesis and should have general consequences for the interpretation of mutation effects on protein stability, as they suggest that residue environments in proteins may be optimized for stabilizing interactions to a remarkable degree of specificity. We also provide evidence that such stabilizing interactions may be detected in sequence alignments, and briefly discuss the implications of this possibility for the derivation of structural information (on native and denatured states) from comparative sequence analyses.     

Rates of phenotypic evolution of ecological characters and sexual traits during the Tanganyikan cichlid adaptive radiation

Theory suggests that sexual traits evolve faster than ecological characters. However, characteristics of a species niche may also influence evolution of sexual traits. Hence, a pending question is whether ecological characters and sexual traits present similar tempo and mode of evolution during periods of rapid ecological divergence, such as adaptive radiation. Here, we use recently developed phylogenetic comparative methods to analyse the temporal dynamics of evolution for ecological and sexual traits in Tanganyikan cichlids. Our results indicate that whereas disparity in ecological characters was concentrated early in the radiation, disparity in sexual traits remained high throughout the radiation. Thus, closely related Tanganyikan cichlids presented higher disparity in sexual traits than ecological characters. Sexual traits were also under stronger selection than ecological characters. In sum, our results suggest that ecological characters and sexual traits present distinct evolutionary patterns, and that sexual traits can evolve faster than ecological characters, even during adaptive radiation.     

Insights from modeling protein evolution with context-dependent mutation and asymmetric amino acid selection

We develop an approximate maximum likelihood method to estimate flanking nucleotide context-dependent mutation rates and amino acid exchange-dependent selection in orthologous protein-coding sequences and use it to analyze genome-wide coding sequence alignments from mammals and yeast. Allowing context-dependent mutation provides a better fit to coding sequence data than simpler (context-independent or CpG "hotspot") models and significantly affects selection parameter estimates. Allowing asymmetric (nonreciprocal) selection on amino acid exchanges gives a better fit than simple dN/dS or symmetric selection models. Relative selection strength estimates from our models show good agreement with independent estimates derived from human disease-causing and engineered mutations. Selection strengths depend on local protein structure, showing expected biophysical trends in helical versus nonhelical regions and increased asymmetry on polar-hydrophobic exchanges with increased burial. The more stringent selection that has previously been observed for highly expressed proteins is primarily concentrated in buried regions, supporting the notion that such proteins are under stronger than average selection for stability. Our analyses indicate that a highly parameterized model of mutation and selection is computationally tractable and is a useful tool for exploring a variety of biological questions concerning protein and coding sequence evolution.     

Scalable method to determine mutations that occur during adaptive evolution of Escherichia coli

Denaturing HPLC was used to determine mutations occurring during the adaptive evolution of Escherichia coli K-12. The strains were evolved over 700 generations on glycerol as the sole carbon source from a sub-optimal to an optimal growth rate. The mutations detected by direct sequencing of amplicons of the glycerol-phosphate regulon repressor (glpR) gene were a synonymous substitution Val20Val in two separately evolved strains. Non-synonymous substitutions, Val119Gly and Gly179Trp, were also observed in each of the two strains. This procedure can be scaled to determine genome-scale sequence variations that have occurred during adaptive evolution.     

Genome architecture drives protein evolution in ciliates

Studies of microbial eukaryotes have been pivotal in the discovery of biological phenomena, including RNA editing, self-splicing RNA, and telomere addition. Here we extend this list by demonstrating that genome architecture, namely the extensive processing of somatic (macronuclear) genomes in some ciliate lineages, is associated with elevated rates of protein evolution. Using newly developed likelihood-based procedures for studying molecular evolution, we investigate 6 genes to compare 1) ciliate protein evolution to that of 3 other clades of eukaryotes (plants, animals, and fungi) and 2) protein evolution in ciliates with extensively processed macronuclear genomes to that of other ciliate lineages. In 5 of the 6 genes, ciliates are estimated to have a higher ratio of nonsynonymous/synonymous substitution rates, consistent with an increase in the rate of protein diversification in ciliates relative to other eukaryotes. Even more striking, there is a significant effect of genome architecture within ciliates as the most divergent proteins are consistently found in those lineages with the most highly processed macronuclear genomes. We propose a model whereby genome architecture-specifically chromosomal processing, amitosis within macronuclei, and epigenetics-allows ciliates to explore protein space in a novel manner. Further, we predict that examination of diverse eukaryotes will reveal additional evidence of the impact of genome architecture on molecular evolution.     

Structure and evolution of the human genes encoding protein C and coagulation factor IX

Human protein C is a vitamin K-dependent plasma protein that serves as a feedback down-regulator of the coagulation cascade by specifically degrading the protein cofactors VIIIa and Va. The protein C precursor consists of the following domains: leader peptide, "gla" region, two epidermal growth factor segments, and the activation peptide/serine protease. Comparison of amino acid sequences reveals that protein C and factor IX are homologous. A comparison of the genes for protein C and factor IX shows that all seven of the introns within the protein coding regions are in identical positions and correspond to protein structure-function domain boundries. However, the base compositions of the two genes (coding and noncoding regions) are remarkably different: approximately 60% guanine + cytosine (G + C) for protein C versus approximately 40% G + C for factor IX. One possible explanation for this phenomenon is that the factor IX gene (located on the X chromosome) has undergone extensive deoxycytosine methylation and subsequent spontaneous deamination mutagenesis, resulting in a net C to thymine (and G to adenine) transition. This would suggest that the protein C gene may represent a more primitive form of the gene duplication precursor.     

Molecular evolution and population genetics of duplicated accessory gland protein genes in Drosophila

To investigate the potential importance of gene duplication in D. melanogaster accessory gland protein (Acp) gene evolution we carried out a computational analysis comparing annotated D. melanogaster Acp genes to the entire D. melanogaster genome. We found that two known Acp genes are actually members of small multigene families. Polymorphism and divergence data from these duplicated genes suggest that in at least four cases, protein divergence between D. melanogaster and D. simulans is a result of directional selection. One putative Acp revealed by our computational analysis shows evidence of a recent selective sweep in a non-African population (but not in an African population). These data support the idea that selection on reproduction-related genes may drive divergence of populations within species, and strengthen the conclusion that Acps may often be under directional selection in Drosophila.     

Effects of metabolic rate on protein evolution

Since the modern evolutionary synthesis was first proposed early in the twentieth century, attention has focused on assessing the relative contribution of mutation versus natural selection on protein evolution. Here we test a model that yields general quantitative predictions on rates of protein evolution by combining principles of individual energetics with Kimura's neutral theory. The model successfully predicts much of the heterogeneity in rates of protein evolution for diverse eukaryotes (i.e. fishes, amphibians, reptiles, birds, mammals) from different thermal environments. Data also show that the ratio of non-synonymous to synonymous nucleotide substitution is independent of body size, and thus presumably of effective population size. These findings indicate that rates of protein evolution are largely controlled by mutation rates, which in turn are strongly influenced by individual metabolic rate.     

Domain tree-based analysis of protein architecture evolution

Understanding the dynamics behind domain architecture evolution is of great importance to unravel the functions of proteins. Complex architectures have been created throughout evolution by rearrangement and duplication events. An interesting question is how many times a particular architecture has been created, a form of convergent evolution or domain architecture reinvention. Previous studies have approached this issue by comparing architectures found in different species. We wanted to achieve a finer-grained analysis by reconstructing protein architectures on complete domain trees. The prevalence of domain architecture reinvention in 96 genomes was investigated with a novel domain tree-based method that uses maximum parsimony for inferring ancestral protein architectures. Domain architectures were taken from Pfam. To ensure robustness, we applied the method to bootstrap trees and only considered results with strong statistical support. We detected multiple origins for 12.4% of the scored architectures. In a much smaller data set, the subset of completely domain-assigned proteins, the figure was 5.6%. These results indicate that domain architecture reinvention is a much more common phenomenon than previously thought. We also determined which domains are most frequent in multiply created architectures and assessed whether specific functions could be attributed to them. However, no strong functional bias was found in architectures with multiple origins.     

Domain rearrangements in protein evolution

Most eukaryotic proteins are multi-domain proteins that are created from fusions of genes, deletions and internal repetitions. An investigation of such evolutionary events requires a method to find the domain architecture from which each protein originates. Therefore, we defined a novel measure, domain distance, which is calculated as the number of domains that differ between two domain architectures. Using this measure the evolutionary events that distinguish a protein from its closest ancestor have been studied and it was found that indels are more common than internal repetition and that the exchange of a domain is rare. Indels and repetitions are common at both the N and C-terminals while they are rare between domains. The evolution of the majority of multi-domain proteins can be explained by the stepwise insertions of single domains, with the exception of repeats that sometimes are duplicated several domains in tandem. We show that domain distances agree with sequence similarity and semantic similarity based on gene ontology annotations. In addition, we demonstrate the use of the domain distance measure to build evolutionary trees. Finally, the evolution of multi-domain proteins is exemplified by a closer study of the evolution of two protein families, non-receptor tyrosine kinases and RhoGEFs.     

The distribution of fitness effects of new beneficial mutations in Pseudomonas fluorescens

Theoretical studies of adaptation emphasize the importance of understanding the distribution of fitness effects (DFE) of new mutations. We report the isolation of 100 adaptive mutants—without the biasing influence of natural selection—from an ancestral genotype whose fitness in the niche occupied by the derived type is extremely low. The fitness of each derived genotype was determined relative to a single reference type and the fitness effects found to conform to a normal distribution. When fitness was measured in a different environment, the rank order changed, but not the shape of the distribution. We argue that, even with detailed knowledge of the genetic architecture underpinning the adaptive types (as is the case here), the DFEs remain unpredictable, and we discuss the possibility that general explanations for the shape of the DFE might not be possible in the absence of organism-specific biological details.     

A quasispecies approach to viral evolution in the context of an adaptive immune system

A deeper understanding of the mechanisms that determine viral evolution in the context of an adaptive immune system is vital for the development of efficient strategies to defeat viral infections. The problem of describing these mechanisms is discussed using the concept of quasispecies. Conditions for both an optimal immune response and for highest viral viability are derived from theoretical models and are supported by empirical data.     

Gradual adaptive changes of a protein facing high salt concentrations

Several experimental techniques were applied to unravel fine molecular details of protein adaptation to high salinity. We compared four homologous enzymes, which suggested a new halo-adaptive state in the process of molecular adaptation to high-salt conditions. Together with comparative functional studies, the structure of malate dehydrogenase from the eubacterium Salinibacter ruber shows that the enzyme shares characteristics of a halo-adapted archaea-bacterial enzyme and of non-halo-adapted enzymes from other eubacterial species. The S. ruber enzyme is active at the high physiological concentrations of KCl but, unlike typical halo-adapted enzymes, remains folded and active at low salt concentrations. Structural aspects of the protein, including acidic residues at the surface, solvent-exposed hydrophobic surface, and buried hydrophobic surface, place it between the typical halo-adapted and non-halo-adapted proteins. The enzyme lacks inter-subunit ion-binding sites often seen in halo-adapted enzymes. These observations permit us to suggest an evolutionary pathway that is highlighted by subtle trade-offs to achieve an optimal compromise among solubility, stability, and catalytic activity.     

Life-history evolution at the molecular level: adaptive amino acid composition of avian vitellogenins

Avian genomes typically encode three distinct vitellogenin (VTG) egg yolk proteins (VTG1, VTG2 and VTG3), which arose by gene duplication prior to the most recent common ancestor of birds. Analysis of VTG sequences from 34 avian species in a phylogenetic framework supported the hypothesis that VTG amino acid composition has co-evolved with embryo incubation time. Embryo incubation time was positively correlated with the proportions of dietary essential amino acids (EAAs) in VTG1 and VTG2, and with the proportion of sulfur-containing amino acids in VTG3. These patterns were seen even when only semi-altricial and/or altricial species were considered, suggesting that the duration of embryo incubation is a major selective factor on the amino acid composition of VTGs, rather than developmental mode alone. The results are consistent with the hypothesis that the level of EAAs provided to the egg represents an adaptation to the loss of amino acids through breakdown over the course of incubation and imply that life-history phenotypes and VTG amino acid composition have co-evolved throughout the evolutionary history of birds.     

Epistasis and frequency dependence influence the fitness of an adaptive mutation in a diversifying lineage

The opportunity for a mutation to invade a population can dramatically vary depending on the context in which this mutation occurs. Such context dependence is difficult to document as it requires the ability to measure how a mutation affects phenotypes and fitness and to manipulate the context in which the mutation occurs. We identified a mutation in a gene encoding a global regulator in one of two ecotypes that diverged from a common ancestor during 1200 generations of experimental evolution. We replaced the ancestral allele by the mutant allele, and vice versa, in several clones isolated during the time course of the evolution experiment, and compared the phenotype and fitness of clones isogenic except for the focal mutation. We show that the fitness and phenotype of the mutation are strongly affected by epistatic interactions between genes in the same genome, as well as by frequency dependent selection resulting from biotic interactions between individuals in the same population. We conclude that amongst the replicate population in which it spread, the mutation we identified is only adaptive when occurring in specific genomes and competing with specific individuals. This study thus demonstrates that the opportunity for an adaptive mutation to spread in an evolutionary lineage can only be understood in the light of its genomic and competitive environments.     

Functional divergence in protein (family) sequence evolution

As widely used today to infer function, the homology search is based on the neutral theory that sites of greatest functional significance are under the strongest selective constraints as well as lowest evolutionary rates, and vice versa. Therefore, site-specific rate changes (or altered selective constraints) are related to functional divergence during protein (family) evolution. In this paper, we review our recent work about this issue. We show a great deal of functional information can be obtained from the evolutionary perspective, which can in turn be used to facilitate high throughput functional assays. The emergence of evolutionary functional genomics is also indicated. The related software DIVERGE can be obtained form http://xgu1.zool.iastate.edu.     

Detecting lineage-specific adaptive evolution of brain-expressed genes in human using rhesus macaque as outgroup

Comparative genetic analysis between human and chimpanzee may detect genetic divergences responsible for human-specific characteristics. Previous studies have identified a series of genes that potentially underwent Darwinian positive selection during human evolution. However, without a closely related species as outgroup, it is difficult to identify human-lineage-specific changes, which is critical in delineating the biological uniqueness of humans. In this study, we conducted phylogeny-based analyses of 2633 human brain-expressed genes using rhesus macaque as the outgroup. We identified 47 candidate genes showing strong evidence of positive selection in the human lineage. Genes with maximal expression in the brain showed a higher evolutionary rate in human than in chimpanzee. We observed that many immune-defense-related genes were under strong positive selection, and this trend was more prominent in chimpanzee than in human. We also demonstrated that rhesus macaque performed much better than mouse as an outgroup in identifying lineage-specific selection in humans.     

Molecular evolution of the mammalian prion protein

Prion protein (PrP) sequences are until now available for only six of the 18 orders of placental mammals. A broader comparison of mammalian prions might help to understand the enigmatic functional and pathogenic properties of this protein. We therefore determined PrP coding sequences in 26 mammalian species to include all placental orders and major subordinal groups. Glycosylation sites, cysteines forming a disulfide bridge, and a hydrophobic transmembrane region are perfectly conserved. Also, the sequences responsible for secondary structure elements, for N- and C-terminal processing of the precursor protein, and for attachment of the glycosyl-phosphatidylinositol membrane anchor are well conserved. The N-terminal region of PrP generally contains five or six repeats of the sequence P(Q/H)GGG(G/-)WGQ, but alleles with two, four, and seven repeats were observed in some species. This suggests, together with the pattern of amino acid replacements in these repeats, the regular occurrence of repeat expansion and contraction. Histidines implicated in copper ion binding and a proline involved in 4-hydroxylation are lacking in some species, which questions their importance for normal functioning of cellular PrP. The finding in certain species of two or seven repeats, and of amino acid substitutions that have been related to human prion diseases, challenges the relevance of such mutations for prion pathology. The gene tree deduced from the PrP sequences largely agrees with the species tree, indicating that no major deviations occurred in the evolution of the prion gene in different placental lineages. In one species, the anteater, a prion pseudogene was present in addition to the active gene.