Cited by 802Evolution and Genetic DynamicsGenomics and Phylogenetic StudiesRNA and protein synthesis mechanisms

RELAX: Detecting Relaxed Selection in a Phylogenetic Framework

Open access
Wertheim, Joel O., Murrell, Ben, Smith, Martin D., Kosakovsky Pond, Sergei L., Scheffler, Konrad·2015·Molecular Biology and Evolution

Abstract

Relaxation of selective strength, manifested as a reduction in the efficiency or intensity of natural selection, can drive evolutionary innovation and presage lineage extinction or loss of function. Mechanisms through which selection can be relaxed range from the removal of an existing selective constraint to a reduction in effective population size. Standard methods for estimating the strength and extent of purifying or positive selection from molecular sequence data are not suitable for detecting relaxed selection, because they lack power and can mistake an increase in the intensity of positive selection for relaxation of both purifying and positive selection. Here, we present a general hypothesis testing framework (RELAX) for detecting relaxed selection in a codon-based phylogenetic framework. Given two subsets of branches in a phylogeny, RELAX can determine whether selective strength was relaxed or intensified in one of these subsets relative to the other. We establish the validity of our test via simulations and show that it can distinguish between increased positive selection and a relaxation of selective strength. We also demonstrate the power of RELAX in a variety of biological scenarios where relaxation of selection has been hypothesized or demonstrated previously. We find that obligate and facultative γ-proteobacteria endosymbionts of insects are under relaxed selection compared with their free-living relatives and obligate endosymbionts are under relaxed selection compared with facultative endosymbionts. Selective strength is also relaxed in asexual Daphnia pulex lineages, compared with sexual lineages. Endogenous, nonfunctional, bornavirus-like elements are found to be under relaxed selection compared with exogenous Borna viruses. Finally, selection on the short-wavelength sensitive, SWS1, opsin genes in echolocating and nonecholocating bats is relaxed only in lineages in which this gene underwent pseudogenization; however, selection on the functional medium/long-wavelength sensitive opsin, M/LWS1, is found to be relaxed in all echolocating bats compared with nonecholocating bats.

AI snapshot

Written from the paper’s full text.

TL;DR

RELAX is a codon-based phylogenetic hypothesis test that detects relaxed or intensified selection by comparing the distribution of nonsynonymous/synonymous rate ratios (ω) between test and reference branch sets. It distinguishes relaxation from intensified positive selection and detects relaxed selection in endosymbionts, bat opsins, endogenous bornaviruses, and asexual *Daphnia pulex*.

Key findings

  • RELAX controls type I error at the nominal level: at *P* = 0.05, 52/1,000 simulated replicates rejected the null, giving an estimated type I error rate of 0.052 (95% CI: 0.040–0.068).
  • In power simulations, RELAX had 100% power to detect relaxed selection at *k* = 0.5 across three ω-distribution scenarios.
  • A single-ω-per-branch-partition approach misidentified intensified positive selection as relaxed selection in 100/100 replicates, whereas RELAX correctly identified intensified selection in 100/100 replicates with mean *k* = 3.79.
  • Empirical tests found relaxed selection in γ-proteobacteria endosymbionts vs. free-living relatives (*k* = 0.30, *P* < 0.0001), primary vs. secondary endosymbionts (*k* = 0.56, *P* < 0.0001), SWS1 pseudogenes in bats (*k* = 0.16, *P* < 0.0001), M/LWS1 in echolocating vs. tree/cave-roosting bats (*k* = 0.21, *P* = 0.0005), endogenous bornavirus elements (*k* = 0.02, *P* < 0.0001), and asexual *D. pulex* mitochondria (*k* = 0.63, *P* = 0.00003).

Methods

  • RELAX extends random effects branch-site (BS-REL) models: branches are partitioned a priori into test (T), reference (R), and optionally unclassified (U) sets; a selection intensity parameter *k* transforms the reference ω distribution to the test distribution, and a likelihood ratio test compares *k* = 1 (null) with a free *k* (alternative) using a 1-df χ² asymptotic distribution.
  • Simulations assessed type I error (1,000 codon alignments over the γ-proteobacteria phylogeny), power (three ω-distribution scenarios, 11 *k* values from 0 to 1, 100 replicates each), and robustness to increased positive selection (100 replicates over a four-taxon phylogeny).
  • Empirical data sets included 69 single-copy orthologs from 50 γ-proteobacteria species, bat SWS1 and M/LWS1 opsin genes, bornavirus nucleoprotein/EBLN-2 elements, and mitochondrial protein-coding regions from 27 *D. pulex* strains plus two *D. melanica* outgroups; alignments used MUSCLE and phylogenies used GTR + Γ4 in PhyML. RELAX is implemented in HyPhy v2.2 and Datamonkey.

