RNA-guided RNA silencing by an Asgard archaeal Argonaute.

Argonaute proteins are the central effectors of RNA-guided RNA silencing pathways in eukaryotes, playing crucial roles in gene repression and defense against viruses and transposons. Eukaryotic Argonautes are subdivided into two clades: AGOs generally facilitate miRNA- or siRNA-mediated silencing, while PIWIs generally facilitate piRNA-mediated silencing. It is currently unclear when and how Argonaute-based RNA silencing mechanisms arose and diverged during the emergence and early evolution of eukaryotes. Here, we show that in Asgard archaea, the closest prokaryotic relatives of eukaryotes, an evolutionary expansion of Argonaute proteins took place. In particular, a deep-branching PIWI protein (HrAgo1) encoded by the genome of the Lokiarchaeon 'Candidatus Harpocratesius repetitus' shares a common origin with eukaryotic PIWI proteins. Contrasting known prokaryotic Argonautes that use single-stranded DNA as guides and/or targets, HrAgo1 mediates RNA-guided RNA cleavage, and facilitates gene silencing when expressed in human cells and supplied with miRNA precursors. A cryo-EM structure of HrAgo1, combined with quantitative single-molecule experiments, reveals that the protein displays structural features and target-binding modes that are a mix of those of eukaryotic AGO and PIWI proteins. Thus, this deep-branching archaeal PIWI may have retained an ancestral molecular architecture that preceded the functional and mechanistic divergence of eukaryotic AGOs and PIWIs.

Description of Asgardarchaeum abyssi gen. nov. spec. nov., a novel species within the class Asgardarchaeia and phylum Asgardarchaeota in accordance with the SeqCode.

Asgardarchaeota, commonly referred to as Asgard archaea, is a candidatus phylum-rank archaeal clade that includes the closest archaeal relatives of eukaryotes. Despite their prevalence in the scientific literature, the name Asgardarchaeota lacks nomenclatural validation. Here, we describe a novel high-quality metagenome-assembled genome (MAG), AB3033_2TS, proposed to serve as the nomenclatural type for the species Asgardarchaeum abyssiTS according to the rules of the SeqCode. Based on protein content and compositional features, we infer that A. abyssi AB3033_2TS is an acetogenic chemoheterotroph, possibly a facultative lithoautotroph, and is adapted to a thermophilic lifestyle. Utilizing genomes from Asgard archaea, TACK, and Euryarchaea, we perform phylogenomic reconstructions using the GTDB archaeal marker genes, the current reference set for taxonomic classification. Calibrating relative evolutionary divergence (RED) values for Asgardarchaeota using established Thermoproteota lineages in the GTDB r207 reference tree, we establish a robust classification and propose Asgardarchaeum as the type genus for the family Asgardarchaeaceae (fam. nov)., the order Asgardarchaeales (ord. nov.), the class Asgardarchaeia (class. nov.), and the phylum Asgardarchaeota (phyl. nov.). This effort aims to preserve taxonomic congruence in the scientific literature.

Microbial Diversity and Open Questions about the Deep Tree of Life.

In this perspective, we explore the transformative impact and inherent limitations of metagenomics and single-cell genomics on our understanding of microbial diversity and their integration into the Tree of Life. We delve into the key challenges associated with incorporating new microbial lineages into the Tree of Life through advanced phylogenomic approaches. Additionally, we shed light on enduring debates surrounding various aspects of the microbial Tree of Life, focusing on recent advances in some of its deepest nodes, such as the roots of bacteria, archaea, and eukaryotes. We also bring forth current limitations in genome recovery and phylogenomic methodology, as well as new avenues of research to uncover additional key microbial lineages and resolve the shape of the Tree of Life.

Inference and reconstruction of the heimdallarchaeial ancestry of eukaryotes.

