The structures of cytochrome c and the rates of molecular evolution

Functional Optimization in Distinct Tissues and Conditions Constrains the Rate of Protein Evolution

Understanding the main determinants of protein evolution is a fundamental challenge in biology. Despite many decades of active research, the molecular and cellular mechanisms underlying the substantial variability of evolutionary rates across cellular proteins are not currently well understood. It also remains unclear how protein molecular function is optimized in the context of multicellular species and why many proteins, such as enzymes, are only moderately efficient on average. Our analysis of genomics and functional datasets reveals in multiple organisms a strong inverse relationship between the optimality of protein molecular function and the rate of protein evolution. Furthermore, we find that highly expressed proteins tend to be substantially more functionally optimized. These results suggest that cellular expression costs lead to more pronounced functional optimization of abundant proteins and that the purifying selection to maintain high levels of functional optimality significantly slows protein evolution. We observe that in multicellular species both the rate of protein evolution and the degree of protein functional efficiency are primarily affected by expression in several distinct cell types and tissues, specifically, in developed neurons with upregulated synaptic processes in animals and in young and fast-growing tissues in plants. Overall, our analysis reveals how various constraints from the molecular, cellular, and species' levels of biological organization jointly affect the rate of protein evolution and the level of protein functional adaptation.

Modeling Substitution Rate Evolution across Lineages and Relaxing the Molecular Clock

Relaxing the molecular clock using models of how substitution rates change across lineages has become essential for addressing evolutionary problems. The diversity of rate evolution models and their implementations are substantial, and studies have demonstrated their impact on divergence time estimates can be as significant as that of calibration information. In this review, we trace the development of rate evolution models from the proposal of the molecular clock concept to the development of sophisticated Bayesian and non-Bayesian methods that handle rate variation in phylogenies. We discuss the various approaches to modeling rate evolution, provide a comprehensive list of available software, and examine the challenges and advancements of the prevalent Bayesian framework, contrasting them to faster non-Bayesian methods. Lastly, we offer insights into potential advancements in the field in the era of big data.

The effect of cytochrome c on Na,K-ATPase

Na,K-ATPase is a crucial enzyme responsible for maintaining Na+, K+-gradients across the cell membrane, which is essential for numerous physiological processes within various organs and tissues. Due to its significance in cellular physiology, inhibiting Na,K-ATPase can have profound physiological consequences. This characteristic makes it a target for various pharmacological applications, and drugs that modulate the pump's activity are thus used in the treatment of several medical conditions. Cytochrome c (Cytc) is a protein with dual functions in the cell. In the mitochondria, it is essential for ATP synthesis and energy production. However, in response to apoptotic stimuli, it is released into the cytosol, where it triggers programmed cell death through the intrinsic apoptosis pathway. Aside from its role in canonical intrinsic apoptosis, Cytc also plays additional roles. For instance, Cytc participates in certain non-apoptotic functions -those which are less well-understood in comparison to its role in apoptosis. Within this in vitro study, we have shown the impact of Cytc on Na,K-ATPase for the first time. Cytc has a biphasic action on Na,K-ATPase, with activation at low concentrations (0.06 ng/ml; 6 ng/ml) and inhibition at high concentration (120 ng/ml). Cytc moreover displays isoform/subunit specificity and regulates the Na+ form of the enzyme, while having no effect on the activity or kinetic parameters of the K+-dependent form of the enzyme. Changing the affinity of p-chloromercuribenzoic acid (PCMB) by Cytc is therefore both a required and sufficient condition for confirming that PCMB and Cytc share the same target, namely the thiol groups of cysteine in Na,K-ATPase.

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?

Unusually Divergent Ubiquitin Genes and Proteins in Plasmodium Species

Ubiquitin is an extraordinarily highly conserved 76 amino acid protein encoded by three different types of gene, where the primary translation products are fusions either of ubiquitin with one of two ribosomal proteins (RPs) or of multiple ubiquitin monomers from head to tail. Here, we investigate the evolution of ubiquitin genes in mammalian malaria parasites (Plasmodium species). The ubiquitin encoded by the RPS27a fusion gene is highly divergent, as previously found in a variety of protists. However, we also find that two other forms of divergent ubiquitin sequence, each previously thought to be extremely rare, have arisen recently during the divergence of Plasmodium subgenera. On two occasions, in two distinct lineages, the ubiquitin encoded by the RPL40 fusion gene has rapidly diverged. In addition, in one of these lineages, the polyubiquitin genes have undergone a single codon insertion, previously considered a unique feature of Rhizaria. There has been disagreement whether the multiple ubiquitin coding repeats within a genome exhibit concerted evolution or undergo a birth-and-death process; the Plasmodium ubiquitin genes show clear signs of concerted evolution, including the spread of this codon insertion to multiple repeats within the polyubiquitin gene.

Correlation between protein abundance and sequence conservation: what do recent experiments say?

Cells evolve in a space of parameter values set by physical and chemical forces. These constraints create associations among cellular properties. A particularly strong association is the negative correlation between the rate of evolution of proteins and their abundance in the cell. Highly expressed proteins evolve slower than lowly expressed ones. Multiple hypotheses have been put forward to explain this relationship, including, for instance, the requirement for higher mRNA stability, misfolding avoidance, and misinteraction avoidance for highly expressed proteins. Here, we review some of these hypotheses, their predictions, and how they are supported to finally discuss recent experiments that have been performed to test these predictions.

Investigating the reliability of molecular estimates of evolutionary time when substitution rates and speciation rates vary

Background

