A new method for estimating nonsynonymous substitutions and its applications to detecting positive selection

An Improved Codon Modeling Approach for Accurate Estimation of the Mutation Bias

Phylogenetic codon models are routinely used to characterize selective regimes in coding sequences. Their parametric design, however, is still a matter of debate, in particular concerning the question of how to account for differing nucleotide frequencies and substitution rates. This problem relates to the fact that nucleotide composition in protein-coding sequences is the result of the interactions between mutation and selection. In particular, because of the structure of the genetic code, the nucleotide composition differs between the three coding positions, with the third position showing a more extreme composition. Yet, phylogenetic codon models do not correctly capture this phenomenon and instead predict that the nucleotide composition should be the same for all three positions. Alternatively, some models allow for different nucleotide rates at the three positions, an approach conflating the effects of mutation and selection on nucleotide composition. In practice, it results in inaccurate estimation of the strength of selection. Conceptually, the problem comes from the fact that phylogenetic codon models do not correctly capture the fixation bias acting against the mutational pressure at the mutation–selection equilibrium. To address this problem and to more accurately identify mutation rates and selection strength, we present an improved codon modeling approach where the fixation rate is not seen as a scalar, but as a tensor. This approach gives an accurate representation of how mutation and selection oppose each other at equilibrium and yields a reliable estimate of the mutational process, while disentangling the mean fixation probabilities prevailing in different mutational directions.

Molecular Evolution in Small Steps under Prevailing Negative Selection: A Nearly Universal Rule of Codon Substitution

The widely accepted view that evolution proceeds in small steps is based on two premises: 1) negative selection acts strongly against large differences and 2) positive selection favors small-step changes. The two premises are not biologically connected and should be evaluated separately. We now extend a previous approach to studying codon evolution in the entire genome. Codon substitution rate is a function of the physicochemical distance between amino acids (AAs), equated with the step size of evolution. Between nine pairs of closely related species of plants, invertebrates, and vertebrates, the evolutionary rate is strongly and negatively correlated with a set of AA distances (Δ_U, scaled to [0, 1]). Δ_U, a composite measure of evolutionary rates across diverse taxa, is influenced by almost all of the 48 physicochemical properties used here. The new analyses reveal a crucial trend hidden from previous studies: Δ_U is strongly correlated with the evolutionary rate (R² > 0.8) only when the genes are predominantly under negative selection. Because most genes in most taxa are strongly constrained by negative selection, Δ_U has indeed appeared to be a nearly universal measure of codon evolution. In conclusion, molecular evolution at the codon level generally takes small steps due to the prevailing negative selection. Whether positive selection may, or may not, follow the small-step rule is addressed in a companion study.

Molecular Evolution in Large Steps-Codon Substitutions under Positive Selection

Molecular evolution is believed to proceed in small steps. The step size can be defined by a distance reflecting physico-chemical disparities between amino acid (AA) pairs that can be exchanged by single 1-bp mutations. We show that AA substitution rates are strongly and negatively correlated with this distance but only when positive selection is relatively weak. We use the McDonald and Kreitman test to separate the influences of positive and negative selection. While negative selection is indeed stronger on AA substitutions generating larger changes in chemical properties of AAs, positive selection operates by different rules. For 65 of the 75 possible pairs, positive selection is comparable in strength regardless of AA distance. However, the ten pairs under the strongest positive selection all exhibit large leaps in chemical properties. Five of the ten pairs are shared between Drosophila and Hominoids, thus hinting at a common but modest biochemical basis of adaptation across taxa. The hypothesis that adaptive changes often take large functional steps will need to be extensively tested. If validated, molecular models will need to better integrate positive and negative selection in the search for adaptive signal.

Amino acid exchangeabilities vary across the tree of life

