The deep archaeal roots of eukaryotes

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.

Endosymbioses Have Shaped the Evolution of Biological Diversity and Complexity Time and Time Again

Life on Earth comprises prokaryotes and a broad assemblage of endosymbioses. The pages of Molecular Biology and Evolution and Genome Biology and Evolution have provided an essential window into how these endosymbiotic interactions have evolved and shaped biological diversity. Here, we provide a current perspective on this knowledge by drawing on decades of revelatory research published in Molecular Biology and Evolution and Genome Biology and Evolution, and insights from the field at large. The accumulated work illustrates how endosymbioses provide hosts with novel phenotypes that allow them to transition between adaptive landscapes to access environmental resources. Such endosymbiotic relationships have shaped and reshaped life on Earth. The early serial establishment of mitochondria and chloroplasts through endosymbioses permitted massive upscaling of cellular energetics, multicellularity, and terrestrial planetary greening. These endosymbioses are also the foundation upon which all later ones are built, including everything from land-plant endosymbioses with fungi and bacteria to nutritional endosymbioses found in invertebrate animals. Common evolutionary mechanisms have shaped this broad range of interactions. Endosymbionts generally experience adaptive and stochastic genome streamlining, the extent of which depends on several key factors (e.g. mode of transmission). Hosts, in contrast, adapt complex mechanisms of resource exchange, cellular integration and regulation, and genetic support mechanisms to prop up degraded symbionts. However, there are significant differences between endosymbiotic interactions not only in how partners have evolved with each other but also in the scope of their influence on biological diversity. These differences are important considerations for predicting how endosymbioses will persist and adapt to a changing planet.

An excavate root for the eukaryote tree of life

Much of the higher-order phylogeny of eukaryotes is well resolved, but the root remains elusive. We assembled a dataset of 183 eukaryotic proteins of archaeal ancestry to test this root. The resulting phylogeny identifies four lineages of eukaryotes currently classified as "Excavata" branching separately at the base of the tree. Thus, Parabasalia appear as the first major branch of eukaryotes followed sequentially by Fornicata, Preaxostyla, and Discoba. All four excavate branch points receive full statistical support from analyses with commonly used evolutionary models, a protein structure partition model that we introduce here, and various controls for deep phylogeny artifacts. The absence of aerobic mitochondria in Parabasalia, Fornicata, and Preaxostyla suggests that modern eukaryotes arose under anoxic conditions, probably much earlier than expected, and without the benefit of mitochondrial respiration.

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.

Diverse viruses of marine archaea discovered using metagenomics

During the past decade, metagenomics became a method of choice for the discovery of novel viruses. However, host assignment for uncultured viruses remains challenging, especially for archaeal viruses, which are grossly undersampled compared to viruses of bacteria and eukaryotes. Here, we assessed the utility of CRISPR spacer targeting, tRNA gene matching and homology searches for viral signature proteins, such as major capsid proteins, for the assignment of archaeal hosts and validated these approaches on metaviromes from Yangshan Harbor (YSH). We report 35 new genomes of viruses which could be confidently assigned to hosts representing diverse lineages of marine archaea. We show that the archaeal YSH virome is highly diverse, with some viruses enriching the previously described virus groups, such as magroviruses of Marine Group II Archaea (Poseidoniales), and others representing novel groups of marine archaeal viruses. Metagenomic recruitment of Tara Oceans datasets on the YSH viral genomes demonstrated the presence of YSH Poseidoniales and Nitrososphaeria viruses in the global oceans, but also revealed the endemic YSH-specific viral lineages. Furthermore, our results highlight the relationship between the soil and marine thaumarchaeal viruses. We propose three new families within the class Caudoviricetes for the classification of the five complete viral genomes predicted to replicate in marine Poseidoniales and Nitrososphaeria, two ecologically important and widespread archaeal groups. This study illustrates the utility of viral metagenomics in exploring the archaeal virome and provides new insights into the diversity, distribution and evolution of marine archaeal viruses.

Cargo Genes of Tn7-Like Transposons Comprise an Enormous Diversity of Defense Systems, Mobile Genetic Elements, and Antibiotic Resistance Genes

Transposition is a major mechanism of horizontal gene mobility in prokaryotes. However, exploration of the genes mobilized by transposons (cargo) is hampered by the difficulty in delineating integrated transposons from their surrounding genetic context. Here, we present a computational approach that allowed us to identify the boundaries of 6,549 Tn7-like transposons. We found that 96% of these transposons carry at least one cargo gene. Delineation of distinct communities in a gene-sharing network demonstrates how transposons function as a conduit of genes between phylogenetically distant hosts. Comparative analysis of the cargo genes reveals significant enrichment of mobile genetic elements (MGEs) nested within Tn7-like transposons, such as insertion sequences and toxin-antitoxin modules, and of genes involved in recombination, anti-MGE defense, and antibiotic resistance. More unexpectedly, cargo also includes genes encoding central carbon metabolism enzymes. Twenty-two Tn7-like transposons carry both an anti-MGE defense system and antibiotic resistance genes, illustrating how bacteria can overcome these combined pressures upon acquisition of a single transposon. This work substantially expands the distribution of Tn7-like transposons, defines their evolutionary relationships, and provides a large-scale functional classification of prokaryotic genes mobilized by transposition.

Ancient Gene Capture and Recent Gene Loss Shape the Evolution of Orthopoxvirus-Host Interaction Genes

The survival of viruses depends on their ability to resist host defenses and, of all animal virus families, the poxviruses have the most antidefense genes. Orthopoxviruses (ORPV), a genus within the subfamily Chordopoxvirinae, infect diverse mammals and include one of the most devastating human pathogens, the now eradicated smallpox virus. ORPV encode ∼200 genes, of which roughly half are directly involved in virus genome replication and expression as well as virion morphogenesis. The remaining ∼100 "accessory" genes are responsible for virus-host interactions, particularly counter-defense of innate immunity. Complete sequences are currently available for several hundred ORPV genomes isolated from a variety of mammalian hosts, providing a rich resource for comparative genomics and reconstruction of ORPV evolution. To identify the provenance and evolutionary trends of the ORPV accessory genes, we constructed clusters including the orthologs of these genes from all chordopoxviruses. Most of the accessory genes were captured in three major waves early in chordopoxvirus evolution, prior to the divergence of ORPV and the sister genus Centapoxvirus from their common ancestor. The capture of these genes from the host was followed by extensive gene duplication, yielding several paralogous gene families. In addition, nine genes were gained during the evolution of ORPV themselves. In contrast, nearly every accessory gene was lost, some on multiple, independent occasions in numerous lineages of ORPV, so that no ORPV retains them all. A variety of functional interactions could be inferred from examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently. IMPORTANCE Orthopoxviruses (ORPV) include smallpox (variola) virus, one of the most devastating human pathogens, and vaccinia virus, comprising the vaccine used for smallpox eradication. Among roughly 200 ORPV genes, about half are essential for genome replication and expression as well as virion morphogenesis, whereas the remaining half consists of accessory genes counteracting the host immune response. We reannotated the accessory genes of ORPV, predicting the functions of uncharacterized genes, and reconstructed the history of their gain and loss during the evolution of ORPV. Most of the accessory genes were acquired in three major waves antedating the origin of ORPV from chordopoxviruses. The evolution of ORPV themselves was dominated by gene loss, with numerous genes lost at the base of each major group of ORPV. Examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently during ORPV evolution allows prediction of different types of functional interactions.

Lauric Acid Is a Potent Biological Control Agent That Damages the Cell Membrane of Phytophthora sojae