An accurate timescale of evolutionary history is essential to testing hypotheses about the influence of historical events and processes, and the timescale for evolution is increasingly derived from analysis of DNA sequences. But variation in the rate of molecular evolution complicates the inference of time from DNA. Evidence is growing for numerous factors, such as life history and habitat, that are linked both to the molecular processes of mutation and fixation and to rates of macroevolutionary diversification. However, the most widely used methods rely on idealised models of rate variation, such as the uncorrelated and autocorrelated clocks, and molecular dating methods are rarely tested against complex models of rate change. One relationship that is not accounted for in molecular dating is the potential for interaction between molecular substitution rates and speciation, a relationship that has been supported by empirical studies in a growing number of taxa. If these relationships are as widespread as current evidence suggests, they may have a significant influence on molecular dates.

Results

We simulate phylogenies and molecular sequences under three different realistic rate variation models—one in which speciation rates and substitution rates both vary but are unlinked, one in which they covary continuously and one punctuated model in which molecular change is concentrated in speciation events, using empirical case studies to parameterise realistic simulations. We test three commonly used “relaxed clock” molecular dating methods against these realistic simulations to explore the degree of error in molecular dates under each model. We find average divergence time inference errors ranging from 12% of node age for the unlinked model when reconstructed under an uncorrelated rate prior using BEAST 2, to up to 91% when sequences evolved under the punctuated model are reconstructed under an autocorrelated prior using PAML.

Conclusions

We demonstrate the potential for substantial errors in molecular dates when both speciation rates and substitution rates vary between lineages. This study highlights the need for tests of molecular dating methods against realistic models of rate variation generated from empirical parameters and known relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02015-8.

Contiguously hydrophobic sequences are functionally significant throughout the human exome

SignificanceProteins rely on the hydrophobic effect to maintain structure and interactions with the environment. Surprisingly, natural selection on amino acid hydrophobicity has not been detected using modern genetic data. Analyses that treat each amino acid separately do not reveal significant results, which we confirm here. However, because the hydrophobic effect becomes more powerful as more hydrophobic molecules are introduced, we tested whether unbroken stretches of hydrophobic amino acids are under selection. Using genetic variant data from across the human genome, we find evidence that selection increases with the length of the unbroken hydrophobic sequence. These results could lead to improvements in a wide range of genomic tools as well as insights into protein-aggregation disease etiology and protein evolutionary history.

Cytochrome c: Using Biological Insight toward Engineering an Optimized Anticancer Biodrug

The heme protein cytochrome c (Cyt c) plays pivotal roles in cellular life and death processes. In the respiratory chain of mitochondria, it serves as an electron transfer protein, contributing to the proliferation of healthy cells. In the cell cytoplasm, it activates intrinsic apoptosis to terminate damaged cells. Insight into these mechanisms and the associated physicochemical properties and biomolecular interactions of Cyt c informs on the anticancer therapeutic potential of the protein, especially in its ability to subvert the current limitations of small molecule-based chemotherapy. In this review, we explore the development of Cyt c as an anticancer drug by identifying cancer types that would be receptive to the cytotoxicity of the protein and factors that can be finetuned to enhance its apoptotic potency. To this end, some information is obtained by characterizing known drugs that operate, in part, by triggering Cyt c induced apoptosis. The application of different smart drug delivery systems is surveyed to highlight important features for maintaining Cyt c stability and activity and improving its specificity for cancer cells and high drug payload release while recognizing the continuing limitations. This work serves to elucidate on the optimization of the strategies to translate Cyt c to the clinical market.

Extracellular domains of transmembrane proteins defy the expression level-evolutionary rate anticorrelation

Highly expressed proteins tend to evolve slowly, a trend known as the expression level-rate of evolution (E-R) anticorrelation. Whereas the reasons for this anticorrelation remain unclear, the most influential hypotheses attribute it to highly expressed proteins being subjected to strong selective pressures to avoid misfolding and/or misinteraction. In accordance with these hypotheses, work in our laboratory has recently shown that extracellular (secreted) proteins lack an E-R anticorrelation (or exhibit a weaker than usual E-R anticorrelation). Extracellular proteins are folded inside the endoplasmic reticulum, where enhanced quality control of folding mechanisms exist, and function in the extracellular space, where misinteraction is unlikely to occur or to produce deleterious effects. Transmembrane proteins contain both intracellular domains (which are folded and function in the cytosol) and extracellular domains (which complete their folding in the endoplasmic reticulum and function in the extracellular space). We thus hypothesized that the extracellular domains of transmembrane proteins should exhibit a weaker E-R anticorrelation than their intracellular domains. Our analyses of human, Saccharomyces and Arabidopsis transmembrane proteins allowed us to confirm our hypothesis. Our results are in agreement with models attributing the E-R anticorrelation to the deleterious effects of misfolding and/or misinteraction.

Interaction of human hemoglobin (HHb) and cytochrome c (Cyt c) with biogenic chloroxine-conjugated silver nanoflowers: spectroscopic and molecular docking approaches

In this research, the biological activity of the antibacterial drug Chloroxine-conjugated biogenic AgNPs (COX-AgNPs) was investigated in simulated physiological conditions (pH = 7.40). Different spectroscopic methods such as UV-visible, fluorescence, and circular dichroism spectroscopic and docking simulation were employed to evaluate the structural changes in the most important blood proteins (human hemoglobin (HHb) and Cytochrome c (Cyt c)) in the presence of COX-AgNPs. The results showed that the COX-AgNPs can bind to HHb and Cyt c and the secondary structure of these proteins remains unchanged, which is crucial in providing insights into the side effects of newly synthesized drugs on their carriers.Communicated by Ramaswamy H. Sarma.

Inferring the deep past from molecular data

