As more archaeal genomes are sequenced, effective research and analysis tools are needed to integrate the diverse information available for any given locus. The feature-rich UCSC Genome Browser, created originally to annotate the human genome, can be applied to any sequenced organism. We have created a UCSC Archaeal Genome Browser, available at http://archaea.ucsc.edu/, currently with 26 archaeal genomes. It displays G/C content, gene and operon annotation from multiple sources, sequence motifs (promoters and Shine-Dalgarno), microarray data, multi-genome alignments and protein conservation across phylogenetic and habitat categories. We encourage submission of new experimental and bioinformatic analysis from contributors. The purpose of this tool is to aid biological discovery and facilitate greater collaboration within the archaeal research community.

In Archaea, previous studies have revealed the presence of multiple intron-containing tRNAs and split tRNAs. The full unexpurgated analysis of archaeal tRNA genes remains a challenging task in the field of bioinformatics, because of the presence of various types of hidden tRNA genes in archaea. Here, we suggested a computational method that searched for widely separated genes encoding tRNA halves to generate suppressive variants of missing tRNAs.

Archaea-infecting viruses are morphologically and genomically among the most diverse entities. Unfortunately, they are also fairly understudied due to a lack of efficient genetic tools. Here, we present a detailed protocol for the CRISPR/Cas-based genome editing of the virus SIRV2 infecting the genus Sulfolobus, which could easily be adapted to other archaeal viruses. This protocol also includes the procedure for endogenous viral protein purification and identification, allowing for assessing the molecular mechanisms behind virus life cycle and virus-host interactions. For complete details on the use and execution of this protocol, please refer to Mayo-Muñoz et al. (2018) and Bhoobalan-Chitty et al. (2019).

As the null hypothesis of genome evolution, population genetic theory suggests that selection strength controls genome size. Through the process of genetic drift, this theory predicts that compact genomes are maintained by strong purifying selection while complex genomes are enabled by weak purifying selection. It offers a unifying framework that explains why prokaryotic genomes are much smaller than their eukaryotic counterparts. However, recent findings suggest that bigger prokaryotic genomes appear to experience stronger purifying selection, indicating that purifying selection may not dominate prokaryotic genome evolution. Since archaeal genomes were underrepresented in those studies, generalization of the conclusions to both archaeal and bacterial genomes may not be warranted. In this study, we revisited this matter by focusing on archaeal and bacterial genomes separately. We found that bigger bacterial genomes indeed experienced stronger purifying selection, but the opposite was observed in archaeal genomes. This new finding would predict an enrichment of noncoding sequences in large archaeal genomes, which was confirmed by an analysis of coding density. In contrast, coding density remained stable regardless of bacterial genome size. In conclusion, this study suggests that purifying selection may play a more important role in archaeal genome evolution than previously hypothesized, indicating that there could be a major difference between the evolutionary regimes of Archaea and Bacteria. IMPORTANCE The evolution of genome complexity is a fundamental question in biology. A hallmark of eukaryotic genome complexity is that larger genomes tend to have more noncoding sequences, which are believed to be minimal in archaeal and bacterial genomes. However, we found that archaeal genomes also possessed this eukaryotic feature while bacterial genomes did not. This could be predicted from our analysis on genetic drift, which showed a relaxation of purifying selection in larger archaeal genomes, also a eukaryotic feature. In contrast, the opposite was evident in bacterial genomes.

Variation in recombination rates across chromosomes has been shown to be a primary force shaping the architecture of genome divergence. In archaea, little is known about variation in recombination across the chromosome or how it shapes genome evolution. We identified significant variations in polymorphism occurring across the chromosomes of ten closely related sympatric strains of the thermoacidophilic archaeon Sulfolobus islandicus. Statistical analyses show that recombination varies across the genome and interacts with selection to define large genomic regions with reduced polymorphism, particularly in the regions surrounding the three origins of replication. Our findings demonstrate how recombination defines the mosaic of variation in this asexually reproducing microorganism and provide insight into the evolutionary origins of genome architecture in this organism from the Archaeal domain.

No abstract available