Sustainable management of plant pathogens is becoming more challenging, and novel solutions are needed. Plant biologically active secondary metabolites are important sources of novel crop protection chemistry. Effective individual compounds of these natural products have the potential to be successful new agrochemicals. In this study, we identified lauric acid (LA) from soybean defense leaf volatiles. LA inhibited the growth of Phytophthora sojae, the causal agent of soybean root rot. It influenced mycelial development, sporangium formation, and zoospore generation and germination by damaging the P. sojae cell membrane. Additionally, we showed that LA and several of its derivatives, such as glycerol monolaurate (GML), had similar biological activities. Both LA and GML were safe to soybean plants when used at less than 0.3 g a.i./plant and could promote soybean growth, implying their potential as eco-friendly biological control agents.

Expanded diversity of Asgard archaea and their relationships with eukaryotes

Asgard is a recently discovered superphylum of archaea that appears to include the closest archaeal relatives of eukaryotes1-5. Debate continues as to whether the archaeal ancestor of eukaryotes belongs within the Asgard superphylum or whether this ancestor is a sister group to all other archaea (that is, a two-domain versus a three-domain tree of life)6-8. Here we present a comparative analysis of 162 complete or nearly complete genomes of Asgard archaea, including 75 metagenome-assembled genomes that-to our knowledge-have not previously been reported. Our results substantially expand the phylogenetic diversity of Asgard and lead us to propose six additional phyla that include a deep branch that we have provisionally named Wukongarchaeota. Our phylogenomic analysis does not resolve unequivocally the evolutionary relationship between eukaryotes and Asgard archaea, but instead-depending on the choice of species and conserved genes used to build the phylogeny-supports either the origin of eukaryotes from within Asgard (as a sister group to the expanded Heimdallarchaeota-Wukongarchaeota branch) or a deeper branch for the eukaryote ancestor within archaea. Our comprehensive protein domain analysis using the 162 Asgard genomes results in a major expansion of the set of eukaryotic signature proteins. The Asgard eukaryotic signature proteins show variable phyletic distributions and domain architectures, which is suggestive of dynamic evolution through horizontal gene transfer, gene loss, gene duplication and domain shuffling. The phylogenomics of the Asgard archaea points to the accumulation of the components of the mobile archaeal 'eukaryome' in the archaeal ancestor of eukaryotes (within or outside Asgard) through extensive horizontal gene transfer.

Thousands of previously unknown phages discovered in whole-community human gut metagenomes

BACKGROUND

Double-stranded DNA bacteriophages (dsDNA phages) play pivotal roles in structuring human gut microbiomes; yet, the gut virome is far from being fully characterized, and additional groups of phages, including highly abundant ones, continue to be discovered by metagenome mining. A multilevel framework for taxonomic classification of viruses was recently adopted, facilitating the classification of phages into evolutionary informative taxonomic units based on hallmark genes. Together with advanced approaches for sequence assembly and powerful methods of sequence analysis, this revised framework offers the opportunity to discover and classify unknown phage taxa in the human gut.

RESULTS

A search of human gut metagenomes for circular contigs encoding phage hallmark genes resulted in the identification of 3738 apparently complete phage genomes that represent 451 putative genera. Several of these phage genera are only distantly related to previously identified phages and are likely to found new families. Two of the candidate families, "Flandersviridae" and "Quimbyviridae", include some of the most common and abundant members of the human gut virome that infect Bacteroides, Parabacteroides, and Prevotella. The third proposed family, "Gratiaviridae," consists of less abundant phages that are distantly related to the families Autographiviridae, Drexlerviridae, and Chaseviridae. Analysis of CRISPR spacers indicates that phages of all three putative families infect bacteria of the phylum Bacteroidetes. Comparative genomic analysis of the three candidate phage families revealed features without precedent in phage genomes. Some "Quimbyviridae" phages possess Diversity-Generating Retroelements (DGRs) that generate hypervariable target genes nested within defense-related genes, whereas the previously known targets of phage-encoded DGRs are structural genes. Several "Flandersviridae" phages encode enzymes of the isoprenoid pathway, a lipid biosynthesis pathway that so far has not been known to be manipulated by phages. The "Gratiaviridae" phages encode a HipA-family protein kinase and glycosyltransferase, suggesting these phages modify the host cell wall, preventing superinfection by other phages. Hundreds of phages in these three and other families are shown to encode catalases and iron-sequestering enzymes that can be predicted to enhance cellular tolerance to reactive oxygen species.

CONCLUSIONS

Analysis of phage genomes identified in whole-community human gut metagenomes resulted in the delineation of at least three new candidate families of Caudovirales and revealed diverse putative mechanisms underlying phage-host interactions in the human gut. Addition of these phylogenetically classified, diverse, and distinct phages to public databases will facilitate taxonomic decomposition and functional characterization of human gut viromes. Video abstract.

Eco-evolutionary feedbacks mediated by bacterial membrane vesicles

Bacterial membrane vesicles (BMVs) are spherical extracellular organelles whose cargo is enclosed by a biological membrane. The cargo can be delivered to distant parts of a given habitat in a protected and concentrated manner. This review presents current knowledge about BMVs in the context of bacterial eco-evolutionary dynamics among different environments and hosts. BMVs may play an important role in establishing and stabilizing bacterial communities in such environments; for example, bacterial populations may benefit from BMVs to delay the negative effect of certain evolutionary trade-offs that can result in deleterious phenotypes. BMVs can also perform ecosystem engineering by serving as detergents, mediators in biochemical cycles, components of different biofilms, substrates for cross-feeding, defense systems against different dangers and enzyme-delivery mechanisms that can change substrate availability. BMVs further contribute to bacteria as mediators in different interactions, with either other bacterial species or their hosts. In short, BMVs extend and deliver phenotypic traits that can have ecological and evolutionary value to both their producers and the ecosystem as a whole.

4-Ethylphenol, A Volatile Organic Compound Produced by Disease-Resistant Soybean, Is a Potential Botanical Agrochemical Against Oomycetes

Oomycetes, represented by Phytophthora, are seriously harmful to agricultural production, resulting in a decline in grain quality and agricultural products and causing great economic losses. Integrated management of oomycete diseases is becoming more challenging, and plant derivatives represent effective alternatives to synthetic chemicals as novel crop protection solutions. Biologically active secondary metabolites are rapidly synthesized and released by plants in response to biotic stress caused by herbivores or insects, as well as pathogens. In this study, we identified groups of volatile organic compounds (VOCs) from soybean plants inoculated with Phytophthora sojae, the causal agent of soybean root rot. 4-Ethylphenol was present among the identified VOCs and was induced in the incompatible interaction between the plants and the pathogen. 4-Ethylphenol inhibited the growth of P. sojae and Phytophthora nicotianae and had toxicity to sporangia formation and zoospore germination by destroying the pathogen cell membrane; it had a good control effect on soybean root rot and tobacco black shank in the safe concentration range. Furthermore, 4-Ethylphenol had a potent antifungal activity against three soil-borne phytopathogenic fungi, Rhizoctonia solani, Fusarium graminearum, and Gaeumannomyces graminis var tritici, and four forma specialis of Fusarium oxysporum, which suggest a potential to be an eco-friendly biological control agent.

Deep phylogeny of cancer drivers and compensatory mutations

