Evolution of DNA specificity in a transcription factor family produced a new gene regulatory module.

Complex gene regulatory networks require transcription factors (TFs) to bind distinct DNA sequences. To understand how novel TF specificity evolves, we combined phylogenetic, biochemical, and biophysical approaches to interrogate how DNA recognition diversified in the steroid hormone receptor (SR) family. After duplication of the ancestral SR, three mutations in one copy radically weakened binding to the ancestral estrogen response element (ERE) and improved binding to a new set of DNA sequences (steroid response elements, SREs). They did so by establishing unfavorable interactions with ERE and abolishing unfavorable interactions with SRE; also required were numerous permissive substitutions, which nonspecifically improved cooperativity and affinity of DNA binding. Our findings indicate that negative determinants of binding play key roles in TFs' DNA selectivity and-with our prior work on the evolution of SR ligand specificity during the same interval-show how a specific new gene regulatory module evolved without interfering with the integrity of the ancestral module.

A hydrophobic ratchet entrenches molecular complexes.

Most proteins assemble into multisubunit complexes1. The persistence of these complexes across evolutionary time is usually explained as the result of natural selection for functional properties that depend on multimerization, such as intersubunit allostery or the capacity to do mechanical work2. In many complexes, however, multimerization does not enable any known function3. An alternative explanation is that multimers could become entrenched if substitutions accumulate that are neutral in multimers but deleterious in monomers; purifying selection would then prevent reversion to the unassembled form, even if assembly per se does not enhance biological function3-7. Here we show that a hydrophobic mutational ratchet systematically entrenches molecular complexes. By applying ancestral protein reconstruction and biochemical assays to the evolution of steroid hormone receptors, we show that an ancient hydrophobic interface, conserved for hundreds of millions of years, is entrenched because exposure of this interface to solvent reduces protein stability and causes aggregation, even though the interface makes no detectable contribution to function. Using structural bioinformatics, we show that a universal mutational propensity drives sites that are buried in multimeric interfaces to accumulate hydrophobic substitutions to levels that are not tolerated in monomers. In a database of hundreds of families of multimers, most show signatures of long-term hydrophobic entrenchment. It is therefore likely that many protein complexes persist because a simple ratchet-like mechanism entrenches them across evolutionary time, even when they are functionally gratuitous.

Origin of complexity in haemoglobin evolution.

Most proteins associate into multimeric complexes with specific architectures1,2, which often have functional properties such as cooperative ligand binding or allosteric regulation3. No detailed knowledge is available about how any multimer and its functions arose during evolution. Here we use ancestral protein reconstruction and biophysical assays to elucidate the origins of vertebrate haemoglobin, a heterotetramer of paralogous α- and ß-subunits that mediates respiratory oxygen transport and exchange by cooperatively binding oxygen with moderate affinity. We show that modern haemoglobin evolved from an ancient monomer and characterize the historical 'missing link' through which the modern tetramer evolved-a noncooperative homodimer with high oxygen affinity that existed before the gene duplication that generated distinct α- and ß-subunits. Reintroducing just two post-duplication historical substitutions into the ancestral protein is sufficient to cause strong tetramerization by creating favourable contacts with more ancient residues on the opposing subunit. These surface substitutions markedly reduce oxygen affinity and even confer cooperativity, because an ancient linkage between the oxygen binding site and the multimerization interface was already an intrinsic feature of the protein's structure. Our findings establish that evolution can produce new complex molecular structures and functions via simple genetic mechanisms that recruit existing biophysical features into higher-level architectures.

Author Correction: Origin of complexity in haemoglobin evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Evolutionary Shortcuts via Multinucleotide Substitutions and Their Impact on Natural Selection Analyses.