There is an expectation that analyses of molecular sequences might be able to distinguish between alternative hypotheses for ancient relationships, but the phylogenetic methods used and types of data analyzed are of critical importance in any attempt to recover historical signal. Here, we discuss some common issues that can influence the topology of trees obtained when using overly simple models to analyze molecular data that often display complicated patterns of sequence heterogeneity. To illustrate our discussion, we have used three examples of inferred relationships which have changed radically as models and methods of analysis have improved. In two of these examples, the sister-group relationship between thermophilic Thermus and mesophilic Deinococcus, and the position of long-branch Microsporidia among eukaryotes, we show that recovering what is now generally considered to be the correct tree is critically dependent on the fit between model and data. In the third example, the position of eukaryotes in the tree of life, the hypothesis that is currently supported by the best available methods is fundamentally different from the classical view of relationships between major cellular domains. Since heterogeneity appears to be pervasive and varied among all molecular sequence data, and even the best available models can still struggle to deal with some problems, the issues we discuss are generally relevant to phylogenetic analyses. It remains essential to maintain a critical attitude to all trees as hypotheses of relationship that may change with more data and better methods.

The Roles of Protein Structure, Taxon Sampling, and Model Complexity in Phylogenomics: A Case Study Focused on Early Animal Divergences

Despite the long history of using protein sequences to infer the tree of life, the potential for different parts of protein structures to retain historical signal remains unclear. We propose that it might be possible to improve analyses of phylogenomic datasets by incorporating information about protein structure. We test this idea using the position of the root of Metazoa (animals) as a model system. We examined the distribution of “strongly decisive” sites (alignment positions that support a specific tree topology) in a dataset comprising >1500 proteins and almost 100 taxa. The proportion of each class of strongly decisive sites in different structural environments was very sensitive to the model used to analyze the data when a limited number of taxa were used but they were stable when taxa were added. As long as enough taxa were analyzed, sites in all structural environments supported the same topology regardless of whether standard tree searches or decisive sites were used to select the optimal tree. However, the use of decisive sites revealed a difference between the support for minority topologies for sites in different structural environments: buried sites and sites in sheet and coil environments exhibited equal support for the minority topologies, whereas solvent-exposed and helix sites had unequal numbers of sites, supporting the minority topologies. This suggests that the relatively slowly evolving buried, sheet, and coil sites are giving an accurate picture of the true species tree and the amount of conflict among gene trees. Taken as a whole, this study indicates that phylogenetic analyses using sites in different structural environments can yield different topologies for the deepest branches in the animal tree of life and that analyzing larger numbers of taxa eliminates this conflict. More broadly, our results highlight the desirability of incorporating information about protein structure into phylogenomic analyses.

Richard Dickerson, Molecular Clocks, and Rates of Protein Evolution

Proteins approximately behave as molecular clocks, accumulating amino acid replacements at a more or less constant rate. Nonetheless, each protein displays a characteristic rate of evolution: whereas some proteins remain largely unaltered over large periods of time, others can rapidly accumulate amino acid replacements. An article by Richard Dickerson, published in the first issue of the Journal of Molecular Evolution (J Mol Evol 1:26-45, 1971), described the first analysis in which the rates of evolution of many proteins were compared, and the differences were interpreted in the light of their function. When comparing the sequences of fibrinopeptides, hemoglobin, and cytochrome c of different species, he observed a linear relationship between the number of amino acid replacements and divergence time. Remarkably, fibrinopeptides had evolved fast, cytochrome c had evolved slowly, and hemoglobin exhibited an intermediate rate of evolution. As the Journal of Molecular Evolution celebrates its 50th anniversary, I highlight this landmark article and reflect on its impact on the field of Molecular Evolution.

On the transformation of MinHash-based uncorrected distances into proper evolutionary distances for phylogenetic inference

Recently developed MinHash-based techniques were proven successful in quickly estimating the level of similarity between large nucleotide sequences. This article discusses their usage and limitations in practice to approximating uncorrected distances between genomes, and transforming these pairwise dissimilarities into proper evolutionary distances. It is notably shown that complex distance measures can be easily approximated using simple transformation formulae based on few parameters. MinHash-based techniques can therefore be very useful for implementing fast yet accurate alignment-free phylogenetic reconstruction procedures from large sets of genomes. This last point of view is assessed with a simulation study using a dedicated bioinformatics tool.

The management of embedded metal fragment patients and the role of chelation Therapy: A workshop of the Department of Veterans Affairs-Walter Reed National Medical Center

Exposure to retained metal fragments from war-related injuries can result in increased systemic metal concentrations, thereby posing potential health risks to target organs far from the site of injury. Given the large number of veterans who have retained fragments and the lack of clear guidance on how to medically manage these individuals, the Department of Veterans Affairs (VA) convened a meeting of chelation experts and clinicians who care for embedded fragment patients to discuss current practices and provide medical management guidance. Based on this group's clinical expertise and review of published literature, the evidence presented suggests that, at least in the case of lead fragments, short-term chelation therapy may be beneficial for embedded fragment patients experiencing acute symptoms associated with metal toxicity; however, in the absence of clinical symptoms or significantly elevated blood lead concentrations (greater than 80 µg/dL), chelation therapy may offer little to no benefit for individuals with retained fragments and pose greater risks due to remobilization of metals stored in bone and other soft tissues. The combination of periodic biomonitoring to assess metal body burden, longitudinal fragment imaging, and selective fragment removal when metal concentrations approach critical injury thresholds offers a more conservative management approach to caring for patients with embedded fragments.

Evolutionary Models for the Diversification of Placental Mammals Across the KPg Boundary

