Recoding of stop codons expands the metabolic potential of two novel Asgardarchaeota lineages

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.

Analysis of nearly 3000 archaeal genomes from terrestrial geothermal springs sheds light on interconnected biogeochemical processes

Terrestrial geothermal springs are physicochemically diverse and host abundant populations of Archaea. However, the diversity, functionality, and geological influences of these Archaea are not well understood. Here we explore the genomic diversity of Archaea in 152 metagenomes from 48 geothermal springs in Tengchong, China, collected from 2016 to 2021. Our dataset is comprised of 2949 archaeal metagenome-assembled genomes spanning 12 phyla and 392 newly identified species, which increases the known species diversity of Archaea by ~48.6%. The structures and potential functions of the archaeal communities are strongly influenced by temperature and pH, with high-temperature acidic and alkaline springs favoring archaeal abundance over Bacteria. Genome-resolved metagenomics and metatranscriptomics provide insights into the potential ecological niches of these Archaea and their potential roles in carbon, sulfur, nitrogen, and hydrogen metabolism. Furthermore, our findings illustrate the interplay of competition and cooperation among Archaea in biogeochemical cycles, possibly arising from overlapping functional niches and metabolic handoffs. Taken together, our study expands the genomic diversity of Archaea inhabiting geothermal springs and provides a foundation for more incisive study of biogeochemical processes mediated by Archaea in geothermal ecosystems.

The Selenoproteome as a Dynamic Response Mechanism to Oxidative Stress in Hydrogenotrophic Methanogenic Communities

Methanogenesis is a critical process in the carbon cycle that is applied industrially in anaerobic digestion and biogas production. While naturally occurring in diverse environments, methanogenesis requires anaerobic and reduced conditions, although varying degrees of oxygen tolerance have been described. Microaeration is suggested as the next step to increase methane production and improve hydrolysis in digestion processes; therefore, a deeper understanding of the methanogenic response to oxygen stress is needed. To explore the drivers of oxygen tolerance in methanogenesis, two parallel enrichments were performed under the addition of H2/CO2 in an environment without reducing agents and in a redox-buffered environment by adding redox mediator 9,10-anthraquinone-2,7-disulfonate disodium. The cellular response to oxidative conditions is mapped using proteomic analysis. The resulting community showed remarkable tolerance to high-redox environments and was unperturbed in its methane production. Next to the expression of pathways to mitigate reactive oxygen species, the higher redox potential environment showed an increased presence of selenocysteine and selenium-associated pathways. By including sulfur-to-selenium mass shifts in a proteomic database search, we provide the first evidence of the dynamic and large-scale incorporation of selenocysteine as a response to oxidative stress in hydrogenotrophic methanogenesis and the presence of a dynamic selenoproteome.

Asgard archaeal selenoproteome reveals a roadmap for the archaea-to-eukaryote transition of selenocysteine incorporation machinery

Selenocysteine (Sec) is encoded by the UGA codon that normally functions as a stop signal and is specifically incorporated into selenoproteins via a unique recoding mechanism. The translational recoding of UGA as Sec is directed by an unusual RNA structure, the SECIS element. Although archaea and eukaryotes adopt similar Sec encoding machinery, the SECIS elements have no similarities to each other with regard to sequence and structure. We analyzed >400 Asgard archaeal genomes to examine the occurrence of both Sec encoding system and selenoproteins in this archaeal superphylum, the closest prokaryotic relatives of eukaryotes. A comprehensive map of Sec utilization trait has been generated, providing the most detailed understanding of the use of this nonstandard amino acid in Asgard archaea so far. By characterizing the selenoproteomes of all organisms, several selenoprotein-rich phyla and species were identified. Most Asgard archaeal selenoprotein genes possess eukaryotic SECIS-like structures with varying degrees of diversity. Moreover, euryarchaeal SECIS elements might originate from Asgard archaeal SECIS elements via lateral gene transfer, indicating a complex and dynamic scenario of the evolution of SECIS element within archaea. Finally, a roadmap for the transition of eukaryotic SECIS elements from archaea was proposed, and selenophosphate synthetase may serve as a potential intermediate for the generation of ancestral eukaryotic SECIS element. Our results offer new insights into a deeper understanding of the evolution of Sec insertion machinery.

Xeno Amino Acids: A Look into Biochemistry as We Do Not Know It