Inference and interpretation of evolutionary processes, in particular of the types and targets of natural selection affecting coding sequences, are critically influenced by the assumptions built into statistical models and tests. If certain aspects of the substitution process (even when they are not of direct interest) are presumed absent or are modeled with too crude of a simplification, estimates of key model parameters can become biased, often systematically, and lead to poor statistical performance. Previous work established that failing to accommodate multinucleotide (or multihit, MH) substitutions strongly biases dN/dS-based inference towards false-positive inferences of diversifying episodic selection, as does failing to model variation in the rate of synonymous substitution (SRV) among sites. Here, we develop an integrated analytical framework and software tools to simultaneously incorporate these sources of evolutionary complexity into selection analyses. We found that both MH and SRV are ubiquitous in empirical alignments, and incorporating them has a strong effect on whether or not positive selection is detected (1.4-fold reduction) and on the distributions of inferred evolutionary rates. With simulation studies, we show that this effect is not attributable to reduced statistical power caused by using a more complex model. After a detailed examination of 21 benchmark alignments and a new high-resolution analysis showing which parts of the alignment provide support for positive selection, we show that MH substitutions occurring along shorter branches in the tree explain a significant fraction of discrepant results in selection detection. Our results add to the growing body of literature which examines decades-old modeling assumptions (including MH) and finds them to be problematic for comparative genomic data analysis. Because multinucleotide substitutions have a significant impact on natural selection detection even at the level of an entire gene, we recommend that selection analyses of this type consider their inclusion as a matter of routine. To facilitate this procedure, we developed, implemented, and benchmarked a simple and well-performing model testing selection detection framework able to screen an alignment for positive selection with two biologically important confounding processes: site-to-site synonymous rate variation, and multinucleotide instantaneous substitutions.

Alternative evolutionary histories in the sequence space of an ancient protein.

To understand why molecular evolution turned out as it did, we must characterize not only the path that evolution followed across the space of possible molecular sequences but also the many alternative trajectories that could have been taken but were not. A large-scale comparison of real and possible histories would establish whether the outcome of evolution represents an optimal state driven by natural selection or the contingent product of historical chance events; it would also reveal how the underlying distribution of functions across sequence space shaped historical evolution. Here we combine ancestral protein reconstruction with deep mutational scanning to characterize alternative histories in the sequence space around an ancient transcription factor, which evolved a novel biological function through well-characterized mechanisms. We find hundreds of alternative protein sequences that use diverse biochemical mechanisms to perform the derived function at least as well as the historical outcome. These alternatives all require prior permissive substitutions that do not enhance the derived function, but not all require the same permissive changes that occurred during history. We find that if evolution had begun from a different starting point within the network of sequences encoding the ancestral function, outcomes with different genetic and biochemical forms would probably have resulted; this contingency arises from the distribution of functional variants in sequence space and epistasis between residues. Our results illuminate the topology of the vast space of possibilities from which history sampled one path, highlighting how the outcome of evolution depends on a serial chain of compounding chance events.

EvoChromo: towards a synthesis of chromatin biology and evolution.

Over the past few years, interest in chromatin and its evolution has grown. To further advance these interests, we organized a workshop with the support of The Company of Biologists to debate the current state of knowledge regarding the origin and evolution of chromatin. This workshop led to prospective views on the development of a new field of research that we term 'EvoChromo'. In this short Spotlight article, we define the breadth and expected impact of this new area of scientific inquiry on our understanding of both chromatin and evolution.

Fitness effects but no temperature-mediated balancing selection at the polymorphic Adh gene of Drosophila melanogaster.

Polymorphism in the alcohol dehydrogenase (ADH) protein of Drosophila melanogaster, like genetic variation in many other enzymes, has long been hypothesized to be maintained by a selective trade-off between thermostability and enzyme activity. Two major Adh variants, named Fast and Slow, are distributed along latitudinal clines on several continents. The balancing selection trade-off hypothesis posits that Fast is favored at high latitudes because it metabolizes alcohol faster, whereas Slow is favored at low latitudes because it is more stable at high temperatures. Here we use biochemical and physiological assays of precisely engineered genetic variants to directly test this hypothesis. As predicted, the Fast protein has higher catalytic activity than Slow, and both the Fast protein and regulatory variants linked to it confer greater ethanol tolerance on transgenic animals. But we found no evidence of a temperature-mediated trade-off: The Fast protein is not less stable or active at high temperatures, and Fast alleles increase ethanol tolerance and survivorship at all temperatures tested. Further, analysis of a population genomic dataset reveals no signature of balancing selection in the Adh gene. These results provide strong evidence against balancing selection driven by a stability/activity trade-off in Adh, and they justify caution about this hypothesis for other enzymes except those for which it has been directly tested. Our findings tentatively suggest that environment-specific selection for the Fast allele, coupled with demographic history, may have produced the observed pattern of Adh variation.

