Detection Of Positive Selection Research Articles

Detection of positive selection eliminating effects of structural constraints in hemagglutinin of H3N2 human influenza A virus

In the evolutionary studies of proteins, the average effect of natural selection operating on amino acid mutations may be examined by comparing the numbers of synonymous (dS) and nonsynonymous (dN) substitutions that have accumulated during the same time period. In this method, destabilizing mutations occurring across protein molecules may interfere with detection of natural selection, particularly positive selection, operating on other mutations. Here an attempt to detect positive selection eliminating effects of structural constraints is demonstrated using hemagglutinin (HA) of H3N2 human influenza A virus as an example. Compatible and incompatible amino acids were inferred at each site from the computational analysis of three-dimensional structure using the thermodynamic stability as an indicator, and natural selection was examined by comparing dS and dN among compatible amino acids. In the analysis of 2701 nucleotide sequences for the entire coding region of HA, the new method identified twice as many positively selected amino acid sites as the ordinary method (16 and 4 sites in the former method without and with correction for multiple testing, respectively, and 8 and 2 sites in the latter method). Positively selected sites were involved in epitopes, receptor-binding pocket, epistasis, and stabilization, which appeared to be biologically reasonable. Nevertheless, there still appeared to be several problems, which may largely render this method conservative. It may be effective to analyze many densely sampled sequences in this method.

Cross-Species Transmission in the Speciation of the Currently Known Murinae-Associated Hantaviruses

To gain more insight into the phylogeny of Dabieshan virus (DBSV), carried by Niviventer confucianus and other Murinae-associated hantaviruses, genome sequences of novel variants of DBSV were recovered from Niviventer rats trapped in the mountainous areas of Wenzhou, China. Genetic analyses show that all known genetic variants of DBSV, including the ones identified in this study, are distinct from other Murinae-associated hantaviruses. DBSV variants show geographic clustering and high intraspecies diversity. The data suggest that DBSV is a distinct species in the genus Hantavirus. Interestingly, DBSV shows the highest sequence identity to Hantaan virus (HTNV), with a >7% difference in the sequences of the N, GPC, and L proteins, while N. confucianus is more closely related to Rattus norvegicus (the host of Seoul virus [SEOV]) than to Apodemus agrarius (the host of HTNV and Saaremaa virus [SAAV]). Further genetic analyses of all known Murinae-associated hantaviruses (both established and tentative species) show that many of them, including DBSV, may have originated from host switching. The estimation of evolutionary rates and divergence time supports the role of cross-species transmission in the evolution of Murinae-associated hantaviruses. The detection of positive selection suggests that genetic drift may contribute to the speciation of Murinae-associated hantaviruses and that adaptation has a role as well.

Evolution of MHC class I genes in the European badger (Meles meles)

The major histocompatibility complex (MHC) plays a central role in the adaptive immune system and provides a good model with which to understand the evolutionary processes underlying functional genes. Trans-species polymorphism and orthology are both commonly found in MHC genes; however, mammalian MHC class I genes tend to cluster by species. Concerted evolution has the potential to homogenize different loci, whereas birth-and-death evolution can lead to the loss of orthologs; both processes result in monophyletic groups within species. Studies investigating the evolution of MHC class I genes have been biased toward a few particular taxa and model species. We present the first study of MHC class I genes in a species from the superfamily Musteloidea. The European badger (Meles meles) exhibits moderate variation in MHC class I sequences when compared to other carnivores. We identified seven putatively functional sequences and nine pseudogenes from genomic (gDNA) and complementary (cDNA) DNA, signifying at least two functional class I loci. We found evidence for separate evolutionary histories of the α1 and α2/α3 domains. In the α1 domain, several sequences from different species were more closely related to each other than to sequences from the same species, resembling orthology or trans-species polymorphism. Balancing selection and probable recombination maintain genetic diversity in the α1 domain, evidenced by the detection of positive selection and a recombination event. By comparison, two recombination breakpoints indicate that the α2/α3 domains have most likely undergone concerted evolution, where recombination has homogenized the α2/α3 domains between genes, leading to species-specific clusters of sequences. Our findings highlight the importance of analyzing MHC domains separately.

The Effects of Alignment Error and Alignment Filtering on the Sitewise Detection of Positive Selection