Deciphering the structure of gene regulatory networks across the tree of life remains one of the major challenges in postgenomic biology. We present a novel ChIP-seq workflow for the archaea using the model organism Halobacterium salinarum sp. NRC-1 and demonstrate its application for mapping the genome-wide binding sites of natively expressed transcription factors. This end-to-end pipeline is the first protocol for ChIP-seq in archaea, with methods and tools for each stage from gene tagging to data analysis and biological discovery. Genome-wide binding sites for transcription factors with many binding sites (TfbD) are identified with sensitivity, while retaining specificity in the identification the smaller regulons (bacteriorhodopsin-activator protein). Chromosomal tagging of target proteins with a compact epitope facilitates a standardized and cost-effective workflow that is compatible with high-throughput immunoprecipitation of natively expressed transcription factors. The Pique package, an open-source bioinformatics method, is presented for identification of binding events. Relative to ChIP-Chip and qPCR, this workflow offers a robust catalog of protein-DNA binding events with improved spatial resolution and significantly decreased cost. While this study focuses on the application of ChIP-seq in H. salinarum sp. NRC-1, our workflow can also be adapted for use in other archaea and bacteria with basic genetic tools.

As the third domain of life, archaea, like the eukarya and bacteria, must have robust DNA replication and repair complexes to ensure genome fidelity. Archaea moreover display a breadth of unique habitats and characteristics, and structural biologists increasingly appreciate these features. As archaea include extremophiles that can withstand diverse environmental stresses, they provide fundamental systems for understanding enzymes and pathways critical to genome integrity and stress responses. Such archaeal extremophiles provide critical data on the periodic table for life as well as on the biochemical, geochemical, and physical limitations to adaptive strategies allowing organisms to thrive under environmental stress relevant to determining the boundaries for life as we know it. Specifically, archaeal enzyme structures have informed the architecture and mechanisms of key DNA repair proteins and complexes. With added abilities to temperature-trap flexible complexes and reveal core domains of transient and dynamic complexes, these structures provide insights into mechanisms of maintaining genome integrity despite extreme environmental stress. The DNA damage response protein structures noted in this review therefore inform the basis for genome integrity in the face of environmental stress, with implications for all domains of life as well as for biomanufacturing, astrobiology, and medicine.

A comprehensive in silico analysis of 71 species representing the different taxonomic classes and physiological genre of the domain Archaea was performed. These organisms differed in their physiological attributes, particularly oxygen tolerance and energy metabolism. We explored the diversity and similarity in the codon usage pattern in the genes and genomes of these organisms, emphasizing on their core cellular pathways. Our thrust was to figure out whether there is any underlying similarity in the design of core pathways within these organisms. Analyses of codon utilization pattern, construction of hierarchical linear models of codon usage, expression pattern and codon pair preference pointed to the fact that, in the archaea there is a trend towards biased use of synonymous codons in the core cellular pathways and the Nc-plots appeared to display the physiological variations present within the different species. Our analyses revealed that aerobic species of archaea possessed a larger degree of freedom in regulating expression levels than could be accounted for by codon usage bias alone. This feature might be a consequence of their enhanced metabolic activities as a result of their adaptation to the relatively O2-rich environment. Species of archaea, which are related from the taxonomical viewpoint, were found to have striking similarities in their ORF structuring pattern. In the anaerobic species of archaea, codon bias was found to be a major determinant of gene expression. We have also detected a significant difference in the codon pair usage pattern between the whole genome and the genes related to vital cellular pathways, and it was not only species-specific but pathway specific too. This hints towards the structuring of ORFs with better decoding accuracy during translation. Finally, a codon-pathway interaction in shaping the codon design of pathways was observed where the transcription pathway exhibited a significantly different coding frequency signature.

Spindle- or lemon-shaped viruses infect archaea in diverse environments. Due to the highly pleomorphic nature of these virions, which can be found with cylindrical tails emanating from the spindle-shaped body, structural studies of these capsids have been challenging. We have determined the atomic structure of the capsid of Sulfolobus monocaudavirus 1, a virus that infects hosts living in nearly boiling acid. A highly hydrophobic protein, likely integrated into the host membrane before the virions assemble, forms 7 strands that slide past each other in both the tails and the spindle body. We observe the discrete steps that occur as the tail tubes expand, and these are due to highly conserved quasiequivalent interactions with neighboring subunits maintained despite significant diameter changes. Our results show how helical assemblies can vary their diameters, becoming nearly spherical to package a larger genome and suggest how all spindle-shaped viruses have evolved from archaeal rod-like viruses.

