Comparative network analysis reveals that tissue specificity and gene function are important factors influencing the mode of expression evolution in Arabidopsis and rice

Microarray experiments have yielded massive amounts of expression information measured under various conditions for the model species Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). Expression compendia grouping multiple experiments make it possible to define correlated gene expression patterns within one species and to study how expression has evolved between species. We developed a robust framework to measure expression context conservation (ECC) and found, by analyzing 4,630 pairs of orthologous Arabidopsis and rice genes, that 77% showed conserved coexpression. Examples of nonconserved ECC categories suggested a link between regulatory evolution and environmental adaptations and included genes involved in signal transduction, response to different abiotic stresses, and hormone stimuli. To identify genomic features that influence expression evolution, we analyzed the relationship between ECC, tissue specificity, and protein evolution. Tissue-specific genes showed higher expression conservation compared with broadly expressed genes but were fast evolving at the protein level. No significant correlation was found between protein and expression evolution, implying that both modes of gene evolution are not strongly coupled in plants. By integration of cis-regulatory elements, many ECC conserved genes were significantly enriched for shared DNA motifs, hinting at the conservation of ancestral regulatory interactions in both model species. Surprisingly, for several tissue-specific genes, patterns of concerted network evolution were observed, unveiling conserved coexpression in the absence of conservation of tissue specificity. These findings demonstrate that orthologs inferred through sequence similarity in many cases do not share similar biological functions and highlight the importance of incorporating expression information when comparing genes across species.     

Protein evolution: structure-function relationships of the oncogene beta-catenin in the evolution of multicellular animals

Beta-catenin functions as a cytoskeletal linker protein in cadherin-mediated adhesion and as a signal mediator in wnt-signal transduction pathways. We use a novel integrative approach, combining evolutionary, genomic, and three-dimensional structural data to analyze and trace the structural and functional evolution of beta-catenin genes. This approach also enabled us to examine the effects of gene duplication on the structure and function of beta-catenin genes in Drosophila, C. elegans, and vertebrates. By sampling a large number of different taxa, we identified both ancestral and derived motifs and residues within the different regions of the beta-catenin proteins. Projecting amino acid substitutions onto the three- dimensional structure established for mouse beta-catenin, we identified specific domains that exhibit loss and gain of selective constraints during beta catenin evolution. Structural changes, changes in the amino acid substitution rate, and the appearance of novel functional domains in beta-catenin can be mapped to specific branches on the metazoan tree. Together, our analyses suggest that a single, beta-catenin gene fulfilled both adhesion and signaling functions in the last common ancestor of metazoans some 700 million years ago. In addition, gene duplications facilitated the evolution of beta-catenins with novel functions and allowed the evolution of multiple, single-function proteins (cell adhesion or wnt-signaling) from the ancestral, dual-function protein. Integrative methods such as those we have applied here, utilizing the 'natural experiments' present in animal diversity, can be employed to identify novel and shared functional motifs and residues in virtually any protein among the proteomes of model systems and humans.     

Phylogenetic--evolutionary approaches to bioinformatics

Phylogenies of organisms are essential to investigating a range of evolutionary questions of interest to researchers in the field of bioinformatics. Phylogenies not only help to define how to study many evolutionary questions, they must also be taken into account when conducting statistical analyses. Here it is shown how phylogenies can be used to investigate variability along the sites of a gene, reconstruct ancestral states of ancient genes and proteins, identify and characterise events of parallel and convergent evolution, find events of gene duplication, analyse predictions from molecular clocks, seek evidence for correlated changes among different parts of the same gene or genome, and test theories of molecular evolution. A table of statistical and phylogenetic methods is presented.     

Functional divergence in the Arabidopsis beta-1,3-glucanase gene family inferred by phylogenetic reconstruction of expression states