Family of neural wiring receptors in bilaterians defined by phylogenetic, biochemical, and structural evidence.

The evolution of complex nervous systems was accompanied by the expansion of numerous protein families, including cell-adhesion molecules, surface receptors, and their ligands. These proteins mediate axonal guidance, synapse targeting, and other neuronal wiring-related functions. Recently, 32 interacting cell surface proteins belonging to two newly defined families of the Ig superfamily (IgSF) in fruit flies were discovered to label different subsets of neurons in the brain and ventral nerve cord. They have been shown to be involved in synaptic targeting and morphogenesis, retrograde signaling, and neuronal survival. Here, we show that these proteins, Dprs and DIPs, are members of a widely distributed family of two- and three-Ig domain molecules with neuronal wiring functions, which we refer to as Wirins. Beginning from a single ancestral Wirin gene in the last common ancestor of Bilateria, numerous gene duplications produced the heterophilic Dprs and DIPs in protostomes, along with two other subfamilies that diversified independently across protostome phyla. In deuterostomes, the ancestral Wirin evolved into the IgLON subfamily of neuronal receptors. We show that IgLONs interact with each other and that their complexes can be broken by mutations designed using homology models based on Dpr and DIP structures. The nematode orthologs ZIG-8 and RIG-5 also form heterophilic and homophilic complexes, and crystal structures reveal numerous apparently ancestral features shared with Dpr-DIP complexes. The evolutionary, biochemical, and structural relationships we demonstrate here provide insights into neural development and the rise of the metazoan nervous system.

Pervasive contingency and entrenchment in a billion years of Hsp90 evolution.

Interactions among mutations within a protein have the potential to make molecular evolution contingent and irreversible, but the extent to which epistasis actually shaped historical evolutionary trajectories is unclear. To address this question, we experimentally measured how the fitness effects of historical sequence substitutions changed during the billion-year evolutionary history of the heat shock protein 90 (Hsp90) ATPase domain beginning from a deep eukaryotic ancestor to modern Saccharomyces cerevisiae We found a pervasive influence of epistasis. Of 98 derived amino acid states that evolved along this lineage, about half compromise fitness when introduced into the reconstructed ancestral Hsp90. And the vast majority of ancestral states reduce fitness when introduced into the extant S. cerevisiae Hsp90. Overall, more than 75% of historical substitutions were contingent on permissive substitutions that rendered the derived state nondeleterious, became entrenched by subsequent restrictive substitutions that made the ancestral state deleterious, or both. This epistasis was primarily caused by specific interactions among sites rather than a general effect on the protein's tolerance to mutation. Our results show that epistasis continually opened and closed windows of mutational opportunity over evolutionary timescales, producing histories and biological states that reflect the transient internal constraints imposed by the protein's fleeting sequence states.

Historical contingency and its biophysical basis in glucocorticoid receptor evolution.

Understanding how chance historical events shape evolutionary processes is a central goal of evolutionary biology. Direct insights into the extent and causes of evolutionary contingency have been limited to experimental systems, because it is difficult to know what happened in the deep past and to characterize other paths that evolution could have followed. Here we combine ancestral protein reconstruction, directed evolution and biophysical analysis to explore alternative 'might-have-been' trajectories during the ancient evolution of a novel protein function. We previously found that the evolution of cortisol specificity in the ancestral glucocorticoid receptor (GR) was contingent on permissive substitutions, which had no apparent effect on receptor function but were necessary for GR to tolerate the large-effect mutations that caused the shift in specificity. Here we show that alternative mutations that could have permitted the historical function-switching substitutions are extremely rare in the ensemble of genotypes accessible to the ancestral GR. In a library of thousands of variants of the ancestral protein, we recovered historical permissive substitutions but no alternative permissive genotypes. Using biophysical analysis, we found that permissive mutations must satisfy at least three physical requirements--they must stabilize specific local elements of the protein structure, maintain the correct energetic balance between functional conformations, and be compatible with the ancestral and derived structures--thus revealing why permissive mutations are rare. These findings demonstrate that GR evolution depended strongly on improbable, non-deterministic events, and this contingency arose from intrinsic biophysical properties of the protein.

