PhyloBot: A Web Portal for Automated Phylogenetics, Ancestral Sequence Reconstruction, and Exploration of Mutational Trajectories

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Molecular advances to study the function, evolution and spectral tuning of arthropod visual opsins

Visual opsins of vertebrates and invertebrates diversified independently and converged to detect ultraviolet to long wavelengths (LW) of green or red light. In both groups, colour vision largely derives from opsin number, expression patterns and changes in amino acids interacting with the chromophore. Functional insights regarding invertebrate opsin evolution have lagged behind those for vertebrates because of the disparity in genomic resources and the lack of robust in vitro systems to characterize spectral sensitivities. Here, we review bioinformatic approaches to identify and model functional variation in opsins as well as recently developed assays to measure spectral phenotypes. In particular, we discuss how transgenic lines, cAMP-spectroscopy and sensitive heterologous expression platforms are starting to decouple genotype-phenotype relationships of LW opsins to complement the classical physiological-behavioural-phylogenetic toolbox of invertebrate visual sensory studies. We illustrate the use of one heterologous method by characterizing novel LW Gq opsins from 10 species, including diurnal and nocturnal Lepidoptera, a terrestrial dragonfly and an aquatic crustacean, expressing them in HEK293T cells, and showing that their maximum absorbance spectra (λmax) range from 518 to 611 nm. We discuss the advantages of molecular approaches for arthropods with complications such as restricted availability, lateral filters, specialized photochemistry and/or electrophysiological constraints. This article is part of the theme issue 'Understanding colour vision: molecular, physiological, neuronal and behavioural studies in arthropods'.

The Unusual Cosubstrate Specificity of NQO2: Conservation Throughout the Amniotes and Implications for Cellular Function

Human Quinone Reductase 2 (NQO2) is a pharmacological target and has appeared in numerous screening efforts as an off-target interactor with kinase-targeted drugs. However the cellular functions of NQO2 are not known. To gain insight into the potential cellular functions of NQO2, we have carried out a detailed evolutionary analysis. One of the most striking characteristics of NQO2 is that it uses conventional dihydronicotinamide cosubstrates, NADH and NADPH, extremely inefficiently, raising questions about an enzymatic function in cells. To characterize the ability of NQO2 to serve as an enzyme, the NQO2 gene was disrupted in HCT116 cells. These NQO2 knockouts along with the parental cells were used to demonstrate that cellular NQO2 is unable to catalyze the activation of the DNA cross-linking reagent, CB1954, without the addition of exogenous dihydronicotinamide riboside (NRH). To find whether the unusual cosubstrate specificity of NQO2 has been conserved in the amniotes, recombinant NQO2 from a reptile, Alligator mississippiensis, and a bird, Anas platyrhynchos, were cloned, purified, and their catalytic activity characterized. Like the mammalian enzymes, the reptile and bird NQO2 were efficient catalysts with the small and synthetic cosubstrate N-benzyl-1,4-dihydronicotinamide but were inefficient in their use of NADH and NADPH. Therefore, the unusual cosubstrate preference of NQO2 appears to be conserved throughout the amniotes; however, we found that NQO2 is not well-conserved in the amphibians. A phylogenetic analysis indicates that NQO1 and NQO2 diverged at the time, approximately 450 MYA, when tetrapods were beginning to evolve.

Resurrected Rubisco suggests uniform carbon isotope signatures over geologic time

The earliest geochemical indicators of microbes-and the enzymes that powered them-extend back ∼3.8 Ga on Earth. Paleobiologists often attempt to understand these indicators by assuming that the behaviors of extant microbes and enzymes are uniform with those of their predecessors. This consistency in behavior seems at odds with our understanding of the inherent variability of living systems. Here, we examine whether a uniformitarian assumption for an enzyme thought to generate carbon isotope indicators of biological activity, RuBisCO, can be corroborated by independently studying the history of changes recorded within RuBisCO's genetic sequences. We resurrected a Precambrian-age RuBisCO by engineering its ancient DNA inside a cyanobacterium genome and measured the engineered organism's fitness and carbon-isotope-discrimination profile. Results indicate that Precambrian uniformitarian assumptions may be warranted but with important caveats. Experimental studies illuminating early innovations are crucial to explore the molecular foundations of life's earliest traces.

Resurrecting Enzymes by Ancestral Sequence Reconstruction

Ancestral Sequence Reconstruction (ASR) allows one to infer the sequences of extinct proteins using the phylogeny of extant proteins. It consists of disclosing the evolutionary history-i.e., the phylogeny-of a protein family of interest and then inferring the sequences of its ancestors-i.e., the nodes in the phylogeny. Assisted by gene synthesis, the selected ancestors can be resurrected in the lab and experimentally characterized. The crucial step to succeed with ASR is starting from a reliable phylogeny. At the same time, it is of the utmost importance to have a clear idea on the evolutionary history of the family under study and the events that influenced it. This allows us to implement ASR with well-defined hypotheses and to apply the appropriate experimental methods. In the last years, ASR has become popular to test hypotheses about the origin of functionalities, changes in activities, understanding physicochemical properties of proteins, among others. In this context, the aim of this chapter is to present the ASR approach applied to the reconstruction of enzymes-i.e., proteins with catalytic roles. The spirit of this contribution is to provide a basic, hands-to-work guide for biochemists and biologists who are unfamiliar with molecular phylogenetics.