Plant beta-1,3-glucanases (beta-1,3-Gs) (E.C. 3.2.1.39) comprise large, highly complex gene families involved in pathogen defense as well as a wide range of normal developmental processes. In spite of previous phylogenetic analyses that classify beta-1,3-Gs by sequence relatedness, the functional evolution of beta-1,3-Gs remains unclear. Here, expression and phylogenetic analyses have been integrated in order to investigate patterns of functional divergence in the Arabidopsis beta-1,3-G gene family. Fifty beta-1,3-G genes were grouped into expression classes through clustering of microarray data, and functions were inferred based on knowledge of coexpressed genes and existing literature. The resulting expression classes were mapped as discrete states onto a phylogenetic tree and parsimony reconstruction of ancestral expression states was performed, providing a model of expression divergence. Results showed a highly nonrandom distribution of developmental expression states in the phylogeny (P = 0.0002) indicating a significant degree of coupling between sequence and developmental expression divergence. A weaker, yet significant level of coupling was found using stress response data, but not using hormone-response or pathogen-response data. According to the model of developmental expression divergence, the ancestral function was most likely involved in cell division and/or cell wall remodeling. The associated expression state is widely distributed in the phylogeny, is retained by over 25% of gene family members, and is consistent with the known functions of beta-1,3-Gs in distantly related species and gene families. Consistent with previous hypotheses, pathogenesis-related (PR) beta-1,3-Gs appear to have evolved from ancestral developmentally regulated beta-1,3-Gs, acquiring PR function through a number of evolutionary events: divergence from the ancestral expression state, acquisition of pathogen/stress-responsive expression patterns, and loss of the C-terminal region including the glycosylphosphatidylinisotol (GPI)-anchoring site thus allowing for extracellular secretion.     

Ancestral inference and the study of codon bias evolution: implications for molecular evolutionary analyses of the Drosophila melanogaster subgroup

Reliable inference of ancestral sequences can be critical to identifying both patterns and causes of molecular evolution. Robustness of ancestral inference is often assumed among closely related species, but tests of this assumption have been limited. Here, we examine the performance of inference methods for data simulated under scenarios of codon bias evolution within the Drosophila melanogaster subgroup. Genome sequence data for multiple, closely related species within this subgroup make it an important system for studying molecular evolutionary genetics. The effects of asymmetric and lineage-specific substitution rates (i.e., varying levels of codon usage bias and departures from equilibrium) on the reliability of ancestral codon usage was investigated. Maximum parsimony inference, which has been widely employed in analyses of Drosophila codon bias evolution, was compared to an approach that attempts to account for uncertainty in ancestral inference by weighting ancestral reconstructions by their posterior probabilities. The latter approach employs maximum likelihood estimation of rate and base composition parameters. For equilibrium and most non-equilibrium scenarios that were investigated, the probabilistic method appears to generate reliable ancestral codon bias inferences for molecular evolutionary studies within the D. melanogaster subgroup. These reconstructions are more reliable than parsimony inference, especially when codon usage is strongly skewed. However, inference biases are considerable for both methods under particular departures from stationarity (i.e., when adaptive evolution is prevalent). Reliability of inference can be sensitive to branch lengths, asymmetry in substitution rates, and the locations and nature of lineage-specific processes within a gene tree. Inference reliability, even among closely related species, can be strongly affected by (potentially unknown) patterns of molecular evolution in lineages ancestral to those of interest.     

Predominant gain of promoter TATA box after gene duplication associated with stress responses

