A framework for evolutionary systems biology

Background  

Many difficult problems in evolutionary genomics are related to mutations that have weak effects on fitness, as the consequences of mutations with large effects are often simple to predict. Current systems biology has accumulated much data on mutations with large effects and can predict the properties of knockout mutants in some systems. However experimental methods are too insensitive to observe small effects.     

The evolutionary dynamics of self-incompatibility systems

Self-incompatible flowering plants reject pollen that expresses the same mating specificity as the pistil (female reproductive tract). In most plant families, pollen and pistil mating specificities segregate as a single locus, the S locus. In at least two self-incompatibility systems, distinct pollen and pistil specificity genes are embedded in an extensive nonrecombining tract. To facilitate consideration of how new S locus specificities arise in systems with distinct pollen and pistil genes, we present a graphical model for the generation of hypotheses. It incorporates the evolutionary principle that nonreciprocal siring success (cross-pollinations between two plants produce seeds in only one direction) tends to favor the rejecting partner. This model suggests that selection within S-allele specificity classes could accelerate the rate of nonsynonymous (amino acid-changing) substitutions, with periodic selective sweeps removing segregating variation within classes. Accelerated substitution within specificity classes could also promote the origin of new S-allele specificities.     

The evolutionary history of quorum-sensing systems in bacteria

Communication among bacterial cells through quorum-sensing (QS) systems is used to regulate ecologically and medically important traits, including virulence to hosts. QS is widespread in bacteria; it has been demonstrated experimentally in diverse phylogenetic groups, and homologs to the implicated genes have been discovered in a large proportion of sequenced bacterial genomes. The widespread distribution of the underlying gene families (LuxI/R and LuxS) raises the questions of how often QS genes have been transferred among bacterial lineages and the extent to which genes in the same QS system exchange partners or coevolve. Phylogenetic analyses of the relevant gene families show that the genes annotated as LuxI/R inducer and receptor elements comprise two families with virtually no homology between them and with one family restricted to the gamma-Proteobacteria and the other more widely distributed. Within bacterial phyla, trees for the LuxS and the two LuxI/R families show broad agreement with the ribosomal RNA tree, suggesting that these systems have been continually present during the evolution of groups such as the Proteobacteria and the Firmicutes. However, lateral transfer can be inferred for some genes (e.g., from Firmicutes to some distantly related lineages for LuxS). In general, the inducer/receptor elements in the LuxI/R systems have evolved together with little exchange of partners, although loss or replacement of partners has occurred in several lineages of gamma-Proteobacteria, the group for which sampling is most intensive in current databases. For instance, in Pseudomonas aeruginosa, a transferred QS system has been incorporated into the pathway of a native one. Gene phylogenies for the main LuxI/R family in Pseudomonas species imply a complex history of lateral transfer, ancestral duplication, and gene loss within the genus.     

Prediction of functional phosphorylation sites by incorporating evolutionary information

Protein phosphorylation is a ubiquitous protein post-translational modification,which plays an important role in cellular signaling systems underlying various physiological and pathological processes.Current in silico methods mainly focused on the prediction of phosphorylation sites,but rare methods considered whether a phosphorylation site is functional or not.Since functional phosphorylation sites are more valuable for further experimental research and a proportion of phosphorylation sites have no direct functional effects,the prediction of functional phosphorylation sites is quite necessary for this research area.Previous studies have shown that functional phosphorylation sites are more conserved than non-functional phosphorylation sites in evolution.Thus,in our method,we developed a web sewer by integrating existing phosphorylation site prediction methods,as well as both absolute and relative evolutionary conservation scores to predict the most likely functional phosphorylation sites.Using our method,we predicted the most likely functional sites of the human,rat and mouse proteomes and built a database for the predicted sites.By the analysis of overall prediction results,we demonstrated that protein phosphorylation plays an important role in all the enriched KEGG pathways.By the analysis of protein-specific prediction results,we demonstrated the usefulness of our method for individual protein studies.Our method would help to characterize the most likely functional phosphorylation sites for further studies in this research area.     

Amphioxus and tunicates as evolutionary model systems

