Selective pressures for accurate altruism targeting: evidence from digital evolution for difficult-to-test aspects of inclusive fitness theory

Inclusive fitness theory predicts that natural selection will favour altruist genes that are more accurate in targeting altruism only to copies of themselves. In this paper, we provide evidence from digital evolution in support of this prediction by competing multiple altruist-targeting mechanisms that vary in their accuracy in determining whether a potential target for altruism carries a copy of the altruist gene. We compete altruism-targeting mechanisms based on (i) kinship (kin targeting), (ii) genetic similarity at a level greater than that expected of kin (similarity targeting), and (iii) perfect knowledge of the presence of an altruist gene (green beard targeting). Natural selection always favoured the most accurate targeting mechanism available. Our investigations also revealed that evolution did not increase the altruism level when all green beard altruists used the same phenotypic marker. The green beard altruism levels stably increased only when mutations that changed the altruism level also changed the marker (e.g. beard colour), such that beard colour reliably indicated the altruism level. For kin- and similarity-targeting mechanisms, we found that evolution was able to stably adjust altruism levels. Our results confirm that natural selection favours altruist genes that are increasingly accurate in targeting altruism to only their copies. Our work also emphasizes that the concept of targeting accuracy must include both the presence of an altruist gene and the level of altruism it produces.     

Compositional patterns in vertebrate genomes: Conservation and change in evolution

Summary The evolution of vertebrate genomes can be investigated by analyzing their regional compositional patterns, namely the compositional distributions of large DNA fragments (in the 30–100-kb size range), of coding sequences, and of their different codon positions. This approach has shown the existence of two evolutionary modes. In the conservative mode, compositional patterns are maintained over long times (many million years), in spite of the accumulation of enormous numbers of base substitutions. In the transitional, or shifting, mode, compositional patterns change into new ones over much shorter times.The conservation of compositional patterns, which has been investigated in mammalian genomes, appears to be due in part to some measure of compositional conservation in the base substitution process, and in part to negative selection acting at regional (isochore) levels in the genome and eliminating deviations from a narrow range of values, presumably corresponding to optimal functional properties. On the other hand, shifts of compositional patterns, such as those that occurred between cold-blooded and warm-blooded vertebrates, appear to be due essentially to both negative and positive selection again operating at the isochore level, largely under the influence of changes in environmental conditions, and possibly taking advantage of mutational biases in the replication/repair enzymes and/or in the enzyme make-up of nucleotide precursor pools. Other events (like translocations and changes in chromosomal structure) also play a role in the transitional mode of genome evolution.The present findings (1) indicate that isochores, which correspond to the DNA segments of individual or contiguous chromatin domains, represent selection units in the vertebrate genome; and (2) shed new light on the selectionist-neutralist controversy.This work was presented at the EMBO Workshop on Evolution (Cambridge, UK, 4–6 July 1988) and at the 16th International Congress of Genetics (Toronto, Canada, 20–27 August 1988)     

分子生物学与进化的新理论

论述了由于分子生物学取得的成就,而形成的生物进化的新理论,并阐述了若干新理论指导生产实践取得的巨大应用成果。     

Thermodynamical interpretation of an adaptive walk on a Mt. Fuji-type fitness landscape: Einstein relation-like formula holds in a stochastic evolution

We have theoretically studied the statistical properties of adaptive walks (or hill-climbing) on a Mt. Fuji-type fitness landscape in the multi-dimensional sequence space through mathematical analysis and computer simulation. The adaptive walk is characterized by the "mutation distance" d as the step-width of the walker and the "population size" N as the number of randomly generated d-fold point mutants to be screened. In addition to the fitness W, we introduced the following quantities analogous to thermodynamical concepts: "free fitness" G(W) is identical with W+T x S(W), where T is the "evolutionary temperature" T infinity square root of d/lnN and S(W) is the entropy as a function of W, and the "evolutionary force" X is identical with d(G(W)/T)/dW, that is caused by the mutation and selection pressure. It is known that a single adaptive walker rapidly climbs on the fitness landscape up to the stationary state where a "mutation-selection-random drift balance" is kept. In our interpretation, the walker tends to the maximal free fitness state, driven by the evolutionary force X. Our major findings are as follows: First, near the stationary point W*, the "climbing rate" J as the expected fitness change per generation is described by J approximately L x X with L approximately V/2, where V is the variance of fitness distribution on a local landscape. This simple relationship is analogous to the well-known Einstein relation in Brownian motion. Second, the "biological information gain" (DeltaG/T) through adaptive walk can be described by combining the Shannon's information gain (DeltaS) and the "fitness information gain" (DeltaW/T).     