Evolutionary biochemistry: revealing the historical and physical causes of protein properties.

The repertoire of proteins and nucleic acids in the living world is determined by evolution; their properties are determined by the laws of physics and chemistry. Explanations of these two kinds of causality - the purviews of evolutionary biology and biochemistry, respectively - are typically pursued in isolation, but many fundamental questions fall squarely at the interface of fields. Here we articulate the paradigm of evolutionary biochemistry, which aims to dissect the physical mechanisms and evolutionary processes by which biological molecules diversified and to reveal how their physical architecture facilitates and constrains their evolution. We show how an integration of evolution with biochemistry moves us towards a more complete understanding of why biological molecules have the properties that they do.

Distal substitutions drive divergent DNA specificity among paralogous transcription factors through subdivision of conformational space.

Many genomes contain families of paralogs--proteins with divergent function that evolved from a common ancestral gene after a duplication event. To understand how paralogous transcription factors evolve divergent DNA specificities, we examined how the glucocorticoid receptor and its paralogs evolved to bind activating response elements [(+)GREs] and negative glucocorticoid response elements (nGREs). We show that binding to nGREs is a property of the glucocorticoid receptor (GR) DNA-binding domain (DBD) not shared by other members of the steroid receptor family. Using phylogenetic, structural, biochemical, and molecular dynamics techniques, we show that the ancestral DBD from which GR and its paralogs evolved was capable of binding both nGRE and (+)GRE sequences because of the ancestral DBD's ability to assume multiple DNA-bound conformations. Subsequent amino acid substitutions in duplicated daughter genes selectively restricted protein conformational space, causing this dual DNA-binding specificity to be selectively enhanced in the GR lineage and lost in all others. Key substitutions that determined the receptors' response element-binding specificity were far from the proteins' DNA-binding interface and interacted epistatically to change the DBD's function through DNA-induced allosteric mechanisms. These amino acid substitutions subdivided both the conformational and functional space of the ancestral DBD among the present-day receptors, allowing a paralogous family of transcription factors to control disparate transcriptional programs despite high sequence identity.

Robustness of Reconstructed Ancestral Protein Functions to Statistical Uncertainty.

Hypotheses about the functions of ancient proteins and the effects of historical mutations on them are often tested using ancestral protein reconstruction (APR)-phylogenetic inference of ancestral sequences followed by synthesis and experimental characterization. Usually, some sequence sites are ambiguously reconstructed, with two or more statistically plausible states. The extent to which the inferred functions and mutational effects are robust to uncertainty about the ancestral sequence has not been studied systematically. To address this issue, we reconstructed ancestral proteins in three domain families that have different functions, architectures, and degrees of uncertainty; we then experimentally characterized the functional robustness of these proteins when uncertainty was incorporated using several approaches, including sampling amino acid states from the posterior distribution at each site and incorporating the alternative amino acid state at every ambiguous site in the sequence into a single "worst plausible case" protein. In every case, qualitative conclusions about the ancestral proteins' functions and the effects of key historical mutations were robust to sequence uncertainty, with similar functions observed even when scores of alternate amino acids were incorporated. There was some variation in quantitative descriptors of function among plausible sequences, suggesting that experimentally characterizing robustness is particularly important when quantitative estimates of ancient biochemical parameters are desired. The worst plausible case method appears to provide an efficient strategy for characterizing the functional robustness of ancestral proteins to large amounts of sequence uncertainty. Sampling from the posterior distribution sometimes produced artifactually nonfunctional proteins for sequences reconstructed with substantial ambiguity.

Evolution of increased complexity in a molecular machine.