TATA box, the core promoter element, exists in a broad range of eukaryotes, and the expression of TATA-containing genes usually responds to various environmental stresses. Hence, the evolution of TATA-box in duplicate genes may provide some clues for the interrelationship among environmental stress, expression differentiation, and duplicate gene preservation. In the present study, we observed that the TATA box is significantly overrepresented in duplicate genes compared with singletons in human, worm, Arabidopsis, and yeast genomes. We then conducted an extensive functional genomic analysis to investigate the evolution of TATA box along over 700 yeast gene family phylogenies. After reconstructing the ancestral TATA-box states (presence or absence), we found that significantly higher numbers of TATA box gain events than loss events had occurred after yeast gene duplications-the overall gain-loss ratio is about 3-4 to 1. Interestingly, these TATA-gain duplicate genes on average have experienced greater expression divergence from the ancestral expression states than their most closely related TATA-less duplicate partners, but only under environmental stress conditions (asymmetric evolution); indeed, under normal physiological conditions, they have similar expression divergence (symmetric evolution). Moreover, we showed that TATA-gain duplicates are enriched in stress-associated functional categories but that is not the case for TATA-ancestral duplicates (those inherited from their ancestors prior to duplication). Together, we conclude that after the gene duplication, gain of the TATA box in duplicate promoters may have played an important role in yeast duplicate preservation by accelerating expression divergence that may facilitate the adaptive evolution of the organism in response to environmental changes.     

CpG site degeneration triggered by the loss of functional constraint created a highly polymorphic macaque drug-metabolizing gene, CYP1A2

Background

Duplicated genes frequently experience asymmetric rates of sequence evolution. Relaxed selective constraints and positive selection have both been invoked to explain the observation that one paralog within a gene-duplicate pair exhibits an accelerated rate of sequence evolution. In the majority of studies where asymmetric divergence has been established, there is no indication as to which gene copy, ancestral or derived, is evolving more rapidly. In this study we investigated the effect of local synteny (gene-neighborhood conservation) and codon usage on the sequence evolution of gene duplicates in the S. cerevisiae genome. We further distinguish the gene duplicates into those that originated from a whole-genome duplication (WGD) event (ohnologs) versus small-scale duplications (SSD) to determine if there exist any differences in their patterns of sequence evolution.

Results

For SSD pairs, the derived copy evolves faster than the ancestral copy. However, there is no relationship between rate asymmetry and synteny conservation (ancestral-like versus derived-like) in ohnologs. mRNA abundance and optimal codon usage as measured by the CAI is lower in the derived SSD copies relative to ancestral paralogs. Moreover, in the case of ohnologs, the faster-evolving copy has lower CAI and lowered expression.

Conclusions

Together, these results suggest that relaxation of selection for codon usage and gene expression contribute to rate asymmetry in the evolution of duplicated genes and that in SSD pairs, the relaxation of selection stems from the loss of ancestral regulatory information in the derived copy.     

Evolutionary constraint in flanking regions of avian genes

An important comprehension from comparative genomic analysis is that sequence conservation beyond neutral expectations is frequently found outside protein-coding regions, indicating important functional roles of noncoding DNA. Understanding the causes of constraint on noncoding sequence evolution forms an important area of research, not least in light of the importance for understanding the evolution of gene expression. We aligned all orthologous genes of chicken and zebra finch together with 5 kb of their upstream and downstream noncoding sequences, to study the evolution of gene flanking sequences in the avian genome. Using ancestral repeats as a neutral reference, we detected significant evolutionary constraint in the 3' flanking region, highest directly after termination (60%) and then gradually decreasing to about 20% 5 kb downstream. Constraint was higher in annotated 3' untranslated regions (UTRs) than in non-UTRs at the same distance from the stop codon and higher in sequences annotated as microRNA (miRNA)-binding sites than in non-miRNA-binding sites within 3' UTRs. Constraint was also higher when estimated for a smaller data set of genes from more closely related songbird species, indicating turnover of functional elements during avian evolution. On the 5' flanking side constraint was readily seen within the first 125 bp immediately upstream of the start codon (34%) and was about 10% for remaining sequence within 5 kb upstream. Analysis of chicken polymorphism data gave further support for the highest constraint directly before and after the translated region. Finally, we found that genes evolving under the highest constraint measured by d(N)/d(S) also had the highest level of constraint in the 3' flanking region. This study broadens the insights into gene flanking sequence evolution by adding new findings from a vertebrate lineage other than mammals.     

