scholarly journals On the interpretation of results returned by branch-site models

Author(s):  
Stephane Guindon

In a recent study, Murrell et al. (2015) compared the performance of several branch-site models of codon evolution. Their interpretation of results published by Lu & Guindon (2014) suggests that the stochastic branch-site model implemented in the software fitmodel is anti-conservative altogether, i.e., positive selection is detected more often than expected when analyzing sequences evolving under a mixture of neutrality and negative selection. I argue here that this presentation of the performance of fitmodel is misleading and should not deter evolutionary biologists from using this approach in exploratory analyses of selection patterns at the molecular level.

2015 ◽  
Author(s):  
Stephane Guindon

In a recent study, Murrell et al. (2015) compared the performance of several branch-site models of codon evolution. Their interpretation of results published by Lu & Guindon (2014) suggests that the stochastic branch-site model implemented in the software fitmodel is anti-conservative altogether, i.e., positive selection is detected more often than expected when analyzing sequences evolving under a mixture of neutrality and negative selection. I argue here that this presentation of the performance of fitmodel is misleading and should not deter evolutionary biologists from using this approach in exploratory analyses of selection patterns at the molecular level.


2019 ◽  
Vol 69 (4) ◽  
pp. 722-738 ◽  
Author(s):  
Christopher T Jones ◽  
Noor Youssef ◽  
Edward Susko ◽  
Joseph P Bielawski

Abstract A central objective in biology is to link adaptive evolution in a gene to structural and/or functional phenotypic novelties. Yet most analytic methods make inferences mainly from either phenotypic data or genetic data alone. A small number of models have been developed to infer correlations between the rate of molecular evolution and changes in a discrete or continuous life history trait. But such correlations are not necessarily evidence of adaptation. Here, we present a novel approach called the phenotype–genotype branch-site model (PG-BSM) designed to detect evidence of adaptive codon evolution associated with discrete-state phenotype evolution. An episode of adaptation is inferred under standard codon substitution models when there is evidence of positive selection in the form of an elevation in the nonsynonymous-to-synonymous rate ratio $\omega$ to a value $\omega > 1$. As it is becoming increasingly clear that $\omega > 1$ can occur without adaptation, the PG-BSM was formulated to infer an instance of adaptive evolution without appealing to evidence of positive selection. The null model makes use of a covarion-like component to account for general heterotachy (i.e., random changes in the evolutionary rate at a site over time). The alternative model employs samples of the phenotypic evolutionary history to test for phenomenological patterns of heterotachy consistent with specific mechanisms of molecular adaptation. These include 1) a persistent increase/decrease in $\omega$ at a site following a change in phenotype (the pattern) consistent with an increase/decrease in the functional importance of the site (the mechanism); and 2) a transient increase in $\omega$ at a site along a branch over which the phenotype changed (the pattern) consistent with a change in the site’s optimal amino acid (the mechanism). Rejection of the null is followed by post hoc analyses to identify sites with strongest evidence for adaptation in association with changes in the phenotype as well as the most likely evolutionary history of the phenotype. Simulation studies based on a novel method for generating mechanistically realistic signatures of molecular adaptation show that the PG-BSM has good statistical properties. Analyses of real alignments show that site patterns identified post hoc are consistent with the specific mechanisms of adaptation included in the alternate model. Further simulation studies show that the covarion-like component of the PG-BSM plays a crucial role in mitigating recently discovered statistical pathologies associated with confounding by accounting for heterotachy-by-any-cause. [Adaptive evolution; branch-site model; confounding; mutation-selection; phenotype–genotype.]


2021 ◽  
Vol 12 ◽  
Author(s):  
Vladimir M. Jovanovic ◽  
Melanie Sarfert ◽  
Carlos S. Reyna-Blanco ◽  
Henrike Indrischek ◽  
Dulce I. Valdivia ◽  
...  

