Hitchhiking under positive Darwinian selection

Positive selection can be inferred from its effect on linked neutral variation. In the restrictive case when there is no recombination, all linked variation is removed. If recombination is present but rare, both deterministic and stochastic models of positive selection show that linked variation hitchhikes to either low or high frequencies. While the frequency distribution of variation can be influenced by a number of evolutionary processes, an excess of derived variants at high frequency is a unique pattern produced by hitchhiking (derived refers to the nonancestral state as determined from an outgroup). We adopt a statistic, H, to measure an excess of high compared to intermediate frequency variants. Only a few high-frequency variants are needed to detect hitchhiking since not many are expected under neutrality. This is of particular utility in regions of low recombination where there is not much variation and in regions of normal or high recombination, where the hitchhiking effect can be limited to a small (<1 kb) region. Application of the H test to published surveys of Drosophila variation reveals an excess of high frequency variants that are likely to have been influenced by positive selection.     

Evolution of snake venom disintegrins by positive Darwinian selection

PII-disintegrins, cysteine-rich polypeptides broadly distributed in the venoms of geographically diverse species of vipers and rattlesnakes, antagonize the adhesive functions of beta(1) and beta(3) integrin receptors. PII-disintegrins evolved in Viperidae by neofunctionalization of disintegrin-like domains of duplicated PIII-snake venom hemorrhagic metalloproteinase (SVMP) genes recruited into the venom proteome before the radiation of the advanced snakes. Minimization of the gene (loss of introns and coding regions) and the protein structures (successive loss of disulfide bonds) underpins the postduplication divergence of disintegrins. However, little is known about the underlying genetic mechanisms that have generated the structural and functional diversity among disintegrins. Phylogenetic inference and maximum likelihood-based codon substitution approaches were used to analyze the evolution of the disintegrin family. The topology of the phylogenetic tree does not parallel that of the species tree. This incongruence is consistent with that expected for a multigene family undergoing a birth-and-death process in which the appearance and disappearance of loci are being driven by selection. Cysteine and buried residues appear to be under strong purifying selection due to their role in maintaining the active conformation of disintegrins. Divergence of disintegrins is strongly influenced by positive Darwinian selection causing accelerated rate of substitution in a substantial proportion of surface-exposed disintegrin residues. Global and lineage-specific sites evolving under diversifying selection were identified. Several sites are located within the integrin-binding loop and the C-terminal tail, two regions that form a conformational functional epitope. Arginine-glycine-aspartic acid (RGD) was inferred to represent the ancestral integrin-recognition motif, which emerged from the subgroup of PIII-SVMPs bearing the RDECD sequence. The most parsimonious nucleotide substitution model required for the emergence of all known disintegrin's integrin inhibitory motifs from an ancestral RGD sequence involves a minimum of three mutations. The adaptive advantage of the emergence of motifs targeting beta(1) integrins and the role of positively selected sites located within nonfunctional disintegrin regions appear to be difficult to rationalize in the context of a predator-prey arms race. Perhaps, this represents a consequence of the neofunctionalization potential of the disintegrin domain, a feature that may underlie its recruitment into the venom proteome followed by its successful transformation into a toxin.     

Evolution of adaptive phenotypic traits without positive Darwinian selection

Recent evidence suggests the frequent occurrence of a simple non-Darwinian (but non-Lamarckian) model for the evolution of adaptive phenotypic traits, here entitled the plasticity-relaxation-mutation (PRM) mechanism. This mechanism involves ancestral phenotypic plasticity followed by specialization in one alternative environment and thus the permanent expression of one alternative phenotype. Once this specialization occurs, purifying selection on the molecular basis of other phenotypes is relaxed. Finally, mutations that permanently eliminate the pathways leading to alternative phenotypes can be fixed by genetic drift. Although the generality of the PRM mechanism is at present unknown, I discuss evidence for its widespread occurrence, including the prevalence of exaptations in evolution, evidence that phenotypic plasticity has preceded adaptation in a number of taxa and evidence that adaptive traits have resulted from loss of alternative developmental pathways. The PRM mechanism can easily explain cases of explosive adaptive radiation, as well as recently reported cases of apparent adaptive evolution over ecological time.     

Molecular evolution and positive Darwinian selection of the chloroplast maturase matK