Deciphering the timing of the placental mammal radiation is a longstanding problem in evolutionary biology, but consensus on the tempo and mode of placental diversification remains elusive. Nevertheless, an accurate timetree is essential for understanding the role of important events in Earth history (e.g., Cretaceous Terrestrial Revolution, KPg mass extinction) in promoting the taxonomic and ecomorphological diversification of Placentalia. Archibald and Deutschman described three competing models for the diversification of placental mammals, which are the Explosive, Long Fuse, and Short Fuse Models. More recently, the Soft Explosive Model and Trans-KPg Model have emerged as additional hypotheses for the placental radiation. Here, we review molecular and paleontological evidence for each of these five models including the identification of general problems that can negatively impact divergence time estimates. The Long Fuse Model has received more support from relaxed clock studies than any of the other models, but this model is not supported by morphological cladistic studies that position Cretaceous eutherians outside of crown Placentalia. At the same time, morphological cladistics has a poor track record of reconstructing higher-level relationships among the orders of placental mammals including the results of new pseudoextinction analyses that we performed on the largest available morphological data set for mammals (4,541 characters). We also examine the strengths and weaknesses of different timetree methods (node dating, tip dating, and fossilized birth-death dating) that may now be applied to estimate the timing of the placental radiation. While new methods such as tip dating are promising, they also have problems that must be addressed if these methods are to effectively discriminate among competing hypotheses for placental diversification. Finally, we discuss the complexities of timetree estimation when the signal of speciation times is impacted by incomplete lineage sorting (ILS) and hybridization. Not accounting for ILS results in dates that are older than speciation events. Hybridization, in turn, can result in dates than are younger or older than speciation dates. Disregarding this potential variation in "gene" history across the genome can distort phylogenetic branch lengths and divergence estimates when multiple unlinked genomic loci are combined together in a timetree analysis.

Molecular function limits divergent protein evolution on planetary timescales

Functional conservation is known to constrain protein evolution. Nevertheless, the long-term divergence patterns of proteins maintaining the same molecular function and the possible limits of this divergence have not been explored in detail. We investigate these fundamental questions by characterizing the divergence between ancient protein orthologs with conserved molecular function. Our results demonstrate that the decline of sequence and structural similarities between such orthologs significantly slows down after ~1-2 billion years of independent evolution. As a result, the sequence and structural similarities between ancient orthologs have not substantially decreased for the past billion years. The effective divergence limit (>25% sequence identity) is not primarily due to protein sites universally conserved in all linages. Instead, less than four amino acid types are accepted, on average, per site across orthologous protein sequences. Our analysis also reveals different divergence patterns for protein sites with experimentally determined small and large fitness effects of mutations.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Weak Chemical Interactions That Drive Protein Evolution: Crowding, Sticking, and Quinary Structure in Folding and Function

In recent years, better instrumentation and greater computing power have enabled the imaging of elusive biomolecule dynamics in cells, driving many advances in understanding the chemical organization of biological systems. The focus of this Review is on interactions in the cell that affect both biomolecular stability and function and modulate them. The same protein or nucleic acid can behave differently depending on the time in the cell cycle, the location in a specific compartment, or the stresses acting on the cell. We describe in detail the crowding, sticking, and quinary structure in the cell and the current methods to quantify them both in vitro and in vivo. Finally, we discuss protein evolution in the cell in light of current biophysical evidence. We describe the factors that drive protein evolution and shape protein interaction networks. These interactions can significantly affect the free energy, ΔG, of marginally stable and low-population proteins and, due to epistasis, direct the evolutionary pathways in an organism. We finally conclude by providing an outlook on experiments to come and the possibility of collaborative evolutionary biology and biophysical efforts.

The Molecular Clock in the Evolution of Protein Structures

The molecular clock hypothesis, which states that substitutions accumulate in protein sequences at a constant rate, plays a fundamental role in molecular evolution but it is violated when selective or mutational processes vary with time. Such violations of the molecular clock have been widely investigated for protein sequences, but not yet for protein structures. Here, we introduce a novel statistical test (Significant Clock Violations) and perform a large scale assessment of the molecular clock in the evolution of both protein sequences and structures in three large superfamilies. After validating our method with computer simulations, we find that clock violations are generally consistent in sequence and structure evolution, but they tend to be larger and more significant in structure evolution. Moreover, changes of function assessed through Gene Ontology and InterPro terms are associated with large and significant clock violations in structure evolution. We found that almost one third of significant clock violations are significant in structure evolution but not in sequence evolution, highlighting the advantage to use structure information for assessing accelerated evolution and gathering hints of positive selection. Clock violations between closely related pairs are frequently significant in sequence evolution, consistent with the observed time dependence of the substitution rate attributed to segregation of neutral and slightly deleterious polymorphisms, but not in structure evolution, suggesting that these substitutions do not affect protein structure although they may affect stability. These results are consistent with the view that natural selection, both negative and positive, constrains more strongly protein structures than protein sequences. Our code for computing clock violations is freely available at https://github.com/ugobas/Molecular_clock.

Interaction of bisphenol A 3, 4-quinone metabolite with human hemoglobin, human serum albumin and cytochrome c in vitro

Since covalent protein-bisphenol A adducts generated by the interaction of protein nucleophiles with bisphenol A quinone affect the physicochemical properties of proteins in functional foods and biological tissues, it has become a hot topic nowadays. Therefore, we investigated the interaction of several different biomacromolecules such as hemoglobin, human serum albumin and cytochrome c with bisphenol A 3, 4-quinone (BPAQ). The effects of binding on changes in biomolecular structure were determined by various spectroscopic methods. BPAQ effects were investigated by using the UV-Vis spectroscopy and the quenching phenomenon from fluorescence emission. It proved that the formation of bio-complex and their aromatic micro-environment was likely to be disturbed with as well. Changes observed in circular dichroism (CD) spectroscopy confirmed the quantitative loss of the alpha-helical structure. Further studies with matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOFMS) and molecular docking indicated combining ratio and binding sites between proteins and BPAQ. The in vitro data of BPAQ-proteins adducts may provide a valuable theoretical basis for the elucidation of the toxicological mechanisms of BPAQ adducts in biological systems and environments.

Six Impossible Things before Breakfast: Assumptions, Models, and Belief in Molecular Dating