Resurrecting Golgi proteins to grasp Golgi ribbon formation and self-association under stress

The Golgi complex is an essential organelle of the eukaryotic exocytic pathway. A subfamily of Golgi matrix proteins, called GRASPs, is central in stress-induced unconventional secretion, Golgi dynamics during mitosis/apoptosis, and Golgi ribbon formation. The Golgi ribbon is vertebrate-specific and correlates with the appearance of two GRASP paralogues and two Golgins (GM130/Golgin45), which form specific GRASP-Golgin pairs. The molecular details of their appearance only in Metazoans are unknown. Moreover, despite new functionalities supported by GRASP paralogy, little is known about their structural and evolutionary differences. Here, we used ancestor sequence reconstruction and biophysical/biochemical approaches to assess the evolution of GRASPs structure/dynamics, fibrillation, and how they started anchoring their Golgin partners. Our data showed that a GRASP ancestor anchored Golgins before gorasp gene duplication in Metazoans. After gene duplication, variations within the GRASP binding pocket determined which paralogue would recruit which Golgin. These interactions are responsible for their specific Golgi location and Golgi ribbon appearance. We also suggest that GRASPs have a long-standing capacity to form supramolecular structures, affecting their participation in stress-induced processes.

Lanza V, Rodríguez-Beltrán J, Galán JC, San Millán A, Cantón R, Coque TM. Evolutionary Pathways and Trajectories in Antibiotic Resistance

Evolution is the hallmark of life. Descriptions of the evolution of microorganisms have provided a wealth of information, but knowledge regarding "what happened" has precluded a deeper understanding of "how" evolution has proceeded, as in the case of antimicrobial resistance. The difficulty in answering the "how" question lies in the multihierarchical dimensions of evolutionary processes, nested in complex networks, encompassing all units of selection, from genes to communities and ecosystems. At the simplest ontological level (as resistance genes), evolution proceeds by random (mutation and drift) and directional (natural selection) processes; however, sequential pathways of adaptive variation can occasionally be observed, and under fixed circumstances (particular fitness landscapes), evolution is predictable. At the highest level (such as that of plasmids, clones, species, microbiotas), the systems' degrees of freedom increase dramatically, related to the variable dispersal, fragmentation, relatedness, or coalescence of bacterial populations, depending on heterogeneous and changing niches and selective gradients in complex environments. Evolutionary trajectories of antibiotic resistance find their way in these changing landscapes subjected to random variations, becoming highly entropic and therefore unpredictable. However, experimental, phylogenetic, and ecogenetic analyses reveal preferential frequented paths (highways) where antibiotic resistance flows and propagates, allowing some understanding of evolutionary dynamics, modeling and designing interventions. Studies on antibiotic resistance have an applied aspect in improving individual health, One Health, and Global Health, as well as an academic value for understanding evolution. Most importantly, they have a heuristic significance as a model to reduce the negative influence of anthropogenic effects on the environment.

An Insight of RuBisCO Evolution through a Multilevel Approach

RuBisCO is the most abundant enzyme on earth; it regulates the organic carbon cycle in the biosphere. Studying its structural evolution will help to develop new strategies of genetic improvement in order to increase food production and mitigate CO₂ emissions. In the present work, we evaluate how the evolution of sequence and structure among isoforms I, II and III of RuBisCO defines their intrinsic flexibility and residue-residue interactions. To do this, we used a multilevel approach based on phylogenetic inferences, multiple sequence alignment, normal mode analysis, and molecular dynamics. Our results show that the three isoforms exhibit greater fluctuation in the loop between αB and βC, and also present a positive correlation with loop 6, an important region for enzymatic activity because it regulates RuBisCO conformational states. Likewise, an increase in the flexibility of the loop structure between αB and βC, as well as Lys330 (form II) and Lys322 (form III) of loop 6, is important to increase photosynthetic efficiency. Thus, the cross-correlation dynamics analysis showed changes in the direction of movement of the secondary structures in the three isoforms. Finally, key amino acid residues related to the flexibility of the RuBisCO structure were indicated, providing important information for its enzymatic engineering.

The beauty of biocatalysis: sustainable synthesis of ingredients in cosmetics

Covering: 2015 up to July 2021The market for cosmetics is consumer driven and the desire for green, sustainable and natural ingredients is increasing. The use of isolated enzymes and whole-cell organisms to synthesise these products is congruent with these values, especially when combined with the use of renewable, recyclable or waste feedstocks. The literature of biocatalysis for the synthesis of ingredients in cosmetics in the past five years is herein reviewed.