Limitations

  • Stated: phylogenetic methods are uninformative about subbranch timescales; if episodic positive selection switches rapidly along a branch, the increased average ω can be indistinguishable from relaxation when ω remains below 1, and this can only be circumvented by sampling more sequences.
  • Stated: point estimates of *k* and ω lack reported sampling variances; high ω values are not reliably estimable due to flat likelihood surfaces; the General Descriptive model should not be used to select partitions for formal RELAX hypothesis testing because it would use the data twice and bias the test.
  • Evident from the text: empirical analyses rely on previously inferred phylogenies and alignments, and the *D. pulex* result conflicts with the recent-transition hypothesis of Tucker et al. (2013), a discrepancy the paper states it cannot resolve.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1093/molbev/msu400

Cites

75 of 75 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Claims the paper makes

31 of 53 cited sources checked · 19 supported · 0 flagged · accuracy not yet measured

  • Premisecitation

    The partitioned lineages test can be misleading because estimates of ! can be increased by processes such as intensification of positive selection when selective pressures vary both along branches and sites.

    “However, this test can be misleading because estimates of ! can be increased by processes such as intensification of positive selection when selective pressures vary both along branches and sites (Murrell et al. 2012).”
    Quote in the paperSubstantial claim1 of 1 sources checked
    • cites Murrell et al. 2012 source

      The Introduction (page unavailable) says tests that assume selection is constant across lineages can miss episodic positive selection, and the Discussion (page unavailable) describes positive selection on some branches being masked by purifying selection on others. These passages support the general warning that varying selective pressures can mislead a test, but they do not say that the partitioned lineages test yields increased ω estimates because positive selection intensifies. (LLM, conf=0.82)

      Experimental — accuracy not yet measured
  • Premisecitation

    Approaches estimating the number of sites under strong purifying selection in one region versus another and directly testing for differences in site-specific ! estimates are relevant primarily for comparing selective pressures across regions within a gene.

    “Relaxed selection across sites in a gene has also been explored by estimating the number of sites under strong purifying selection in one region versus another (Wertheim and Worobey 2009) and by directly testing for differences in sitespecific !estimates (Pond et al. 2006); however, such approaches are relevant primarily for comparing selective pressures across regions within a gene.”
    Quote in the paperSubstantial claim2 of 2 sources checked
    • cites Wertheim and Worobey 2009 source
    • cites Pond et al. 2006 source
  • Premisecitation

    RELAX builds upon the random effects branch-site (BS-REL) models, which were developed to overcome statistical and methodological restrictions imposed by the popular branch-site class of evolutionary models.

    “It builds upon the random effects branch-site (BS-REL) models (Kosakovsky Pond et al. 2011), which were developed to overcome statistical and methodological restrictions imposed by the popular”
    Quote in the paperSubstantial claim2 of 2 sources checked
    • cites Kosakovsky Pond et al. 2011 source
    • cites Yang and Nielsen 2002 source

      The Introduction (page not available) describes this paper’s branch-site models, and the Theory section states assumptions they impose. The copy does not address RELAX or BS-REL, so it cannot establish that RELAX builds on BS-REL or that BS-REL was developed to overcome restrictions of these models. (LLM, conf=0.78)

      Experimental — accuracy not yet measured
  • Premisecitation

    In a BS-REL framework, the ! ratio along each branch at a given site is drawn from a discrete distribution of ! values independently of other branches in the tree.

    “In a BS-REL framework, the !ratio along each branch at a given site is drawn from a discrete distribution of !values independently of other branches in the tree, affording RELAX a flexible and powerful mechanism for handling heterogeneous selective pressures.”
    Quote in the paperSubstantial claim1 of 1 sources checked
    • cites Kosakovsky Pond et al. 2011 source
  • Premisecitation

    The sum over !-configurations can be computed efficiently using Felsenstein's pruning algorithm.

    “We have previously shown (Kosakovsky Pond et al. 2011) that this sum can be computed efficiently using Felsenstein's pruning algorithm (Felsenstein 1981) with the transition probability matrix operating on branch b (with branch length t b ) given by a weighted sum of the three transition probability matrices:”
    Quote in the paperSubstantial claim1 of 2 sources checked
    • cites Kosakovsky Pond et al. 2011 source
    • cites Felsenstein 1981 source
  • May combine two claims. 31
  • Classified as probably not a substantial claim. 34
  • May combine two claims.

    Premisecitation

    Standard methods for estimating the strength and extent of purifying or positive selection from molecular sequence data are not suitable for detecting relaxed selection because they lack power and can mistake an increase in the intensity of positive selection for relaxation of both purifying and positive selection.

    “Standard methods for estimating the strength and extent of purifying or positive selection from molecular sequence data are not suitable for detecting relaxed selection, because they lack power and can mistake an increase in the intensity of positive selection for relaxation of both purifying and positive selection.”
    • Also stated in: § Introduction · plain
      “However, this test can be misleading because estimates of ! can be increased by processes such as intensification of positive selection when selective pressures vary both along branches and sites (Murrell et al. 2012).”
    Quote in the paperSubstantial claim1 of 1 sources checked
    • cites Murrell et al. 2012 source

      Split into 2 atoms — 0 confirmed, 2 partially supported, 0 disputed, 0 misattributed, 0 unverifiable. Expand the claim to see per-atom verdicts.

      Experimental — accuracy not yet measured
  • Classified as probably not a substantial claim.

    Premisecitation

    The most obvious approach to look for evidence of relaxed selection has been to identify loss of trait or pseudogenization.

    “The most obvious approach has been to identify loss of trait or pseudogenization (e.g., Dainat et al. [2012]).”
    Quote in the paperSubstantial claim1 of 1 sources checked
    • cites Dainat et al. [2012] source
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisecitation

    Relaxed selection has been implicated as both a driver of speciation and a precursor to extinction.

    “Relaxed selection has also been implicated as both a driver of speciation (Templeton 2008) and a precursor to extinction (Moran 1996).”
    Quote in the paperNot a substantial claim1 of 2 sources checked
    • cites Templeton 2008 source
    • cites Moran 1996 source
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Variation in the ! distribution across the phylogeny can be studied using the BS-REL model.

    “Such variation across the phylogeny can be studied using the BS-REL model (Kosakovsky Pond et al. 2011).”
    Quote in the paperNot a substantial claim1 of 1 sources checked
    • cites Kosakovsky Pond et al. 2011 source
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    The introduction of branch-site models addressed evolutionary rate variation not only over sites but also over branches.

    “The introduction of branch-site models addressed evolutionary rate variation not only over sites but also over branches (Yang and Nielsen 2002).”
    Quote in the paperNot a substantial claim1 of 1 sources checked
    • cites Yang and Nielsen 2002 source
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisecitation

    Random effects and mixed effects branch-site models allow detection of individual branches on which only a subset of sites are under positive selection or, conversely, individual sites which are under positive selection on only a subset of branches.

    “Our random effects (Kosakovsky Pond et al. 2011) and mixed effects (Murrell et al. 2012) branch-site models allow detection of individual branches on which only a subset of sites are under positive selection or, conversely, individual sites which are under positive selection on only a subset of branches.”
    Quote in the paperNot a substantial claim2 of 2 sources checked
    • cites Kosakovsky Pond et al. 2011 source

      Split into 2 atoms — 1 confirmed, 1 partially supported, 0 disputed, 0 misattributed, 0 unverifiable. Expand the claim to see per-atom verdicts.

      Experimental — accuracy not yet measured
    • cites Murrell et al. 2012 source
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    An ! distribution can be thought of as an evolutionary fingerprint characterizing the evolutionary pattern of the gene, where a small subset of sites can have very large ! values even when the majority of sites evolved under strong purifying selection.

    “An !distribution can be thought of as an evolutionary fingerprint characterizing the evolutionary pattern of the gene, where a small subset of sites can have very large !values even when the majority of sites evolved under strong purifying selection (Kosakovsky Pond, Scheffler, et al. 2010).”
    Quote in the paperNot a substantial claim1 of 1 sources checked
    • cites Kosakovsky Pond, Scheffler, et al. 2010 source

      Split into 2 atoms — 0 confirmed, 0 partially supported, 0 disputed, 0 misattributed, 2 unverifiable. Expand the claim to see per-atom verdicts.

      Experimental — accuracy not yet measured
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    The introduction of site-to-site rate variation allowed detection of individual sites evolving under positive selection by classifying sites into different rate categories and estimating a separate ! value specifically for the sites in the positive selection category.

    “The introduction of site-to-site rate variation allowed detection of individual sites evolving under positive selection, by classifying sites into different rate categories and estimating a separate !value specifically for the sites in the positive selection category (Nielsen and Yang 1998).”
    Quote in the paperNot a substantial claim1 of 1 sources checked
    • cites Nielsen and Yang 1998 source
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    The first codon models estimated parameters representing an averaged nonsynonymous-synonymous rate ratio.

    “The first codon models estimated parameters representing an averaged nonsynonymous-synonymous rate ratio.”
    Quote in the paperNot a substantial claim1 of 2 sources checked
    • cites Goldman and Yang 1994 source
    • cites Muse and Gaut 1994 source

      The Probabilistic Model for Codon Evolution section says the model’s parameters determine synonymous and nonsynonymous substitution rates. But Parameter Interpretation describes the estimated parameters as products of instantaneous rates and time, not as an averaged nonsynonymous-to-synonymous rate ratio; its averaging over codons concerns expected substitution counts. (LLM, conf=0.82)

      Experimental — accuracy not yet measured
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    It became conventional to think in terms of the distribution of ! values over sites instead of a single gene-wide !.

    “Therefore, it became conventional to think in terms of the distribution of !values over sites instead of a single gene-wide !(Yang et al. 2000).”
    Quote in the paperNot a substantial claim1 of 1 sources checked
    • cites Yang et al. 2000 source

      The Theory section models site-to-site variation in the dN/dS ratio using a statistical distribution, and the Results compare those models with a one-ratio model. This supports the conceptual contrast in the claim, but the paper does not establish that thinking in distributions became conventional across the field. (LLM, conf=0.82)

      Experimental — accuracy not yet measured
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Branch-level resolution vastly increased the detectable instances of positive selection.

    “This development vastly increased the detectable instances of positive selection (Murrell et al. 2012).”
    Quote in the paperNot a substantial claim1 of 1 sources checked
    • cites Murrell et al. 2012 source
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisecitation

    Traces of relaxed selection have been documented in bacterial endosymbionts, including fixation of deleterious alleles and overall elevated rates of nonsynonymous and radical amino acid replacements.

    “Accordingly, traces of relaxed selection have been documented in bacterial endosymbionts, including fixation of deleterious alleles (e.g., loss of DNA replication proofreading enzymes) and overall elevated rates of nonsynonymous and radical amino acid replacements (Wernegreen 2011).”
    Quote in the paperNot a substantial claim1 of 1 sources checked
    • cites Wernegreen 2011 source

      The Introduction (page unavailable) reports genome-wide rate acceleration concentrated at nonsynonymous sites, and the Discussion reports elevated radical replacements and evidence consistent with fixation of slightly deleterious mutations. The copy does not document the cited example of lost DNA-replication proofreading enzymes. The Discussion also says relaxed selection may contribute but is unlikely to explain the genome-wide pattern on its own; it favors genetic drift and linkage. (LLM, conf=0.84)

      Experimental — accuracy not yet measured
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Insects can also carry facultative (secondary) endosymbionts whose fitness effects on their hosts can be variable and complex.

    “Insects can also carry facultative (secondary) endosymbionts whose fitness effects on their hosts can be variable and complex (Oliver et al. 2006(Oliver et al. , 2009;;Degnan et al. 2009;Burke et al. 2010).”
    Quote in the paperNot a substantial claim2 of 3 sources checked
    • cites Oliver et al. 2006 source
    • cites Oliver et al. 2009 source
    • cites Degnan et al. 2009 source
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Many insect species have obligate (primary) bacterial endosymbionts that perform crucial metabolic processes for the host.

    “Many insect species have obligate (primary) bacterial endosymbionts that perform crucial metabolic processes for the host (Moran and Baumann 2000;McCutcheon et al. 2009;Vogel and Moran 2011).”
    Quote in the paperNot a substantial claim1 of 3 sources checked
    • cites Moran and Baumann 2000 source
    • cites McCutcheon et al. 2009 source
    • cites Vogel and Moran 2011 source
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Borna disease viruses may have integrated into mammalian genomes over 65 Ma.

    “Borna disease viruses may have integrated into mammalian genomes over 65 Ma (Horie et al. 2013).”
    Quote in the paperNot a substantial claim1 of 1 sources checked
    • cites Horie et al. 2013 source

      The ‘Paleovirology of bornaviruses’ section says that bornavirus-like elements in several mammalian genomes appear to have been established approximately 65 Ma, and that bornavirus infection in mammals may go back to 65 Ma. It does not establish integration over 65 Ma ago. (LLM, conf=0.85)

      Experimental — accuracy not yet measured
  • Classified as probably not a substantial claim.

    Premisecitation

    Tucker et al. (2013) examined the genomes of sexual and asexual D. pulex lineages and found evidence for more recent transitions to asexuality, on the order of decades.

    “Recently, Tucker et al. (2013) examined the genomes of sexual and asexual D. pulex lineages and found evidence for more recent transitions to asexuality, on the order of decades.”
    Quote in the paperSubstantial claim1 of 1 sources checked
    • cites Tucker et al. (2013) source
  • May combine two claims.

    Premisecitation

    Asexual animal species are rare and typically do not persist for long periods of time.

    “however, asexual animal species are rare and typically do not persist for long periods of time (Barton 2009).”
    Quote in the paperSubstantial claim1 of 1 sources checked
    • cites Barton 2009 source

      The paper states that asexual species have a “relatively short life,” supporting the persistence part of the claim (section and page not provided in our copy). Its introduction says sexual reproduction dominates eukaryotic biology, but does not specifically establish that asexual animal species are rare. (LLM, conf=0.83)

      Experimental — accuracy not yet measured
  • May combine two claims.

    Findingown data

    RELAX correctly identified intensified selection along test branches in 100/100 replicates, with a mean selection intensity parameter k of 3.79.

    “RELAX, by avoiding the unrealistic single-!assumption, correctly identified intensified selection along test branches in 100/100 replicates, with a mean selection intensity parameter (k) of 3.79.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    In endosymbionts, !T,1 representing the strongest purifying selection shifted up by an order of magnitude; !T,2 representing weakly purifying selection shifted from 0.3 to 1.0; and !T,3 representing positive selection shifted down by an order of magnitude.

    “Specifically, !T;1 , representing the strongest purifying selection, shifted up by an order of magnitude; !T;2 , representing weakly purifying selection shifted from 0.3 to 1.0; and !T;3 , representing positive selection, shifted down by an order of magnitude.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Relative to free-living bacteria, there was sufficient signal to detect relaxation in both primary endosymbionts (k = 0.28; P < 0.0001) and secondary endosymbionts (k = 0.61; P < 0.0001) independently.

    “Relative to free-living bacteria, there was sufficient signal to detect relaxation in both primary endosymbionts (k = 0.28; P-value < 0.0001) and secondary endosymbionts (k = 0.61; Pvalue < 0.0001) independently.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The backbone of the phylogeny is under intense selection (k_b > 1), as are terminal branches leading to free-living -proteobacteria, whereas the endosymbiont lineages tend to be more relaxed (k_b < 1).

    “The backbone of the phylogeny is under intense selection (k b 41 shown as red in fig.6), as are terminal branches leading to free-living -proteobacteria, whereas the endosymbiont lineages tend to be more relaxed (k b < 1 shown as blue/purple in fig.6).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    RELAX is well behaved on simulated data and has power to detect relaxed selection via simulation and in four biological systems predicted to exhibit changes in the intensity of selection on a subset of branches.

    “We demonstrated that RELAX is well behaved on simulated data and established its power to detect relaxed selection via simulation and in four biological systems predicted to exhibit changes in the intensity of selection on a subset of branches.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    We inferred 78% of the sites evolved under an omega of 0.82, with the remaining 22% evolving under omega = 1.2.

    “We inferred 78% of the sites evolved under an ! of 0.82, with the remaining 22% evolving under !¼ 1:2.”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim. May combine two claims.

    Premisereasoning

    A higher ! within a group is interpreted as evidence for relaxation of selective constraints, as those constraints tend to reduce the rate of nonsynonymous substitutions relative to that of synonymous substitutions, thereby keeping ! values low.

    “A higher !within a group is interpreted as evidence for relaxation of selective constraints, as those constraints tend to reduce the rate of nonsynonymous substitutions relative to that of synonymous substitutions, thereby keeping !values low.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premisenone given

    Identifying loss of trait or pseudogenization detects the terminal effects of relaxation but is not capable of detecting subtler changes in the intensity of selection that may not drive a gene to obsolescence or cause a loss of function.

    “This approach detects the terminal effects of relaxation, but it is not capable of detecting subtler changes in the intensity of selection that may not drive a gene to obsolescence or cause a loss of function.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Study design and conductown data

    RELAX's validity is established via simulations, and it can distinguish between increased positive selection and a relaxation of selective strength.

    “We establish the validity of our test via simulations and show that it can distinguish between increased positive selection and a relaxation of selective strength.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S4 (based).

    Premisenone given

    When rate variation over time is not modeled, ! values for specific sites or site categories represent time-averages over all branches in the phylogeny.

    “When rate variation over time is not modeled, !values for specific sites or site categories represent time-averages over all branches in the phylogeny.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Premisereasoning

    In branch-site models with multiple ! categories represented by different proportions of sites, relaxation can cause the ! values inferred for the selection categories to move toward 1 and/or the proportions of sites belonging to the different categories to change in such a way that more sites are assigned to categories with ! values closer to 1.

    “In the context of branch-site models with multiple !categories represented by different proportions of sites, this trend can have two different effects: The ! values inferred for the selection categories can move toward 1 (fig.1) and/or the proportions of sites belonging to the different categories can Detecting Relaxed Selection .doi:10.1093/molbev/msu400MBE change in such a way that more sites”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisereasoning

    Gene-wide average ! values were typically below 1 because most sites in protein-coding genes are strongly conserved.

    “Such gene-wide average !values were typically below 1, because most sites in protein-coding genes are strongly conserved.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisereasoning

    Site-to-site rate variation led to the detection of sites for which ! is well above 1, even in genes where the average ! value estimated across all sites is below 1.

    “This approach led to the detection of sites for which ! is well above 1, even in genes where the average !value estimated across all sites is below 1.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Interpretationreasoning

    When the introduction of k produces a statistically significant improvement in model fit, selection on test branches is intensified (k > 1) or relaxed (k < 1) compared with background branches.

    “When the introduction of k produces a statistically significant improvement in model fit, we conclude that selection on test branches is intensified (k 4 1) or relaxed (k < 1) compared with background branches.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisereasoning

    The change from phylogeny-level resolution to branch-level resolution yielded higher branch-specific ! values than phylogeny-wide ! values.

    “The change from phylogeny-level resolution to branch-level resolution yielded higher branch-specific !values than phylogeny-wide !values, mirroring the effect obtained when site-to-site rate variation was introduced.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Premisenone given

    When selection is relaxed, ! values smaller than 1 increase toward 1, whereas ! values above 1 decrease.

    “When selection is relaxed, the smaller !values increase toward 1, whereas !values above 1 decrease.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Study design and conductown data

    RELAX's performance is validated using simulations as well as four biological examples where relaxation would be expected according to evolutionary theory.

    “We validate the performance of RELAX using simulations as well as four biological examples where relaxation would be expected according to evolutionary theory: 1) endosymbiosis in -proteobacteria, 2) color vision genes in bats, 3) viral endogenization in bornaviruses, and 4) asexual reproduction in crustaceans.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationreasoning

    If the null model is rejected by the LRT, selection on the test branches is intensified (k > 1) or relaxed (k < 1) compared with the reference branches.

    “If the null model is rejected by the LRT, it indicates that selection on the test branches is intensified (k 4 1) or relaxed (k < 1) compared with the reference branches.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S2 (checkable), S3 (load-bearing), S4 (based).

    Interpretationnone given

    The simulated scenarios capture only a tiny fraction of the breadth under which selection may actually operate.

    “Importantly, these simulated scenarios capture only a tiny fraction of the breadth under which selection may actually operate.”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim.

    Study design and conductnone given

    Three sets of simulations were performed to assess the validity of RELAX.

    “Three sets of simulations were performed to assess the validity of RELAX.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing), S4 (based).

    Interpretationnone given

    In empirical data sets, the signal for relaxed selection is expected to be driven primarily by sites evolving under purifying rather than positive selection, which are by far the majority.

    “In empirical data sets, we expect the signal for relaxed selection to be driven primarily by sites evolving under purifying rather than positive selection, which are by far the majority.”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim.

    Study design and conductown data

    The robustness of RELAX to the effect of increasing only the selection intensity of positive selection in a subset of branches (model misspecification) was investigated.

    “We also investigated the robustness of RELAX to the effect of increasing only the selection intensity of positive selection in a subset of branches (i.e., model misspecification).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premiseown data

    High ! values have a flat likelihood surface, making larger empirical estimates for this parameter unreliable.

    “We enforced this constraint due to the flat likelihood surface associated with high !values, which makes the larger empirical estimate for this parameter unreliable.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisecitation

    Natural selection is less efficient in this setting, and Muller's Ratchet can lead to the accumulation of slightly deleterious mutations.

    “Natural selection is less efficient in this setting, and Muller's Ratchet can lead to the accumulation of slightly deleterious mutations (Muller 1964).”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites Muller 1964 source
  • May combine two claims.

    Premisecitation

    Obligate endosymbionts have arisen multiple times among the -proteobacteria and are passed vertically from mother to child, thereby introducing convergent features in the bacterial genomes: gene loss, smaller genomes, lower GC content, and reduced thermal stability of proteins.

    “These obligate endosymbionts have arisen multiple times among the -proteobacteria (Husn ık et al. 2011) and are passed vertically from mother to child, thereby introducing convergent features in the bacterial genomes: gene loss, smaller genomes, lower GC content, and reduced thermal stability of proteins (Wernegreen 2002).”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites Husn ık et al. 2011

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites Wernegreen 2002 source
  • Classified as probably not a substantial claim. May combine two claims.

    Interpretationreasoning

    Using the General Descriptive model to select partitions for formal RELAX hypothesis testing would amount to using the data twice and substantially bias the test result.

    “Furthermore, the General Descriptive model should not be used to select partitions for formal RELAX hypothesis testing, because it would amount to using the data twice and substantially bias the test result.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premisecitation

    Phylogenetic inference that does not account for variable substitution rates arising from skewed GC content found in endosymbionts tends to force their lineages together in a phylogenetic tree, whereas nonhomogeneous models of molecular evolution suggest a larger number of independent origins of endosymbionts from free-living bacteria.

    “Phylogenetic inference that does not account for variable substitution rates arising from skewed GC content found in endosymbionts tends to force their lineages together in a phylogenetic tree, whereas nonhomogeneous models of molecular evolution suggest a larger number of independent origins of endosymbionts from free-living bacteria (Husn ık et al. 2011).”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites Husn ık et al. 2011

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Interpretationcitation

    The pattern of stronger relaxation in primary than secondary endosymbionts was expected because primary endosymbionts likely experience tighter bottlenecks during vertical transmission and their genomes show greater evidence of extreme evolutionary features (e.g., gene loss and AT bias).

    “This pattern was expected since primary endosymbionts likely experience tighter bottlenecks during vertical transmission and their genomes show greater evidence of extreme evolutionary features (e.g., gene loss and AT bias; McCutcheon and Moran 2012).”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites McCutcheon and Moran 2012 source
  • May combine two claims.

    Study design and conductreasoning

    RELAX determines the strength of evidence that k is different from 1 (reported as P-values) but not the sampling variances of the estimates.

    “For all of the point estimates of k and ! reported below, it should be borne in mind that our statistical methodology allows us to determine the strength of evidence that k is different from 1 (reported as P-values) but not the sampling variances of the estimates.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisecitation

    Zhao, Rossiter, et al. (2009) inferred lineage-specific omega values and consistent patterns of purifying selection (omega ranging between 0.0 and 0.48).

    “Zhao, Rossiter, et al. (2009) inferred lineage-specific !values and consistent patterns of purifying selection (! ranging between 0.0 and 0.48).”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites Zhao, Rossiter, et al. (2009)

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    The loss of vision in animals that reside in total darkness is a classic case of trait loss caused by relaxed selection.

    “The loss of vision in animals that reside in total darkness (e.g., cave-dwelling fish) is a classic case of trait loss caused by relaxed selection (Darwin 1859).”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites Darwin 1859

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Interpretationown data

    The relaxation of selection in all echolocating species over tens of millions of years has not led to the formation of M/LWS1 pseudogenes in echolocating species, suggesting some functionality of this gene may have been retained.

    “This relaxation of selection in all echolocating species, over tens of millions of years, has not led to the formation of M/LWS1 pseudogenes in the echolocating species, which suggests some functionality of this gene may have been retained.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Interpretationreasoning

    One would expect an excess of mildly deleterious mutations accumulating in mitochondrial genomes due to relaxed selection in asexual lineages.

    “Therefore, one would expect an excess of mildly deleterious mutations accumulating in mitochondrial genomes due to relaxed selection in asexual lineages.”
    Quote in the paperNot a substantial claim
    • Rests on an assumption: After the transition to asexuality, there is no longer segregation between the nuclear and mitochondrial genomes, and the mitochondrial and nuclear genomes are in complete linkage disequilibrium.
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    This bias can be corrected by using nonhomogeneous nucleotide substitution models in phylogenetic inference.

    “This bias can be corrected by using nonhomogeneous nucleotide substitution models in phylogenetic inference (Galtier andGouy 1995, 1998).”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites Galtier andGouy 1995, 1998

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Probably not a substantial claim: may not pass S4 (based). May combine two claims.

    Study design and conductnone given

    A caveat of phylogenetic methods in general is that sequence data are uninformative regarding changes that happen at subbranch timescales and consequently that changes in selective pressure over short timescales are difficult to characterize.

    “A caveat of phylogenetic methods in general is that sequence data are uninformative regarding changes that happen at subbranch timescales and consequently that changes in selective pressure over short timescales are difficult to characterize.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Premisecitation

    Mitochondria are not part of the sexually reproducing genome, but they should experience the selective effects of a transition from sexuality to asexuality.

    “Mitochondria are not part of the sexually reproducing genome, but they should experience the selective effects of a transition from sexuality to asexuality (Paland and Lynch 2006).”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites Paland and Lynch 2006 source
  • Classified as probably not a substantial claim.

    Premisecitation

    Sexual reproduction is costly from an evolutionary perspective.

    “From an evolutionary perspective, sexual reproduction is costly (Maynard Smith 1978);”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites Maynard Smith 1978

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Interpretationreasoning

    These selective constraints likely arise from the necessity to encode genes for complex biological functions in compact genomes.

    “These selective constraints likely arise from the necessity to encode genes for complex biological functions in compact genomes.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Interpretationreasoning

    The small deviation from neutrality may be driven by the components of these branches preceding viral endogenization.

    “The small deviation from neutrality may be driven by the components of these branches preceding viral endogenization.”
    Quote in the paperNot a substantial claim
    • Rests on an assumption: If the endogenous viral elements were not subject to any level of purifying selection at the amino acid level, one would expect all the sites to evolve under an omega equal to 1.