Enhancer evolution in chordates: Lessons from functional analyses of cephalochordate cis-regulatory modules

Chordates comprise three major groups, cephalochordates (amphioxus), tunicates (urochordates), and vertebrates. Since cephalochordates were the early branching group, comparisons between amphioxus and other chordates help us to speculate about ancestral chordates. Here, I summarize accumulating data from functional studies analyzing amphioxus cis-regulatory modules (CRMs) in model systems of other chordate groups, such as mice, chickens, clawed frogs, fish, and ascidians. Conservatism and variability of CRM functions illustrate how gene regulatory networks have evolved in chordates. Amphioxus CRMs, which correspond to CRMs deeply conserved among animal phyla, govern reporter gene expression in conserved expression domains of the putative target gene in host animals. In addition, some CRMs located in similar genomic regions (intron, upstream, or downstream) also possess conserved activity, even though their sequences are divergent. These conservative CRM functions imply ancestral genomic structures and gene regulatory networks in chordates. However, interestingly, if expression patterns of amphioxus genes do not correspond to those of orthologs of experimental models, some amphioxus CRMs recapitulate expression patterns of amphioxus genes, but not those of endogenous genes, suggesting that these amphioxus CRMs are close to the ancestral states of chordate CRMs, while vertebrates/tunicates innovated new CRMs to reconstruct gene regulatory networks subsequent to the divergence of the cephalochordates. Alternatively, amphioxus CRMs may have secondarily lost ancestral CRM activity and evolved independently. These data help to solve fundamental questions of chordate evolution, such as neural crest cells, placodes, a forebrain/midbrain, and genome duplication. Experimental validation is crucial to verify CRM functions and evolution.     

Duplicated genes evolve slower than singletons despite the initial rate increase

Background

Gene duplication is an important mechanism that can lead to the emergence of new functions during evolution. The impact of duplication on the mode of gene evolution has been the subject of several theoretical and empirical comparative-genomic studies. It has been shown that, shortly after the duplication, genes seem to experience a considerable relaxation of purifying selection.

Results

Here we demonstrate two opposite effects of gene duplication on evolutionary rates. Sequence comparisons between paralogs show that, in accord with previous observations, a substantial acceleration in the evolution of paralogs occurs after duplication, presumably due to relaxation of purifying selection. The effect of gene duplication on evolutionary rate was also assessed by sequence comparison between orthologs that have paralogs (duplicates) and those that do not (singletons). It is shown that, in eukaryotes, duplicates, on average, evolve significantly slower than singletons. Eukaryotic ortholog evolutionary rates for duplicates are also negatively correlated with the number of paralogs per gene and the strength of selection between paralogs. A tally of annotated gene functions shows that duplicates tend to be enriched for proteins with known functions, particularly those involved in signaling and related cellular processes; by contrast, singletons include an over-abundance of poorly characterized proteins.

Conclusions

These results suggest that whether or not a gene duplicate is retained by selection depends critically on the pre-existing functional utility of the protein encoded by the ancestral singleton. Duplicates of genes of a higher biological import, which are subject to strong functional constraints on the sequence, are retained relatively more often. Thus, the evolutionary trajectory of duplicated genes appears to be determined by two opposing trends, namely, the post-duplication rate acceleration and the generally slow evolutionary rate owing to the high level of functional constraints.

     

Gene ordering in partitive clustering using microarray expressions

