Widespread positive selection in synonymous sites of mammalian genes

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.     

Fungal Smn and Spf30 homologues are mainly present in filamentous fungi and genomes with many introns: implications for spinal muscular atrophy

Spinal muscular atrophy is an important rare genetic disease characterized by the loss of motor neurons, where the main gene responsible is smn1. Orthologous genes have only been characterized in a single fungal genome: Schizosaccharomyces pombe. We have searched for putative SMN orthologues in publically available fungal genomes, finding that they are predominately present in filamentous fungi. SMN binding partners and the SPF30 SMN paralogue, which are all involved in mRNA splicing, were found to be present in a similar but non-identical subset of fungal genomes. The Saccharomycces cerevisiae yeast genome contains neither smn1 orthologues nor paralogues and it has been suggested that this might be related to the low number of introns in this yeast. Here we have tested this hypothesis by looking at other fungal genomes. Significantly, we find that fungal genomes with high numbers of introns also possess an SMN orthologue or at least its paralogue, SPF30.     

Adaptive protein evolution and regulatory divergence in Drosophila

Two recent studies demonstrated a positive correlation between divergence in gene expression and protein sequence in Drosophila. This correlation could be driven by positive selection or variation in functional constraint. To distinguish between these alternatives, we compared patterns of molecular evolution for 1,862 genes with two previously reported estimates of expression divergence in Drosophila. We found a slight negative trend (nonsignificant) between positive selection on protein sequence and divergence in expression levels between Drosophila melanogaster and Drosophila simulans. Conversely, shifts in expression patterns during Drosophila development showed a positive association with adaptive protein evolution, though as before the relationship was weak and not significant. Overall, we found no strong evidence for an increase in the incidence of positive selection on protein-coding regions in genes with divergent expression in Drosophila, suggesting that the previously reported positive association between protein and regulatory divergence primarily reflects variation in functional constraint.     

Inferring natural selection on fine-scale chromatin organization in yeast

Despite its potential role in the evolution of complex phenotypes, the detection of negative (purifying) and positive selection on noncoding regulatory sequence has been elusive because of the inherent difficulty in predicting the functional consequences of mutations on noncoding sequence. Because the functioning of regulatory sequence depends upon both chromatin configuration and cis-regulatory factor binding, we investigate the idea that the functional conservation of regulatory regions should be associated with the conservation of sequence-dependent bending properties of DNA that determine its affinity for the nucleosome. Recent advances in the computational prediction of sequence-dependent affinity to nucleosomes provide an opportunity to distinguish between neutral and nonneutral evolution of fine-scale chromatin organization. Here, a statistical test is presented for detecting evolutionary conservation and/or adaptive evolution of nucleosome affinity from interspecies comparisons of DNA sequences. Local nucleosome affinities of homologous sequences were calculated using 2 recently published methods. A randomization test was applied to sites of mutation to evaluate the similarity of DNA-nucleosome affinity between several closely related species of Saccharomyces yeast. For most of the genes we analyzed, the conservation of local nucleosome affinity was detected at a few distinct locations in the upstream noncoding region. Our results also demonstrate that different patterns of chromatin evolution have shaped DNA-nucleosome interaction at the core promoters of TATA-containing and TATA-less genes and that elevated purifying selection has maintained low affinity for nucleosome in the core promoters of the latter group. Across the entire yeast genome, DNA-nucleosome interaction was also discovered to be significantly more conserved in TATA-less genes compared with TATA-containing genes.     

Host-symbiont conflicts: positive selection on an outer membrane protein of parasitic but not mutualistic Rickettsiaceae

