Evolutionary significance of gene expression divergence

Recent large-scale studies of evolutionary changes in gene expression among mammalian species have led to the proposal that gene expression divergence may be neutral with respect to organismic fitness. Here, we employ a comparative analysis of mammalian gene sequence divergence and gene expression divergence to test the hypothesis that the evolution of gene expression is predominantly neutral. Two models of neutral gene expression evolution are considered: 1-purely neutral evolution (i.e., no selective constraint) of gene expression levels and patterns and 2-neutral evolution accompanied by selective constraint. With respect to purely neutral evolution, levels of change in gene expression between human-mouse orthologs are correlated with levels of gene sequence divergence that are determined largely by purifying selection. In contrast, evolutionary changes of tissue-specific gene expression profiles do not show such a correlation with sequence divergence. However, divergence of both gene expression levels and profiles are significantly lower for orthologous human-mouse gene pairs than for pairs of randomly chosen human and mouse genes. These data clearly point to the action of selective constraint on gene expression divergence and are inconsistent with the purely neutral model; however, there is likely to be a neutral component in evolution of gene expression, particularly, in tissues where the expression of a given gene is low and functionally irrelevant. The model of neutral evolution with selective constraint predicts a regular, clock-like accumulation of gene expression divergence. However, relative rate tests of the divergence among human-mouse-rat orthologous gene sets reveal clock-like evolution for gene sequence divergence, and to a lesser extent for gene expression level divergence, but not for the divergence of tissue-specific gene expression profiles. Taken together, these results indicate that gene expression divergence is subject to the effects of purifying selective constraint and suggest that it might also be substantially influenced by positive Darwinian selection.     

Conservation and coevolution in the scale-free human gene coexpression network

The role of natural selection in biology is well appreciated. Recently, however, a critical role for physical principles of network self-organization in biological systems has been revealed. Here, we employ a systems level view of genome-scale sequence and expression data to examine the interplay between these two sources of order, natural selection and physical self-organization, in the evolution of human gene regulation. The topology of a human gene coexpression network, derived from tissue-specific expression profiles, shows scale-free properties that imply evolutionary self-organization via preferential node attachment. Genes with numerous coexpressed partners (the hubs of the coexpression network) evolve more slowly on average than genes with fewer coexpressed partners, and genes that are coexpressed show similar rates of evolution. Thus, the strength of selective constraints on gene sequences is affected by the topology of the gene coexpression network. This connection is strong for the coding regions and 3' untranslated regions (UTRs), but the 5' UTRs appear to evolve under a different regime. Surprisingly, we found no connection between the rate of gene sequence divergence and the extent of gene expression profile divergence between human and mouse. This suggests that distinct modes of natural selection might govern sequence versus expression divergence, and we propose a model, based on rapid, adaptation-driven divergence and convergent evolution of gene expression patterns, for how natural selection could influence gene expression divergence.     

Molecular evolution and network-level analysis of the N-glycosylation metabolic pathway across primates

N-glycosylation is one of the most important forms of protein modification, serving key biological functions in multicellular organisms. N-glycans at the cell surface mediate the interaction between cells and the surrounding matrix and may act as pathogen receptors, making the genes responsible for their synthesis good candidates to show signatures of adaptation to different pathogen environments. Here, we study the forces that shaped the evolution of the genes involved in the synthesis of the N-glycans during the divergence of primates within the framework of their functional network. We have found that, despite their function of producing glycan repertoires capable of evading rapidly evolving pathogens, genes involved in the synthesis of the glycans are highly conserved, and no signals of positive selection have been detected within the time of divergence of primates. This suggests strong functional constraints as the main force driving their evolution. We studied the strength of the purifying selection acting on the genes in relation to the network structure considering the position of each gene along the pathway, its connectivity, and the rates of evolution in neighboring genes. We found a strong and highly significant negative correlation between the strength of purifying selection and the connectivity of each gene, indicating that genes encoding for highly connected enzymes evolve slower and thus are subject to stronger selective constraints. This result confirms that network topology does shape the evolution of the genes and that the connectivity within metabolic pathways and networks plays a major role in constraining evolutionary rates.     