Would another origin of life resemble Earth's biochemical use of amino acids? Here, we review current knowledge at three levels: (1) Could other classes of chemical structure serve as building blocks for biopolymer structure and catalysis? Amino acids now seem both readily available to, and a plausible chemical attractor for, life as we do not know it. Amino acids thus remain important and tractable targets for astrobiological research. (2) If amino acids are used, would we expect the same L-alpha-structural subclass used by life? Despite numerous ideas, it is not clear why life favors L-enantiomers. It seems clearer, however, why life on Earth uses the shortest possible (alpha-) amino acid backbone, and why each carries only one side chain. However, assertions that other backbones are physicochemically impossible have relaxed into arguments that they are disadvantageous. (3) Would we expect a similar set of side chains to those within the genetic code? Many plausible alternatives exist. Furthermore, evidence exists for both evolutionary advantage and physicochemical constraint as explanatory factors for those encoded by life. Overall, as focus shifts from amino acids as a chemical class to specific side chains used by post-LUCA biology, the probable role of physicochemical constraint diminishes relative to that of biological evolution. Exciting opportunities now present themselves for laboratory work and computing to explore how changing the amino acid alphabet alters the universe of protein folds. Near-term milestones include: (a) expanding evidence about amino acids as attractors within chemical evolution; (b) extending characterization of other backbones relative to biological proteins; and (c) merging computing and laboratory explorations of structures and functions unlocked by xeno peptides.

Ribosomal Protein Cluster Organization in Asgard Archaea

It has been proposed that the superphylum of Asgard Archaea may represent a historical link between the Archaea and Eukarya. Following the discovery of the Archaea, it was soon appreciated that archaeal ribosomes were more similar to those of Eukarya rather than Bacteria. Coupled with other eukaryotic-like features, it has been suggested that the Asgard Archaea may be directly linked to eukaryotes. However, the genomes of Bacteria and non-Asgard Archaea generally organize ribosome-related genes into clusters that likely function as operons. In contrast, eukaryotes typically do not employ an operon strategy. To gain further insight into conservation of the r-protein genes, the genome order of conserved ribosomal protein (r-protein) coding genes was identified in 17 Asgard genomes (thirteen complete genomes and four genomes with less than 20 contigs) and compared with those found previously in non-Asgard archaeal and bacterial genomes. A universal core of two clusters of 14 and 4 cooccurring r-proteins, respectively, was identified in both the Asgard and non-Asgard Archaea. The equivalent genes in the E. coli version of the cluster are found in the S10 and spc operons. The large cluster of 14 r-protein genes (uS19-uL22-uS3-uL29-uS17 from the S10 operon and uL14-uL24-uL5-uS14-uS8-uL6-uL18-uS5-uL30-uL15 from the spc operon) occurs as a complete set in the genomes of thirteen Asgard genomes (five Lokiarchaeotes, three Heimdallarchaeotes, one Odinarchaeote, and four Thorarchaeotes). Four less conserved clusters with partial bacterial equivalents were found in the Asgard. These were the L30e (str operon in Bacteria) cluster, the L18e (alpha operon in Bacteria) cluster, the S24e-S27ae-rpoE1 cluster, and the L31e, L12..L1 cluster. Finally, a new cluster referred to as L7ae was identified. In many cases, r-protein gene clusters/operons are less conserved in their organization in the Asgard group than in other Archaea. If this is generally true for nonribosomal gene clusters, the results may have implications for the history of genome organization. In particular, there may have been an early transition to or from the operon approach to genome organization. Other nonribosomal cellular features may support different relationships. For this reason, it may be important to consider ribosome features separately.

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.

Update of the Pyrrolysyl-tRNA Synthetase/tRNAPyl Pair and Derivatives for Genetic Code Expansion

The cotranslational incorporation of pyrrolysine (Pyl), the 22nd proteinogenic amino acid, into proteins in response to the UAG stop codon represents an outstanding example of natural genetic code expansion. Genetic encoding of Pyl is conducted by the pyrrolysyl-tRNA synthetase (PylRS) and its cognate tRNA, tRNA^Pyl. Owing to the high tolerance of PylRS toward diverse amino acid substrates and great orthogonality in various model organisms, the PylRS/tRNA^Pyl-derived pairs are ideal for genetic code expansion to insert noncanonical amino acids (ncAAs) into proteins of interest. Since the discovery of cellular components involved in the biosynthesis and genetic encoding of Pyl, synthetic biologists have been enthusiastic about engineering PylRS/tRNA^Pyl-derived pairs to rewrite the genetic code of living cells. Recently, considerable progress has been made in understanding the molecular phylogeny, biochemical properties, and structural features of the PylRS/tRNA^Pyl pair, guiding its further engineering and optimization. In this review, we cover the basic and updated knowledge of the PylRS/tRNA^Pyl pair's unique characteristics that make it an outstanding tool for reprogramming the genetic code. In addition, we summarize the recent efforts to create efficient and (mutually) orthogonal PylRS/tRNA^Pyl-derived pairs for incorporation of diverse ncAAs by genome mining, rational design, and advanced directed evolution methods.

