Estimates of positive Darwinian selection are inflated by errors in sequencing, annotation, and alignment.
Ontology highlight
ABSTRACT:
SUBMITTER: Schneider A
PROVIDER: S-EPMC2817407 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Current estimates of diversifying positive selection rely on first having an accurate multiple sequence alignment. Simulation studies have shown that under biologically plausible conditions, relying on a single estimate of the alignment from commonly used alignment software can lead to unacceptably high false-positive rates in detecting diversifying positive selection. We present a novel statistical method that eliminates excess false positives resulting from alignment error by jointly estimating the degree of positive selection and the alignment under an evolutionary model. Our model treats both substitutions and insertions/deletions as sequence changes on a tree and allows site heterogeneity in the substitution process. We conduct inference starting from unaligned sequence data by integrating over all alignments. This approach naturally accounts for ambiguous alignments without requiring ambiguously aligned sites to be identified and removed prior to analysis. We take a Bayesian approach and conduct inference using Markov chain Monte Carlo to integrate over all alignments on a fixed evolutionary tree topology. We introduce a Bayesian version of the branch-site test and assess the evidence for positive selection using Bayes factors. We compare two models of differing dimensionality using a simple alternative to reversible-jump methods. We also describe a more accurate method of estimating the Bayes factor using Rao-Blackwellization. We then show using simulated data that jointly estimating the alignment and the presence of positive selection solves the problem with excessive false positives from erroneous alignments and has nearly the same power to detect positive selection as when the true alignment is known. We also show that samples taken from the posterior alignment distribution using the software BAli-Phy have substantially lower alignment error compared with MUSCLE, MAFFT, PRANK, and FSA alignments.
| S-EPMC4155473 | biostudies-literature
Project description:Genomic imprinting is an epigenetic phenomenon where autosomal genes display uniparental expression depending on whether they are maternally or paternally inherited. Genomic imprinting can arise from parental conflicts over resource allocation to the offspring, which could drive imprinted loci to evolve by positive selection. We investigate whether positive selection is associated with genomic imprinting in the inbreeding species Arabidopsis thaliana. Our analysis of 140 genes regulated by genomic imprinting in the A. thaliana seed endosperm demonstrates they are evolving more rapidly than expected. To investigate whether positive selection drives this evolutionary acceleration, we identified orthologs of each imprinted gene across 34 plant species and elucidated their evolutionary trajectories. Increased positive selection was sought by comparing its incidence among imprinted genes with nonimprinted controls. Strikingly, we find a statistically significant enrichment of imprinted paternally expressed genes (iPEGs) evolving under positive selection, 50.6% of the total, but no such enrichment for positive selection among imprinted maternally expressed genes (iMEGs). This suggests that maternally- and paternally expressed imprinted genes are subject to different selective pressures. Almost all positively selected amino acids were fixed across 80 sequenced A. thaliana accessions, suggestive of selective sweeps in the A. thaliana lineage. The imprinted genes under positive selection are involved in processes important for seed development including auxin biosynthesis and epigenetic regulation. Our findings support a genomic imprinting model for plants where positive selection can affect paternally expressed genes due to continued conflict with maternal sporophyte tissues, even when parental conflict is reduced in predominantly inbreeding species.
| S-EPMC6526901 | biostudies-literature
Project description:Comparisons of intron-exon structures between homologous genes in different eukaryotic species have revealed substantial variation in the number of introns. These observations imply that, in each case, an intron presence-absence polymorphism must have existed in the past. Such a polymorphism, created by a recent intron-loss mutation, is reported here in a eukaryotic organism. This gene structure, detected in the jingwei (jgw) gene, segregates at high frequency (77%) in natural populations of Drosophila teissieri and is associated with a marked change in mRNA levels. Furthermore, the intron loss does not result from a mRNA-mediated mechanism as is usually proposed, but from a partial deletion at the DNA level that also results in the addition of four new amino acids to the JGW protein. Population genetic analyses of the pattern of nucleotide variation surrounding the intron polymorphism indicate the action of positive Darwinian selection on the intron-absent variant. Forward simulations suggest that the intensity of this selection is weak to moderate, roughly equal to the selection intensity on most replacement mutations in Drosophila.
| S-EPMC123031 | biostudies-literature
Project description:Trace amine-associated receptors (TAARs) in mammals recently have been shown to function as olfactory receptors. We have delineated the taar gene family in jawless, cartilaginous, and bony fish (zero, 2, and >100 genes, respectively). We conclude that taar genes are evolutionary much younger than the related OR and ORA/V1R olfactory receptor families, which are present already in lamprey, a jawless vertebrate. The 2 cartilaginous fish genes appear to be ancestral for 2 taar classes, each with mammalian and bony fish (teleost) representatives. Unexpectedly, a whole new clade, class III, of taar genes originated even later, within the teleost lineage. Taar genes from all 3 classes are expressed in subsets of zebrafish olfactory receptor neurons, supporting their function as olfactory receptors. The highly conserved TAAR1 (shark, mammalian, and teleost orthologs) is not expressed in the olfactory epithelium and may constitute the sole remnant of a primordial, nonolfactory function of this family. Class III comprises three-fourths of all teleost taar genes and is characterized by the complete loss of the aminergic ligand-binding motif, stringently conserved in the other 2 classes. Two independent intron gains in class III taar genes represent extraordinary evolutionary dynamics, considering the virtual absence of intron gains during vertebrate evolution. The d(N)/d(S) analysis suggests both minimal global negative selection and an unparalleled degree of local positive selection as another hallmark of class III genes. The accelerated evolution of class III teleost taar genes conceivably might mark the birth of another olfactory receptor gene family.
| S-EPMC2657432 | biostudies-literature
Project description:The relative rates of nucleotide substitution at synonymous and nonsynonymous sites within protein-coding regions have been widely used to infer the action of natural selection from comparative sequence data. It is known, however, that mutational and repair biases can affect rates of evolution at both synonymous and nonsynonymous sites. More importantly, it is also known that synonymous sites are particularly prone to the effects of nucleotide bias. This means that nucleotide biases may affect the calculated ratio of substitution rates at synonymous and nonsynonymous sites. Using a large data set of animal mitochondrial sequences, we demonstrate that this is, in fact, the case. Highly biased nucleotide sequences are characterized by significantly elevated dN/dS ratios, but only when the nucleotide frequencies are not taken into account. When the analysis is repeated taking the nucleotide frequencies at each codon position into account, such elevated ratios disappear. These results suggest that the recently reported differences in dN/dS ratios between vertebrate and invertebrate mitochondrial sequences could be explained by variations in mitochondrial nucleotide frequencies rather than the effects of positive Darwinian selection.
| S-EPMC2582980 | biostudies-literature
Project description:Rapid evolution driven by positive Darwinian selection is a recurrent theme in male reproductive protein evolution. In contrast, positive selection has never been demonstrated for female reproductive proteins. Here, we perform phylogeny-based tests on three female mammalian fertilization proteins and demonstrate positive selection promoting their divergence. Two of these female fertilization proteins, the zona pellucida glycoproteins ZP2 and ZP3, are part of the mammalian egg coat. Several sites identified in ZP3 as likely to be under positive selection are located in a region previously demonstrated to be involved in species-specific sperm-egg interaction, suggesting the selective pressure is related to male-female interaction. The results provide long-sought evidence for two evolutionary hypotheses: sperm competition and sexual conflict.
| S-EPMC30168 | biostudies-literature
Project description:Members of the Family Bufonidae, true toads, are famous for their endogenously synthesized cardioactive steroids that serve as defensive toxins. Evolution of resistance to these toxins is not understood. We sequenced a key region of the toxin's binding site in the Na(+)/K(+) ATPase for relevant taxa representing Hyloidea (including bufonids), Ranoidea and Archaeobatrachia and tested for positive selection in a phylogenetic context. Bufonidae were distinct from other Hyloidea at 4-6 of 12 sites and, with one exception, had a homologous amino acid sequence. Melanophryniscus stelzneri had a distinct sequence, consistent with other independent evidence for a differentiated toxin. Tests within Bufonidae detected positive selection within the binding region, providing, to our knowledge, the first evidence of this type for positive selection within Amphibia. There was no evidence for positive selection on Bufonidae or M. stelzneri lineages. Sequence change in Leptodactylus ocellatus, a leptodactylid predator of Bufonidae, provides a molecular basis for predator resistance possibly associated with gene duplication.
| S-EPMC2781935 | biostudies-literature
Project description: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.
| S-EPMC2478720 | biostudies-literature
Project description:BACKGROUND:Molecular evolution is a very active field of research, with several complementary approaches, including dN/dS, HON90, MM01, and others. Each has documented strengths and weaknesses, and no one approach provides a clear picture of how natural selection works at the molecular level. The purpose of this work is to present a simple new method that uses quantitative amino acid properties to identify and characterize directional selection in proteins. METHODS:Inferred amino acid replacements are viewed through the prism of a single physicochemical property to determine the amount and direction of change caused by each replacement. This allows the calculation of the probability that the mean change in the single property associated with the amino acid replacements is equal to zero (H0: ? = 0; i.e., no net change) using a simple two-tailed t-test. RESULTS:Example data from calanoid and cyclopoid copepod cytochrome oxidase subunit I sequence pairs are presented to demonstrate how directional selection may be linked to major shifts in adaptive zones, and that convergent evolution at the whole organism level may be the result of convergent protein adaptations. CONCLUSIONS:Rather than replace previous methods, this new method further complements existing methods to provide a holistic glimpse of how natural selection shapes protein structure and function over evolutionary time.
| S-EPMC3849801 | biostudies-literature
Project description: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. RESULTS: We have used genome-wide strain-specific single nucleotide polymorphism data from 64 strains of budding yeast (Saccharomyces cerevisiae or Saccharomyces paradoxus) to determine whether adaptive positive selection is correlated with protein regions showing propensity for different classes of structure conformation. Data from phylogenetic and population genetic analysis of 3,746 gene alignments consistently shows a significantly higher degree of positive Darwinian selection in intrinsically disordered regions of proteins compared to regions of alpha helix, beta sheet or tertiary structure. Evidence of positive selection is significantly enriched in classes of proteins whose functions and molecular mechanisms can be coupled to adaptive processes and these classes tend to have a higher average content of intrinsically unstructured protein regions. CONCLUSIONS: We suggest that intrinsically disordered protein regions may be important for the production and maintenance of genetic variation with adaptive potential and that they may thus be of central significance for the evolvability of the organism or cell in which they occur.
| S-EPMC3218827 | biostudies-other