56 of 56
  • The authors state this as an assumption.

    Site likelihoods are independent.

    “under the standard assumption of independence”
  • The authors state this as an assumption.

    Nucleotide substitution rates follow the general time reversible (GTR) form.

    “assumed to follow the general time reversible [GTR] form”
  • The paper relies on this without support.

    The phylogeny used in the analyses is accurate.

    “when provided with an accurate phylogeny”
  • The paper relies on this without support.

    Once a gene has lost functionality, purifying selection is expected to abate.

    “Once a gene has lost functionality, purifying selection is expected to abate”
  • The paper relies on this without support.

    The SWS1 gene has been lost in cave-roosting fruit bats and HDC echolocating species, but functional SWS1 genes have been maintained in tree-roosting and LDC echolocating species.

    “The SWS1 gene has been lost in caveroosting fruit bats and HDC echolocating species, but functional SWS1 genes have been maintained in tree-roosting and LDC echolocating species.”
  • The paper relies on this without support.

    The M/LWS1 gene has been retained across all documented bat species.

    “A second opsin gene involved in green-red vision, medium/long wavelength sensitive (M/LWS1), has been retained across all documented bat species.”
  • The paper relies on this without support.

    The Borna disease virus has experienced more recent endogenization events, such as integration into the 13-lined ground squirrel genome.

    “The virus has also experienced more recent endogenization events, such as the integration into the 13-lined ground squirrel genome.”
  • The paper relies on this without support.

    After the transition to asexuality, there is no longer segregation between the nuclear and mitochondrial genomes, and the mitochondrial and nuclear genomes are in complete linkage disequilibrium.

    “After the transition to asexuality, there is no longer segregation between the nuclear and mitochondrial genomes, and the mitochondrial and nuclear genomes are in complete linkage disequilibrium.”
  • The paper relies on this without support.

    If the endogenous viral elements were not subject to any level of purifying selection at the amino acid level, one would expect all the sites to evolve under an omega equal to 1.

    “If the endogenous viral elements were not subject to any level of purifying selection at the amino acid level, one would expect all the sites to evolve under an !equal to 1.”