Integrating Phylogenetics With Intron Positions Illuminates the Origin of the Complex Spliceosome

Eukaryotic genes are characterized by the presence of introns that are removed from pre-mRNA by a spliceosome. This ribonucleoprotein complex is comprised of multiple RNA molecules and over a hundred proteins, which makes it one of the most complex molecular machines that originated during the prokaryote-to-eukaryote transition. Previous works have established that these introns and the spliceosomal core originated from self-splicing introns in prokaryotes. Yet, how the spliceosomal core expanded by recruiting many additional proteins remains largely elusive. In this study, we use phylogenetic analyses to infer the evolutionary history of 145 proteins that we could trace back to the spliceosome in the last eukaryotic common ancestor. We found that an overabundance of proteins derived from ribosome-related processes was added to the prokaryote-derived core. Extensive duplications of these proteins substantially increased the complexity of the emerging spliceosome. By comparing the intron positions between spliceosomal paralogs, we infer that most spliceosomal complexity postdates the spread of introns through the proto-eukaryotic genome. The reconstruction of early spliceosomal evolution provides insight into the driving forces behind the emergence of complexes with many proteins during eukaryogenesis.

Actin cytoskeleton and complex cell architecture in an Asgard archaeon

Asgard archaea are considered to be the closest known relatives of eukaryotes. Their genomes contain hundreds of eukaryotic signature proteins (ESPs), which inspired hypotheses on the evolution of the eukaryotic cell^1-3. A role of ESPs in the formation of an elaborate cytoskeleton and complex cellular structures has been postulated^4-6, but never visualized. Here we describe a highly enriched culture of 'Candidatus Lokiarchaeum ossiferum', a member of the Asgard phylum, which thrives anaerobically at 20 °C on organic carbon sources. It divides every 7-14 days, reaches cell densities of up to 5 × 10⁷ cells per ml and has a significantly larger genome compared with the single previously cultivated Asgard strain⁷. ESPs represent 5% of its protein-coding genes, including four actin homologues. We imaged the enrichment culture using cryo-electron tomography, identifying 'Ca. L. ossiferum' cells on the basis of characteristic expansion segments of their ribosomes. Cells exhibited coccoid cell bodies and a network of branched protrusions with frequent constrictions. The cell envelope consists of a single membrane and complex surface structures. A long-range cytoskeleton extends throughout the cell bodies, protrusions and constrictions. The twisted double-stranded architecture of the filaments is consistent with F-actin. Immunostaining indicates that the filaments comprise Lokiactin-one of the most highly conserved ESPs in Asgard archaea. We propose that a complex actin-based cytoskeleton predated the emergence of the first eukaryotes and was a crucial feature in the evolution of the Asgard phylum by scaffolding elaborate cellular structures.

The expanding Asgard archaea invoke novel insights into Tree of Life and eukaryogenesis

The division of organisms on the Tree of Life into either a three-domain (3D) tree or a two-domain (2D) tree has been disputed for a long time. Ever since the discovery of Archaea by Carl Woese in 1977 using 16S ribosomal RNA sequence as the evolutionary marker, there has been a great advance in our knowledge of not only the growing diversity of Archaea but also the evolutionary relationships between different lineages of living organisms. Here, we present this perspective to summarize the progress of archaeal diversity and changing notion of the Tree of Life. Meanwhile, we provide the latest progress in genomics/physiology-based discovery of Asgard archaeal lineages as the closest relative of Eukaryotes. Furthermore, we propose three major directions for future research on exploring the "next one" closest Eukaryote relative, deciphering the function of archaeal eukaryotic signature proteins and eukaryogenesis from both genomic and physiological aspects, and understanding the roles of horizontal gene transfer, viruses, and mobile elements in eukaryogenesis.

The selenophosphate synthetase family: A review

Selenophosphate synthetases use selenium and ATP to synthesize selenophosphate. This is required for biological utilization of selenium, most notably for the synthesis of the non-canonical amino acid selenocysteine (Sec). Therefore, selenophosphate synthetases underlie all functions of selenoproteins, which include redox homeostasis, protein quality control, hormone regulation, metabolism, and many others. This protein family comprises two groups, SelD/SPS2 and SPS1. The SelD/SPS2 group represent true selenophosphate synthetases, enzymes central to selenium metabolism which are present in all Sec-utilizing organisms across the tree of life. Notably, many SelD/SPS2 proteins contain Sec as catalytic residue in their N-terminal flexible selenium-binding loop, while others replace it with cysteine (Cys). The SPS1 group comprises proteins originated through gene duplications of SelD/SPS2 in metazoa in which the Sec/Cys-dependent catalysis was disrupted. SPS1 proteins do not synthesize selenophosphate and are not required for Sec synthesis. They have essential regulatory functions related to redox homeostasis and pyridoxal phosphate, which affect signaling pathways for growth and differentiation. In this review, we summarize the knowledge about the selenophosphate synthetase family acquired through decades of research, encompassing their structure, mechanism, function, and evolution.

