Parallel adaptations to high temperatures in the Archaean eon

Phylogenomic analyses of ochrophytes (stramenopiles) with an emphasis on neglected lineages

Ochrophyta is a photosynthetic lineage that crowns the phylogenetic tree of stramenopiles, one of the major eukaryotic supergroups. Due to their ecological impact as a major primary producer, ochrophytes are relatively well-studied compared to the rest of the stramenopiles, yet their evolutionary relationships remain poorly understood. This is in part due to a number of missing lineages in large-scale multigene analyses, and an apparently rapid radiation leading to many short internodes between ochrophyte subgroups in the tree. These short internodes are also found across deep-branching lineages of stramenopiles with limited phylogenetic signal, leaving many relationships controversial overall. We have addressed this issue with other deep-branching stramenopiles recently, and now examine whether contentious relationships within the ochrophytes may be resolved with the help of filling in missing lineages in an updated phylogenomic dataset of ochrophytes, along with exploring various gene filtering criteria to identify the most phylogenetically informative genes. We generated ten new transcriptomes from various culture collections and a single-cell isolation from an environmental sample, added these to an existing phylogenomic dataset, and examined the effects of selecting genes with high phylogenetic signal or low phylogenetic noise. For some previously contentious relationships, we find a variety of analyses and gene filtering criteria consistently unite previously unstable groupings with strong statistical support. For example, we recovered a robust grouping of Eustigmatophyceae with Raphidophyceae-Phaeophyceae-Xanthophyceae while Olisthodiscophyceae formed a sister-lineage to Pinguiophyceae. Selecting genes with high phylogenetic signal or data quality recovered more stable topologies. Overall, we find that adding under-represented groups across different lineages is still crucial in resolving phylogenetic relationships, and discrete gene properties affect lineages of stramenopiles differently. This is something which may be explored to further our understanding of the molecular evolution of stramenopiles.

Structure and Activity of Reconstructed Pseudo-Ancestral Cyclotides

Cyclotides are cyclic peptides that are promising scaffolds for the design of drug candidates and chemical tools. However, despite there being hundreds of reported cyclotides, drug design studies have commonly focussed on a select few prototypic examples. Here, we explored whether ancestral sequence reconstruction could be used to generate new cyclotides for further optimization. We show that the reconstructed 'pseudo-ancestral' sequences, named Ancy-m (for the ancestral cyclotide of the Möbius sub-family) and Ancy-b (for the bracelet sub-family), have well-defined structures like their extant members, comprising the core structural feature of a cyclic cystine knot. This motif underpins efforts to re-engineer cyclotides for agrochemical and therapeutic applications. We further show that the reconstructed sequences are resistant to temperatures approaching boiling, bind to phosphatidyl-ethanolamine lipid bilayers at micromolar affinity, and inhibit the growth of insect cells at inhibitory concentrations in the micromolar range. Interestingly, the Ancy-b cyclotide had a higher oxidative folding yield than its comparator cyclotide cyO2, which belongs to the bracelet cyclotide subfamily known to be notoriously difficult to fold. Overall, this study provides new cyclotide sequences not yet found naturally that could be valuable starting points for the understanding of cyclotide evolution and for further optimization as drug leads.

The nature of the last universal common ancestor and its impact on the early Earth system

The nature of the last universal common ancestor (LUCA), its age and its impact on the Earth system have been the subject of vigorous debate across diverse disciplines, often based on disparate data and methods. Age estimates for LUCA are usually based on the fossil record, varying with every reinterpretation. The nature of LUCA's metabolism has proven equally contentious, with some attributing all core metabolisms to LUCA, whereas others reconstruct a simpler life form dependent on geochemistry. Here we infer that LUCA lived ~4.2 Ga (4.09-4.33 Ga) through divergence time analysis of pre-LUCA gene duplicates, calibrated using microbial fossils and isotope records under a new cross-bracing implementation. Phylogenetic reconciliation suggests that LUCA had a genome of at least 2.5 Mb (2.49-2.99 Mb), encoding around 2,600 proteins, comparable to modern prokaryotes. Our results suggest LUCA was a prokaryote-grade anaerobic acetogen that possessed an early immune system. Although LUCA is sometimes perceived as living in isolation, we infer LUCA to have been part of an established ecological system. The metabolism of LUCA would have provided a niche for other microbial community members and hydrogen recycling by atmospheric photochemistry could have supported a modestly productive early ecosystem.

Persistent homology reveals strong phylogenetic signal in 3D protein structures

Changes that occur in proteins over time provide a phylogenetic signal that can be used to decipher their evolutionary history and the relationships between organisms. Sequence comparison is the most common way to access this phylogenetic signal, while those based on 3D structure comparisons are still in their infancy. In this study, we propose an effective approach based on Persistent Homology Theory (PH) to extract the phylogenetic information contained in protein structures. PH provides efficient and robust algorithms for extracting and comparing geometric features from noisy datasets at different spatial resolutions. PH has a growing number of applications in the life sciences, including the study of proteins (e.g. classification, folding). However, it has never been used to study the phylogenetic signal they may contain. Here, using 518 protein families, representing 22,940 protein sequences and structures, from 10 major taxonomic groups, we show that distances calculated with PH from protein structures correlate strongly with phylogenetic distances calculated from protein sequences, at both small and large evolutionary scales. We test several methods for calculating PH distances and propose some refinements to improve their relevance for addressing evolutionary questions. This work opens up new perspectives in evolutionary biology by proposing an efficient way to access the phylogenetic signal contained in protein structures, as well as future developments of topological analysis in the life sciences.

Evolution of Consciousness

The natural evolution of consciousness in different animal species mandates that conscious experiences are causally potent in order to confer any advantage in the struggle for survival. Any endeavor to construct a physical theory of consciousness based on emergence within the framework of classical physics, however, leads to causally impotent conscious experiences in direct contradiction to evolutionary theory since epiphenomenal consciousness cannot evolve through natural selection. Here, we review recent theoretical advances in describing sentience and free will as fundamental aspects of reality granted by quantum physical laws. Modern quantum information theory considers quantum states as a physical resource that endows quantum systems with the capacity to perform physical tasks that are classically impossible. Reductive identification of conscious experiences with the quantum information comprised in quantum brain states allows for causally potent consciousness that is capable of performing genuine choices for future courses of physical action. The consequent evolution of brain cortical networks contributes to increased computational power, memory capacity, and cognitive intelligence of the living organisms.

A moderately thermophilic origin of a novel family of marine group II euryarchaeota from deep ocean

Marine group II (MGII) is the most abundant planktonic heterotrophic archaea in the ocean. The evolutionary history of MGII archaea is elusive. In this study, 13 new MGII metagenome-assembled genomes were recovered from surface to the hadal zone in Challenger Deep of the Mariana Trench; four of them from the deep ocean represent a novel group. The optimal growth temperature (OGT) of the common ancestor of MGII has been estimated to be at about 60°C and OGTs of MGIIc, MGIIb, and MGIIa at 47°C-50ºC, 37°C-44ºC, and 30°C-37ºC, respectively, suggesting the adaptation of these species to different temperatures during evolution. The estimated OGT range of MGIIc was supported by experimental measurements of cloned β-galactosidase that showed optimal enzyme activity around 50°C. These results indicate that MGIIc may have originated from a common ancestor that lived in warm or even hot marine environment, such as hydrothermal vents.

Early Divergence and Gene Exchange Highways in the Evolutionary History of Mesoaciditogales