When detecting positive selection in proteins, the prevalence of errors resulting from misalignment and the ability of alignment filters to mitigate such errors are not well understood, but filters are commonly applied to try to avoid false positive results. Focusing on the sitewise detection of positive selection across a wide range of divergence levels and indel rates, we performed simulation experiments to quantify the false positives and false negatives introduced by alignment error and the ability of alignment filters to improve performance. We found that some aligners led to many false positives, whereas others resulted in very few. False negatives were a problem for all aligners, increasing with sequence divergence. Of the aligners tested, PRANK's codon-based alignments consistently performed the best and ClustalW performed the worst. Of the filters tested, GUIDANCE performed the best and Gblocks performed the worst. Although some filters showed good ability to reduce the error rates from ClustalW and MAFFT alignments, none were found to substantially improve the performance of PRANK alignments under most conditions. Our results revealed distinct trends in error rates and power levels for aligners and filters within a biologically plausible parameter space. With the best aligner, a low false positive rate was maintained even with extremely divergent indel-prone sequences. Controls using the true alignment and an optimal filtering method suggested that performance improvements could be gained by improving aligners or filters to reduce the prevalence of false negatives, especially at higher divergence levels and indel rates.

Lack of Low Frequency Variants Masks Patterns of Non-Neutral Evolution following Domestication

Detecting artificial selection in the genome of domesticated species can not only shed light on human history but can also be beneficial to future breeding strategies. Evidence for selection has been documented in domesticated species including maize and rice, but few studies have to date detected signals of artificial selection in the Sorghum bicolor genome. Based on evidence that domesticated S. bicolor and its wild relatives show significant differences in endosperm structure and quality, we sequenced three candidate seed storage protein (kafirin) loci and three candidate starch biosynthesis loci to test whether these genes show non-neutral evolution resulting from the domestication process. We found strong evidence of non-neutral selection at the starch synthase IIa gene, while both starch branching enzyme I and the beta kafirin gene showed weaker evidence of non-neutral selection. We argue that the power to detect consistent signals of non-neutral selection in our dataset is confounded by the absence of low frequency variants at four of the six candidate genes. A future challenge in the detection of positive selection associated with domestication in sorghum is to develop models that can accommodate for skewed frequency spectrums.

Targeted resequencing of a genomic region influencing tameness and aggression reveals multiple signals of positive selection

The identification of the causative genetic variants in quantitative trait loci (QTL) influencing phenotypic traits is challenging, especially in crosses between outbred strains. We have previously identified several QTL influencing tameness and aggression in a cross between two lines of wild-derived, outbred rats (Rattus norvegicus) selected for their behavior towards humans. Here, we use targeted sequence capture and massively parallel sequencing of all genes in the strongest QTL in the founder animals of the cross. We identify many novel sequence variants, several of which are potentially functionally relevant. The QTL contains several regions where either the tame or the aggressive founders contain no sequence variation, and two regions where alternative haplotypes are fixed between the founders. A re-analysis of the QTL signal showed that the causative site is likely to be fixed among the tame founder animals, but that several causative alleles may segregate among the aggressive founder animals. Using a formal test for the detection of positive selection, we find 10 putative positively selected regions, some of which are close to genes known to influence behavior. Together, these results show that the QTL is probably not caused by a single selected site, but may instead represent the joint effects of several sites that were targets of polygenic selection.

Back to the sea twice: identifying candidate plant genes for molecular evolution to marine life

BackgroundSeagrasses are a polyphyletic group of monocotyledonous angiosperms that have adapted to a completely submerged lifestyle in marine waters. Here, we exploit two collections of expressed sequence tags (ESTs) of two wide-spread and ecologically important seagrass species, the Mediterranean seagrass Posidonia oceanica (L.) Delile and the eelgrass Zostera marina L., which have independently evolved from aquatic ancestors. This replicated, yet independent evolutionary history facilitates the identification of traits that may have evolved in parallel and are possible instrumental candidates for adaptation to a marine habitat.ResultsIn our study, we provide the first quantitative perspective on molecular adaptations in two seagrass species. By constructing orthologous gene clusters shared between two seagrasses (Z. marina and P. oceanica) and eight distantly related terrestrial angiosperm species, 51 genes could be identified with detection of positive selection along the seagrass branches of the phylogenetic tree. Characterization of these positively selected genes using KEGG pathways and the Gene Ontology uncovered that these genes are mostly involved in translation, metabolism, and photosynthesis.ConclusionsThese results provide first insights into which seagrass genes have diverged from their terrestrial counterparts via an initial aquatic stage characteristic of the order and to the derived fully-marine stage characteristic of seagrasses. We discuss how adaptive changes in these processes may have contributed to the evolution towards an aquatic and marine existence.