Driver mutations (DM) are the genetic impetus for most cancers. The DM are assumed to be deleterious in species evolution, being eliminated by purifying selection unless compensated by other mutations. We present deep phylogenies for 84 cancer driver genes and investigate the prevalence of 434 DM across gene-species trees. The DM are rare in species evolution, and 181 are completely absent, validating their negative fitness effect. The DM are more common in unicellular than in multicellular eukaryotes, suggesting a link between these mutations and cell proliferation control. 18 DM appear as the ancestral state in one or more major clades, including 3 among mammals. We identify within-gene, compensatory mutations for 98 DM and infer likely interactions between the DM and compensatory sites in protein structures. These findings elucidate the evolutionary status of DM and are expected to advance the understanding of the functions and evolution of oncogenes and tumor suppressors.

Rochman et al. present deep phylogenies for 84 cancer driver genes and examine the prevalence of driver mutations across gene-species trees. Their results show that driver mutations are rare in species evolution and give insight into the evolution of driver mutations and oncogenes.

A Conserved Glycoside Hydrolase Family 7 Cellobiohydrolase PsGH7a of Phytophthora sojae Is Required for Full Virulence on Soybean

Phytopathogens deploy glycoside hydrolases (GHs) to disintegrate plant cell walls for nutrition and invasion. However, the pathogenic mechanisms of the majority of GHs in virulence remain unknown, especially in oomycetes. In this study, a Phytophthora sojae gene encodes a GH7 family cellobiohydrolase, named PsGH7a, was identified. PsGH7a was highly induced during the cyst germination and infection stages. PsGH7a is conserved in oomycetes, and shares a high amino acid sequence identity (>85%) within Phytophthora genus. The recombinant PsGH7a catalyzes the hydrolysis of β-1,4-glucan and avicel, which represent the major components of cellulose in plant cell wall. The mutation of catalytic residue Glu236 to alanine resulted in a lower catalytic activity. In addition, the PsGH7a promotes Phytophthora invasion, while the mutant can not. Notably, PsGH7a protein triggers hypersensitive cell death in diverse plants. PsGH7a knockout mutants were generated via CRISPR/Cas9 system, to investigate its biological function. Compared to wild-type strain P6497, the mutants showed reduced virulence on susceptible soybean, indicates PsGH7a is indispensable to P. sojae virulence.

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.

A New Family of DNA Viruses Causing Disease in Crustaceans from Diverse Aquatic Biomes

Recent genomic and metagenomic studies have led to a dramatic expansion of the known diversity of nucleocytoplasmic large DNA viruses (NCLDVs) of eukaryotes, which include giant viruses of protists and important pathogens of vertebrates, such as poxviruses. However, the characterization of viruses from nonmodel hosts still lags behind. We sequenced the complete genomes of three viruses infecting crustaceans, the Caribbean spiny lobster, demon shrimp, and European shore crab. These viruses have the smallest genomes among the known NCLDVs, with losses of many core genes, some of which are shared with iridoviruses. The deterioration of the transcription apparatus is compatible with microscopic and ultrastructural observations indicating that these viruses replicate in the nucleus of infected cells rather than in the cytoplasm. Phylogenomic analysis indicates that these viruses are sufficiently distinct from all other NCLDVs to justify the creation of a separate family, for which we propose the name “Mininucleoviridae” (i.e., small viruses reproducing in the cell nucleus).

Panulirus argus virus 1 (PaV1) is the only known virus infecting the Caribbean spiny lobster (Panulirus argus) from the Caribbean Sea. Recently, related viruses, Dikerogammarus haemobaphes virus 1 (DhV1) and Carcinus maenas virus 1 (CmV1), have been detected in the demon shrimp (Dikerogammarus haemobaphes) and the European shore crab (Carcinus maenas), respectively, from sites in the United Kingdom. The virion morphology of these crustacean viruses is similar to that of iridoviruses. However, unlike iridoviruses and other nucleocytoplasmic large DNA viruses (NCLDVs), these viruses complete their morphogenesis in the host cell nucleus rather than in the cytoplasm. To date, these crustacean viruses have remained unclassified due to a lack of genomic data. Using an Illumina MiSeq sequencer, we sequenced the complete genomes of PaV1, CmV1, and DhV1. Comparative genome analysis shows that these crustacean virus genomes encode the 10 hallmark proteins previously described for the NCLDVs of eukaryotes, strongly suggesting that they are members of this group. With a size range of 70 to 74 kb, these are the smallest NCLDV genomes identified to date. Extensive gene loss, divergence of gene sequences, and the accumulation of low-complexity sequences reflect the extreme degradation of the genomes of these “minimal” NCLDVs rather than any direct relationship with the NCLDV ancestor. Phylogenomic analysis supports the classification of these crustacean viruses as a distinct family, “Mininucleoviridae,” within the pitho-irido-Marseille branch of the NCLDVs.

Adhesions Assemble!-Autoinhibition as a Major Regulatory Mechanism of Integrin-Mediated Adhesion

The advent of cell-cell and cell-extracellular adhesion enabled cells to interact in a coherent manner, forming larger structures and giving rise to the development of tissues, organs and complex multicellular life forms. The development of such organisms required tight regulation of dynamic adhesive structures by signaling pathways that coordinate cell attachment. Integrin-mediated adhesion to the extracellular matrix provides cells with support, survival signals and context-dependent cues that enable cells to run different cellular programs. One mysterious aspect of the process is how hundreds of proteins assemble seemingly spontaneously onto the activated integrin. An emerging concept is that adhesion assembly is regulated by autoinhibition of key proteins, a highly dynamic event that is modulated by a variety of signaling events. By enabling precise control of the activation state of proteins, autoinhibition enables localization of inactive proteins and the formation of pre-complexes. In response to the correct signals, these proteins become active and interact with other proteins, ultimately leading to development of cell-matrix junctions. Autoinhibition of key components of such adhesion complexes—including core components integrin, talin, vinculin, and FAK and important peripheral regulators such as RIAM, Src, and DLC1—leads to a view that the majority of proteins involved in complex assembly might be regulated by intramolecular interactions. Autoinhibition is relieved via multiple different signals including post-translation modification and proteolysis. More recently, mechanical forces have been shown to stabilize and increase the lifetimes of active conformations, identifying autoinhibition as a means of encoding mechanosensitivity. The complexity and scope for nuanced adhesion dynamics facilitated via autoinhibition provides numerous points of regulation. In this review, we discuss what is known about this mode of regulation and how it leads to rapid and tightly controlled assembly and disassembly of cell-matrix adhesion.

Unique structure and function of viral rhodopsins

Recently, two groups of rhodopsin genes were identified in large double-stranded DNA viruses. The structure and function of viral rhodopsins are unknown. We present functional characterization and high-resolution structure of an Organic Lake Phycodnavirus rhodopsin II (OLPVRII) of group 2. It forms a pentamer, with a symmetrical, bottle-like central channel with the narrow vestibule in the cytoplasmic part covered by a ring of 5 arginines, whereas 5 phenylalanines form a hydrophobic barrier in its exit. The proton donor E42 is placed in the helix B. The structure is unique among the known rhodopsins. Structural and functional data and molecular dynamics suggest that OLPVRII might be a light-gated pentameric ion channel analogous to pentameric ligand-gated ion channels, however, future patch clamp experiments should prove this directly. The data shed light on a fundamentally distinct branch of rhodopsins and may contribute to the understanding of virus-host interactions in ecologically important marine protists.

Genome and Environmental Activity of a Chrysochromulina parva Virus and Its Virophages