Amazonian soil microbial communities are known to be affected by the forest-to-pasture conversion, but the identity and metabolic potential of most of their organisms remain poorly characterized. To contribute to the understanding of these communities, here we describe metagenome-assembled genomes (MAGs) recovered from 12 forest and pasture soil metagenomes of the Brazilian Eastern Amazon. We obtained 11 forest and 30 pasture MAGs (≥50% of completeness and ≤10 % of contamination), distributed among two archaeal and 11 bacterial phyla. The taxonomic classification results suggest that most MAGs may represent potential novel microbial taxa. MAGs selected for further evaluation included members of Acidobacteriota, Actinobacteriota, Desulfobacterota_B, Desulfobacterota_F, Dormibacterota, Eremiobacterota, Halobacteriota, Proteobacteria, and Thermoproteota, thus revealing their roles in carbohydrate degradation and mercury detoxification as well as in the sulphur, nitrogen, and methane cycles. A methane-producing Archaea of the genus Methanosarcina was almost exclusively recovered from pasture soils, which can be linked to a sink-to-source shift after the forest-to-pasture conversion. The novel MAGs constitute an important resource to help us unravel the yet-unknown microbial diversity in Amazonian soils and its functional potential and, consequently, the responses of these microorganisms to land-use change.

The consensus genome sequence of a new member of the family Fuselloviridae designated as SMF1 (Sulfolobales Mexican fusellovirus 1) is presented. The complete circular genome was recovered from a metagenomic study of a Mexican hot spring. SMF1 exhibits an exceptional coding strand bias and a reduced set of fuselloviral core genes.

While biofilms formed by bacteria have received great attention due to their importance in pathogenesis, much less research has been focused on the biofilms formed by archaea. It has been known that extracellular filaments in archaea, such as type IV pili, hami, and cannulae, play a part in the formation of archaeal biofilms. We have used cryo-electron microscopy to determine the atomic structure of a previously uncharacterized class of archaeal surface filaments from hyperthermophilic Pyrobaculum calidifontis. These filaments, which we call archaeal bundling pili (ABP), assemble into highly ordered bipolar bundles. The bipolar nature of these bundles most likely arises from the association of filaments from at least two different cells. The component protein, AbpA, shows homology, both at the sequence and structural level, to the bacterial protein TasA, a major component of the extracellular matrix in bacterial biofilms, contributing to biofilm stability. We show that AbpA forms very stable filaments in a manner similar to the donor-strand exchange of bacterial TasA fibers and chaperone-usher pathway pili where a β-strand from one subunit is incorporated into a β-sheet of the next subunit. Our results reveal likely mechanistic similarities and evolutionary connection between bacterial and archaeal biofilms, and suggest that there could be many other archaeal surface filaments that are as yet uncharacterized.

Archaea are highly diverse and represent a primary life domain, but the majority of them remain uncultured. Currently, 16S rRNA phylogeny is widely used in archaeal taxonomy and diversity surveys. However, highly conserved sequence of 16S rRNA possibly results in generation of chimera in the amplicons and metagenome-assembled genomes (MAGs) and therefore limits its application. The newly developed phylogenomic approach has overcome these flaws, but it demands high-quality MAGs and intensive computation. In this study, we investigated the use of the archaeal transcription termination factor aCPSF1 in archaeal classification and diversity surveys. The phylogenetic analysis of 1,964 aCPSF1 orthologs retrieved from the available archaeal (meta)genomes resulted in convergent clustering patterns with those of archaeal phylogenomics and 16S rRNA phylogeny. The aCPSF1 phylogeny also displayed comparable clustering with the methanoarchaeal McrABG phylogeny and the haloarchaeal phylogenomics. Normalization of 779 aCPSF1 sequences including 261 from cultured archaeal species yielded a taxonomic ranking system with higher resolutions than that obtained with 16S rRNA for genus and species. Using the aCPSF1 taxonomy, 144 unclassified archaea in NCBI database were identified to various taxonomic ranks. Moreover, aCPSF1- and 16S rRNA-based surveys of the archaeal diversity in a sample from a South China Sea cold seep produced similar results. Our results demonstrate that aCPSF1 is an alternative archaeal phylogenetic marker, which exhibits higher resolution than 16S rRNA, and is more readily usable than phylogenomics in the taxonomic study of archaea. IMPORTANCE Archaea represent a unique type of prokaryote, which inhabit in various environments including extreme environments, and so define the boundary of biosphere, and play pivotal ecological roles, particularly in extreme environments. Since their discovery over 40 years ago, environmental archaea have been widely investigated using the 16S rRNA sequence comparison, and the recently developed phylogenomic approach because the majority of archaea are recalcitrant to laboratory cultivation. However, the highly conserved sequence of 16S rRNA and intensive bioinformatic computation of phylogenomics limit their applications in archaeal species delineation and diversity investigations. aCPSF1 is a ubiquitously distributed and vertically inherited transcription termination factor in archaea. In this study, we developed an aCPSF1-based archaeal taxonomic system which exhibits congruent phylogenic clustering patterns with archaeal phylogenomics and higher resolution than 16S rRNA in distinguishing archaea at lower taxonomic ranks. Therefore, aCPSF1 is a new phylogenetic marker in the taxonomic and diversity studies of archaea.