The Rickettsiaceae is a family of intracellular bacterial symbionts that includes both vertically transmitted parasites that spread by manipulating the reproduction of their host (Wolbachia in arthropods) and horizontally transmitted parasites (represented by Cowdria ruminantium), and mutualists (Wolbachia pipientis in nematode worms). We have investigated the nature of natural selection acting on an outer membrane protein, the wsp gene in Wolbachia and its homologue map1 in Cowdria, thought likely to be involved in host-parasite interactions in these bacteria. The ratio of nonsynonymous to synonymous substitution rates (d(N)/d(S)) at individual amino acid sites or at lineages within the gene's phylogeny was estimated using maximum likelihood models of codon substitution. The first hypothesis we tested was that this protein is under positive selection in the parasitic but not in the mutualistic Rickettsiaceae. This hypothesis was supported as positive selection and was detected in Cowdria and arthropod Wolbachia sequence evolution but not in the evolution of Wolbachia sequences from nematodes. Furthermore, this selection was concentrated outside the transmembrane region of the protein and, therefore, in the regions of the protein that may interact with the host. The second hypothesis tested was that positive selection would be stronger in the strains of arthropod Wolbachia that distort the host sex ratio than in those that induce cytoplasmic incompatibility. However, we found no support for this hypothesis. In conclusion, our results are consistent with the hypothesis that antagonistic coevolution causes faster evolution of surface protein sequences in parasites than in mutualists. Confirmation of this conclusion awaits the replication of these results both in additional genes and across more bacterial taxa. The regions of the wsp and map1 genes we identified as likely to be involved in host-parasite arms races should be examined in future studies of parasite virulence and host immune responses, and during the design of vaccines.     

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.     

Epistasis Creates Invariant Sites and Modulates the Rate of Molecular Evolution

Invariant sites are a common feature of amino acid sequence evolution. The presence of invariant sites is frequently attributed to the need to preserve function through site-specific conservation of amino acid residues. Amino acid substitution models without a provision for invariant sites often fit the data significantly worse than those that allow for an excess of invariant sites beyond those predicted by models that only incorporate rate variation among sites (e.g., a Gamma distribution). An alternative is epistasis between sites to preserve residue interactions that can create invariant sites. Through computer-simulated sequence evolution, we evaluated the relative effects of site-specific preferences and site-site couplings in the generation of invariant sites and the modulation of the rate of molecular evolution. In an analysis of ten major families of protein domains with diverse sequence and functional properties, we find that the negative selection imposed by epistasis creates many more invariant sites than site-specific residue preferences alone. Further, epistasis plays an increasingly larger role in creating invariant sites over longer evolutionary periods. Epistasis also dictates rates of domain evolution over time by exerting significant additional purifying selection to preserve site couplings. These patterns illuminate the mechanistic role of epistasis in the processes underlying observed site invariance and evolutionary rates.     

Sequence Fingerprints of MicroRNA Conservation

It is known that the conservation of protein-coding genes is associated with their sequences both various species, such as animals and plants. However, the association between microRNA (miRNA) conservation and their sequences in various species remains unexplored. Here we report the association of miRNA conservation with its sequence features, such as base content and cleavage sites, suggesting that miRNA sequences contain the fingerprints for miRNA conservation. More interestingly, different species show different and even opposite patterns between miRNA conservation and sequence features. For example, mammalian miRNAs show a positive/negative correlation between conservation and AU/GC content, whereas plant miRNAs show a negative/positive correlation between conservation and AU/GC content. Further analysis puts forward the hypothesis that the introns of protein-coding genes may be a main driving force for the origin and evolution of mammalian miRNAs. At the 5′ end, conserved miRNAs have a preference for base U, while less-conserved miRNAs have a preference for a non-U base in mammals. This difference does not exist in insects and plants, in which both conserved miRNAs and less-conserved miRNAs have a preference for base U at the 5′ end. We further revealed that the non-U preference at the 5′ end of less-conserved mammalian miRNAs is associated with miRNA function diversity, which may have evolved from the pressure of a highly sophisticated environmental stimulus the mammals encountered during evolution. These results indicated that miRNA sequences contain the fingerprints for conservation, and these fingerprints vary according to species. More importantly, the results suggest that although species share common mechanisms by which miRNAs originate and evolve, mammals may develop a novel mechanism for miRNA origin and evolution. In addition, the fingerprint found in this study can be predictor of miRNA conservation, and the findings are helpful in achieving a clearer understanding of miRNA function and evolution.     