Alignment-Integrated Reconstruction of Ancestral Sequences Improves Accuracy

Ancestral sequence reconstruction (ASR) uses an alignment of extant protein sequences, a phylogeny describing the history of the protein family and a model of the molecular-evolutionary process to infer the sequences of ancient proteins, allowing researchers to directly investigate the impact of sequence evolution on protein structure and function. Like all statistical inferences, ASR can be sensitive to violations of its underlying assumptions. Previous studies have shown that, whereas phylogenetic uncertainty has only a very weak impact on ASR accuracy, uncertainty in the protein sequence alignment can more strongly affect inferred ancestral sequences. Here, we show that errors in sequence alignment can produce errors in ASR across a range of realistic and simplified evolutionary scenarios. Importantly, sequence reconstruction errors can lead to errors in estimates of structural and functional properties of ancestral proteins, potentially undermining the reliability of analyses relying on ASR. We introduce an alignment-integrated ASR approach that combines information from many different sequence alignments. We show that integrating alignment uncertainty improves ASR accuracy and the accuracy of downstream structural and functional inferences, often performing as well as highly accurate structure-guided alignment. Given the growing evidence that sequence alignment errors can impact the reliability of ASR studies, we recommend that future studies incorporate approaches to mitigate the impact of alignment uncertainty. Probabilistic modeling of insertion and deletion events has the potential to radically improve ASR accuracy when the model reflects the true underlying evolutionary history, but further studies are required to thoroughly evaluate the reliability of these approaches under realistic conditions.

Fully Automated Ancestral Sequence Reconstruction using FireProtASR

Protein evolution and protein engineering techniques are of great interest in basic science and industrial applications such as pharmacology, medicine, or biotechnology. Ancestral sequence reconstruction (ASR) is a powerful technique for probing evolutionary relationships and engineering robust proteins with good thermostability and broad substrate specificity. The following protocol describes the setting up and execution of an automated FireProt^ASR workflow using a dedicated web site. The service allows for inference of ancestral proteins automatically, from a single protein sequence. Once a protein sequence is submitted, the server will build a dataset of homology sequences, perform a multiple sequence alignment (MSA), build a phylogenetic tree, and reconstruct ancestral nodes. The protocol is also highly flexible and allows for multiple forms of input, advanced settings, and the ability to start jobs from: (i) a single sequence, (ii) a set of homologous sequences, (iii) an MSA, and (iv) a phylogenetic tree. This approach automates all necessary steps and offers a way for novices with limited exposure to ASR techniques to improve the properties of a protein of interest. The technique can even be used to introduce catalytic promiscuity into an enzyme. A web server for accessing the fully automated workflow is freely accessible at https://loschmidt.chemi.muni.cz/fireprotasr/. © 2021 Wiley Periodicals LLC. Basic Protocol: ASR using the Web Server FireProt^ASR.

Ancestral sequence reconstruction - An underused approach to understand the evolution of gene function in plants?

Whilst substantial research effort has been placed on understanding the interactions of plant proteins with their molecular partners, relatively few studies in plants - by contrast to work in other organisms - address how these interactions evolve. It is thought that ancestral proteins were more promiscuous than modern proteins and that specificity often evolved following gene duplication and subsequent functional refining. However, ancestral protein resurrection studies have found that some modern proteins have evolved de novo from ancestors lacking those functions. Intriguingly, the new interactions evolved as a consequence of just a few mutations and, as such, acquisition of new functions appears to be neither difficult nor rare, however, only a few of them are incorporated into biological processes before they are lost to subsequent mutations. Here, we detail the approach of ancestral sequence reconstruction (ASR), providing a primer to reconstruct the sequence of an ancestral gene. We will present case studies from a range of different eukaryotes before discussing the few instances where ancestral reconstructions have been used in plants. As ASR is used to dig into the remote evolutionary past, we will also present some alternative genetic approaches to investigate molecular evolution on shorter timescales. We argue that the study of plant secondary metabolism is particularly well suited for ancestral reconstruction studies. Indeed, its ancient evolutionary roots and highly diverse landscape provide an ideal context in which to address the focal issue around the emergence of evolutionary novelties and how this affects the chemical diversification of plant metabolism.

Computational design of enzymes for biotechnological applications

Enzymes are the natural catalysts that execute biochemical reactions upholding life. Their natural effectiveness has been fine-tuned as a result of millions of years of natural evolution. Such catalytic effectiveness has prompted the use of biocatalysts from multiple sources on different applications, including the industrial production of goods (food and beverages, detergents, textile, and pharmaceutics), environmental protection, and biomedical applications. Natural enzymes often need to be improved by protein engineering to optimize their function in non-native environments. Recent technological advances have greatly facilitated this process by providing the experimental approaches of directed evolution or by enabling computer-assisted applications. Directed evolution mimics the natural selection process in a highly accelerated fashion at the expense of arduous laboratory work and economic resources. Theoretical methods provide predictions and represent an attractive complement to such experiments by waiving their inherent costs. Computational techniques can be used to engineer enzymatic reactivity, substrate specificity and ligand binding, access pathways and ligand transport, and global properties like protein stability, solubility, and flexibility. Theoretical approaches can also identify hotspots on the protein sequence for mutagenesis and predict suitable alternatives for selected positions with expected outcomes. This review covers the latest advances in computational methods for enzyme engineering and presents many successful case studies.