Confidence in molecular dating analyses has grown with the increasing sophistication of the methods. Some problematic cases where molecular dates disagreed with paleontological estimates appear to have been resolved with a growing agreement between molecules and fossils. But we cannot relax just yet. The growing analytical sophistication of many molecular dating methods relies on an increasingly large number of assumptions about evolutionary history and processes. Many of these assumptions are based on statistical tractability rather than being informed by improved understanding of molecular evolution, yet changing the assumptions can influence molecular dates. How can we tell if the answers we get are driven more by the assumptions we make than by the molecular data being analyzed?

N-glycoproteins exhibit a positive expression level-evolutionary rate correlation

The different proteins of any proteome evolve at enormously different rates. One of the primary factors influencing rates of protein evolution is expression level, with highly expressed proteins tending to evolve at slow rates. This phenomenon, known as the expression level-evolutionary rate (E-R) anticorrelation, has been attributed to the abundance-dependent deleterious effects of misfolding or misinteraction. We have recently shown that secreted proteins either lack an E-R anticorrelation or exhibit a significantly reduced E-R anticorrelation. This effect may be due to the strict quality control to which secreted proteins are subject in the endoplasmic reticulum (which is expected to reduce the rate of misfolding and its deleterious effects) or to their extracellular location (expected to reduce the rate of misinteraction and its deleterious effects). Among secreted proteins, N-glycosylated ones are under particularly strong quality control. Here, we investigate how N-linked glycosylation affects the E-R anticorrelation. Strikingly, we observe a positive E-R correlation among N-glycosylated proteins. That is, N-glycoproteins that are highly expressed evolve at faster rates than lowly expressed N-glycoproteins, in contrast to what is observed among intracellular proteins.

The molecular clock and evolutionary timescales

The molecular clock provides a valuable means of estimating evolutionary timescales from genetic and biochemical data. Proposed in the early 1960s, it was first applied to amino acid sequences and immunological measures of genetic distances between species. The molecular clock has undergone considerable development over the years, and it retains profound relevance in the genomic era. In this mini-review, we describe the history of the molecular clock, its impact on evolutionary theory, the challenges brought by evidence of evolutionary rate variation among species, and the statistical models that have been developed to account for these heterogeneous rates of genetic change. We explain how the molecular clock can be used to infer rates and timescales of evolution, and we list some of the key findings that have been obtained when molecular clocks have been applied to genomic data. Despite the numerous challenges that it has faced over the decades, the molecular clock continues to offer the most effective method of resolving the details of the evolutionary timescale of the Tree of Life.

Cytochrome b5 reductase is the component from neuronal synaptic plasma membrane vesicles that generates superoxide anion upon stimulation by cytochrome c

In this work, we measured the effect of cytochrome c on the NADH-dependent superoxide anion production by synaptic plasma membrane vesicles from rat brain. In these membranes, the cytochrome c stimulated NADH-dependent superoxide anion production was inhibited by antibodies against cytochrome b5 reductase linking the production to this enzyme. Measurement of the superoxide anion radical generated by purified recombinant soluble and membrane cytochrome b5 reductase corroborates the production of the radical by different enzyme isoforms. In the presence of cytochrome c, a burst of superoxide anion as well as the reduction of cytochrome c by cytochrome b5 reductase was measured. Complex formation between both proteins suggests that cytochrome b5 reductase is one of the major partners of cytochrome c upon its release from mitochondria to the cytosol during apoptosis. Superoxide anion production and cytochrome c reduction are the consequences of the stimulated NADH consumption by cytochrome b5 reductase upon complex formation with cytochrome c and suggest a major role of this enzyme as an anti-apoptotic protein during cell death.

Position Matters: Network Centrality Considerably Impacts Rates of Protein Evolution in the Human Protein-Protein Interaction Network

The proteins of any organism evolve at disparate rates. A long list of factors affecting rates of protein evolution have been identified. However, the relative importance of each factor in determining rates of protein evolution remains unresolved. The prevailing view is that evolutionary rates are dominantly determined by gene expression, and that other factors such as network centrality have only a marginal effect, if any. However, this view is largely based on analyses in yeasts, and accurately measuring the importance of the determinants of rates of protein evolution is complicated by the fact that the different factors are often correlated with each other, and by the relatively poor quality of available functional genomics data sets. Here, we use correlation, partial correlation and principal component regression analyses to measure the contributions of several factors to the variability of the rates of evolution of human proteins. For this purpose, we analyzed the entire human protein–protein interaction data set and the human signal transduction network—a network data set of exceptionally high quality, obtained by manual curation, which is expected to be virtually free from false positives. In contrast with the prevailing view, we observe that network centrality (measured as the number of physical and nonphysical interactions, betweenness, and closeness) has a considerable impact on rates of protein evolution. Surprisingly, the impact of centrality on rates of protein evolution seems to be comparable, or even superior according to some analyses, to that of gene expression. Our observations seem to be independent of potentially confounding factors and from the limitations (biases and errors) of interactomic data sets.

Disulfide Bonds Enable Accelerated Protein Evolution

The different proteins of any proteome evolve at enormously different rates. What factors contribute to this variability, and to what extent, is still a largely open question. We hypothesized that disulfide bonds, by increasing protein stability, should make proteins' structures relatively independent of their amino acid sequences, thus acting as buffers of deleterious mutations and enabling accelerated sequence evolution. In agreement with this hypothesis, we observed that membrane proteins with disulfide bonds evolved 88% faster than those without disulfide bonds, and that extracellular proteins with disulfide bonds evolved 49% faster than those without disulfide bonds. In addition, genes encoding proteins with disulfide bonds exhibit an increased likelihood of showing signatures of positive selection. Multivariate analyses indicate that the trend is independent of a number of potentially confounding factors. The effect, however, is not observed among the longest proteins, which can become stabilized by mechanisms other than disulfide bonds.

