Ancestral Sequence Research Articles

The species coalescent indicates possible bat and pangolin origins of the COVID-19 pandemic

A consensus species tree is reconstructed from 11 gene trees for human, bat, and pangolin beta coronaviruses from samples taken early in the pandemic (prior to April 1, 2020). Using coalescent theory, the shallow (short branches relative to the hosts) consensus species tree provides evidence of recent gene flow events between bat and pangolin beta coronaviruses predating the zoonotic transfer to humans. The consensus species tree was also used to reconstruct the ancestral sequence of human SARS-CoV-2, which was 2 nucleotides different from the Wuhan sequence. The time to most recent common ancestor was estimated to be Dec 8, 2019 with a bat origin. Some human, bat, and pangolin coronavirus lineages found in China are phylogenetically distinct, a rare example of a class II phylogeography pattern (Avise et al. in Ann Rev Eco Syst 18:489–422, 1987). The consensus species tree is a product of evolutionary factors, providing evidence of repeated zoonotic transfers between bat and pangolin as a reservoir for future zoonotic transfers to humans.

Deciphering the evolution of flavin-dependent monooxygenase stereoselectivity using ancestral sequence reconstruction

Controlling the selectivity of a reaction is critical for target-oriented synthesis. Accessing complementary selectivity profiles enables divergent synthetic strategies, but is challenging to achieve in biocatalytic reactions given enzymes' innate preferences of a single selectivity. Thus, it is critical to understand the structural features that control selectivity in biocatalytic reactions to achieve tunable selectivity. Here, we investigate the structural features that control the stereoselectivity in an oxidative dearomatization reaction that is key to making azaphilone natural products. Crystal structures of enantiocomplementary biocatalysts guided the development of multiple hypotheses centered on the structural features that control the stereochemical outcome of the reaction; however, in many cases, direct substitutions of active site residues in natural proteins led to inactive enzymes. Ancestral sequence reconstruction(ASR) and resurrection were employed as an alternative strategy to probe the impact of each residue on the stereochemical outcome of the dearomatization reaction. These studies suggest that two mechanisms are active in controlling the stereochemical outcome of the oxidative dearomatization reaction: one involving multiple active site residues in AzaH and the other dominated by a single Phe to Tyr switch in TropB and AfoD. Moreover, this study suggests that the flavin-dependent monooxygenases (FDMOs) adopt simple and flexible strategies to control stereoselectivity, which has led to stereocomplementary azaphilone natural products produced by fungi. This paradigm of combining ASR and resurrection with mutational and computational studies showcases sets of tools for understanding enzyme mechanisms and provides a solid foundation for future protein engineering efforts.

Fortuitously compatible protein surfaces primed allosteric control in cyanobacterial photoprotection

Highly specific interactions between proteins are a fundamental prerequisite for life, but how they evolve remains an unsolved problem. In particular, interactions between initially unrelated proteins require that they evolve matching surfaces. It is unclear whether such surface compatibilities can only be built by selection in small incremental steps, or whether they can also emerge fortuitously. Here, we used molecular phylogenetics, ancestral sequence reconstruction and biophysical characterization of resurrected proteins to retrace the evolution of an allosteric interaction between two proteins that act in the cyanobacterial photoprotection system. We show that this interaction between the orange carotenoid protein (OCP) and its unrelated regulator, the fluorescence recovery protein (FRP), evolved when a precursor of FRP was horizontally acquired by cyanobacteria. FRP’s precursors could already interact with and regulate OCP even before these proteins first encountered each other in an ancestral cyanobacterium. The OCP–FRP interaction exploits an ancient dimer interface in OCP, which also predates the recruitment of FRP into the photoprotection system. Together, our work shows how evolution can fashion complex regulatory systems easily out of pre-existing components.

Immunoanalytical Approach for Detecting and Identifying Ancestral Peptide Biomarkers in Early Earth Analogue Environments.

