A neutral model of transcriptome evolution

Microarray technologies allow the identification of large numbers of expression differences within and between species. Although environmental and physiological stimuli are clearly responsible for changes in the expression levels of many genes, it is not known whether the majority of changes of gene expression fixed during evolution between species and between various tissues within a species are caused by Darwinian selection or by stochastic processes. We find the following: (1) expression differences between species accumulate approximately linearly with time; (2) gene expression variation among individuals within a species correlates positively with expression divergence between species; (3) rates of expression divergence between species do not differ significantly between intact genes and expressed pseudogenes; (4) expression differences between brain regions within a species have accumulated approximately linearly with time since these regions emerged during evolution. These results suggest that the majority of expression differences observed between species are selectively neutral or nearly neutral and likely to be of little or no functional significance. Therefore, the identification of gene expression differences between species fixed by selection should be based on null hypotheses assuming functional neutrality. Furthermore, it may be possible to apply a molecular clock based on expression differences to infer the evolutionary history of tissues.     

Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome

Despite growing evidence of rapid evolution in protein coding genes, the contribution of positive selection to intra- and interspecific differences in protein coding regions of the genome is unclear. We attempted to see if genes coding for secreted proteins and genes with narrow expression, specifically those preferentially expressed in the mammary gland, have diverged at a faster rate between domestic cattle (Bos taurus) and humans (Homo sapiens) than other genes and whether positive selection is responsible. Using a large data set, we identified groups of genes based on secretion and expression patterns and compared them for the rate of nonsynonymous (dN) and synonymous (dS) substitutions per site and the number of radical (Dr) and conservative (Dc) amino acid substitutions. We found evidence of rapid evolution in genes with narrow expression, especially for those expressed in the liver and mammary gland and for genes coding for secreted proteins. We compared common human polymorphism data with human-cattle divergence and found that genes with high evolutionary rates in human-cattle divergence also had a large number of common human polymorphisms. This argues against positive selection causing rapid divergence in these groups of genes. In most cases dN/dS ratios were lower in human-cattle divergence than in common human polymorphism presumably due to differences in the effectiveness of purifying selection between long-term divergence and short-term polymorphism.     

Statistical properties of neutral evolution

Neutral evolution is the simplest model of molecular evolution and thus it is most amenable to a comprehensive theoretical investigation. In this paper, we characterize the statistical properties of neutral evolution of proteins under the requirement that the native state remains thermodynamically stable, and compare them to the ones of Kimura's model of neutral evolution. Our study is based on the Structurally Constrained Neutral (SCN) model which we recently proposed. We show that, in the SCN model, the substitution rate decreases as longer time intervals are considered. Fluctuations from one branch of the evolutionary tree to another are strong, leading to a non-Poissonian statistics for the substitution process. Such strong fluctuations are in part due to the fact that neutral substitution rates for individual residues are strongly correlated for most residue pairs. Interestingly, structurally conserved residues, characterized by a much below average substitution rate, are also much less correlated to other residues and evolve in a much more regular way. Our results can improve methods aimed at distinguishing between neutral and adaptive substitutions as well as methods for computing the expected number of substitutions occurred since the divergence of two protein sequences. In particular, we compute the minimal sequence similarity below which no information about the evolutionary divergence of the compared sequences can be obtained.     

Thermodynamics of neutral protein evolution

Naturally evolving proteins gradually accumulate mutations while continuing to fold to stable structures. This process of neutral evolution is an important mode of genetic change and forms the basis for the molecular clock. We present a mathematical theory that predicts the number of accumulated mutations, the index of dispersion, and the distribution of stabilities in an evolving protein population from knowledge of the stability effects (delta deltaG values) for single mutations. Our theory quantitatively describes how neutral evolution leads to marginally stable proteins and provides formulas for calculating how fluctuations in stability can overdisperse the molecular clock. It also shows that the structural influences on the rate of sequence evolution observed in earlier simulations can be calculated using just the single-mutation delta deltaG values. We consider both the case when the product of the population size and mutation rate is small and the case when this product is large, and show that in the latter case the proteins evolve excess mutational robustness that is manifested by extra stability and an increase in the rate of sequence evolution. All our theoretical predictions are confirmed by simulations with lattice proteins. Our work provides a mathematical foundation for understanding how protein biophysics shapes the process of evolution.     