The placement of a nonhyperthermophilic order Mesoaciditogales as the earliest branching clade within the Thermotogota phylum challenges the prevailing hypothesis that the last common ancestor of Thermotogota was a hyperthermophile. Yet, given the long branch leading to the only two Mesoaciditogales described to date, the phylogenetic position of the order may be due to the long branch attraction artifact. By testing various models and applying data recoding in phylogenetic reconstructions, we observed that early branching of Mesoaciditogales within Thermotogota is strongly supported by the conserved marker genes assumed to be vertically inherited. However, based on the taxonomic content of 1,181 gene families and a phylogenetic analysis of 721 gene family trees, we also found that a substantial number of Mesoaciditogales genes are more closely related to species from the order Petrotogales. These genes contribute to coenzyme transport and metabolism, fatty acid biosynthesis, genes known to respond to heat and cold stressors, and include many genes of unknown functions. The Petrotogales comprise moderately thermophilic and mesophilic species with similar temperature tolerances to that of Mesoaciditogales. Our findings hint at extensive horizontal gene transfer (HGT) between, or parallel independent gene gains by, the two ecologically similar lineages and suggest that the exchanged genes may be important for adaptation to comparable temperature niches.

Compositionally Constrained Sites Drive Long-Branch Attraction

Accurate phylogenies are fundamental to our understanding of the pattern and process of evolution. Yet, phylogenies at deep evolutionary timescales, with correspondingly long branches, have been fraught with controversy resulting from conflicting estimates from models with varying complexity and goodness of fit. Analyses of historical as well as current empirical datasets, such as alignments including Microsporidia, Nematoda, or Platyhelminthes, have demonstrated that inadequate modeling of across-site compositional heterogeneity, which is the result of biochemical constraints that lead to varying patterns of accepted amino acids along sequences, can lead to erroneous topologies that are strongly supported. Unfortunately, models that adequately account for across-site compositional heterogeneity remain computationally challenging or intractable for an increasing fraction of contemporary datasets. Here, we introduce "compositional constraint analysis," a method to investigate the effect of site-specific constraints on amino acid composition on phylogenetic inference. We show that more constrained sites with lower diversity and less constrained sites with higher diversity exhibit ostensibly conflicting signals under models ignoring across-site compositional heterogeneity that lead to long-branch attraction artifacts and demonstrate that more complex models accounting for across-site compositional heterogeneity can ameliorate this bias. We present CAT-posterior mean site frequencies (PMSF), a pipeline for diagnosing and resolving phylogenetic bias resulting from inadequate modeling of across-site compositional heterogeneity based on the CAT model. CAT-PMSF is robust against long-branch attraction in all alignments we have examined. We suggest using CAT-PMSF when convergence of the CAT model cannot be assured. We find evidence that compositionally constrained sites are driving long-branch attraction in two metazoan datasets and recover evidence for Porifera as the sister group to all other animals. [Animal phylogeny; cross-site heterogeneity; long-branch attraction; phylogenomics.].

High speciation rate of niche specialists in hot springs

Ecological and evolutionary processes simultaneously regulate microbial diversity, but the evolutionary processes and their driving forces remain largely unexplored. Here we investigated the ecological and evolutionary characteristics of microbiota in hot springs spanning a broad temperature range (54.8-80 °C) by sequencing the 16S rRNA genes. Our results demonstrated that niche specialists and niche generalists are embedded in a complex interaction of ecological and evolutionary dynamics. On the thermal tolerance niche axis, thermal (T) sensitive (at a specific temperature) versus T-resistant (at least in five temperatures) species were characterized by different niche breadth, community abundance and dispersal potential, consequently differing in potential evolutionary trajectory. The niche-specialized T-sensitive species experienced strong temperature barriers, leading to completely species shift and high fitness but low abundant communities at each temperature ("home niche"), and such trade-offs thus reinforced peak performance, as evidenced by high speciation across temperatures and increasing diversification potential with temperature. In contrast, T-resistant species are advantageous of niche expansion but with poor local performance, as shown by wide niche breadth with high extinction, indicating these niche generalists are "jack-of-all-trades, master-of-none". Despite of such differences, the T-sensitive and T-resistant species are evolutionarily interacted. Specifically, the continuous transition from T-sensitive to T-resistant species insured the exclusion probability of T-resistant species at a relatively constant level across temperatures. The co-evolution and co-adaptation of T-sensitive and T-resistant species were in line with the red queen theory. Collectively, our findings demonstrate that high speciation of niche specialists could alleviate the environmental-filtering-induced negative effect on diversity.

Temperature-Promoted Giant Unilamellar Vesicle (GUV) Aggregation: A Way of Multicellular Formation

The evolution of unicellular to multicellular life is considered to be an important step in the origin of life, and it is crucial to study the influence of environmental factors on this process through cell models in the laboratory. In this paper, we used giant unilamellar vesicles (GUVs) as a cell model to investigate the relationship between environmental temperature changes and the evolution of unicellular to multicellular life. The zeta potential of GUVs and the conformation of the headgroup of phospholipid molecules at different temperatures were examined using phase analysis light scattering (PALS) and attenuated total reflection-Fourier transform infrared spectroscopy (ATR-FTIR), respectively. In addition, the effect of increasing temperature on the aggregation of GUVs was further investigated in ionic solutions, and the possible mechanisms involved were explored. The results showed that increasing temperature reduced the repulsive forces between cells models and promoted their aggregation. This study could effectively contribute to our understanding of the evolution of primitive unicellular to multicellular life.

Genomic Legacies of Ancient Adaptation Illuminate GC-Content Evolution in Bacteria

Bacterial evolution is characterized by strong purifying selection as well as rapid adaptive evolution in changing environments. In this context, the genomic GC content (genomic GC) varies greatly but presents some level of phylogenetic stability, making it challenging to explain based on current hypotheses. To illuminate the evolutionary mechanisms of the genomic GC, we analyzed the base composition and functional inventory of 11,083 representative genomes. A phylogenetically constrained bimodal distribution of the genomic GC, which mainly originated from parallel divergences in the early evolution, was demonstrated. Such variation of the genomic GC can be well explained by DNA replication and repair (DRR), in which multiple pathways correlate with the genomic GC. Furthermore, the biased conservation of various stress-related genes, especially the DRR-related ones, implies distinct adaptive processes in the ancestral lineages of high- or low-GC clades which are likely induced by major environmental changes. Our findings support that the mutational biases resulting from these legacies of ancient adaptation have changed the course of adaptive evolution and generated great variation in the genomic GC. This highlights the importance of indirect effects of natural selection, which indicates a new model for bacterial evolution. IMPORTANCE GC content has been shown to be an important factor in microbial ecology and evolution, and the genomic GC of bacteria can be characterized by great intergenomic heterogeneity, high intragenomic homogeneity, and strong phylogenetic inertia, as well as being associated with the environment. Current hypotheses concerning direct selection or mutational biases cannot well explain these features simultaneously. Our findings of the genomic GC showing that ancient adaptations have transformed the DRR system and that the resulting mutational biases further contributed to a bimodal distribution of it offer a more reasonable scenario for the mechanism. This would imply that, when thinking about the evolution of life, diverse processes of adaptation exist, and combined effects of natural selection should be considered.

Ancestral Sequence Reconstruction of the Ribosomal Protein uS8 and Reduction of Amino Acid Usage to a Smaller Alphabet

Understanding the origin and early evolution of proteins is important for unveiling how the RNA world developed into an RNA-protein world. Because the composition of organic molecules in the Earth's primitive environment was plausibly not as diverse as today, the number of different amino acids used in early protein synthesis is likely to be substantially less than the current 20 proteinogenic residues. In this study, we have explored the thermal stability and RNA binding of ancestral variants of the ribosomal protein uS8 constructed from a reduced-alphabet of amino acids. First, we built a phylogenetic tree based on the amino acid sequences of uS8 from multiple extant organisms and used the tree to infer two plausible amino acid sequences corresponding to the last bacterial common ancestor of uS8. Both ancestral proteins were thermally stable and bound to an RNA fragment. By eliminating individual amino acid letters and monitoring thermal stability and RNA binding in the resulting proteins, we reduced the size of the amino acid set constituting one of the ancestral proteins, eventually finding that convergent sequences consisting of 15- or 14-amino acid alphabets still folded into stable structures that bound to the RNA fragment. Furthermore, a simplified variant reconstructed from a 13-amino-acid alphabet retained affinity for the RNA fragment, although it lost conformational stability. Collectively, RNA-binding activity may be achieved with a subset of the current 20 amino acids, raising the possibility of a simpler composition of RNA-binding proteins in the earliest stage of protein evolution.