Significant Selective Constraint at 4-Fold Degenerate Sites in the Avian Genome and Its Consequence for Detection of Positive Selection

A major conclusion from comparative genomics is that many sequences that do not code for proteins are conserved beyond neutral expectations, indicating that they evolve under the influence of purifying selection and are likely to have functional roles. Due to the degeneracy of the genetic code, synonymous sites within protein-coding genes have previously been seen as “silent” with respect to function and thereby invisible to selection. However, there are indications that synonymous sites of vertebrate genomes are also subject to selection and this is not necessarily because of potential codon bias. We used divergence in ancestral repeats as a neutral reference to estimate the constraint on 4-fold degenerate sites of avian genes in a whole-genome approach. In the pairwise comparison of chicken and zebra finch, constraint was estimated at 24–32%. Based on three-species alignments of chicken, turkey, and zebra finch, lineage-specific estimates of constraint were 43%, 29%, and 24%, respectively. The finding of significant constraint at 4-fold degenerate sites from data on interspecific divergence was replicated in an analysis of intraspecific diversity in the chicken genome. These observations corroborate recent data from mammalian genomes and call for a reappraisal of the use of synonymous substitution rates as neutral standards in molecular evolutionary analysis, for example, in the use of the well-known dN/dS ratio and in inferences on positive selection. We show by simulations that the rate of false positives in the detection of positively selected genes and sites increases several-fold at the levels of constraint at 4-fold degenerate sites found in this study.

Genome scan in the mosquito Aedes rusticus: population structure and detection of positive selection after insecticide treatment

Identification of genes involved in local adaptation is particularly challenging for species functioning as a network of interconnected populations undergoing frequent extinctions-recolonizations, because populations are submitted to contrasted evolutionary pressures. Using amplified fragment length polymorphism markers, population genetic structure of the mosquito Aedes rusticus was analysed in five geographical areas of the French Rhône-Alpes region. We included a number of sites that were treated with the bio-insecticide Bacillus thuringiensis israelensis (Bti) for more than 15 years. Analysis of molecular variance revealed that most of the genetic variability was found within populations (96%), with no significant variation among geographical areas, although variation among populations within areas (4%) was significant. The global genetic differentiation index F(ST) was low (0.0366 +/- 0.167). However, pairwise F(ST) values were significant and no isolation-by-distance at the regional level was observed, suggesting a metapopulation structure in this species. Bti-treatment had no effect on genetic structure and on within-population genetic diversity. Potential signatures of positive selection associated with Bti-treatment were detected for five loci, even though toxicological bioassays performed on field-collected larvae showed no significant difference in mortality between Bti-treated and nontreated sites. The difficulty of detecting moderate resistance in field-collected larvae together with possible differential persistence of toxins in the environment may explain our inability to detect a toxicological response to Bti in treated sites. The evidence for positive selection occurring at several genomic regions suggests a first step towards Bti resistance in natural mosquito populations treated with this bio-insecticide. Furthermore, this signal was detectable using genomic tools before any toxicological evidence for resistance could be identified.

Response to Yang et al.: Problems with Bayesian methods of detecting positive selection at the DNA sequence level

It is unfortunate that Yang et al. (1) misrepresent the contents and conclusions of our recent study (2). Our study was motivated by the recent publication of papers reporting the detection of positive selection in many genes of the human lineage by using the branch-site method (BSM) and others, where the ratio (ω) of nonsynonymous to synonymous nucleotide substitutions is used as a measure of positive selection (e.g., ref. 3). Because the number of nucleotide substitutions per gene (n) in the human lineage after separation from the chimpanzee is very small (n ≈ 4) and BSM depends on the large sample theory, we examined the theoretical basis of BSM in comparison with the small-sample method (SSM) using Fisher's exact test (4), which gives accurate statistical tests for small samples.