One important question in evolutionary biology concerns the origin of vertebrates from invertebrates. The current consensus is that the proximate ancestor of vertebrates was an invertebrate chordate. Today, the invertebrate chordates comprise cephalochordates (amphioxus) and tunicates (each a subphylum in the phylum Chordata, which also includes the vertebrate subphylum). It was widely accepted that, within the chordates, tunicates represent the sister group of a clade of cephalochordates plus vertebrates. However, recent studies suggest that the evolutionary positions of tunicates and cephalochordates should be reversed, the implications of which are considered here. We also review the two major groups of invertebrate chordates and compare relative advantages (and disadvantages) of each as model systems for elucidating the origin of the vertebrates.     

Mechanical sensitivity reveals evolutionary dynamics of mechanical systems

A classic question in evolutionary biology is how form–function relationships promote or limit diversification. Mechanical metrics, such as kinematic transmission (KT) in linkage systems, are useful tools for examining the evolution of form and function in a comparative context. The convergence of disparate systems on equivalent metric values (mechanical equivalence) has been highlighted as a source of potential morphological diversity under the assumption that morphology can evolve with minimal impact on function. However, this assumption does not account for mechanical sensitivity—the sensitivity of the metric to morphological changes in individual components of a structure. We examined the diversification of a four-bar linkage system in mantis shrimp (Stomatopoda), and found evidence for both mechanical equivalence and differential mechanical sensitivity. KT exhibited variable correlations with individual linkage components, highlighting the components that influence KT evolution, and the components that are free to evolve independently from KT and thereby contribute to the observed pattern of mechanical equivalence. Determining the mechanical sensitivity in a system leads to a deeper understanding of both functional convergence and morphological diversification. This study illustrates the importance of multi-level analyses in delineating the factors that limit and promote diversification in form–function systems.     

The evolutionary genetics of sexual systems in flowering plants.

Population genetic studies of the evolution of breeding systems in flowering plants are reviewed. The selective advantage of a gene's increasing the selfing rate is stressed. In the evolution of outbreeding mechanisms, some strong disadvantage to selfing must therefore be acting; it is suggested that this disadvantage is inbreeding depression. Populations with no absolute barrier to selfing, and with intermediate levels of self-fertilization, appear to be the most likely starting state for the evolution of outbreeding mechanisms. There is some evidence for inbreeding depression in such populations. The evolution of distyly and dioecy are considered in some detail. An explanation for the existence of supergenes controlling these systems is proposed. The breakdown of distyly and tristyly are also considered. The evolution of recombination rates in selfing and outcrossing species is examined briefly.     

Novel evolutionary relationship among four fish model systems

Knowledge of the correct phylogenetic relationships among animals is crucial for the valid interpretation of evolutionary trends in biology. Zebrafish, medaka, pufferfish and cichilds are fish models for development, genomics and comparative genetics studies, although their phylogenetic relationships have not been tested rigorously. The results of phylogenomic analysis based on 20 nuclear protein-coding genes confirmed the basal placement of zebrafish in the fish phylogeny but revealed an unexpected relationship among the other three species, contrary to traditionally held systematic views based on morphology. Our analyses show that medaka (Beloniformes) and cichlids (Perciformes) appear to be more closely related to each other than either of them is to pufferfish (Tetraodontiformes), suggesting that a re-interpretation of some findings in comparative biology might be required. In addition, phylogenomic analyses show that fish typically have more copies of nuclear genes than land vertebrates, supporting the fish-specific genome duplication hypothesis.     

Plant self-incompatibility systems: a molecular evolutionary perspective

Incompatibility recognition systems preventing self-fertilization have evolved several times in independent lineages of Angiosperm plants, and three main model systems are well characterized at the molecular level [the gametophytic self-incompatibility (SI) systems of Solanaceae, Rosaceae and Anthirrhinum, the very different system of poppy, and the system in Brassicaceae with sporophytic control of pollen SI reactions]. In two of these systems, the genes encoding both components of pollen-pistil recognition are now known, showing clearly that these two proteins are distinct, that is, SI is a lock-and-key mechanism. Here, we review recent findings in the three well-studied systems in the light of these results and analyse their implications for understanding polymorphism and coevolution of the two SI genes, in the context of a tightly linked genome region.     

Engineering microbial systems to explore ecological and evolutionary dynamics