In the ongoing debates about eukaryogenesis-the series of evolutionary events leading to the emergence of the eukaryotic cell from prokaryotic ancestors-members of the Asgard archaea play a key part as the closest archaeal relatives of eukaryotes1. However, the nature and phylogenetic identity of the last common ancestor of Asgard archaea and eukaryotes remain unresolved2-4. Here we analyse distinct phylogenetic marker datasets of an expanded genomic sampling of Asgard archaea and evaluate competing evolutionary scenarios using state-of-the-art phylogenomic approaches. We find that eukaryotes are placed, with high confidence, as a well-nested clade within Asgard archaea and as a sister lineage to Hodarchaeales, a newly proposed order within Heimdallarchaeia. Using sophisticated gene tree and species tree reconciliation approaches, we show that analogous to the evolution of eukaryotic genomes, genome evolution in Asgard archaea involved significantly more gene duplication and fewer gene loss events compared with other archaea. Finally, we infer that the last common ancestor of Asgard archaea was probably a thermophilic chemolithotroph and that the lineage from which eukaryotes evolved adapted to mesophilic conditions and acquired the genetic potential to support a heterotrophic lifestyle. Our work provides key insights into the prokaryote-to-eukaryote transition and a platform for better understanding the emergence of cellular complexity in eukaryotic cells.

Genome Evolution of a Symbiont Population for Pathogen Defense in Honeybees.

The honeybee gut microbiome is thought to be important for bee health, but the role of the individual members is poorly understood. Here, we present closed genomes and associated mobilomes of 102 Apilactobacillus kunkeei isolates obtained from the honey crop (foregut) of honeybees sampled from beehives in Helsingborg in the south of Sweden and from the islands Gotland and Åland in the Baltic Sea. Each beehive contained a unique composition of isolates and repeated sampling of similar isolates from two beehives in Helsingborg suggests that the bacterial community is stably maintained across bee generations during the summer months. The sampled bacterial population contained an open pan-genome structure with a high genomic density of transposons. A subset of strains affiliated with phylogroup A inhibited growth of the bee pathogen Melissococcus plutonius, all of which contained a 19.5 kb plasmid for the synthesis of the antimicrobial compound kunkecin A, while a subset of phylogroups B and C strains contained a 32.9 kb plasmid for the synthesis of a putative polyketide antibiotic. This study suggests that the mobile gene pool of A. kunkeei plays a key role in pathogen defense in honeybees, providing new insights into the evolutionary dynamics of defensive symbiont populations.

Asgard archaea shed light on the evolutionary origins of the eukaryotic ubiquitin-ESCRT machinery.

The ESCRT machinery, comprising of multiple proteins and subcomplexes, is crucial for membrane remodelling in eukaryotic cells, in processes that include ubiquitin-mediated multivesicular body formation, membrane repair, cytokinetic abscission, and virus exit from host cells. This ESCRT system appears to have simpler, ancient origins, since many archaeal species possess homologues of ESCRT-III and Vps4, the components that execute the final membrane scission reaction, where they have been shown to play roles in cytokinesis, extracellular vesicle formation and viral egress. Remarkably, metagenome assemblies of Asgard archaea, the closest known living relatives of eukaryotes, were recently shown to encode homologues of the entire cascade involved in ubiquitin-mediated membrane remodelling, including ubiquitin itself, components of the ESCRT-I and ESCRT-II subcomplexes, and ESCRT-III and Vps4. Here, we explore the phylogeny, structure, and biochemistry of Asgard homologues of the ESCRT machinery and the associated ubiquitylation system. We provide evidence for the ESCRT-I and ESCRT-II subcomplexes being involved in ubiquitin-directed recruitment of ESCRT-III, as it is in eukaryotes. Taken together, our analyses suggest a pre-eukaryotic origin for the ubiquitin-coupled ESCRT system and a likely path of ESCRT evolution via a series of gene duplication and diversification events.

A closed Candidatus Odinarchaeum chromosome exposes Asgard archaeal viruses.

