On the dynamics of evolutionary discontinuities

Intraspecific competition and selection in polymorphic populations are considered as gradual, yet highly nonlinear, processes which possibly lead to structurally unstable states in species and ecosystems. Discontinuous genetic and morphological changes may thus occur, which are interpreted as evolutionary “catastrophes,” that is, instantaneous transitions between distant, secularly stable evolutionary branches. It is argued that evolutionary gradualism and saltationism represent complementary modes of macroevolution corresponding respectively to the local and global topographical structure of adaptive surfaces.     

Myostatin三维结构模建及分子进化分析

Myostatin(MST)为肌肉生长负调节因子,其功能受抑制可导致肌肉量增加.对MST核酸序列进行序列比对,构建进化树;采用同源模建方法首次模建MST成熟肽生物活性二聚体的四级结构,并预测MST与其受体ActRIIB的相互作用模式.进化树将肌肉生长抑制素基因(MSTN)分成4个亚家族:哺乳动物MSTN,鸟类MSTN以及鱼类MSTN 1和2.MST受纯化选择作用,在不同物种的直系同源基因具有较高的刚源性,其中哺乳动物、鸟类MST C端活性肽氨基酸序列高度保守.表明哺乳动物、鸟类MST的结构、功能类似,且信号传导路径可能一致;而鱼类MST的调控机制可能存在较大差异.MST结构及其表面静电势和疏水氨基酸分布表明静电力和疏水相互作用在MST与其受体结合过程中可能起到十分重要的作用.     

Mutation-absorption model of the enzyme