Some giant viruses are ecological agents that are predicted to be involved in the top-down control of single-celled eukaryotic algae populations in aquatic ecosystems. Despite an increased interest in giant viruses since the discovery and characterization of Mimivirus and other viral giants, little is known about their physiology and ecology. In this study, we characterized the genome and functional potential of a giant virus that infects the freshwater haptophyte Chrysochromulina parva, originally isolated from Lake Ontario. This virus, CpV-BQ2, is a member of the nucleo-cytoplasmic large DNA virus (NCLDV) group and possesses a 437 kb genome encoding 503 ORFs with a GC content of 25%. Phylogenetic analyses of core NCLDV genes place CpV-BQ2 amongst the emerging group of algae-infecting Mimiviruses informally referred to as the “extended Mimiviridae,” making it the first virus of this group to be isolated from a freshwater ecosystem. During genome analyses, we also captured and described the genomes of three distinct virophages that co-occurred with CpV-BQ2 and likely exploit CpV for their own replication. These virophages belong to the polinton-like viruses (PLV) group and encompass 19–23 predicted genes, including all of the core PLV genes as well as several genes implicated in genome modifications. We used the CpV-BQ2 and virophage reference sequences to recruit reads from available environmental metatranscriptomic data to estimate their activity in fresh waters. We observed moderate recruitment of both virus and virophage transcripts in samples obtained during Microcystis aeruginosa blooms in Lake Erie and Lake Tai, China in 2013, with a spike in activity in one sample. Virophage transcript abundance for two of the three isolates strongly correlated with that of the CpV-BQ2. Together, the results highlight the importance of giant viruses in the environment and establish a foundation for future research on the physiology and ecology CpV-BQ2 as a model system for algal Mimivirus dynamics in freshwaters.

Comparative genomics and evolution of trans-activating RNAs in Class 2 CRISPR-Cas systems

Trans-activating CRISPR (tracr) RNA is a distinct RNA species that interacts with the CRISPR (cr) RNA to form the dual guide (g) RNA in type II and subtype V-B CRISPR-Cas systems. The tracrRNA-crRNA interaction is essential for pre-crRNA processing as well as target recognition and cleavage. The tracrRNA consists of an antirepeat, which forms an imperfect hybrid with the repeat in the crRNA, and a distal region containing a Rho-independent terminator. Exhaustive comparative analysis of the sequences and predicted structures of the Class 2 CRISPR guide RNAs shows that all these guide RNAs share distinct structural features, in particular, the nexus stem-loop that separates the repeat-antirepeat hybrid from the distal portion of the tracrRNA and the conserved GU pair at that end of the hybrid. These structural constraints might ensure full exposure of the spacer for target recognition. Reconstruction of tracrRNA evolution for 4 tight bacterial groups demonstrates random drift of repeat-antirepeat complementarity within a window of hybrid stability that is, apparently, maintained by selection. An evolutionary scenario is proposed whereby tracrRNAs evolved on multiple occasions, via rearrangement of a CRISPR array to form the antirepeat in different locations with respect to the array. A functional tracrRNA would form if, in the new location, the antirepeat is flanked by sequences that meet the minimal requirements for a promoter and a Rho-independent terminator. Alternatively, or additionally, the antirepeat sequence could be occasionally 'reset' by recombination with a repeat, restoring the functionality of tracrRNAs that drift beyond the required minimal hybrid stability.

Virus Genomes from Deep Sea Sediments Expand the Ocean Megavirome and Support Independent Origins of Viral Gigantism

Genomics and evolution of giant viruses are two of the most vigorously developing areas of virus research. Lately, metagenomics has become the main source of new virus genomes. Here we describe a metagenomic analysis of the genomes of large and giant viruses from deep sea sediments. The assembled new virus genomes substantially expand the known diversity of the nucleocytoplasmic large DNA viruses of eukaryotes. The results support the concept of independent evolution of giant viruses from smaller ancestors in different virus branches.

The nucleocytoplasmic large DNA viruses (NCLDV) of eukaryotes (proposed order, “Megavirales”) include the families Poxviridae, Asfarviridae, Iridoviridae, Ascoviridae, Phycodnaviridae, Marseilleviridae, and Mimiviridae, as well as still unclassified pithoviruses, pandoraviruses, molliviruses, and faustoviruses. Several of these virus groups include giant viruses, with genome and particle sizes exceeding those of many bacterial and archaeal cells. We explored the diversity of the NCLDV in deep sea sediments from the Loki’s Castle hydrothermal vent area. Using metagenomics, we reconstructed 23 high-quality genomic bins of novel NCLDV, 15 of which are related to pithoviruses, 5 to marseilleviruses, 1 to iridoviruses, and 2 to klosneuviruses. Some of the identified pithovirus-like and marseillevirus-like genomes belong to deep branches in the phylogenetic tree of core NCLDV genes, substantially expanding the diversity and phylogenetic depth of the respective groups. The discovered viruses, including putative giant members of the family Marseilleviridae, have a broad range of apparent genome sizes, in agreement with the multiple, independent origins of gigantism in different branches of the NCLDV. Phylogenomic analysis reaffirms the monophyly of the pithovirus-iridovirus-marseillevirus branch of the NCLDV. Similarly to other giant viruses, the pithovirus-like viruses from Loki’s Castle encode translation systems components. Phylogenetic analysis of these genes indicates a greater bacterial contribution than had been detected previously. Genome comparison suggests extensive gene exchange between members of the pithovirus-like viruses and Mimiviridae. Further exploration of the genomic diversity of Megavirales in additional sediment samples is expected to yield new insights into the evolution of giant viruses and the composition of the ocean megavirome.

Catalytic linkage between caspase activity and proteostasis in Archaea

The model haloarchaeon, Haloferax volcanii possess an extremely high, and highly specific, basal caspase activity in exponentially growing cells that closely resembles caspase-4. This activity is specifically inhibited by the pan-caspase inhibitor, z-VAD-FMK, and has no cross-reactivity with other known protease families. Although it is one of the dominant cellular proteolytic activities in exponentially growing H. volcanii cells, the interactive cellular roles remain unknown and the protein(s) responsible for this activity remain elusive. Here, biochemical purification and in situ trapping with caspase targeted covalent inhibitors combined with genome-enabled proteomics, structural analysis, targeted gene knockouts and treatment with canavanine demonstrated a catalytic linkage between caspase activity and thermosomes, proteasomes and cdc48b, a cell division protein and proteasomal degradation facilitating ATPase, as part of an 'interactase' of stress-related protein complexes with an established link to the unfolded protein response (UPR). Our findings provide novel cellular and biochemical context for the observed caspase activity in Archaea and add new insight to understanding the role of this activity, implicating their possible role in the establishment of protein stress and ER associated degradation pathways in Eukarya.

Extreme Deviations from Expected Evolutionary Rates in Archaeal Protein Families

Origin of new biological functions is a complex phenomenon ranging from single-nucleotide substitutions to the gain of new genes via horizontal gene transfer or duplication. Neofunctionalization and subfunctionalization of proteins is often attributed to the emergence of paralogs that are subject to relaxed purifying selection or positive selection and thus evolve at accelerated rates. Such phenomena potentially could be detected as anomalies in the phylogenies of the respective gene families. We developed a computational pipeline to search for such anomalies in 1,834 orthologous clusters of archaeal genes, focusing on lineage-specific subfamilies that significantly deviate from the expected rate of evolution. Multiple potential cases of neofunctionalization and subfunctionalization were identified, including some ancient, house-keeping gene families, such as ribosomal protein S10, general transcription factor TFIIB and chaperone Hsp20. As expected, many cases of apparent acceleration of evolution are associated with lineage-specific gene duplication. On other occasions, long branches in phylogenetic trees correspond to horizontal gene transfer across long evolutionary distances. Significant deceleration of evolution is less common than acceleration, and the underlying causes are not well understood; functional shifts accompanied by increased constraints could be involved. Many gene families appear to be “highly evolvable,” that is, include both long and short branches. Even in the absence of precise functional predictions, this approach allows one to select targets for experimentation in search of new biology.

