Relative contributions of structural designability and functional diversity in molecular evolution of duplicates

Analysis of increasingly saturated sequence databases have shown that gene family sizes are highly skewed with many families being small and few containing many, far-diverged homologs. Additionally, recently published results have identified a structural determinant of mutational plasticity: designability that correlates strongly with gene family size. In this paper, we explore the possible links between the two observations, exploring the possible effect of designability on duplication and divergence. We show that designability has an inverse of expected relationship with strength of selection. More designable domains that should have more mutational plasticity evolve slower. However, we also present evidence that recently duplicated genes have variable probability of locus fixation correlated with strength of selection. As expected, paralogs under stronger evolutionary pressure have a lower failure rate. Finally, we show that probability of pseudogene formation from gene duplication can be directly tied to designability and functional flexibility of the family. We present evidence that gene families with higher designability have diverged farther because of lower probability of pseudogenization. Additionally, mutational plasticity may play an integral role by influencing pseudogenization rate. Either way, we show that considering the failure rate of duplications is integral in understanding the determinants and dynamics of molecular evolution.     

Role of pleiotropy in the evolution of a cryptic developmental variation in Caenorhabditis elegans

Robust biological systems are expected to accumulate cryptic genetic variation that does not affect the system output in standard conditions yet may play an evolutionary role once phenotypically expressed under a strong perturbation. Genetic variation that is cryptic relative to a robust trait may accumulate neutrally as it does not change the phenotype, yet it could also evolve under selection if it affects traits related to fitness in addition to its cryptic effect. Cryptic variation affecting the vulval intercellular signaling network was previously uncovered among wild isolates of Caenorhabditis elegans. Using a quantitative genetic approach, we identify a non-synonymous polymorphism of the previously uncharacterized nath-10 gene that affects the vulval phenotype when the system is sensitized with different mutations, but not in wild-type strains. nath-10 is an essential protein acetyltransferase gene and the homolog of human NAT10. The nath-10 polymorphism also presents non-cryptic effects on life history traits. The nath-10 allele carried by the N2 reference strain leads to a subtle increase in the egg laying rate and in the total number of sperm, a trait affecting the trade-off between fertility and minimal generation time in hermaphrodite individuals. We show that this allele appeared during early laboratory culture of N2, which allowed us to test whether it may have evolved under selection in this novel environment. The derived allele indeed strongly outcompetes the ancestral allele in laboratory conditions. In conclusion, we identified the molecular nature of a cryptic genetic variation and characterized its evolutionary history. These results show that cryptic genetic variation does not necessarily accumulate neutrally at the whole-organism level, but may evolve through selection for pleiotropic effects that alter fitness. In addition, cultivation in the laboratory has led to adaptive evolution of the reference strain N2 to the laboratory environment, which may modify other phenotypes of interest.     

The role of developmental plasticity in evolutionary innovation

Explaining the origins of novel traits is central to evolutionary biology. Longstanding theory suggests that developmental plasticity, the ability of an individual to modify its development in response to environmental conditions, might facilitate the evolution of novel traits. Yet whether and how such developmental flexibility promotes innovations that persist over evolutionary time remains unclear. Here, we examine three distinct ways by which developmental plasticity can promote evolutionary innovation. First, we show how the process of genetic accommodation provides a feasible and possibly common avenue by which environmentally induced phenotypes can become subject to heritable modification. Second, we posit that the developmental underpinnings of plasticity increase the degrees of freedom by which environmental and genetic factors influence ontogeny, thereby diversifying targets for evolutionary processes to act on and increasing opportunities for the construction of novel, functional and potentially adaptive phenotypes. Finally, we examine the developmental genetic architectures of environment-dependent trait expression, and highlight their specific implications for the evolutionary origin of novel traits. We critically review the empirical evidence supporting each of these processes, and propose future experiments and tests that would further illuminate the interplay between environmental factors, condition-dependent development, and the initiation and elaboration of novel phenotypes.     

An evolutionary model for protein-coding regions with conserved RNA structure