Different amino acid pairs have drastically different relative exchangeabilities (REs), and accounting for this variation is an important and common practice in inferring phylogenies, testing selection, and predicting mutational effects, among other analyses. In all such endeavors, REs have been generally considered invariant among species; this assumption, however, has not been scrutinized. Using maximum likelihood to analyze 180 genome sequences, we estimated REs from 90 clades representing all three domains of life, and found numerous instances of substantial between-clade differences in REs. REs show more differences between orthologous proteins of different clades than unrelated proteins of the same clade, suggesting that REs are genome-wide, clade-specific features, probably a result of proteome-wide evolutionary changes in the physicochemical environments of amino acid residues. The discovery of among-clade RE variations cautions against assuming constant REs in various analyses and demonstrates a higher-than-expected complexity in mechanisms of proteome evolution.

Transcriptome analyses provide insights into the phylogeny and adaptive evolution of the mangrove fern genus Acrostichum

The mangrove fern genus Acrostichum grows in the extremely unstable marine intertidal zone under harsh conditions, such as high salt concentrations, tidal rhythms and long-term climate changes. To explore the phylogenetic relationships and molecular mechanisms underlying adaptations in this genus, we sequenced the transcriptomes of two species of Acrostichum, A. aureum and A. speciosum, as well as a species in the sister genus, Ceratopteris thalictroides. We obtained 47,517, 36,420 and 60,823 unigenes for the three ferns, of which 24.39-45.63% were annotated using public databases. The estimated divergence time revealed that Acrostichum adapted to the coastal region during the late Cretaceous, whereas the two mangrove ferns from the Indo West-Pacific (IWP) area diverged more recently. Two methods (the modified branch-site model and the K_h method) were used to identify several positively selected genes, which may contribute to differential adaptation of the two Acrostichum species to different light and salt conditions. Our study provides abundant transcriptome data and new insights into the evolution and adaptations of mangrove ferns in the inhospitable intertidal zone.

PSP: rapid identification of orthologous coding genes under positive selection across multiple closely related prokaryotic genomes

Background

With genomic sequences of many closely related bacterial strains made available by deep sequencing, it is now possible to investigate trends in prokaryotic microevolution. Positive selection is a sub-process of microevolution, in which a particular mutation is favored, causing the allele frequency to continuously shift in one direction. Wide scanning of prokaryotic genomes has shown that positive selection at the molecular level is much more frequent than expected. Genes with significant positive selection may play key roles in bacterial adaption to different environmental pressures. However, selection pressure analyses are computationally intensive and awkward to configure.

Results

Here we describe an open access web server, which is designated as PSP (Positive Selection analysis for Prokaryotic genomes) for performing evolutionary analysis on orthologous coding genes, specially designed for rapid comparison of dozens of closely related prokaryotic genomes. Remarkably, PSP facilitates functional exploration at the multiple levels by assignments and enrichments of KO, GO or COG terms. To illustrate this user-friendly tool, we analyzed Escherichia coli and Bacillus cereus genomes and found that several genes, which play key roles in human infection and antibiotic resistance, show significant evidence of positive selection. PSP is freely available to all users without any login requirement at: http://db-mml.sjtu.edu.cn/PSP/.

Conclusions

PSP ultimately allows researchers to do genome-scale analysis for evolutionary selection across multiple prokaryotic genomes rapidly and easily, and identify the genes undergoing positive selection, which may play key roles in the interactions of host-pathogen and/or environmental adaptation.

Transcriptome sequencing and annotation for the Jamaican fruit bat (Artibeus jamaicensis)

