Evidence Of Positive Darwinian Selection Research Articles

Improved inference of site-specific positive selection under a generalized parametric codon model when there are multinucleotide mutations and multiple nonsynonymous rates

BackgroundAn excess of nonsynonymous substitutions, over neutrality, is considered evidence of positive Darwinian selection. Inference for proteins often relies on estimation of the nonsynonymous to synonymous ratio (ω = dN/dS) within a codon model. However, to ease computational difficulties, ω is typically estimated assuming an idealized substitution process where (i) all nonsynonymous substitutions have the same rate (regardless of impact on organism fitness) and (ii) instantaneous double and triple (DT) nucleotide mutations have zero probability (despite evidence that they can occur). It follows that estimates of ω represent an imperfect summary of the intensity of selection, and that tests based on the ω > 1 threshold could be negatively impacted.ResultsWe developed a general-purpose parametric (GPP) modelling framework for codons. This novel approach allows specification of all possible instantaneous codon substitutions, including multiple nonsynonymous rates (MNRs) and instantaneous DT nucleotide changes. Existing codon models are specified as special cases of the GPP model. We use GPP models to implement likelihood ratio tests for ω > 1 that accommodate MNRs and DT mutations. Through both simulation and real data analysis, we find that failure to model MNRs and DT mutations reduces power in some cases and inflates false positives in others. False positives under traditional M2a and M8 models were very sensitive to DT changes. This was exacerbated by the choice of frequency parameterization (GY vs. MG), with rates sometimes > 90% under MG. By including MNRs and DT mutations, accuracy and power was greatly improved under the GPP framework. However, we also find that over-parameterized models can perform less well, and this can contribute to degraded performance of LRTs.ConclusionsWe suggest GPP models should be used alongside traditional codon models. Further, all codon models should be deployed within an experimental design that includes (i) assessing robustness to model assumptions, and (ii) investigation of non-standard behaviour of MLEs. As the goal of every analysis is to avoid false conclusions, more work is needed on model selection methods that consider both the increase in fit engendered by a model parameter and the degree to which that parameter is affected by un-modelled evolutionary processes.

Positive selection adaptation of two-domain arginine kinase (AK) from cold seep Vesicomyidae clams.

Arginine kinase (AK) is an important member of Phosphagen kinases which engage in energy metabolism process, and AKs from cold seep clams may develop an effective mechanism to adapt a special habitat (e.g. low temperature). Three Vesicomyidae clams and seven Veneridae clams (belong to the same Order Veneroida) were chosen to analyze the evolution of two-domain AKs. In the present study, ten two-domain AKs were identified and Neighbor-joining tree showed that AKs were divided into two groups. Branch-site model indicated that two-domain AKs were subjected to strong positive selection (ω2a = 17.5058). 16 positively selective sites were detected and five of them showed posterior probabilities of 0.95 or more. Comparative analysis found that domain 2 might be suffered from more evolutionary selection pressure than domain 1, as most positively sites were located at domain 2. Residue Pro (positively selective site) (587P in ApAK) in domain 2 from all Vesicomyidae AKs might participate in change of the synergism and in the function of its cold-adapted characteristics. In conclusion, our studies provide evidence of positive Darwinian selection in the two-domain AKs family of Vesicomyidae clams, and may contribute to a better understanding of its adaptation mechanisms to cold seep habitats.

Evolution of methicillin-resistant Staphylococcus aureus: Evidence of positive selection in a penicillin-binding protein (PBP) 2a coding gene mecA

Methicillin-resistant Staphylococcus aureus (S. aureus) (MRSA) represents more and more S. aureus infections. MecA, the novel coding gene of penicillin-binding protein (PBP) 2a of MRSA, is the key resistance factor of β-lactam, but little is known about the evolution of this gene. Given the crucial role of mecA in S. aureus physiology and β-lactam resistance, the selective forces may contribute to adaptation of the bacteria to the special environments such as its host or antibiotics. To understand the evolution of this gene, we screened GenBank database and analyzed mecA of 249 S. aureus strains. Twenty-nine unique alleles with 26 unique amino acid sequences were identified. Phylogenetic analysis showed three main groups of mecA in the global S. aureus strains. Analysis of these alleles using codon-substitution models (M8, M3, and M2a) and likelihood ratio tests (LRTs) of the codeML package and a random-effects likelihood (REL) method of HyPhy package for the site-specific ratio of nonsynonymous to synonymous substitution rates suggested that fourteen sites in the allosteric domain of PBP2a have been subjected to strong positive selection pressure. Mutations of two positive selection sites (N146K and E239K) were reported to be essential for ceftaroline- or L-695, 256-resistant. Further study indicated that the positive selection pressure might be more likely related to the host's inflammatory or immune response during S. aureus infection. Our studies provide the first evidence of positive Darwinian selection in the mecA of S. aureus, contributing to a better understanding of the adaptive mechanism of this bacterium.