Gradual changes in function of proteins in response to single changes in primary structure are often observed to occur and are a necessary condition for evolution by variation and natural selection at the protein level. A probabilistic (entropy theory_ analysis of the effect of changes in primary structure on three-dimensional shape and function shows that such gradualism is based on the presence of a control system in the molecule involving a definite general form of structure-function degeneracy. The assumptions of the analysis are that primary structure determines tertiary structure (or a thermal distribution of tertiary configurations and allosteric forms), tertiary structure determines function (characterized by rate and other parameters), and that certain features of tertiary structure may be specialized for particular functions. The main conclusion is that embodied in the molecule is a subsystem which serves as a buffer, absorbing mutation or other forms of genetic variation and expressing these as graceful variations in features of the shape critical for function. This buffer system may be realized by numerical redundancy of amino acids or other mechanisms which increase the redundancy of weak interactions responsible for folding, utilization of amino acids having a greater number of analogs with redundant features, or local and global structural formats which allow for more effective utilization of redundancy. The mutation-absorption model has implications for the interpretation of structure-function relations in biology, the topology of the adaptive landscape, the interpretation of isoenzymes and allozymes, the relationship between selection and neutralism in evolution, and the relation between the complexity of and energy required by biological systems and the effectiveness of evolutionary optimization.     

Genes that determine flower color: the role of regulatory changes in the evolution of phenotypic adaptations

A central goal of evolutionary genetics is to trace the causal pathway between mutations at particular genes and adaptation at the phenotypic level. The proximate objective is to identify adaptations through the analysis of molecular sequence data from specific candidate genes or their regulatory elements. In this paper, we consider the molecular evolution of floral color in the morning glory genus (Ipomoea) as a model for relating molecular and phenotypic evolution. To begin, flower color variation usually conforms to simple Mendelian transmission, thus facilitating genetic and molecular analyses. Population genetic studies of flower color polymorphisms in the common morning glory (Ipomoea purpurea) have shown that some morphs are subject to complex patterns of selection. Striking differences in floral color and morphology are also associated with speciation in the genus Ipomoea. The molecular bases for these adaptive shifts can be dissected because the biosynthetic pathways that determine floral pigmentation are well understood and many of the genes of flavonoid biosynthesis have been isolated and extensively studied. We present a comparative analysis of the level of gene expression in Ipomoea for several key genes in flavonoid biosynthesis. Specifically we ask: how frequently are adaptive shifts in flower color phenotypes associated with changes in regulation of gene expression versus mutations in structural genes? The results of this study show that most species differences in this crucial phenotype are associated with changes in the regulation of gene expression.     

Physicochemical evolution and molecular adaptation of the cetacean and artiodactyl cytochrome b proteins

Cetaceans have most likely experienced metabolic shifts since evolutionarily diverging from their terrestrial ancestors, shifts that may be reflected in the proteins such as cytochrome b that are responsible for metabolic efficiency. However, accepted statistical methods for detecting molecular adaptation are largely biased against even moderately conservative proteins because the primary criterion involves a comparison of nonsynonymous and synonymous substitution rates (dN/dS); they do not allow for the possibility that adaptation may come in the form of very few amino acid changes. We apply the MM01 model to the possible molecular adaptation of cytochrome b among cetaceans because it does not rely on a dN/dS ratio, instead evaluating positive selection in terms of the amino acid properties that comprise protein phenotypes that selection at the molecular level may act upon. We also apply the codon-degeneracy model (CDM), which focuses on evaluating overall patterns of nucleotide substitution in terms of base exchange, codon position, and synonymy to estimate the overall effect of selection. Using these relatively new models, we characterize the molecular adaptation that has occurred in the cetacean cytochrome b protein by comparing revealed amino acid replacement patterns to those found among artiodactyls, the modern terrestrial mammals found to be most closely related to cetaceans. Our findings suggest that several regions of the cetacean cytochrome b protein have experienced molecular adaptation. Also, these adaptations are spatially associated with domain structure, protein function, and the structure and function of the cytochrome bc(1) complex and its constituents. We also have found a general correlation between the results of the analytical software programs TreeSAAP (which implements the MM01 model) and CDM (which implements the codon-degeneracy model).     

Cancer Reviews 2021 series: Cancer evolution: Darwin and beyond

Clinical and laboratory studies over recent decades have established branched evolution as a feature of cancer. However, while grounded in somatic selection, several lines of evidence suggest a Darwinian model alone is insufficient to fully explain cancer evolution. First, the role of macroevolutionary events in tumour initiation and progression contradicts Darwin''s central thesis of gradualism. Whole‐genome doubling, chromosomal chromoplexy and chromothripsis represent examples of single catastrophic events which can drive tumour evolution. Second, neutral evolution can play a role in some tumours, indicating that selection is not always driving evolution. Third, increasing appreciation of the role of the ageing soma has led to recent generalised theories of age‐dependent carcinogenesis. Here, we review these concepts and others, which collectively argue for a model of cancer evolution which extends beyond Darwin. We also highlight clinical opportunities which can be grasped through targeting cancer vulnerabilities arising from non‐Darwinian patterns of evolution.     

Subtle evolutionary changes in the distribution of N-glycosylation sequons in the HIV-1 envelope glycoprotein 120

Many viruses are known to undergo rapid evolutionary changes under selective pressures. The HIV-1 envelope glycoprotein 120 (gp120) shows extreme selection for NXS/T sequons, the potential sites of N-glycosylation. Although the average number of sequons in gp120 appears to be relatively stable in the recent past, even slight changes in the distribution of sequons may potentially play crucial roles in protein interaction and viral infection. This study tracked the prevalence and distribution of NXS/T sequons in gp120 over a period of 29 years (from 1981 to 2009). The gp120 showed location specific distribution of sequons with higher density in the outer domain of the molecule. The NXT sequon density decreased in the outer domain (despite the increase in the sequon specific amino acid threonine), but increased in the inner domain. By contrast, the NXS sequon density increased specifically in the outer domain. Related changes were also seen in the distribution probabilities of sequons in two domains. The results indicate that the gp120, chiefly in subtype B, is redistributing NXS/T sequons within the molecule with specific selection for NXS sequons. The subtle evolution of sequons in gp120 may have implications in viral resistance and infection.     

Analysis of metabolic network based on conservation of molecular structure

Recent work has revealed much about chemical reactions inside hundreds of organisms as well as universal characteristics of metabolic networks, which shed light on the evolution of the networks. However, characteristics of individual metabolites have been neglected. For example, some carbohydrates have structures that are decomposed into small molecules by metabolic reactions, but coenzymes such as ATP are mostly preserved. Such differences in metabolite characteristics are important for understanding the universal characteristics of metabolic networks. To quantify the structure conservation of metabolites, we defined the "structure conservation index" (SCI) for each metabolite as the fraction of metabolite atoms restored to their original positions through metabolic reactions. As expected, coenzymes and coenzyme-like metabolites that have reaction loops in the network show a higher SCI. Using the index, we found that the sum of metabolic fluxes is negatively correlated with the structure preservation of metabolite. Also, we found that each reaction path around high SCI metabolites changes independently, while changes in reaction paths involving low SCI metabolites coincide through evolution processes. These correlations may provide a clue to universal properties of metabolic networks.     

Equilibrium state of molecular breeding

We investigate the equilibrium state of the model of Peng, et al. for molecular breeding. In the model, a population of DNA sequences is successively culled by removing the sequences with the lowest binding affinity to a particular target sequence. The remaining sequences are then amplified to restore the original population size, undergoing some degree of point-substitution of nucleotides in the process. Working in the infinite population size limit, we derive an equation for the equilibrium distribution of binding affinity, here modeled by the number of matches to the target sequence. The equation is then solved approximately in the limit of large sequence length, in the three regimes of strong, intermediate and weak selection. The approximate solutions are verified via comparison to exact numerical results.     

Population structure, levels of selection, and the evolution of intracellular symbionts

Many obligately intracellular symbionts exhibit a characteristic set of genetic changes that include an increase in substitution rates, loss of many genes, and apparent destabilization of many proteins and structural RNAs. Authors have suggested that these changes are due to increased mutation rates, or, more commonly, decreased effective population size due to population bottlenecks at the symbiont or, perhaps, host level. I propose that the increase in substitution rates and accumulation of deleterious mutations is a consequence of the population structure imposed on the endosymbionts by strict host association, loss of horizontal transmission and potentially conflicting levels of selection. I analyze a population genetic model of endosymbiont evolution, and demonstrate that substitution rates will increase, and the effect of those substitutions on endosymbiont fitness will become more deleterious as horizontal transmission among hosts decreases. Additionally, I find that there is a critical level of horizontal transmission below which natural selection cannot effectively purge deleterious mutations, leading to an expected loss of fitness over time. This critical level varies across loci with the degree of correlation between host and endosymbiont fitness, and may help explain differential retention and loss of certain genes.     

Extreme differences in charge changes during protein evolution

Summary The maintenance of a proper distribution of charged amino acid residues might be expected to be an important factor in protein evolution. We therefore compared the inferred changes in charge during the evolution of 43 protein families with the changes expected on the basis of random base substitutions. It was found that certain proteins, like the eye lens crystallins and most histones, display an extreme avoidance of changes in charge. Other proteins, like phospholipase A2 and ferredoxin, apparently have sustained more charged replacements than expected, suggesting a positive selection for changes in charge. Depending on function and structure of a protein, charged residues apparently can be important targets for selective forces in protein evolution. It appears that actual biased codon usage tends to decrease the proportion of charged amino acid replacements. The influence of nonrandomness of mutations is more equivocal. Genes that use the mitochondrial instead of the universal code lower the probability that charge changes will occur in the encoded proteins.     

Towards an alternative evolution model

Lamarck and Darwin agreed on the inconstancy of species and on the exclusive gradualism of evolution (nature does not jump). Darwinism, revived as neo-Darwinism, was almost generally accepted from about 1930 till 1960. In the sixties the evolutionary importance of selection has been called in question by the neutralists. The traditional conception of the gene is disarranged by recent molecular-biological findings. Owing to the increasing confusion about the concept of genotype, this concept is reconsidered. The idea of the genotype as a cluster of genes is replaced by a cybernetical interpretation of the genotype. As nature does jump, exclusive gradualism is dismissed. Saltatory evolution is a natural phenomenon, provoked by a sudden collapse of the thresholds which resist against evolution. The fossil record and the taxonomic system call for a macromutational interpretation. As Lamarck and Darwin overlooked the resistance of evolutionary thresholds, an alternative evolution model is needed, the first to be constructed on a palaeontological and taxonomic basis.Progress emerges from error, far more than from chaos. (Gray, 1963)     

Estimation of evolutionary distances from protein spatial structures

New equations are derived to estimate the number of amino acid substitutions per site between two homologous proteins from the root mean square (RMS) deviation between two spatial structures and from the fraction of identical residues between two sequences. The equations are based on evolutionary models, analyzing predominantly structural changes and not sequence changes. Evolution of spatial structure is treated as a diffusion in an elastic force field. Diffusion accounts for structural changes caused by amino acid substitutions, and elastic force reflects selection, which preserves protein fold. Obtained equations are supported by analysis of protein spatial structures. Received: 21 September 1995 / Accepted: 19 May 1997     

载脂蛋白基因家族的密码子空间分析——核苷酸替换的非随机进化选择

ＤＮＡ分子进化中,对核苷酸替换的选择可呈选择中性或选择倾向性。为研究载脂蛋白基因进化过程中对核苷酸变化的选择方式,本文建立了基因的密码子空间分析方法。密码子空间是由密码子３个位置上核苷酸出现机率所组成的矩阵。对该空间中核苷酸分布的非随机性度量可以反映进化过程中核苷酸替换的选择方式。应用该法,我们发现载脂蛋白基因密码子空间第一及第三位的核苷酸分布呈高度非随机性。进一步研究表明：这种核苷酸的非随机分布可能与腺苷酸、胸苷酸对密码子位置的非中性选择有关。此外,还研究了同义密码子的选择使用与分支种系发生的关系。结果显示：载脂蛋白分子演化中存在着同义密码子使用的分子进化钟。这些研究提示密码子空间中核苷酸替换的非随机选择可能是载脂蛋白基因进化的一种特征。     

Quantitative genetic variation in Daphnia: temporal changes in genetic architecture

Nonadditive genetic variation and genetic disequilibrium are two important factors that influence the evolutionary trajectory of natural populations. We assayed quantitative genetic variation in a temporary-pond-dwelling population of Daphnia pulex over a full season to examine the role of nonadditive genetic variation and genetic disequilibrium in determining the short-term evolutionary trajectory of a cyclic parthenogen. Quantitative traits were influenced by three factors: (1) clonal selection significantly changed the population mean phenotype during the course of the growing season; (2) sexual reproduction and recombination led to significant changes in life-history trait means and the levels of expressed genetic variation, implying the presence of substantial nonadditive genetic variation and genetic disequilibrium; and (3) Egg-bank effects were found to be an important component of the realized year-to-year change. Additionally, we examined the impact of genetic disequilibria induced by clonal selection on the genetic (co)variance structure with a common principal components model. Clonal selection caused significant changes in the (co)variance structure that were eliminated by a single bout of random mating, suggesting that a build-up of disequilibria was the primary source of changes in the (co)variance structure. The results of this study highlight the complexity of natural selection operating on populations that undergo alternating phases of sexual and asexual reproduction.     

Structures of in vitro evolved binding sites on neocarzinostatin scaffold reveal unanticipated evolutionary pathways

We have recently applied in vitro evolution methods to create in Neocarzinostatin a new binding site for a target molecule unrelated to its natural ligand. The main objective of this work was to solve the structure of some of the selected binders in complex with the target molecule: testosterone. Three proteins (1a.15, 3.24 and 4.1) were chosen as representative members of sequence families that came out of the selection process within different randomization schemes. In order to evaluate ligand-induced conformational adaptation, we also determined the structure of one of the proteins (3.24) in the free and complexed forms. Surprisingly, all these mutants bind not one but two molecules of testosterone in two very different ways. The 3.24 structure revealed that the protein spontaneously evolved in the system to bind two ligand molecules in one single binding crevice. These two binding sites are formed by substituted as well as by non-variable side-chains. The comparison with the free structure shows that only limited structural changes are observed upon ligand binding. The X-ray structures of the complex formed by 1a.15 and 4.1 Neocarzinostatin mutants revealed that the two variants form very similar dimers. These dimers were observed neither for the uncomplexed variants nor for wild-type Neocarzinostatin but were shown here to be induced by ligand binding. Comparison of the three complexed forms clearly suggests that these unanticipated structural responses resulted from the molecular arrangement used for the selection experiments.     

Compensatory Evolution of a WW Domain Variant Lacking the Strictly Conserved Trp Residue

Replacement of conserved amino acid residues during evolution of proteins can lead to divergence and the formation of new families with novel functions, but is often deleterious to both protein structure and function. Using the WW domain, we experimentally examined whether and to what degree second-site mutations can compensate for the reduction of function and loss of structure that accompany substitution of a strictly conserved amino acid residue. The W17F mutant of the WW domain, with substitution of the most strictly conserved Trp residue, is known to lack a specific three-dimensional structure and shows reduced binding affinity in comparison to the wild type. To obtain second-site revertants, we performed a selection experiment based on the proline-rich peptide (PY ligand) binding affinity using the W17F mutant as the initial sequence. After selection by ribosome display, we were able to select revertants that exhibited a maximum ninefold higher affinity to the PY ligand than the W17F mutant and showed an even better affinity than the wild type. In addition, we found that the functional restoration resulted in increased binding specificity in selected revertants, and the structures were more compact, with increased amounts of secondary structure, in comparison to the W17F mutant. Our results suggest that the defective structure and function of the proteins caused by mutations in highly conserved residues occurring through divergent evolution not only can be restored but can be further improved by compensatory mutations.     

TRANSLATING BETWEEN MICROEVOLUTIONARY PROCESS AND MACROEVOLUTIONARY PATTERNS: THE CORRELATION STRUCTURE OF INTERSPECIFIC DATA

As species evolve along a phylogenetic tree, we expect closely related species to retain some phenotypic similarities due to their shared evolutionary histories. The amount of expected similarity depends both on the hierarchical phylogenetic structure, and on the specific magnitude and types of evolutionary changes that accumulate during each generation. In this study, we show how models of microevolutionary change can be translated into the resulting macroevolutionary patterns. We illustrate how the structure of phenotypic covariances expected in interspecific measurements can be derived, and how this structure depends on the microevolutionary forces guiding phenotypic change at each generation. We then explore the covariance structure expected from several simple microevolutionary models of phenotypic evolution, including various combinations of random genetic drift, directional selection, stabilizing selection, and environmental change, as well as models of punctuated or burst-like evolution. We find that stabilizing selection leads to patterns of exponential decrease of between species covariance with phylogenetic distance. This is different from the usual linear patterns of decrease assumed in most comparative and systematic methods. Nevertheless, linear patterns of decrease can result from many processes in addition to random genetic drift, such as directional and fluctuating selection as well as modes of punctuated change. Our framework can be used to develop methods for (1) phylogenetic reconstruction; (2) inference of the evolutionary process from comparative data; and (3) conducting or evaluating statistical analyses of comparative data while taking phylogenetic history into account.