Do essential genes evolve slowly?

Approximately two thirds of all knockouts of individual mouse genes give rise to viable fertile mice. These genes have thus been termed 'non-essential' in contrast to 'essential' genes, the knockouts of which result in death or infertility. Although non-essential genes are likely to be under selection that favours sequence conservation [1], it is predicted that they are less subject to such stabilising selection than essential genes, and hence evolve faster [2]. We have addressed this issue by analysing the molecular evolution of 108 non-essential and 67 essential genes that have been sequenced in both mouse and rat. On preliminary analysis, the non-essential genes appeared to be faster evolving than the essential ones. We found, however, that the non-essential class contains a disproportionate number of immune-system genes that may be under directional selection (that is, selection favouring change) because of host-parasite coevolution. After correction for this bias, we found that the rate at which genes evolve does not correlate with the severity of the knockout phenotype. This was corroborated by the finding that, whereas neuron-specific genes have significantly lower rates of change than other genes, essential and non-essential neuronal genes have comparable rates of evolution. Our findings most probably reflect strong selection acting against even very subtle deleterious phenotypes, and indicate that the putative involvement of directional selection in host-parasite coevolution and gene expression within the nervous system explains much more of the variance in rates of gene evolution than does the knockout phenotype.     

Eight new mtDNA sequences of glass sponges reveal an extensive usage of + 1 frameshifting in mitochondrial translation

Three previously studied mitochondrial genomes of glass sponges (phylum Porifera, class Hexactinellida) contained single nucleotide insertions in protein coding genes inferred as sites of + 1 translational frameshifting. To investigate the distribution and evolution of these sites and to help elucidate the mechanism of frameshifting, we determined eight new complete or nearly complete mtDNA sequences from glass sponges and examined individual mitochondrial genes from three others. We found nine new instances of single nucleotide insertions in these sequences and analyzed them both comparatively and phylogenetically. The base insertions appear to have been gained and lost repeatedly in hexactinellid mt protein genes, suggesting no functional significance for the frameshifting sites. A high degree of sequence conservation, the presence of unusual tRNAs, and a distinct pattern of codon usage suggest the “out-of-frame pairing” model of translational frameshifting. Additionally, we provide evidence that relaxed selection pressure on glass sponge mtDNA – possibly a result of their low growth rates and deep-water lifestyle – has allowed frameshift insertions to be tolerated for hundreds of millions of years. Our study provides the first example of a phylogenetically diverse and extensive usage of translational frameshifting in animal mitochondrial coding sequences.     

Evolutionary replacement of components in a salamander pheromone signaling complex: more evidence for phenotypic-molecular decoupling

In this article we explore the evolutionary history of a functional complex at the molecular level in plethodontid salamanders. The complex consists of a proteinaceous courtship pheromone, a pheromone-producing gland on the male's chin, and a set of behaviors for delivering the pheromone to the female. Long-term evolutionary stasis is the defining feature of this complex at both the morphological and behavioral levels. However, our previous assessment of the pheromone gene, plethodontid receptivity factor (PRF), revealed rapid evolution at the molecular level despite stasis at higher levels of organization. Analysis of a second pheromone gene, sodefrin precursor-like factor (SPF), now indicates that evolutionary decoupling in this complex is pervasive. The evolutionary profiles of SPF and PRF are remarkably similar in that: (a) both genes exhibit high levels of sequence diversity both within and across taxa, (b) genetic diversity has been driven by strong positive selection, and (c) the genes have evolved heterogeneously in different salamander lineages. The composition of the pheromone signal as a whole, however, has experienced an extraordinary evolutionary transition. Whereas SPF has been retained throughout the 100 MY radiation of salamanders, PRF has only recently been recruited to a pheromone function (27 million years ago). When SPF and PRF coexist in the same clade, they show contrasting patterns of evolution. When one shows rapid evolution driven by positive selection, the other shows neutral divergence restrained by purifying selection. In one clade, the origin and subsequent rapid evolution of PRF appear to have interfered with the evolution and persistence of SPF, leading to a pattern of evolutionary replacement. Overall, these two pheromone genes provide a revealing window on the dynamics that drive the evolution of multiple traits in a signaling complex.     