9 of 9

Cites

Flagged first, then by citations

2 of the 77 references couldn’t be matched to a known paper.

Research works printed but not linked: 1

  1. Maynard Smith, 1978, The evolution of sex

Cited by

Sorted by citations

Authors
Wertheim, J. O., Murrell, B., Smith, M. D., Kosakovsky Pond, S. L., & Scheffler, K.
Affiliations
  • University of California San Diego (US)
  • Stellenbosch University (ZA)
Corresponding author
  • Wertheim, Joel O.
  • Scheffler, Konrad
Venue
Molecular Biology and Evolution
Published
2014-12-23
Volume
32
Issue
3
Pages
820-832
ISSN
  • 0737-4038 (ISSN-L)
  • 1537-1719
OpenAlex
W2028098940
Cited by
778
References
77
Citations by year
2026: 113 · 2025: 85 · 2024: 119 · 2023: 86 · 2022: 101 · 2021: 79 · 2020: 56 · 2019: 54 · 2018: 32 · 2017: 30 · 2016: 18 · 2015: 5
FWCI
11.57
Type
journal article
Language
en
Open access
Locations
  • Molecular Biology and Evolution (publishedVersion)
  • PubMed Central (publishedVersion)
  • PubMed (publishedVersion)
  • eScholarship (California Digital Library) (submittedVersion)