Evolution of CRISPR-associated endonucleases as inferred from resurrected proteins

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated Cas9 is an effector protein that targets invading DNA and plays a major role in the prokaryotic adaptive immune system. Although Streptococcus pyogenes CRISPR-Cas9 has been widely studied and repurposed for applications including genome editing, its origin and evolution are poorly understood. Here, we investigate the evolution of Cas9 from resurrected ancient nucleases (anCas) in extinct firmicutes species that last lived 2.6 billion years before the present. We demonstrate that these ancient forms were much more flexible in their guide RNA and protospacer-adjacent motif requirements compared with modern-day Cas9 enzymes. Furthermore, anCas portrays a gradual palaeoenzymatic adaptation from nickase to double-strand break activity, exhibits high levels of activity with both single-stranded DNA and single-stranded RNA targets and is capable of editing activity in human cells. Prediction and characterization of anCas with a resurrected protein approach uncovers an evolutionary trajectory leading to functionally flexible ancient enzymes.

Proteome-wide 3D structure prediction provides insights into the ancestral metabolism of ancient archaea and bacteria

Ancestral metabolism has remained controversial due to a lack of evidence beyond sequence-based reconstructions. Although prebiotic chemists have provided hints that metabolism might originate from non-enzymatic protometabolic pathways, gaps between ancestral reconstruction and prebiotic processes mean there is much that is still unknown. Here, we apply proteome-wide 3D structure predictions and comparisons to investigate ancestorial metabolism of ancient bacteria and archaea, to provide information beyond sequence as a bridge to the prebiotic processes. We compare representative bacterial and archaeal strains, which reveal surprisingly similar physiological and metabolic characteristics via microbiological and biophysical experiments. Pairwise comparison of protein structures identify the conserved metabolic modules in bacteria and archaea, despite interference from overly variable sequences. The conserved modules (for example, middle of glycolysis, partial TCA, proton/sulfur respiration, building block biosynthesis) constitute the basic functions that possibly existed in the archaeal-bacterial common ancestor, which are remarkably consistent with the experimentally confirmed protometabolic pathways. These structure-based findings provide a new perspective to reconstructing the ancestral metabolism and understanding its origin, which suggests high-throughput protein 3D structure prediction is a promising approach, deserving broader application in future ancestral exploration.

Engineering functional thermostable proteins using ancestral sequence reconstruction

Natural proteins are often only slightly more stable in the native state than the denatured state, and an increase in environmental temperature can easily shift the balance towards unfolding. Therefore, the engineering of proteins to improve protein stability is an area of intensive research. Thermostable proteins are required to withstand industrial process conditions, for increased shelf-life of protein therapeutics, for developing robust 'biobricks' for synthetic biology applications, and for research purposes (e.g. structure determination). In addition, thermostability buffers the often destabilizing effects of mutations introduced to improve other properties. Rational design approaches to engineering thermostability require structural information, but even with advanced computational methods, it is challenging to predict or parameterize all the relevant structural factors with sufficient precision to anticipate the results of a given mutation. Directed evolution is an alternative when structures are unavailable but requires extensive screening of mutant libraries. Recently however, bioinspired approaches based on phylogenetic analyses have shown great promise. Leveraging the rapid expansion in sequence data and bioinformatic tools, ancestral sequence reconstruction (ASR) can generate highly stable folds for novel applications in industrial chemistry, medicine, and synthetic biology. This review provides an overview of the factors important for successful inference of thermostable proteins by ASR and what it can reveal about the determinants of stability in proteins.

Resolving the History of Life on Earth by Seeking Life As We Know It on Mars

An origin of Earth life on Mars would resolve significant inconsistencies between the inferred history of life and Earth's geologic history. Life as we know it utilizes amino acids, nucleic acids, and lipids for the metabolic, informational, and compartment-forming subsystems of a cell. Such building blocks may have formed simultaneously from cyanosulfidic chemical precursors in a planetary surface scenario involving ultraviolet light, wet-dry cycling, and volcanism. On the inferred water world of early Earth, such an origin would have been limited to volcanic island hotspots. A cyanosulfidic origin of life could have taken place on Mars via photoredox chemistry, facilitated by orders-of-magnitude more sub-aerial crust than early Earth, and an earlier transition to oxidative conditions that could have been involved in final fixation of the genetic code. Meteoritic bombardment may have generated transient habitable environments and ejected and transferred life to Earth. Ongoing and future missions to Mars offer an unprecedented opportunity to confirm or refute evidence consistent with a cyanosulfidic origin of life on Mars, search for evidence of ancient life, and constrain the evolution of Mars' oxidation state over time. We should seek to prove or refute a martian origin for life on Earth alongside other possibilities.

The hyperthermophilic archaeon Thermococcus kodakarensis is resistant to pervasive negative supercoiling activity of DNA gyrase

In all cells, DNA topoisomerases dynamically regulate DNA supercoiling allowing essential DNA processes such as transcription and replication to occur. How this complex system emerged in the course of evolution is poorly understood. Intriguingly, a single horizontal gene transfer event led to the successful establishment of bacterial gyrase in Archaea, but its emergent function remains a mystery. To better understand the challenges associated with the establishment of pervasive negative supercoiling activity, we expressed the gyrase of the bacterium Thermotoga maritima in a naïve archaeon Thermococcus kodakarensis which naturally has positively supercoiled DNA. We found that the gyrase was catalytically active in T. kodakarensis leading to strong negative supercoiling of plasmid DNA which was stably maintained over at least eighty generations. An increased sensitivity of gyrase-expressing T. kodakarensis to ciprofloxacin suggested that gyrase also modulated chromosomal topology. Accordingly, global transcriptome analyses revealed large scale gene expression deregulation and identified a subset of genes responding to the negative supercoiling activity of gyrase. Surprisingly, the artificially introduced dominant negative supercoiling activity did not have a measurable effect on T. kodakarensis growth rate. Our data suggest that gyrase can become established in Thermococcales archaea without critically interfering with DNA transaction processes.

An Overview of 7α- and 7β-Hydroxysteroid Dehydrogenases: Structure, Specificity and Practical Application

7α-Hydroxysteroid dehydrogenase and 7β-hydroxysteroid dehydrogenase are key enzymes involved in bile acid metabolism. They catalyze the epimerization of a hydroxyl group through 7-keto bile acid intermediates. Basic research of the two enzymes has focused on exploring new enzymes and the structure-function relationship. The application research focused on the in vitro biosynthesis of bile acid drugs and the exploration and improvement of their catalytic ability based on molecular engineering. This article summarized the primary and advanced structural characteristics, specificities, biochemical properties, and applications of the two enzymes. The emphasis is also given to obtaining of novel 7α-hydroxysteroid dehydrogenase and 7β-hydroxysteroid dehydrogenase that are thermally stable and active in the presence of organic solvents, high substrate concentration, and extreme pH values. To achieve these goals, enzyme redesigning based on protein engineering and genomics may be the most useful approaches.

Introduction-First Billion Years: Habitability