Evolution of Gene Regulatory Networks by Fluctuating Selection and Intrinsic Constraints

Various characteristics of complex gene regulatory networks (GRNs) have been discovered during the last decade, e.g., redundancy, exponential indegree distributions, scale-free outdegree distributions, mutational robustness, and evolvability. Although progress has been made in this field, it is not well understood whether these characteristics are the direct products of selection or those of other evolutionary forces such as mutational biases and biophysical constraints. To elucidate the causal factors that promoted the evolution of complex GRNs, we examined the effect of fluctuating environmental selection and some intrinsic constraining factors on GRN evolution by using an individual-based model. We found that the evolution of complex GRNs is remarkably promoted by fixation of beneficial gene duplications under unpredictably fluctuating environmental conditions and that some internal factors inherent in organisms, such as mutational bias, gene expression costs, and constraints on expression dynamics, are also important for the evolution of GRNs. The results indicate that various biological properties observed in GRNs could evolve as a result of not only adaptation to unpredictable environmental changes but also non-adaptive processes owing to the properties of the organisms themselves. Our study emphasizes that evolutionary models considering such intrinsic constraining factors should be used as null models to analyze the effect of selection on GRN evolution.     

A selection index for gene expression evolution and its application to the divergence between humans and chimpanzees

The importance of gene regulation in animal evolution is a matter of long-standing interest, but measuring the impact of selection on gene expression has proven a challenge. Here, we propose a selection index of gene expression as a straightforward method for assessing the mode and strength of selection operating on gene expression levels. The index is based on the widely used McDonald-Kreitman test and requires the estimation of four quantities: the within-species and between-species expression variances as well as the sequence heterozygosity and divergence of neutrally evolving sequences. We apply the method to data from human and chimpanzee lymphoblastoid cell lines and show that gene expression is in general under strong stabilizing selection. We also demonstrate how the same framework can be used to estimate the proportion of adaptive gene expression evolution.     

Genes under weaker stabilizing selection increase network evolvability and rapid regulatory adaptation to an environmental shift

Regulatory networks play a central role in the modulation of gene expression, the control of cellular differentiation, and the emergence of complex phenotypes. Regulatory networks could constrain or facilitate evolutionary adaptation in gene expression levels. Here, we model the adaptation of regulatory networks and gene expression levels to a shift in the environment that alters the optimal expression level of a single gene. Our analyses show signatures of natural selection on regulatory networks that both constrain and facilitate rapid evolution of gene expression level towards new optima. The analyses are interpreted from the standpoint of neutral expectations and illustrate the challenge to making inferences about network adaptation. Furthermore, we examine the consequence of variable stabilizing selection across genes on the strength and direction of interactions in regulatory networks and in their subsequent adaptation. We observe that directional selection on a highly constrained gene previously under strong stabilizing selection was more efficient when the gene was embedded within a network of partners under relaxed stabilizing selection pressure. The observation leads to the expectation that evolutionarily resilient regulatory networks will contain optimal ratios of genes whose expression is under weak and strong stabilizing selection. Altogether, our results suggest that the variable strengths of stabilizing selection across genes within regulatory networks might itself contribute to the long‐term adaptation of complex phenotypes.     

Mistranslation-induced protein misfolding as a dominant constraint on coding-sequence evolution

Strikingly consistent correlations between rates of coding-sequence evolution and gene expression levels are apparent across taxa, but the biological causes behind the selective pressures on coding-sequence evolution remain controversial. Here, we demonstrate conserved patterns of simple covariation between sequence evolution, codon usage, and mRNA level in E. coli, yeast, worm, fly, mouse, and human that suggest that all observed trends stem largely from a unified underlying selective pressure. In metazoans, these trends are strongest in tissues composed of neurons, whose structure and lifetime confer extreme sensitivity to protein misfolding. We propose, and demonstrate using a molecular-level evolutionary simulation, that selection against toxicity of misfolded proteins generated by ribosome errors suffices to create all of the observed covariation. The mechanistic model of molecular evolution that emerges yields testable biochemical predictions, calls into question the use of nonsynonymous-to-synonymous substitution ratios (Ka/Ks) to detect functional selection, and suggests how mistranslation may contribute to neurodegenerative disease.     