Toward a molecular taxonomy

Summary The concepts of a super information source and ensemble averaging are used to estimate the amount of information stored in protein andt-RNA sequences. Specifically applied to cytochrome c and hemoglobins, information measures analogous to those found to be highly significant for DNA pair frequency data (D₂ vs.R) by Gatlin (1968) prove to be extremely highly correlated with organism complexity. Super source stability and the possible taxonomic utility of the extraordinary clusterings obtained are discussed. A restriction on the construction of ancestral sequences and a possible handle on homology are also detailed.     

On the molecular evolutionary clock

Summary The conceptual framework surrounding the origin of the molecular evolutionary clock and circumstances of this origin are described. In regard to the quest for the best available molecular clocks, a return to protein clocks is conditionally recommended. On the basis of recent data and certain considerations, it is pointed out that the realm of neutrality in evolution is probably less extensive than is now commonly thought, in the three distinct senses of the term neutrality—neutrality as nonfunctionality of mutations, neutrality as equifunctionality of mutations, and neutrality as a mode of fixation of mutations. The possibility is raised that complex sets of interacting components forming a system that is bounded with respect to its environment may quite generally display an intrinsic trend to a quasi-clockwise evolutionary behavior.     

Dating the early evolution of plants: detection and molecular clock analyses of orthologs

Orthologs generally are under selective pressure against loss of function, while paralogs usually accumulate mutations and finally die or deviate in terms of function or regulation. Most ortholog detection methods contaminate the resulting datasets with a substantial amount of paralogs. Therefore we aimed to implement a straightforward method that allows the detection of ortholog clusters with a reduced amount of paralogs from completely sequenced genomes. The described cross-species expansion of the reciprocal best BLAST hit method is a time-effective method for ortholog detection, which results in 68% truly orthologous clusters and the procedure specifically enriches single-copy orthologs. The detection of true orthologs can provide a phylogenetic toolkit to better understand evolutionary processes. In a study across six photosynthetic eukaryotes, nuclear genes of putative mitochondrial origin were shown to be over-represented among single copy orthologs. These orthologs are involved in fundamental biological processes like amino acid metabolism or translation. Molecular clock analyses based on this dataset yielded divergence time estimates for the red/green algae (1,142 MYA), green algae/land plant (725 MYA), mosses/seed plant (496 MYA), gymno-/angiosperm (385 MYA) and monocotyledons/core eudicotyledons (301 MYA) divergence times. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The theory of everything and molecular evolution

The theory of everything is discussed in relationship to early bacterial molecular evolution. The emphasis is on time, space (or location at the molecular level), the universal construction kit (elements contained in periodic table) and change per units of time that were necessary for early bacterial molecular evolution to occur.     

Phycobilisomes linker family in cyanobacterial genomes: divergence and evolution

Cyanobacteria are the oldest life form making important contributions to global CO2 fixation on the Earth. Phycobilisomes (PBSs) are the major light harvesting systems of most cyanobacteria species. Recent availability of the whole genome database of cyanobacteria provides us a global and further view on the complex structural PBSs. A PBSs linker family is crucial in structure and function of major light-harvesting PBSs complexes. Linker polypeptides are considered to have the same ancestor with other phycobiliproteins (PBPs), and might have been diverged and evolved under particularly selective forces together. In this paper, a total of 192 putative linkers including 167 putative PBSs-associated linker genes and 25 Ferredoxin-NADP oxidoreductase (FNR) genes were detected through whole genome analysis of all 25 cyanobacterial genomes (20 finished and 5 in draft state). We compared the PBSs linker family of cyanobacteria in terms of gene structure, chromosome location, conservation domain, and polymorphic variants, and discussed the features and functions of the PBSs linker family. Most of PBSs-associated linkers in PBSs linker family are assembled into gene clusters with PBPs. A phylogenetic analysis based on protein data demonstrates a possibility of six classes of the linker family in cyanobacteria. Emergence, divergence, and disappearance of PBSs linkers among cyanobacterial species were due to speciation, gene duplication, gene transfer, or gene loss, and acclimation to various environmental selective pressures especially light.     