The physical processes active during the first billion years (FBY) of Earth's history, such as accretion, differentiation, and impact cratering, provide constraints on the initial conditions that were conducive to the formation and establishment of life on Earth. This motivated the Lunar and Planetary Institute's FBY topical initiative, which was a four-part conference series intended to look at each of these physical processes to study the basic structure and composition of our Solar System that was set during the FBY. The FBY Habitability conference, held in September 2019, was the last in this series and was intended to synthesize the initiative; specifically, to further our understanding of the origins of life, planetary and environmental habitability, and the search for life beyond Earth. The conference included discussions of planetary habitability and the potential emergence of life on bodies within our Solar System, as well as extrasolar systems by applying our knowledge of the Solar System's FBY, and in particular Earth's early history. To introduce this Special Collection, which resulted from work discussed at the conference, we provide a review of the main themes and a synopsis of the FBY Habitability conference.

The Mutational Robustness of the Genetic Code and Codon Usage in Environmental Context: A Non-Extremophilic Preference?

The genetic code was evolved, to some extent, to minimize the effects of mutations. The effects of mutations depend on the amino acid repertoire, the structure of the genetic code and frequencies of amino acids in proteomes. The amino acid compositions of proteins and corresponding codon usages are still under selection, which allows us to ask what kind of environment the standard genetic code is adapted to. Using simple computational models and comprehensive datasets comprising genomic and environmental data from all three domains of Life, we estimate the expected severity of non-synonymous genomic mutations in proteins, measured by the change in amino acid physicochemical properties. We show that the fidelity in these physicochemical properties is expected to deteriorate with extremophilic codon usages, especially in thermophiles. These findings suggest that the genetic code performs better under non-extremophilic conditions, which not only explains the low substitution rates encountered in halophiles and thermophiles but the revealed relationship between the genetic code and habitat allows us to ponder on earlier phases in the history of Life.

Plausibility of Early Life in a Relatively Wide Temperature Range: Clues from Simulated Metabolic Network Expansion

The debate on the temperature of the environment where life originated is still inconclusive. Metabolic reactions constitute the basis of life, and may be a window to the world where early life was born. Temperature is an important parameter of reaction thermodynamics, which determines whether metabolic reactions can proceed. In this study, the scale of the prebiotic metabolic network at different temperatures was examined by a thermodynamically constrained network expansion simulation. It was found that temperature has limited influence on the scale of the simulated metabolic networks, implying that early life may have occurred in a relatively wide temperature range.

The Molecular Determinants of Thermoadaptation: Methanococcales as a Case Study

Previous reports have shown that environmental temperature impacts proteome evolution in Bacteria and Archaea. However, it is unknown whether thermoadaptation mainly occurs via the sequential accumulation of substitutions, massive horizontal gene transfers, or both. Measuring the real contribution of amino acid substitution to thermoadaptation is challenging, because of confounding environmental and genetic factors (e.g., pH, salinity, genomic G + C content) that also affect proteome evolution. Here, using Methanococcales, a major archaeal lineage, as a study model, we show that optimal growth temperature is the major factor affecting variations in amino acid frequencies of proteomes. By combining phylogenomic and ancestral sequence reconstruction approaches, we disclose a sequential substitutional scheme in which lysine plays a central role by fine tuning the pool of arginine, serine, threonine, glutamine, and asparagine, whose frequencies are strongly correlated with optimal growth temperature. Finally, we show that colonization to new thermal niches is not associated with high amounts of horizontal gene transfers. Altogether, although the acquisition of a few key proteins through horizontal gene transfer may have favored thermoadaptation in Methanococcales, our findings support sequential amino acid substitutions as the main factor driving thermoadaptation.

Phylogenetic analysis of mutational robustness based on codon usage supports that the standard genetic code does not prefer extreme environments

The mutational robustness of the genetic code is rarely discussed in the context of biological diversity, such as codon usage and related factors, often considered as independent of the actual organism's proteome. Here we put the living beings back to picture and use distortion as a metric of mutational robustness. Distortion estimates the expected severities of non-synonymous mutations measuring it by amino acid physicochemical properties and weighting for codon usage. Using the biological variance of codon frequencies, we interpret the mutational robustness of the standard genetic code with regards to their corresponding environments and genomic compositions (GC-content). Employing phylogenetic analyses, we show that coding fidelity in physicochemical properties can deteriorate with codon usages adapted to extreme environments and these putative effects are not the artefacts of phylogenetic bias. High temperature environments select for codon usages with decreased mutational robustness of hydrophobic, volumetric, and isoelectric properties. Selection at high saline concentrations also leads to reduced fidelity in polar and isoelectric patterns. These show that the genetic code performs best with mesophilic codon usages, strengthening the view that LUCA or its ancestors preferred lower temperature environments. Taxonomic implications, such as rooting the tree of life, are also discussed.

Resurrection of Ancestral Malate Dehydrogenases Reveals the Evolutionary History of Halobacterial Proteins : Deciphering Gene Trajectories and Changes in Biochemical Properties

Extreme halophilic Archaea thrive in high salt, where, through proteomic adaptation, they cope with the strong osmolarity and extreme ionic conditions of their environment. In spite of wide fundamental interest, however, studies providing insights into this adaptation are scarce, because of practical difficulties inherent to the purification and characterization of halophilic enzymes. In this work, we describe the evolutionary history of malate dehydrogenases (MalDH) within Halobacteria (a class of the Euryarchaeota phylum). We resurrected nine ancestors along the inferred halobacterial MalDH phylogeny, including the Last Common Ancestral MalDH of Halobacteria (LCAHa) and compared their biochemical properties with those of five modern halobacterial MalDHs. We monitored the stability of these various MalDHs, their oligomeric states and enzymatic properties, as a function of concentration for different salts in the solvent. We found that a variety of evolutionary processes such as amino acid replacement, gene duplication, loss of MalDH gene and replacement owing to horizontal transfer resulted in significant differences in solubility, stability and catalytic properties between these enzymes in the three Halobacteriales, Haloferacales and Natrialbales orders since the LCAHa MalDH.We also showed how a stability trade-off might favor the emergence of new properties during adaptation to diverse environmental conditions. Altogether, our results suggest a new view of halophilic protein adaptation in Archaea.

In search of the RNA world on Mars

Advances in origins of life research and prebiotic chemistry suggest that life as we know it may have emerged from an earlier RNA World. However, it has been difficult to reconcile the conditions used in laboratory experiments with real-world geochemical environments that may have existed on the early Earth and hosted the origin(s) of life. This challenge is due to geologic resurfacing and recycling that have erased the overwhelming majority of the Earth's prebiotic history. We therefore propose that Mars, a planet frozen in time, comprised of many surfaces that have remained relatively unchanged since their formation > 4 Gya, is the best alternative to search for environments consistent with geochemical requirements imposed by the RNA world. In this study, we synthesize in situ and orbital observations of Mars and modeling of its early atmosphere into solutions containing a range of pHs and concentrations of prebiotically relevant metals (Fe²⁺ , Mg²⁺ , and Mn²⁺ ) spanning various candidate aqueous environments. We then experimentally determine RNA degradation kinetics due to metal-catalyzed hydrolysis (cleavage) and evaluate whether early Mars could have been permissive toward the accumulation of long-lived RNA polymers. Our results indicate that a Mg²⁺ -rich basalt sourcing metals to a slightly acidic (pH 5.4) environment mediates the slowest rates of RNA cleavage, though geologic evidence and basalt weathering models suggest aquifers on Mars would be near neutral (pH ~ 7). Moreover, the early onset of oxidizing conditions on Mars has major consequences regarding the availability of oxygen-sensitive metals (i.e., Fe²⁺ and Mn²⁺ ) due to increased RNA degradation rates and precipitation. Overall, (a) low pH decreases RNA cleavage at high metal concentrations; (b) acidic to neutral pH environments with Fe²⁺ or Mn²⁺ cleave more RNA than Mg²⁺ ; and (c) alkaline environments with Mg²⁺ dramatically cleaves more RNA while precipitates were observed for Fe²⁺ and Mn²⁺ .