It is not clear whether matK evolves under Darwinian selection. In this study, 70 plant groups, representing 2,279 species at various evolutionary levels, were used to illustrate the molecular adaptation and evolutionary dynamics of gene divergence in matKs. Selective influences were investigated using standard dN/dS ratio methods. Analyses revealed the presence of positive selection in matKs of 32 plant groups. More positively selected sites were detected in the N-terminal region than in the RT domain and domain X of matK. Moreover, removing amino acid sites that are under positive selection has a significant effect on the bootstrap values of phylogenetic reconstruction. Our results suggest that the rapidly evolving matK evolves under positive selection in some lineages of land plants. Several regions of matK have experienced molecular adaptation, which fine-tunes maturase performance.     

The influence of phylogenetic uncertainty on the detection of positive Darwinian selection

The power of maximum likelihood tests of positive selection on protein-coding genes depends heavily on detecting and accounting for potential biases in the studied data set. Although the influence of transition:transversion and codon biases have been investigated in detail, little is known about how inaccuracy in the phylogeny used during the calculations affects the performance of these tests. In this study, 3 empirical data sets are analyzed using sets of simulated topologies corresponding to low, intermediate, and high levels of phylogenetic uncertainty. The detection of positive selection was largely unaffected by errors in the underlying phylogeny. However, the number of sites identified as being under positive selection tended to be overestimated.     

Mitochondrial targeting adaptation of the hominoid-specific glutamate dehydrogenase driven by positive Darwinian selection

Many new gene copies emerged by gene duplication in hominoids, but little is known with respect to their functional evolution. Glutamate dehydrogenase (GLUD) is an enzyme central to the glutamate and energy metabolism of the cell. In addition to the single, GLUD-encoding gene present in all mammals (GLUD1), humans and apes acquired a second GLUD gene (GLUD2) through retroduplication of GLUD1, which codes for an enzyme with unique, potentially brain-adapted properties. Here we show that whereas the GLUD1 parental protein localizes to mitochondria and the cytoplasm, GLUD2 is specifically targeted to mitochondria. Using evolutionary analysis and resurrected ancestral protein variants, we demonstrate that the enhanced mitochondrial targeting specificity of GLUD2 is due to a single positively selected glutamic acid-to-lysine substitution, which was fixed in the N-terminal mitochondrial targeting sequence (MTS) of GLUD2 soon after the duplication event in the hominoid ancestor approximately 18-25 million years ago. This MTS substitution arose in parallel with two crucial adaptive amino acid changes in the enzyme and likely contributed to the functional adaptation of GLUD2 to the glutamate metabolism of the hominoid brain and other tissues. We suggest that rapid, selectively driven subcellular adaptation, as exemplified by GLUD2, represents a common route underlying the emergence of new gene functions.     

Role of mutational bias and natural selection on genome-wide nucleotide bias in prokaryotic organisms

Correlations between genomic GC contents and amino acid frequencies were studied in the homologous sequences of 12 eubacterial genomes. Results show that amino acids encoded by GC-rich codons increases significantly with genomic GC contents, whereas opposite trend was observed in case of amino acids encoded by GC-poor codons. Further studies show all the amino acids do not change in the predicted direction according to their genomic GC pressure, suggesting that protein evolution is not entirely dictated by their nucleotide frequencies. Amino acid substitution matrix calculated among hydrophobic, amphipathic and hydrophilic amino acid groups' shows that amphipathic and hydrophilic amino acids are more frequently substituted by hydrophobic amino acids than from hydrophobic to hydrophilic or amphipathic amino acids. This indicates that nucleotide bias induces a directional changes in proteome composition in such a way that underwent strong changes in hydropathy values. In fact, significant increases in hydrophobicity values have also been observed with the increase of genomic GC contents. Correlations between GC contents and amino acid compositions in three different predicted protein secondary structures show that hydropathy values increases significantly with GC contents in aperiodic and helix structures whereas strand structure remains insensitive with the genomic GC levels. The relative importance of mutation and selection on the evolution of proteins have been discussed on the basis of these results.     

Proteome-wide evidence for enhanced positive Darwinian selection within intrinsically disordered regions in proteins

Background  