Reconstructed evolutionary history of the yeast septins Cdc11 and Shs1

Septins are GTP-binding proteins conserved across metazoans. They can polymerize into extended filaments and, hence, are considered a component of the cytoskeleton. The number of individual septins varies across the tree of life—yeast (Saccharomyces cerevisiae) has seven distinct subunits, a nematode (Caenorhabditis elegans) has two, and humans have 13. However, the overall geometric unit (an apolar hetero-octameric protomer and filaments assembled there from) has been conserved. To understand septin evolutionary variation, we focused on a related pair of yeast subunits (Cdc11 and Shs1) that appear to have arisen from gene duplication within the fungal clade. Either Cdc11 or Shs1 occupies the terminal position within a hetero-octamer, yet Cdc11 is essential for septin function and cell viability, whereas Shs1 is not. To discern the molecular basis of this divergence, we utilized ancestral gene reconstruction to predict, synthesize, and experimentally examine the most recent common ancestor (“Anc.11-S”) of Cdc11 and Shs1. Anc.11-S was able to occupy the terminal position within an octamer, just like the modern subunits. Although Anc.11-S supplied many of the known functions of Cdc11, it was unable to replace the distinct function(s) of Shs1. To further evaluate the history of Shs1, additional intermediates along a proposed trajectory from Anc.11-S to yeast Shs1 were generated and tested. We demonstrate that multiple events contributed to the current properties of Shs1: (1) loss of Shs1–Shs1 self-association early after duplication, (2) co-evolution of heterotypic Cdc11–Shs1 interaction between neighboring hetero-octamers, and (3) eventual repurposing and acquisition of novel function(s) for its C-terminal extension domain. Thus, a pair of duplicated proteins, despite constraints imposed by assembly into a highly conserved multi-subunit structure, could evolve new functionality via a complex evolutionary pathway.

Reconstructing the evolutionary history of nitrogenases: Evidence for ancestral molybdenum-cofactor utilization

The nitrogenase metalloenzyme family, essential for supplying fixed nitrogen to the biosphere, is one of life's key biogeochemical innovations. The three forms of nitrogenase differ in their metal dependence, each binding either a FeMo-, FeV-, or FeFe-cofactor where the reduction of dinitrogen takes place. The history of nitrogenase metal dependence has been of particular interest due to the possible implication that ancient marine metal availabilities have significantly constrained nitrogenase evolution over geologic time. Here, we reconstructed the evolutionary history of nitrogenases, and combined phylogenetic reconstruction, ancestral sequence inference, and structural homology modeling to evaluate the potential metal dependence of ancient nitrogenases. We find that active-site sequence features can reliably distinguish extant Mo-nitrogenases from V- and Fe-nitrogenases and that inferred ancestral sequences at the deepest nodes of the phylogeny suggest these ancient proteins most resemble modern Mo-nitrogenases. Taxa representing early-branching nitrogenase lineages lack one or more biosynthetic nifE and nifN genes that both contribute to the assembly of the FeMo-cofactor in studied organisms, suggesting that early Mo-nitrogenases may have utilized an alternate and/or simplified pathway for cofactor biosynthesis. Our results underscore the profound impacts that protein-level innovations likely had on shaping global biogeochemical cycles throughout the Precambrian, in contrast to organism-level innovations that characterize the Phanerozoic Eon.

Iroki: automatic customization and visualization of phylogenetic trees

Phylogenetic trees are an important analytical tool for evaluating community diversity and evolutionary history. In the case of microorganisms, the decreasing cost of sequencing has enabled researchers to generate ever-larger sequence datasets, which in turn have begun to fill gaps in the evolutionary history of microbial groups. However, phylogenetic analyses of these types of datasets create complex trees that can be challenging to interpret. Scientific inferences made by visual inspection of phylogenetic trees can be simplified and enhanced by customizing various parts of the tree. Yet, manual customization is time-consuming and error prone, and programs designed to assist in batch tree customization often require programming experience or complicated file formats for annotation. Iroki, a user-friendly web interface for tree visualization, addresses these issues by providing automatic customization of large trees based on metadata contained in tab-separated text files. Iroki’s utility for exploring biological and ecological trends in sequencing data was demonstrated through a variety of microbial ecology applications in which trees with hundreds to thousands of leaf nodes were customized according to extensive collections of metadata. The Iroki web application and documentation are available at https://www.iroki.net or through the VIROME portal http://virome.dbi.udel.edu. Iroki’s source code is released under the MIT license and is available at https://github.com/mooreryan/iroki.

