Selectionism and neutralism in molecular evolution

Genes under positive selection in Mycobacterium tuberculosis

We employed an evolutionary genomics approach to detect genes under lineage-specific positive selection for the two closely related Mycobacterium tuberculosis strains, the virulent H37Rv and the avirulent H37Ra, with the clinical isolate CDC1551 as the outgroup. We found six H37Rv-specific and six H37Ra-specific positively selected genes, among which the former comprised a flavoprotein, a RNA polymerase sigma factor SigM, two PPE family proteins, as well as two hypothetical proteins, while the latter consisted of a dehydrogenase, a (3R)-hydroxyacyl-ACP dehydratase subunit HadA, a PPE family protein, and three PE-PGRS family proteins. Obviously, the PE/PPE/PE-PGRS family proteins were the main targets of positive selection. The functional discussion of our findings implied that those positively selected genes were highly involved in antigen variations and immune evasions of Mycobacterium tuberculosis.

The role of the effective population size in compensatory evolution

The impact of the effective population size (N_e) on the efficacy of selection has been the focus of many theoretical and empirical studies over the recent years. Yet, the effect of N_e on evolution under epistatic fitness interactions is not well understood. In this study, we compare selective constraints at independently evolving (unpaired) and coevolving (paired) sites in orthologous transfer RNAs (tRNA molecules for vertebrate and drosophilid species pairs of different N_e. We show that patterns of nucleotide variation for the two classes of sites are explained well by Kimura's one- and two-locus models of sequence evolution under mutational pressure. We find that constraints in orthologous tRNAs increase with increasing N_e of the investigated species pair. Thereby, the effect of N_e on the efficacy of selection is stronger at unpaired sites than at paired sites. Furthermore, we identify a “core” set of tRNAs with high structural similarity to tRNAs from all major kingdoms of life and a “peripheral” set with lower similarity. We observe that tRNAs in the former set are subject to higher constraints and less prone to the effect of N_e, whereas constraints in tRNAs of the latter set show a large influence of N_e. Finally, we are able to demonstrate that constraints are relaxed in X-linked drosophilid tRNAs compared with autosomal tRNAs and suggest that N_e is responsible for this difference. The observed effects of N_e are consistent with the hypothesis that evolution of most tRNAs is governed by slightly to moderately deleterious mutations (i.e., |N_es| ≤ 5).

Cloning and characterization of boron transporters in Brassica napus

Six full-length cDNA encoding boron transporters (BOR) were isolated from Brassica napus (AACC) by rapid amplification of cDNA ends (RACE). The phylogenic analysis revealed that the six BORs were the orthologues of AtBOR1, which formed companying with the triplication and allotetra-ploidization process of B. napus, and were divided into three groups in B. napus. Each group was comprised of two members, one of which was originated from Brassica rapa (AA) and the other from Brassica oleracea (CC). Based on the phylogenetic relationships, the six genes were named as BnBOR1;1a, BnBOR1;1c, BnBOR1;2a, BnBOR1;2c, BnBOR1;3a and BnBOR1;3c, respectively. The deduced BnBOR1 s had extensive similarity with other plant BORs, with the identity of 74-96.8% in amino acid sequence. The BnBOR1;3a and BnBOR1;3c resembled AtBOR1 in number and positions of the 11 introns, but the others only have 9 introns. After the gene duplication, there was evidence of purifying selection under a divergent selective pressure. The expression patterns of the six BnBOR1 s were detected by semi-quantitative RT-PCR. The BnBOR1;3a and BnBOR1;3c showed a ubiquitous expression in all of the investigated tissues, whereas the other four genes showed similar tissue-specific expression profile. Unlike the non-transcriptional regulation of AtBOR1, the expression of BnBOR1;1c and BnBOR1;2a were obviously induced by boron deficiency. This study suggested that the BOR1 s had undergone a divergent expression pattern in the genome of B. napus after that the B. napus diverged from Arabidopsis thaliana.

Revisiting the yeast PPR proteins--application of an Iterative Hidden Markov Model algorithm reveals new members of the rapidly evolving family

Pentatricopeptide repeat (PPR) proteins are the largest known RNA-binding protein family, and are found in all eukaryotes, being particularly abundant in higher plants. PPR proteins localize mostly to mitochondria and chloroplasts, and many were shown to modulate organellar genome expression on the posttranscriptional level. Although the genomes of land plants encode hundreds of PPR proteins, only a few have been identified in Fungi and Metazoa. As the current PPR motif profiles are built mainly on the basis of the predominant plant sequences, they are unlikely to be optimal for detecting fungal and animal members of the family, and many putative PPR proteins in these genomes may remain undetected. In order to verify this hypothesis, we designed a hidden Markov model-based bioinformatic tool called Supervised Clustering-based Iterative Phylogenetic Hidden Markov Model algorithm for the Evaluation of tandem Repeat motif families (SCIPHER) using sequence data from orthologous clusters from available yeast genomes. This approach allowed us to assign 12 new proteins in Saccharomyces cerevisiae to the PPR family. Similarly, in other yeast species, we obtained a 5-fold increase in the detection of PPR motifs, compared with the previous tools. All the newly identified S. cerevisiae PPR proteins localize in the mitochondrion and are a part of the RNA processing interaction network. Furthermore, the yeast PPR proteins seem to undergo an accelerated divergent evolution. Analysis of single and double amino acid substitutions in the Dmr1 protein of S. cerevisiae suggests that cooperative interactions between motifs and pseudoreversion could be the force driving this rapid evolution.

Expression level, cellular compartment and metabolic network position all influence the average selective constraint on mammalian enzymes

BACKGROUND

A gene's position in regulatory, protein interaction or metabolic networks can be predictive of the strength of purifying selection acting on it, but these relationships are neither universal nor invariably strong. Following work in bacteria, fungi and invertebrate animals, we explore the relationship between selective constraint and metabolic function in mammals.

RESULTS

We measure the association between selective constraint, estimated by the ratio of nonsynonymous (Ka) to synonymous (Ks) substitutions, and several, primarily metabolic, measures of gene function. We find significant differences between the selective constraints acting on enzyme-coding genes from different cellular compartments, with the nucleus showing higher constraint than genes from either the cytoplasm or the mitochondria. Among metabolic genes, the centrality of an enzyme in the metabolic network is significantly correlated with Ka/Ks. In contrast to yeasts, gene expression magnitude does not appear to be the primary predictor of selective constraint in these organisms.

CONCLUSIONS

Our results imply that the relationship between selective constraint and enzyme centrality is complex: the strength of selective constraint acting on mammalian genes is quite variable and does not appear to exclusively follow patterns seen in other organisms.

Pharmacogenetics in medico-legal context

Medico-legal autopsy is the primary method in determining the cause and manner of death when the death is suspected to be unnatural. In some of these autopsies, the death remains ambiguous, even after a complete autopsy including histological investigation and toxicological screenings. In cases where there are no morphological abnormalities, medico-legal genetics may offer additional means to provide knowledge of possible genetic mutations, which may have initiated the process or predisposed the individual to stress risk conditions leading to death. One class of ambiguous deaths consists of drug-related deaths where the interpretation of the toxicological results are not clear. In such situations post mortem genotyping and the analysis of metabolite rations may provide an insight to the findings. A few cases demonstrating the potential strength of pharmacogenetics in medico-legal context has been published. However, there is a paramount need for serious scientific studies before the field of post mortem pharmacogenetics can be utilized in routine medico-legal analyses casework and brought routinely into courtroom.

Community transcriptomics reveals universal patterns of protein sequence conservation in natural microbial communities

Background

Combined metagenomic and metatranscriptomic datasets make it possible to study the molecular evolution of diverse microbial species recovered from their native habitats. The link between gene expression level and sequence conservation was examined using shotgun pyrosequencing of microbial community DNA and RNA from diverse marine environments, and from forest soil.

Results

Across all samples, expressed genes with transcripts in the RNA sample were significantly more conserved than non-expressed gene sets relative to best matches in reference databases. This discrepancy, observed for many diverse individual genomes and across entire communities, coincided with a shift in amino acid usage between these gene fractions. Expressed genes trended toward GC-enriched amino acids, consistent with a hypothesis of higher levels of functional constraint in this gene pool. Highly expressed genes were significantly more likely to fall within an orthologous gene set shared between closely related taxa (core genes). However, non-core genes, when expressed above the level of detection, were, on average, significantly more highly expressed than core genes based on transcript abundance normalized to gene abundance. Finally, expressed genes showed broad similarities in function across samples, being relatively enriched in genes of energy metabolism and underrepresented by genes of cell growth.

Conclusions

These patterns support the hypothesis, predicated on studies of model organisms, that gene expression level is a primary correlate of evolutionary rate across diverse microbial taxa from natural environments. Despite their complexity, meta-omic datasets can reveal broad evolutionary patterns across taxonomically, functionally, and environmentally diverse communities.

Noise driven evolutionary waves

Adaptation in spatially extended populations entails the propagation of evolutionary novelties across habitat ranges. Driven by natural selection, beneficial mutations sweep through the population in a "wave of advance". The standard model for these traveling waves, due to R. Fisher and A. Kolmogorov, plays an important role in many scientific areas besides evolution, such as ecology, epidemiology, chemical kinetics, and recently even in particle physics. Here, we extend the Fisher-Kolmogorov model to account for mutations that confer an increase in the density of the population, for instance as a result of an improved metabolic efficiency. We show that these mutations invade by the action of random genetic drift, even if the mutations are slightly deleterious. The ensuing class of noise-driven waves are characterized by a wave speed that decreases with increasing population sizes, contrary to conventional Fisher-Kolmogorov waves. When a trade-off exists between density and growth rate, an evolutionary optimal population density can be predicted. Our simulations and analytical results show that genetic drift in conjunction with spatial structure promotes the economical use of limited resources. The simplicity of our model, which lacks any complex interactions between individuals, suggests that noise-induced pattern formation may arise in many complex biological systems including evolution.

Correlated evolution of nearby residues in Drosophilid proteins

Here we investigate the correlations between coding sequence substitutions as a function of their separation along the protein sequence. We consider both substitutions between the reference genomes of several Drosophilids as well as polymorphisms in a population sample of Zimbabwean Drosophila melanogaster. We find that amino acid substitutions are “clustered” along the protein sequence, that is, the frequency of additional substitutions is strongly enhanced within ≈10 residues of a first such substitution. No such clustering is observed for synonymous substitutions, supporting a “correlation length” associated with selection on proteins as the causative mechanism. Clustering is stronger between substitutions that arose in the same lineage than it is between substitutions that arose in different lineages. We consider several possible origins of clustering, concluding that epistasis (interactions between amino acids within a protein that affect function) and positional heterogeneity in the strength of purifying selection are primarily responsible. The role of epistasis is directly supported by the tendency of nearby substitutions that arose on the same lineage to preserve the total charge of the residues within the correlation length and by the preferential cosegregation of neighboring derived alleles in our population sample. We interpret the observed length scale of clustering as a statistical reflection of the functional locality (or modularity) of proteins: amino acids that are near each other on the protein backbone are more likely to contribute to, and collaborate toward, a common subfunction.

Genes are templates for proteins, yet evolutionary studies of genes and proteins often bear little resemblance. Analyses of gene evolution typically treat each codon independently, quantifying gene evolution by summing over the constituent codons. In contrast, studies of protein evolution generally incorporate protein structure and interactions between amino acids explicitly. We investigate correlations in the evolution of codons as a function of their distance from each other along the protein coding sequence. This approach is motivated by the expectation that codons near each other in sequence often encode amino acids belonging to the same functional unit. Consequently, these amino acids are more likely to interact and/or experience similar selective regimes, introducing correlation between the evolution of the underlying codons. We find codon evolution in Drosophilids to be correlated over a characteristic length scale of ≈10 codons. Specifically, the presence of a non-synonymous substitution substantially increases the probability of further such substitutions nearby, particularly within that lineage. Further analysis suggests both functional interactions between amino acids and correlation in the strength of selection contribute to this effect. These findings are relevant for understanding the relative importance of different modes of selection, and particularly the role of epistasis, in gene and protein evolution.

Nonadaptive processes in primate and human evolution

Evolutionary biology has tended to focus on adaptive evolution by positive selection as the primum mobile of evolutionary trajectories in species while underestimating the importance of nonadaptive evolutionary processes. In this review, I describe evidence that suggests that primate and human evolution has been strongly influenced by nonadaptive processes, particularly random genetic drift and mutation. This is evidenced by three fundamental effects: a relative relaxation of selective constraints (i.e., purifying selection), a relative increase in the fixation of slightly deleterious mutations, and a general reduction in the efficacy of positive selection. These effects are observed in protein-coding, regulatory regions, and in gene expression data, as well as in an augmentation of fixation of large-scale mutations, including duplicated genes, mobile genetic elements, and nuclear mitochondrial DNA. The evidence suggests a general population-level explanation such as a reduction in effective population size (N(e)). This would have tipped the balance between the evolutionary forces of natural selection and random genetic drift toward genetic drift for variants having small selective effects. After describing these proximate effects, I describe the potential consequences of these effects for primate and human evolution. For example, an increase in the fixation of slightly deleterious mutations could potentially have led to an increase in the fixation rate of compensatory mutations that act to suppress the effects of slightly deleterious substitutions. The potential consequences of compensatory evolution for the evolution of novel gene functions and in potentially confounding the detection of positively selected genes are explored. The consequences of the passive accumulation of large-scale genomic mutations by genetic drift are unclear, though evidence suggests that new gene copies as well as insertions of transposable elements into genes can potentially lead to adaptive phenotypes. Finally, because a decrease in selective constraint at the genetic level is expected to have effects at the morphological level, I review studies that compare rates of morphological change in various mammalian and island populations where N(e) is reduced. Furthermore, I discuss evidence that suggests that craniofacial morphology in the Homo lineage has shifted from an evolutionary rate constrained by purifying selection toward a neutral evolutionary rate.

Fixation of a deleterious allele under mutation pressure and finite selection intensity

The mean fixation time of a deleterious mutant allele is studied beyond the diffusion approximation. As in Kimura's classical work [M. Kimura, Proc. Natl. Acad. Sci. USA. 77, 522 (1980)], that was motivated by the problem of fixation in the presence of amorphic or hypermorphic mutations, we consider a diallelic model at a single locus comprising a wild-type A and a mutant allele A' produced irreversibly from A at small uniform rate v. The relative fitnesses of the mutant homozygotes A'A', mutant heterozygotes A'A and wild-type homozygotes AA are 1-s, 1-h and 1, respectively, where it is assumed that v<<s. Here, we employ a WKB theory and directly treat the underlying Markov chain (formulated as a birth-death process) obeyed by the allele frequency (whose dynamics is prescribed by the Moran model). Importantly, this approach allows to accurately account for effects of large fluctuations. After a general description of the theory, we focus on the case of a deleterious mutant allele (i.e. s>0) and discuss three situations: when the mutant is (i) completely dominant (s=h); (ii) completely recessive (h=0), and (iii) semi-dominant (h=s/2). Our theoretical predictions for the mean fixation time and the quasi-stationary distribution of the mutant population in the coexistence state, are shown to be in excellent agreement with numerical simulations. Furthermore, when s is finite, we demonstrate that our results are superior to those of the diffusion theory, while the latter is shown to be an accurate approximation only when N(e)s(2)<<1, where N(e) is the effective population size.

Models of RNA virus evolution and their roles in vaccine design

Viruses are fast evolving pathogens that continuously adapt to the highly variable environments they live and reproduce in. Strategies devoted to inhibit virus replication and to control their spread among hosts need to cope with these extremely heterogeneous populations and with their potential to avoid medical interventions. Computational techniques such as phylogenetic methods have broadened our picture of viral evolution both in time and space, and mathematical modeling has contributed substantially to our progress in unraveling the dynamics of virus replication, fitness, and virulence. Integration of multiple computational and mathematical approaches with experimental data can help to predict the behavior of viral pathogens and to anticipate their escape dynamics. This piece of information plays a critical role in some aspects of vaccine development, such as viral strain selection for vaccinations or rational attenuation of viruses. Here we review several aspects of viral evolution that can be addressed quantitatively, and we discuss computational methods that have the potential to improve vaccine design.

Understanding infectious agents from an in silico perspective

Knowledge of infectious diseases now emerging from genomic, proteomic, epidemiological and clinical data can provide insights into the mechanisms of immune function, disease pathogenesis and epidemiology. Here, we describe how considerable advances in computational methods of data mining, mathematical modeling in epidemiology and simulation have been used to enhance our understanding of infectious agents and discuss their impact on the discovery of new therapeutics and controlling their spread.

Transition from positive to neutral in mutation fixation along with continuing rising fitness in thermal adaptive evolution

It remains to be determined experimentally whether increasing fitness is related to positive selection, while stationary fitness is related to neutral evolution. Long-term laboratory evolution in Escherichia coli was performed under conditions of thermal stress under defined laboratory conditions. The complete cell growth data showed common continuous fitness recovery to every 2°C or 4°C stepwise temperature upshift, finally resulting in an evolved E. coli strain with an improved upper temperature limit as high as 45.9°C after 523 days of serial transfer, equivalent to 7,560 generations, in minimal medium. Two-phase fitness dynamics, a rapid growth recovery phase followed by a gradual increasing growth phase, was clearly observed at diverse temperatures throughout the entire evolutionary process. Whole-genome sequence analysis revealed the transition from positive to neutral in mutation fixation, accompanied with a considerable escalation of spontaneous substitution rate in the late fitness recovery phase. It suggested that continually increasing fitness not always resulted in the reduction of genetic diversity due to the sequential takeovers by fit mutants, but caused the accumulation of a considerable number of mutations that facilitated the neutral evolution.

The detailed results of a two-year in vitro thermal adaptive evolution experiment are described. A laboratory-evolved E. coli strain with an improved upper temperature limit, as high as 45.9°C, was acquired after 523 days of serial transfer, equivalent to 7,560 generations, in nutrient-limited medium. The complete daily records of cell growth exhibited universal two-phase fitness dynamics, a rapid growth recovery phase followed by a gradual increasing growth phase, throughout the entire evolutionary process. Genome-sequence analysis not only showed considerable escalation of the spontaneous substitution rate, but also revealed the transition from positive to nearly neutral in mutation fixation. Particularly, even with the rising fitness of bacterial cells, neutrality was observed in molecular evolution. These observations suggested that a discrete evolutionary mode occurred in the continuous evolutionary route, linking Darwinian adaptive selection with Kimura's neutral evolution. Such transition from beneficial to neutral path may be adopted as an evolutionary strategy robust to rigorous environmental changes.

Robustness and evolvability in the functional anatomy of a PER-ARNT-SIM (PAS) domain

The robustness of proteins against point mutations implies that only a small subset of residues determines functional properties. We test this prediction using photoactive yellow protein (PYP), a 125-residue prototype of the PER-ARNT-SIM (PAS) domain superfamily of signaling proteins. PAS domains are defined by a small number of conserved residues of unknown function. We report high-throughput biophysical measurements on a complete Ala scan set of purified PYP mutants. The dataset of 1,193 values on active site properties, functional kinetics, stability, and production level reveals that 124 mutants retain the characteristic photocycle of PYP, but that the majority of substitutions significantly alter functional properties. Only 35% of substitutions that strongly affect function are located at the active site. Unexpectedly, most PAS-conserved residues are required for maintaining protein production. PAS domain activation often involves conformational changes in α-helices linked to the PAS core. However, the mechanism of transmission and kinetic regulation of allosteric structural changes from the PAS domain to these helices is not clear. The Ala scan data reveal interactions governing allosteric switching in PYP. The photocycle kinetics is significantly altered by substitutions at 58 positions and spans a 3,000-fold range. Nine residues that dock the N-terminal α-helices of PYP to its PAS core regulate signaling kinetics. Ile39 and Asn43 are identified as part of a mechanism for regulating allosteric switching that is conserved among PAS domains. These results show that PYP combines robustness with a high degree of evolvability and imply production level as an important factor in protein evolution.

The neutral theory of molecular evolution in the genomic era

The neutral theory of molecular evolution has been widely accepted and is the guiding principle for studying evolutionary genomics and the molecular basis of phenotypic evolution. Recent data on genomic evolution are generally consistent with the neutral theory. However, many recently published papers claim the detection of positive Darwinian selection via the use of new statistical methods. Examination of these methods has shown that their theoretical bases are not well established and often result in high rates of false-positive and false-negative results. When the deficiencies of these statistical methods are rectified, the results become largely consistent with the neutral theory. At present, genome-wide analyses of natural selection consist of collections of single-locus analyses. However, because phenotypic evolution is controlled by the interaction of many genes, the study of natural selection ought to take such interactions into account. Experimental studies of evolution will also be crucial.

Evolutionary history of Arabidopsis thaliana aminoacyl-tRNA synthetase dual-targeted proteins

Aminoacyl-transfer RNA (tRNA) synthetases (aaRS) are key players in translation and act early in protein synthesis by mediating the attachment of amino acids to their cognate tRNA molecules. In plants, protein synthesis may occur in three subcellular compartments (cytosol, mitochondria, and chloroplasts), which requires multiple versions of the protein to be correctly delivered to its proper destination. The organellar aaRS are nuclear encoded and equipped with targeting information at the N-terminal sequence, which enables them to be specifically translocated to their final location. Most of the aaRS families present organellar proteins that are dual targeted to mitochondria and chloroplasts. Here, we examine the dual targeting behavior of aaRS from an evolutionary perspective. Our results show that Arabidopsis thaliana aaRS sequences are a result of a horizontal gene transfer event from bacteria. However, there is no evident bias indicating one single ancestor (Cyanobacteria or Proteobacteria). The dual-targeted aaRS phylogenetic relationship was characterized into two different categories (paralogs and homologs) depending on the state recovered for both dual-targeted and cytosolic proteins. Taken together, our results suggest that the dual-targeted condition is a gain-of-function derived from gene duplication. Selection may have maintained the original function in at least one of the copies as the additional copies diverged.

Different phylogenomic approaches to resolve the evolutionary relationships among model fish species

Comparative genomics holds the promise to magnify the information obtained from individual genome sequencing projects, revealing common features conserved across genomes and identifying lineage-specific characteristics. To implement such a comparative approach, a robust phylogenetic framework is required to accurately reconstruct evolution at the genome level. Among vertebrate taxa, teleosts represent the second best characterized group, with high-quality draft genome sequences for five model species (Danio rerio, Gasterosteus aculeatus, Oryzias latipes, Takifugu rubripes, and Tetraodon nigroviridis), and several others are in the finishing lane. However, the relationships among the acanthomorph teleost model fishes remain an unresolved taxonomic issue. Here, a genomic region spanning over 1.2 million base pairs was sequenced in the teleost fish Dicentrarchus labrax. Together with genomic data available for the above fish models, the new sequence was used to identify unique orthologous genomic regions shared across all target taxa. Different strategies were applied to produce robust multiple gene and genomic alignments spanning from 11,802 to 186,474 amino acid/nucleotide positions. Ten data sets were analyzed according to Bayesian inference, maximum likelihood, maximum parsimony, and neighbor joining methods. Extensive analyses were performed to explore the influence of several factors (e.g., alignment methodology, substitution model, data set partitions, and long-branch attraction) on the tree topology. Although a general consensus was observed for a closer relationship between G. aculeatus (Gasterosteidae) and Di. labrax (Moronidae) with the atherinomorph O. latipes (Beloniformes) sister taxon of this clade, with the tetraodontiform group Ta. rubripes and Te. nigroviridis (Tetraodontiformes) representing a more distantly related taxon among acanthomorph model fish species, conflicting results were obtained between data sets and methods, especially with respect to the choice of alignment methodology applied to noncoding parts of the genomic region under study. This may limit the use of intergenic/noncoding sequences in phylogenomics until more robust alignment algorithms are developed.

Inference for the initial stage of domain shuffling: tracing the evolutionary fate of the PIPSL retrogene in hominoids

Domain shuffling has provided extraordinarily diverse functions to proteins. Nevertheless, how newly combined domains are coordinated to create novel functions remains a fundamental question of genetic and phenotypic evolution. Previously, we reported a unique mechanism of gene creation, whereby new combinations of functional domains are assembled from distinct genes at the RNA level, reverse transcribed, and integrated into the genome by the L1 retrotransposon. The novel gene PIPSL, created by the fusion of phosphatidylinositol-4-phosphate 5-kinase (PIP5K1A) and 26S proteasome subunit (S5a/PSMD4) genes, is specifically transcribed in human and chimpanzee testes. We present the first evidence for the translation of PIPSL in humans. The human PIPSL locus showed a low nucleotide diversity within 11 populations (125 individuals) compared with other genomic regions such as introns and overall chromosomes. It was equivalent to the average for coding sequences or exons from other genes, suggesting that human PIPSL has some function and is conserved among modern populations. Two linked amino acid-altering single-nucleotide polymorphisms were found in the PIPSL kinase domain of non-African populations. They are positioned in the vicinity of the substrate-binding cavity of the parental PIP5K1A protein and change the charge of both residues. The relatively rapid expansion of this haplotype might indicate a selective advantage for it in modern humans. We determined the evolutionary fate of PIPSL domains created by domain shuffling. During hominoid diversification, the S5a-derived domain was retained in all lineages, whereas the ubiquitin-interacting motif (UIM) 1 in the domain experienced critical amino acid replacements at an early stage, being conserved under subsequent high levels of nonsynonymous substitutions to UIM2 and other domains, suggesting that adaptive evolution diversified these functional compartments. Conversely, the PIP5K1A-derived domain is degenerated in gibbons and gorillas. These observations provide a possible scheme of domain shuffling in which the combined parental domains are not tightly linked in the novel chimeric protein, allowing for changes in their functional roles, leading to their fine-tuning. Selective pressure toward a novel function initially acted on one domain, whereas the other experienced a nearly neutral state. Over time, the latter also gained a new function or was degenerated.

Evolvability and Speed of Evolutionary Algorithms in Light of Recent Developments in Biology

Biological and artificial evolutionary systems exhibit varying degrees of evolvability and different rates of evolution. Such quantities can be affected by various factors. Here, we review some evolutionary mechanisms and discuss new developments in biology that can potentially improve evolvability or accelerate evolution in artificial systems. Biological notions are discussed to the degree they correspond to notions in Evolutionary Computation. We hope that the findings put forward here can be used to design computational models of evolution that produce significant gains in evolvability and evolutionary speed.

New insights into molecular evolution: prospects from the Barcode of Life Initiative (BOLI)

Geographic and temporal patterns of morphological and behavioral diversifications among species stimulated Darwin to propose a mechanism for evolutionary change through natural selection. Scientific developments have revealed an even more fundamental level of biological complexity: sequence variation in DNA. While genome projects yield spectacular insights into molecular evolution, they have targeted only a few species. In contrast, the Barcode of Life Initiative (BOLI) proposes a horizontal approach to genomics, examining short, standardized genome segments across the sweep of eukaryotic life, all 10 million species. BOLI will extend our understanding of evolution and speciation in varied ways. It will facilitate quantification of biological diversity by disclosing cryptic species and enabling a rapid survey of taxon diversity in groups that have hitherto received scant morphological examination. It will facilitate assignment of life history stages to known species and provide a first estimate of species ages. It will also reveal key features of the mitochondrial genome, because the evolutionary properties of barcodes relate to those in the mitochondrial genome as a whole, acting to flag taxonomic groups or species with unusual nucleotide composition or evolutionary rates. The growing volume of barcode records has revealed that sequence variability within species is generally much lower than divergence among species (barcoding gap), a pattern that occurs in diverse lineages, suggesting a pervasive evolutionary process. Low variability may reflect recurrent selective sweeps of favored mitochondrial variants propagating as single linkage units across species. If this hypothesis is substantiated, the implications are significant, particularly for our understanding of molecular evolution of mitochondrial DNA and its relationship with species delineation.

Degeneration of an intracellular ion channel in the primate lineage by relaxation of selective constraints

Ion channel genes are highly conserved and are rarely degenerated in the primate lineage leading to humans. So far, the only well-characterized ion channel known to be degenerated in primates is the plasma membrane transient receptor potential channel TRPC2, possibly due to changes in the pheromone signaling. Here, by analyzing the sequence data from ten primate species, we have determined the degeneration process of the TPC3 gene that encodes a member of the two-pore channel (TPC) family recently implicated in Ca(2+) release by nicotinic acid adenine dinucleotide phosphate from intracellular acidic stores in animals. We show that degeneration of TPC3 likely began in the common ancestors of Apes and Old World monkeys through a conserved inactivating mutation, followed by additional deleterious mutations resulting in the generation of a TPC3 pseudogene in the descendant catarrhine lineage. Located at a chromosome recombination hot spot, catarrhine TPC3 pseudogenes underwent a series of lineage-specific rearrangements, including exon deletion and duplication. In contrast, we identify near full-length TPC3 sequences in New World monkeys and Prosimians and show that the gene is subjected to strong purifying selection and therefore likely functional. Our data provide the first evidence for relaxed functional constraints for an intracellular ion channel in primates and shed novel insights into the evolution and regulation of Ca(2+) signaling in the primate lineage.

Lineage pattern, trans-species polymorphism, and selection pressure among the major lineages of feline MHC-DRB peptide-binding region

The long-term evolution of major histocompatibility complex (MHC) involves the birth-and-death process and independent divergence of loci during episodes punctuated by natural selection. Here, we investigated the molecular signatures of natural selection at exon-2 of MHC class II DRB gene which includes a part of the peptide-binding region (PBR) in seven of eight putative extant Felidae lineages. The DRB alleles in felids can be mainly divided into five lineages. Signatures of trans-species polymorphism among major allelic lineages indicate that balancing selection has maintained the MHC polymorphism for a long evolutionary time. Analysis based on maximum likelihood models of codon substitution revealed overall purifying selection acting on the feline DRB. Sites that have undergone positive selection and those that are under divergent selective pressure among lineages were detected and found to fall within the putative PBR. This study increased our understanding of the nature of selective forces acting on DRB during feline radiation.

A zinc-binding site by negative selection induces metallodrug susceptibility in an essential chaperonin

GroES is an indispensable chaperonin virtually found throughout all life forms. Consequently, mutations of this protein must be critically scrutinized by natural selection. Nevertheless, the homolog from a potentially virulent gastric pathogen, Helicobacter pylori, strikingly features a histidine/cysteine-rich C terminus that shares no significant homology with other family members. Additionally, three more (H45, C51, and C53) are uniquely present in its apical domain. The statistical analyses show that these residues may have originated from negative selection, presumably driven by either dependent or independent amino acid mutations. In the absence of the C-terminal metal-binding domain, the mutant protein still exhibits a substantial capacity for zinc binding in vivo. The biochemical properties of site-directed mutants indicate that H45, C51, and C53 make up an oxidation-sensitive zinc-binding site that may donate the bound metal to a zinc acceptor. Of interest, bismuth antiulcer drugs strongly bind at this site (K(d) of approximately 7 x 10(-26) M), replacing the bound zinc and consequently inducing the disruption of the quaternary structure. Because biological features by negative selection are usually inert to change during evolution, this study sheds light on a promising field whereby medicines can be designed or improved to specifically target the residues that uniquely evolved in pathogenic proteins so as to retard the emergence of drug resistance.

Gene duplication in the genome of parasitic Giardia lamblia

Background

Giardia are a group of widespread intestinal protozoan parasites in a number of vertebrates. Much evidence from G. lamblia indicated they might be the most primitive extant eukaryotes. When and how such a group of the earliest branching unicellular eukaryotes developed the ability to successfully parasitize the latest branching higher eukaryotes (vertebrates) is an intriguing question. Gene duplication has long been thought to be the most common mechanism in the production of primary resources for the origin of evolutionary novelties. In order to parse the evolutionary trajectory of Giardia parasitic lifestyle, here we carried out a genome-wide analysis about gene duplication patterns in G. lamblia.

Results

Although genomic comparison showed that in G. lamblia the contents of many fundamental biologic pathways are simplified and the whole genome is very compact, in our study 40% of its genes were identified as duplicated genes. Evolutionary distance analyses of these duplicated genes indicated two rounds of large scale duplication events had occurred in G. lamblia genome. Functional annotation of them further showed that the majority of recent duplicated genes are VSPs (Variant-specific Surface Proteins), which are essential for the successful parasitic life of Giardia in hosts. Based on evolutionary comparison with their hosts, it was found that the rapid expansion of VSPs in G. lamblia is consistent with the evolutionary radiation of placental mammals.

Conclusions

Based on the genome-wide analysis of duplicated genes in G. lamblia, we found that gene duplication was essential for the origin and evolution of Giardia parasitic lifestyle. The recent expansion of VSPs uniquely occurring in G. lamblia is consistent with the increment of its hosts. Therefore we proposed a hypothesis that the increment of Giradia hosts might be the driving force for the rapid expansion of VSPs.

Local adaptation maintains clinal variation in melanin-based coloration of European barn owls (Tyto alba)

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin-based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish-brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.

Lineage-specific differences in the amino acid substitution process

In Darwinian evolution, mutations occur approximately at random in a gene, turned into amino acid mutations by the genetic code. Some mutations are fixed to become substitutions and some are eliminated from the population. Partitioning pairs of closely related species with complete genome sequences by average population size of each pair, we looked at the substitution matrices generated for these partitions and compared the substitution patterns between species. We estimated a population genetic model that relates the relative fixation probabilities of different types of mutations to the selective pressure and population size. Parameterizations of the average and distribution of selective pressures for different amino acid substitution types in different population size comparisons were generated with a Bayesian framework. We found that partitions in population size as well as in substitution type are required to explain the substitution data. Selection coefficients were found to decrease with increasingly radical amino acid substitution and with increasing effective population size. To further explore the role of underlying processes in amino acid substitution, we analyzed embryophyte (plant) gene families from TAED (The Adaptive Evolution Database), where solved structures for at least one member exist in the Protein Data Bank. Using PAML, we assigned branches to three categories: strong negative selection, moderate negative selection/neutrality, and positive diversifying selection. Focusing on the first and third categories, we identified sites changing along gene family lineages and observed the spatial patterns of substitution. Selective sweeps were expected to create primary sequence clustering under positive diversifying selection. Co-evolution through direct physical interaction was expected to cause tertiary structural clustering. Under both positive and negative selection, the substitution patterns were found to be nonrandom. Under positive diversifying selection, significant independent signals were found for primary and tertiary sequence clustering, suggesting roles for both selective sweeps and direct physical interaction. Under strong negative selection, the signals were not found to be independent. All together, a complex interplay of population genetic and protein thermodynamics forces is suggested.

How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings

Major histocompatibility complex (MHC) genes have been put forward as a model for studying how genetic diversity is maintained in wild populations. Pathogen-mediated selection (PMS) is believed to generate the extraordinary levels of MHC diversity observed. However, establishing the relative importance of the three proposed mechanisms of PMS (heterozygote advantage, rare-allele advantage and fluctuating selection) has proved extremely difficult. Studies have attempted to differentiate between mechanisms of PMS using two approaches: (i) comparing MHC diversity with that expected under neutrality and (ii) relating MHC diversity to pathogen regime. Here, we show that in many cases the same predictions arise from the different mechanisms under these approaches, and that most studies that have inferred one mechanism of selection have not fully considered the alternative explanations. We argue that, while it may be possible to demonstrate that particular mechanisms of PMS are occurring, resolving their relative importance within a system is probably impossible. A more realistic target is to continue to demonstrate when and where the different mechanisms of PMS occur, with the aim of determining their relative importance across systems. We put forward what we believe to be the most promising approaches that will allow us to progress towards achieving this.

Immunogenicity characterization of two ancient wheat α-gliadin peptides related to coeliac disease

The immunogenic potential of α-gliadin protein from two ancient wheats was studied with reference to coeliac disease. To this aim we investigated Graziella Ra® and Kamut® (the latter is considered an ancient relative of modern durum wheat) in comparison to four durum wheat accessions (Senatore Cappelli, Flaminio, Grazia and Svevo). ELISA and Western Blot analyses - carried out by two monoclonal antibodies raised against the α-gliadin peptides p31-49 (LGQQQPFPQQPYPQPQPF) and p56-75 (LQLQPFPQPQLPYPQPQLPY) containing a core region (underlined) reported to be toxic for coeliac patients - always showed an antibody-antigen positive reaction. For all accessions, an α-gliadin gene has also been cloned and sequenced. Deduced amino acid sequences constantly showed the toxic motifs. In conclusion, we strongly recommend that coeliac patients should avoid consuming Graziella Ra® or Kamut®. In fact their α-gliadin not only is as toxic as one of the other wheat accessions, but also occurs in greater amount, which is in line with the higher level of proteins in ancient wheats when compared to modern varieties.

Polymorphism of the ovine keratin-associated protein 1-4 gene (KRTAP1-4)

Keratin-associated proteins (KAPs) are one of the main structural components of the wool fibre and form a semi-rigid matrix in which the keratin intermediate filaments are embedded. Variation in the KAP genes may affect the structure of KAPs and hence wool characteristics. In this study, we used PCR-SSCP to analyse ovine KRTAP1-4 (previously B2D), a gene encoding a member of the KAP1-x family. Nine different PCR-SSCP patterns were detected in the 320 sheep that were analysed. Either one or a combination of two patterns was observed for each sheep, which was consistent with these sheep being either homozygous or heterozygous for this gene. DNA sequencing revealed that these patterns represent nine different DNA sequences. All of these sequences were unique, but shared a high homology with the published ovine KRTAP1-4 sequence, suggesting that these sequences represent allelic variants of KRTAP1-4. There were a total of 14 single nucleotide polymorphisms (SNPs) identified and these SNPs tended to be clustered in two regions. Of the 13 SNPs found in the coding region, nine were non-synonymous SNPs and would result in amino acid changes. The variation detected here may have an impact on the structure of KAP1-4 and hence affect wool traits.

Anticancer drug discovery in the future: an evolutionary perspective

Identification of agents that are pharmacologically active against human cancer has depended largely on the screening of natural products and their analogs. Many anticancer drugs have been discovered fortuitously through random investigation of organisms; indeed, serendipity remains important in anticancer drug discovery. Although it is broadly accepted that cancers comprise an evolutionary microcosm, this idea has not been advanced to understand and control carcinogenic progression. Here, we address anticancer drug discovery from an evolutionary perspective and present a series of case studies that demonstrate that the rate of anticancer drug discovery can be increased greatly by targeted screening of natural compounds from ancient species.

A universal nonmonotonic relationship between gene compactness and expression levels in multicellular eukaryotes

Analysis of gene architecture and expression levels of four organisms, Homo sapiens, Caenorhabditis elegans, Drosophila melanogaster, and Arabidopsis thaliana, reveals a surprising, nonmonotonic, universal relationship between expression level and gene compactness. With increasing expression level, the genes tend at first to become longer but, from a certain level of expression, they become more and more compact, resulting in an approximate bell-shaped dependence. There are two leading hypotheses to explain the compactness of highly expressed genes. The selection hypothesis predicts that gene compactness is predominantly driven by the level of expression, whereas the genomic design hypothesis predicts that expression breadth across tissues is the driving force. We observed the connection between gene expression breadth in humans and gene compactness to be significantly weaker than the connection between expression level and compactness, a result that is compatible with the selection hypothesis but not the genome design hypothesis. The initial gene elongation with increasing expression level could be explained, at least in part, by accumulation of regulatory elements enhancing expression, in particular, in introns. This explanation is compatible with the observed positive correlation between intron density and expression level of a gene. Conversely, the trend toward increasing compactness for highly expressed genes could be caused by selection for minimization of energy and time expenditure during transcription and splicing and for increased fidelity of transcription, splicing, and/or translation that is likely to be particularly critical for highly expressed genes. Regardless of the exact nature of the forces that shape the gene architecture, we present evidence that, at least, in animals, coding and noncoding parts of genes show similar architectonic trends.

Neutral evolution on mammalian protein surfaces

Because of their low effective population sizes, natural selection is expected to have reduced effectiveness in organisms such as mammals. By comparing the amino acid substitution rates between mammalian protein surfaces and interiors, it was found that almost a third of the proteins surveyed failed to reject the null hypothesis of neutral substitutions among surface residues. Proteins with such partly neutral evolution nonetheless have no fewer protein interactions than do other proteins. I suggest that natural selection can function to preserve protein interactions without requiring strict conservation of the individual residue contacts that impart those interactions.

Evolution of plant MADS box transcription factors: evidence for shifts in selection associated with early angiosperm diversification and concerted gene duplications

Phylogenomic analyses show that gene and genome duplication events have led to the diversification of transcription factor gene families throughout the evolutionary history of land plants and that gene duplications have played an important role in shaping regulatory networks influencing key phenotypic characters including floral development and flowering time. A molecular evolutionary investigation of the mode and tempo of selection acting on the angiosperm MADS box AP1/SQUA, AP3/PI, AG/AGL11, and SEP gene subfamilies revealed site-specific patterns of shifting evolutionary constraint throughout angiosperm history. Specific positions in the four canonical MADS box gene regions, especially K domains and C-terminal regions of all four of these MADS box gene subfamilies exhibited clade-specific shifts in selective constraint following concerted duplication events. Moreover, the frequency of site-specific shifts in constraint was correlated with gene duplications and early angiosperm diversification. We hypothesize that coevolution among interacting MADS box proteins may be responsible for simultaneous increases in the ratio of nonsynonymous to synonymous substitutions (d(N)/d(S) = omega) early in angiosperm history and following concerted duplication events.

In the light of directed evolution: pathways of adaptive protein evolution

Directed evolution is a widely-used engineering strategy for improving the stabilities or biochemical functions of proteins by repeated rounds of mutation and selection. These experiments offer empirical lessons about how proteins evolve in the face of clearly-defined laboratory selection pressures. Directed evolution has revealed that single amino acid mutations can enhance properties such as catalytic activity or stability and that adaptation can often occur through pathways consisting of sequential beneficial mutations. When there are no single mutations that improve a particular protein property experiments always find a wealth of mutations that are neutral with respect to the laboratory-defined measure of fitness. These neutral mutations can open new adaptive pathways by at least 2 different mechanisms. Functionally-neutral mutations can enhance a protein's stability, thereby increasing its tolerance for subsequent functionally beneficial but destabilizing mutations. They can also lead to changes in "promiscuous" functions that are not currently under selective pressure, but can subsequently become the starting points for the adaptive evolution of new functions. These lessons about the coupling between adaptive and neutral protein evolution in the laboratory offer insight into the evolution of proteins in nature.

Genomic patterns of adaptive divergence between chromosomally differentiated sunflower species

Understanding the genetic mechanisms of speciation and basis of species differences is among the most important challenges in evolutionary biology. Two questions of particular interest are what roles divergent selection and chromosomal differentiation play in these processes. A number of recently proposed theories argue that chromosomal rearrangements can facilitate the development and maintenance of reproductive isolation and species differences by suppressing recombination within rearranged regions. Reduced recombination permits the accumulation of alleles contributing to isolation and adaptive differentiation and protects existing differences from the homogenizing effects of introgression between incipient species. Here, we examine patterns of genetic diversity and divergence in rearranged versus collinear regions in two widespread, extensively hybridizing sunflower species, Helianthus annuus and Helianthus petiolaris, using sequence data from 77 loci distributed throughout the genomes of the two species. We find weak evidence for increased genetic divergence near chromosomal break points but not within rearranged regions overall. We find no evidence for increased rates of adaptive divergence on rearranged chromosomes; in fact, collinear chromosomes show a far greater excess of fixed amino acid differences between the two species. A comparison with a third sunflower species indicates that much of the nonsynonymous divergence between H. annuus and H. petiolaris probably occurred during or soon after their formation. Our results suggest a limited role for chromosomal rearrangements in genetic divergence, but they do document substantial adaptive divergence and provide further evidence of how species integrity and genetic identity can be maintained at many loci in the face of extensive hybridization and gene flow.

Reliabilities of identifying positive selection by the branch-site and the site-prediction methods

Natural selection operating in protein-coding genes is often studied by examining the ratio (omega) of the rates of nonsynonymous to synonymous nucleotide substitution. The branch-site method (BSM) based on a likelihood ratio test is one of such tests to detect positive selection for a predetermined branch of a phylogenetic tree. However, because the number of nucleotide substitutions involved is often very small, we conducted a computer simulation to examine the reliability of BSM in comparison with the small-sample method (SSM) based on Fisher's exact test. The results indicate that BSM often generates false positives compared with SSM when the number of nucleotide substitutions is approximately 80 or smaller. Because the omega value is also used for predicting positively selected sites, we examined the reliabilities of the site-prediction methods, using nucleotide sequence data for the dim-light and color vision genes in vertebrates. The results showed that the site-prediction methods have a low probability of identifying functional changes of amino acids experimentally determined and often falsely identify other sites where amino acid substitutions are unlikely to be important. This low rate of predictability occurs because most of the current statistical methods are designed to identify codon sites with high omega values, which may not have anything to do with functional changes. The codon sites showing functional changes generally do not show a high omega value. To understand adaptive evolution, some form of experimental confirmation is necessary.

Inferring human stem cell behaviour from epigenetic drift

Although conventional experimental manipulations are impractical, it may be possible to infer human stem cell fates by 'reading' histories recorded within their genomes. Genomes are almost perfect copies of copies, and ancestries may be surreptitiously recorded by replication errors that inevitably accumulate. The greater the number of divisions, the greater the number of replication errors ('a molecular clock hypothesis'). Mutations rarely occur during a lifetime, but DNA methylation patterns are also copied after DNA replication and measurably drift with ageing at certain CpG sites in mitotic tissues, such as the colon. Such passenger methylation pattern variation may effectively function as 'epigenetic' somatic cell mitotic clocks. Replication errors can only accumulate in long-lived stem cell lineages, so methylation pattern drift largely records stem cell behaviour. How methylation patterns may encode stem cell ancestries is illustrated with two types of small reproductive units--colon crypt niches with continuous genealogies, and hair follicles with punctuated genealogies. Potentially, the genealogy of any human cell may be inferred by 'reading' its genome.

From fitness landscapes to seascapes: non-equilibrium dynamics of selection and adaptation

Evolution is a quest for innovation. Organisms adapt to changing natural selection by evolving new phenotypes. Can we read this dynamics in their genomes? Not every mutation under positive selection responds to a change in selection: beneficial changes also occur at evolutionary equilibrium, repairing previous deleterious changes and restoring existing functions. Adaptation, by contrast, is viewed here as a non-equilibrium phenomenon: the genomic response to time-dependent selection. Our approach extends the static concept of fitness landscapes to dynamic fitness seascapes. It shows that adaptation requires a surplus of beneficial substitutions over deleterious ones. Here, we focus on the evolution of yeast and Drosophila genomes, providing examples where adaptive evolution can and cannot be inferred, despite the presence of positive selection.

Quantification of adaptive evolution of genes expressed in avian brain and the population size effect on the efficacy of selection

Whether protein evolution is mainly due to fixation of beneficial alleles by positive selection or to random genetic drift has remained a contentious issue over the years. Here, we use two genomewide polymorphism data sets collected from chicken populations, together with divergence data from >5,000 chicken-zebra finch gene orthologs expressed in brain, to assess the amount of adaptive evolution in protein-coding genes of birds. First, we show that estimates of the fixation index (FI, the ratio of fixed nonsynonymous-to-synonymous changes over the ratio of the corresponding polymorphisms) are highly dependent on the character of the underlying data sets. Second, by using polymorphism data from high-frequency alleles, to avoid the confounding effect of slightly deleterious mutations segregating at low frequency, we estimate that about 20% of amino acid changes have been brought to fixation through positive selection during avian evolution. This estimate is intermediate to that obtained in humans (lower) and flies as well as bacteria (higher), and is consistent with population genetics theory that stipulates a positive relationship between the efficiency of selection and the effective population size. Further, by comparing the FIs for common and all alleles, we estimate that approximately 20% of nonsynonymous variation segregating in chicken populations represent slightly deleterious mutations, which is less than in Drosophila. Overall, these results highlight the link between the effective population size and positive as well as negative selection.

Nucleotide variation in the NCED3 region of Arabidopsis thaliana and its association study with abscisic acid content under drought stress

Drought tolerance is a comprehensive quantitative trait that is being understood further at the molecular genetic level. Abscisic acid (ABA) is the main drought-induced hormone that regulates the expression of many genes related to drought responses. 9-cis-epoxycarotenoid dioxygenase (NCED3) is thought to be a key enzyme in ABA biosynthesis. In this paper, we measured the ABA content increase under drought stress, and sequenced and compared the sequence of AtNCED3 among 22 Arabidopsis thaliana accessions. The results showed that the fold of ABA content increase under drought stress was highly variable among these accessions. High density single nucleotide polymorphism (SNP) and insertion/deletion (indel) were found in the AtNCED3 region, on average one SNP per 87.4 bp and one indel per 502 bp. Nucleotide diversity was significantly lower in the coding region than that in non-coding regions. The results of an association study with anova analysis suggested that the 274th site (P<-->S) and the 327th site (P<-->R) amino acid variations might be the cause of ABA content increase of 163av accession under drought stress.

Non-random pre-transcriptional evolution in HIV-1. A refutation of the foundational conditions for neutral evolution

The complete base sequence of HIV-1 virus and GP120 ENV gene were analyzed to establish their distance to the expected neutral random sequence. An especial methodology was devised to achieve this aim. Analyses included: a) proportion of dinucleotides (signatures); b) homogeneity in the distribution of dinucleotides and bases (isochores) by dividing both segments in ten and three sub-segments, respectively; c) probability of runs of bases and No-bases according to the Bose-Einstein distribution. The analyses showed a huge deviation from the random distribution expected from neutral evolution and neutral-neighbor influence of nucleotide sites. The most significant result is the tremendous lack of CG dinucleotides (p < 10^-50 ), a selective trait of eukaryote and not of single stranded RNA virus genomes. Results not only refute neutral evolution and neutral neighbor influence, but also strongly indicate that any base at any nucleotide site correlates with all the viral genome or sub-segments. These results suggest that evolution of HIV-1 is pan-selective rather than neutral or nearly neutral.

Neutralism and selectionism: a network-based reconciliation

Neutralism and selectionism are extremes of an explanatory spectrum for understanding patterns of molecular evolution and the emergence of evolutionary innovation. Although recent genome-scale data from protein-coding genes argue against neutralism, molecular engineering and protein evolution data argue that neutral mutations and mutational robustness are important for evolutionary innovation. Here I propose a reconciliation in which neutral mutations prepare the ground for later evolutionary adaptation. Key to this perspective is an explicit understanding of molecular phenotypes that has only become accessible in recent years.

The evolution of animal chemosensory receptor gene repertoires: roles of chance and necessity

Chemosensory receptors are essential for the survival of organisms that range from bacteria to mammals. Recent studies have shown that the numbers of functional chemosensory receptor genes and pseudogenes vary enormously among the genomes of different animal species. Although much of the variation can be explained by the adaptation of organisms to different environments, it has become clear that a substantial portion is generated by genomic drift, a random process of gene duplication and deletion. Genomic drift also generates a substantial amount of copy-number variation in chemosensory receptor genes within species. It seems that mutation by gene duplication and inactivation has important roles in both the adaptive and non-adaptive evolution of chemosensation.