The use of evolutionary patterns in protein annotation

With genomic data skyrocketing, their biological interpretation remains a serious challenge. Diverse computational methods address this problem by pointing to the existence of recurrent patterns among sequence, structure, and function. These patterns emerge naturally from evolutionary variation, natural selection, and divergence--the defining features of biological systems--and they identify molecular events and shapes that underlie specificity of function and allosteric communication. Here we review these methods, and the patterns they identify in case studies and in proteome-wide applications, to infer and rationally redesign function.     

Can Evolution Explain Insanity?

I distinguish three evolutionary explanations of mental illness: first, breakdowns in evolved computational systems; second, evolved systems performing their evolutionary function in a novel environment; third, evolved personality structures. I concentrate on the second and third explanations, as these are distinctive of an evolutionary psychopathology, with progressively less credulity in the light of the empirical evidence. General morals are drawn for evolutionary psychiatry.     

Phylogeny as a guide to structure and function of membrane transport proteins

Protein phylogeny, based on primary amino acid sequence relatedness, reflects the evolutionary process and therefore provides a guide to structure, mechanism and function. Any two proteins that are related by common descent are expected to exhibit similar structures and functions to a degree proportional to the degree of their sequence similarity; but two independently evolving proteins should not. This principle provides the impetus to define protein phylogenetic relationships and interrelate families when possible. In this mini-review, we summarize the computational approaches and criteria we use to establish common evolutionary origin. We apply these tools to define distant superfamily relationships between several previously recognized transport protein families. In some cases, available structural and functional data are evaluated in order to substantiate our claim that molecular phylogeny provides a reliable guide to protein structure and function.     

Stabilizing selection of protein function and distribution of selection coefficient among sites

In this study, I take a new approach to modeling the evolutionary constraint of protein sequence, introducing the stabilizing selection of protein function into the nearly-neutral theory. In other words, protein function under stabilizing selection generates the evolutionary conservation at the sequence level. With the help of random mutational effects of nucleotides on protein function, I have derived the distribution of selection coefficient among sites, called the S-distribution whose parameters have clear biological interpretations. Moreover, I have studied the inverse relationship between the evolutionary rate and the effective population size, showing that the number of molecular phenotypes of protein function, i.e., independent components in the fitness of the organism, may play a key role for the molecular clock under the nearly-neutral theory. These results are helpful for having a better understanding of the underlying evolutionary mechanism of protein sequences, as well as human disease-related mutations.     

Adaptive evolution and functional redesign of core metabolic proteins in snakes

Background

Adaptive evolutionary episodes in core metabolic proteins are uncommon, and are even more rarely linked to major macroevolutionary shifts.

Methodology/Principal Findings

We conducted extensive molecular evolutionary analyses on snake mitochondrial proteins and discovered multiple lines of evidence suggesting that the proteins at the core of aerobic metabolism in snakes have undergone remarkably large episodic bursts of adaptive change. We show that snake mitochondrial proteins experienced unprecedented levels of positive selection, coevolution, convergence, and reversion at functionally critical residues. We examined Cytochrome C oxidase subunit I (COI) in detail, and show that it experienced extensive modification of normally conserved residues involved in proton transport and delivery of electrons and oxygen. Thus, adaptive changes likely altered the flow of protons and other aspects of function in CO, thereby influencing fundamental characteristics of aerobic metabolism. We refer to these processes as “evolutionary redesign” because of the magnitude of the episodic bursts and the degree to which they affected core functional residues.

Conclusions/Significance

The evolutionary redesign of snake COI coincided with adaptive bursts in other mitochondrial proteins and substantial changes in mitochondrial genome structure. It also generally coincided with or preceded major shifts in ecological niche and the evolution of extensive physiological adaptations related to lung reduction, large prey consumption, and venom evolution. The parallel timing of these major evolutionary events suggests that evolutionary redesign of metabolic and mitochondrial function may be related to, or underlie, the extreme changes in physiological and metabolic efficiency, flexibility, and innovation observed in snake evolution.     