Ancestral Archaea Expanded the Genetic Code with Pyrrolysine

The pyrrolysyl-tRNA synthetase (PylRS) facilitates the co-translational installation of the 22^nd amino acid pyrrolysine. Owing to its tolerance for diverse amino acid substrates, and its orthogonality in multiple organisms, PylRS has emerged as a major route to install noncanonical amino acids into proteins in living cells. Recently, a novel class of PylRS enzymes was identified in a subset of methanogenic archaea. Enzymes within this class (ΔPylSn) lack the N-terminal tRNA-binding domain that is widely conserved amongst PylRS enzymes, yet remain highly active and orthogonal in bacteria and eukaryotes. In this study, we use biochemical and in vivo UAG-readthrough assays to characterize the aminoacylation efficiency and substrate spectrum of a ΔPylSn class PylRS from the archaeon Ca. Methanomethylophilus alvus. We show that, compared to the full-length enzyme from Methanosarcina mazei, the Ca. M. alvus PylRS displays reduced aminoacylation efficiency, but an expanded amino acid substrate spectrum. To gain insight into the evolution of ΔPylSn enzymes, we performed molecular phylogeny using 156 PylRS and 105 tRNA^Pyl sequences from diverse anaerobic archaea and bacteria. This analysis suggests that the PylRS•tRNA^Pyl pair diverged before the evolution of the three domains of life, placing an early limit on the evolution of the Pyl-decoding trait. Furthermore, our results document the co-evolutionary history of PylRS and tRNA^Pyl and reveal the emergence of tRNA^Pyl sequences with unique A73 and U73 discriminator bases. The orthogonality of these tRNA^Pyl species with the more common G73-containing tRNA^Pyl will enable future efforts to engineer PylRS systems for further genetic code expansion.

Unconventional genetic code systems in archaea

Archaea constitute the third domain of life, distinct from bacteria and eukaryotes given their ability to tolerate extreme environments. To survive these harsh conditions, certain archaeal lineages possess unique genetic code systems to encode either selenocysteine or pyrrolysine, rare amino acids not found in all organisms. Furthermore, archaea utilize alternate tRNA-dependent pathways to biosynthesize and incorporate members of the 20 canonical amino acids. Recent discoveries of new archaeal species have revealed the co-occurrence of these genetic code systems within a single lineage. This review discusses the diverse genetic code systems of archaea, while detailing the associated biochemical elements and molecular mechanisms.

Three families of Asgard archaeal viruses identified in metagenome-assembled genomes

Asgardarchaeota harbour many eukaryotic signature proteins and are widely considered to represent the closest archaeal relatives of eukaryotes. Whether similarities between Asgard archaea and eukaryotes extend to their viromes remains unknown. Here we present 20 metagenome-assembled genomes of Asgardarchaeota from deep-sea sediments of the basin off the Shimokita Peninsula, Japan. By combining a CRISPR spacer search of metagenomic sequences with phylogenomic analysis, we identify three family-level groups of viruses associated with Asgard archaea. The first group, verdandiviruses, includes tailed viruses of the class Caudoviricetes (realm Duplodnaviria); the second, skuldviruses, consists of viruses with predicted icosahedral capsids of the realm Varidnaviria; and the third group, wyrdviruses, is related to spindle-shaped viruses previously identified in other archaea. More than 90% of the proteins encoded by these viruses of Asgard archaea show no sequence similarity to proteins encoded by other known viruses. Nevertheless, all three proposed families consist of viruses typical of prokaryotes, providing no indication of specific evolutionary relationships between viruses infecting Asgard archaea and eukaryotes. Verdandiviruses and skuldviruses are likely to be lytic, whereas wyrdviruses potentially establish chronic infection and are released without host cell lysis. All three groups of viruses are predicted to play important roles in controlling Asgard archaea populations in deep-sea ecosystems.

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.

Widespread stop-codon recoding in bacteriophages may regulate translation of lytic genes