Gene Conversion Among Paralogs Results in Moderate False Detection of Positive Selection Using Likelihood Methods

Previous studies have shown that recombination between allelic sequences can cause likelihood-based methods for detecting positive selection to produce many false-positive results. In this article, we use simulations to study the impact of nonallelic gene conversion on the specificity of PAML to detect positive selection among gene duplicates. Our results show that, as expected, gene conversion leads to higher rates of false-positive results, although only moderately. These rates increase with the genetic distance between sequences, the length of converted tracts, and when no outgroup sequences are included in the analysis. We also find that branch-site models will incorrectly identify unconverted sequences as the targets of positive selection when their close paralogs are converted. Bayesian prediction of sites undergoing adaptive evolution implemented in PAML is affected by conversion, albeit in a less straightforward way. Our work suggests that particular attention should be devoted to the evolutionary analysis of recent duplicates that may have experienced gene conversion because they may provide false signals of positive selection. Fortunately, these results also imply that those cases most susceptible to false-positive results--i.e., high divergence between paralogs, long conversion tracts--are also the cases where detecting gene conversion is the easiest.

Statistical Comparison of Nucleotide, Amino Acid, and Codon Substitution Models for Evolutionary Analysis of Protein-Coding Sequences

Statistical models for the evolution of molecular sequences play an important role in the study of evolutionary processes. For the evolutionary analysis of protein-coding sequences, 3 types of evolutionary models are available: 1) nucleotide, 2) amino acid, and 3) codon substitution models. Selecting appropriate models can greatly improve the estimation of phylogenies and divergence times and the detection of positive selection. Although much attention has been paid to the comparisons among the same types of models, relatively little attention has been paid to the comparisons among the different types of models. Additionally, because such models have different data structures, comparison of those models using conventional model selection criteria such as Akaike information criterion (AIC) or Bayesian information criterion (BIC) is not straightforward. Here, we suggest new procedures to convert models of the above-mentioned 3 types to 64-dimensional models with nucleotide triplet substitution. These conversion procedures render it possible to statistically compare the models of these 3 types by using AIC or BIC. By analyzing divergent and conserved interspecific mammalian sequences and intraspecific human population data, we show the superiority of the codon substitution models and discuss the advantages and disadvantages of the models of the 3 types.

Genetic analysis of foot-and-mouth disease virus serotype A of Indian origin and detection of positive selection and recombination in leader protease-and capsid-coding regions.

The leader protease (Lpro) and capsid-coding sequences (P1) constitute approximately 3 kb of the foot-and-mouth disease virus (FMDV). We studied the phylogenetic relationship of 46 FMDV serotype A isolates of Indian origin collected during the period 1968–2005 and also eight vaccine strains using the neighbour-joining tree and Bayesian tree methods. The viruses were categorized under three major groups — Asian, Euro-South American and European. The Indian isolates formed a distinct genetic group among the Asian isolates. The Indian isolates were further classified into different genetic subgroups (<5% divergence). Post-1995 isolates were divided into two subgroups while a few isolates which originated in the year 2005 from Andhra Pradesh formed a separate group. These isolates were closely related to the isolates of the 1970s. The FMDV isolates seem to undergo reverse mutation or convergent evolution wherein sequences identical to the ancestors are present in the isolates in circulation. The eight vaccine strains included in the study were not related to each other and belonged to different genetic groups. Recombination was detected in the Lpro region in one isolate (A IND 20/82) and in the VP1 coding 1D region in another isolate (A RAJ 21/96). Positive selection was identified at aa positions 23 in the Lpro (P<0.05; 0.046*) and at aa 171 in the capsid protein VP1 (P<0.01; 0.003**).

Positive Selection in Prion Protein

The prion diseases, such as Creutzfeldt-Jakob disease of humans and bovine spongiform encephalopathy, involve the aberrant metabolism and accumulation of prion protein PrP. There are three contradictory hypotheses about evolution of prion protein gene PRNP. Population genetic studies have proposed that PRNP could be under balancing selection, strong purifying selection, or mainly positive selection. We made use of the maximum likelihood tests for detection of positive selection at the amino acid level and present availability of PRNP coding sequences to contribute to these disagreements. Positive selection could occur at amino acids residing in active sites, and at amino acids involved in protein-protein interactions. Thus we tested a hypothesis that positive selection at the amino acid level in PrP might have taken place in human and related species from the superordinal group Euarchonta, as well as in bovine and related species from the superordinal clade Laurasiatheria. Our study and the present experimental evidences indicate that positive selection at the amino acid level might have taken place in the PrP signal sequences and conformationally plastic PrP regions, as well as at the protein X binding sites.