Several mass spectrometry and spectroscopic techniques have been used in the search for molecular biomarkers on Mars. A major constraint is their capability to detect and identify large and complex compounds such as peptides or other biopolymers. Multiplex immunoassays can detect these compounds, but antibodies must be produced for a large number of sequence-dependent molecular targets. Ancestral Sequence Reconstruction (ASR) followed by protein "resurrection" in the lab can help to narrow the selection of targets. Herein, we propose an immunoanalytical method to identify ancient and universally conserved protein/peptide sequences as targets for identifying ancestral biomarkers in nature. We have developed, tested, and validated this approach by producing antibodies to eight previously described ancestral resurrected proteins (three β-lactamases, three thioredoxins, one Elongation Factor Tu, and one RuBisCO, all of them theoretically dated as Precambrian), and used them as a proxy to search for any potential feature of them that could be present in current natural environments. By fluorescent sandwich microarray immunoassays (FSMI), we have detected positive immunoreactions with antibodies to the oldest β-lactamase and thioredoxin proteins (ca. 4 Ga) in samples from a hydrothermal environment. Fine epitope mapping and inhibitory immunoassays allowed the identification of well-conserved epitope peptide sequences that resulted from ASR and were present in the sample. We corroborated these results by metagenomic sequencing and found several genes encoding analogue proteins with significant matches to the peptide epitopes identified with the antibodies. The results demonstrated that peptides inferred from ASR studies have true counterpart analogues in Nature, which validates and strengthens the well-known ASR/protein resurrection technique and our immunoanalytical approach for investigating ancient environments and metabolisms on Earth and elsewhere.

Mutation-driven parallel evolution in emergence of ACE2-utilizing sarbecoviruses.

Mutation and recombination are two major genetic mechanisms that drive the evolution of viruses. They both exert an interplay during virus evolution, in which mutations provide a first ancestral source of genetic diversity for subsequent recombination. Sarbecoviruses are a group of evolutionarily related β-coronaviruses including human severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 and a trove of related animal viruses called SARS-like CoVs (SL-CoVs). This group of members either use or not use angiotensin-converting enzyme 2 (ACE2) as their entry receptor, which has been linked to the properties of their spike protein receptor binding domains (RBDs). This raises an outstanding question regarding how ACE2 binding originated within sarbecoviruses. Using a combination of analyses of phylogenies, ancestral sequences, structures, functions and molecular dynamics, we provide evidence in favor of an evolutionary scenario, in which three distinct ancestral RBDs independently developed the ACE2 binding trait via parallel amino acid mutations. In this process, evolutionary intermediate RBDs might be firstly formed through loop extensions to offer key functional residues accompanying point mutations to remove energetically unfavorable interactions and to change the dynamics of the functional loops, all required for ACE2 binding. Subsequent optimization in the context of evolutionary intermediates led to the independent emergence of ACE2-binding RBDs in the SARS-CoV and SARS-CoV-2 clades of Asian origin and the clade comprising SL-CoVs of European and African descent. These findings will help enhance our understanding of mutation-driven evolution of sarbecoviruses in their early history.

Evidence for an early green/red photocycle that precedes the diversification of GAF domain photoreceptor cyanobacteriochromes

Phytochromes are linear tetrapyrrole-binding photoreceptors in eukaryotes and bacteria, primarily responding to red and far-red light signals reversibly. Among the GAF domain-based phytochrome superfamily, cyanobacteria-specific cyanobacteriochromes show various optical properties covering the entire visible region. It is unknown what physiological demands drove the evolution of cyanobacteriochromes in cyanobacteria. Here, we utilize ancestral sequence reconstruction and biochemical verification to show that the resurrected ancestral cyanobacteriochrome proteins reversibly respond to green and red light signals. pH titration analyses indicate that the deprotonation of the bound phycocyanobilin chromophore is crucial to perceive green light. The ancestral cyanobacteriochromes show only modest thermal reversion to the green light-absorbing form, suggesting that they evolved to sense the incident green/red light ratio. Many cyanobacteria can utilize green light for photosynthesis using phycobilisome light-harvesting complexes. The green/red sensing cyanobacteriochromes may have allowed better acclimation to changing light environments by rearranging the absorption capacity of the phycobilisome through chromatic acclimation.Graphical abstract