Correlated evolutionary rates across genomic compartments in Annonaceae

The molecular clock hypothesis is an important concept in biology. Deviations from a constant rate of nucleotide substitution have been found widely among lineages, genomes, genes and individual sites. Phylogenetic research can accommodate for these differences in applying specific models of evolution. Lineage-specific rate heterogeneity however can generate bi- or multimodal distributions of substitution rates across the branches of a tree and this may mislead phylogenetic inferences with currently available models. The plant family Annonaceae is an excellent case to study lineage-specific rate heterogeneity. The two major sister subfamilies, Annonoideae and Malmeoideae, have shown great discrepancies in branch lengths. We used high-throughput sequencing data of 72 genes, 99 spacers and 16 introns from 24 chloroplast genomes and nuclear ribosomal DNA of 23 species to study the molecular rate of evolution in Annonaceae. In all analyses, longer branch lengths and/or higher substitution rates were found for the Annonoideae compared to the Malmeoideae. The Annonaceae had wide variability in chloroplast length, ranging from minimal 175,684bp to 201,723 for Annonoideae and minimal 152,357 to 170,985bp in Malmeoideae, mostly reflecting variation in inverted-repeat length. The Annonoideae showed a higher GC-content in the conserved parts of the chloroplast genome and higher omega (d_N/d_S)-ratios than the Malmeoideae, which could indicate less stringent purifying selection, a pattern that has been found in groups with small population sizes. This study generates new insights into the processes causing lineage-specific rate heterogeneity, which could lead to improved phylogenetic methods.

Substitutional Analysis of Orthologous Protein Families Using BLOCKS

Orthologous proteins, form due to divergence of parental sequence, perform similar function under different environmental and biological conditions. Amino acid changes at locus specific positions form hetero-pairs whose role in BLOCK evolution is yet to be understood. We involve eight protein BLOCKs of known divergence rate to gain insight into the role of hetero-pairs in evolution. Our procedure APBEST uses BLOCK-FASTA file to extract BLOCK specific evolutionary parameters such as dominantly used hetero-pair (D), usage of hetero-pairs (E), non-conservative to conservative substitution ratio (R), maximally-diverse residue (MDR), residue (RD) and class (CD) specific diversity. All these parameters show BLOCK specific variation. Conservative nature of D points towards restoration of function of BLOCK. While E sets the upper-limit of usage of hereto-pairs, strong correlation of R with divergence-rate indicates that the later is directly dependent on non-conservative substitutions. The observation that MDR, measure of positional diversity, occupy very limited positions in BLOCK indicates accommodation of diversity is positionally restricted. Overall, the study extract observed hetero-pair related quantitative and multi-parametric details of BLOCK, which finds application in evolutionary biology.

The impacts of drift and selection on genomic evolution in insects

Genomes evolve through a combination of mutation, drift, and selection, all of which act heterogeneously across genes and lineages. This leads to differences in branch-length patterns among gene trees. Genes that yield trees with the same branch-length patterns can be grouped together into clusters. Here, we propose a novel phylogenetic approach to explain the factors that influence the number and distribution of these gene-tree clusters. We apply our method to a genomic dataset from insects, an ancient and diverse group of organisms. We find some evidence that when drift is the dominant evolutionary process, each cluster tends to contain a large number of fast-evolving genes. In contrast, strong negative selection leads to many distinct clusters, each of which contains only a few slow-evolving genes. Our work, although preliminary in nature, illustrates the use of phylogenetic methods to shed light on the factors driving rate variation in genomic evolution.

Essentiality Is a Strong Determinant of Protein Rates of Evolution during Mutation Accumulation Experiments in Escherichia coli

The Neutral Theory of Molecular Evolution is considered the most powerful theory to understand the evolutionary behavior of proteins. One of the main predictions of this theory is that essential proteins should evolve slower than dispensable ones owing to increased selective constraints. Comparison of genomes of different species, however, has revealed only small differences between the rates of evolution of essential and nonessential proteins. In some analyses, these differences vanish once confounding factors are controlled for, whereas in other cases essentiality seems to have an independent, albeit small, effect. It has been argued that comparing relatively distant genomes may entail a number of limitations. For instance, many of the genes that are dispensable in controlled lab conditions may be essential in some of the conditions faced in nature. Moreover, essentiality can change during evolution, and rates of protein evolution are simultaneously shaped by a variety of factors, whose individual effects are difficult to isolate. Here, we conducted two parallel mutation accumulation experiments in Escherichia coli, during 5,500–5,750 generations, and compared the genomes at different points of the experiments. Our approach (a short-term experiment, under highly controlled conditions) enabled us to overcome many of the limitations of previous studies. We observed that essential proteins evolved substantially slower than nonessential ones during our experiments. Strikingly, rates of protein evolution were only moderately affected by expression level and protein length.

Significance of positive selection and gene duplication in adaptive evolution: in memory of Austin L. Hughes

In this minireview, we highlight the contributions of the late Austin L. Hughes to two areas of molecular evolution: the role of positive (Darwinian) selection, and the impact of gene duplications during genome evolution.

Revisiting a Classic Study of the Molecular Clock

A constant rate of molecular evolution among homologous proteins and across lineages is known as the molecular clock. This concept has been useful for estimating divergence times. Here, we revisit a study by Richard Dickerson (J Mol Evol 1:26-45, 1971), wherein he provided striking visual evidence for a constant rate of amino acid changes among various evolutionary branch points. Dickerson's study is commonly cited as support of the molecular clock and a figure from it is often reproduced in textbooks. Since its publication, however, there have been updates made to dates of common ancestors based on the fossil record that should be considered. Additionally, collecting the accession numbers and carefully outlining Dickerson's methods serves as a resource to students of the molecular clock hypothesis.

Cytochrome c adducts with PCB quinoid metabolites