Exploring short k-mer profiles in cells and mobile elements from Archaea highlights the major influence of both the ecological niche and evolutionary history

BACKGROUND

K-mer-based methods have greatly advanced in recent years, largely driven by the realization of their biological significance and by the advent of next-generation sequencing. Their speed and their independence from the annotation process are major advantages. Their utility in the study of the mobilome has recently emerged and they seem a priori adapted to the patchy gene distribution and the lack of universal marker genes of viruses and plasmids. To provide a framework for the interpretation of results from k-mer based methods applied to archaea or their mobilome, we analyzed the 5-mer DNA profiles of close to 600 archaeal cells, viruses and plasmids. Archaea is one of the three domains of life. Archaea seem enriched in extremophiles and are associated with a high diversity of viral and plasmid families, many of which are specific to this domain. We explored the dataset structure by multivariate and statistical analyses, seeking to identify the underlying factors.

RESULTS

For cells, the 5-mer profiles were inconsistent with the phylogeny of archaea. At a finer taxonomic level, the influence of the taxonomy and the environmental constraints on 5-mer profiles was very strong. These two factors were interdependent to a significant extent, and the respective weights of their contributions varied according to the clade. A convergent adaptation was observed for the class Halobacteria, for which a strong 5-mer signature was identified. For mobile elements, coevolution with the host had a clear influence on their 5-mer profile. This enabled us to identify one previously known and one new case of recent host transfer based on the atypical composition of the mobile elements involved. Beyond the effect of coevolution, extrachromosomal elements strikingly retain the specific imprint of their own viral or plasmid taxonomic family in their 5-mer profile.

CONCLUSION

This specific imprint confirms that the evolution of extrachromosomal elements is driven by multiple parameters and is not restricted to host adaptation. In addition, we detected only recent host transfer events, suggesting the fast evolution of short k-mer profiles. This calls for caution when using k-mers for host prediction, metagenomic binning or phylogenetic reconstruction.

Evolution of frustrated and stabilising contacts in reconstructed ancient proteins

Energetic properties of a protein are a major determinant of its evolutionary fitness. Using a reconstruction algorithm, dating the reconstructed proteins and calculating the interaction network between their amino acids through a coevolutionary approach, we studied how the interactions that stabilise 890 proteins, belonging to five families, evolved for billions of years. In particular, we focused our attention on the network of most strongly attractive contacts and on that of poorly optimised, frustrated contacts. Our results support the idea that the cluster of most attractive interactions extends its size along evolutionary time, but from the data, we cannot conclude that protein stability or that the degree of frustration tends always to decrease.

A Place for Viruses on the Tree of Life

Viruses are ubiquitous. They infect almost every species and are probably the most abundant biological entities on the planet, yet they are excluded from the Tree of Life (ToL). However, there can be no doubt that viruses play a significant role in evolution, the force that facilitates all life on Earth. Conceptually, viruses are regarded by many as non-living entities that hijack living cells in order to propagate. A strict separation between living and non-living entities places viruses far from the ToL, but this may be theoretically unsound. Advances in sequencing technology and comparative genomics have expanded our understanding of the evolutionary relationships between viruses and cellular organisms. Genomic and metagenomic data have revealed that co-evolution between viral and cellular genomes involves frequent horizontal gene transfer and the occasional co-option of novel functions over evolutionary time. From the giant, ameba-infecting marine viruses to the tiny Porcine circovirus harboring only two genes, viruses and their cellular hosts are ecologically and evolutionarily intertwined. When deciding how, if, and where viruses should be placed on the ToL, we should remember that the Tree functions best as a model of biological evolution on Earth, and it is important that models themselves evolve with our increasing understanding of biological systems.

Bayesian genome scale modelling identifies thermal determinants of yeast metabolism

The molecular basis of how temperature affects cell metabolism has been a long-standing question in biology, where the main obstacles are the lack of high-quality data and methods to associate temperature effects on the function of individual proteins as well as to combine them at a systems level. Here we develop and apply a Bayesian modeling approach to resolve the temperature effects in genome scale metabolic models (GEM). The approach minimizes uncertainties in enzymatic thermal parameters and greatly improves the predictive strength of the GEMs. The resulting temperature constrained yeast GEM uncovers enzymes that limit growth at superoptimal temperatures, and squalene epoxidase (ERG1) is predicted to be the most rate limiting. By replacing this single key enzyme with an ortholog from a thermotolerant yeast strain, we obtain a thermotolerant strain that outgrows the wild type, demonstrating the critical role of sterol metabolism in yeast thermosensitivity. Therefore, apart from identifying thermal determinants of cell metabolism and enabling the design of thermotolerant strains, our Bayesian GEM approach facilitates modelling of complex biological systems in the absence of high-quality data and therefore shows promise for becoming a standard tool for genome scale modeling.

While temperature impacts the function of all cellular components, it’s hard to rule out how the temperature dependence of cell phenotypes emerged from the dependence of individual components. Here, the authors develop a Bayesian genome scale modelling approach to identify thermal determinants of yeast metabolism.

Bedrock radioactivity influences the rate and spectrum of mutation

All organisms on Earth are exposed to low doses of natural radioactivity but some habitats are more radioactive than others. Yet, documenting the influence of natural radioactivity on the evolution of biodiversity is challenging. Here, we addressed whether organisms living in naturally more radioactive habitats accumulate more mutations across generations using 14 species of waterlice living in subterranean habitats with contrasted levels of radioactivity. We found that the mitochondrial and nuclear mutation rates across a waterlouse species’ genome increased on average by 60% and 30%, respectively, when radioactivity increased by a factor of three. We also found a positive correlation between the level of radioactivity and the probability of G to T (and complementary C to A) mutations, a hallmark of oxidative stress. We conclude that even low doses of natural bedrock radioactivity influence the mutation rate possibly through the accumulation of oxidative damage, in particular in the mitochondrial genome.

Revisiting the Relationships Between Genomic G + C Content, RNA Secondary Structures, and Optimal Growth Temperature

Over twenty years ago Galtier and Lobry published a manuscript entitled "Relationships between Genomic G + C Content, RNA Secondary Structure, and Optimal Growth Temperature" in the Journal of Molecular Evolution that showcased the lack of a relationship between genomic G + C content and optimal growth temperature (OGT) in a set of about 200 prokaryotes. Galtier and Lobry also assessed the relationship between RNA secondary structures (rRNA stems, tRNAs) and OGT, and in this case a clear relationship emerged. Increasing structured RNA G + C content (particularly in regions that are double-stranded) correlates with increased OGT. Both of these fundamental relationships have withstood test of many additional sequences and spawned a variety of different applications that include prediction of OGT from rRNA sequence and computational ncRNA identification approaches. In this work, I present the motivation behind Galtier and Lobry's original paper and the larger questions addressed by the work, how these questions have evolved over the last two decades, and the impact of Galtier and Lobry's manuscript in fields beyond these questions.

The LUCA and its complex virome

The last universal cellular ancestor (LUCA) is the most recent population of organisms from which all cellular life on Earth descends. The reconstruction of the genome and phenotype of the LUCA is a major challenge in evolutionary biology. Given that all life forms are associated with viruses and/or other mobile genetic elements, there is no doubt that the LUCA was a host to viruses. Here, by projecting back in time using the extant distribution of viruses across the two primary domains of life, bacteria and archaea, and tracing the evolutionary histories of some key virus genes, we attempt a reconstruction of the LUCA virome. Even a conservative version of this reconstruction suggests a remarkably complex virome that already included the main groups of extant viruses of bacteria and archaea. We further present evidence of extensive virus evolution antedating the LUCA. The presence of a highly complex virome implies the substantial genomic and pan-genomic complexity of the LUCA itself.