Modelling 'evo-devo' with RNA

The folding of RNA sequences into secondary structures is a simple yet biophysically grounded model of a genotype-phenotype map. Its computational and mathematical analysis has uncovered a surprisingly rich statistical structure characterized by shape space covering, neutral networks and plastogenetic congruence. I review these concepts and discuss their evolutionary implications.     

Experimental evolution, loss-of-function mutations, and "the first rule of adaptive evolution"

Adaptive evolution can cause a species to gain, lose, or modify a function; therefore, it is of basic interest to determine whether any of these modes dominates the evolutionary process under particular circumstances. Because mutation occurs at the molecular level, it is necessary to examine the molecular changes produced by the underlying mutation in order to assess whether a given adaptation is best considered as a gain, loss, or modification of function. Although that was once impossible, the advance of molecular biology in the past half century has made it feasible. In this paper, I review molecular changes underlying some adaptations, with a particular emphasis on evolutionary experiments with microbes conducted over the past four decades. I show that by far the most common adaptive changes seen in those examples are due to the loss or modification of a pre-existing molecular function, and I discuss the possible reasons for the prominence of such mutations.     

Bacterial mechanosensitive channels: Experiment and theory

Since their discovery in Escherichia coli some 20 years ago, studies of bacterial mechanosensitive (MS) ion channels have been at the forefront of the MS channel research field. Two major events greatly advanced the research on bacterial MS channels: (i) cloning of MscL and MscS, the MS channels of Large and Small conductance, and (ii) solving their 3D crystal structure. These events enabled further experimental studies employing EPR and FRET spectroscopy in addition to patch clamp and molecular biological techniques that have successfully been used in characterization of the structure and function of bacterial MS channels. In parallel with the experimental studies computational modelling has been applied to elucidate the molecular dynamics of MscL and MscS, which has significantly contributed to our understanding of basic physical principles of the mechanosensory transduction in living organisms.     

Analysis of the overdispersed clock in the short-term evolution of hepatitis C virus: Using the E1/E2 gene sequences to infer infection dates in a single source outbreak

The assumption of a molecular clock for dating events from sequence information is often frustrated by the presence of heterogeneity among evolutionary rates due, among other factors, to positively selected sites. In this work, our goal is to explore methods to estimate infection dates from sequence analysis. One such method, based on site stripping for clock detection, was proposed to unravel the clocklike molecular evolution in sequences showing high variability of evolutionary rates and in the presence of positive selection. Other alternatives imply accommodating heterogeneity in evolutionary rates at various levels, without eliminating any information from the data. Here we present the analysis of a data set of hepatitis C virus (HCV) sequences from 24 patients infected by a single individual with known dates of infection. We first used a simple criterion of relative substitution rate for site removal prior to a regression analysis. Time was regressed on maximum likelihood pairwise evolutionary distances between the sequences sampled from the source individual and infected patients. We show that it is indeed the fastest evolving sites that disturb the molecular clock and that these sites correspond to positively selected codons. The high computational efficiency of the regression analysis allowed us to compare the site-stripping scheme with random removal of sites. We demonstrate that removing the fast-evolving sites significantly increases the accuracy of estimation of infection times based on a single substitution rate. However, the time-of-infection estimations improved substantially when a more sophisticated and computationally demanding Bayesian method was used. This method was used with the same data set but keeping all the sequence positions in the analysis. Consequently, despite the distortion introduced by positive selection on evolutionary rates, it is possible to obtain quite accurate estimates of infection dates, a result of especial relevance for molecular epidemiology studies.     

Effects of Mitochondrial DNA Rate Variation on Reconstruction of Pleistocene Demographic History in a Social Avian Species,Pomatostomus superciliosus