Extremely high genetic diversity in a single tumor points to prevalence of non-Darwinian cell evolution

The prevailing view that the evolution of cells in a tumor is driven by Darwinian selection has never been rigorously tested. Because selection greatly affects the level of intratumor genetic diversity, it is important to assess whether intratumor evolution follows the Darwinian or the non-Darwinian mode of evolution. To provide the statistical power, many regions in a single tumor need to be sampled and analyzed much more extensively than has been attempted in previous intratumor studies. Here, from a hepatocellular carcinoma (HCC) tumor, we evaluated multiregional samples from the tumor, using either whole-exome sequencing (WES) (n = 23 samples) or genotyping (n = 286) under both the infinite-site and infinite-allele models of population genetics. In addition to the many single-nucleotide variations (SNVs) present in all samples, there were 35 "polymorphic" SNVs among samples. High genetic diversity was evident as the 23 WES samples defined 20 unique cell clones. With all 286 samples genotyped, clonal diversity agreed well with the non-Darwinian model with no evidence of positive Darwinian selection. Under the non-Darwinian model, MALL (the number of coding region mutations in the entire tumor) was estimated to be greater than 100 million in this tumor. DNA sequences reveal local diversities in small patches of cells and validate the estimation. In contrast, the genetic diversity under a Darwinian model would generally be orders of magnitude smaller. Because the level of genetic diversity will have implications on therapeutic resistance, non-Darwinian evolution should be heeded in cancer treatments even for microscopic tumors.

Genome-wide detection of selection and other evolutionary forces.

As is well known, pathogenic microbes evolve rapidly to escape from the host immune system and antibiotics. Genetic variations among microbial populations occur frequently during the long-term pathogen-host evolutionary arms race, and individual mutation beneficial for the fitness can be fixed preferentially. Many recent comparative genomics studies have pointed out the importance of selective forces in the molecular evolution of bacterial pathogens. The public availability of large-scale next-generation sequencing data and many state-of-the-art statistical methods of molecular evolution enable us to scan genome-wide alignments for evidence of positive Darwinian selection, recombination, and other evolutionary forces operating on the coding regions. In this chapter, we describe an integrative analysis pipeline and its application to tracking featured evolutionary trajectories on the genome of an animal pathogen. The evolutionary analysis of the protein-coding part of the genomes will provide a wide spectrum of genetic variations that play potential roles in adaptive evolution of bacteria.

Phylogenomics of strongylocentrotid sea urchins