Understanding the adaptive changes that alter the function of proteins during evolution is an important question for biology and medicine. The increasing number of completely sequenced genomes from closely related organisms, as well as individuals within species, facilitates systematic detection of recent selection events by means of comparative genomics.     

Unrecognized fine-scale recombination can mimic the effects of adaptive radiation

Gene sequences can undergo accelerated nucleotide changes and rapid diversification. The rapid sequence changes can then potentially lead to phylogenetic incongruence. Recently, Bodilis et al. (2011) observed artificial phylogenetic incongruence using the Pseudomonas surface protein gene oprF, and hypothesized that it was the result of a long-branch attraction artifact ultimately caused by adaptive radiation. In this study, an alternative hypothesis, namely fine-scale recombination, was tested on the same dataset. The results reveal that regions in oprF are of different evolutionary origins, and the mosaic gene structure resulted in confounding phylogenetic signals. These findings demonstrate that unrecognized fine-scale recombination can confound the phylogenetic interpretation and emphasize the limitation of using whole genes as the unit of phylogenetic analysis.     

Reciprocal relationship between stem-loop potential and substitution density in retroviral quasispecies under positive Darwinian selection

Nucleic acids have the potential to form in trastrand stem-loops if complementary bases are suitably located. Computer analyses of poliovirus and retroviral RNAs have revealed a reciprocal relationship between statistically significant stem-loop potential and sequence variability. The statistically significant stem-loop potential of a nucleic acid segment has been defined as a function of the difference between the folding energy of the natural segment (FONS) and the mean folding energy of a set of randomized (shuffled) versions of the natural segment (FORS-M). Since FONS is dependent on both base composition and base order, whereas FORS-M is solely dependent on base composition (a genomic characteristic), it follows that statistically significant stem-loop potential (FORS-D) is a function of base order (a local characteristic). In retroviral genomes, as in all DNA genomes studied, positive FORS-D values are widely distributed. Thus there have been pressures on base order both to encode specific functions and to encode stem-loops. As in the case of DNA genomes under positive Darwinian selection pressure, in HIV-1 specific function appears to dominate in rapidly evolving regions. Here high sequence variability, expressed as substitution density (not indel density), is associated with negative FORS-D values (impaired base-order-dependent stem-loop potential). This suggests that in these regions HIV-1 genomes are under positive selection pressure by host defenses. The general function of stem-loops is recombination. This is a vital process if, from among members of viral quasispecies, functional genomes are to be salvaged. Thus, for rapidly evolving RNA genomes, it is as important to conserve base-order-dependent stem-loop potential as to conserve other functions.     

Rapid evolution by positive Darwinian selection in the extracellular domain of the abundant lymphocyte protein CD45 in primates

CD45, encoded by PTPRC in humans, is the most abundantly expressed protein on the surface of many lymphocytes. We investigated whether the extracellular region of CD45 was under positive selection in Old World primates, and whether there was differential selection across this region, particularly on exons that were involved in alternative splicing and those that were not alternatively spliced. The results show extraordinarily strong and consistent positive Darwinian selection on the extracellular part of CD45 throughout the evolution of Old World monkeys, apes and humans. Positive selection is concentrated in exons 9 and 14, which code for the previously neglected linker and fibronectin III domains. These exons have a high rate of evolution at nonsynonymous sites that is roughly twice as high as that of the intronic rate in this gene. In contrast, alternatively spliced exons 4-6, which code for the variable domains, are under weaker positive selection and are evolving more slowly than the intronic rate. These data provide a striking example of positive selection in a well-known gene that should provide an impetus for further functional studies to elucidate its species-specific function.     

Generic Darwinian selection in catalytic protocell assemblies

To satisfy the minimal requirements for life, an information carrying molecular structure must be able to convert resources into building blocks and also be able to adapt to or modify its environment to enhance its own proliferation. Furthermore, new copies of itself must have variable fitness such that evolution is possible. In practical terms, a minimal protocell should be characterized by a strong coupling between its metabolism and genetic subsystem, which is made possible by the container. There is still no general agreement on how such a complex system might have been naturally selected for in a prebiotic environment. However, the historical details are not important for our investigations as they are related to assembling and evolution of protocells in the laboratory. Here, we study three different minimal protocell models of increasing complexity, all of them incorporating the coupling between a 'genetic template', a container and, eventually, a toy metabolism. We show that for any local growth law associated with template self-replication, the overall temporal evolution of all protocell's components follows an exponential growth (efficient or uninhibited autocatalysis). Thus, such a system attains exponential growth through coordinated catalytic growth of its component subsystems, independent of the replication efficiency of the involved subsystems. As exponential growth implies the survival of the fittest in a competitive environment, these results suggest that protocell assemblies could be efficient vehicles in terms of evolving through Darwinian selection.     