Ancestral Sequence Reconstruction Enhances Gene Mining Efforts for Industrial Ene Reductases by Expanding Enzyme Panels with Thermostable Catalysts

Ancestral Sequence Reconstruction Enhances Gene Mining Efforts for Industrial Ene Reductases by Expanding Enzyme Panels with Thermostable Catalysts

Big data mining, rational modification, and ancestral sequence reconstruction inferred multiple xylose isomerases for biorefinery.

The isomerization of xylose to xylulose is considered the most promising approach to initiate xylose bioconversion. Here, phylogeny-guided big data mining, rational modification, and ancestral sequence reconstruction strategies were implemented to explore new active xylose isomerases (XIs) for Saccharomyces cerevisiae. Significantly, 13 new active XIs for S. cerevisiae were mined or artificially created. Moreover, the importance of the amino-terminal fragment for maintaining basic XI activity was demonstrated. With the mined XIs, four efficient xylose-utilizing S. cerevisiae were constructed and evolved, among which the strain S. cerevisiae CRD5HS contributed to ethanol titers as high as 85.95 and 94.76 g/liter from pretreated corn stover and corn cob, respectively, without detoxifying or washing pretreated biomass. Potential genetic targets obtained from adaptive laboratory evolution were further analyzed by sequencing the high-performance strains. The combined XI mining methods described here provide practical references for mining other scarce and valuable enzymes.

Biophysical models for studying patterns of RuBisCO small subunit evolution.

For much of Earth's history, key enzymes have dominated the conversion of energy into the production of biomass. This history extends hundreds of millions of years into the past, but the available record provides few direct indicators of the emergence and evolution of enzymes that shaped our planet. A good example is Ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO), which is responsible for the fixation of CO2 in conjunction with oxygenic photosynthesis and is estimated to be the most abundant enzyme on Earth. RuBisCO form-I is the only form composed of a large and small subunit. The biochemical role that the small subunit plays is unclear, and the circumstances of its coevolution with the large subunit or with drastic changes to the Earth's surface environment such as the concentrations of CO2 and O2 gases are also unknown. Here, we explore how the small subunit-mediated RuBisCO form-I transition by generating a phylogenetic tree of proteins, reconstructing its ancestral sequences, and analyzing the structure as well as function of the ancestral large-small subunit RuBisCO complexes. From ancient and extant structures, we identify the general attributes that define the interaction between the two RuBisCO subunits and inspect how these interactions have evolved with reference to the ancestral large subunit RuBisCO forms and the partition and diffusion of CO2 and O2 at different relative concentrations of the two gasses. Our results are consistent with the possibility that the small subunit was an evolutionary response to the global rise of atmospheric O2 over 2 billion years ago, and that it impacted large subunit CO2 diffusion. Our research highlights how a protein subunit that emerged billions of years ago in response to a change in the planet's atmosphere impacts the function of an essential protein complex.

Experimental reconstruction of membrane transporter evolution.

Cellular homeostasis relies on solutes crossing cell membranes against their concentration gradients. Membrane transporters catalyze this energetically unfavorable movement using energy from the downhill movement of ions. Often, related transporters couple to either Na+ or H+ gradients; it can be difficult to glean the mechanistic reasons between such differences because the ability to use ions with such different physicochemical properties might require many changes in the transporter. Since evolutionary adaptations to environmental constraints such as salt concentrations and pH likely gave rise to the switch between sodium- and proton-motive forces, we hypothesized that phylogenetic reconstruction of ion-coupling may be a method to study such complex changes. We used phylogenetic analysis to reconstruct the evolution of prokaryotic glutamate transporters and deduce stepwise amino acid sequence changes resulting in the emergence of H+-coupled transporters from Na+-coupled ancestors. The analysis suggests that in addition to changing protonatable residues, introducing H+-dependence increases the hydrophobicity of previously polar regions proximal to the substrate and Na+ binding sites. To contextualize these changes, we used ancestral sequence reconstruction (ASR) to infer transporter sequences before and after the emergence of H+-coupling. We expressed and purified the inferred ancestral transporters, enabling their biochemical characterization and high-resolution cryo-EM structures. Surprisingly, the ancestral transporter before the emergence of H+-coupling shows Na+-independent substrate binding even though its Na+-binding pockets are structurally near-identical to Na+-coupled homologs. The emergence of H+-coupling preserves the protein architecture while subtly restructuring these critical regions. Our results suggest that changes in local structure and complex global allosteric properties underlie the diversification of the ion-motive forces used within the glutamate transporter family. The approach we developed might be broadly applicable to dissecting complex allosteric properties of membrane proteins.