Here we present a model of nucleotide substitution in protein-coding regions that also encode the formation of conserved RNA structures. In such regions, apparent evolutionary context dependencies exist, both between nucleotides occupying the same codon and between nucleotides forming a base pair in the RNA structure. The overlap of these fundamental dependencies is sufficient to cause "contagious" context dependencies which cascade across many nucleotide sites. Such large-scale dependencies challenge the use of traditional phylogenetic models in evolutionary inference because they explicitly assume evolutionary independence between short nucleotide tuples. In our model we address this by replacing context dependencies within codons by annotation-specific heterogeneity in the substitution process. Through a general procedure, we fragment the alignment into sets of short nucleotide tuples based on both the protein coding and the structural annotation. These individual tuples are assumed to evolve independently, and the different tuple sets are assigned different annotation-specific substitution models shared between their members. This allows us to build a composite model of the substitution process from components of traditional phylogenetic models. We applied this to a data set of full-genome sequences from the hepatitis C virus where five RNA structures are mapped within the coding region. This allowed us to partition the effects of selection on different structural elements and to test various hypotheses concerning the relation of these effects. Of particular interest, we found evidence of a functional role of loop and bulge regions, as these were shown to evolve according to a different and more constrained selective regime than the nonpairing regions outside the RNA structures. Other potential applications of the model include comparative RNA structure prediction in coding regions and RNA virus phylogenetics.     

Structural mapping of protein interactions reveals differences in evolutionary pressures correlated to mRNA level and protein abundance

Genome-wide studies in Saccharomyces cerevisiae concluded that the dominant determinant of protein evolutionary rates is expression level: highly expressed proteins generally evolve most slowly. To determine how this constraint affects the evolution of protein interactions, we directly measure evolutionary rates of protein interface, surface, and core residues by structurally mapping domain interactions to yeast genomes. We find that mRNA level and protein abundance, though correlated, report on pressures affecting regions of proteins differently. Pressures proportional to mRNA level slow evolutionary rates of all structural regions and reduce the variability in rate differences between interfaces and other surfaces. In contrast, the evolutionary rate variation within a domain is much less correlated to protein abundance. Distinct pressures may be associated primarily with the cost (mRNA level) and functional (protein abundance) benefit of protein production. Interfaces of proteins with low mRNA levels may have higher evolutionary flexibility and could constitute the raw material for new functions.     

The vertebrate social behavior network: evolutionary themes and variations

Based on a wide variety of data, it is now clear that birds and teleost (bony) fish possess a core "social behavior network" within the basal forebrain and midbrain that is homologous to the social behavior network of mammals. The nodes of this network are reciprocally connected, contain receptors for sex steroid hormones, and are involved in multiple forms of social behavior. Other hodological features and neuropeptide distributions are likewise very similar across taxa. This evolutionary conservation represents a boon for experiments on phenotypic behavioral variation, as the extraordinary social diversity of teleost fish and songbirds can now be used to generate broadly relevant insights into issues of brain function that are not particularly tractable in other vertebrate groups. Two such lines of research are presented here, each of which addresses functional variation within the network as it relates to divergent patterns of social behavior. In the first set of experiments, we have used a sexually polymorphic fish to demonstrate that natural selection can operate independently on hypothalamic neuroendocrine functions that are relevant for (1) gonadal regulation and (2) sex-typical behavioral modulation. In the second set of experiments, we have exploited the diversity of avian social organizations and ecologies to isolate species-typical group size as a quasi-independent variable. These experiments have shown that specific areas and peptidergic components of the social behavior network possess functional properties that evolve in parallel with divergence and convergence in sociality.     

Inverse relationship between evolutionary rate and age of mammalian genes

A large number of genes is shared by all living organisms, whereas many others are unique to some specific lineages, indicating their different times of origin. The availability of a growing number of eukaryotic genomes allows us to estimate which mammalian genes are novel genes and, approximately, when they arose. In this article, we classify human genes into four different age groups and estimate evolutionary rates in human and mouse orthologs. We show that older genes tend to evolve more slowly than newer ones; that is, proteins that arose earlier in evolution currently have a larger proportion of sites subjected to negative selection. Interestingly, this property is maintained when a fraction of the fastest-evolving genes is excluded or when only genes belonging to a given functional class are considered. One way to explain this relationship is by assuming that genes maintain their functional constraints along all their evolutionary history, but the nature of more recent evolutionary innovations is such that the functional constraints operating on them are increasingly weaker. Alternatively, our results would also be consistent with a scenario in which the functional constraints acting on a gene would not need to be constant through evolution. Instead, starting from weak functional constraints near the time of origin of a gene-as supported by mechanisms proposed for the origin of orphan genes-there would be a gradual increase in selective pressures with time, resulting in fewer accepted mutations in older versus more novel genes.     

In silico identification of functional protein interfaces