Rates of Evolution of Hominoid Seminal Proteins are Correlated with Function and Expression, Rather than Mating System

In species where females mate with multiple males during a single ovulatory cycle, sperm competition is hypothesized to increase the rate of adaptive evolution of proteins expressed in male reproductive tissues through recurrent selective sweeps (positive selection). The hominoids, comprising apes and humans, are a group of closely related primates with extensive variation in mating behaviors and predicted levels of sperm competition. Since previous studies of individual male reproductive genes have shown evidence of positive selection, we estimated rates of evolution of a comprehensive set of proteins expressed in ejaculated semen. Our results show that these proteins in hominoids do not have elevated rates of nonsynonymous substitutions (Ka) compared with a control dataset of nonreproductive genes. Species with greater sperm competition do not have faster rates of seminal protein evolution. Although at these broad levels our hypotheses were not confirmed, further analyses indicate specific patterns of molecular evolution. Namely, the Ka of seminal genes is more strongly correlated with measures of tissue specificity than nonreproductive genes, suggesting that the former may more readily adapt to tissue-specific functions. Proteins expressed from the seminal vesicles evolve more rapidly than those from other male reproductive tissues. Also, several gene ontology categories show elevated rates of protein evolution, not seen in the control data set. While the generalization that male reproductive genes evolve rapidly in hominoids is an oversimplification, a subset of proteins can be identified that are likely targets for adaptive evolution driven by sexual selection.     

The evolution of social philopatry in female primates

The transition from solitary life to sociality is considered one of the major transitions in evolution. In primates, this transition is currently not well understood. Traditional verbal models appear insufficient to unravel the complex interplay of environmental and demographic factors involved in the evolution of primate sociality, and recent phylogenetic reconstructions have produced conflicting results. We therefore analyze a theoretical model for the evolution of female social philopatry that sheds new light on the question why most primates live in groups. In individual-based simulations, we study the evolution of dispersal strategies of both resident females and their offspring. The model reveals that social philopatry can evolve through kin selection, even if retention of offspring is costly in terms of within-group resource competition and provides no direct benefits. Our model supports the role of predator avoidance as a selective pressure for group-living in primates, but it also suggests that a second benefit of group-living, communal resource defense, might be required to trigger the evolution of sizable groups. Lastly, our model reveals that seemingly small differences in demographic parameters can have profound effects on primate social evolution.     

Mouse models and the evolutionary developmental biology of the skull

Understanding development is relevant to understanding evolutionbecause developmental processes structure the expression ofphenotypic variation upon which natural selection acts. Advancesin developmental biology are fueling a new synthesis of developmentaland evolutionary biology, but it remains unclear how to usedevelopmental information that largely derives from a few modelorganisms to test hypotheses about the evolutionary developmentalbiology of taxa such as humans and other primates that havenot been or are not amenable to direct study through experimentaldevelopmental biology. In this article, we discuss how and whenmodel organisms like mice are useful for studying the evolutionarydevelopmental biology of even rather distantly related and morphologicallydifferent groups like primates. A productive approach is tofocus on processes that are likely to play key roles in producingevolutionarily significant phenotypic variation across a largephylogenetic range. We illustrate this approach by applyingthe analysis of craniofacial variation in mouse mutant modelsto primate and human evolution.     

Differential Evolution of MAGE Genes Based on Expression Pattern and Selection Pressure