Conservation of Structure and Immune Antagonist Functions of Filoviral VP35 Homologs Present in Microbat Genomes

Non-retroviral integrated RNA viral sequences (NIRVs) potentially encoding ∼280 amino acid homologs to filovirus VP35 proteins are present across the Myotis genus of bats. These are estimated to have been maintained for ∼18 million years, indicating their co-option. To address the reasons for co-option, 16 Myotis VP35s were characterized in comparison to VP35s from the extant filoviruses Ebola virus and Marburg virus, in which VP35s play critical roles in immune evasion and RNA synthesis. The Myotis VP35s demonstrated a conserved suppression of innate immune signaling, albeit with reduced potency, in either human or Myotis cells. Their attenuation reflects a lack of dsRNA binding that in the filoviral VP35s correlates with potent suppression of interferon responses. Despite divergent function, evolution has preserved in Myotis the structure of the filoviral VP35s, indicating that this structure is critical for co-opted function, possibly as a regulator of innate immune signaling.

Enamel proteome shows that Gigantopithecus was an early diverging pongine

Gigantopithecus blacki was a giant hominid that inhabited densely forested environments of Southeast Asia during the Pleistocene¹. Its evolutionary relationships to other great ape species, and their divergence during the Middle and Late Miocene (16-5.3 Mya), remains disputed^2,3. Hypotheses regarding relationships between Gigantopithecus and extinct and extant hominids are difficult to substantiate because of its highly derived dentognathic morphology and the absence of cranial and post-cranial remains^1,3-6. Therefore, proposed hypotheses on the phylogenetic position of Gigantopithecus among hominids have been wide-ranging, but none have received independent molecular validation. We retrieved dental enamel proteome sequences from a 1.9 million years (Mya) old Gigantopithecus blacki molar found in Chuifeng Cave, China^7,8. The thermal age of these protein sequences is approximately five times older than any previously published mammalian proteome or genome. We demonstrate that Gigantopithecus is a sister clade to orangutans (genus Pongo) with a common ancestor about 10-12 Mya, implying that the Gigantopithecus divergence from Pongo is part of the Miocene radiation of great apes. Additionally, we hypothesize that the expression of alpha-2-HS-glycoprotein (AHSG), which has not been observed in enamel proteomes previously, had a role in the biomineralization of the thick enamel crowns that characterize the large molars in the genus^9,10. The survival of an Early Pleistocene dental enamel proteome in the subtropics further expands the scope of palaeoproteomic analysis into geographic areas and time periods previously considered incompatible with genetic preservation.

Evolutionary Analysis of JAZ Proteins in Plants: An Approach in Search of the Ancestral Sequence

Jasmonates are phytohormones that regulate development, metabolism and immunity. Signal transduction is critical to activate jasmonate responses, but the evolution of some key regulators such as jasmonate-ZIM domain (JAZ) repressors is not clear. Here, we identified 1065 JAZ sequence proteins in 66 lower and higher plants and analyzed their evolution by bioinformatics methods. We found that the TIFY and Jas domains are highly conserved along the evolutionary scale. Furthermore, the canonical degron sequence LPIAR(R/K) of the Jas domain is conserved in lower and higher plants. It is noteworthy that degron sequences showed a large number of alternatives from gymnosperms to dicots. In addition, ethylene-responsive element binding factor-associated amphiphilic repression (EAR) motifs are displayed in all plant lineages from liverworts to angiosperms. However, the cryptic MYC2-interacting domain (CMID) domain appeared in angiosperms for the first time. The phylogenetic analysis performed using the Maximum Likelihood method indicated that JAZ ortholog proteins are grouped according to their similarity and plant lineage. Moreover, ancestral JAZ sequences were constructed by PhyloBot software and showed specific changes in the TIFY and Jas domains during evolution from liverworts to dicots. Finally, we propose a model for the evolution of the ancestral sequences of the main eight JAZ protein subgroups. These findings contribute to the understanding of the JAZ family origin and expansion in land plants.

Assessing circovirus gene flow in multiple spill-over events

The establishment of viral pathogens in new host environments following spillover events probably requires adaptive changes within both the new host and pathogen. After many generations, signals for ancient cross-species transmission may become lost and a strictly host-adapted phylogeny may mimic true co-divergence while the virus may retain an inherent ability to jump host species. The mechanistic basis for such processes remains poorly understood. To study the dynamics of virus-host co-divergence and the arbitrary chances of spillover in various reservoir hosts with equal ecological opportunity, we examined structural constraints of capsid protein in extant populations of Beak and feather disease virus (BFDV) during known spillover events. By assessing reservoir-based genotype stratification, we identified co-divergence defying signatures in the evolution BFDV which highlighted primordial processes of cryptic host adaptation and competing forces of host co-divergence and cross-species transmission. We demonstrate that, despite extensive surface plasticity gathered over a longer span of evolution, structural constraints of the capsid protein allow opportunistic host switching in host-adapted populations. This study provides new insights into how small populations of endangered psittacine species may face multidirectional forces of infection from reservoirs with apparently co-diverging genotypes.