Rapid evolution of a primate sperm protein: relaxation of functional constraint or positive Darwinian selection?

Protamines are arginine-rich proteins that replace histones and bind sperm DNA during spermatogenesis in vertebrates. Previous studies have shown that protamine exons evolve faster than does the protamine intron. It has been suggested that this is a result of a relaxation of functional constraint. However, a more likely explanation is that the evolutionary rate of exons has been accelerated by positive Darwinian selection, because introns are generally believed to evolve in a neutral fashion. Therefore, we examined the possibility that positive selection has been acting on the protamine genes of three groups of placental mammals: primates (hominoids and Old World monkeys), rodents (mice, rats, and guinea pigs), and pecoran ruminants (deer and bovids). We found that the nucleotide substitution rate at nonsynonymous sites is significantly higher than the rate at synonymous and intron sites for protamine P1 of hominoids and Old World monkeys. This result suggests that positive selection has been operating on protamine P1 of these species. In contrast, no clear-cut evidence of positive selection was found for protamine P1 of ruminants and rodents or protamine P2 of primates. The agent of positive selection on primate protamine P1 remains unknown, though sperm competition is a possibility. Further investigations on the function and intraspecific polymorphism of this protein are needed in order to identify the selection agent.     

A small post-translocation energy bias aids nucleotide selection in T7 RNA polymerase transcription

The RNA polymerase (RNAP) of bacteriophage T7 is a single subunit enzyme that can transcribe DNA to RNA in the absence of additional protein factors. In this work, we present a model of T7 RNAP translocation during elongation. Based on structural information and experimental data from single-molecule force measurements, we show that a small component of facilitated translocation or power stroke coexists with the Brownian-ratchet-driven motions, and plays a crucial role in nucleotide selection at pre-insertion. The facilitated translocation is carried out by the conserved Tyr⁶³⁹ that moves its side chain into the active site, pushing aside the 3′-end of the RNA, and forming a locally stabilized post-translocation intermediate. Pre-insertion of an incoming nucleotide into this stabilized intermediate state ensures that Tyr⁶³⁹ closely participates in selecting correct nucleotides. A similar translocation mechanism has been suggested for multi-subunit RNAPs involving the bridge-helix bending. Nevertheless, the bent bridge-helix sterically prohibits nucleotide binding in the post-transolocation intermediate analog; moreover, the analog is not stabilized unless an inhibitory protein factor binds to the enzyme. Using our scheme, we also compared the efficiencies of different strategies for nucleotide selection, and examined effects of facilitated translocation on forward tracking.     

Interplay of Darwinian and frequency-dependent selection in the host-associated microbial populations

In order to analyze the microevolutionary processes in host-associated microorganisms, we simulated the dynamics of rhizobia populations composed of a parental strain and its mutants possessing the altered fitness within "plant-soil" system. The population dynamics was presented as a series of cycles (each one involves "soil-->rhizosphere-->nodules-->soil" succession) described using recurrent equations. For representing the selection and mutation pressures, we used a universal approach based on calculating the shifts in the genetic ratios of competing bacterial genotypes within the particular habitats and across several habitats. Analysis of the model demonstrated that a balanced polymorphism may be established in rhizobia population: mutants with an improved fitness do not supplant completely the parental strain while mutants with a decreased fitness may be maintained stably. This polymorphism is caused by a rescue of low-fitted genotypes via negative frequency-dependent selection (FDS) that is implemented during inoculation of nodules and balances the Darwinian selection that occurs during multiplication or extinction of bacteria at different habitats. The most diverse populations are formed if the rhizobia are equally successful in soil and nodules, while a marked preference for any of these habitats results in the decrease of diversity. Our simulation suggests that FDS can maintain the mutualistic rhizobia-legume interactions under the stress conditions deleterious for surviving the bacterial strains capable for intensive N2 fixation. Genetic consequences of releasing the modified rhizobia strains may be addressed using the presented model.     