Asgard archaea have recently been identified as the closest archaeal relatives of eukaryotes. Their ecology, and particularly their virome, remain enigmatic. We reassembled and closed the chromosome of Candidatus Odinarchaeum yellowstonii LCB_4, through long-range PCR, revealing CRISPR spacers targeting viral contigs. We found related viruses in the genomes of diverse prokaryotes from geothermal environments, including other Asgard archaea. These viruses open research avenues into the ecology and evolution of Asgard archaea.

Chlamydial contribution to anaerobic metabolism during eukaryotic evolution.

The origin of eukaryotes is a major open question in evolutionary biology. Multiple hypotheses posit that eukaryotes likely evolved from a syntrophic relationship between an archaeon and an alphaproteobacterium based on H2 exchange. However, there are no strong indications that modern eukaryotic H2 metabolism originated from archaea or alphaproteobacteria. Here, we present evidence for the origin of H2 metabolism genes in eukaryotes from an ancestor of the Anoxychlamydiales-a group of anaerobic chlamydiae, newly described here, from marine sediments. Among Chlamydiae, these bacteria uniquely encode genes for H2 metabolism and other anaerobiosis-associated pathways. Phylogenetic analyses of several components of H2 metabolism reveal that Anoxychlamydiales homologs are the closest relatives to eukaryotic sequences. We propose that an ancestor of the Anoxychlamydiales contributed these key genes during the evolution of eukaryotes, supporting a mosaic evolutionary origin of eukaryotic metabolism.

Marine Sediments Illuminate Chlamydiae Diversity and Evolution.

The bacterial phylum Chlamydiae is so far composed of obligate symbionts of eukaryotic hosts. Well known for Chlamydiaceae, pathogens of humans and other animals, Chlamydiae also include so-called environmental lineages that primarily infect microbial eukaryotes. Environmental surveys indicate that Chlamydiae are found in a wider range of environments than anticipated previously. However, the vast majority of this chlamydial diversity has been underexplored, biasing our current understanding of their biology, ecological importance, and evolution. Here, we report that previously undetected and active chlamydial lineages dominate microbial communities in deep anoxic marine sediments taken from the Arctic Mid-Ocean Ridge. Reaching relative abundances of up to 43% of the bacterial community, and a maximum diversity of 163 different species-level taxonomic units, these Chlamydiae represent important community members. Using genome-resolved metagenomics, we reconstructed 24 draft chlamydial genomes, expanding by over a third the known genomic diversity in this phylum. Phylogenomic analyses revealed several novel clades across the phylum, including a previously unknown sister lineage of the Chlamydiaceae, providing new insights into the origin of pathogenicity in this family. We were unable to identify putative eukaryotic hosts for these marine sediment chlamydiae, despite identifying genomic features that may be indicative of host-association. The high abundance and genomic diversity of Chlamydiae in these anoxic marine sediments indicate that some members could play an important, and thus far overlooked, ecological role in such environments and may indicate alternate lifestyle strategies.

Origin and Evolution of the Bartonella Gene Transfer Agent.

Gene transfer agents (GTAs) are domesticated bacteriophages that have evolved into molecular machines for the transfer of bacterial DNA. Despite their widespread nature and their biological implications, the mechanisms and selective forces that drive the emergence of GTAs are still poorly understood. Two GTAs have been identified in the Alphaproteobacteria: the RcGTA, which is widely distributed in a broad range of species; and the BaGTA, which has a restricted host range that includes vector-borne intracellular bacteria of the genus Bartonella. The RcGTA packages chromosomal DNA randomly, whereas the BaGTA particles contain a relatively higher fraction of genes for host interaction factors that are amplified from a nearby phage-derived origin of replication. In this study, we compare the BaGTA genes with homologous bacteriophage genes identified in the genomes of Bartonella species and close relatives. Unlike the BaGTA, the prophage genes are neither present in all species, nor inserted into homologous genomic sites. Phylogenetic inferences and substitution frequency analyses confirm codivergence of the BaGTA with the host genome, as opposed to multiple integration and recombination events in the prophages. Furthermore, the organization of segments flanking the BaGTA differs from that of the prophages by a few rearrangement events, which have abolished the normal coordination between phage genome replication and phage gene expression. Based on the results of our comparative analysis, we propose a model for how a prophage may be transformed into a GTA that transfers amplified bacterial DNA segments.