Understanding the evolution of cold sensing and menthol activation in TRPM8.

Transient receptor potential melastatin 8 (TRPM8) is a polymodally gated ion channel activated by a variety of stimuli including menthol, cold, and pH. These chemical and physical stimuli ultimately regulate the ion conductance pore. While the menthol binding site and ion conduction pathway gate have been identified, the physics that underlies the mechanism of activation among different stimuli remains elusive. Particularly, the temperature modulation, how TRPM8 converts small changes in temperature into a conformational change, ultimately resulting in signal transduction is not well understood. Here to shed light into polymodal regulation of TRPM8, we perform a comparative analysis of conformational dynamics and functional characterizations of ancestral and extant human TRPM8. We use ancestral sequence reconstruction, a computational molecular evolution tool, to explore that ability to disentangle TRPM8 activation modes, which presumably emerged distinctly over evolutionary time. Based on the extant sequences of TRPM8, ancestral reconstruction was used to determine the DNA sequences of the last common primate, mammalian, and vertebrate ancestors. Resurrection and experimental analysis of these ancestral TRPM8 proteins show them to be functional in heterologous mammalian cells. Whole-cell patch-clamp electrophysiology and calcium mobilization assays confirm our hypothesis that the cold and menthol responses differ in terms of emergence and sensitivity between the ancestral constructs suggesting a new method to dissect polymodal protein behavior. Furthermore, comparative conformational analysis of ancestral and extant TRMP8 through the Dynamic Flexibility and Dynamic Coupling analysis also further support this hypothesis and provide dynamical finger prints for menthol and cold sensation.

Novavax NVX-COV2373 triggers neutralization of Omicron sub-lineages

The SARS-CoV-2 Omicron (B.1.1.529) Variant of Concern (VOC) and its sub-lineages (including BA.2, BA.4, BA.5, BA.2.12.1) contain spike mutations that confer high level resistance to neutralizing antibodies induced by vaccination with ancestral spike or infection with previously circulating variants. The NVX-CoV2373 vaccine, a protein nanoparticle vaccine containing the ancestral spike sequence, has value in countries with constrained cold-chain requirements. Here we report neutralizing titers following two or three doses of NVX-CoV2373. We show that after two doses, Omicron sub-lineages BA.1 and BA.4/BA.5 were resistant to neutralization by 72% (21/29) and 59% (17/29) of samples respectively. However, after a third dose of NVX-CoV2373, we observed high titers against Omicron BA.1 (GMT: 1,197) and BA.4/BA.5 (GMT: 582), with responses similar in magnitude to those triggered by three doses of an mRNA vaccine. These data are of particular relevance as BA.4/BA.5 is dominating in multiple locations, and highlight the potential utility of the NVX-CoV2373 vaccine as a booster in resource-limited environments.

Topiary: Pruning the manual labor from ancestral sequence reconstruction.