Ancestral reconstruction reveals mechanisms of ERK regulatory evolution

Protein kinases are crucial to coordinate cellular decisions and therefore their activities are strictly regulated. Previously we used ancestral reconstruction to determine how CMGC group kinase specificity evolved (Howard et al., 2014). In the present study, we reconstructed ancestral kinases to study the evolution of regulation, from the inferred ancestor of CDKs and MAPKs, to modern ERKs. Kinases switched from high to low autophosphorylation activity at the transition to the inferred ancestor of ERKs 1 and 2. Two synergistic amino acid changes were sufficient to induce this change: shortening of the β3-αC loop and mutation of the gatekeeper residue. Restoring these two mutations to their inferred ancestral state led to a loss of dependence of modern ERKs 1 and 2 on the upstream activating kinase MEK in human cells. Our results shed light on the evolutionary mechanisms that led to the tight regulation of a kinase that is central in development and disease.

Exploring the sequence, function, and evolutionary space of protein superfamilies using sequence similarity networks and phylogenetic reconstructions

Integrative computational methods can facilitate the discovery of new protein functions and enzymatic reactions by enabling the observation and investigation of complex sequence-structure-function and evolutionary relationships within protein superfamilies. Here, we highlight the use of sequence similarity networks (SSNs) and phylogenetic reconstructions to map the functional divergence and evolutionary history of protein superfamilies. We exemplify this approach using the nitroreductase (NTR) flavoenzyme superfamily, demonstrating that SSN investigations can provide a rapid and effective means to classify groups of proteins, expose sequence similarity relationships across the global scale of a protein superfamily, and efficiently support detailed phylogenetic analyses. Integration of such approaches with systematic experimental characterization will expand our understanding of the functional diversity of enzymes, their evolution, and their associated physiological roles.

Crystal structure of phosphoribulokinase from Synechococcus sp. strain PCC 6301

Phosphoribulokinase (PRK) catalyses the ATP-dependent phosphorylation of ribulose 5-phosphate to give ribulose 1,5-bisphosphate. Regulation of this reaction in response to light controls carbon fixation during photosynthesis. Here, the crystal structure of PRK from the cyanobacterium Synechococcus sp. strain PCC 6301 is presented. The enzyme is dimeric and has an α/β-fold with an 18-stranded β-sheet at its core. Interestingly, a disulfide bond is found between Cys40 and the P-loop residue Cys18, revealing the structural basis for the redox inactivation of PRK activity. A second disulfide bond appears to rigidify the dimer interface and may thereby contribute to regulation by the adaptor protein CP12 and glyceraldehyde-3-phosphate dehydrogenase.

Sequence selection by FitSS4ASR alleviates ancestral sequence reconstruction as exemplified for geranylgeranylglyceryl phosphate synthase

For evolutionary studies, but also for protein engineering, ancestral sequence reconstruction (ASR) has become an indispensable tool. The first step of every ASR protocol is the preparation of a representative sequence set containing at most a few hundred recent homologs whose composition determines decisively the outcome of a reconstruction. A common approach for sequence selection consists of several rounds of manual recompilation that is driven by embedded phylogenetic analyses of the varied sequence sets. For ASR of a geranylgeranylglyceryl phosphate synthase, we additionally utilized FitSS4ASR, which replaces this time-consuming protocol with an efficient and more rational approach. FitSS4ASR applies orthogonal filters to a set of homologs to eliminate outlier sequences and those bearing only a weak phylogenetic signal. To demonstrate the usefulness of FitSS4ASR, we determined experimentally the oligomerization state of eight predecessors, which is a delicate and taxon-specific property. Corresponding ancestors deduced in a manual approach and by means of FitSS4ASR had the same dimeric or hexameric conformation; this concordance testifies to the efficiency of FitSS4ASR for sequence selection. FitSS4ASR-based results of two other ASR experiments were added to the Supporting Information. Program and documentation are available at https://gitlab.bioinf.ur.de/hek61586/FitSS4ASR.

Diversification of DNA binding specificities enabled SREBP transcription regulators to expand the repertoire of cellular functions that they govern in fungi

