On the evolution of scale-free topologies with a gene regulatory network model

A novel approach to generating scale-free network topologies is introduced, based on an existing artificial gene regulatory network model. From this model, different interaction networks can be extracted, based on an activation threshold. By using an evolutionary computation approach, the model is allowed to evolve, in order to reach specific network statistical measures. The results obtained show that, when the model uses a duplication and divergence initialisation, such as seen in nature, the resulting regulation networks not only are closer in topology to scale-free networks, but also require only a few evolutionary cycles to achieve a satisfactory error value.     

The evolutionary potential of the Drosophila sex determination gene network

The evolution of sex determination mechanisms is known to be relatively rapid, though recent evidence indicates that certain parts of the mechanism may be more highly conserved. These characteristics establish the sex determination mechanism as a good candidate for the theoretical study of gene network evolution, particularly of networks involved in development. We investigate the short-term evolutionary potential of the sex determination mechanism in Drosophila melanogaster with the aid of a synchronous logical model. We introduce general theoretical concepts such as a network-specific form of mutation, and a notion of functional equivalence between networks. We apply this theoretical framework to the sex determination mechanism and compare it to a population of random networks, enabling us to find features both general to sex determination networks, and particular to the Drosophila network. In general, sex determination networks exist within large sets of functionally equivalent networks all of which satisfy the sex determination task. These large sets are in turn composed of subsets which are mutationally related, suggesting a high degree of flexibility is available without compromising the core functionality. Two particular characteristics of the Drosophila network are found: (a) a parsimonious use of gene interactions, and (b) the network structure can produce a relatively large number of dynamical pattern variations through single network mutations.     

Expression diversity and evolutionary dynamics of rice duplicate genes

Duplicate genes are believed to be a major source of new gene functions over evolutionary time. In order to evaluate the evolutionary dynamics of rice duplicate genes, formed principally by paleoployploidization prior to the speciation of the Poaceae family, we have employed a public microarray dataset including 155 gene expression omnibus sample plates and bioinformatics tools. At least 57.4% of old ~70 million years ago (MYA) duplicate gene pairs exhibit divergences in expression over the given experimental set, whereas at least 50.9% of young ~7.7-MYA duplicate gene pairs were shown to be divergent. When grouping the rice duplicate genes according to functional categories, we noted a striking and significant enrichment of divergent duplicate metabolism-associated genes, as compared to that observed in non-divergent duplicate genes. While both non-synonymous substitution (Ka) and synonymous substitution (Ks) values between non- and divergent duplicate gene pairs evidenced significant differences, the Ka/Ks values between them exhibited no significant differences. Interestingly, the average numbers of conserved motifs of the duplicate gene pairs revealed a pattern of decline along with an increase in expression diversity, partially supporting the subfunctionalization model with degenerative complementation in regulatory motifs. Duplicate gene pairs with high local similarity (HLS) segments, which might be formed via conversion between rice paleologs, evidenced higher expression correlations than were observed in the gene pairs without the HLS segments; this probably resulted in an increased likelihood of gene conversion in promoters of the gene pairs harboring HLS segments. More than 60% of the rice gene families exhibited similar high expression diversity between members as compared to that of randomly selected gene pairs. These findings are likely reflective of the evolutionary dynamics of rice duplicate genes for gene retention. An erratum to this article can be found at     

Evolutionary dynamics of duplicated genes in plants

Gene duplication, arising from region-specific duplication or genome-wide polyploidization, is a prominent feature in plant genome evolution. Understanding the mechanisms generating duplicate gene copies and the subsequent dynamics among gene duplicates is vital because these investigations shed light on regional and genome-wide aspects of evolutionary forces shaping intra- and interspecific genome contents, evolutionary relationships, and interactions. This review discusses recent gene duplication analyses in plants, focusing on the molecular and evolutionary dynamics occurring at three different timescales following duplication: (1). initial establishment and persistence of cytotypes, (2). interactions among duplicate gene copies, and (3). longer term differentiation between duplicated genes. These relative time points are presented in terms of their potential adaptive significance and impact on plant evolutionary genomics research.     

Comparative analysis of essential genes and nonessential genes in Escherichia coli K12