Molecular clock in neutral protein evolution

Background

A frequent observation in molecular evolution is that amino-acid substitution rates show an index of dispersion (that is, ratio of variance to mean) substantially larger than one. This observation has been termed the overdispersed molecular clock. On the basis of in silico protein-evolution experiments, Bastolla and coworkers recently proposed an explanation for this observation: Proteins drift in neutral space, and can temporarily get trapped in regions of substantially reduced neutrality. In these regions, substitution rates are suppressed, which results in an overall substitution process that is not Poissonian. However, the simulation method of Bastolla et al. is representative only for cases in which the product of mutation rate μ and population size N_e is small. How the substitution process behaves when μN_e is large is not known.

Results

Here, I study the behavior of the molecular clock in in silico protein evolution as a function of mutation rate and population size. I find that the index of dispersion decays with increasing μN_e, and approaches 1 for large μN_e . This observation can be explained with the selective pressure for mutational robustness, which is effective when μN_e is large. This pressure keeps the population out of low-neutrality traps, and thus steadies the ticking of the molecular clock.

Conclusions

The molecular clock in neutral protein evolution can fall into two distinct regimes, a strongly overdispersed one for small μN_e, and a mostly Poissonian one for large μN_e. The former is relevant for the majority of organisms in the plant and animal kingdom, and the latter may be relevant for RNA viruses.

     

Cytocatalytic evolution in plants

Low fertility and vigor, and requirements of ecological niches distinct from diploids are not universal characteristics of autoploids. InClaytonia, Hedyotis, Oldenlandia, and other genera, occurrence and frequency of both polyploids and aneuploids within species populations suggest a greater role of such mutations in the evolution of vascular plants than heretofore presumed.     

Adaptive evolution on neutral networks

We study the evolution of large but finite asexual populations evolving in fitness landscapes in which all mutations are either neutral or strongly deleterious. We demonstrate that despite the absence of higher fitness genotypes, adaptation takes place as regions with more advantageous distributions of neutral genotypes are discovered. Since these discoveries are typically rare events, the population dynamics can be subdivided into separate epochs, with rapid transitions between them. Within one epoch, the average fitness in the population is approximately constant. The transitions between epochs, however, are generally accompanied by a significant increase in the average fitness. We verify our theoretical considerations with two analytically tractable bitstring models.     

Misconceptions and false expectations in neutral evolution

Neutral evolution results from random recurrent mutation and genetic drift. A small part of random evolution, that which is related to protein or DNA polymorphisms, is the subject of the Neutral Theory of Evolution. One of the foundations of this theory is the demonstration that the mutation rate (m) is equal to the substitution rate. Since both rates are independent of population size, they are independent of drift, which is dependent upon population size. Neutralists have erroneously equated the substitution rate with the fixation rate, despite the fact that they are antithetical conceptions. The neutralists then applied the random walk stochastic model to justify alleles or bases that were fixated or eliminated. In this model, once the allele or base frequencies reach the monomorphic states (values of 1.0 or 0.0), the absorbing barriers, they can no longer return to the polymorphic state. This operates in a pure mathematical model. If recurrent mutation occurs (as in biotic real systems) fixation and elimination are impossible. A population of bacteria in which m = 10(-8) base mutation (or substitution)/site/generation and the reproduction rate is 1000 cell cycle/year should replace all its genome bases in approximately 100,000 years. The expected situation for all sites is polymorphism for the four bases rather than monomorphism at 1.0 or 0.0 frequencies. If fixation and elimination of a base for more than 500,000 years are impossible, then most of the neutral theory is untenable. A new complete neutral model, which allows for recurrent substitutions, is proposed here based on recurrent mutation or substitution and drift alone. The model fits a binomial or Poisson distribution and not a geometric one, as does neutral theory.     