Mitochondrial sequence data is often used to reconstruct the demographic history of Pleistocene populations in an effort to understand how species have responded to past climate change events. However, departures from neutral equilibrium conditions can confound evolutionary inference in species with structured populations or those that have experienced periods of population expansion or decline. Selection can affect patterns of mitochondrial DNA variation and variable mutation rates among mitochondrial genes can compromise inferences drawn from single markers. We investigated the contribution of these factors to patterns of mitochondrial variation and estimates of time to most recent common ancestor (T_MRCA) for two clades in a co-operatively breeding avian species, the white-browed babbler Pomatostomus superciliosus. Both the protein-coding ND3 gene and hypervariable domain I control region sequences showed departures from neutral expectations within the superciliosus clade, and a two-fold difference in T_MRCA estimates. Bayesian phylogenetic analysis provided evidence of departure from a strict clock model of molecular evolution in domain I, leading to an over-estimation of T_MRCA for the superciliosus clade at this marker. Our results suggest mitochondrial studies that attempt to reconstruct Pleistocene demographic histories should rigorously evaluate data for departures from neutral equilibrium expectations, including variation in evolutionary rates across multiple markers. Failure to do so can lead to serious errors in the estimation of evolutionary parameters and subsequent demographic inferences concerning the role of climate as a driver of evolutionary change. These effects may be especially pronounced in species with complex social structures occupying heterogeneous environments. We propose that environmentally driven differences in social structure may explain observed differences in evolutionary rate of domain I sequences, resulting from longer than expected retention times for matriarchal lineages in the superciliosus clade.     

Evolutionary convergence of alternative splicing in ion channels

In Drosophila melanogaster and humans, members of three different ion-channel gene families share tandem exon duplications, which are alternatively spliced. In this article, I demonstrate that the duplication events that give rise to these mutually exclusive exons are unlikely to be ancestral but have probably occurred independently in different lineages. These events provide remarkable examples of evolutionary convergence in alternative splicing. The result has important implications for the analysis of regulation of alternative splicing using comparative genomics and our understanding of molecular evolution.     

The duplication of an eight-residue helical stretch in Staphylococcal nuclease is not helical: a model for evolutionary change

A common method of evolutionary change is gene duplication, followed by other events that lead to new function, decoration of folds, oligomerization, or other changes. As part of a study on the potential for evolutionary change created by duplicated sequences, we have carried out a crystallographic study on a mutant of Staphylococcal nuclease in which residues 55-62 have been duplicated in a wild-type variant termed PHS. In the parental protein (PHS) these residues form the first two turns of a helix running from residue 54 to 68 (hereafter designated as helix I). The crystal structure of the mutant is very similar to that of the parental, with helix I being unaltered. The duplicated residues are accommodated by expanding an existing loop N-terminal to helix I. In addition, circular dichroism (CD) studies have been carried out on a parental peptide containing helix I with six flanking residues at each terminus (residues 48-74) and on the same peptide expanded by the duplication, as a function of 2,2,2-trifluoroethanol (TFE) concentration. Each peptide possesses only modest helical propensity in solution. Our data, which is different from what was observed in T4 lysozyme, show that the conformation of the duplicated sequence is determined by a balance of sequential and longer-range effects. Thus duplicating sequence need not mean duplicating structure. Proteins 2000;40:465-472.     

Molecular evolution of the antiretroviral <Emphasis Type="Italic">TRIM5</Emphasis> gene