Enzyme-coding genes as molecular clocks: The molecular evolution of animal alpha-amylases

Summary We constructed a cDNA library for the beetle,Tribolium castaneum. This library was screened using a cloned amylase gene fromDrosophila melanogaster as a molecular probe. Beetle amylase cDNA clones were isolated from this bank, and the nucleotide sequence was obtained for a cDNA clone with a coding capacity for 228 amino acids. Both the nucleotide sequence and predicted amino acid sequence were compared to our recent results forD. melanogaster alpha-amylases, along with published sequences for other alpha-amylases. The results show that animal alpha-amylases are highly conserved over their entire length. A borader comparison, which includes plant and microbial alpha-amylase sequences, indicates that parts of the gene are conserved between prokaryotes, plants, and animals. We discuss the potential importance of this and other enzyme-coding genes for the construction of molecular phylogenies and for the study of the general question of molecular clocks in evolution.     

Nonuniformity of nucleotide substitution rates in molecular evolution: Computer simulation and analysis of 5S ribosomal RNA sequences

Summary The effects of temporal (among different branches of a phylogeny) and spatial (among different nucleotide sites within a gene) nonuniformities of nucleotide substitution rates on the construction of phylogenetic trees from nucleotide sequences are addressed. Spatial nonuniformity may be estimated by using Shannon's (1948) entropy formula to measure the Relative Nucleotide Variability (RNV) at each nucleotide site in an aligned set of sequences; this is demonstrated by a comparative analysis of 5S rRNAs. New methods of constructing phylogenetic trees are proposed that augment the Unweighted Pair-Group Using Arithmetic Averages (UPGMA) algorithm by estimating and compensating for both spatial and temporal nonuniformity in substitution rates. These methods are evaluated by computer simulations of 5S rRNA evolution that include both kinds of nonuniformities. It was found that the proposed Reference Ratio Method improved both the ability to reconstruct the correct topology of a tree and also the estimation of branch lengths as compared to UPGMA. A previous method (Farris et al. 1970; Klotz et al. 1979; Li 1981) was found to be less successful in reconstructing topologies when there is high probability of multiple mutations at some sites. Phylogenetic analyses of 5S rRNA sequences support the endosymbiotic origins of both chloroplasts and mitochondria, even though the latter exhibit an accelerated rate of nucleotide substitution. Phylogenetic trees also reveal an adaptive radiation within the eubacteria and another within the eukaryotes for the origins of most major phyla within each group during the Precambrian era.     

Functional and morphological stasis during molecular evolution

The evolutionary distance between two sets of proteins was estimated using the techniques of Miyata and Yasunaga (1980) and Kimura (1980). Human beta 2-microglobulin was compared with the homologous murine molecule, while human and equine alpha-globin were similarly treated. It was found that a large amount of molecular evolution has occurred in beta 2-microglobulin since its divergence from the common ancestor of mice and humans. Kimura's estimate of evolutionary distance, K, is 0.353, while those of Miyata and Yasunaga are KS = 0.708 and KA = 0.171. The respective values for human and equine alpha-globin are 0.152, 0.293, and 0.084. In spite of this molecular evolution, it is shown that murine beta 2-microglobulin can effect the expression of HLA class I antigens on the surface of human-mouse hybrid cells and that the tertiary structures of human and equine deoxyhemoglobin are nearly identical. These observations are discussed in the light of Kimura's theory of neutral allelic drift.     

The evolution of protamine P1 genes in dasyurid marsupials

We report the complete DNA sequences of the protamine P1 gene and flanking regions for 13 species of the marsupial family Dasyuridae. The structure of the protamine locus is conserved in dasyurids and consists of two exons (of lengths 142–151 and 47 bp) separated by an intron (208–240 bp). A key feature of the dasyund intron is a 38–40 by duplication found in all species examined to date. This duplication apparently predates the radiation of modern dasyurid lineages and may be homologous to a similar feature in the marsupial mole (Notoryctes). Sequences from a species of Planigale demonstrate that this genus is unique among marsupials in possessing cysteine residues in its protamine P1 molecules. Cysteines may provide enhanced chemical stability for condensed sperm nuclei, a physiological feature that would converge on the common eutherian pattern. Phylogenetic analysis of the protamine genes yields a tree that is largely congruent with previous molecular systematic studies in two areas: (1) There are three main dasyurid lineages corresponding to the Sminthopsinae, Dasyurinae, and Phascogalinae; (2) Dasyurinae and Phascogalinae are sister groups. This study is the first estimate of dasyurid relationships based on a nuclear DNA sequence.Correspondence to: J.D. Retief     