Evolutionary insights into Trm112-methyltransferase holoenzymes involved in translation between archaea and eukaryotes

Protein synthesis is a complex and highly coordinated process requiring many different protein factors as well as various types of nucleic acids. All translation machinery components require multiple maturation events to be functional. These include post-transcriptional and post-translational modification steps and methylations are the most frequent among these events. In eukaryotes, Trm112, a small protein (COG2835) conserved in all three domains of life, interacts and activates four methyltransferases (Bud23, Trm9, Trm11 and Mtq2) that target different components of the translation machinery (rRNA, tRNAs, release factors). To clarify the function of Trm112 in archaea, we have characterized functionally and structurally its interaction network using Haloferax volcanii as model system. This led us to unravel that methyltransferases are also privileged Trm112 partners in archaea and that this Trm112 network is much more complex than anticipated from eukaryotic studies. Interestingly, among the identified enzymes, some are functionally orthologous to eukaryotic Trm112 partners, emphasizing again the similarity between eukaryotic and archaeal translation machineries. Other partners display some similarities with bacterial methyltransferases, suggesting that Trm112 is a general partner for methyltransferases in all living organisms.

Archaea-First and the Co-Evolutionary Diversification of Domains of Life

The origins and evolution of the Archaea, Bacteria, and Eukarya remain controversial. Phylogenomic-wide studies of molecular features that are evolutionarily conserved, such as protein structural domains, suggest Archaea is the first domain of life to diversify from a stem line of descent. This line embodies the last universal common ancestor of cellular life. Here, we propose that ancestors of Euryarchaeota co-evolved with those of Bacteria prior to the diversification of Eukarya. This co-evolutionary scenario is supported by comparative genomic and phylogenomic analyses of the distributions of fold families of domains in the proteomes of free-living organisms, which show horizontal gene recruitments and informational process homologies. It also benefits from the molecular study of cell physiologies responsible for membrane phospholipids, methanogenesis, methane oxidation, cell division, gas vesicles, and the cell wall. Our theory however challenges popular cell fusion and two-domain of life scenarios derived from sequence analysis, demanding phylogenetic reconciliation. Also see the video abstract here: https://youtu.be/9yVWn_Q9faY.

The Rhizome of Lokiarchaeota Illustrates the Mosaicity of Archaeal Genomes

Genome remodeling and exchange of sequences are widespread in the prokaryotic world and mosaic genomes challenge the classification of prokaryotes, which cannot be properly achieved in terms of a single gene or group of genes. Here, we studied individually the gene collection of the archaic microorganism Lokiarchaeum sp., suggested as an archaeal host close to the emergence of the eukaryotes. The network or rhizome of all Lokiarchaeum sp. genes revealed that the genomic repertoire is mainly composed of genes from archaeal (∼36%) and bacterial origin (∼28%), distantly followed by components of eukaryotic origin (∼2%). Thirty-three percent of genes were unique to this species (ORFans). The mosaicity of archaea was also supported by studying Methanomassiliicoccus luminyensis, an archaea from the gut, in which 67% of the genomic repertoire arised from archaea and 22% from bacteria. Our results illustrate the intricate evolutionary relationships of the archaeal genome repertoire and highlight the rhizome-like processes of evolution in archaea, their mosaicity, and chimeric origin composed of different domains of life, questioning the reality of a tree of life.

Systematic prediction of genes functionally linked to CRISPR-Cas systems by gene neighborhood analysis

The CRISPR-Cas systems of bacterial and archaeal adaptive immunity consist of direct repeat arrays separated by unique spacers and multiple CRISPR-associated (cas) genes encoding proteins that mediate all stages of the CRISPR response. In addition to the relatively small set of core cas genes that are typically present in all CRISPR-Cas systems of a given (sub)type and are essential for the defense function, numerous genes occur in CRISPR-cas loci only sporadically. Some of these have been shown to perform various ancillary roles in CRISPR response, but the functional relevance of most remains unknown. We developed a computational strategy for systematically detecting genes that are likely to be functionally linked to CRISPR-Cas. The approach is based on a "CRISPRicity" metric that measures the strength of CRISPR association for all protein-coding genes from sequenced bacterial and archaeal genomes. Uncharacterized genes with CRISPRicity values comparable to those of cas genes are considered candidate CRISPR-linked genes. We describe additional criteria to predict functionally relevance for genes in the candidate set and identify 79 genes as strong candidates for functional association with CRISPR-Cas systems. A substantial majority of these CRISPR-linked genes reside in type III CRISPR-cas loci, which implies exceptional functional versatility of type III systems. Numerous candidate CRISPR-linked genes encode integral membrane proteins suggestive of tight membrane association of CRISPR-Cas systems, whereas many others encode proteins implicated in various signal transduction pathways. These predictions provide ample material for improving annotation of CRISPR-cas loci and experimental characterization of previously unsuspected aspects of CRISPR-Cas system functionality.

Comparative in Silico Analysis of Ferric Reduction Oxidase (FRO) Genes Expression Patterns in Response to Abiotic Stresses, Metal and Hormone Applications

The ferric reduction oxidase (FRO) gene family is involved in various biological processes widely found in plants and may play an essential role in metal homeostasis, tolerance and intricate signaling networks in response to a number of abiotic stresses. Our study describes the identification, characterization and evolutionary relationships of FRO genes families. Here, total 50 FRO genes in Plantae and 15 ‘FRO like’ genes in non-Plantae were retrieved from 16 different species. The entire FRO genes have been divided into seven clades according to close similarity in biological and functional behavior. Three conserved domains were common in FRO genes while in two FROs sub genome have an extra NADPH-Ox domain, separating the function of plant FROs. OsFRO1 and OsFRO7 genes were expressed constitutively in rice plant. Real-time RT-PCR analysis demonstrated that the expression of OsFRO1 was high in flag leaf, and OsFRO7 gene expression was maximum in leaf blade and flag leaf. Both genes showed vigorous expressions level in response to different abiotic and hormones treatments. Moreover, the expression of both genes was also substantial under heavy metal stresses. OsFRO1 gene expression was triggered following 6 h under Zn, Pb, Co and Ni treatments, whereas OsFRO7 gene expression under Fe, Pb and Ni after 12 h, Zn and Cr after 6 h, and Mn and Co after 3 h treatments. These findings suggest the possible involvement of both the genes under abiotic and metal stress and the regulation of phytohormones. Therefore, our current work may provide the foundation for further functional characterization of rice FRO genes family.

Vast diversity of prokaryotic virus genomes encoding double jelly-roll major capsid proteins uncovered by genomic and metagenomic sequence analysis

BACKGROUND

Analysis of metagenomic sequences has become the principal approach for the study of the diversity of viruses. Many recent, extensive metagenomic studies on several classes of viruses have dramatically expanded the visible part of the virosphere, showing that previously undetected viruses, or those that have been considered rare, actually are important components of the global virome.

RESULTS