The Genome Taxonomy Database (GTDB; https://gtdb.ecogenomic.org) provides a phylogenetically consistent and rank normalized genome-based taxonomy for prokaryotic genomes sourced from the NCBI Assembly database. GTDB R06-RS202 spans 254 090 bacterial and 4316 archaeal genomes, a 270% increase since the introduction of the GTDB in November, 2017. These genomes are organized into 45 555 bacterial and 2339 archaeal species clusters which is a 200% increase since the integration of species clusters into the GTDB in June, 2019. Here, we explore prokaryotic diversity from the perspective of the GTDB and highlight the importance of metagenome-assembled genomes in expanding available genomic representation. We also discuss improvements to the GTDB website which allow tracking of taxonomic changes, easy assessment of genome assembly quality, and identification of genomes assembled from type material or used as species representatives. Methodological updates and policy changes made since the inception of the GTDB are then described along with the procedure used to update species clusters in the GTDB. We conclude with a discussion on the use of average nucleotide identities as a pragmatic approach for delineating prokaryotic species.

The genomes of all organisms throughout the tree of life are compacted and organized in chromatin by association of chromatin proteins. Eukaryotic genomes encode histones, which are assembled on the genome into octamers, yielding nucleosomes. Post-translational modifications of the histones, which occur mostly on their N-terminal tails, define the functional state of chromatin. Like eukaryotes, most archaeal genomes encode histones, which are believed to be involved in the compaction and organization of their genomes. Instead of discrete multimers, in vivo data suggest assembly of "nucleosomes" of variable size, consisting of multiples of dimers, which are able to induce repression of transcription. Based on these data and a model derived from X-ray crystallography, it was recently proposed that archaeal histones assemble on DNA into "endless" hypernucleosomes. In this review, we discuss the amino acid determinants of hypernucleosome formation and highlight differences with the canonical eukaryotic octamer. We identify archaeal histones differing from the consensus, which are expected to be unable to assemble into hypernucleosomes. Finally, we identify atypical archaeal histones with short N- or C-terminal extensions and C-terminal tails similar to the tails of eukaryotic histones, which are subject to post-translational modification. Based on the expected characteristics of these archaeal histones, we discuss possibilities of involvement of histones in archaeal transcription regulation.

The domain Archaea has historically been divided into two phyla, the Crenarchaeota and Euryarchaeota. Although regarded as members of the Crenarchaeota based on small subunit rRNA phylogeny, environmental genomics and efforts for cultivation have recently revealed two novel phyla/divisions in the Archaea; the 'Thaumarchaeota' and 'Korarchaeota'. Here, we show the genome sequence of Candidatus 'Caldiarchaeum subterraneum' that represents an uncultivated crenarchaeotic group. A composite genome was reconstructed from a metagenomic library previously prepared from a microbial mat at a geothermal water stream of a sub-surface gold mine. The genome was found to be clearly distinct from those of the known phyla/divisions, Crenarchaeota (hyperthermophiles), Euryarchaeota, Thaumarchaeota and Korarchaeota. The unique traits suggest that this crenarchaeotic group can be considered as a novel archaeal phylum/division. Moreover, C. subterraneum harbors an ubiquitin-like protein modifier system consisting of Ub, E1, E2 and small Zn RING finger family protein with structural motifs specific to eukaryotic system proteins, a system clearly distinct from the prokaryote-type system recently identified in Haloferax and Mycobacterium. The presence of such a eukaryote-type system is unprecedented in prokaryotes, and indicates that a prototype of the eukaryotic protein modifier system is present in the Archaea.

Methylmercury (MeHg) is a microbially produced neurotoxin derived from inorganic mercury (Hg), which accumulation in rice represents a major health concern to humans. However, the microbial control of MeHg dynamics in the environment remains elusive. Here, leveraging three rice paddy fields with distinct concentrations of Hg (Total Hg (THg): 0.21-513 mg kg-1 dry wt. soil; MeHg: 1.21-6.82 ng g-1 dry wt. soil), we resorted to metagenomics to determine the microbial determinants involved in MeHg production under contrasted contamination settings. We show that Hg methylating Archaea, along with methane-cycling genes, were enriched in severely contaminated paddy soils. Metagenome-resolved Genomes of novel putative Hg methylators belonging to Nitrospinota (UBA7883), with poorly resolved taxonomy despite high completeness, showed evidence of facultative anaerobic metabolism and adaptations to fluctuating redox potential. Furthermore, we found evidence of environmental filtering effects that influenced the phylogenies of not only hgcA genes under different THg concentrations, but also of two housekeeping genes, rpoB and glnA, highlighting the need for further experimental validation of whether THg drives the evolution of hgcAB. Finally, assessment of the genomic environment surrounding hgcAB suggests that this gene pair may be regulated by an archaeal toxin-antitoxin (TA) system, instead of the more frequently found arsR-like genes in bacterial methylators. This suggests the presence of distinct hgcAB regulation systems in bacteria and archaea. Our results support the emerging role of Archaea in MeHg cycling under mining-impacted environments and shed light on the differential control of the expression of genes involved in MeHg formation between Archaea and Bacteria.

Oceanic viruses that infect bacteria, or phages, are known to modulate host diversity, metabolisms, and biogeochemical cycling, while the viruses that infect marine Archaea remain understudied despite the critical ecosystem roles played by their hosts. Here we introduce "MArVD", for Metagenomic Archaeal Virus Detector, an annotation tool designed to identify putative archaeal virus contigs in metagenomic datasets. MArVD is made publicly available through the online iVirus analytical platform. Benchmarking analysis of MArVD showed it to be >99% accurate and 100% sensitive in identifying the 127 known archaeal viruses among the 12,499 viruses in the VirSorter curated dataset. Application of MArVD to 10 viral metagenomes from two depth profiles in the Eastern Tropical North Pacific (ETNP) oxygen minimum zone revealed 43 new putative archaeal virus genomes and large genome fragments ranging in size from 10 to 31 kb. Network-based classifications, which were consistent with marker gene phylogenies where available, suggested that these putative archaeal virus contigs represented six novel candidate genera. Ecological analyses, via fragment recruitment and ordination, revealed that the diversity and relative abundances of these putative archaeal viruses were correlated with oxygen concentration and temperature along two OMZ-spanning depth profiles, presumably due to structuring of the host Archaea community. Peak viral diversity and abundances were found in surface waters, where Thermoplasmata 16S rRNA genes are prevalent, suggesting these archaea as hosts in the surface habitats. Together these findings provide a baseline for identifying archaeal viruses in sequence datasets, and an initial picture of the ecology of such viruses in non-extreme environments.

Commonly used 16S rRNA gene primers do not detect the full range of archaeal diversity present in the vertebrate gut. As a result, several questions regarding the archaeal component of the gut microbiota remain, including which Archaea are host-associated, the specificities of such associations and the major factors influencing archaeal diversity. Using 16S rRNA gene amplicon sequencing with primers that specifically target Archaea, we obtained sufficient sequence data from 185 gastrointestinal samples collected from 110 vertebrate species that span five taxonomic classes (Mammalia, Aves, Reptilia, Amphibia and Actinopterygii), of which the majority were wild. We provide evidence for previously undescribed Archaea-host associations, including Bathyarchaeia and Methanothermobacter, the latter of which was prevalent among Aves and relatively abundant in species with higher body temperatures, although this association could not be decoupled from host phylogeny. Host phylogeny explained archaeal diversity more strongly than diet, while specific taxa were associated with both factors, and cophylogeny was significant and strongest for mammalian herbivores. Methanobacteria was the only class predicted to be present in the last common ancestors of mammals and all host species. Further analysis indicated that Archaea-Bacteria interactions have a limited effect on archaeal diversity. These findings expand our current understanding of Archaea-vertebrate associations.