A major goal of biological research is to provide a mechanistic understanding of diverse biological processes. To this end, synthetic biology offers a powerful approach, whereby biological questions can be addressed in a well-defined framework. By constructing simple gene circuits, such studies have generated new insights into the design principles of gene regulatory networks. Recently, this strategy has been applied to analyze ecological and evolutionary questions, where population-level interactions are critical. Here, we highlight recent development of such systems and discuss how they were used to address problems in ecology and evolutionary biology. As illustrated by these examples, synthetic ecosystems provide a unique platform to study ecological and evolutionary phenomena that are challenging to study in their natural contexts.     

Prediction of binding hot spot residues by using structural and evolutionary parameters

In this work, we present a method for predicting hot spot residues by using a set of structural and evolutionary parameters. Unlike previous studies, we use a set of parameters which do not depend on the structure of the protein in complex, so that the predictor can also be used when the interface region is unknown. Despite the fact that no information concerning proteins in complex is used for prediction, the application of the method to a compiled dataset described in the literature achieved a performance of 60.4%, as measured by F-Measure, corresponding to a recall of 78.1% and a precision of 49.5%. This result is higher than those reported by previous studies using the same data set.     

Prediction of effluent pollution level in activated sludge systems

This paper deals with theoretical and practical problems of modeling activated sludge wastewater treatment processes. A substrate removal model suitable for predicting the concentration of poorly biodegradable substances in purified water, based on batch laboratory experiments, is presented.     

Fish species diversity in southern African estuarine systems: an evolutionary perspective

Synopsis The ichthyofauna of southern African estuaries consists primarily of juvenile marine species that use these habitats as nursery areas. The abundance and biomass of fishes in estuarine systems are typically high but species diversity is generally low, with only a few taxa dominating the community. This relatively low species diversity is attributed to the fact that estuaries in the region are unpredictable environments which lack any degree of permanence and are dominated by mobile marine eurytopes. Although stenotopes, represented mainly by small resident species from the families Gobiidae, Clinidae and Syngnathidae, are present in southern African estuaries, little speciation appears to have occurred. A possible reason for this lack of speciation, apart from the seasonal and annual variability of the abiotic environment, is that the lifetime of individual systems is usually limited to a few thousand years. In addition, fishes utilising southern African estuaries need to remain flexible (eurytopic) in their responses to an external environment which is unlikely to become more stable in the future. Thus the lack of permanence and fluctuating nature of southern African estuaries on both a spatial and temporal scale, together with the dominance of eurytopes in these systems, does not favour the evolution of new species. A preliminary examination of the available literature indicates that a detailed review of estuarine ichthyofaunal communities on a global basis will probably mirror the trends outlined above, and reveal a domination of these dynamic ecosystems by eurytopic taxa with low speciation potential.     

Genetic variation inLycopersicon pimpinellifolium: Evidence of evolutionary change in mating systems

L. pimpinellifolium is a highly heterogeneous species, exhibiting pronounced trends from one end of its linear distribution to the other in nearly every studied genetic locus. Drastic differences between populations were also detected in genetic variability and rates of outcrossing. Highly significant positive correlations exist in every possible comparison between flower size, degree of stigma exsertion, heterozygosity, and allelic polymorphism. The hypothesis most compatible with observations proposes that the very uniform, highly self-pollinated biotypes originated from the more primitive, more variable, facultatively allogamous forms.Dedicated with affection, gratitude and respect to the memory ofTheodosius Dobzhansky.     

Physiological and evolutionary studies of NAP systems in Shewanella piezotolerans WP3

Most of the Shewanella species contain two periplasmic nitrate reductases (NAP-α and NAP-β), which is a unique feature of this genus. In the present study, the physiological function and evolutionary relationship of the two NAP systems were studied in the deep-sea bacterium Shewanella piezotolerans WP3. Both of the WP3 nap gene clusters: nap-α (napD1A1B1C) and nap-β (napD2A2B2) were shown to be involved in nitrate respiration. Phylogenetic analyses suggest that NAP-β originated earlier than NAP-α. Tetraheme cytochromes NapC and CymA were found to be the major electron deliver proteins, and CymA also served as a sole electron transporter towards nitrite reductase. Interestingly, a ΔnapA2 mutant with the single functional NAP-α system showed better growth than the wild-type strain, when grown in nitrate medium, and it had a selective advantage to the wild-type strain. On the basis of these results, we proposed the evolution direction of nitrate respiration system in Shewanella: from a single NAP-β to NAP-β and NAP-α both, followed by the evolution to a single NAP-α. Moreover, the data presented here will be very useful for the designed engineering of Shewanella for more efficient respiring capabilities for environmental bioremediation.     