We investigated the provenance of viruses related to tail-less bacteriophages of the family Tectiviridae by searching genomic and metagenomics sequence databases for distant homologs of the tectivirus-like Double Jelly-Roll major capsid proteins (DJR MCP). These searches resulted in the identification of numerous genomes of virus-like elements that are similar in size to tectiviruses (10-15 kilobases) and have diverse gene compositions. By comparison of the gene repertoires, the DJR MCP-encoding genomes were classified into 6 distinct groups that can be predicted to differ in reproduction strategies and host ranges. Only the DJR MCP gene that is present by design is shared by all these genomes, and most also encode a predicted DNA-packaging ATPase; the rest of the genes are present only in subgroups of this unexpectedly diverse collection of DJR MCP-encoding genomes. Only a minority encode a DNA polymerase which is a hallmark of the family Tectiviridae and the putative family "Autolykiviridae". Notably, one of the identified putative DJR MCP viruses encodes a homolog of Cas1 endonuclease, the integrase involved in CRISPR-Cas adaptation and integration of transposon-like elements called casposons. This is the first detected occurrence of Cas1 in a virus. Many of the identified elements are individual contigs flanked by inverted or direct repeats and appear to represent complete, extrachromosomal viral genomes, whereas others are flanked by bacterial genes and thus can be considered as proviruses. These contigs come from metagenomes of widely different environments, some dominated by archaea and others by bacteria, suggesting that collectively, the DJR MCP-encoding elements have a broad host range among prokaryotes.

CONCLUSIONS

The findings reported here greatly expand the known host range of (putative) viruses of bacteria and archaea that encode a DJR MCP. They also demonstrate the extreme diversity of genome architectures in these viruses that encode no universal proteins other than the capsid protein that was used as the marker for their identification. From a supposedly minor group of bacterial and archaeal viruses, these viruses are emerging as a substantial component of the prokaryotic virome.

Giant viruses with an expanded complement of translation system components

The discovery of giant viruses blurred the sharp division between viruses and cellular life. Giant virus genomes encode proteins considered as signatures of cellular organisms, particularly translation system components, prompting hypotheses that these viruses derived from a fourth domain of cellular life. Here we report the discovery of a group of giant viruses (Klosneuviruses) in metagenomic data. Compared with other giant viruses, the Klosneuviruses encode an expanded translation machinery, including aminoacyl transfer RNA synthetases with specificities for all 20 amino acids. Notwithstanding the prevalence of translation system components, comprehensive phylogenomic analysis of these genes indicates that Klosneuviruses did not evolve from a cellular ancestor but rather are derived from a much smaller virus through extensive gain of host genes.

Proposed Role for KaiC-Like ATPases as Major Signal Transduction Hubs in Archaea

All organisms must adapt to ever-changing environmental conditions and accordingly have evolved diverse signal transduction systems. In bacteria, the most abundant networks are built around the two-component signal transduction systems that include histidine kinases and receiver domains. In contrast, eukaryotic signal transduction is dominated by serine/threonine/tyrosine protein kinases. Both of these systems are also found in archaea, but they are not as common and diversified as their bacterial and eukaryotic counterparts, suggesting the possibility that archaea have evolved other, still uncharacterized signal transduction networks. Here we propose a role for KaiC family ATPases, known to be key components of the circadian clock in cyanobacteria, in archaeal signal transduction. The KaiC family is notably expanded in most archaeal genomes, and although most of these ATPases remain poorly characterized, members of the KaiC family have been shown to control archaellum assembly and have been found to be a stable component of the gas vesicle system in Halobacteria Computational analyses described here suggest that KaiC-like ATPases and their homologues with inactivated ATPase domains are involved in many other archaeal signal transduction pathways and comprise major hubs of complex regulatory networks. We predict numerous input and output domains that are linked to KaiC-like proteins, including putative homologues of eukaryotic DEATH domains that could function as adapters in archaeal signaling networks. We further address the relationships of the archaeal family of KaiC homologues to the bona fide KaiC of cyanobacteria and implications for the existence of a KaiC-based circadian clock apparatus in archaea.IMPORTANCE Little is currently known about signal transduction pathways in Archaea Recent studies indicate that KaiC-like ATPases, known as key components of the circadian clock apparatus in cyanobacteria, are involved in the regulation of archaellum assembly and, likely, type IV pili and the gas vesicle system in Archaea We performed comprehensive comparative genomic analyses of the KaiC family. A vast protein interaction network was revealed, with KaiC family proteins as hubs for numerous input and output components, many of which are shared with two-component signal transduction systems. Putative KaiC-based signal transduction systems are predicted to regulate the activities of membrane-associated complexes and individual proteins, such as signal recognition particle and membrane transporters, and also could be important for oxidative stress response regulation. KaiC-centered signal transduction networks are predicted to play major roles in archaeal physiology, and this work is expected to stimulate their experimental characterization.

Discovery of an expansive bacteriophage family that includes the most abundant viruses from the human gut

Metagenomic sequence analysis is rapidly becoming the primary source of virus discovery ^1–3. A substantial majority of the currently available virus genomes comes from metagenomics, and some of these represent extremely abundant viruses even if never grown in the laboratory. A particularly striking case of a virus discovered via metagenomics is crAssphage, which is by far the most abundant human-associated virus known, comprising up to 90% of the sequences in the gut virome ⁴. Over 80% of the predicted proteins encoded in the approximately 100 kilobase crAssphage genome showed no significant similarity to available protein sequences, precluding classification of this virus and hampering further study. Here we combine comprehensive search of genomic and metagenomic databases with sensitive methods for protein sequence analysis to identify an expansive, diverse group of bacteriophages related to crAssphage and predict the functions of the majority of phage proteins, in particular, those that comprise the structural, replication and expression modules. Most if not all of the crAss-like phages appear to be associated with diverse bacteria from the phylum Bacteroidetes, which includes some of the most abundant bacteria in the human gut microbiome and are also common in various other habitats. These findings provide for experimental characterization of the most abundant but poorly understood members of the human-associated virome.

A proteome view of structural, functional, and taxonomic characteristics of major protein domain clusters

Proteome-scale bioinformatics research is increasingly conducted as the number of completely sequenced genomes increases, but analysis of protein domains (PDs) usually relies on similarity in their amino acid sequences and/or three-dimensional structures. Here, we present results from a bi-clustering analysis on presence/absence data for 6,580 unique PDs in 2,134 species with a sequenced genome, thus covering a complete set of proteins, for the three superkingdoms of life, Bacteria, Archaea, and Eukarya. Our analysis revealed eight distinctive PD clusters, which, following an analysis of enrichment of Gene Ontology functions and CATH classification of protein structures, were shown to exhibit structural and functional properties that are taxa-characteristic. For examples, the largest cluster is ubiquitous in all three superkingdoms, constituting a set of 1,472 persistent domains created early in evolution and retained in living organisms and characterized by basic cellular functions and ancient structural architectures, while an Archaea and Eukarya bi-superkingdom cluster suggests its PDs may have existed in the ancestor of the two superkingdoms, and others are single superkingdom- or taxa (e.g. Fungi)-specific. These results contribute to increase our appreciation of PD diversity and our knowledge of how PDs are used in species, yielding implications on species evolution.

Sequence complementarity at the ribosomal Peptidyl Transferase Centre implies self-replicating origin

A feasible scenario for the emergence of life requires the spontaneous materialization and sustainability of a proto-ribosome that could have catalysed the formation of the first peptides. Models of proto-ribosomes were derived from the ribosomal Peptidyl Transferase Centre (PTC) region, but the poor prebiotic copying abilities give rise to the question of their mode of replication. Here, complementarity is demonstrated in bacterial ribosomes, between nucleotides that constitute the two halves of the PTC cavity. The complementarity corroborates the dimeric nature of the proto-ribosome and is likely to underlie the symmetry of the PTC region. Furthermore, it indicates a simple and efficient replication mode; the strand of each monomer could have acted as a template for the synthesis of its counterpart, forming a self-replicating ribozyme.

Characterization of Reconstructed Ancestral Proteins Suggests a Change in Temperature of the Ancient Biosphere