Starting from publicly-accessible datasets, we have utilized comparative and phylogenetic genome analyses to characterize the evolution of the human MAGE gene family. Our characterization of genomic structures in representative genomes of primates, rodents, carnivora, and macroscelidea indicates that both Type I and Type II MAGE genes have undergone lineage-specific evolution. The restricted expression pattern in germ cells of Type I MAGE orthologs is observed throughout evolutionary history. Unlike Type II MAGEs that have conserved promoter sequences, Type I MAGEs lack promoter conservation, suggesting that epigenetic regulation is a central mechanism for controlling their expression. Codon analysis shows that Type I but not Type II MAGE genes have been under positive selection. The combination of genomic and expression analysis suggests that Type 1 MAGE promoters and genes continue to evolve in the hominin lineage, perhaps towards functional diversification or acquiring additional specific functions, and that selection pressure at codon level is associated with expression spectrum.     

Host-parasite interactions and the evolution of gene expression

Interactions between hosts and parasites provide an ongoing source of selection that promotes the evolution of a variety of features in the interacting species. Here, we use a genetically explicit mathematical model to explore how patterns of gene expression evolve at genetic loci responsible for host resistance and parasite infection. Our results reveal the striking yet intuitive conclusion that gene expression should evolve along very different trajectories in the two interacting species. Specifically, host resistance loci should frequently evolve to co-express alleles, whereas parasite infection loci should evolve to express only a single allele. This result arises because hosts that co-express resistance alleles are able to recognize and clear a greater diversity of parasite genotypes. By the same token, parasites that co-express antigen or elicitor alleles are more likely to be recognized and cleared by the host, and this favours the expression of only a single allele. Our model provides testable predictions that can help interpret accumulating data on expression levels for genes relevant to host-parasite interactions.     

The evolution of genomic imprinting via variance minimization: an evolutionary genetic model

A small number of mammalian loci exhibit genomic imprinting, in which only one copy of a gene is expressed while the other is silenced. At some such loci, the maternally inherited allele is inactivated; others show paternal inactivation. Several hypotheses have been put forward to explain how this genetic system could have evolved in the face of the selective advantages of diploidy. In this study, we examine the variance-minimization hypothesis, which proposes that imprinting arose through selection for reduced variation in levels of gene expression. We present an evolutionary genetic model incorporating both this selection pressure and deleterious mutations to elucidate the conditions under which imprinting could evolve. Our analysis implies that additional mechanisms such as genetic drift are required for imprinting to evolve from an initial nonimprinting state. Other predictions of this hypothesis do not appear to fit the available data as well as predictions for two alternative hypotheses, genetic conflict and the ovarian time bomb. On the basis of this evidence, we conclude that the variance-minimization hypothesis appears less adequate to explain the evolution of genomic imprinting.     

Comparative studies of gene expression and the evolution of gene regulation

The hypothesis that differences in gene regulation have an important role in speciation and adaptation is more than 40 years old. With the advent of new sequencing technologies, we are able to characterize and study gene expression levels and associated regulatory mechanisms in a large number of individuals and species at an unprecedented resolution and scale. We have thus gained new insights into the evolutionary pressures that shape gene expression levels and have developed an appreciation for the relative importance of evolutionary changes in different regulatory genetic and epigenetic mechanisms. The current challenge is to link gene regulatory changes to adaptive evolution of complex phenotypes. Here we mainly focus on comparative studies in primates and how they are complemented by studies in model organisms.     

The evolution of sex-biased genes and sex-biased gene expression

Differences between males and females in the optimal phenotype that is favoured by selection can be resolved by the evolution of differential gene expression in the two sexes. Microarray experiments have shown that such sex-biased gene expression is widespread across organisms and genomes. Sex-biased genes show unusually rapid sequence evolution, are often labile in their pattern of expression, and are non-randomly distributed in the genome. Here we discuss the characteristics and expression of sex-biased genes, and the selective forces that shape this previously unappreciated source of phenotypic diversity. Sex-biased gene expression has implications beyond just evolutionary biology, including for medical genetics.