The Sterol Regulatory Element Binding Proteins (SREBPs) are basic-helix-loop-helix transcription regulators that control the expression of sterol biosynthesis genes in higher eukaryotes and some fungi. Surprisingly, SREBPs do not regulate sterol biosynthesis in the ascomycete yeasts (Saccharomycotina) as this role was handed off to an unrelated transcription regulator in this clade. The SREBPs, nonetheless, expanded in fungi such as the ascomycete yeasts Candida spp., raising questions about their role and evolution in these organisms. Here we report that the fungal SREBPs diversified their DNA binding preferences concomitantly with an expansion in function. We establish that several branches of fungal SREBPs preferentially bind non-palindromic DNA sequences, in contrast to the palindromic DNA motifs recognized by most basic-helix-loop-helix proteins (including SREBPs) in higher eukaryotes. Reconstruction and biochemical characterization of the likely ancestor protein suggest that an intrinsic DNA binding promiscuity in the family was resolved by alternative mechanisms in different branches of fungal SREBPs. Furthermore, we show that two SREBPs in the human commensal yeast Candida albicans drive a transcriptional cascade that inhibits a morphological switch under anaerobic conditions. Preventing this morphological transition enhances C. albicans colonization of the mammalian intestine, the fungus' natural niche. Thus, our results illustrate how diversification in DNA binding preferences enabled the functional expansion of a family of eukaryotic transcription regulators.

Botulinum Neurotoxin Diversity from a Gene-Centered View

Botulinum neurotoxins (BoNTs) rank amongst the most potent toxins known. The factors responsible for the emergence of the many known and yet unknown BoNT variants remain elusive. It also remains unclear why anaerobic bacteria that are widely distributed in our environment and normally do not pose a threat to humans, produce such deadly toxins. Even the possibility of accidental toxicity to humans has not been excluded. Here, I review the notion that BoNTs may have specifically evolved to target vertebrates. Considering the extremely complex molecular architecture of the toxins, which enables them to reach the bloodstream, to recognize and enter neurons, and to block neurotransmitter release, it seems highly unlikely that BoNT toxicity to vertebrates is a coincidence. The carcass⁻maggot cycle provides a plausible explanation for a natural role of the toxins: to enable mass reproduction of bacteria, spores, and toxins, using toxin-unaffected invertebrates, such as fly maggots, as the vectors. There is no clear correlation between toxigenicity and a selective advantage of clostridia in their natural habitat. Possibly, non-toxigenic strains profit from carcasses resulting from the action of toxigenic strains. Alternatively, a gene-centered view of toxin evolution would also explain this observation. Toxin-coding mobile genetic elements may have evolved as selfish genes, promoting their own propagation, similar to commensal viruses, using clostridia and other bacteria as the host. Research addressing the role of BoNTs in nature and the origin of toxin variability goes hand in hand with the identification of new toxin variants and the design of improved toxin variants for medical applications. These research directions may also reveal yet unknown natural antidotes against these extremely potent neurotoxins.

Convergent evolution of tertiary structure in rhodopsin visual proteins from vertebrates and box jellyfish

Box jellyfish and vertebrates are separated by >500 million years of evolution yet have structurally analogous lens eyes that employ rhodopsin photopigments for vision. All opsins possess a negatively charged residue-the counterion-to maintain visible-light sensitivity and facilitate photoisomerization of their retinaldehyde chromophore. In vertebrate rhodopsins, the molecular evolution of the counterion position-from a highly conserved distal location in the second extracellular loop (E181) to a proximal location in the third transmembrane helix (E113)-is established as a key driver of higher fidelity photoreception. Here, we use computational biology and heterologous action spectroscopy to determine whether the appearance of the advanced visual apparatus in box jellyfish was also accompanied by changes in the opsin tertiary structure. We found that the counterion in an opsin from the lens eye of the box jellyfish Carybdea rastonii (JellyOp) has also moved to a unique proximal location within the transmembrane bundle-E94 in TM2. Furthermore, we reveal that this Schiff base/counterion system includes an additional positive charge-R186-that has coevolved with E94 to functionally separate E94 and E181 in the chromophore-binding pocket of JellyOp. By engineering this pocket-neutralizing R186 and E94, or swapping E94 with the vertebrate counterion E113-we can recreate versions of the invertebrate and vertebrate counterion systems, respectively, supporting a relatively similar overall architecture in this region of animal opsins. In summary, our data establish the third only counterion site in animal opsins and reveal convergent evolution of tertiary structure in opsins from distantly related species with advanced visual systems.

Resurrecting ancestral genes in bacteria to interpret ancient biosignatures

Two datasets, the geologic record and the genetic content of extant organisms, provide complementary insights into the history of how key molecular components have shaped or driven global environmental and macroevolutionary trends. Changes in global physico-chemical modes over time are thought to be a consistent feature of this relationship between Earth and life, as life is thought to have been optimizing protein functions for the entirety of its approximately 3.8 billion years of history on the Earth. Organismal survival depends on how well critical genetic and metabolic components can adapt to their environments, reflecting an ability to optimize efficiently to changing conditions. The geologic record provides an array of biologically independent indicators of macroscale atmospheric and oceanic composition, but provides little in the way of the exact behaviour of the molecular components that influenced the compositions of these reservoirs. By reconstructing sequences of proteins that might have been present in ancient organisms, we can downselect to a subset of possible sequences that may have been optimized to these ancient environmental conditions. How can one use modern life to reconstruct ancestral behaviours? Configurations of ancient sequences can be inferred from the diversity of extant sequences, and then resurrected in the laboratory to ascertain their biochemical attributes. One way to augment sequence-based, single-gene methods to obtain a richer and more reliable picture of the deep past, is to resurrect inferred ancestral protein sequences in living organisms, where their phenotypes can be exposed in a complex molecular-systems context, and then to link consequences of those phenotypes to biosignatures that were preserved in the independent historical repository of the geological record. As a first step beyond single-molecule reconstruction to the study of functional molecular systems, we present here the ancestral sequence reconstruction of the beta-carbonic anhydrase protein. We assess how carbonic anhydrase proteins meet our selection criteria for reconstructing ancient biosignatures in the laboratory, which we term palaeophenotype reconstruction.This article is part of the themed issue 'Reconceptualizing the origins of life'.