Evolution and postglacial colonization of Seewis hantavirus with Sorex araneus in Finland.

Hantaviruses have co-existed with their hosts for millions of years. Seewis virus (SWSV), a soricomorph-borne hantavirus, is widespread in Eurasia, ranging from Central Siberia to Western Europe. To gain insight into the phylogeography and evolutionary history of SWSV in Finland, lung tissue samples of 225 common shrews (Sorex araneus) trapped from different parts of Finland were screened for the presence of SWSV RNA. Forty-two of the samples were positive. Partial small (S), medium (M) and large (L) segments of the virus were sequenced, and analyzed together with all SWSV sequences available in Genbank. The phylogenetic analysis of the partial S-segment sequences suggested that all Finnish SWSV strains shared their most recent common ancestor with the Eastern European strains, while the L-segment suggested multiple introductions. The difference between the L- and S-segment phylogenies implied that reassortment events play a role in the evolution of SWSV. Of the Finnish strains, variants from Eastern Finland occupied the root position in the phylogeny, and had the highest genetic diversity, supporting the hypothesis that SWSV reached Finland first form the east. During the spread in Finland, the virus has formed three separate lineages, identified here by correlation analysis of genetic versus geographic distance combined with median-joining network analysis. These results support the hypothesis that Finnish SWSV recolonized Finland with its host, the common shrew, from east after the last ice age 12,000-8000years ago, and then subsequently spread along emerging land bridges towards west or north with the migration and population expansion of its host.

Determinism and Contingency Shape Metabolic Complementation in an Endosymbiotic Consortium.

Bacterial endosymbionts and their insect hosts establish an intimate metabolic relationship. Bacteria offer a variety of essential nutrients to their hosts, whereas insect cells provide the necessary sources of matter and energy to their tiny metabolic allies. These nutritional complementations sustain themselves on a diversity of metabolite exchanges between the cell host and the reduced yet highly specialized bacterial metabolism-which, for instance, overproduces a small set of essential amino acids and vitamins. A well-known case of metabolic complementation is provided by the cedar aphid Cinara cedri that harbors two co-primary endosymbionts, Buchnera aphidicola BCc and Ca. Serratia symbiotica SCc, and in which some metabolic pathways are partitioned between different partners. Here we present a genome-scale metabolic network (GEM) for the bacterial consortium from the cedar aphid iBSCc. The analysis of this GEM allows us the confirmation of cases of metabolic complementation previously described by genome analysis (i.e., tryptophan and biotin biosynthesis) and the redefinition of an event of metabolic pathway sharing between the two endosymbionts, namely the biosynthesis of tetrahydrofolate. In silico knock-out experiments with iBSCc showed that the consortium metabolism is a highly integrated yet fragile network. We also have explored the evolutionary pathways leading to the emergence of metabolic complementation between reduced metabolisms starting from individual, complete networks. Our results suggest that, during the establishment of metabolic complementation in endosymbionts, adaptive evolution is significant in the case of tryptophan biosynthesis, whereas vitamin production pathways seem to adopt suboptimal solutions.

Switches in Genomic GC Content Drive Shifts of Optimal Codons under Sustained Selection on Synonymous Sites.

The major codon preference model suggests that codons read by tRNAs in high concentrations are preferentially utilized in highly expressed genes. However, the identity of the optimal codons differs between species although the forces driving such changes are poorly understood. We suggest that these questions can be tackled by placing codon usage studies in a phylogenetic framework and that bacterial genomes with extreme nucleotide composition biases provide informative model systems. Switches in the background substitution biases from GC to AT have occurred in Gardnerella vaginalis (GC = 32%), and from AT to GC in Lactobacillus delbrueckii (GC = 62%) and Lactobacillus fermentum (GC = 63%). We show that despite the large effects on codon usage patterns by these switches, all three species evolve under selection on synonymous sites. In G. vaginalis, the dramatic codon frequency changes coincide with shifts of optimal codons. In contrast, the optimal codons have not shifted in the two Lactobacillus genomes despite an increased fraction of GC-ending codons. We suggest that all three species are in different phases of an on-going shift of optimal codons, and attribute the difference to a stronger background substitution bias and/or longer time since the switch in G. vaginalis. We show that comparative and correlative methods for optimal codon identification yield conflicting results for genomes in flux and discuss possible reasons for the mispredictions. We conclude that switches in the direction of the background substitution biases can drive major shifts in codon preference patterns even under sustained selection on synonymous codon sites.