Anagenetic evolution in island plants

Aim  Plants in islands have often evolved through adaptive radiation, providing the classical model of evolution of closely related species each with strikingly different morphological and ecological features and with low levels of genetic divergence. We emphasize the importance of an alternative (anagenetic) model of evolution, whereby a single island endemic evolves from a progenitor and slowly builds up genetic variation through time.
Location  Continental and oceanic islands.
Methods  We surveyed 2640 endemic angiosperm species in 13 island systems of the world, both oceanic and continental, for anagenetic and cladogenetic patterns of speciation. Genetic data were evaluated from a progenitor and derivative species pair in Ullung Island, Korea, and Japan.
Results  We show that the anagenetic model of evolution is much more important in oceanic islands than previously believed, accounting for levels of endemic specific diversity from 7% in the Hawaiian Islands to 88% in Ullung Island, Korea, with a mean for all islands of 25%. Examination of an anagenetically derived endemic species in Ullung Island reveals genetic (amplified fragment length polymorphism) variation equal or nearly equal to that of its continental progenitor.
Main conclusions  We hypothesize that, during anagenetic speciation, initial founder populations proliferate, and then accumulate genetic variation slowly through time by mutation and recombination in a relatively uniform environment, with drift and/or selection yielding genetic and morphological divergence sufficient for the recognition of new species. Low-elevation islands with low habitat heterogeneity are highly correlated with high levels of anagenetic evolution, allowing prediction of levels of the two models of evolution from these data alone. Both anagenetic and adaptive radiation models of speciation are needed to explain the observed levels of specific and genetic diversity in oceanic islands.     

Subgenome evolution in allotetraploid plants

Polyploidization is a well-known speciation and adaptation mechanism. Traces of former polyploidization events were discovered within many genomes, and especially in plants. Allopolyploidization by interspecific hybridization between two species is common. Among hybrid plants, many are domesticated species of agricultural interest and many of their genomes and of their presumptive parents have been sequenced. Hybrid genomes remain challenging to analyse because of the presence of multiple subgenomes. The genomes of hybrids often undergo rearrangement and degradation over time. Based on 10 hybrid plant genomes from six different genera, with hybridization dating from 10,000 to 5 million years ago, we assessed subgenome degradation, subgenomic intermixing and biased subgenome fractionation. The restructuring of hybrid genomes does not proceed proportionally with the age of the hybrid. The oldest hybrids in our data set display completely different fates: whereas the subgenomes of the tobacco plant Nicotiana benthamiana are in an advanced stage of degradation, the subgenomes of quinoa (Chenopodium quinoa) are exceptionally well conserved by structure and sequence. We observed statistically significant biased subgenome fractionation in seven out of 10 hybrids, which had different ages and subgenomic intermixing levels. Hence, we conclude that no correlation exists between biased fractionation and subgenome intermixing. Lastly, domestication may encourage or hinder subgenome intermixing, depending on the evolutionary context. In summary, comparative analysis of hybrid genomes and their presumptive parents allowed us to determine commonalities and differences between their evolutionary fates. In order to facilitate the future analysis of further hybrid genomes, we automated the analysis steps within manticore , which is publicly available at https://github.com/MatteoSchiavinato/manticore.git .     

Exploring phenotype space through neutral evolution

RNA secondary-structure folding algorithms predict the existence of connected networks of RNA sequences with identical secondary structures. Fitness landscapes that are based on the mapping between RNA sequence and RNA secondary structure hence have many neutral paths. A neutral walk on these fitness landscapes gives access to a virtually unlimited number of secondary structures that are a single point mutation from the neutral path. This shows that neutral evolution explores phenotype space and can play a role in adaptation. Received: 23 December 1995 / Accepted: 17 March 1996