Proteins perform many of their biological roles through protein-protein, protein-DNA or protein-ligand interfaces. The identification of the amino acids comprising these interfaces often enhances our understanding of the biological function of the proteins. Many methods for the detection of functional interfaces have been developed, and large-scale analyses have provided assessments of their accuracy. Among them are those that consider the size of the protein interface, its amino acid composition and its physicochemical and geometrical properties. Other methods to this effect use statistical potential functions of pairwise interactions, and evolutionary information. The rationale of the evolutionary approach is that functional and structural constraints impose selective pressure; hence, biologically important interfaces often evolve at a slower pace than do other external regions of the protein. Recently, an algorithm, Rate4Site, and a web-server, ConSurf (http://consurf.tau.ac.il/), for the identification of functional interfaces based on the evolutionary relations among homologous proteins as reflected in phylogenetic trees, were developed in our laboratory. The explicit use of the tree topology and branch lengths makes the method remarkably accurate and sensitive. Here we demonstrate its potency in the identification of the functional interfaces of a hypothetical protein, the structure of which was determined as part of the international structural genomics effort. Finally, we propose to combine complementary procedures, in order to enhance the overall performance of methods for the identification of functional interfaces in proteins.     

Imprint of evolutionary conservation and protein structure variation on the binding function of protein tyrosine kinases

MOTIVATION: According to the models of divergent molecular evolution, the evolvability of new protein function may depend on the induction of new phenotypic traits by a small number of mutations of the binding site residues. Evolutionary relationships between protein kinases are often employed to infer inhibitor binding profiles from sequence analysis. However, protein kinases binding profiles may display inhibitor selectivity within a given kinase subfamily, while exhibiting cross-activity between kinases that are phylogenetically remote from the prime target. The emerging insights into kinase function and evolution combined with a rapidly growing number of publically available crystal structures of protein kinases complexes have motivated structural bioinformatics analysis of sequence-structure relationships in determining the binding function of protein tyrosine kinases. RESULTS: In silico profiling of Imatinib mesylate and PD-173955 kinase inhibitors with protein tyrosine kinases is conducted on kinome scale by using evolutionary analysis and fingerprinting inhibitor-protein interactions with the panel of all publically available protein tyrosine kinases crystal structures. We have found that sequence plasticity of the binding site residues alone may not be sufficient to enable protein tyrosine kinases to readily evolve novel binding activities with inhibitors. While evolutionary signal derived solely from the tyrosine kinase sequence conservation can not be readily translated into the ligand binding phenotype, the proposed structural bioinformatics analysis can discriminate a functionally relevant kinase binding signal from a simple phylogenetic relationship. The results of this work reveal that protein conformational diversity is intimately linked with sequence plasticity of the binding site residues in achieving functional adaptability of protein kinases towards specific drug binding. This study offers a plausible molecular rationale to the experimental binding profiles of the studied kinase inhibitors and provides a theoretical basis for constructing functionally relevant kinase binding trees.     

Mitochondrial DNA as a marker of molecular diversity: a reappraisal

Over the last three decades, mitochondrial DNA has been the most popular marker of molecular diversity, for a combination of technical ease-of-use considerations, and supposed biological and evolutionary properties of clonality, near-neutrality and clock-like nature of its substitution rate. Reviewing recent literature on the subject, we argue that mitochondrial DNA is not always clonal, far from neutrally evolving and certainly not clock-like, questioning its relevance as a witness of recent species and population history. We critically evaluate the usage of mitochondrial DNA for species delineation and identification. Finally, we note the great potential of accumulating mtDNA data for evolutionary and functional analysis of the mitochondrial genome.     

Contact Density Affects Protein Evolutionary Rate from Bacteria to Animals

The density of contacts or the fraction of buried sites in a protein structure is thought to be related to a protein’s designability, and genes encoding more designable proteins should evolve faster than other genes. Several recent studies have tested this hypothesis but have found conflicting results. Here, we investigate how a gene’s evolutionary rate is affected by its protein’s contact density, considering the four species Escherichia coli, Saccharomyces cerevisiae, Drosophila melanogaster, and Homo sapiens. We find for all four species that contact density correlates positively with evolutionary rate, and that these correlations do not seem to be confounded by gene expression level. The strength of this signal, however, varies widely among species. We also study the effect of contact density on domain evolution in multidomain proteins and find that a domain’s contact density influences the domain’s evolutionary rate. Within the same protein, a domain with higher contact density tends to evolve faster than a domain with lower contact density. Our study provides evidence that contact density can increase evolutionary rates, and that it acts similarly on the level of entire proteins and of individual protein domains. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Protein subcellular relocalization: a new perspective on the origin of novel genes

Gene duplication is considered to be the most important evolutionary process for generating novel genes. However, the mechanisms involved in the evolution of such genetic innovations remain unclear. There is compelling evidence to suggest that changing the subcellular location of a protein can also alter its function, and that diversity in subcellular targeting within gene families is common. Here, we introduce the idea that protein subcellular relocalization might be an important evolutionary mechanism for the origins of new genes.     

Prediction of site-specific amino acid distributions and limits of divergent evolutionary changes in protein sequences

We derive an analytic expression for site-specific stationary distributions of amino acids from the structurally constrained neutral (SCN) model of protein evolution with conservation of folding stability. The stationary distributions that we obtain have a Boltzmann-like shape, and their effective temperature parameter, measuring the limit of divergent evolutionary changes at a given site, can be predicted from a site-specific topological property, the principal eigenvector of the contact matrix of the native conformation of the protein. These analytic results, obtained without free parameters, are compared with simulations of the SCN model and with the site-specific amino acid distributions obtained from the Protein Data Bank. These results also provide new insights into how the topology of a protein fold influences its designability, i.e., the number of sequences compatible with that fold. The dependence of the effective temperature on the principal eigenvector decreases for longer proteins, as a possible consequence of the fact that selection for thermodynamic stability becomes weaker in this case.     

Eco‐evolutionary feedbacks between prey densities and linkage disequilibrium in the predator maintain diversity

Diversity occurs at multiple scales. Within a single population, there is diversity in genotypes and phenotypes. At a larger scale, within ecological communities, there is diversity in species. A number of studies have investigated how diversity at these two scales influence each other through what has been termed eco‐evolutionary feedbacks. Here we study a three‐species ecological module called apparent competition, in which the predator is evolving in a trait that determines its interaction with two prey species. Unlike previous studies on apparent competition, which employed evolutionary frameworks with very simple genetics, we study an eco‐evolutionary model in which the predator's trait is determined by two recombining diallelic loci, so that its mean and variance can evolve, as well as associations (linkage disequilibrium) between the loci. We ask how eco‐evolutionary feedbacks with these two loci affect the coexistence of the prey species and the maintenance of polymorphisms within the predator species. We uncover a novel eco‐evolutionary feedback between the prey densities and the linkage disequilibrium between the predator's loci. Through a stability analysis, we demonstrate how these feedbacks affect polymorphisms at both loci and, among others, may generate stable cycling.     

Odonata (dragonflies and damselflies) as a bridge between ecology and evolutionary genomics

Odonata (dragonflies and damselflies) present an unparalleled insect model to integrate evolutionary genomics with ecology for the study of insect evolution. Key features of Odonata include their ancient phylogenetic position, extensive phenotypic and ecological diversity, several unique evolutionary innovations, ease of study in the wild and usefulness as bioindicators for freshwater ecosystems worldwide. In this review, we synthesize studies on the evolution, ecology and physiology of odonates, highlighting those areas where the integration of ecology with genomics would yield significant insights into the evolutionary processes that would not be gained easily by working on other animal groups. We argue that the unique features of this group combined with their complex life cycle, flight behaviour, diversity in ecological niches and their sensitivity to anthropogenic change make odonates a promising and fruitful taxon for genomics focused research. Future areas of research that deserve increased attention are also briefly outlined.     

Beyond the evolutionary theory of ageing, from functional genomics to evo-gero

By the mid 1970s, the mechanisms by which ageing can evolve had a secure theoretical basis in population genetics. Here, we discuss how subsequent evolutionary work has focussed on testing and extending this theory, and on attempting to integrate it with other emerging facets of the biology of ageing, such as genetic studies of long-lived mutants and of phenotypic plasticity in ageing, such as in response to nutritional status. We also describe how functional genomic studies are providing new insights into the evolutionary forces shaping genome evolution and lifespan control. Future challenges include understanding the biochemistry of longevity and how its failure generates ageing and associated diseases, and the determination of the genetic basis of lifespan evolution and the great plasticity that it displays.     

Heterostyly accelerates diversification via reduced extinction in primroses

The exceptional species diversity of flowering plants, exceeding that of their sister group more than 250-fold, is especially evident in floral innovations, interactions with pollinators and sexual systems. Multiple theories, emphasizing flower–pollinator interactions, genetic effects of mating systems or high evolvability, predict that floral evolution profoundly affects angiosperm diversification. However, consequences for speciation and extinction dynamics remain poorly understood. Here, we investigate trajectories of species diversification focusing on heterostyly, a remarkable floral syndrome where outcrossing is enforced via cross-compatible floral morphs differing in placement of their respective sexual organs. Heterostyly evolved at least 20 times independently in angiosperms. Using Darwin''s model for heterostyly, the primrose family, we show that heterostyly accelerates species diversification via decreasing extinction rates rather than increasing speciation rates, probably owing to avoidance of the negative genetic effects of selfing. However, impact of heterostyly appears to differ over short and long evolutionary time-scales: the accelerating effect of heterostyly on lineage diversification is manifest only over long evolutionary time-scales, whereas recent losses of heterostyly may prompt ephemeral bursts of speciation. Our results suggest that temporal or clade-specific conditions may ultimately determine the net effects of specific traits on patterns of species diversification.