Darwinian balancing selection: Predation counters sexual selection in a wild insect

The potential viability costs of sexually selected traits are central to hypotheses about the evolution of exaggerated traits. Estimates of these costs in nature can come from selection analyses using multiple components of fitness during the same time frame. For a population of tree crickets (Oecanthus nigricornis: Gryllidae), we analyzed viability and sexual selection on male traits by comparing Oecanthus prey of a solitary wasp to those that survived, and comparing mating individuals to solitary males. We measured forewing width (sexually size dimorphic and used for singing), head width, pronotum length, and size of hind jumping legs as potential targets of selection. Supporting the hypothesis that sexually selected traits have viability costs, we found that significant directional sexual selection for wider heads was opposed by significant viability selection for narrower heads. Nonlinear selection revealed that individuals with wide heads and small legs were most attractive, but individuals with narrow heads, large legs, and intermediate pronotum length were most likely to survive. Successful mating may put males at greater risk of predation, especially if copulation per se is risky. Such balancing selection in tree crickets may have constrained the evolution of sexual dimorphism in head size—a condition seen in other gryllids and orthopterans.     

POTION: an end-to-end pipeline for positive Darwinian selection detection in genome-scale data through phylogenetic comparison of protein-coding genes

Background

Detection of genes evolving under positive Darwinian evolution in genome-scale data is nowadays a prevailing strategy in comparative genomics studies to identify genes potentially involved in adaptation processes. Despite the large number of studies aiming to detect and contextualize such gene sets, there is virtually no software available to perform this task in a general, automatic, large-scale and reliable manner. This certainly occurs due to the computational challenges involved in this task, such as the appropriate modeling of data under analysis, the computation time to perform several of the required steps when dealing with genome-scale data and the highly error-prone nature of the sequence and alignment data structures needed for genome-wide positive selection detection.

Results

We present POTION, an open source, modular and end-to-end software for genome-scale detection of positive Darwinian selection in groups of homologous coding sequences. Our software represents a key step towards genome-scale, automated detection of positive selection, from predicted coding sequences and their homology relationships to high-quality groups of positively selected genes. POTION reduces false positives through several sophisticated sequence and group filters based on numeric, phylogenetic, quality and conservation criteria to remove spurious data and through multiple hypothesis corrections, and considerably reduces computation time thanks to a parallelized design. Our software achieved a high classification performance when used to evaluate a curated dataset of Trypanosoma brucei paralogs previously surveyed for positive selection. When used to analyze predicted groups of homologous genes of 19 strains of Mycobacterium tuberculosis as a case study we demonstrated the filters implemented in POTION to remove sources of errors that commonly inflate errors in positive selection detection. A thorough literature review found no other software similar to POTION in terms of customization, scale and automation.

Conclusion

To the best of our knowledge, POTION is the first tool to allow users to construct and check hypotheses regarding the occurrence of site-based evidence of positive selection in non-curated, genome-scale data within a feasible time frame and with no human intervention after initial configuration. POTION is available at http://www.lmb.cnptia.embrapa.br/share/POTION/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1765-0) contains supplementary material, which is available to authorized users.     

Dntf-2r,a young Drosophila retroposed gene with specific male expression under positive Darwinian selection

A direct approach to investigating new gene origination is to examine recently evolved genes. We report a new gene in the Drosophila melanogaster subgroup, Drosophila nuclear transport factor-2-related (Dntf-2r). Its sequence features and phylogenetic distribution indicate that Dntf-2r is a retroposed functional gene and originated in the common ancestor of D. melanogaster, D. simulans, D. sechellia, and D. mauritiana, within the past 3-12 million years (MY). Dntf-2r evolved more rapidly than the parental gene, under positive Darwinian selection as revealed by the McDonald-Kreitman test and other evolutionary analyses. Comparative expression analysis shows that Dntf-2r is male specific whereas the parental gene, Dntf-2, is widely expressed in D. melanogaster. In agreement with its new expression pattern, the Dntf-2r putative promoter sequence is similar to the late testis promoter of beta2-tubulin. We discuss the possibility that the action of positive selection in Dntf-2r is related to its putative male-specific functions.