Many cellular processes are carried out by molecular 'machines'-assemblies of multiple differentiated proteins that physically interact to execute biological functions. Despite much speculation, strong evidence of the mechanisms by which these assemblies evolved is lacking. Here we use ancestral gene resurrection and manipulative genetic experiments to determine how the complexity of an essential molecular machine--the hexameric transmembrane ring of the eukaryotic V-ATPase proton pump--increased hundreds of millions of years ago. We show that the ring of Fungi, which is composed of three paralogous proteins, evolved from a more ancient two-paralogue complex because of a gene duplication that was followed by loss in each daughter copy of specific interfaces by which it interacts with other ring proteins. These losses were complementary, so both copies became obligate components with restricted spatial roles in the complex. Reintroducing a single historical mutation from each paralogue lineage into the resurrected ancestral proteins is sufficient to recapitulate their asymmetric degeneration and trigger the requirement for the more elaborate three-component ring. Our experiments show that increased complexity in an essential molecular machine evolved because of simple, high-probability evolutionary processes, without the apparent evolution of novel functions. They point to a plausible mechanism for the evolution of complexity in other multi-paralogue protein complexes.

Thermodynamic system drift in protein evolution.

Proteins from thermophiles are generally more thermostable than their mesophilic homologs, but little is known about the evolutionary process driving these differences. Here we attempt to understand how the diverse thermostabilities of bacterial ribonuclease H1 (RNH) proteins evolved. RNH proteins from Thermus thermophilus (ttRNH) and Escherichia coli (ecRNH) share similar structures but differ in melting temperature (T(m)) by 20 °C. ttRNH's greater stability is caused in part by the presence of residual structure in the unfolded state, which results in a low heat capacity of unfolding (ΔCp) relative to ecRNH. We first characterized RNH proteins from a variety of extant bacteria and found that Tm correlates with the species' growth temperatures, consistent with environmental selection for stability. We then used ancestral sequence reconstruction to statistically infer evolutionary intermediates along lineages leading to ecRNH and ttRNH from their common ancestor, which existed approximately 3 billion years ago. Finally, we synthesized and experimentally characterized these intermediates. The shared ancestor has a melting temperature between those of ttRNH and ecRNH; the T(m)s of intermediate ancestors along the ttRNH lineage increased gradually over time, while the ecRNH lineage exhibited an abrupt drop in Tm followed by relatively little change. To determine whether the underlying mechanisms for thermostability correlate with the changes in T(m), we measured the thermodynamic basis for stabilization--ΔCp and other thermodynamic parameters--for each of the ancestors. We observed that, while the T(m) changes smoothly, the mechanistic basis for stability fluctuates over evolutionary time. Thus, even while overall stability appears to be strongly driven by selection, the proteins explored a wide variety of mechanisms of stabilization, a phenomenon we call "thermodynamic system drift." This suggests that even on lineages with strong selection to increase stability, proteins have wide latitude to explore sequence space, generating biophysical diversity and potentially opening new evolutionary pathways.

Vestigialization of an allosteric switch: genetic and structural mechanisms for the evolution of constitutive activity in a steroid hormone receptor.

An important goal in molecular evolution is to understand the genetic and physical mechanisms by which protein functions evolve and, in turn, to characterize how a protein's physical architecture influences its evolution. Here we dissect the mechanisms for an evolutionary shift in function in the mollusk ortholog of the steroid hormone receptors (SRs), a family of biologically essential transcription factors. In vertebrates, the activity of SRs allosterically depends on binding a hormonal ligand; in mollusks, however, the SR ortholog (called ER, because of high sequence similarity to vertebrate estrogen receptors) activates transcription in the absence of ligand and does not respond to steroid hormones. To understand how this shift in regulation evolved, we combined evolutionary, structural, and functional analyses. We first determined the X-ray crystal structure of the ER of the Pacific oyster Crassostrea gigas (CgER), and found that its ligand pocket is filled with bulky residues that prevent ligand occupancy. To understand the genetic basis for the evolution of mollusk ERs' unique functions, we resurrected an ancient SR progenitor and characterized the effect of historical amino acid replacements on its functions. We found that reintroducing just two ancient replacements from the lineage leading to mollusk ERs recapitulates the evolution of full constitutive activity and the loss of ligand activation. These substitutions stabilize interactions among key helices, causing the allosteric switch to become "stuck" in the active conformation and making activation independent of ligand binding. Subsequent changes filled the ligand pocket without further affecting activity; by degrading the allosteric switch, these substitutions vestigialized elements of the protein's architecture required for ligand regulation and made reversal to the ancestral function more complex. These findings show how the physical architecture of allostery enabled a few large-effect mutations to trigger a profound evolutionary change in the protein's function and shaped the genetics of evolutionary reversibility.