Human testis-specific genes are under relaxed negative selection

Recent studies have suggested that selective forces and constraints acting on genes varied during human evolution depending on the organ in which they are expressed. To gain insight into the evolution of organ determined negative selection forces, we compared the non-synonymous SNP diversity of genes expressed in different organs. Based on a HAPMAP dataset, we determined for each SNP its frequency in 11 human populations and, in each case, predicted whether or not the change it produces is deleterious. We have shown that, for all organs under study, SNPs predicted to be deleterious are present at a significantly lower frequency than SNPs predicted to be tolerated. However, testis-specific genes contain a higher proportion of deleterious SNPs than other organs. This study shows that negative selection is acting on the whole human genome, but that the action of negative selection is relaxed on testis-specific genes. This result adds to and expands the hypothesis of a recent evolutionary change in the human male reproductive system and its behavior.     

Rapid Evolution of Plethodontid Modulating Factor, a Hypervariable Salamander Courtship Pheromone, is Driven by Positive Selection

Sexual communication in plethodontid salamanders is mediated by a proteinaceous pheromone that a male delivers to a female during courtship, boosting her receptivity. The pheromone consists of three proteins from three unrelated protein families. These proteins are among a small group of pheromones known to affect female receptivity in vertebrates. Previously, we showed that the genes of two of these proteins (PRF and SPF) are prone to incessant evolution driven by positive selection, presumably as a consequence of coevolution with female receptors. In this report, we focus on the evolution of the third pheromone protein gene family, plethodontid modulating factor (PMF), to determine whether it shows the same pattern of diversification. We used RT-PCR in mental gland cDNA to survey PMF sequences from three genera of plethodontid salamanders (27 spp.) to measure rates of evolution, level of gene diversification, modes of selection, and types of amino acid substitution. Like PRF and SPF, PMF is produced by a multigene family characterized by gene duplication and high levels of polymorphism. PMF evolution is rapid, incessant, and driven by positive selection. PMF is more extreme in these dimensions than both PRF and SPF. Nestled within this extraordinary variation, however, is a signature of purifying selection, acting to preserve important structural and biochemical features of the PMF protein (i.e., secretion signal, cysteine residues, and pI). Although a pattern of persistent diversification exists at the molecular level, the morphological and behavioral aspects of the pheromone delivery system show evolutionary stasis over millions of years.     

High intron sequence conservation across three mammalian orders suggests functional constraints

Several studies have demonstrated high levels of sequence conservation in noncoding DNA compared between two species (e.g., human and mouse), and interpreted this conservation as evidence for functional constraints. If this interpretation is correct, it suggests the existence of a hidden class of abundant regulatory elements. However, much of the noncoding sequence conserved between two species may result from chance or from small-scale heterogeneity in mutation rates. Stronger inferences are expected from sequence comparisons using more than two taxa, and by testing for spatial patterns of conservation in addition to primary sequence similarity. We used a Bayesian local alignment method to compare approximately 10 kb of intron sequence from nine genes in a pairwise manner between human, whale, and seal to test whether the degree and pattern of conservation is consistent with neutral divergence. Comparison of the three sets of conserved gapless pairwise blocks revealed the following patterns: The proportion of identical intron nucleotides averaged 47% in pairwise comparisons and 28% across the three taxa. Proportions of conserved sequence were similar in unique sequence and general mammalian repetitive elements. We simulated sequence evolution under a neutral model using published estimates of substitution rate heterogeneity for noncoding DNA and found pairwise identity at 33% and three-taxon identity at 16% of nucleotide sites. Spatial patterns of primary sequence conservation were also nonrandomly distributed within introns. Overall, segments of intron sequence closer to flanking exons were significantly more conserved than interior intron sequence. This level of intron sequence conservation is above that expected by chance and strongly suggests that intron sequences are playing a larger functional role in gene regulation than previously realized.     