In 2004, the first report of TRIM5α as a cellular antiretroviral factor triggered intense interest among virologists, particularly because some primate orthologs of TRIM5α have activity against HIV. Since that time, a complex and eventful evolutionary history of the TRIM5 locus has emerged. A review of the TRIM5 literature constitutes a veritable compendium of evolutionary phenomena, including elevated rates of nonsynonymous substitution, divergence in subdomains due to short insertions and deletions, expansions and contractions in gene copy number, pseudogenization, balanced polymorphism, trans-species polymorphism, convergent evolution, and the acquisition of new domains by exon capture. Unlike most genes, whose history is dominated by long periods of purifying selection interspersed with rare instances of genetic innovation, analysis of restriction factor loci is likely to be complicated by the unpredictable and more-or-less constant influence of positive selection. In this regard, the molecular evolution and population genetics of restriction factor loci most closely resemble patterns that have been documented for immunity genes, such as class I and II MHC genes, whose products interact directly with microbial targets. While the antiretroviral activity encoded by TRIM5 provides plausible mechanistic hypotheses for these unusual evolutionary observations, evolutionary analyses have reciprocated by providing significant insights into the structure and function of the TRIM5α protein. Many of the lessons learned from TRIM5 should be applicable to the study of other restriction factor loci, and molecular evolutionary analysis could facilitate the discovery of new antiviral factors, particularly among the many TRIM genes whose functions remain as yet unidentified.     

The evolution of gene regulatory networks controlling Arabidopsis thaliana L. trichome development

Background

The variation in structure and function of gene regulatory networks (GRNs) participating in organisms development is a key for understanding species-specific evolutionary strategies. Even the tiniest modification of developmental GRN might result in a substantial change of a complex morphogenetic pattern. Great variety of trichomes and their accessibility makes them a useful model for studying the molecular processes of cell fate determination, cell cycle control and cellular morphogenesis. Nowadays, a large number of genes regulating the morphogenesis of A. thaliana trichomes are described. Here we aimed at a study the evolution of the GRN defining the trichome formation, and evaluation its importance in other developmental processes.

Results

In study of the evolution of trichomes formation GRN we combined classical phylogenetic analysis with information on the GRN topology and composition in major plants taxa. This approach allowed us to estimate both times of evolutionary emergence of the GRN components which are mainly proteins, and the relative rate of their molecular evolution. Various simplifications of protein structure (based on the position of amino acid residues in protein globula, secondary structure type, and structural disorder) allowed us to demonstrate the evolutionary associations between changes in protein globules and speciations/duplications events. We discussed their potential involvement in protein-protein interactions and GRN function.

Conclusions

We hypothesize that the divergence and/or the specialization of the trichome-forming GRN is linked to the emergence of plant taxa. Information about the structural targets of the protein evolution in the GRN may predict switching points in gene networks functioning in course of evolution. We also propose a list of candidate genes responsible for the development of trichomes in a wide range of plant species.

     

MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment

With its theoretical basis firmly established in molecular evolutionary and population genetics, the comparative DNA and protein sequence analysis plays a central role in reconstructing the evolutionary histories of species and multigene families, estimating rates of molecular evolution, and inferring the nature and extent of selective forces shaping the evolution of genes and genomes. The scope of these investigations has now expanded greatly owing to the development of high-throughput sequencing techniques and novel statistical and computational methods. These methods require easy-to-use computer programs. One such effort has been to produce Molecular Evolutionary Genetics Analysis (MEGA) software, with its focus on facilitating the exploration and analysis of the DNA and protein sequence variation from an evolutionary perspective. Currently in its third major release, MEGA3 contains facilities for automatic and manual sequence alignment, web-based mining of databases, inference of the phylogenetic trees, estimation of evolutionary distances and testing evolutionary hypotheses. This paper provides an overview of the statistical methods, computational tools, and visual exploration modules for data input and the results obtainable in MEGA.     

Contrasting recombination patterns and demographic histories of the plant pathogen Ralstonia solanacearum inferred from MLSA

We used multilocus sequence analysis (MLSA) on a worldwide collection of the plant pathogenic Ralstonia solanacearum (Betaproteobacteria) to retrace its complex evolutionary history. Using genetic imprints left during R. solanacearum evolution, we were able to delineate distinct evolutionary complex displaying contrasting dynamics. Among the phylotypes already described (I, IIA, IIB, III, IV), eight groups of strains with distinct evolutionary patterns, named clades, were identified. From our recombination analysis, we identified 21 recombination events that occurred within and across these lineages. Although appearing the most divergent and ancestral phylotype, phylotype IV was inferred as a gene donor for the majority of the recombination events that we detected. Whereas this phylotype apparently fuelled the species diversity, ongoing diversification was mainly detected within phylotype I, IIA and III. These three groups presented a recent expanding population structure, a high level of homologous recombination and evidences of long-distance migrations. Factors such as adaptation to a specific host or intense trading of infected crops may have promoted this diversification. Whether R. solanacearum lineages will eventually evolve in distinct species remains an open question. The intensification of cropping and increase of geographical dispersion may favour situations of phylotype sympatry and promote higher exchange of key factors for host adaptation from their common genetic pool.     