Funders
  • International AIDS Vaccine Initiative
  • National Research Foundation
  • National Institutes of Health (AI090970, AI-36214)
  • Center for AIDS Research, University of Washington
  • University of California, San Diego
  • National Heart, Lung, and Blood Institute (U54HL108460)
  • National Institute on Drug Abuse (DA-034978, DP1DA034978)
  • National Institute of General Medical Sciences (GM093939, R01GM093939, U01GM110749)
  • National Institute of Allergy and Infectious Diseases (AI090970, AI110181, AI-100665, K24AI100665, K01AI110181, P30AI036214, U19AI090970)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Evolution and Genetic Dynamics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA and protein synthesis mechanisms
Keywords
relaxed selection, purifying selection, positive selection, endosymbionts, Daphnia pulex, asexual lineages, opsin genes, pseudogenization, echolocating bats
MeSH terms
Animals, Chiroptera, Codon, Computer Simulation, Daphnia, Models, Genetic, Phylogeny, Selection, Genetic, Sequence Analysis, DNA, Evolution, Molecular, Computational Biology, Gammaproteobacteria, Bornaviridae, Opsins

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Claims read by deepseek-v4.1-flash under definition version 3 (extract-v3) on 2026-11-02.

Classified by jev-1.13-20260917 (gold-guideline-v0.2) on 2026-11-02.

Cited sources are checked by gpt-6-sol.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review