Genes can be classified as essential or nonessential based on their indispensability for a living organism. Previous researches have suggested that essential genes evolve more slowly than nonessential genes and the impact of gene dispensability on a gene’s evolutionary rate is not as strong as expected. However, findings have not been consistent and evidence is controversial regarding the relationship between the gene indispensability and the rate of gene evolution. Understanding how different classes of genes evolve is essential for a full understanding of evolutionary biology, and may have medical relevance in the design of new antibacterial agents. We therefore performed an investigation into the properties of essential and nonessential genes. Analysis of evolutionary conservation, protein length distribution and amino acid usage between essential and nonessential genes in Escherichia coli K12 demonstrated that essential genes are relatively preserved throughout the bacterial kingdom when compared to nonessential genes. Furthermore, results show that essential genes, compared to nonessential genes, have a significantly higher proportion of large (>534 amino acids) and small proteins (<139 amino acids) relative to medium-sized proteins. The pattern of amino acids usage shows a similar trend for essential and nonessential genes, although some notable exceptions are observed. These findings help to clarify our understanding of the evolutionary mechanisms of essential and nonessential genes, relevant to the study of mutagenesis and possibly allowing prediction of gene properties in other poorly understood organisms.     

When does gene flow facilitate evolutionary rescue?

Experimental and theoretical studies have highlighted the impact of gene flow on the probability of evolutionary rescue in structured habitats. Mathematical modeling and simulations of evolutionary rescue in spatially or otherwise structured populations showed that intermediate migration rates can often maximize the probability of rescue in gradually or abruptly deteriorating habitats. These theoretical results corroborate the positive effect of gene flow on evolutionary rescue that has been identified in experimental yeast populations. The observations that gene flow can facilitate adaptation are in seeming conflict with traditional population genetics results that show that gene flow usually hampers (local) adaptation. Identifying conditions for when gene flow facilitates survival chances of populations rather than reducing them remains a key unresolved theoretical question. We here present a simple analytically tractable model for evolutionary rescue in a two-deme model with gene flow. Our main result is a simple condition for when migration facilitates evolutionary rescue, as opposed as no migration. We further investigate the roles of asymmetries in gene flow and/or carrying capacities, and the effects of density regulation and local growth rates on evolutionary rescue.     

An ancient frame-shifting event in the highly conserved KPNA gene family has undergone extensive compensation by natural selection in vertebrates

One of the prevailing arguments in evolutionary theory is that the duplicates of genes can acquire novel functionality. This is because only one of the paralogs need maintain the ancestral function, leaving room for natural experimentation due to a respite in purifying selection. Although many duplicates can subsequently become disabled by nullifying mutations, a few may also go on to diverge along a novel evolutionary trajectory. Here, evidence is provided that demonstrates how this scenario may not always be true. Rather, in the case of the highly conserved KPNA importin family, an initial relaxation in selection induced a frameshift that was later suppressed and heavily compensated for as part of a reparative and optimizing process. Despite a resulting divergence, there remains a distinct preservation of both sequence and functionality among the paralogs. This would indicate that duplicates can be retained by selection for reasons related to their redundant functionality. It also shows that, even when positive selection is inferred in duplicate genes, this may be of a compensatory nature rather than one representing any biochemical innovation. Generally, this development would perhaps be a more common outcome for gene duplication than is currently maintained.     

Accelerated Evolutionary Rate May Be Responsible for the Emergence of Lineage-Specific Genes in Ascomycota

The evolutionary origin of “orphan” genes, genes that lack sequence similarity to any known gene, remains a mystery. One suggestion has been that most orphan genes evolve rapidly so that similarity to other genes cannot be traced after a certain evolutionary distance. This can be tested by examining the divergence rates of genes with different degrees of lineage specificity. Here the lineage specificity (LS) of a gene describes the phylogenetic distribution of that gene’s orthologues in related species. Highly lineage-specific genes will be distributed in fewer species in a phylogeny. In this study, we have used the complete genomes of seven ascomycotan fungi and two animals to define several levels of LS, such as Eukaryotes-core, Ascomycota-core, Euascomycetes-specific, Hemiascomycetes-specific, Aspergillus-specific, and Saccharomyces-specific. We compare the rates of gene evolution in groups of higher LS to those in groups with lower LS. Molecular evolutionary analyses indicate an increase in nonsynonymous nucleotide substitution rates in genes with higher LS. Several analyses suggest that LS is correlated with the evolutionary rate of the gene. This correlation is stronger than those of a number of other factors that have been proposed as predictors of a gene’s evolutionary rate, including the expression level of genes, gene essentiality or dispensability, and the number of protein-protein interactions. The accelerated evolutionary rates of genes with higher LS may reflect the influence of selection and adaptive divergence during the emergence of orphan genes. These analyses suggest that accelerated rates of gene evolution may be responsible for the emergence of apparently orphan genes. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Martin Kreitman]     