Molecular Evolution of Immune Genes in the Malaria Mosquito Anopheles gambiae

BackgroundAs pathogens that circumvent the host immune response are favoured by selection, so are host alleles that reduce parasite load. Such evolutionary processes leave their signature on the genes involved. Deciphering modes of selection operating on immune genes might reveal the nature of host-pathogen interactions and factors that govern susceptibility in host populations. Such understanding would have important public health implications.Methodology/FindingsWe analyzed polymorphisms in four mosquito immune genes (SP14D1, GNBP, defensin, and gambicin) to decipher selection effects, presumably mediated by pathogens. Using samples of Anopheles arabiensis, An. quadriannulatus and four An. gambiae populations, as well as published sequences from other Culicidae, we contrasted patterns of polymorphisms between different functional units of the same gene within and between populations. Our results revealed selection signatures operating on different time scales. At the most recent time scale, within-population diversity revealed purifying selection. Between populations and between species variation revealed reduced differentiation (GNBP and gambicin) at coding vs. noncoding- regions, consistent with balancing selection. McDonald-Kreitman tests between An. quadriannulatus and both sibling species revealed higher fixation rate of synonymous than nonsynonymous substitutions (GNBP) in accordance with frequency dependent balancing selection. At the longest time scale (>100 my), PAML analysis using distant Culicid taxa revealed positive selection at one codon in gambicin. Patterns of genetic variation were independent of exposure to human pathogens.Significance and ConclusionsPurifying selection is the most common form of selection operating on immune genes as it was detected on a contemporary time scale on all genes. Selection for “hypervariability” was not detected, but negative balancing selection, detected at a recent evolutionary time scale between sibling species may be rather common. Detection of positive selection at the deepest evolutionary time scale suggests that it occurs infrequently, possibly in association with speciation events. Our results provided no evidence to support the hypothesis that selection was mediated by pathogens that are transmitted to humans.

Models of coding sequence evolution

Probabilistic models of sequence evolution are in widespread use in phylogenetics and molecular sequence evolution. These models have become increasingly sophisticated and combined with statistical model comparison techniques have helped to shed light on how genes and proteins evolve. Models of codon evolution have been particularly useful, because, in addition to providing a significant improvement in model realism for protein-coding sequences, codon models can also be designed to test hypotheses about the selective pressures that shape the evolution of the sequences. Such models typically assume a phylogeny and can be used to identify sites or lineages that have evolved adaptively. Recently some of the key assumptions that underlie phylogenetic tests of selection have been questioned, such as the assumption that the rate of synonymous changes is constant across sites or that a single phylogenetic tree can be assumed at all sites for recombining sequences. While some of these issues have been addressed through the development of novel methods, others remain as caveats that need to be considered on a case-by-case basis. Here, we outline the theory of codon models and their application to the detection of positive selection. We review some of the more recent developments that have improved their power and utility, laying a foundation for further advances in the modeling of coding sequence evolution.

Detection of positive selection in the major capsid protein VP60 of the rabbit haemorrhagic disease virus (RHDV)

Mutations were analysed in the major capsid protein VP60 of the rabbit haemorrhagic disease virus (RHDV), a calicivirus responsible for high mortality rates in both wild and domestic European rabbits (Oryctolagus cuniculus). Likelihood of positive selection was estimated using the PAML software applied to 43 non-identical complete sequences of the major capsid protein. Three codons showed signs of positive selection (with posterior probabilities over 95%), one of them is located in the region containing the major antigenic determinants (region E). The presence of positively selected codons (PSCs) in other regions may suggest the existence of other antigenic regions on the major capsid protein that stimulate protective immune responses. At all the 3 PSCs, variation contributes to putative N-glycosylation sites of the protein. An N-glycosylation site is deleted in the non-pathogenic strain RCV. Some of the substitutions at PSCs may alter the polarity and the charge of the protein with possible implications in the protein structure and host interaction. The detection of PSCs should allow a better understanding of the interaction between RHDV and the rabbit immune system.