The genome of Rhizobiales bacteria in predatory ants reveals urease gene functions but no genes for nitrogen fixation.

Gut-associated microbiota of ants include Rhizobiales bacteria with affiliation to the genus Bartonella. These bacteria may enable the ants to fix atmospheric nitrogen, but no genomes have been sequenced yet to test the hypothesis. Sequence reads from a member of the Rhizobiales were identified in the data collected in a genome project of the ant Harpegnathos saltator. We present an analysis of the closed 1.86 Mb genome of the ant-associated bacterium, for which we suggest the species name Candidatus Tokpelaia hoelldoblerii. A phylogenetic analysis reveals a relationship to Bartonella and Brucella, which infect mammals. Novel gene acquisitions include a gene for a putative extracellular protein of more than 6,000 amino acids secreted by the type I secretion system, which may be involved in attachment to the gut epithelium. No genes for nitrogen fixation could be identified, but genes for a multi-subunit urease protein complex are present in the genome. The urease genes are also present in Brucella, which has a fecal-oral transmission pathway, but not in Bartonella, which use blood-borne transmission pathways. We hypothesize that the gain and loss of the urease function is related to transmission strategies and lifestyle changes in the host-associated members of the Rhizobiales.

Functionally Structured Genomes in Lactobacillus kunkeei Colonizing the Honey Crop and Food Products of Honeybees and Stingless Bees.

Lactobacillus kunkeei is the most abundant bacterial species in the honey crop and food products of honeybees. The 16 S rRNA genes of strains isolated from different bee species are nearly identical in sequence and therefore inadequate as markers for studies of coevolutionary patterns. Here, we have compared the 1.5 Mb genomes of ten L. kunkeei strains isolated from all recognized Apis species and another two strains from Meliponini species. A gene flux analysis, including previously sequenced Lactobacillus species as outgroups, indicated the influence of reductive evolution. The genome architecture is unique in that vertically inherited core genes are located near the terminus of replication, whereas genes for secreted proteins and putative host-adaptive traits are located near the origin of replication. We suggest that these features have resulted from a genome-wide loss of genes, with integrations of novel genes mostly occurring in regions flanking the origin of replication. The phylogenetic analyses showed that the bacterial topology was incongruent with the host topology, and that strains of the same microcluster have recombined frequently across the host species barriers, arguing against codiversification. Multiple genotypes were recovered in the individual hosts and transfers of mobile elements could be demonstrated for strains isolated from the same host species. Unlike other bacteria with small genomes, short generation times and multiple rRNA operons suggest that L. kunkeei evolves under selection for rapid growth in its natural growth habitat. The results provide an extended framework for reductive genome evolution and functional genome organization in bacteria.

Extensive intra-phylotype diversity in lactobacilli and bifidobacteria from the honeybee gut.