Gene loss and evolutionary rates following whole-genome duplication in teleost fishes

Teleost fishes provide the first unambiguous support for ancient whole-genome duplication in an animal lineage. Studies in yeast or plants have shown that the effects of such duplications can be mediated by a complex pattern of gene retention and changes in evolutionary pressure. To explore such patterns in fishes, we have determined by phylogenetic analysis the evolutionary origin of 675 Tetraodon duplicated genes assigned to chromosomes, using additional data from other species of actinopterygian fishes. The subset of genes, which was retained in double after the genome duplication, is enriched in development, signaling, behavior, and regulation functional categories. The evolutionary rate of duplicate fish genes appears to be determined by 3 forces: 1) fish proteins evolve faster than mammalian orthologs; 2) the genes kept in double after genome duplication represent the subset under strongest purifying selection; and 3) following duplication, there is an asymmetric acceleration of evolutionary rate in one of the paralogs. These results show that similar mechanisms are at work in fishes as in yeast or plants and provide a framework for future investigation of the consequences of duplication in fishes and other animals.     

以发育模块重复征用和整合为基础的进化创新机制

进化新征的起源和分化是进化发育生物学研究的核心问题。通过对多细胞生物早期发育调控机制的比较分析,发现亲缘关系较远的生物所共有的一些形态特征受保守的发育调控程序调节（深同源性）。许多创新性状的发生是基于对预先存在的基因或发育调控模块的重复利用和整合。发育基因调控网络在结构和功能上高度模块化,因此不仅可以通过模块拆分和重复征用改变发育程式,而且也增强了调控网络自身的进化力。研究基因调控网络和发育系统的进化动态将有助于更深入地认识生物演化过程中创新性状发生和表型进化的分子机制。     

Coordination of cell differentiation and migration in mathematical models of caudal embryonic axis extension

Vertebrate embryos display a predominant head-to-tail body axis whose formation is associated with the progressive development of post-cranial structures from a pool of caudal undifferentiated cells. This involves the maintenance of active FGF signaling in this caudal region as a consequence of the restricted production of the secreted factor FGF8. FGF8 is transcribed specifically in the caudal precursor region and is down-regulated as cells differentiate and the embryo extends caudally. We are interested in understanding the progressive down-regulation of FGF8 and its coordination with the caudal movement of cells which is also known to be FGF-signaling dependent. Our study is performed using mathematical modeling and computer simulations. We use an individual-based hybrid model as well as a caricature continuous model for the simulation of experimental observations (ours and those known from the literature) in order to examine possible mechanisms that drive differentiation and cell movement during the axis elongation. Using these models we have identified a possible gene regulatory network involving self-repression of a caudal morphogen coupled to directional domain movement that may account for progressive down-regulation of FGF8 and conservation of the FGF8 domain of expression. Furthermore, we have shown that chemotaxis driven by molecules, such as FGF8 secreted in the stem zone, could underlie the migration of the caudal precursor zone and, therefore, embryonic axis extension. These mechanisms may also be at play in other developmental processes displaying a similar mode of axis extension coupled to cell differentiation.     

Of chicken wings and frog legs: a smorgasbord of evolutionary variation in mechanisms of tetrapod limb development

The tetrapod limb, which has served as a paradigm for the study of development and morphological evolution, is becoming a paradigm for developmental evolution as well. In its origin and diversification, the tetrapod limb has undergone a great deal of remodeling. These morphological changes and other evolutionary phenomena have produced variation in mechanisms of tetrapod limb development. Here, we review that variation in the four major clades of limbed tetrapods. Comparisons in a phylogenetic context reveal details of development and evolution that otherwise may have been unclear. Such details include apparent differences in the mechanisms of dorsal-ventral patterning and limb identity specification between mouse and chick and mechanistic novelties in amniotes, anurans, and urodeles. As we gain a better understanding of the details of limb development, further differences among taxa will be revealed. The use of appropriate comparative techniques in a phylogenetic context thus sheds light on evolutionary transitions in limb morphology and the generality of developmental models across species and is therefore important to both evolutionary and developmental biologists.     

Transitions and transversions in evolutionary descent: An approach to understanding