Understanding the evolution of ancestral life, and especially the ability of some organisms to flourish in the variable environments experienced in Earth’s early biosphere, requires knowledge of the characteristics and the environment of these ancestral organisms. Information about early life and environmental conditions has been obtained from fossil records and geological surveys. Recent advances in phylogenetic analysis, and an increasing number of protein sequences available in public databases, have made it possible to infer ancestral protein sequences possessed by ancient organisms. However, the in silico studies that assess the ancestral base content of ribosomal RNAs, the frequency of each amino acid in ancestral proteins, and estimate the environmental temperatures of ancient organisms, show conflicting results. The characterization of ancestral proteins reconstructed in vitro suggests that ancient organisms had very thermally stable proteins, and therefore were thermophilic or hyperthermophilic. Experimental data supports the idea that only thermophilic ancestors survived the catastrophic increase in temperature of the biosphere that was likely associated with meteorite impacts during the early history of Earth. In addition, by expanding the timescale and including more ancestral proteins for reconstruction, it appears as though the Earth’s surface temperature gradually decreased over time, from Archean to present.

Symbiosis in eukaryotic evolution

Fifty years ago, Lynn Margulis, inspiring in early twentieth-century ideas that put forward a symbiotic origin for some eukaryotic organelles, proposed a unified theory for the origin of the eukaryotic cell based on symbiosis as evolutionary mechanism. Margulis was profoundly aware of the importance of symbiosis in the natural microbial world and anticipated the evolutionary significance that integrated cooperative interactions might have as mechanism to increase cellular complexity. Today, we have started fully appreciating the vast extent of microbial diversity and the importance of syntrophic metabolic cooperation in natural ecosystems, especially in sediments and microbial mats. Also, not only the symbiogenetic origin of mitochondria and chloroplasts has been clearly demonstrated, but improvement in phylogenomic methods combined with recent discoveries of archaeal lineages more closely related to eukaryotes further support the symbiogenetic origin of the eukaryotic cell. Margulis left us in legacy the idea of 'eukaryogenesis by symbiogenesis'. Although this has been largely verified, when, where, and specifically how eukaryotic cells evolved are yet unclear. Here, we shortly review current knowledge about symbiotic interactions in the microbial world and their evolutionary impact, the status of eukaryogenetic models and the current challenges and perspectives ahead to reconstruct the evolutionary path to eukaryotes.

Ancient, highly conserved proteins from a LUCA with complex cell biology provide evidence in support of the nuclear compartment commonality (NuCom) hypothesis

The nuclear compartment commonality (NuCom) hypothesis posits a complex last common ancestor (LUCA) with membranous compartments including a nuclear membrane. Such a LUCA then evolved to produce two nucleated lineages of the tree of life: the Planctomycetes-Verrucomicrobia-Chlamydia superphylum (PVC) within the Bacteria, and the Eukarya. We propose that a group of ancient essential protokaryotic signature proteins (PSPs) originating in LUCA were incorporated into ancestors of PVC Bacteria and Eukarya. Tubulins, ubiquitin system enzymes and sterol-synthesizing enzymes are consistent with early origins of these features shared between the PVC superphylum and Eukarya.

The mother of all actins?

New insights into the structure of filaments made of crenactin, a homolog of actin found in archaea, shed light on how the cytoskeleton might have evolved.

Could a Proto-Ribosome Emerge Spontaneously in the Prebiotic World?

An indispensable prerequisite for establishing a scenario of life emerging by natural processes is the requirement that the first simple proto-molecules could have had a realistic probability of self-assembly from random molecular polymers in the prebiotic world. The vestige of the proto-ribosome, which is believed to be still embedded in the contemporary ribosome, is used to assess the feasibility of such spontaneous emergence. Three concentric structural elements of different magnitudes, having a dimeric nature derived from the symmetrical region of the ribosomal large subunit, were suggested to constitute the vestige of the proto-ribosome. It is assumed to have materialized spontaneously in the prebiotic world, catalyzing non-coded peptide bond formation and simple elongation. Probabilistic and energetic considerations are applied in order to evaluate the suitability of the three contenders for being the initial proto-ribosome. The analysis points to the simplest proto-ribosome, comprised of a dimer of tRNA-like molecules presently embedded in the core of the symmetrical region, as the only one having a realistic statistical likelihood of spontaneous emergence from random RNA chains. Hence it offers a feasible starting point for a continuous evolutionary path from the prebiotic matter, through natural processes, into the intricate modern translation system.

Quest for Ancestors of Eukaryal Cells Based on Phylogenetic Analyses of Aminoacyl-tRNA Synthetases

The three-domain phylogenetic system of life has been challenged, particularly with regard to the position of Eukarya. The recent increase of known genome sequences has allowed phylogenetic analyses of all extant organisms using concatenated sequence alignment of universally conserved genes; these data supported the two-domain hypothesis, which place eukaryal species as ingroups of the Domain Archaea. However, the origin of Eukarya is complicated: the closest archaeal species to Eukarya differs in single-gene phylogenetic analyses depending on the genes. In this report, we performed molecular phylogenetic analyses of 23 aminoacyl-tRNA synthetases (ARS). Cytoplasmic ARSs in 12 trees showed a monophyletic Eukaryotic branch. One ARS originated from TACK superphylum. One ARS originated from Euryarchaeota and three originated from DPANN superphylum. Four ARSs originated from different bacterial species. The other 8 cytoplasmic ARSs were split into two or three groups in respective trees, which suggested that the cytoplasmic ARSs were replaced by secondary ARSs, and the original ARSs have been lost during evolution of Eukarya. In these trees, one original cytoplasmic ARS was derived from Euryarchaeota and three were derived from DPANN superphylum. Our results strongly support the two-domain hypothesis. We discovered that rampant-independent lateral gene transfers from several archaeal species of DPANN superphylum have contributed to the formation of Eukaryal cells. Based on our phylogenetic analyses, we proposed a model for the establishment of Eukarya.

Origin of eukaryotes from within archaea, archaeal eukaryome and bursts of gene gain: eukaryogenesis just made easier?

The origin of eukaryotes is a fundamental, forbidding evolutionary puzzle. Comparative genomic analysis clearly shows that the last eukaryotic common ancestor (LECA) possessed most of the signature complex features of modern eukaryotic cells, in particular the mitochondria, the endomembrane system including the nucleus, an advanced cytoskeleton and the ubiquitin network. Numerous duplications of ancestral genes, e.g. DNA polymerases, RNA polymerases and proteasome subunits, also can be traced back to the LECA. Thus, the LECA was not a primitive organism and its emergence must have resulted from extensive evolution towards cellular complexity. However, the scenario of eukaryogenesis, and in particular the relationship between endosymbiosis and the origin of eukaryotes, is far from being clear. Four recent developments provide new clues to the likely routes of eukaryogenesis. First, evolutionary reconstructions suggest complex ancestors for most of the major groups of archaea, with the subsequent evolution dominated by gene loss. Second, homologues of signature eukaryotic proteins, such as actin and tubulin that form the core of the cytoskeleton or the ubiquitin system, have been detected in diverse archaea. The discovery of this ‘dispersed eukaryome’ implies that the archaeal ancestor of eukaryotes was a complex cell that might have been capable of a primitive form of phagocytosis and thus conducive to endosymbiont capture. Third, phylogenomic analyses converge on the origin of most eukaryotic genes of archaeal descent from within the archaeal evolutionary tree, specifically, the TACK superphylum. Fourth, evidence has been presented that the origin of the major archaeal phyla involved massive acquisition of bacterial genes. Taken together, these findings make the symbiogenetic scenario for the origin of eukaryotes considerably more plausible and the origin of the organizational complexity of eukaryotic cells more readily explainable than they appeared until recently.