Polychlorinated biphenyls (PCBs) are a group of 209 individual congeners widely used as industrial chemicals. PCBs are found as by-products in dye and paint manufacture and are legacy, ubiquitous, and persistent as human and environmental contaminants. PCBs with fewer chlorine atoms may be metabolized to hydroxy- and dihydroxy-metabolites and further oxidized to quinoid metabolites both in vitro and in vivo. Specifically, quinoid metabolites may form adducts on nucleophilic sites within cells. We hypothesized that the PCB-quinones covalently bind to cytochrome c and, thereby, cause defects in the function of cytochrome c. In this study, synthetic PCB quinones, 2-(4'-chlorophenyl)-1,4-benzoquinone (PCB3-pQ), 4-4'-chlorophenyl)-1,2-benzoquinone (PCB3-oQ), 2-(3', 5'-dichlorophenyl)-1,4-benzoquinone, 2-(3',4', 5'-trichlorophenyl)-1,4-benzoquinone, and 2-(4'-chlorophenyl)-3,6-dichloro-1,4-benzoquinone, were incubated with cytochrome c, and adducts were detected by liquid chromatography-mass spectrometry (LC-MS) and matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI TOF). Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) was employed to separate the adducted proteins, while trypsin digestion and liquid chromatography-tandem mass spectrometry (LC-MS/MS) were applied to identify the amino acid binding sites on cytochrome c. Conformation change of cytochrome c after binding with PCB3-pQ was investigated by SYBYL-X simulation and cytochrome c function was examined. We found that more than one molecule of PCB-quinone may bind to one molecule of cytochrome c. Lysine and glutamic acid were identified as the predominant binding sites. Software simulation showed conformation changes of adducted cytochrome c. Additionally, cross-linking of cytochrome c was observed on the SDS-PAGE gel. Cytochrome c was found to lose its function as electron acceptor after incubation with PCB quinones. These data provide evidence that the covalent binding of PCB quinone metabolites to cytochrome c may be included among the toxic effects of PCBs.

New thoughts on an old riddle: What determines genetic diversity within and between species?

The question of what determines genetic diversity has long remained unsolved by the modern evolutionary theory (MET). However, it has not deterred researchers from producing interpretations of genetic diversity by using MET. We examine the two observations of genetic diversity made in the 1960s that contributed to the development of MET. The interpretations of these observations by MET are widely known to be inadequate. We review the recent progress of an alternative framework, the maximum genetic diversity (MGD) hypothesis, that uses axioms and natural selection to explain the vast majority of genetic diversity as being at equilibrium that is largely determined by organismal complexity. The MGD hypothesis absorbs the proven virtues of MET and considers its assumptions relevant only to a much more limited scope. This new synthesis has accounted for the overlooked phenomenon of progression towards higher complexity, and more importantly, been instrumental in directing productive research.

Fifty years of co-evolution and beyond: integrating co-evolution from molecules to species

Fifty years after Ehrlich and Raven's seminal paper, the idea of co-evolution continues to grow as a key concept in our understanding of organic evolution. This concept has not only provided a compelling synthesis between evolutionary biology and community ecology, but has also inspired research that extends beyond its original scope. In this article, we identify unresolved questions about the co-evolutionary process and advocate for the integration of co-evolutionary research from molecular to interspecific interactions. We address two basic questions: (i) What is co-evolution and how common is it? (ii) What is the unit of co-evolution? Both questions aim to explore the heart of the co-evolutionary process. Despite the claim that co-evolution is ubiquitous, we argue that there is in fact little evidence to support the view that reciprocal natural selection and coadaptation are common in nature. We also challenge the traditional view that co-evolution only occurs between traits of interacting species. Co-evolution has the potential to explain evolutionary processes and patterns that result from intra- and intermolecular biochemical interactions within cells, intergenomic interactions (e.g. nuclear-cytoplasmic) within species, as well as intergenomic interactions mediated by phenotypic traits between species. Research that bridges across these levels of organization will help to advance our understanding of the importance of the co-evolutionary processes in shaping the diversity of life on Earth.

Expression of temperature-sensitive ion channel TRPM8 in sperm cells correlates with vertebrate evolution

Transient Receptor Potential cation channel, subfamily Melastatin, member 8 (TRPM8) is involved in detection of cold temperature, different noxious compounds and in execution of thermo- as well as chemo-sensitive responses at cellular levels. Here we explored the molecular evolution of TRPM8 by analyzing sequences from various species. We elucidate that several regions of TRPM8 had different levels of selection pressure but the 4th-5th transmembrane regions remain highly conserved. Analysis of synteny suggests that since vertebrate origin, TRPM8 gene is linked with SPP2, a bone morphogen. TRPM8, especially the N-terminal region of it, seems to be highly variable in human population. We found 16,656 TRPM8 variants in 1092 human genomes with top variations being SNPs, insertions and deletions. A total of 692 missense mutations are also mapped to human TRPM8 protein of which 509 seem to be delateroiours in nature as supported by Polyphen V2, SIFT and Grantham deviation score. Using a highly specific antibody, we demonstrate that TRPM8 is expressed endogenously in the testis of rat and sperm cells of different vertebrates ranging from fish to higher mammals. We hypothesize that TRPM8 had emerged during vertebrate evolution (ca 450 MYA). We propose that expression of TRPM8 in sperm cell and its role in regulating sperm function are important factors that have guided its molecular evolution, and that these understandings may have medical importance.

Mammalian genome evolution is governed by multiple pacemakers