The Jamaican fruit bat (Artibeus jamaicensis) is one of the most common bats in the tropical Americas. It is thought to be a potential reservoir host of Tacaribe virus, an arenavirus closely related to the South American hemorrhagic fever viruses. We performed transcriptome sequencing and annotation from lung, kidney and spleen tissues using 454 and Illumina platforms to develop this species as an animal model. More than 100,000 contigs were assembled, with 25,000 genes that were functionally annotated. Of the remaining unannotated contigs, 80% were found within bat genomes or transcriptomes. Annotated genes are involved in a broad range of activities ranging from cellular metabolism to genome regulation through ncRNAs. Reciprocal BLAST best hits yielded 8,785 sequences that are orthologous to mouse, rat, cattle, horse and human. Species tree analysis of sequences from 2,378 loci was used to achieve 95% bootstrap support for the placement of bat as sister to the clade containing horse, dog, and cattle. Through substitution rate estimation between bat and human, 32 genes were identified with evidence for positive selection. We also identified 466 immune-related genes, which may be useful for studying Tacaribe virus infection of this species. The Jamaican fruit bat transcriptome dataset is a resource that should provide additional candidate markers for studying bat evolution and ecology, and tools for analysis of the host response and pathology of disease.

Identification of HCV protease inhibitor resistance mutations by selection pressure-based method

A major challenge to successful antiviral therapy is the emergence of drug-resistant viruses. Recent studies have developed several automated analyses of HIV sequence polymorphism based on calculations of selection pressure (K_a/K_s) to predict drug resistance mutations. Similar resistance analysis programs for HCV inhibitors are not currently available. Taking advantage of the recently available sequence data of patient HCV samples from a Phase II clinical study of protease inhibitor boceprevir, we calculated the selection pressure for all codons in the HCV protease region (amino acid 1–181) to identify potential resistance mutations. The correlation between mutations was also calculated to evaluate linkage between any two mutations. Using this approach, we identified previously known major resistant mutations, including a recently reported mutation V55A. In addition, a novel mutation V158I was identified, and we further confirmed its resistance to boceprevir in protease enzyme and replicon assay. We also extended the approach to analyze potential interactions between individual mutations and identified three pairs of correlated changes. Our data suggests that selection pressure-based analysis and correlation mapping could provide useful tools to analyze large amount of sequencing data from clinical samples and to identify new drug resistance mutations as well as their linkage and correlations.

Highly conserved regimes of neighbor-base-dependent mutation generated the background primary-structural heterogeneities along vertebrate chromosomes

The content of guanine+cytosine varies markedly along the chromosomes of homeotherms and great effort has been devoted to studying this heterogeneity and its biological implications. Already before the DNA-sequencing era, however, it was established that the dinucleotides in the DNA of mammals in particular, and of most organisms in general, show striking over- and under-representations that cannot be explained by the base composition. Here we show that in the coding regions of vertebrates both GC content and codon occurrences are strongly correlated with such "motif preferences" even though we quantify the latter using an index that is not affected by the base composition, codon usage, and protein-sequence encoding. These correlations are likely to be the result of the long-term shaping of the primary structure of genic and non-genic DNA by a regime of mutation of which central features have been maintained by natural selection. We find indeed that these preferences are conserved in vertebrates even more rigidly than codon occurrences and we show that the occurrence-preference correlations are stronger in intronic and non-genic DNA, with the R(2)s reaching 99% when GC content is approximately 0.5. The mutation regime appears to be characterized by rates that depend markedly on the bases present at the site preceding and at that following each mutating site, because when we estimate such rates of neighbor-base-dependent mutation (NBDM) from substitutions retrieved from alignments of coding, intronic, and non-genic mammalian DNA sorted and grouped by GC content, they suffice to simulate DNA sequences in which motif occurrences and preferences as well as the correlations of motif preferences with GC content and with motif occurrences, are very similar to the mammalian ones. The best fit, however, is obtained with NBDM regimes lacking strand effects, which indicates that over the long term NBDM switches strands in the germline as one would expect for effects due to loosely contained background transcription. Finally, we show that human coding regions are less mutable under the estimated NBDM regimes than under matched context-independent mutation and that this entails marked differences between the spectra of amino-acid mutations that either mutation regime should generate. In the Discussion we examine the mechanisms likely to underlie NBDM heterogeneity along chromosomes and propose that it reflects how the diversity and activity of lesion-bypass polymerases (LBPs) track the landscapes of scheduled and non-scheduled genome repair, replication, and transcription during the cell cycle. We conclude that the primary structure of vertebrate genic DNA at and below the trinucleotide level has been governed over the long term by highly conserved regimes of NBDM which should be under direct natural selection because they alter drastically missense-mutation rates and hence the somatic and the germline mutational loads. Therefore, the non-coding DNA of vertebrates may have been shaped by NBDM only epiphenomenally, with non-genic DNA being affected mainly when found in the proximity of genes.