Why do evolutionary systems stick to the edge of chaos

Summary  The long-term behaviour of dynamic systems can be classified in two different regimes, regular or chaotic, depending on the values of the control parameters, which are kept constant during the time evolution. Starting from slightly different initial conditions, a regular system converges to the same final trajectory, whereas a chaotic system follows two distinct trajectories exponentially diverging from each other. In spite of these differences, regular and chaotic systems share a common property: both arrive exponentially fast to their final destiny, becoming trapped there. In both cases one has finite transient times. This is not a profitable property in what concerns evolutionary strategies, where the eternal search for new forms, better than the current one, is imperative. That is why evolutionary dynamic systems tend to tune themselves in very particular situations in between regular and chaotic regimes. These particular situations present eternal transients, and the system actually never reaches its final destiny, preserving diversity. This feature allows the system to visit other regions of the space of possibilities, not only the tiny region covered by its final attractor.     

Driving developmental and evolutionary change: A systems biology view

Embryonic development is underpinned by ～50 core processes that drive morphogenesis, growth, patterning and differentiation, and each is the functional output of a complex molecular network. Processes are thus the natural and parsimonious link between genotype and phenotype and the obvious focus for any discussion of biological change. Here, the implications of this approach are explored. One is that many features of developmental change can be modeled as mathematical graphs, or sets of connected triplets of the general form <noun><verb><noun>. In these, the verbs (edges) are the outputs of the processes that drive change and the nouns (nodes) are the time-dependent states of biological entities (from molecules to tissues). Such graphs help unpick the multi-level complexity of developmental phenomena and may help suggest new experiments. Another comes from analyzing the effect of mutation that lead to tinkering with the dynamic properties of these processes and to congenital abnormalities; if these changes are both inherited and advantageous, they become evolutionary modifications. In this context, protein networks often represents what classical evolutionary genetics sees as genes, and the realization that traits reflect the output processes of complex networks, particularly for growth, patterning and pigmentation, rather than anything simpler clarifies some problems that the evolutionary synthesis of the 1950s has found hard to solve. In the wider context, most processes are used many times in development and cooperate to produce tissue modules (bones, branching duct systems, muscles etc.). Their underlying generative networks can thus be thought of as genomic modules or subroutines.     

Biogeographical and evolutionary importance of the European high mountain systems

Background

Menopause is a seemingly maladaptive life-history trait that is found in many long-lived mammals. There are two competing evolutionary hypotheses for this phenomenon; in the adaptive view of menopause, the cessation of reproduction may increase the fitness of older females; in the non-adaptive view, menopause may be explained by physiological deterioration with age. The decline and eventual cessation of reproduction has been documented in a number of mammalian species, however the evolutionary cause of this trait is unknown.

Results

We examined a unique 30-year time series of killer whales, tracking the reproductive performance of individuals through time. Killer whales are extremely long-lived, and may have the longest documented post-reproductive lifespan of any mammal, including humans. We found no strong support for either of the adaptive hypotheses of menopause; there was little support for the presence of post-reproductive females benefitting their daughter's reproductive performance (interbirth interval and reproductive lifespan of daughters), or the number of mature recruits to the population. Oldest mothers (> 35) did appear to have a small positive impact on calf survival, suggesting that females may gain experience with age. There was mixed support for the grandmother hypothesis – grandoffspring survival probabilities were not influenced by living grandmothers, but grandmothers may positively influence survival of juveniles at a critical life stage.

Conclusion

Although existing data do not allow us to examine evolutionary tradeoffs between survival and reproduction for this species, we were able to examine the effect of maternal age on offspring survival. Our results are consistent with similar studies of other mammals – oldest mothers appear to be better mothers, producing calves with higher survival rates. Studies of juvenile survival in humans have reported positive benefits of grandmothers on newly weaned infants; our results indicate that 3-year old killer whales may experience a positive benefit from helpful grandmothers. While our research provides little support for menopause evolving to provide fitness benefits to mothers or grandmothers, our work supports previous research showing that menopause and long post-reproductive lifespans are not a human phenomenon.