Cofactors as Molecular Fossils To Trace the Origin and Evolution of Proteins

Due to their early origin and extreme conservation, cofactors are valuable molecular fossils for tracing the origin and evolution of proteins. First, as the order of protein folds binding with cofactors roughly coincides with protein-fold chronology, cofactors are considered to have facilitated the origin of primitive proteins by selecting them from pools of random amino acid sequences. Second, in the subsequent evolution of proteins, cofactors still played an important role. More interestingly, as metallic cofactors evolved with geochemical variations, some geochemical events left imprints in the chronology of protein architecture; this provides further evidence supporting the coevolution of biochemistry and geochemistry. In this paper, we attempt to review the molecular fossils used in tracing the origin and evolution of proteins, with a special focus on cofactors.

Structural studies of geranylgeranylglyceryl phosphate synthase, a prenyltransferase found in thermophilic Euryarchaeota

Archaea are uniquely adapted to thrive in harsh environments, and one of these adaptations involves the archaeal membrane lipids, which are characterized by their isoprenoid alkyl chains connected via ether linkages to glycerol 1-phosphate. The membrane lipids of the thermophilic and acidophilic euryarchaeota Thermoplasma volcanium are exclusively glycerol dibiphytanyl glycerol tetraethers. The first committed step in the biosynthetic pathway of these archaeal lipids is the formation of the ether linkage between glycerol 1-phosphate and geranylgeranyl diphosphate, and is catalyzed by the enzyme geranylgeranylglyceryl phosphate synthase (GGGPS). The 1.72 Å resolution crystal structure of GGGPS from T. volcanium (TvGGGPS) in complex with glycerol and sulfate is reported here. The crystal structure reveals TvGGGPS to be a dimer, which is consistent with the absence of the aromatic anchor residue in helix α5a that is required for hexamerization in other GGGPS homologs; the hexameric quaternary structure in GGGPS is thought to provide thermostability. A phylogenetic analysis of the Euryarchaeota and a parallel ancestral state reconstruction investigated the relationship between optimal growth temperature and the ancestral sequences. The presence of an aromatic anchor residue is not explained by temperature as an ecological parameter. An examination of the active site of the TvGGGPS dimer revealed that it may be able to accommodate longer isoprenoid substrates, supporting an alternative pathway of isoprenoid membrane-lipid synthesis.

Multidomain ribosomal protein trees and the planctobacterial origin of neomura (eukaryotes, archaebacteria)

Palaeontologically, eubacteria are > 3× older than neomura (eukaryotes, archaebacteria). Cell biology contrasts ancestral eubacterial murein peptidoglycan walls and derived neomuran N-linked glycoprotein coats/walls. Misinterpreting long stems connecting clade neomura to eubacteria on ribosomal sequence trees (plus misinterpreted protein paralogue trees) obscured this historical pattern. Universal multiprotein ribosomal protein (RP) trees, more accurate than rRNA trees, are taxonomically undersampled. To reduce contradictions with genically richer eukaryote trees and improve eubacterial phylogeny, we constructed site-heterogeneous and maximum-likelihood universal three-domain, two-domain, and single-domain trees for 143 eukaryotes (branching now congruent with 187-protein trees), 60 archaebacteria, and 151 taxonomically representative eubacteria, using 51 and 26 RPs. Site-heterogeneous trees greatly improve eubacterial phylogeny and higher classification, e.g. showing gracilicute monophyly, that many 'rDNA-phyla' belong in Proteobacteria, and reveal robust new phyla Synthermota and Aquithermota. Monoderm Posibacteria and Mollicutes (two separate wall losses) are both polyphyletic: multiple outer membrane losses in Endobacteria occurred separately from Actinobacteria; neither phylum is related to Chloroflexi, the most divergent prokaryotes, which originated photosynthesis (new model proposed). RP trees support an eozoan root for eukaryotes and are consistent with archaebacteria being their sisters and rooted between Filarchaeota (=Proteoarchaeota, including 'Asgardia') and Euryarchaeota sensu-lato (including ultrasimplified 'DPANN' whose long branches often distort trees). Two-domain trees group eukaryotes within Planctobacteria, and archaebacteria with Planctobacteria/Sphingobacteria. Integrated molecular/palaeontological evidence favours negibacterial ancestors for neomura and all life. Unique presence of key pre-neomuran characters favours Planctobacteria only as ancestral to neomura, which apparently arose by coevolutionary repercussions (explained here in detail, including RP replacement) of simultaneous outer membrane and murein loss. Planctobacterial C-1 methanotrophic enzymes are likely ancestral to archaebacterial methanogenesis and β-propeller-α-solenoid proteins to eukaryotic vesicle coats, nuclear-pore-complexes, and intraciliary transport. Planctobacterial chaperone-independent 4/5-protofilament microtubules and MamK actin-ancestors prepared for eukaryote intracellular motility, mitosis, cytokinesis, and phagocytosis. We refute numerous wrong ideas about the universal tree.

Multidomain ribosomal protein trees and the planctobacterial origin of neomura (eukaryotes, archaebacteria)

Palaeontologically, eubacteria are > 3× older than neomura (eukaryotes, archaebacteria). Cell biology contrasts ancestral eubacterial murein peptidoglycan walls and derived neomuran N-linked glycoprotein coats/walls. Misinterpreting long stems connecting clade neomura to eubacteria on ribosomal sequence trees (plus misinterpreted protein paralogue trees) obscured this historical pattern. Universal multiprotein ribosomal protein (RP) trees, more accurate than rRNA trees, are taxonomically undersampled. To reduce contradictions with genically richer eukaryote trees and improve eubacterial phylogeny, we constructed site-heterogeneous and maximum-likelihood universal three-domain, two-domain, and single-domain trees for 143 eukaryotes (branching now congruent with 187-protein trees), 60 archaebacteria, and 151 taxonomically representative eubacteria, using 51 and 26 RPs. Site-heterogeneous trees greatly improve eubacterial phylogeny and higher classification, e.g. showing gracilicute monophyly, that many 'rDNA-phyla' belong in Proteobacteria, and reveal robust new phyla Synthermota and Aquithermota. Monoderm Posibacteria and Mollicutes (two separate wall losses) are both polyphyletic: multiple outer membrane losses in Endobacteria occurred separately from Actinobacteria; neither phylum is related to Chloroflexi, the most divergent prokaryotes, which originated photosynthesis (new model proposed). RP trees support an eozoan root for eukaryotes and are consistent with archaebacteria being their sisters and rooted between Filarchaeota (=Proteoarchaeota, including 'Asgardia') and Euryarchaeota sensu-lato (including ultrasimplified 'DPANN' whose long branches often distort trees). Two-domain trees group eukaryotes within Planctobacteria, and archaebacteria with Planctobacteria/Sphingobacteria. Integrated molecular/palaeontological evidence favours negibacterial ancestors for neomura and all life. Unique presence of key pre-neomuran characters favours Planctobacteria only as ancestral to neomura, which apparently arose by coevolutionary repercussions (explained here in detail, including RP replacement) of simultaneous outer membrane and murein loss. Planctobacterial C-1 methanotrophic enzymes are likely ancestral to archaebacterial methanogenesis and β-propeller-α-solenoid proteins to eukaryotic vesicle coats, nuclear-pore-complexes, and intraciliary transport. Planctobacterial chaperone-independent 4/5-protofilament microtubules and MamK actin-ancestors prepared for eukaryote intracellular motility, mitosis, cytokinesis, and phagocytosis. We refute numerous wrong ideas about the universal tree.

Phylogenetic Analyses of Sites in Different Protein Structural Environments Result in Distinct Placements of the Metazoan Root