Null mutations in human and mouse orthologs frequently result in different phenotypes

One-to-one orthologous genes of relatively closely related species are widely assumed to have similar functions and cause similar phenotypes when deleted from the genome. Although this assumption is the foundation of comparative genomics and the basis for the use of model organisms to study human biology and disease, its validity is known only from anecdotes rather than from systematic examination. Comparing documented phenotypes of null mutations in humans and mice, we find that >20% of human essential genes have nonessential mouse orthologs. These changes of gene essentiality appear to be associated with adaptive evolution at the protein-sequence, but not gene-expression, level. Proteins localized to the vacuole, a cellular compartment for waste management, are highly enriched among essentiality-changing genes. It is probable that the evolution of the prolonged life history in humans required enhanced waste management for proper cellular function until the time of reproduction, which rendered these vacuole proteins essential and generated selective pressures for their improvement. If our gene sample represents the entire genome, our results would mean frequent changes of phenotypic effects of one-to-one orthologous genes even between relatively closely related species, a possibility that should be considered in comparative genomic studies and in making cross-species inferences of gene function and phenotypic effect.

The new mutation theory of phenotypic evolution

Recent studies of developmental biology have shown that the genes controlling phenotypic characters expressed in the early stage of development are highly conserved and that recent evolutionary changes have occurred primarily in the characters expressed in later stages of development. Even the genes controlling the latter characters are generally conserved, but there is a large component of neutral or nearly neutral genetic variation within and between closely related species. Phenotypic evolution occurs primarily by mutation of genes that interact with one another in the developmental process. The enormous amount of phenotypic diversity among different phyla or classes of organisms is a product of accumulation of novel mutations and their conservation that have facilitated adaptation to different environments. Novel mutations may be incorporated into the genome by natural selection (elimination of preexisting genotypes) or by random processes such as genetic and genomic drift. However, once the mutations are incorporated into the genome, they may generate developmental constraints that will affect the future direction of phenotypic evolution. It appears that the driving force of phenotypic evolution is mutation, and natural selection is of secondary importance.

Adaptive evolution in humans revealed by the negative correlation between the polymorphism and fixation phases of evolution

The selective forces acting on amino acid substitutions may be different in the two phases of molecular evolution: polymorphism and fixation. Negative selection and genetic drift may dominate the first phase, whereas positive selection may become much more significant in the second phase. However, the conventional dichotomy of synonymous vs. nonsynonymous changes does not offer the resolution needed to study the dynamics of these two phases. Following previously published methods, we separated amino acid changes into 75 elementary types (1-bp substitution between their respective codons). The likelihood of each type of amino acid change becoming polymorphic (PI, which stands for "polymorphic index"), relative to synonymous changes, can then be calculated. Similarly, the likelihood of fixation (FI, for "fixation index"), conditional on common polymorphisms, is also calculated. Using Perlegen and HapMap data on human polymorphisms and the chimpanzee sequences as the outgroup, we compared the evolutionary dynamics of the 75 elementary changes in the two phases. We found a strong "L-shaped" negative correlation (P < 0.001) between FI and PI. Only those changes with low PIs show FI > 1, which is often a signature of adaptive evolution. These patterns suggest that negative and positive selection operate more effectively on the same set of amino acid changes and that approximately 10-13% of amino acid substitutions between humans and chimpanzee may be adaptive.