A central step in the analysis of gene expression data is the identification of groups of genes that exhibit similar expression patterns. Clustering and ordering the genes using gene expression data into homogeneous groups was shown to be useful in functional annotation, tissue classification, regulatory motif identification, and other applications. Although there is a rich literature on gene ordering in hierarchical clustering framework for gene expression analysis, there is no work addressing and evaluating the importance of gene ordering in partitive clustering framework, to the best knowledge of the authors. Outside the framework of hierarchical clustering, different gene ordering algorithms are applied on the whole data set, and the domain of partitive clustering is still unexplored with gene ordering approaches. A new hybrid method is proposed for ordering genes in each of the clusters obtained from partitive clustering solution, using microarray gene expressions.Two existing algorithms for optimally ordering cities in travelling salesman problem (TSP), namely, FRAG_GALK and Concorde, are hybridized individually with self organizing MAP to show the importance of gene ordering in partitive clustering framework. We validated our hybrid approach using yeast and fibroblast data and showed that our approach improves the result quality of partitive clustering solution, by identifying subclusters within big clusters, grouping functionally correlated genes within clusters, minimization of summation of gene expression distances, and the maximization of biological gene ordering using MIPS categorization. Moreover, the new hybrid approach, finds comparable or sometimes superior biological gene order in less computation time than those obtained by optimal leaf ordering in hierarchical clustering solution.     

Gene duplication, the evolution of novel gene functions, and detecting functional divergence of duplicates in silico

Duplication of genes increases the amount of genetic material on which evolution can work and has been considered of major importance for the development of biological novelties or to explain important transitions that have occurred during biological evolution. Recently, much research has been devoted to the study of the evolutionary and functional divergence of duplicated genes. Since the majority of genes are part of gene families, there is considerable interest in predicting differences in function between duplicates and assessing the functional redundancy of genes within gene families. In this review, we discuss the strengths and limitations of both older and novel approaches to investigate the evolution of duplicated genes in silico.     

Gene duplication is an evolutionary mechanism for expanding spectral diversity in the long-wavelength photopigments of butterflies

Butterfly long-wavelength (L) photopigments are interesting for comparative studies of adaptive evolution because of the tremendous phenotypic variation that exists in their wavelength of peak absorbance (lambda(max) value). Here we present a comprehensive survey of L photopigment variation by measuring lambda(max) in 12 nymphalid and 1 riodinid species using epi-microspectrophotometry. Together with previous data, we find that L photopigment lambda(max) varies from 510-565 nm in 22 nymphalids, with an even broader 505- to 600-nm range in riodinids. We then surveyed the L opsin genes for which lambda(max) values are available as well as from related taxa and found 2 instances of L opsin gene duplication within nymphalids, in Hermeuptychia hermes and Amathusia phidippus, and 1 instance within riodinids, in the metalmark butterfly Apodemia mormo. Using maximum parsimony and maximum likelihood ancestral state reconstructions to map the evolution of spectral shifts within the L photopigments of nymphalids, we estimate the ancestral pigment had a lambda(max) = 540 nm +/- 10 nm standard error and that blueshifts in wavelength have occurred at least 4 times within the family. We used ancestral state reconstructions to investigate the importance of several amino acid substitutions (Ile17Met, Ala64Ser, Asn70Ser, and Ser137Ala) previously shown to have evolved under positive selection that are correlated with blue spectral shifts. These reconstructions suggest that the Ala64Ser substitution has indeed occurred along the newly identified blueshifted L photopigment lineages. Substitutions at the other 3 sites may also be involved in the functional diversification of L photopigments. Our data strongly suggest that there are limits to the evolution of L photopigment spectral shifts among species with only one L opsin gene and that opsin gene duplication broadens the potential range of lambda(max) values.     

Drosophila duplicate genes evolve new functions on the fly

Gene duplication is thought to play a key role in phenotypic innovation. While several processes have been hypothesized to drive the retention and functional evolution of duplicate genes, their genomic contributions have never been determined. We recently developed the first genome-wide method to classify these processes by comparing distances between expression profiles of duplicate genes and their ancestral single-copy orthologs. Application of our approach to spatial gene expression profiles in two Drosophila species revealed that a majority of young duplicate genes possess new functions, and that new functions are acquired rapidly—often within a few million years. Surprisingly, new functions tend to arise in younger copies of duplicate gene pairs. Moreover, we found that young duplicates are often specifically expressed in testes, whereas old duplicates are broadly expressed across several tissues, providing strong support for the hypothetical “out-of-testes” origin of new genes. In this Extra View, I discuss our findings in the context of theoretical predictions about gene duplication, with a particular emphasis on the importance of natural selection in the evolution of novel phenotypes.     