BackgroundStrongylocentrotid sea urchins have a long tradition as model organisms for studying many fundamental processes in biology including fertilization, embryology, development and genome regulation but the phylogenetic relationships of the group remain largely unresolved. Although the differing isolating mechanisms of vicariance and rapidly evolving gamete recognition proteins have been proposed, a stable and robust phylogeny is unavailable.ResultsWe used a phylogenomic approach with mitochondrial and nuclear genes taking advantage of the whole-genome sequencing of nine species in the group to establish a stable (i.e. concordance in tree topology among multiple lies of evidence) and robust (i.e. high nodal support) phylogenetic hypothesis for the family Strongylocentrotidae. We generated eight draft mitochondrial genome assemblies and obtained 13 complete mitochondrial genes for each species. Consistent with previous studies, mitochondrial sequences failed to provide a reliable phylogeny. In contrast, we obtained a very well-supported phylogeny from 2301 nuclear genes without evidence of positive Darwinian selection both from the majority of most-likely gene trees and the concatenated fourfold degenerate sites: ((P. depressus, (M. nudus, M. franciscanus), (H. pulcherrimus, (S. purpuratus, (S. fragilis, (S. pallidus, (S. droebachiensis, S. intermedius)). This phylogeny was consistent with a single invasion of deep-water environments followed by a holarctic expansion by Strongylocentrotus. Divergence times for each species estimated with reference to the divergence times between the two major clades of the group suggest a correspondence in the timing with the opening of the Bering Strait and the invasion of the holarctic regions.ConclusionsNuclear genome data contains phylogenetic signal informative for understanding the evolutionary history of this group. However, mitochondrial genome data does not. Vicariance can explain major patterns observed in the phylogeny. Other isolating mechanisms are appropriate to explore in this system to help explain divergence patterns not well supported by vicariance, such as the effects of rapidly evolving gamete recognition proteins on isolating populations. Our findings of a stable and robust phylogeny, with the increase in mitochondrial and nuclear comparative genomic data, provide a system in which we can enhance our understanding of molecular evolution and adaptation in this group of sea urchins.

The first evidence of positive selection in peptidoglycan recognition protein (PGRP) genes of Crassostrea gigas

The oyster Crassostrea gigas is thought to have developed effective immunity to potentially harmful pathogens while under continuous exposure to marine microorganisms; however, the evolutionary mechanisms by which such immunity developed has not been understood. To understand the evolution of immunity, we characterized the family of peptidoglycan recognition proteins in the oyster (CgPGRPs). These proteins are crucial pattern recognition receptors for peptidoglycans (PGNs) and thereby, for activating the innate immune response of host. Herein, we identify seven new CgPGRP genes. Phylogenetic analysis of the seven new and five previously reported CgPRGP genes reveals that the CgPRGP gene family can be clustered into two groups, CgPRGPS and CgPRGPL. Moreover, the CgPRGPS group can be further divided into five subgroups. A codon-substitution model and three likelihood ratio tests (LRTs) suggest that seven sites in the CgPGRP family of genes have been subjected to strong positive selection (ω = 3.035–4.143), Three dimensional modeling revealed that these sites are found primarily at the periphery of coils and α-helices rather than in β-strands, perhaps allowing PGRP to adapt to, and recognize, variability of PGN structure. In conclusion, our studies provide the first evidence of positive Darwinian selection in the CgPGRP family, contributing to a better understanding of the adaptive mechanism of host-pathogens interaction in marine mollusks.

Near Neutrality

The nearly neutral theory represents a development of Kimura's neutral theory of molecular evolution that makes testable predictions that go beyond a mere null model. Recent evidence has strongly supported several of these predictions, including the prediction that slightly deleterious variants will accumulate in a species that has undergone a severe bottleneck or in cases where recombination is reduced or absent. Because bottlenecks often occur in speciation and slightly deleterious mutations in coding regions will usually be nonsynonymous, we should expect that the ratio of nonsynonymous to synonymous fixed differences between species should often exceed the ratio of nonsynonymous to synonymous polymorphisms within species. Many data support this prediction, although they have often been wrongly interpreted as evidence for positive Darwinian selection. The use of conceptually flawed tests for positive selection has become widespread in recent years, seriously harming the quest for an understanding of genome evolution. When properly analyzed, many (probably most) claimed cases of positive selection will turn out to involve the fixation of slightly deleterious mutations by genetic drift in bottlenecked populations. Slightly deleterious variants are a transient feature of evolution in the long term, but they have substantially affected contemporary species, including our own.

Evidence for positive Darwinian selection on the hepcidin gene of Perciform and Pleuronectiform fishes

Hepcidin is a small cysteine-rich peptide that plays an important role in antimicrobial activity and in maintaining iron homeostasis in vertebrates. Here we report on the underlying mechanism that maintains high sequence diversities among the hepcidin-like variants of perciform and pleuronectiform fishes. In contrast to mammals, maximum likelihood-based codon substitution analyses revealed that positive Darwinian selection (nonsynonymous to synonymous substitution, omega > 1) is the likely cause of accelerated rate of amino acid substitutions in the hepcidin mature peptide region of these fishes. Comparison of models incorporating positive selection (omega > 1) at certain sites with models not incorporating positive selection (omega < 1) failed to reject (p = 0) the evidence of positive selection among the codon sites of percifom and pleuronectiform hepcidin. The adaptive evolution of this peptide in perciform and pleuronectiform fishes might be directed by pathogens when the host is exposed to new habitats/environments.

Looking for Darwin in all the wrong places: the misguided quest for positive selection at the nucleotide sequence level

Recent years have seen an explosion of interest in evidence for positive Darwinian selection at the molecular level. This quest has been hampered by the use of statistical methods that fail adequately to rule out alternative hypotheses, particularly the relaxation of purifying selection and the effects of population bottlenecks, during which the effectiveness of purifying selection is reduced. A further problem has been the assumption that positive selection will generally involve repeated amino-acid changes to a single protein. This model was derived from the case of the vertebrate major histocompatibility complex (MHC), but the MHC proteins are unusual in being involved in protein-protein recognition and in a co-evolutionary process of pathogens. There is no reason to suppose that repeated amino-acid changes to a single protein are involved in selectively advantageous phenotypes in general. Rather adaptive phenotypes are much more likely to result from other causes, including single amino-acid changes; deletion or silencing of genes or changes in the pattern of gene expression.

Evidence for Positive Darwinian Selection of Vip Gene in Bacillus thuringiensis

Vegetative insecticidal proteins (VIPs), produced during the vegetative stage of their growth in Bacillus thuringiensis, are a group of insecticidal proteins and represent the second generation of insecticidal trans-genes that will complement the novel d-endotoxins in future. Fewer structural and functional relationships of Vip proteins are known in comparison with those of d-endotoxins. In this study, both the maximum-likelihood methods and the maximum parsimony based sliding window analysis were used to evaluate the molecular evolution of Vip proteins. As a result, strong evidence was found that Vip proteins are subject to the high rates of positive selection, and 16 sites are identified to be under positive selection using the Bayes Empirical Bayesian method. Interestingly, all these positively selected sites are located from site-705 to site-809 in the C-terminus of the Vip proteins. Most of these sites are exposed and clustered in the loop regions when mapped onto its computational predicted secondary tertiary and a part of the tertiary structure. It has been postulated that the high divergence in the C-terminal of Vip proteins may not result from the lack of functional constraints, but rather from the rapid mutation to adapt their targeted insects, driven by positive selection. The potential positive selection pressures may be an attempt to adapt for the “arm race” between Vip proteins and the targeted insects, or to enlarge their target's host range. Sites identified to be under positive selection may be related to the insect host range, which may shed a light on the investigation of the Vip proteins’ structural and functional relationships.

Detecting positive darwinian selection in brain-expressed genes during human evolution

To understand the genetic basis that underlies the phenotypic divergence between human and non-human primates, we screened a total of 7176 protein-coding genes expressed in the human brain and compared them with the chimpanzee orthologs to identify genes that show evidence of rapid evolution in the human lineage. Our results showed that the nonsynonymous/synonymous substitution (Ka/Ks) ratio for genes expressed in the brain of human and chimpanzee is 0.3854, suggesting that the brain-expressed genes are under functional constraint. The X-linked human brain-expressed genes evolved more rapidly than autosomal ones. We further dissected the molecular evolutionary patterns of 34 candidate genes by sequencing representative primate species to identify lineage-specific adaptive evolution. Fifteen out of the 34 candidate genes showed evidence of positive Darwinian selection in human and/or chimpanzee lineages. These genes are predicted to play diverse functional roles in embryonic development, spermatogenesis and male fertility, signal transduction, sensory nociception, and neural function. This study together with others demonstrated the usefulness and power of phylogenetic comparison of multiple closely related species in detecting lineage-specific adaptive evolution, and the identification of the positively selected brain-expressed genes may add new knowledge to the understanding of molecular mechanism of human origin.

Adaptive evolution after duplication of penaeidin antimicrobial peptides

Penaeidin antimicrobial peptides in penaeid shrimps are an important component of their innate immune system that provides immunity against infection caused by several gram-positive bacteria and filamentous fungal species. Despite the knowledge on the identification and characterization of these peptides in penaeid shrimps, little is known about the evolutionary pattern of these peptides and the underlying genetic mechanisms that maintain high sequence diversities in the penaeidin gene family. Based on the phylogenetic analyses and maximum likelihood-based codon substitution analyses, here we present the convincing evidence that multiple copies of penaeidins have evolved by gene duplication, and positive Darwinian selection (adaptive evolution) is the likely cause of accelerated rate of amino acid substitutions among these duplicated genes. While the average ratio of non-synonymous to synonymous substitutions (ω) for the entire coding region of both active domains is 0.9805, few codon sites showed significantly higher ω (3.73). The likelihood ratio tests that compare models incorporating positive selection (ω > 1) at certain codon sites with models not incorporating positive selection (ω < 1), failed to reject ( p = 0) the evidence of positive Darwinian selection. The rapid adaptive evolution of this gene family might be directed by the pathogens and the faster rate of amino acid substitutions in the N-terminal proline-rich and C-terminal cysteine-rich domains could be due to their direct involvement in the protection against pathogens. When the host expose to different habitats/environment an accelerated rate of amino acid substitutions in both the active domains may also be expected.

Genomic organization and evolution of the ULBP genes in cattle

BackgroundThe cattle UL16-binding protein 1 (ULBP1) and ULBP2 genes encode members of the MHC Class I superfamily that have homology to the human ULBP genes. Human ULBP1 and ULBP2 interact with the NKG2D receptor to activate effector cells in the immune system. The human cytomegalovirus UL16 protein is known to disrupt the ULBP-NKG2D interaction, thereby subverting natural killer cell-mediated responses. Previous Southern blotting experiments identified evidence of increased ULBP copy number within the genomes of ruminant artiodactyls. On the basis of these observations we hypothesized that the cattle ULBPs evolved by duplication and sequence divergence to produce a sufficient number and diversity of ULBP molecules to deliver an immune activation signal in the presence of immunogenic peptides. Given the importance of the ULBPs in antiviral immunity in other species, our goal was to determine the copy number and genomic organization of the ULBP genes in the cattle genome.ResultsSequencing of cattle bacterial artificial chromosome genomic inserts resulted in the identification of 30 cattle ULBP loci existing in two gene clusters. Evidence of extensive segmental duplication and approximately 14 Kbp of novel repetitive sequences were identified within the major cluster. Ten ULBPs are predicted to be expressed at the cell surface. Substitution analysis revealed 11 outwardly directed residues in the predicted extracellular domains that show evidence of positive Darwinian selection. These positively selected residues have only one residue that overlaps with those proposed to interact with NKG2D, thus suggesting the interaction with molecules other than NKG2D.ConclusionThe ULBP loci in the cattle genome apparently arose by gene duplication and subsequent sequence divergence. Substitution analysis of the ULBP proteins provided convincing evidence for positive selection on extracellular residues that may interact with peptide ligands. These results support our hypothesis that the cattle ULBPs evolved under adaptive diversifying selection to avoid interaction with a UL16-like molecule whilst preserving the NKG2D binding site. The large number of ULBPs in cattle, their extensive diversification, and the high prevalence of bovine herpesvirus infections make this gene family a compelling target for studies of antiviral immunity.

Sequence divergence at the putative flowering time locus COL1 in Brassicaceae

An insertion/deletion polymorphism ( Ind2) in the Brassica nigra CONSTANS LIKE 1 ( Bni COL1) gene was previously found to be associated with variation in flowering time. In the present study we examine the inter-specific divergence of COL1 in the family Brassicaceae. Analysis of codon substitution models did not reveal evidence of positive Darwinian selection, but comparisons of the COL1 gene in different species revealed a surprising number of indels. A total of 24 indels were found in the 650 bp of the middle variable region of the gene. This high number of indels could reflect a lack of constraint on length of this region of the protein, or the effect of positive selection. The number of indels was close to that expected in non-coding DNA, but the indels were longer in COL1 than those observed in non-coding regions. Reconstruction of indel evolution indicated that most indels resulted from deletions rather than insertions. The Ind2 indel that has shown association with flowering time in Brassica nigra exhibited a remarkable distribution in the Brassicaceae family, indicating that the polymorphism may have persisted more than ten million years. Considering presumed historic populations sizes of Brassicaceae species, such a long persistence time seems unlikely for a neutral polymorphism.

Evidence of Positive Darwinian Selection in Putative Meningococcal Vaccine Antigens

Meningococcal meningitidis is a life-threatening disease. In Europe and the United States the majority of cases are caused by virulent meningococcal strains belonging to serogroup B. Presently there is no effective vaccine against serogroup B strains, as traditional vaccine antigens such as polysaccharide capsules are unusable as they lead to autoimmunity. The year 2000 saw the publication of the complete genome of Neisseria meningitidis MC58, a virulent serogroup B bacterium. Working in conjunction with the sequencing project, researchers endeavored to locate highly conserved membrane-associated proteins that elicit an immune response. It is hoped that these proteins will provide a basis for novel vaccines against serogroup B strains. A number of potential vaccine antigens have been located and are presently in phase I clinical trials. Recently many reports pertaining to the evidence of positive Darwinian selection in membrane proteins of pathogens have been reported. This study utilized in silico methods to test for evidence of historical positive Darwinian selection in seven such vaccine candidates. We found that two of these proteins show signatures of adaptive evolution, while the remaining proteins show evidence of strong purifying selection. This has significant implications for the design of a vaccine against serogroup B strains, as it has been shown that vaccines that target epitopes that are under strong purifying selection are better than those that target variable epitopes.

Evidence of Positive Darwinian Selection in Omp85, a Highly ConservedBacterial Outer Membrane Protein Essential for Cell Viability

Omp85 is a highly conserved outer membrane protein found in all gram-negative bacteria. It is essential for bacterial cell viability and plays an integral function in the positioning and folding of other outer membrane proteins into the bacterial outer membrane. We have employed a maximum likelihood and a maximum parsimony approach to detect evidence of positive Darwinian selection in Omp85 homologues from 10 delta-proteobacteria and have identified 14 amino acid sites that show evidence of being under the influence of adaptive evolution. Interestingly all sites bar one are concentrated within surface loops of the protein that most likely interact with host immune response or the surrounding environment. Alternatively amino acids within membrane-spanning regions of the protein are found to be under purifying selection most likely as a result of structural constraints.

Mitochondrial DNA haplotype frequencies in natural and experimental populations of Drosophila subobscura.

The evolution of Drosophila subobscura mitochondrial DNA has been studied in experimental populations, founded with flies from a natural population from Esporles (Majorca, Balearic Islands, Spain). This population, like other European ones, is characterized by the presence of two very common (>96%) mitochondrial haplotypes (called I and II) and rare and endemic haplotypes that appear at very low frequencies. There is no statistical evidence of positive Darwinian selection acting on the mitochondrial DNA variants according to Tajima's neutrality test. Two experimental populations, with one replicate each, were established with flies having a heterogeneous nuclear genetic background, which was representative of the composition of the natural population. Both populations were started with the two most frequent mitochondrial haplotypes, but at different initial frequencies. After 13 to 16 generations, haplotype II reached fixation in three cages and its frequency was 0.89 by generation 25 in the fourth cage. Random drift can be rejected as the force responsible for the observed changes in haplotype frequencies. There is not only statistical evidence of a linear trend favoring a mtDNA (haploid) fitness effect, but also of a significant nonlinear deviation that could be due to a nuclear component.

Sequence evolution of the porB gene of Neisseria gonorrhoeae and Neisseria meningitidis: evidence of positive Darwinian selection.

Protein 1 (PI) is a major porin of Neisseria gonorrhoeae and Neisseria meningitidis and is encoded by a single locus, porB. Alleles of the porB locus of N. gonorrhoeae are assigned to two homology groups, PI(A) and PI(B), on the basis of immunological and structural similarity. In a like manner, alleles of the porB locus of the closely related bacterium, N. meningitidis, are allocated into class 2 and class 3 homology groups. An individual strain of N. gonorrhoeae or N. meningitidis expresses either one or other of these porin homology groups but never both, and the antigenic reactions of these highly diverse outer membrane proteins form part of the N. gonorrhoeae and N. meningitidis serotyping schemes. A comparison of the number of synonymous and nonsynonymous substitutions per site between the two most divergent alleles of each of these four groups of porB alleles shows that PI(A) alleles have accumulated significantly more nonsynonymous substitutions per site than synonymous substitutions. In contrast the distribution of synonymous and nonsynonymous substitutions between alleles of class 2 and class 3 porins are not significantly different from random. We localize the regions of the PI(A) alleles with an excess of amino acid changes to the surface-exposed loops of these outer membrane proteins and suggest that positive Darwinian selection for diversity, driven by the human immune system, can most easily explain the allelic polymorphism and the pattern of synonymous and nonsynonymous substitutions.