Ancestral sequence reconstruction (ASR) is a powerful tool to study the evolution of proteins and thus gain deep insight into the relationships among protein sequence, structure, and function. A major barrier to its broad use is the complexity of the task: it requires multiple software packages, complex file manipulations, and expert phylogenetic knowledge. Here we introduce topiary, a software pipeline that aims to overcome this barrier. To use topiary, users prepare a spreadsheet with a handful of sequences. Topiary then: (1) Infers the taxonomic scope for the ASR study and finds relevant sequences by BLAST; (2) Does taxonomically informed sequence quality control and redundancy reduction; (3) Constructs a multiple sequence alignment; (4) Generates a maximum-likelihood gene tree; (5) Reconciles the gene tree to the species tree; (6) Reconstructs ancestral amino acid sequences; and (7) Determines branch supports. The pipeline returns annotated evolutionary trees, spreadsheets with sequences, and graphical summaries of ancestor quality. This is achieved by integrating modern phylogenetics software (Muscle5, RAxML-NG, GeneRax, and PastML) with online databases (NCBI and the Open Tree of Life). In this paper, we introduce non-expert readers to the steps required for ASR, describe the specific design choices made in topiary, provide a detailed protocol for users, and then validate the pipeline using datasets from a broad collection of protein families. Topiary is freely available for download: https://github.com/harmslab/topiary.

Reconstruction of hundreds of reference ancestral genomes across the eukaryotic kingdom

Ancestral sequence reconstruction is a fundamental aspect of molecular evolution studies and can trace small-scale sequence modifications through the evolution of genomes and species. In contrast, fine-grained reconstructions of ancestral genome organizations are still in their infancy, limiting our ability to draw comprehensive views of genome and karyotype evolution. Here we reconstruct the detailed gene contents and organizations of 624 ancestral vertebrate, plant, fungi, metazoan and protist genomes, 183 of which are near-complete chromosomal gene order reconstructions. Reconstructed ancestral genomes are similar to their descendants in terms of gene content as expected and agree precisely with reference cytogenetic and in silico reconstructions when available. By comparing successive ancestral genomes along the phylogenetic tree, we estimate the intra- and interchromosomal rearrangement history of all major vertebrate clades at high resolution. This freely available resource introduces the possibility to follow evolutionary processes at genomic scales in chronological order, across multiple clades and without relying on a single extant species as reference.

Investigation of the Molecular Evolution of Treg Suppression Mechanisms Indicates a Convergent Origin.

Regulatory T cell (Treg) suppression of conventional T cells is a central mechanism that ensures immune system homeostasis. The exact time point of Treg emergence is still disputed. Furthermore, the time of Treg-mediated suppression mechanisms’ emergence has not been identified. It is not yet known whether Treg suppression mechanisms diverged from a single pathway or converged from several sources. We investigated the evolutionary history of Treg suppression pathways using various phylogenetic analysis tools. To ensure the conservation of function for investigated proteins, we augmented our study using nonhomology-based methods to predict protein functions among various investigated species and mined the literature for experimental evidence of functional convergence. Our results indicate that a minority of Treg suppressor mechanisms could be homologs of ancient conserved pathways. For example, CD73, an enzymatic pathway known to play an essential role in invertebrates, is highly conserved between invertebrates and vertebrates, with no evidence of positive selection (w = 0.48, p-value < 0.00001). Our findings indicate that Tregs utilize homologs of proteins that diverged in early vertebrates. However, our findings do not exclude the possibility of a more evolutionary pattern following the duplication degeneration−complementation (DDC) model. Ancestral sequence reconstruction showed that Treg suppression mechanism proteins do not belong to one family; rather, their emergence seems to follow a convergent evolutionary pattern.

Ancestral reconstruction of the MotA stator subunit reveals that conserved residues far from the pore are required to drive flagellar motility.

The bacterial flagellar motor (BFM) is a rotary nanomachine powered by the translocation of ions across the inner membrane through the stator complex. The stator complex consists of two membrane proteins: MotA and MotB (in H+-powered motors), or PomA and PomB (in Na+-powered motors). In this study, we used ancestral sequence reconstruction (ASR) to probe which residues of MotA correlate with function and may have been conserved to preserve motor function. We reconstructed 10 ancestral sequences of MotA and found four of them were motile in combination with contemporary Escherichia coli MotB and in combination with our previously published functional ancestral MotBs. Sequence comparison between wild-type (WT) E. coli MotA and MotA-ASRs revealed 30 critical residues across multiple domains of MotA that were conserved among all motile stator units. These conserved residues included pore-facing, cytoplasm-facing, and MotA-MotA intermolecular facing sites. Overall, this work demonstrates the role of ASR in assessing conserved variable residues in a subunit of a molecular complex.