Phylogenomics, the use of large datasets to examine phylogeny, has revolutionized the study of evolutionary relationships. However, genome-scale data have not been able to resolve all relationships in the tree of life; this could reflect, at least in part, the poor-fit of the models used to analyze heterogeneous datasets. Some of the heterogeneity may reflect the different patterns of selection on proteins based on their structures. To test that hypothesis, we developed a pipeline to divide phylogenomic protein datasets into subsets based on secondary structure and relative solvent accessibility. We then tested whether amino acids in different structural environments had distinct signals for the topology of the deepest branches in the metazoan tree. We focused on a dataset that appeared to have a mixture of signals and we found that the most striking difference in phylogenetic signal reflected relative solvent accessibility. Analyses of exposed sites (residues located on the surface of proteins) yielded a tree that placed ctenophores sister to all other animals whereas sites buried inside proteins yielded a tree with a sponge+ctenophore clade. These differences in phylogenetic signal were not ameliorated when we conducted analyses using a set of maximum-likelihood profile mixture models. These models are very similar to the Bayesian CAT model, which has been used in many analyses of deep metazoan phylogeny. In contrast, analyses conducted after recoding amino acids to limit the impact of deviations from compositional stationarity increased the congruence in the estimates of phylogeny for exposed and buried sites; after recoding amino acid trees estimated using the exposed and buried site both supported placement of ctenophores sister to all other animals. Although the central conclusion of our analyses is that sites in different structural environments yield distinct trees when analyzed using models of protein evolution, our amino acid recoding analyses also have implications for metazoan evolution. Specifically, our results add to the evidence that ctenophores are the sister group of all other animals and they further suggest that the placozoa+cnidaria clade found in some other studies deserves more attention. Taken as a whole, these results provide striking evidence that it is necessary to achieve a better understanding of the constraints due to protein structure to improve phylogenetic estimation.

The Evolution of Reverse Gyrase Suggests a Nonhyperthermophilic Last Universal Common Ancestor

Reverse gyrase (RG) is the only protein found ubiquitously in hyperthermophilic organisms, but absent from mesophiles. As such, its simple presence or absence allows us to deduce information about the optimal growth temperature of long-extinct organisms, even as far as the last universal common ancestor of extant life (LUCA). The growth environment and gene content of the LUCA has long been a source of debate in which RG often features. In an attempt to settle this debate, we carried out an exhaustive search for RG proteins, generating the largest RG data set to date. Comprising 376 sequences, our data set allows for phylogenetic reconstructions of RG with unprecedented size and detail. These RG phylogenies are strikingly different from those of universal proteins inferred to be present in the LUCA, even when using the same set of species. Unlike such proteins, RG does not form monophyletic archaeal and bacterial clades, suggesting RG emergence after the formation of these domains, and/or significant horizontal gene transfer. Additionally, the branch lengths separating archaeal and bacterial groups are very short, inconsistent with the tempo of evolution from the time of the LUCA. Despite this, phylogenies limited to archaeal RG resolve most archaeal phyla, suggesting predominantly vertical evolution since the time of the last archaeal ancestor. In contrast, bacterial RG indicates emergence after the last bacterial ancestor followed by significant horizontal transfer. Taken together, these results suggest a nonhyperthermophilic LUCA and bacterial ancestor, with hyperthermophily emerging early in the evolution of the archaeal and bacterial domains.

The Archean atmosphere

The atmosphere of the Archean eon-one-third of Earth's history-is important for understanding the evolution of our planet and Earth-like exoplanets. New geological proxies combined with models constrain atmospheric composition. They imply surface O2 levels <10-6 times present, N2 levels that were similar to today or possibly a few times lower, and CO2 and CH4 levels ranging ~10 to 2500 and 102 to 104 times modern amounts, respectively. The greenhouse gas concentrations were sufficient to offset a fainter Sun. Climate moderation by the carbon cycle suggests average surface temperatures between 0° and 40°C, consistent with occasional glaciations. Isotopic mass fractionation of atmospheric xenon through the Archean until atmospheric oxygenation is best explained by drag of xenon ions by hydrogen escaping rapidly into space. These data imply that substantial loss of hydrogen oxidized the Earth. Despite these advances, detailed understanding of the coevolving solid Earth, biosphere, and atmosphere remains elusive, however.

Chemistry of Abiotic Nucleotide Synthesis

The chemistry of abiotic nucleotide synthesis of RNA and DNA in the context of their prebiotic origins on early earth is a continuing challenge. How did (or how can) the nucleotides form and assemble from the small molecule inventories and under conditions that prevailed on early earth 3.5-4 billion years ago? This review provides a background and up-to-date progress that will allow the reader to judge where the field stands currently and what remains to be achieved. We start with a brief primer on the biological synthesis of nucleotides, followed by an extensive focus on the prebiotic formation of the components of nucleotides-either via the synthesis of ribose and the canonical nucleobases and then joining them together or by building both the conjoined sugar and nucleobase, part-by-part-toward the ultimate goal of forming RNA and DNA by polymerization. The review will emphasize that there are-and will continue to be-many more questions than answers from the synthetic, mechanistic, and analytical perspectives. We wrap up the review with a cautionary note in this context about coming to conclusions as to whether the problem of chemistry of prebiotic nucleotide synthesis has been solved.

Viral Emerging Pathogen Evolution

The emergence of emerging and reemerging infectious diseases throughout history has been the product of the interaction of infectious agents, immunity, selective pressure, and environmental factors. The origin of the viral ancestors remains controversial, the debate remains whether the viruses existed before their host cells, evolved as molecular structures that adapted to the intracellular development or of the cells that infected or housed them as hosts. The nature of the viruses could propose a hypothetical evolutionary model based on the way they parasitize cells and the way they fuse with their host cells. So far there are no “universal” viral genes that have not allowed the hypotheses of the “common viral ancestor” to be conclusive. Insects have a special role in understanding viral evolution and how they became pathogens, suggesting a coevolutionary scenario of colonization that would initially allow the survival of viruses and later an ecological niche required for the cycle of virus replication. The coevolution of viruses and defense systems of hosts represents a fundamental role in the evolution of both viruses and host cells, a kind of arms race between viruses and hosts, and the deescalation of that war could explain the coexistence of viruses and hosts. Viruses are playing an important role in the understanding of cell evolution and have earned the right to be included in the tree of life. The aim of this chapter is to review the viral evolution and the insects as host.

Insights into the ecological roles and evolution of methyl-coenzyme M reductase-containing hot spring Archaea

Several recent studies have shown the presence of genes for the key enzyme associated with archaeal methane/alkane metabolism, methyl-coenzyme M reductase (Mcr), in metagenome-assembled genomes (MAGs) divergent to existing archaeal lineages. Here, we study the mcr-containing archaeal MAGs from several hot springs, which reveal further expansion in the diversity of archaeal organisms performing methane/alkane metabolism. Significantly, an MAG basal to organisms from the phylum Thaumarchaeota that contains mcr genes, but not those for ammonia oxidation or aerobic metabolism, is identified. Together, our phylogenetic analyses and ancestral state reconstructions suggest a mostly vertical evolution of mcrABG genes among methanogens and methanotrophs, along with frequent horizontal gene transfer of mcr genes between alkanotrophs. Analysis of all mcr-containing archaeal MAGs/genomes suggests a hydrothermal origin for these microorganisms based on optimal growth temperature predictions. These results also suggest methane/alkane oxidation or methanogenesis at high temperature likely existed in a common archaeal ancestor.

Methane metabolism by some lineages of Archaea contributes to the cycling of carbon on Earth. Here, the authors show high diversity of methyl-coenzyme M reductase (Mcr), a key enzyme associated with archaeal methane/alkane metabolism, in hot spring Archaea, and investigate their ecological roles and evolution.