Gene regulatory factors (GRFs), such as transcription factors, co-factors and histone-modifying enzymes, play many important roles in modifying gene expression in biological processes. They have also been proposed to underlie speciation and adaptation. To investigate potential contributions of GRFs to primate evolution, we analyzed GRF genes in 27 publicly available primate genomes. Genes coding for zinc finger (ZNF) proteins, especially ZNFs with a Krüppel-associated box (KRAB) domain were the most abundant TFs in all genomes. Gene numbers per TF family differed between all species. To detect signs of positive selection in GRF genes we investigated more than 3,000 human GRFs with their more than 70,000 orthologs in 26 non-human primates. We implemented two independent tests for positive selection, the branch-site-model of the PAML suite and aBSREL of the HyPhy suite, focusing on the human and great ape branch. Our workflow included rigorous procedures to reduce the number of false positives: excluding distantly similar orthologs, manual corrections of alignments, and considering only genes and sites detected by both tests for positive selection. Furthermore, we verified the candidate sites for selection by investigating their variation within human and non-human great ape population data. In order to approximately assign a date to positively selected sites in the human lineage, we analyzed archaic human genomes. Our work revealed with high confidence five GRFs that have been positively selected on the human lineage and one GRF that has been positively selected on the great ape lineage. These GRFs are scattered on different chromosomes and have been previously linked to diverse functions. For some of them a role in speciation and/or adaptation can be proposed based on the expression pattern or association with human diseases, but it seems that they all contributed independently to human evolution. Four of the positively selected GRFs are KRAB-ZNF proteins, that induce changes in target genes co-expression and/or through arms race with transposable elements. Since each positively selected GRF contains several sites with evidence for positive selection, we suggest that these GRFs participated pleiotropically to phenotypic adaptations in humans.


Author(s):  
Joshua H T Potter ◽  
Kalina T J Davies ◽  
Laurel R Yohe ◽  
Miluska K R Sanchez ◽  
Edgardo M Rengifo ◽  
...  

Abstract Dietary adaptation is a major feature of phenotypic and ecological diversification, yet the genetic basis of dietary shifts is poorly understood. Among mammals, Neotropical leaf-nosed bats (family Phyllostomidae) show unmatched diversity in diet; from a putative insectivorous ancestor, phyllostomids have radiated to specialize on diverse food sources, including blood, nectar, and fruit. To assess whether dietary diversification in this group was accompanied by molecular adaptations for changing metabolic demands, we sequenced 89 transcriptomes across 58 species, and combined these with published data to compare ∼13,000 protein coding genes across 66 species. We tested for positive selection on focal lineages, including those inferred to have undergone dietary shifts. Unexpectedly, we found a broad signature of positive selection in the ancestral phyllostomid branch, spanning genes implicated in the metabolism of all major macronutrients, yet few positively selected genes at the inferred switch to plantivory. Branches corresponding to blood- and nectar-based diets showed selection in loci underpinning nitrogenous waste excretion and glycolysis, respectively. Intriguingly, patterns of selection in metabolism genes were mirrored by those in loci implicated in craniofacial remodelling, a trait previously linked to phyllostomid dietary specialisation. Finally, using simulations, we show that the widely-used branch-site model is likely to be misspecified, with the implication that it is too conservative and probably under-reports true cases of positive selection. Our findings point to a complex picture of adaptive radiation, in which the evolution of new dietary specialisations has been facilitated by early adaptations combined with the generation of new genetic variation.


2019 ◽  
Author(s):  
Qingjian Chen ◽  
Ao Lan ◽  
Xu Shen ◽  
Chung-I Wu

AbstractThe widely accepted view that evolution proceeds in small steps is based on two premises: i) negative selection acts strongly against large differences (Kimura 1983); and ii) positive selection favors small-step changes (Fisher 1930). The two premises are not biologically connected and should be evaluated separately. We now extend the approach of Tang et al. (2004) to codon evolution for the entire genome. Codon substitution rate is a function of the physico-chemical distance between amino acids (AAs), equated with the step size of evolution. This step size depends on a large number of physico-chemical properties as 46 of the 48 properties examined affect the rate. Between 9 pairs of closely-related species of plants, invertebrates and vertebrates, the evolutionary rate is indeed strongly andnegativelycorrelated with the AA distance (ΔU, scaled to [0, 1]). While the analyses corroborate the published results that relied on partial genomes, there is an important difference: ΔUis strongly correlated with the evolutionary rate (R2> 0.8) only when the genes are under predominant negative selection. Nevertheless, since most genes in most taxa are strongly constrained by negative selection, ΔUwould appear to be a nearly-universal measure of codon evolution. In conclusion, the driving force of the small-step evolution at the codon level is negative selection. The unanswered question of whether positive selection may, or may not, follow the small-step rule will be addressed in a companion study (Chen, et al. 2019).


2019 ◽  
Vol 11 (10) ◽  
pp. 2702-2712 ◽  
Author(s):  
Qingjian Chen ◽  
Ao Lan ◽  
Xu Shen ◽  
Chung-I Wu