Summary In this paper I lay a quantitative theoretical groundwork for understanding the proportions of the possible types of base substitutions observed between 12 genes sharing a common ancestor and isolated from extant species. The experimentally observed types of base substitution between two sequenced genes do not give a direct measure of the types of base substitutions that occur during evolutionary descent. However, by use of a statistical assemblage of these observations, we can recover, without the assumption of parsimony, the conditional base substitution probabilities that determine this descent. Three methods—direct count, regression, and informational entropy maximization—are described by which these probabilities can be estimated from experimental data. The methods are complementary in that each is most useful for somewhat different types of experimental data. These methods are used to study the ratio of transversions to transitions during gene divergence. Though this ratio is not constant during divergence, it does approach a stable limiting value that in principle can vary from zero, corresponding to 100% transition differences, to infinity, corresponding to 0% transition differences. In practice the limiting ratio tends to hover around a value of two, which is expected on a random basis. However, base substitution pathways that are very nonrandom also may lead to a limiting ratio of exactly two, so that such a value is not diagnostic for random pathways. The limiting ratio can be directly calculated from a knowledge of the twelve conditional probabilities for each type of base substitution, or from a knowledge of the equilibrium base composition of the DNAs compared. An expression is given for this calculation. Fifteen years ago Jean Derancourt, Andrew Lebor and Emile Zuckerkandl (1967), analyzing the amino acid sequence of globin chains coded by nuclear genes, made the original observation that the proportion of transition differences decreases with increasing evolutionary time. Recently Brown et al. (1982) and Brown and Simpson (1982) have reported a decrease in the observed proportion of transition differences in mitochondrial DNA with increasing evolutionary divergence. The conditions that must be satisfied for this type of behavior to occur at stable base composition and with stable base substitution probabilities are defined. Multiple substitutionsper se do not lead to a decrease in transition differences with increasing evolutionary divergence.     

Evolving homeostatic tissue using genetic algorithms

Multicellular organisms maintain form and function through a multitude of homeostatic mechanisms. The details of these mechanisms are in many cases unknown, and so are their evolutionary origin and their link to development. In order to illuminate these issues we have investigated the evolution of structural homeostasis in the simplest of cases, a tissue formed by a mono-layer of cells. To this end, we made use of a 3-dimensional hybrid cellular automaton, an individual-based model in which the behaviour of each cell depends on its local environment. Using an evolutionary algorithm (EA) we evolved cell signalling networks, both with a fixed and an incremental fitness evaluation, which give rise to and maintain a mono-layer tissue structure. Analysis of the solutions provided by the EA shows that the two evaluation methods gives rise to different types of solutions to the problem of homeostasis. The fixed method leads to almost optimal solutions, where the tissue relies on a high rate of cell turnover, while the solutions from the incremental scheme behave in a more conservative manner, only dividing when necessary. In order to test the robustness of the solutions we subjected them to environmental stress, by wounding the tissue, and to genetic stress, by introducing mutations. The results show that the robustness very much depends on the mechanism responsible for maintaining homeostasis. The two evolved cell types analysed present contrasting mechanisms by which tissue homeostasis can be maintained. This compares well to different tissue types found in multicellular organisms. For example the epithelial cells lining the colon in humans are shed at a considerable rate, while in other tissue types, which are not as exposed, the conservative type of homeostatic mechanism is normally found. These results will hopefully shed light on how multicellular organisms have evolved homeostatic mechanisms and what might occur when these mechanisms fail, as in the case of cancer.     

Phylogeny, taxonomy, and evolution of the endothelin receptor gene family

A gene phylogeny provides the natural historical order to classify genes and to understand their functional, structural, and genomic diversity. The gene family of endothelin receptors (EDNR) is responsible for many key physiological and developmental processes of tetrapods and teleosts. This study provides a well-defined gene phylogeny for the EDNR family, which is used to classify its members and to assess their evolution. The EDNR phylogeny supports the recognition of the EDNRA, EDNRB, and EDNRC subfamilies, as well as more lineage-specific duplicates of teleosts and the African clawed frog. The duplications for these nominal genes are related to the various whole-genome amplifications of vertebrates, jawed vertebrates, fishes, and frog. The EDNR phylogeny also identifies several gene losses, including that of EDNRC from placental and marsupial (therian) mammals. When coupled with structural and biochemical information, site-specific analyses of evolutionary rate shifts reveal two distinct patterns of potential functional changes at the sequence level between therian versus non-therian EDNRA and EDNRB (i.e., between groups without and with EDNRC). An analysis of linkage maps and tetrapod synteny further suggests that the loss of therian EDNRC may be related to a chromosomal deletion in its common ancestor.