BACKGROUND: In the honeybee Apis mellifera, the bacterial gut community is consistently colonized by eight distinct phylotypes of bacteria. Managed bee colonies are of considerable economic interest and it is therefore important to elucidate the diversity and role of this microbiota in the honeybee. In this study, we have sequenced the genomes of eleven strains of lactobacilli and bifidobacteria isolated from the honey crop of the honeybee A. mellifera. RESULTS: Single gene phylogenies confirmed that the isolated strains represent the diversity of lactobacilli and bifidobacteria in the gut, as previously identified by 16S rRNA gene sequencing. Core genome phylogenies of the lactobacilli and bifidobacteria further indicated extensive divergence between strains classified as the same phylotype. Phylotype-specific protein families included unique surface proteins. Within phylotypes, we found a remarkably high level of gene content diversity. Carbohydrate metabolism and transport functions contributed up to 45% of the accessory genes, with some genomes having a higher content of genes encoding phosphotransferase systems for the uptake of carbohydrates than any previously sequenced genome. These genes were often located in highly variable genomic segments that also contained genes for enzymes involved in the degradation and modification of sugar residues. Strain-specific gene clusters for the biosynthesis of exopolysaccharides were identified in two phylotypes. The dynamics of these segments contrasted with low recombination frequencies and conserved gene order structures for the core genes. Hits for CRISPR spacers were almost exclusively found within phylotypes, suggesting that the phylotypes are associated with distinct phage populations. CONCLUSIONS: The honeybee gut microbiota has been described as consisting of a modest number of phylotypes; however, the genomes sequenced in the current study demonstrated a very high level of gene content diversity within all three described phylotypes of lactobacilli and bifidobacteria, particularly in terms of metabolic functions and surface structures, where many features were strain-specific. Together, these results indicate niche differentiation within phylotypes, suggesting that the honeybee gut microbiota is more complex than previously thought.

The Gypsy Database (GyDB) of mobile genetic elements: release 2.0.

This article introduces the second release of the Gypsy Database of Mobile Genetic Elements (GyDB 2.0): a research project devoted to the evolutionary dynamics of viruses and transposable elements based on their phylogenetic classification (per lineage and protein domain). The Gypsy Database (GyDB) is a long-term project that is continuously progressing, and that owing to the high molecular diversity of mobile elements requires to be completed in several stages. GyDB 2.0 has been powered with a wiki to allow other researchers participate in the project. The current database stage and scope are long terminal repeats (LTR) retroelements and relatives. GyDB 2.0 is an update based on the analysis of Ty3/Gypsy, Retroviridae, Ty1/Copia and Bel/Pao LTR retroelements and the Caulimoviridae pararetroviruses of plants. Among other features, in terms of the aforementioned topics, this update adds: (i) a variety of descriptions and reviews distributed in multiple web pages; (ii) protein-based phylogenies, where phylogenetic levels are assigned to distinct classified elements; (iii) a collection of multiple alignments, lineage-specific hidden Markov models and consensus sequences, called GyDB collection; (iv) updated RefSeq databases and BLAST and HMM servers to facilitate sequence characterization of new LTR retroelement and caulimovirus queries; and (v) a bibliographic server. GyDB 2.0 is available at http://gydb.org.

Immunity and other defenses in pea aphids, Acyrthosiphon pisum.

BACKGROUND: Recent genomic analyses of arthropod defense mechanisms suggest conservation of key elements underlying responses to pathogens, parasites and stresses. At the center of pathogen-induced immune responses are signaling pathways triggered by the recognition of fungal, bacterial and viral signatures. These pathways result in the production of response molecules, such as antimicrobial peptides and lysozymes, which degrade or destroy invaders. Using the recently sequenced genome of the pea aphid (Acyrthosiphon pisum), we conducted the first extensive annotation of the immune and stress gene repertoire of a hemipterous insect, which is phylogenetically distantly related to previously characterized insects models. RESULTS: Strikingly, pea aphids appear to be missing genes present in insect genomes characterized to date and thought critical for recognition, signaling and killing of microbes. In line with results of gene annotation, experimental analyses designed to characterize immune response through the isolation of RNA transcripts and proteins from immune-challenged pea aphids uncovered few immune-related products. Gene expression studies, however, indicated some expression of immune and stress-related genes. CONCLUSIONS: The absence of genes suspected to be essential for the insect immune response suggests that the traditional view of insect immunity may not be as broadly applicable as once thought. The limitations of the aphid immune system may be representative of a broad range of insects, or may be aphid specific. We suggest that several aspects of the aphid life style, such as their association with microbial symbionts, could facilitate survival without strong immune protection.