A second estrogen receptor from Japanese lamprey (Lethenteron japonicum) does not have activities for estrogen binding and transcription.

Estrogens regulate many physiological responses in vertebrates by binding to the estrogen receptor (ER), a ligand-activated transcription factor. To understand the evolution of vertebrate ERs and to investigate how estrogen acts in a jawless vertebrate, we used degenerate primer sets and PCR to isolate DNA fragments encoding two distinct ER subtypes, Esr1a and Esr1b from the Japanese lamprey, Lethenteron japonicum. Phylogenetic analysis indicates that these two ERs are the result of lineage-specific gene duplication within the jawless fishes, different from the previous duplication event of Esr1 (ERα) and Esr2 (ERß) within the jawed vertebrates. Reporter gene assays show that lamprey Esr1a displays both constitutive and estrogen-dependent activation of gene transcription. Domain swapping experiments indicate that constitutive activity resides in the A/B domain of lamprey Esr1a. Unexpectedly, lamprey Esr1b does not bind estradiol and is not stimulated by other estrogens, androgens or corticosteroids. A 3D model of lamprey Esr1b suggests that although estradiol fits into the steroid binding site, some stabilizing contacts between the ligand and side chains that are found in human Esr1 and Esr2 are missing in lamprey Esr1b.

Biophysical mechanisms for large-effect mutations in the evolution of steroid hormone receptors.

The genetic and biophysical mechanisms by which new protein functions evolve is a central question in evolutionary biology, biochemistry, and biophysics. Of particular interest is whether major shifts in protein function are caused by a few mutations of large effect and, if they are, the mechanisms that mediate these changes. Here we combine ancestral protein reconstruction with genetic manipulation and explicit studies of protein structure and dynamics to dissect an ancient and discrete shift in ligand specificity in the steroid receptors, a family of biologically essential hormone-controlled transcription factors. We previously found that the ancestor of the entire steroid receptor family was highly specific for estrogens, but its immediate phylogenetic descendant was sensitive only to androgens, progestogens, and corticosteroids. Here we show that this shift in function was driven primarily by two historical amino acid changes, which caused a ∼70,000-fold change in the ancestral protein's specificity. These replacements subtly changed the chemistry of two amino acids, but they dramatically reduced estrogen sensitivity by introducing an excess of interaction partners into the receptor/estrogen complex, inducing a frustrated ensemble of suboptimal hydrogen bond networks unique to estrogens. This work shows how the protein's architecture and dynamics shaped its evolution, amplifying a few biochemically subtle mutations into major shifts in the energetics and function of the protein.

An epistatic ratchet constrains the direction of glucocorticoid receptor evolution.

The extent to which evolution is reversible has long fascinated biologists. Most previous work on the reversibility of morphological and life-history evolution has been indecisive, because of uncertainty and bias in the methods used to infer ancestral states for such characters. Further, despite theoretical work on the factors that could contribute to irreversibility, there is little empirical evidence on its causes, because sufficient understanding of the mechanistic basis for the evolution of new or ancestral phenotypes is seldom available. By studying the reversibility of evolutionary changes in protein structure and function, these limitations can be overcome. Here we show, using the evolution of hormone specificity in the vertebrate glucocorticoid receptor as a case-study, that the evolutionary path by which this protein acquired its new function soon became inaccessible to reverse exploration. Using ancestral gene reconstruction, protein engineering and X-ray crystallography, we demonstrate that five subsequent 'restrictive' mutations, which optimized the new specificity of the glucocorticoid receptor, also destabilized elements of the protein structure that were required to support the ancestral conformation. Unless these ratchet-like epistatic substitutions are restored to their ancestral states, reversing the key function-switching mutations yields a non-functional protein. Reversing the restrictive substitutions first, however, does nothing to enhance the ancestral function. Our findings indicate that even if selection for the ancestral function were imposed, direct reversal would be extremely unlikely, suggesting an important role for historical contingency in protein evolution.