Expression patterns of the rogue Hox genes Hox3/zen and fushi tarazu in the apterygote insect Thermobia domestica

Many embryonic patterning genes are remarkably conserved between vertebrates and invertebrates, and the Hox genes are paradigmatic examples of this conservation. Yet even Hox genes can change dramatically in evolution. Two genes in particular--Hox3 and fushi tarazu--lost their ancestral roles as homeotic genes and play very different developmental roles in the fruit fly Drosophila melanogaster. The Drosophila Hox3 homologs zerknullt and bicoid act in extraembryonic tissues and in establishment of the anteroposterior axis, respectively, whereas fushi tarazu acts in segmentation and neurogenesis. It would be valuable to know what mechanisms allowed Hox3 and ftz to abandon their ancestral roles as homeotic genes and take on new roles. To explore the evolutionary transition of these genes, we analyzed their expression in a primitive insect, the firebrat Thermobia domestica. The expression patterns seem to represent a stage of evolution intermediate between the ancestral state seen in basal arthropods and the derived expression patterns in Drosophila. These expression data help us to narrow the period in which the gene transitions took place. Hox3 appears to have evolved directly into zen within the insects, whereas ftz seems to have adopted the expression patterns of a segmentation and neurogenesis gene earlier in the mandibulate arthropods.     

Bayesian estimation of ancestral character states on phylogenies

Biologists frequently attempt to infer the character states at ancestral nodes of a phylogeny from the distribution of traits observed in contemporary organisms. Because phylogenies are normally inferences from data, it is desirable to account for the uncertainty in estimates of the tree and its branch lengths when making inferences about ancestral states or other comparative parameters. Here we present a general Bayesian approach for testing comparative hypotheses across statistically justified samples of phylogenies, focusing on the specific issue of reconstructing ancestral states. The method uses Markov chain Monte Carlo techniques for sampling phylogenetic trees and for investigating the parameters of a statistical model of trait evolution. We describe how to combine information about the uncertainty of the phylogeny with uncertainty in the estimate of the ancestral state. Our approach does not constrain the sample of trees only to those that contain the ancestral node or nodes of interest, and we show how to reconstruct ancestral states of uncertain nodes using a most-recent-common-ancestor approach. We illustrate the methods with data on ribonuclease evolution in the Artiodactyla. Software implementing the methods (BayesMultiState) is available from the authors.     

LIKELIHOOD OF ANCESTOR STATES IN ADAPTIVE RADIATION

Theories of ecological diversification make predictions about the timing and ordering of character state changes through history. These theories are testable by “reconstructing” ancestor states using phylogenetic trees and measurements of contemporary species. Here we use maximum likelihood to estimate and evaluate the accuracy of ancestor reconstructions. We present likelihoods of discrete ancestor states and derive probability distributions for continuous ancestral traits. The methods are applied to several examples: diets of ancestral Darwin's finches; origin of inquilinism in gall wasps; microhabitat partitioning and body size evolution in scrubwrens; digestive enzyme evolution in artiodactyl mammals; origin of a sexually selected male trait, the sword, in platies and swordtails; and evolution of specialization in Anolis lizards. When changes between discrete character states are rare, the maximum-likelihood results are similar to parsimony estimates. In this case the accuracy of estimates is often high, with the exception of some nodes deep in the tree. If change is frequent then reconstructions are highly uncertain, especially of distant ancestors. Ancestor states for continuous traits are typically highly uncertain. We conclude that measures of uncertainty are useful and should always be provided, despite simplistic assumptions about the probabilistic models that underlie them. If uncertainty is too high, reconstruction should be abandoned in favor of approaches that fit different models of trait evolution to species data and phylogenetic trees, taking into account the range of ancestor states permitted by the data.