蕨类植物叶绿体rps12基因的分子进化研究

以68种蕨类植物和2种石松类植物的rps12基因为对象,在系统发育背景下,结合最大似然法,使用HyPhy和PAML软件对该基因进行进化速率和适应性进化研究。结果显示:位于IR区的外显子2~3,其替换率明显降低,rps12基因编码序列的替换率也随之降低,且rps12基因密码子第3位的GC含量明显升高;在蕨类植物的进化过程中,3′-rps12更倾向定位于IR区,以保持较低的替换率;rps12基因编码的123个氨基酸位点中,共检测到4个正选择位点和116个负选择位点。研究结果表明基因序列进入到IR区后,显示出降低的替换率;强烈的负选择压力表明RPS12蛋白的高度保守性以及rps12基因的功能和结构已经趋于稳定。     

Patterns of molecular evolution and predicted function in thaumatin-like proteins of Populus trichocarpa

Some pathogenesis-related proteins (PR proteins) are subject to positive selection, while others are under negative selection. Here, we report the patterns of molecular evolution in thaumatin-like protein (TLP, PR5 protein) genes of Populus trichocarpa. Signs of positive selection were found in 20 out of 55 Populus TLPs using the likelihood ratio test and ML-based Bayesian methods. Due to the connection between the acidic cleft and the antifungal activity, the secondary structure and three-dimensional structure analyses predicted antifungal activity β-1,3-glucanase activities in these TLPs. Moreover, the coincidence with variable basic sites in the acidic cleft and positively selected sites suggested that fungal diseases may have been the main environmental stress that drove rapid adaptive evolution in Populus.     

Regulation of neural determination by evolutionarily conserved signals: anti-BMP factors and what next?

The evolutionary conservation of Chordin/bone morphogenetic protein (BMP) signaling supports the hypothesis of dorsal-ventral axis inversion of vertebrates and invertebrates, and implies that the invention of a central nervous system occurred only once during animal evolution. This hypothesis is further strengthened by recent findings of the conservation of downstream genes and modifier genes of neural induction. On the other hand, in contrast with such gross conservation, recent data suggest that the requirement for some signals in neural determination may differ even within the vertebrate subphylum.     

Evolutionary Constraint and Disease Associations of Post-Translational Modification Sites in Human Genomes

Interpreting the impact of human genome variation on phenotype is challenging. The functional effect of protein-coding variants is often predicted using sequence conservation and population frequency data, however other factors are likely relevant. We hypothesized that variants in protein post-translational modification (PTM) sites contribute to phenotype variation and disease. We analyzed fraction of rare variants and non-synonymous to synonymous variant ratio (Ka/Ks) in 7,500 human genomes and found a significant negative selection signal in PTM regions independent of six factors, including conservation, codon usage, and GC-content, that is widely distributed across tissue-specific genes and function classes. PTM regions are also enriched in known disease mutations, suggesting that PTM variation is more likely deleterious. PTM constraint also affects flanking sequence around modified residues and increases around clustered sites, indicating presence of functionally important short linear motifs. Using target site motifs of 124 kinases, we predict that at least ∼180,000 motif-breaker amino acid residues that disrupt PTM sites when substituted, and highlight kinase motifs that show specific negative selection and enrichment of disease mutations. We provide this dataset with corresponding hypothesized mechanisms as a community resource. As an example of our integrative approach, we propose that PTPN11 variants in Noonan syndrome aberrantly activate the protein by disrupting an uncharacterized cluster of phosphorylation sites. Further, as PTMs are molecular switches that are modulated by drugs, we study mutated binding sites of PTM enzymes in disease genes and define a drug-disease network containing 413 novel predicted disease-gene links.