Diversity and Evolution of Type IV pili Systems in Archaea

Many surface structures in archaea including various types of pili and the archaellum (archaeal flagellum) are homologous to bacterial type IV pili systems (T4P). The T4P consist of multiple proteins, often with poorly conserved sequences, complicating their identification in sequenced genomes. Here we report a comprehensive census of T4P encoded in archaeal genomes using sensitive methods for protein sequence comparison. This analysis confidently identifies as T4P components about 5000 archaeal gene products, 56% of which are currently annotated as hypothetical in public databases. Combining results of this analysis with a comprehensive comparison of genomic neighborhoods of the T4P, we present models of organization of 10 most abundant variants of archaeal T4P. In addition to the differentiation between major and minor pilins, these models include extra components, such as S-layer proteins, adhesins and other membrane and intracellular proteins. For most of these systems, dedicated major pilin families are identified including numerous stand alone major pilin genes of the PilA family. Evidence is presented that secretion ATPases of the T4P and cognate TadC proteins can interact with different pilin sets. Modular evolution of T4P results in combinatorial variability of these systems. Potential regulatory or modulating proteins for the T4P are identified including KaiC family ATPases, vWA domain-containing proteins and the associated MoxR/GvpN ATPase, TFIIB homologs and multiple unrelated transcription regulators some of which are associated specific T4P. Phylogenomic analysis suggests that at least one T4P system was present in the last common ancestor of the extant archaea. Multiple cases of horizontal transfer and lineage-specific duplication of T4P loci were detected. Generally, the T4P of the archaeal TACK superphylum are more diverse and evolve notably faster than those of euryarchaea. The abundance and enormous diversity of T4P in hyperthermophilic archaea present a major enigma. Apparently, fundamental aspects of the biology of hyperthermophiles remain to be elucidated.

Chimeric Filoviruses for Identification and Characterization of Monoclonal Antibodies

Recent experiments suggest that some glycoprotein (GP)-specific monoclonal antibodies (MAbs) can protect experimental animals against the filovirus Ebola virus (EBOV). There is a need for isolation of MAbs capable of neutralizing multiple filoviruses. Antibody neutralization assays for filoviruses frequently use surrogate systems such as the rhabdovirus vesicular stomatitis Indiana virus (VSV), lentiviruses or gammaretroviruses with their envelope proteins replaced with EBOV GP or pseudotyped with EBOV GP. It is optimal for both screening and in-depth characterization of newly identified neutralizing MAbs to generate recombinant filoviruses that express a reporter fluorescent protein in order to more easily monitor and quantify the infection. Our study showed that unlike neutralization-sensitive chimeric VSV, authentic filoviruses are highly resistant to neutralization by MAbs. We used reverse genetics techniques to replace EBOV GP with its counterpart from the heterologous filoviruses Bundibugyo virus (BDBV), Sudan virus, and even Marburg virus and Lloviu virus, which belong to the heterologous genera in the filovirus family. This work resulted in generation of multiple chimeric filoviruses, demonstrating the ability of filoviruses to tolerate swapping of the envelope protein. The sensitivity of chimeric filoviruses to neutralizing MAbs was similar to that of authentic biologically derived filoviruses with the same GP. Moreover, disabling the expression of the secreted GP (sGP) resulted in an increased susceptibility of an engineered virus to the BDBV52 MAb isolated from a BDBV survivor, suggesting a role for sGP in evasion of antibody neutralization in the context of a human filovirus infection.

IMPORTANCE

The study demonstrated that chimeric rhabdoviruses in which G protein is replaced with filovirus GP, widely used as surrogate targets for characterization of filovirus neutralizing antibodies, do not accurately predict the ability of antibodies to neutralize authentic filoviruses, which appeared to be resistant to neutralization. However, a recombinant EBOV expressing a fluorescent protein tolerated swapping of GP with counterparts from heterologous filoviruses, allowing high-throughput screening of B cell lines to isolate MAbs of any filovirus specificity. Human MAb BDBV52, which was isolated from a survivor of BDBV infection, was capable of partially neutralizing a chimeric EBOV carrying BDBV GP in which expression of sGP was disabled. In contrast, the parental virus expressing sGP was resistant to the MAb. Thus, the ability of filoviruses to tolerate swapping of GP can be used for identification of neutralizing MAbs specific to any filovirus and for the characterization of MAb specificity and mechanism of action.

Evolution of replication machines

The machines that decode and regulate genetic information require the translation, transcription and replication pathways essential to all living cells. Thus, it might be expected that all cells share the same basic machinery for these pathways that were inherited from the primordial ancestor cell from which they evolved. A clear example of this is found in the translation machinery that converts RNA sequence to protein. The translation process requires numerous structural and catalytic RNAs and proteins, the central factors of which are homologous in all three domains of life, bacteria, archaea and eukarya. Likewise, the central actor in transcription, RNA polymerase, shows homology among the catalytic subunits in bacteria, archaea and eukarya. In contrast, while some "gears" of the genome replication machinery are homologous in all domains of life, most components of the replication machine appear to be unrelated between bacteria and those of archaea and eukarya. This review will compare and contrast the central proteins of the "replisome" machines that duplicate DNA in bacteria, archaea and eukarya, with an eye to understanding the issues surrounding the evolution of the DNA replication apparatus.

Archaeal ancestors of eukaryotes: not so elusive any more

The origin of eukaryotes is one of the hardest problems in evolutionary biology and sometimes raises the ominous specter of irreducible complexity. Reconstruction of the gene repertoire of the last eukaryotic common ancestor (LECA) has revealed a highly complex organism with a variety of advanced features but no detectable evolutionary intermediates to explain their origin. Recently, however, genome analysis of diverse archaea led to the discovery of apparent ancestral versions of several signature eukaryotic systems, such as the actin cytoskeleton and the ubiquitin network, that are scattered among archaea. These findings inspired the hypothesis that the archaeal ancestor of eukaryotes was an unusually complex form with an elaborate intracellular organization. The latest striking discovery made by deep metagenomic sequencing vindicates this hypothesis by showing that in phylogenetic trees eukaryotes fall within a newly identified archaeal group, the Lokiarchaeota, which combine several eukaryotic signatures previously identified in different archaea. The discovery of complex archaea that are the closest living relatives of eukaryotes is most compatible with the symbiogenetic scenario for eukaryogenesis.

Evolution of plant δ(1)-pyrroline-5-carboxylate reductases from phylogenetic and structural perspectives

Proline plays a crucial role in cell growth and stress responses, and its accumulation is essential for the tolerance of adverse environmental conditions in plants. Two routes are used to biosynthesize proline in plants. The main route uses glutamate as a precursor, while in the other route proline is derived from ornithine. The terminal step of both pathways, the conversion of δ(1)-pyrroline-5-carboxylate (P5C) to L-proline, is catalyzed by P5C reductase (P5CR) using NADH or NADPH as a cofactor. Since P5CRs are important housekeeping enzymes, they are conserved across all domains of life and appear to be relatively unaffected throughout evolution. However, global analysis of these enzymes unveiled significant functional diversity in the preference for cofactors (NADPH vs. NADH), variation in metal dependence and the differences in the oligomeric state. In our study we investigated evolutionary patterns through phylogenetic and structural analysis of P5CR representatives from all kingdoms of life, with emphasis on the plant species. We also attempted to correlate local sequence/structure variation among the functionally and structurally characterized members of the family.