Abstract 2623: Extant sequence reconstruction is a new method for evaluating and improving ancestral sequence reconstructions to limit the butterfly effect

Abstract 2623: Extant sequence reconstruction is a new method for evaluating and improving ancestral sequence reconstructions to limit the butterfly effect

Experimental and computational analysis of the ancestry of an evolutionary young enzyme from histidine biosynthesis.

The conservation of fold and chemistry of the enzymes associated with histidine biosynthesis suggests that this pathway evolved prior to the diversification of Bacteria, Archaea, and Eukaryotes. The only exception is the histidinol phosphate phosphatase (HolPase). So far, non-homologous HolPases that possess distinct folds and belong to three different protein superfamilies have been identified in various phylogenetic clades. However, their evolution has remained unknown to date. Here, we analyzed the evolutionary history of the HolPase from γ-Proteobacteria (HisB-N). It has been argued that HisB-N and its closest homologue d-glycero-d-manno-heptose-1,7-bisphosphate 7-phosphatase (GmhB) have emerged from the same promiscuous ancestral phosphatase. GmhB variants catalyze the hydrolysis of the anomeric d-glycero-d-manno-heptose-1,7-bisphosphate (αHBP or βHBP) with a strong preference for one anomer (αGmhB or βGmhB). We found that HisB-N from Escherichia coli shows promiscuous activity for βHBP but not αHBP, while βGmhB from Crassaminicella sp. shows promiscuous activity for HolP. Accordingly, a combined phylogenetic tree of αGmhBs, βGmhBs, and HisB-N sequences revealed that HisB-Ns form a compact subcluster derived from βGmhBs. Ancestral sequence reconstruction and in vitro analysis revealed a promiscuous HolPase activity in the resurrected enzymes prior to functional divergence of the successors. The following increase in catalytic efficiency of the HolP turnover is reflected in the shape and electrostatics of the active site predicted by AlphaFold. An analysis of the phylogenetic tree led to a revised evolutionary model that proposes the horizontal gene transfer of a promiscuous βGmhB from δ- to γ-Proteobacteria where it evolved to the modern HisB-N.

Editorial Retraction.

In the Research Article “Gradual emergence followed by exponential spread of the SARS-CoV-2 Omicron variant in Africa”(1), we reported data from retrospective characterization of viral genomes of putative ancestors of the SARS-CoV-2 Omicron variant from western Africa months before the first detection of Omicron. After several social media posts suggested that these putative early Omicron ancestor sequences may have been false positives, we reanalyzed our data and the residual samples. We found a mixture of different SARS-CoV-2 genomic fragments contaminating some of the samples and sequence data on which we based our analysis. The residual samples are now exhausted, and the reconstruction of evolutionary intermediates cannot be replicated. Therefore, we are retracting our Research Article. The epidemiological data are not called into question and will be made available.

Engineering the Thermostability of a d-Carbamoylase Based on Ancestral Sequence Reconstruction for the Efficient Synthesis of d-Tryptophan.

Employing ancestral sequence reconstruction and consensus sequence analysis, the thermostability of a novel d-carbamoylase derived from Nitratireductor indicus (NiHyuC) was engineered through greedy-oriented iterative combinatorial mutagenesis. A mutant S202P/E208D/R277L (M4Th3) was obtained with significantly elevated thermostability. M4Th3 has a half-life of 36.5 h at 40 °C, about 28.5 times of 1.3 h of its parent M4. For the reaction at 40 °C, M4Th3 can catalyze 10 mM N-carbamoyl-d-tryptophan to produce d-tryptophan with a conversion ratio of 96.4% after 12 h, which is significantly higher than 64.1% of M4. MD simulation reveals that new hydrogen bonds emerging from E208D on the surface can increase the hydrophobicity of the protein, leading to improved stability. More importantly, R277L could contribute to enhanced interface stability of homodimeric M4. This study provides a thermostable d-carbamoylase for the "hydantoinase process", which has potential in the industrial synthesis of optically pure natural and non-natural amino acids.