Evolutionary Modes in Protein Observable Space: The Case of Thioredoxins

In this article, we investigated the structural and dynamical evolutionary behaviour of a set of ten thioredoxin proteins as formed by three extant forms and seven resurrected ones in laboratory. Starting from the crystallographic structures, we performed all-atom molecular dynamics simulations and compare the trajectories in terms of structural and dynamical properties. Interestingly, the structural properties related to the protein density (i.e. the number of residues divided by the excluded molecular volume) well describe the protein evolutionary behaviour. Our results also suggest that the changes in sequence as occurred during the evolution have affected the protein essential motions, allowing us to discriminate between ancient and extant proteins in terms of their dynamical behaviour. Such results are yet more evident when the bacterial, archaeal and eukaryotic thioredoxins are separately analysed.

A Shift in Aggregation Avoidance Strategy Marks a Long-Term Direction to Protein Evolution

To detect a direction to evolution, without the pitfalls of reconstructing ancestral states, we need to compare "more evolved" to "less evolved" entities. But because all extant species have the same common ancestor, none are chronologically more evolved than any other. However, different gene families were born at different times, allowing us to compare young protein-coding genes to those that are older and hence have been evolving for longer. To be retained during evolution, a protein must not only have a function, but must also avoid toxic dysfunction such as protein aggregation. There is conflict between the two requirements: hydrophobic amino acids form the cores of protein folds, but also promote aggregation. Young genes avoid strongly hydrophobic amino acids, which is presumably the simplest solution to the aggregation problem. Here we show that young genes' few hydrophobic residues are clustered near one another along the primary sequence, presumably to assist folding. The higher aggregation risk created by the higher hydrophobicity of older genes is counteracted by more subtle effects in the ordering of the amino acids, including a reduction in the clustering of hydrophobic residues until they eventually become more interspersed than if distributed randomly. This interspersion has previously been reported to be a general property of proteins, but here we find that it is restricted to old genes. Quantitatively, the index of dispersion delineates a gradual trend, i.e., a decrease in the clustering of hydrophobic amino acids over billions of years.

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.

Enhancement and Inhibitory Activities of Minerals for Alanine Oligopeptide Elongation Under Hydrothermal Conditions

In a previous study, we have showed that the elongation of an alanine oligopeptide [L-alanyl-L-alanyl-L-alanyl-L-alanine ((Ala)₄)] to higher oligopeptides is enhanced by calcite and dolomite at 275°C, using a mineral-mediated hydrothermal flow reactor system. However, a problem during the use of hydrothermal flow reactor system was that some of the minerals, such as clay, could not be tested due to their clogging in the reactor. In this article, we attempted to analyze the scope of enhancement for the formation of L-alanyl-L-alanyl-L-alanyl-L-alanyl-L-alanine ((Ala)₅) and higher oligopeptides with different minerals including clay minerals for the elongation of alanine oligopeptide at 175°C. First, carbonate minerals and some clay minerals showed an enhancement of the formation of (Ala)₅ from (Ala)₄. On the contrary, volcanic products showed strong inhibitory activities. According to the pH dependence on the (Ala)₄ elongations, we confirmed that most enhancement and inhibitory activities are due to the pH influence on the elongation of (Ala)₄. However, the enhancement of montmorillonite (Tsukinuno), sphalerite, apatite, tourmaline, calcite (Nitto Funka), and the inhibitory activities by volcanic ash (Shinmoedake), volcanic ash (Sakurajima), dickite, and pyrophillite are not simply due to the pH change in the presence of these minerals. The difference found between the previous and present studies suggests that the interaction kinetics of the aqueous phase with the mineral phase is also an important factor for the elongation of (Ala)₄. These data imply that the environments with pH near neutral to weak alkaline and with minerals might have been useful for the accumulation of oligopeptides in hydrothermal conditions.

A series of new E. coli-Thermococcus shuttle vectors compatible with previously existing vectors

Hyperthermophilic microorganisms are an important asset in the toolkits of biotechnologists, biochemists and evolutionary biologists. The anaerobic archaeon, Thermococcus kodakarensis, has become one of the most useful hyperthermophilic model species, not least due to its natural competence and genetic tractability. Despite this, the range of genetic tools available for T. kodakarensis remains limited. Using sequencing and phylogenetic analyses, we determined that the rolling-circle replication origin of the cryptic mini-plasmid pTP2 from T. prieurii is suitable for plasmid replication in T. kodakarensis. Based on this replication origin, we present a novel series of replicative E. coli–T. kodakarensis shuttle vectors. These shuttle vectors have been constructed with three different selectable markers, allowing selection in a range of T. kodakarensis backgrounds. Moreover, these pTP2-derived plasmids are compatible with the single-existing E. coli–T. kodakarensis shuttle vector, pLC70. We show that both pTP2-derived and pLC70-derived plasmids replicate faithfully while cohabitating in T. kodakarensis cells. These plasmids open the door for new areas of research in plasmid segregation, DNA replication and gene expression.

Metabolic shift at the class level sheds light on adaptation of methanogens to oxidative environments

Methanogens have long been considered strictly anaerobic and oxygen-sensitive microorganisms, but their ability to survive oxygen stress has also been documented. Indeed, methanogens have been found in oxidative environments, and antioxidant genes have been detected in their genomes. How methanogens adapt to oxidative environments, however, remain poorly understood. Here, we systematically predicted and annotated antioxidant features from representative genomes across six well-established methanogen orders. Based on functional gene content involved in production of reactive oxygen species, Hierarchical Clustering analyses grouped methanogens into two distinct clusters, corresponding to the Class I and II methanogens, respectively. Comparative genomics suggested a systematic shift in metabolisms across the two classes, resulting in an enrichment of antioxidant features in the Class II. Moreover, meta-analysis of 16 S rRNA gene sequences obtained from EnvDB indicated that members of Class II were more frequently recovered from microaerophilic and even oxic environments than the Class I members. Phylogenomic analysis suggested that the Class I and II methanogens might have evolved before and around the Great Oxygenation Event, respectively. The enrichment of antioxidant features in the Class II methanogens may have played a key role in the adaption of this group to oxidative environments today and historically.

Proteome Evolution of Deep-Sea Hydrothermal Vent Alvinellid Polychaetes Supports the Ancestry of Thermophily and Subsequent Adaptation to Cold in Some Lineages

Temperature, perhaps more than any other environmental factor, is likely to influence the evolution of all organisms. It is also a very interesting factor to understand how genomes are shaped by selection over evolutionary timescales, as it potentially affects the whole genome. Among thermophilic prokaryotes, temperature affects both codon usage and protein composition to increase the stability of the transcriptional/translational machinery, and the resulting proteins need to be functional at high temperatures. Among eukaryotes less is known about genome evolution, and the tube-dwelling worms of the family Alvinellidae represent an excellent opportunity to test hypotheses about the emergence of thermophily in ectothermic metazoans. The Alvinellidae are a group of worms that experience varying thermal regimes, presumably having evolved into these niches over evolutionary times. Here we analyzed 423 putative orthologous loci derived from 6 alvinellid species including the thermophilic Alvinella pompejana and Paralvinella sulfincola. This comparative approach allowed us to assess amino acid composition, codon usage, divergence, direction of residue changes and the strength of selection along the alvinellid phylogeny, and to design a new eukaryotic thermophilic criterion based on significant differences in the residue composition of proteins. Contrary to expectations, the alvinellid ancestor of all present-day species seems to have been thermophilic, a trait subsequently maintained by purifying selection in lineages that still inhabit higher temperature environments. In contrast, lineages currently living in colder habitats likely evolved under selective relaxation, with some degree of positive selection for low-temperature adaptation at the protein level.