The molecular evolution of pancreatic ribonuclease

Summary The primary structures of pancreatic ribonucleases from 26 species (18 artiodactyls, horse, whale, 5 rodents and turtle) are known. Several species contain identical ribonucleases (cow/bison; sheep/goat), other species show polymorphism (arabian camel) or the presence of two structural gene loci (guinea pig pancreas contains two ribonucleases that differ at 31 positions). 26 different sequences (including the ribonuclease from bovine seminal plasma which is paralogous to the pancreatic ribonucleases) were used to construct a most parsimonious tree. A second tree that most closely approximates current biological opinion requires 402 whereas the most parsimonious tree requires 389 nucleotide substitutions. The artiodactyl part of the most parsimonious tree conforms quite well with the biological one of this order, except for the position of the giraffe which is placed with the pronghorn. Other parts of the most parsimonious tree agree less with the biological tree, probably as a result of the occurrence of many parallel and back substitutions. Bovine seminal ribonuclease was found to be the result of a gene duplication which occurred before the divergence of the true ruminants, but after the divergence of this group from the cameloids.The evolutionary rate of ribonuclease was found to be 390, 3.0 and 11 nucleotide substitutions per 10⁹ yrs per ribonuclease gene, codon and covarion respectively. However, there is much variation in evolutionary rate in different taxa. Values ranging from about 100 (in the bovidae) to about 700 (in the rodents) nucleotide substitutions per 10⁹ yrs per gene were found.A method for counting parallel and back mutations is presented. The 389 nucleotide substitutions in the most parsimonious tree occur at 88 codon positions; 154 of them are the result of parallel and back mutations. Parallel evolution to a similar structure, including the presence of 2 sites with carbohydrate, was demonstrated in an extensive region at the surface of pig and guinea pig ribonuclease B. The presence of carbohydrate probably is important in a number of species. A correlation between the presence of heavily glycosidated ribonucleases and coecal digestion was observed. Hypothetical sequences of ancestral ungulate ribonucleases contain many recognition sites for carbohydrate attachment; this suggests that herbivores with coecal digestion might have preceded the true ruminants in mammalian evolution.     

A model for the evolution of paralog families in genomes

We introduce and analyse a simple probabilistic model of genome evolution. It is based on three fundamental evolutionary events: gene loss, duplication and accumulated change. This is motivated by previous works which consisted in fitting the available genomic data into, what is called paralog distributions. This formalism is described by a system of infinite number of linear equations. We show that this system generates a semigroup of linear operators on the space l ¹. We prove that size distribution of paralogous gene families in a genome converges to the equilibrium as time goes to infinity. Moreover we show that when probabilities of gene removal and duplication are close to each other, then the resulting distribution is close to logarithmic distribution. Some empirical results for yeast genomes are presented.     

Improved production and function of llama heavy chain antibody fragments by molecular evolution

The aim of this study was to improve production level of llama heavy chain antibody fragments (V_HH) in Saccharomyces cerevisiae while retaining functional characteristics. For this purpose, the DNA shuffling technique was used on llama V_HH fragments specific for the azo-dye reactive red-6. In the DNA shuffling process, three parental llama V_HH with high amino acid sequence identity with significant differences in production and functional characteristics were used. From these parental sequences, a S. cerevisiae library was created and 16 antigen specific shuffled V_HH fragments were selected. We found that these shuffled V_HH fragments were, (i) unique in sequence; (ii) composed of two or three parental sequences; (iii) in three V_HHs point mutations occurred; and (iv) antigen specificity was not changed. The four highest producers in the yeast S. cerevisiae were selected and production, affinity, and antigen binding at 90°C were compared with parental V_HHs. One shuffled V_HH was enhanced both in production (3.4-fold) and affinity (four-fold). A second shuffled V_HH displayed increased production (1.9-fold), and improved stability (2.4-fold) in antigen binding at 90°C. Structural analysis suggested that improved antigen binding is associated with the A24→V24 substitution, which reduces the size of the hydrophobic pit at the llama V_HH surface. We demonstrate that it is possible to improve desired characteristics of the same V_HH fragment simultaneously using DNA shuffling. Finally, this is one of the first examples of DNA shuffling improving temperature stability of an antibody fragment.