Hitchhiking Both Ways: Effect of Two Interfering Selective Sweeps on Linked Neutral Variation

The neutral polymorphism pattern in the vicinity of a selective sweep can be altered by both stochastic and deterministic factors. Here, we focus on the impact of another selective sweep in the region of influence of a first one. We study the signature left on neutral polymorphism by positive selection at two closely linked loci, when both beneficial mutations reach fixation. We show that, depending on the timing of selective sweeps and on their selection coefficients, the two hitchhiking effects can interfere with each other, leading to less reduction in heterozygosity than a single selective sweep of the same magnitude and more importantly to an excess of intermediate-frequency variants relative to neutrality under some parameter values. This pattern can be sustained and potentially alter the detection of positive selection, including by provoking spurious detection of balancing selection. In situations where positive selection is suspected a priori at several closely linked loci, the polymorphism pattern in the region may also be informative about their selective histories.

Improved methods for detecting selection by mutation analysis of Ig V region sequences

Statistical methods based on the relative frequency of replacement mutations in B lymphocyte Ig V region sequences have been widely used to detect the forces of selection that shape the B cell repertoire. However, current methods produce an unexpectedly high frequency of false positives and are sensitive to intrinsic biases of somatic hypermutation that can give the appearance of selection. The new statistical test proposed here provides a better trade-off between sensitivity and specificity compared with previous approaches. The low specificity of existing methods was shown in silico to result from an interaction between the effects of positive and negative selection. False detection of positive selection was confirmed in vivo through a re-analysis of published sequence data from diffuse large B cell lymphomas, highlighting the need for re-analysis of some existing studies. The sensitivity of the proposed method to detect selection was validated using new Ig transgenic mouse models in which positive selection was expected to be a significant force, as well as with a simulation-based approach. Previous concerns that intrinsic biases of somatic hypermutation could give the appearance of selection were addressed by extending the current mutation models to more fully account for the impact of microsequence on relative mutability and to include transition bias. High specificity was confirmed using a large set of non-productively rearranged Ig sequences. These results show that selection can be detected in vivo with high specificity using the new method proposed here, allowing greater insight into the existence and direction of antigen-driven selection.

Whole-Gene Positive Selection, Elevated Synonymous Substitution Rates, Duplication, and Indel Evolution of the Chloroplast clpP1 Gene

BackgroundSynonymous DNA substitution rates in the plant chloroplast genome are generally relatively slow and lineage dependent. Non-synonymous rates are usually even slower due to purifying selection acting on the genes. Positive selection is expected to speed up non-synonymous substitution rates, whereas synonymous rates are expected to be unaffected. Until recently, positive selection has seldom been observed in chloroplast genes, and large-scale structural rearrangements leading to gene duplications are hitherto supposed to be rare.Methodology/Principle FindingsWe found high substitution rates in the exons of the plastid clpP1 gene in Oenothera (the Evening Primrose family) and three separate lineages in the tribe Sileneae (Caryophyllaceae, the Carnation family). Introns have been lost in some of the lineages, but where present, the intron sequences have substitution rates similar to those found in other introns of their genomes. The elevated substitution rates of clpP1 are associated with statistically significant whole-gene positive selection in three branches of the phylogeny. In two of the lineages we found multiple copies of the gene. Neighboring genes present in the duplicated fragments do not show signs of elevated substitution rates or positive selection. Although non-synonymous substitutions account for most of the increase in substitution rates, synonymous rates are also markedly elevated in some lineages. Whereas plant clpP1 genes experiencing negative (purifying) selection are characterized by having very conserved lengths, genes under positive selection often have large insertions of more or less repetitive amino acid sequence motifs.Conclusions/SignificanceWe found positive selection of the clpP1 gene in various plant lineages to correlated with repeated duplication of the clpP1 gene and surrounding regions, repetitive amino acid sequences, and increase in synonymous substitution rates. The present study sheds light on the controversial issue of whether negative or positive selection is to be expected after gene duplications by providing evidence for the latter alternative. The observed increase in synonymous substitution rates in some of the lineages indicates that the detection of positive selection may be obscured under such circumstances. Future studies are required to explore the functional significance of the large inserted repeated amino acid motifs, as well as the possibility that synonymous substitution rates may be affected by positive selection.