Constraining the timing of the Great Oxidation Event within the Rubisco phylogenetic tree

Ribulose 1,5-bisphosphate (RuBP) carboxylase/oxygenase (RuBisCO, or Rubisco) catalyzes a key reaction by which inorganic carbon is converted into organic carbon in the metabolism of many aerobic and anaerobic organisms. Across the broader Rubisco protein family, homologs exhibit diverse biochemical characteristics and metabolic functions, but the evolutionary origins of this diversity are unclear. Evidence of the timing of Rubisco family emergence and diversification of its different forms has been obscured by a meager paleontological record of early Earth biota, their subcellular physiology and metabolic components. Here, we use computational models to reconstruct a Rubisco family phylogenetic tree, ancestral amino acid sequences at branching points on the tree, and protein structures for several key ancestors. Analysis of historic substitutions with respect to their structural locations shows that there were distinct periods of amino acid substitution enrichment above background levels near and within its oxygen-sensitive active site and subunit interfaces over the divergence between Form III (associated with anoxia) and Form I (associated with oxia) groups in its evolutionary history. One possible interpretation is that these periods of substitutional enrichment are coincident with oxidative stress exerted by the rise of oxygenic photosynthesis in the Precambrian era. Our interpretation implies that the periods of Rubisco substitutional enrichment inferred near the transition from anaerobic Form III to aerobic Form I ancestral sequences predate the acquisition of Rubisco by fully derived cyanobacterial (i.e., dual photosystem-bearing, oxygen-evolving) clades. The partitioning of extant lineages at high clade levels within our Rubisco phylogeny indicates that horizontal transfer of Rubisco is a relatively infrequent event. Therefore, it is possible that the mutational enrichment periods between the Form III and Form I common ancestral sequences correspond to the adaptation of key oxygen-sensitive components of Rubisco prior to, or coincident with, the Great Oxidation Event.

Exploring the past and the future of protein evolution with ancestral sequence reconstruction: the 'retro' approach to protein engineering

A central goal in molecular evolution is to understand the ways in which genes and proteins evolve in response to changing environments. In the absence of intact DNA from fossils, ancestral sequence reconstruction (ASR) can be used to infer the evolutionary precursors of extant proteins. To date, ancestral proteins belonging to eubacteria, archaea, yeast and vertebrates have been inferred that have been hypothesized to date from between several million to over 3 billion years ago. ASR has yielded insights into the early history of life on Earth and the evolution of proteins and macromolecular complexes. Recently, however, ASR has developed from a tool for testing hypotheses about protein evolution to a useful means for designing novel proteins. The strength of this approach lies in the ability to infer ancestral sequences encoding proteins that have desirable properties compared with contemporary forms, particularly thermostability and broad substrate range, making them good starting points for laboratory evolution. Developments in technologies for DNA sequencing and synthesis and computational phylogenetic analysis have led to an escalation in the number of ancient proteins resurrected in the last decade and greatly facilitated the use of ASR in the burgeoning field of synthetic biology. However, the primary challenge of ASR remains in accurately inferring ancestral states, despite the uncertainty arising from evolutionary models, incomplete sequences and limited phylogenetic trees. This review will focus, firstly, on the use of ASR to uncover links between sequence and phenotype and, secondly, on the practical application of ASR in protein engineering.

Reconstructing Ancient Proteins to Understand the Causes of Structure and Function

A central goal in biochemistry is to explain the causes of protein sequence, structure, and function. Mainstream approaches seek to rationalize sequence and structure in terms of their effects on function and to identify function's underlying determinants by comparing related proteins to each other. Although productive, both strategies suffer from intrinsic limitations that have left important aspects of many proteins unexplained. These limits can be overcome by reconstructing ancient proteins, experimentally characterizing their properties, and retracing their evolution through time. This approach has proven to be a powerful means for discovering how historical changes in sequence produced the functions, structures, and other physical/chemical characteristics of modern proteins. It has also illuminated whether protein features evolved because of functional optimization, historical constraint, or blind chance. Here we review recent studies employing ancestral protein reconstruction and show how they have produced new knowledge not only of molecular evolutionary processes but also of the underlying determinants of modern proteins' physical, chemical, and biological properties.