Abstract The widely accepted view that evolution proceeds in small steps is based on two premises: 1) negative selection acts strongly against large differences and 2) positive selection favors small-step changes. The two premises are not biologically connected and should be evaluated separately. We now extend a previous approach to studying codon evolution in the entire genome. Codon substitution rate is a function of the physicochemical distance between amino acids (AAs), equated with the step size of evolution. Between nine pairs of closely related species of plants, invertebrates, and vertebrates, the evolutionary rate is strongly and negatively correlated with a set of AA distances (ΔU, scaled to [0, 1]). ΔU, a composite measure of evolutionary rates across diverse taxa, is influenced by almost all of the 48 physicochemical properties used here. The new analyses reveal a crucial trend hidden from previous studies: ΔU is strongly correlated with the evolutionary rate (R2 > 0.8) only when the genes are predominantly under negative selection. Because most genes in most taxa are strongly constrained by negative selection, ΔU has indeed appeared to be a nearly universal measure of codon evolution. In conclusion, molecular evolution at the codon level generally takes small steps due to the prevailing negative selection. Whether positive selection may, or may not, follow the small-step rule is addressed in a companion study.


Insects ◽  
2021 ◽  
Vol 12 (7) ◽  
pp. 656
Author(s):  
Xiao-Dong Xu ◽  
Jia-Yin Guan ◽  
Zi-Yi Zhang ◽  
Yu-Rou Cao ◽  
Yin-Yin Cai ◽  
...  

We determined 15 complete and two nearly complete mitogenomes of Heptageniidae belonging to three subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae) and six genera (Afronurus, Epeorus, Leucrocuta, Maccaffertium, Stenacron, and Stenonema). Species of Rhithrogeninae and Ecdyonurinae had the same gene rearrangement of CR-I-M-Q-M-ND2, whereas a novel gene rearrangement of CR-I-M-Q-NCR-ND2 was found in Heptageniinae. Non-coding regions (NCRs) of 25–47 bp located between trnA and trnR were observed in all mayflies of Heptageniidae, which may be a synapomorphy for Heptageniidae. Both the BI and ML phylogenetic analyses supported the monophyly of Heptageniidae and its subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae). The phylogenetic results combined with gene rearrangements and NCR locations confirmed the relationship of the subfamilies as (Heptageniinae + (Rhithrogeninae + Ecdyonurinae)). To assess the effects of low-temperature stress on Heptageniidae species from Ottawa, Canada, we found 27 positive selection sites in eight protein-coding genes (PCGs) using the branch-site model. The selection pressure analyses suggested that mitochondrial PCGs underwent positive selection to meet the energy requirements under low-temperature stress.


PeerJ ◽  
2017 ◽  
Vol 5 ◽  
pp. e3145 ◽  
Author(s):  
Edson Ishengoma ◽  
Morris Agaba ◽  
Douglas R. Cavener

BackgroundThe capacity of visually oriented species to perceive and respond to visual signal is integral to their evolutionary success. Giraffes are closely related to okapi, but the two species have broad range of phenotypic differences including their visual capacities. Vision studies rank giraffe’s visual acuity higher than all other artiodactyls despite sharing similar vision ecological determinants with many of them. The extent to which the giraffe’s unique visual capacity and its difference with okapi is reflected by changes in their vision genes is not understood.MethodsThe recent availability of giraffe and okapi genomes provided opportunity to identify giraffe and okapi vision genes. Multiple strategies were employed to identify thirty-six candidate mammalian vision genes in giraffe and okapi genomes. Quantification of selection pressure was performed by a combination of branch-site tests of positive selection and clade models of selection divergence through comparing giraffe and okapi vision genes and orthologous sequences from other mammals.ResultsSignatures of selection were identified in key genes that could potentially underlie giraffe and okapi visual adaptations. Importantly, some genes that contribute to optical transparency of the eye and those that are critical in light signaling pathway were found to show signatures of adaptive evolution or selection divergence. Comparison between giraffe and other ruminants identifies significant selection divergence inCRYAAandOPN1LW. Significant selection divergence was identified inSAGwhile positive selection was detected inLUMwhen okapi is compared with ruminants and other mammals. Sequence analysis ofOPN1LWshowed that at least one of the sites known to affect spectral sensitivity of the red pigment is uniquely divergent between giraffe and other ruminants.DiscussionBy taking a systemic approach to gene function in vision, the results provide the first molecular clues associated with giraffe and okapi vision adaptations. At least some of the genes that exhibit signature of selection may reflect adaptive response to differences in giraffe and okapi habitat. We hypothesize that requirement for long distance vision associated with predation and communication with conspecifics likely played an important role in the adaptive pressure on giraffe vision genes.


Sign in / Sign up

Export Citation Format

Share Document