Bacteriophages (phages) are obligate parasites that use host bacterial translation machinery to produce viral proteins. However, some phages have alternative genetic codes with reassigned stop codons that are predicted to be incompatible with bacterial translation systems. We analysed 9,422 phage genomes and found that stop-codon recoding has evolved in diverse clades of phages that infect bacteria present in both human and animal gut microbiota. Recoded stop codons are particularly over-represented in phage structural and lysis genes. We propose that recoded stop codons might function to prevent premature production of late-stage proteins. Stop-codon recoding has evolved several times in closely related lineages, which suggests that adaptive recoding can occur over very short evolutionary timescales.

Unique mobile elements and scalable gene flow at the prokaryote-eukaryote boundary revealed by circularized Asgard archaea genomes

Eukaryotic genomes are known to have garnered innovations from both archaeal and bacterial domains but the sequence of events that led to the complex gene repertoire of eukaryotes is largely unresolved. Here, through the enrichment of hydrothermal vent microorganisms, we recovered two circularized genomes of Heimdallarchaeum species that belong to an Asgard archaea clade phylogenetically closest to eukaryotes. These genomes reveal diverse mobile elements, including an integrative viral genome that bidirectionally replicates in a circular form and aloposons, transposons that encode the 5,000 amino acid-sized proteins Otus and Ephialtes. Heimdallaechaeal mobile elements have garnered various genes from bacteria and bacteriophages, likely playing a role in shuffling functions across domains. The number of archaea- and bacteria-related genes follow strikingly different scaling laws in Asgard archaea, exhibiting a genome size-dependent ratio and a functional division resembling the bacteria- and archaea-derived gene repertoire across eukaryotes. Bacterial gene import has thus likely been a continuous process unaltered by eukaryogenesis and scaled up through genome expansion. Our data further highlight the importance of viewing eukaryogenesis in a pan-Asgard context, which led to the proposal of a conceptual framework, that is, the Heimdall nucleation–decentralized innovation–hierarchical import model that accounts for the emergence of eukaryotic complexity.

The recovery of two circularized genomes of the Heimdallarchaeum species from hydrothermal vent enrichment cultures reveals that these Asgard archaea carry diverse mobile genetic elements, such as an integrative viral genome and aloposons. These mobile genetic elements contain several bacteria- and phage-derived genes, modulating the shuffling of information between bacteria and archaea, and potentially influencing eukaryogenesis.

Indirect Routes to Aminoacyl-tRNA: The Diversity of Prokaryotic Cysteine Encoding Systems

Universally present aminoacyl-tRNA synthetases (aaRSs) stringently recognize their cognate tRNAs and acylate them with one of the proteinogenic amino acids. However, some organisms possess aaRSs that deviate from the accurate translation of the genetic code and exhibit relaxed specificity toward their tRNA and/or amino acid substrates. Typically, these aaRSs are part of an indirect pathway in which multiple enzymes participate in the formation of the correct aminoacyl-tRNA product. The indirect cysteine (Cys)-tRNA pathway, originally thought to be restricted to methanogenic archaea, uses the unique O-phosphoseryl-tRNA synthetase (SepRS), which acylates the non-proteinogenic amino acid O-phosphoserine (Sep) onto tRNACys. Together with Sep-tRNA:Cys-tRNA synthase (SepCysS) and the adapter protein SepCysE, SepRS forms a transsulfursome complex responsible for shuttling Sep-tRNACys to SepCysS for conversion of the tRNA-bound Sep to Cys. Here, we report a comprehensive bioinformatic analysis of the diversity of indirect Cys encoding systems. These systems are present in more diverse groups of bacteria and archaea than previously known. Given the occurrence and distribution of some genes consistently flanking SepRS, it is likely that this gene was part of an ancient operon that suffered a gradual loss of its original components. Newly identified bacterial SepRS sequences strengthen the suggestion that this lineage of enzymes may not rely on the m1G37 identity determinant in tRNA. Some bacterial SepRSs possess an N-terminal fusion resembling a threonyl-tRNA synthetase editing domain, which interestingly is frequently observed in the vicinity of archaeal SepCysS genes. We also found several highly degenerate SepRS genes that likely have altered amino acid specificity. Cross-analysis of selenocysteine (Sec)-utilizing traits confirmed the co-occurrence of SepCysE and the Sec-utilizing machinery in archaea, but also identified an unusual O-phosphoseryl-tRNASec kinase fusion with an archaeal Sec elongation factor in some lineages, where it may serve in place of SepCysE to prevent crosstalk between the two minor aminoacylation systems. These results shed new light on the variations in SepRS and SepCysS enzymes that may reflect adaptation to lifestyle and habitat, and provide new information on the evolution of the genetic code.