分子系统地理学及其应用

分子系统地理学是20世纪70年代中期伴随着对线粒体DNA的认识而开始酝酿发展的,本文回顾了该学科自诞生以来的发展简史,阐述了其研究内容,着重从3个方面介绍了该领域在其它相关学科中的研究进展,总结了该学科存在的问题,并提出重点应解决的3个问题：1）将细胞器基因与遗传信息更为丰富的核基因相结合;2）合生理论在非平衡种群中的发展应用;3）提高对分化时间估计的精确性。     

Docking protein domains in contact space

Background  

Many biological processes involve the physical interaction between protein domains. Understanding these functional associations requires knowledge of the molecular structure. Experimental investigations though present considerable difficulties and there is therefore a need for accurate and reliable computational methods. In this paper we present a novel method that seeks to dock protein domains using a contact map representation. Rather than providing a full three dimensional model of the complex, the method predicts contacting residues across the interface. We use a scoring function that combines structural, physicochemical and evolutionary information, where each potential residue contact is assigned a value according to the scoring function and the hypothesis is that the real configuration of contacts is the one that maximizes the score. The search is performed with a simulated annealing algorithm directly in contact space.     

Evolutionary Sequence Modeling for Discovery of Peptide Hormones

There are currently a large number of “orphan” G-protein-coupled receptors (GPCRs) whose endogenous ligands (peptide hormones) are unknown. Identification of these peptide hormones is a difficult and important problem. We describe a computational framework that models spatial structure along the genomic sequence simultaneously with the temporal evolutionary path structure across species and show how such models can be used to discover new functional molecules, in particular peptide hormones, via cross-genomic sequence comparisons. The computational framework incorporates a priori high-level knowledge of structural and evolutionary constraints into a hierarchical grammar of evolutionary probabilistic models. This computational method was used for identifying novel prohormones and the processed peptide sites by producing sequence alignments across many species at the functional-element level. Experimental results with an initial implementation of the algorithm were used to identify potential prohormones by comparing the human and non-human proteins in the Swiss-Prot database of known annotated proteins. In this proof of concept, we identified 45 out of 54 prohormones with only 44 false positives. The comparison of known and hypothetical human and mouse proteins resulted in the identification of a novel putative prohormone with at least four potential neuropeptides. Finally, in order to validate the computational methodology, we present the basic molecular biological characterization of the novel putative peptide hormone, including its identification and regional localization in the brain. This species comparison, HMM-based computational approach succeeded in identifying a previously undiscovered neuropeptide from whole genome protein sequences. This novel putative peptide hormone is found in discreet brain regions as well as other organs. The success of this approach will have a great impact on our understanding of GPCRs and associated pathways and help to identify new targets for drug development.     

Unique features of the pre-T-cell receptor alpha-chain: not just a surrogate

The pre-T-cell receptor (pre-TCR) has a crucial role in the normal development of alphabeta T cells. Different views have emerged concerning the structure and function of the pre-TCR. This molecular complex can be viewed as a variant of the alphabeta-TCR in which the pre-TCR alpha-chain that is covalently associated with the TCR beta-chain is a 'surrogate' TCR alpha-chain. Alternatively, the unique structure of the pre-TCR might be associated with a unique function, owing to evolutionary selection of a pre-TCR alpha-chain that has different capabilities from the TCR alpha-chain. As described here, I consider that experimental evidence favours the latter view.