Genomic evolution is shaped by a dynamic combination of mutation, selection and genetic drift. These processes lead to evolutionary rate variation across loci and among lineages. In turn, interactions between these two forms of rate variation can produce residual effects, whereby the pattern of among-lineage rate heterogeneity varies across loci. The nature of rate variation is encapsulated in the pacemaker models of genome evolution, which differ in the degree of importance assigned to residual effects: none (Universal Pacemaker), some (Multiple Pacemaker) or total (Degenerate Multiple Pacemaker). Here we use a phylogenetic method to partition the rate variation across loci, allowing comparison of these pacemaker models. Our analysis of 431 genes from 29 mammalian taxa reveals that rate variation across these genes can be explained by 13 pacemakers, consistent with the Multiple Pacemaker model. We find no evidence that these pacemakers correspond to gene function. Our results have important consequences for understanding the factors driving genomic evolution and for molecular-clock analyses.

AVAILABILITY AND IMPLEMENTATION

ClockstaR-G is freely available for download from github (https://github.com/sebastianduchene/clockstarg).

The changing face of the molecular evolutionary clock

The molecular clock has played an important role in biological research, both as a description of the evolutionary process and as a tool for inferring evolutionary timescales. Genomic data have provided valuable insights into the molecular clock, allowing the patterns and causes of evolutionary rate variation to be characterized in increasing detail. I explain how genome sequences offer exciting opportunities for estimating the timescale of the Tree of Life. I describe the different approaches that have been used to deal with the computational and statistical challenges encountered in molecular clock analyses of genomic data. Finally, I offer a perspective on the future of molecular clocks, highlighting some of the key limitations and the most promising research directions.

Molecular phylogenetics before sequences: oligonucleotide catalogs as k-mer spectra

From 1971 to 1985, Carl Woese and colleagues generated oligonucleotide catalogs of 16S/18S rRNAs from more than 400 organisms. Using these incomplete and imperfect data, Carl and his colleagues developed unprecedented insights into the structure, function, and evolution of the large RNA components of the translational apparatus. They recognized a third domain of life, revealed the phylogenetic backbone of bacteria (and its limitations), delineated taxa, and explored the tempo and mode of microbial evolution. For these discoveries to have stood the test of time, oligonucleotide catalogs must carry significant phylogenetic signal; they thus bear re-examination in view of the current interest in alignment-free phylogenetics based on k-mers. Here we consider the aims, successes, and limitations of this early phase of molecular phylogenetics. We computationally generate oligonucleotide sets (e-catalogs) from 16S/18S rRNA sequences, calculate pairwise distances between them based on D₂ statistics, compute distance trees, and compare their performance against alignment-based and k-mer trees. Although the catalogs themselves were superseded by full-length sequences, this stage in the development of computational molecular biology remains instructive for us today.

Pseudogenes and their composers: delving in the 'debris' of human genome

Pseudogenes, the nonfunctional homologs of functional genes and thus exemplified as 'genomic fossils' provide intriguing snapshots of the evolutionary history of human genome. These defunct copies generally arise by retrotransposition or duplication followed by various genetic disablements. In this study, focusing on human pseudogenes and their functional homologues we describe their characteristic features and relevance to protein sequence evolution. We recapitulate that pseudogenes harbor disease-causing degenerative sequence variations in conjunction with the immense disease gene association of their progenitors. Furthermore, we also discuss the issue of functional resurrection and the potentiality observed in some pseudogenes to regulate their functional counterparts.

A scan for human-specific relaxation of negative selection reveals unexpected polymorphism in proteasome genes

Environmental or genomic changes during evolution can relax negative selection pressure on specific loci, permitting high frequency polymorphisms at previously conserved sites. Here, we jointly analyze population genomic and comparative genomic data to search for functional processes showing relaxed negative selection specifically in the human lineage, whereas remaining evolutionarily conserved in other mammals. Consistent with previous studies, we find that olfactory receptor genes display such a signature of relaxation in humans. Intriguingly, proteasome genes also show a prominent signal of human-specific relaxation: multiple proteasome subunits, including four members of the catalytic core particle, contain high frequency nonsynonymous polymorphisms at sites conserved across mammals. Chimpanzee proteasome genes do not display a similar trend. Human proteasome genes also bear no evidence of recent positive or balancing selection. These results suggest human-specific relaxation of negative selection in proteasome subunits; the exact biological causes, however, remain unknown.

Conservation of tubulin-binding sequences in TRPV1 throughout evolution

Background

Transient Receptor Potential Vanilloid sub type 1 (TRPV1), commonly known as capsaicin receptor can detect multiple stimuli ranging from noxious compounds, low pH, temperature as well as electromagnetic wave at different ranges. In addition, this receptor is involved in multiple physiological and sensory processes. Therefore, functions of TRPV1 have direct influences on adaptation and further evolution also. Availability of various eukaryotic genomic sequences in public domain facilitates us in studying the molecular evolution of TRPV1 protein and the respective conservation of certain domains, motifs and interacting regions that are functionally important.

Methodology and Principal Findings

Using statistical and bioinformatics tools, our analysis reveals that TRPV1 has evolved about ∼420 million years ago (MYA). Our analysis reveals that specific regions, domains and motifs of TRPV1 has gone through different selection pressure and thus have different levels of conservation. We found that among all, TRP box is the most conserved and thus have functional significance. Our results also indicate that the tubulin binding sequences (TBS) have evolutionary significance as these stretch sequences are more conserved than many other essential regions of TRPV1. The overall distribution of positively charged residues within the TBS motifs is conserved throughout evolution. In silico analysis reveals that the TBS-1 and TBS-2 of TRPV1 can form helical structures and may play important role in TRPV1 function.

Conclusions and Significance

Our analysis identifies the regions of TRPV1, which are important for structure – function relationship. This analysis indicates that tubulin binding sequence-1 (TBS-1) near the TRP-box forms a potential helix and the tubulin interactions with TRPV1 via TBS-1 have evolutionary significance. This interaction may be required for the proper channel function and regulation and may also have significance in the context of Taxol®-induced neuropathy.