Metagenome sequencing and 107 microbial genomes from seamount sediments along the Yap and Mariana trenches.

Microbes in the sediments across a series of seamounts along the island arc of the Yap and Mariana trenches were investigated by metagenome. In this study, we reconstructed 107 metagenome-assembled genomes (MAGs), including 100 bacteria and 7 archaea. All the MAGs exhibited >75% completeness and <10% contamination, with 26 MAGs being classified as 'nearly complete' (completeness >90%), while 50 falling within 80-90% range and 31 between 75-80% complete. Phylogenomic analysis revealed that 86% (n = 92) of these MAGs represented new taxa at different taxonomical levels. The species composition of these MAGs was most consistent with the previous reports, with the most abundant phyla being Proteobacteria (n = 39), Methylomirabilota (n = 27), and Nitrospirota (n = 7). These draft genomes provided novel data on species diversity and function in the seamount microbial community, which will provide reference data for extensive comparative genomic studies across crucial phylogenetic groups worldwide.

Intermolecular Gene Conversion for the Equalization of Genome Copies in the Polyploid Haloarchaeon Haloferax volcanii: Identification of Important Proteins.

The model haloarchaeon Haloferax volcanii is polyploid with about 20 copies of its major chromosome. Recently it has been described that highly efficient intermolecular gene conversion operates in H. volcanii to equalize the chromosomal copies. In the current study, 24 genes were selected that encode proteins with orthologs involved in gene conversion or homologous recombination in archaea, bacteria, or eukaryotes. Single gene deletion strains of 22 genes and a control gene were constructed in two parent strains for a gene conversion assay; only radA and radB were shown to be essential. Protoplast fusions were used to generate strains that were heterozygous for the gene HVO_2528, encoding an enzyme for carotinoid biosynthesis. It was revealed that a lack of six of the proteins did not influence the efficiency of gene conversion, while sixteen mutants had severe gene conversion defects. Notably, lack of paralogous proteins of gene families had very different effects, e.g., mutant Δrad25b had no phenotype, while mutants Δrad25a, Δrad25c, and Δrad25d were highly compromised. Generation of a quadruple rad25 and a triple sph deletion strain also indicated that the paralogs have different functions, in contrast to sph2 and sph4, which cannot be deleted simultaneously. There was no correlation between the severity of the phenotypes and the respective transcript levels under non-stressed conditions, indicating that gene expression has to be induced at the onset of gene conversion. Phylogenetic trees of the protein families Rad3/25, MutL/S, and Sph/SMC/Rad50 were generated to unravel the history of the paralogous proteins of H. volcanii. Taken together, unselected intermolecular gene conversion in H. volcanii involves at least 16 different proteins, the molecular roles of which can be studied in detail in future projects.

Phenotype-driven assessment of the ancestral trajectory of sulfur biooxidation in the thermoacidophilic archaea Sulfolobaceae.

Certain members of the family Sulfolobaceae represent the only archaea known to oxidize elemental sulfur, and their evolutionary history provides a framework to understand the development of chemolithotrophic growth by sulfur oxidation. Here, we evaluate the sulfur oxidation phenotype of Sulfolobaceae species and leverage comparative genomic and transcriptomic analysis to identify the key genes linked to sulfur oxidation. Metabolic engineering of the obligate heterotroph Sulfolobus acidocaldarius revealed that the known cytoplasmic components of sulfur oxidation alone are not sufficient to drive prolific sulfur oxidation. Imaging analysis showed that Sulfolobaceae species maintain proximity to the sulfur surface but do not necessarily contact the substrate directly. This indicates that a soluble form of sulfur must be transported to initiate cytoplasmic sulfur oxidation. Conservation patterns and transcriptomic response implicate an extracellular tetrathionate hydrolase and putative thiosulfate transporter in a newly proposed mechanism of sulfur acquisition in the Sulfolobaceae.IMPORTANCESulfur is one of the most abundant elements on earth (2.9% by mass), so it makes sense that the earliest biology found a way to use sulfur to create and sustain life. However, beyond evolutionary significance, sulfur and the molecules it comprises have important technological significance, not only in chemicals such as sulfuric acid and in pyritic ores containing critical metals but also as a waste product from oil and gas production. The thermoacidophilic Sulfolobaceae are unique among the archaea as sulfur oxidizers. The trajectory for how sulfur biooxidation arose and evolved can be traced using experimental and bioinformatic analyses of the available genomic data set. Such analysis can also inform the process by which extracellular sulfur is acquired and transported by thermoacidophilic archaea, a phenomenon that is critical to these microorganisms but has yet to be elucidated.

Asgard archaea modulate potential methanogenesis substrates in wetland soil.

The roles of Asgard archaea in eukaryogenesis and marine biogeochemical cycles are well studied, yet their contributions in soil ecosystems remain unknown. Of particular interest are Asgard archaeal contributions to methane cycling in wetland soils. To investigate this, we reconstructed two complete genomes for soil-associated Atabeyarchaeia, a new Asgard lineage, and a complete genome of Freyarchaeia, and predicted their metabolism in situ. Metatranscriptomics reveals expression of genes for [NiFe]-hydrogenases, pyruvate oxidation and carbon fixation via the Wood-Ljungdahl pathway. Also expressed are genes encoding enzymes for amino acid metabolism, anaerobic aldehyde oxidation, hydrogen peroxide detoxification and carbohydrate breakdown to acetate and formate. Overall, soil-associated Asgard archaea are predicted to include non-methanogenic acetogens, highlighting their potential role in carbon cycling in terrestrial environments.

The parasitic lifestyle of an archaeal symbiont.

DPANN archaea are a diverse group of microorganisms characterised by small cells and reduced genomes. To date, all cultivated DPANN archaea are ectosymbionts that require direct cell contact with an archaeal host species for growth and survival. However, these interactions and their impact on the host species are poorly understood. Here, we show that a DPANN archaeon (Candidatus Nanohaloarchaeum antarcticus) engages in parasitic interactions with its host (Halorubrum lacusprofundi) that result in host cell lysis. During these interactions, the nanohaloarchaeon appears to enter, or be engulfed by, the host cell. Our results provide experimental evidence for a predatory-like lifestyle of an archaeon, suggesting that at least some DPANN archaea may have roles in controlling host populations and their ecology.

Microdiversity in marine pelagic ammonia-oxidizing archaeal populations in a Mediterranean long-read metagenome.

The knowledge of the different population-level processes operating within a species, and the genetic variability of the individual prokaryotic genomes, is key to understanding the adaptability of microbial populations. Here, we characterized the flexible genome of ammonia-oxidizing archaeal (AOA) populations using a metagenomic recruitment approach and long-read (PacBio HiFi) metagenomic sequencing. In the lower photic zone of the western Mediterranean Sea (75 m deep), the genomes Nitrosopelagicus brevis CN25 and Nitrosopumilus catalinensis SPOT1 had the highest recruitment values among available complete AOA genomes. They were used to analyse the diversity of flexible genes (variable from strain to strain) by examining the long-reads located within the flexible genomic islands (fGIs) identified by their under-recruitment. Both AOA genomes had a large fGI involved in the glycosylation of exposed structures, highly variable, and rich in glycosyltransferases. N. brevis had two fGIs related to the transport of phosphorus and ammonium respectively. N. catalinensis had fGIs involved in phosphorus transportation and metal uptake. A fGI5 previously reported as 'unassigned function' in N. brevis could be associated with defense. These findings demonstrate that the microdiversity of marine microbe populations, including AOA, can be effectively characterized using an approach that incorporates third-generation sequencing metagenomics.

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.

Borg extrachromosomal elements of methane-oxidizing archaea have conserved and expressed genetic repertoires.

Borgs are huge extrachromosomal elements (ECE) of anaerobic methane-consuming "Candidatus Methanoperedens" archaea. Here, we used nanopore sequencing to validate published complete genomes curated from short reads and to reconstruct new genomes. 13 complete and four near-complete linear genomes share 40 genes that define a largely syntenous genome backbone. We use these conserved genes to identify new Borgs from peatland soil and to delineate Borg phylogeny, revealing two major clades. Remarkably, Borg genes encoding nanowire-like electron-transferring cytochromes and cell surface proteins are more highly expressed than those of host Methanoperedens, indicating that Borgs augment the Methanoperedens activity in situ. We reconstructed the first complete 4.00 Mbp genome for a Methanoperedens that is inferred to be a Borg host and predicted its methylation motifs, which differ from pervasive TC and CC methylation motifs of the Borgs. Thus, methylation may enable Methanoperedens to distinguish their genomes from those of Borgs. Very high Borg to Methanoperedens ratios and structural predictions suggest that Borgs may be capable of encapsulation. The findings clearly define Borgs as a distinct class of ECE with shared genomic signatures, establish their diversification from a common ancestor with genetic inheritance, and raise the possibility of periodic existence outside of host cells.

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.

Halobacterium yunchengense sp. nov., Natronomonas amylolytica sp. nov., Halorientalis halophila sp. nov., Halobellus salinisoli sp. nov., halophilic archaea isolated from a saline lake and inland saline soil.

Four halophilic archaeal strains YCN1T, YCN58T, LT38T, and LT62T were isolated from Yuncheng Salt Lake (Shanxi, China) and Tarim Basin (Xinjiang, China), respectively. Phylogenetic and phylogenomic analyses showed that these four strains tightly cluster with related species of Halobacterium, Natronomonas, Halorientalis, and Halobellus, respectively. The AAI, ANI, and dDDH values between these four strains and their related species of respective genera were lower than the proposed threshold values for species delineation. Strains YCN1T, YCN58T, LT38T, and LT62T could be differentiated from the current species of Halobacterium, Natronomonas, Halorientalis, and Halobellus, respectively, based on the comparison of diverse phenotypic characteristics. The polar lipid profiles of these four strains were closely similar to those of respective relatives within the genera Halobacterium, Natronomonas, Halorientalis, and Halobellus, respectively. The phenotypic, phylogenetic, and genome-based analyses indicated that strains YCN1T, YCN58T, LT38T, and LT62T represent respective novel species within the genera Halobacterium, Natronomonas, Halorentalis, and Halobellus, for which the names Halobacterium yunchengense sp. nov., Natronomonas amylolytica sp. nov., Halorientalis halophila sp. nov., and Halobellus salinisoli sp. nov. are proposed, respectively.

How Do Thermophiles Organize Their Genomes?

All cells must maintain the structural and functional integrity of the genome under a wide range of environments. High temperatures pose a formidable challenge to cells by denaturing the DNA double helix, causing chemical damage to DNA, and increasing the random thermal motion of chromosomes. Thermophiles, predominantly classified as bacteria or archaea, exhibit an exceptional capacity to mitigate these detrimental effects and prosper under extreme thermal conditions, with some species tolerating temperatures higher than 100°C. Their genomes are mainly characterized by the presence of reverse gyrase, a unique topoisomerase that introduces positive supercoils into DNA. This enzyme has been suggested to maintain the genome integrity of thermophiles by limiting DNA melting and mediating DNA repair. Previous studies provided significant insights into the mechanisms by which NAPs, histones, SMC superfamily proteins, and polyamines affect the 3D genomes of thermophiles across different scales. Here, I discuss current knowledge of the genome organization in thermophiles and pertinent research questions for future investigations.

A non-methanogenic archaeon within the order Methanocellales.

Serpentinization, a geochemical process found on modern and ancient Earth, provides an ultra-reducing environment that can support microbial methanogenesis and acetogenesis. Several groups of archaea, such as the order Methanocellales, are characterized by their ability to produce methane. Here, we generate metagenomic sequences from serpentinized springs in The Cedars, California, and construct a circularized metagenome-assembled genome of a Methanocellales archaeon, termed Met12, that lacks essential methanogenesis genes. The genome includes genes for an acetyl-CoA pathway, but lacks genes encoding methanogenesis enzymes such as methyl-coenzyme M reductase, heterodisulfide reductases and hydrogenases. In situ transcriptomic analyses reveal high expression of a multi-heme c-type cytochrome, and heterologous expression of this protein in a model bacterium demonstrates that it is capable of accepting electrons. Our results suggest that Met12, within the order Methanocellales, is not a methanogen but a CO2-reducing, electron-fueled acetogen without electron bifurcation.

The biosynthesis mechanism of bacterioruberin in halophilic archaea revealed by genome and transcriptome analysis.

Halophilic archaea are promising microbial cell factories for bacterioruberin (BR) production. BR is a natural product with multi-bioactivities, allowing potential application in many fields. In the previous work, a haloarchaeon Halorubrum sp. HRM-150 with a high proportion of BR (about 85%) was isolated, but the low yield impeded its large-scale production. This work figured out BR synthesis characteristics and mechanisms, and proposed strategies for yield improvement. First, glucose (10 g/L) and tryptone (15 g/L) were tested to be better sources for BR production. Besides, the combination of glucose and starch achieved the diauxic growth, and the biomass and BR productivity increased by 85% and 54% than using glucose. Additionally, this work first proposed the BR synthesis pattern, which differs from that of other carotenoids. As a structural component of cell membranes, the BR synthesis is highly coupled with growth, which was most active in the logarithm phase. Meanwhile, the osmotic down shock at the logarithm phase could increase the BR productivity without sacrificing the biomass. Moreover, the de-novo pathway for BR synthesis with a key gene of lyeJ, and its competitive pathways (notably tetraether lipids and retinal) were revealed through genome, transcriptome, and osmotic down shock. Therefore, the BR yield is expected to be improved through mutant construction, such as the overexpression of key gene lyeJ and the knockout of competitive genes, which need to be further explored. The findings will contribute to a better understanding of the metabolism mechanism in haloarchaea and the development of haloarchaea as microbial cell factories. IMPORTANCE: Recent studies have revealed that halophilic microorganism is a promising microbial factory for the next-generation industrialization. Among them, halophilic archaea are advantageous as microbial factories due to their low contamination risk and low freshwater consumption. The halophilic archaea usually accumulate long chain C50 carotenoids, which are barely found in other organisms. Bacterioruberin (BR), the major C50 carotenoid, has multi-bioactivities, allowing potential application in food, cosmetic, and biomedical industries. However, the low yield impedes its large-scale application. This work figured out the BR synthesis characteristics and mechanism, and proposed several strategies for BR yield improvement, encouraging halophilic archaea to function as microbial factories for BR production. Meanwhile, the archaea have special evolutionary status and unique characteristics in taxonomy, the revelation of BR biosynthesis mechanism is beneficial for a better understanding of archaea.

EzBioCloud: a genome-driven database and platform for microbiome identification and discovery.

With the continued evolution of DNA sequencing technologies, the role of genome sequence data has become more integral in the classification and identification of Bacteria and Archaea. Six years after introducing EzBioCloud, an integrated platform representing the taxonomic hierarchy of Bacteria and Archaea through quality-controlled 16S rRNA gene and genome sequences, we present an updated version, that further refines and expands its capabilities. The current update recognizes the growing need for accurate taxonomic information as defining a species increasingly relies on genome sequence comparisons. We also incorporated an advanced strategy for addressing underrepresented or less studied lineages, bolstering the comprehensiveness and accuracy of our database. Our rigorous quality control protocols remain, where whole-genome assemblies from the NCBI Assembly Database undergo stringent screening to remove low-quality sequence data. These are then passed through our enhanced identification bioinformatics pipeline which initiates a 16S rRNA gene similarity search and then calculates the average nucleotide identity (ANI). For genome sequences lacking a 16S rRNA sequence and without a closely related genomic representative for ANI calculation, we apply a different ANI approach using bacterial core genes for improved taxonomic placement (core gene ANI, cgANI). Because of the increase in genome sequences available in NCBI and our newly introduced cgANI method, EzBioCloud now encompasses a total of 109 835 species, of which 21 964 have validly published names. 47 896 are candidate species identified either through 16S rRNA sequence similarity (phylotypes) or through whole genome ANI (genomospecies), and the remaining 39 975 were positioned in the taxonomic tree by cgANI (species clusters). Our EzBioCloud database is accessible at www.ezbiocloud.net/db.

Minimal and hybrid hydrogenases are active from archaea.

Microbial hydrogen (H2) cycling underpins the diversity and functionality of diverse anoxic ecosystems. Among the three evolutionarily distinct hydrogenase superfamilies responsible, [FeFe] hydrogenases were thought to be restricted to bacteria and eukaryotes. Here, we show that anaerobic archaea encode diverse, active, and ancient lineages of [FeFe] hydrogenases through combining analysis of existing and new genomes with extensive biochemical experiments. [FeFe] hydrogenases are encoded by genomes of nine archaeal phyla and expressed by H2-producing Asgard archaeon cultures. We report an ultraminimal hydrogenase in DPANN archaea that binds the catalytic H-cluster and produces H2. Moreover, we identify and characterize remarkable hybrid complexes formed through the fusion of [FeFe] and [NiFe] hydrogenases in ten other archaeal orders. Phylogenetic analysis and structural modeling suggest a deep evolutionary history of hybrid hydrogenases. These findings reveal new metabolic adaptations of archaea, streamlined H2 catalysts for biotechnological development, and a surprisingly intertwined evolutionary history between the two major H2-metabolizing enzymes.

A genome catalogue of mercury-methylating bacteria and archaea from sediments of a boreal river facing human disturbances.

Methyl mercury, a toxic compound, is produced by anaerobic microbes and magnifies in aquatic food webs, affecting the health of animals and humans. The exploration of mercury methylators based on genomes is still limited, especially in the context of river ecosystems. To address this knowledge gap, we developed a genome catalogue of potential mercury-methylating microorganisms. This was based on the presence of hgcAB from the sediments of a river affected by two run-of-river hydroelectric dams, logging activities and a wildfire. Through the use of genome-resolved metagenomics, we discovered a unique and diverse group of mercury methylators. These were dominated by members of the metabolically versatile Bacteroidota and were particularly rich in microbes that ferment butyrate. By comparing the diversity and abundance of mercury methylators between sites subjected to different disturbances, we found that ongoing disturbances, such as the input of organic matter related to logging activities, were particularly conducive to the establishment of a mercury-methylating niche. Finally, to gain a deeper understanding of the environmental factors that shape the diversity of mercury methylators, we compared the mercury-methylating genome catalogue with the broader microbial community. The results suggest that mercury methylators respond to environmental conditions in a manner similar to the overall microbial community. Therefore, it is crucial to interpret the diversity and abundance of mercury methylators within their specific ecological context.

Diverse and unconventional methanogens, methanotrophs, and methylotrophs in metagenome-assembled genomes from subsurface sediments of the Slate River floodplain, Crested Butte, CO, USA.

We use metagenome-assembled genomes (MAGs) to understand single-carbon (C1) compound-cycling-particularly methane-cycling-microorganisms in montane riparian floodplain sediments. We generated 1,233 MAGs (>50% completeness and <10% contamination) from 50- to 150-cm depth below the sediment surface capturing the transition between oxic, unsaturated sediments and anoxic, saturated sediments in the Slate River (SR) floodplain (Crested Butte, CO, USA). We recovered genomes of putative methanogens, methanotrophs, and methylotrophs (n = 57). Methanogens, found only in deep, anoxic depths at SR, originate from three different clades (Methanoregulaceae, Methanotrichaceae, and Methanomassiliicoccales), each with a different methanogenesis pathway; putative methanotrophic MAGs originate from within the Archaea (Candidatus Methanoperedens) in anoxic depths and uncultured bacteria (Ca. Binatia) in oxic depths. Genomes for canonical aerobic methanotrophs were not recovered. Ca. Methanoperedens were exceptionally abundant (~1,400× coverage, >50% abundance in the MAG library) in one sample that also contained aceticlastic methanogens, indicating a potential C1/methane-cycling hotspot. Ca. Methylomirabilis MAGs from SR encode pathways for methylotrophy but do not harbor methane monooxygenase or nitrogen reduction genes. Comparative genomic analysis supports that one clade within the Ca. Methylomirabilis genus is not methanotrophic. The genetic potential for methylotrophy was widespread, with over 10% and 19% of SR MAGs encoding a methanol dehydrogenase or substrate-specific methyltransferase, respectively. MAGs from uncultured Thermoplasmata archaea in the Ca. Gimiplasmatales (UBA10834) contain pathways that may allow for anaerobic methylotrophic acetogenesis. Overall, MAGs from SR floodplain sediments reveal a potential for methane production and consumption in the system and a robust potential for methylotrophy.IMPORTANCEThe cycling of carbon by microorganisms in subsurface environments is of particular relevance in the face of global climate change. Riparian floodplain sediments contain high organic carbon that can be degraded into C1 compounds such as methane, methanol, and methylamines, the fate of which depends on the microbial metabolisms present as well as the hydrological conditions and availability of oxygen. In the present study, we generated over 1,000 MAGs from subsurface sediments from a montane river floodplain and recovered genomes for microorganisms that are capable of producing and consuming methane and other C1 compounds, highlighting a robust potential for C1 cycling in subsurface sediments both with and without oxygen. Archaea from the Ca. Methanoperedens genus were exceptionally abundant in one sample, indicating a potential C1/methane-cycling hotspot in the Slate River floodplain system.

Investigating the nature of prokaryotic genomic island locations within a genome.

Horizontal gene transfer (HGT) is a powerful evolutionary force that considerably shapes the structure of prokaryotic genomes and is associated with genomic islands (GIs). A GI is a DNA segment composed of transferred genes that can be found within a prokaryotic genome, obtained through HGT. Much research has focused on detecting GIs in genomes, but here we pursue a new course, which is identifying possible preferred locations of GIs in the prokaryotic genome. Here, we identify the locations of the GIs within prokaryotic genomes to examine patterns in those locations. Prokaryotic GIs were analyzed according to the genome structure that they are located in, whether it be a circular or a linear genome. The analytical investigations employed are: (1) studying the GI locations in relation to the origin of replication (oriC); (2) exploring the distances between GIs; and (3) determining the distribution of GIs across the genomes. For each of the investigations, the analysis was performed on all of the GIs in the data set. Moreover, to void bias caused by the distribution of the genomes represented, the GIs in one genome from each species and the GIs of the most frequent species are also analyzed. Overall, the results showed that there are preferred sites for the GIs in the genome. In the linear genomes, these sites are usually located in the oriC region and terminus region, while in the circular genomes, they are located solely in the terminus region. These results also showed that the distance distribution between the GIs is almost exponential, which proves that GIs have preferred sites within genomes. The oriC and termniuns are preferred sites for the GIs and a possible natural explanation for this could be connected to the content of the oriC region. Moreover, the content of the GIs in terms of its protein families was studied and the results demonstrated that the majority of frequent protein families are close to identical in each section.

Fusion/fission protein family identification in Archaea.

The majority of newly discovered archaeal lineages remain without a cultivated representative, but scarce experimental data from the cultivated organisms show that they harbor distinct functional repertoires. To unveil the ecological as well as evolutionary impact of Archaea from metagenomics, new computational methods need to be developed, followed by in-depth analysis. Among them is the genome-wide protein fusion screening performed here. Natural fusions and fissions of genes not only contribute to microbial evolution but also complicate the correct identification and functional annotation of sequences. The products of these processes can be defined as fusion (or composite) proteins, the ones consisting of two or more domains originally encoded by different genes and split proteins, and the ones originating from the separation of a gene in two (fission). Fusion identifications are required for proper phylogenetic reconstructions and metabolic pathway completeness assessments, while mappings between fused and unfused proteins can fill some of the existing gaps in metabolic models. In the archaeal genome-wide screening, more than 1,900 fusion/fission protein clusters were identified, belonging to both newly sequenced and well-studied lineages. These protein families are mainly associated with different types of metabolism, genetic, and cellular processes. Moreover, 162 of the identified fusion/fission protein families are archaeal specific, having no identified fused homolog within the bacterial domain. Our approach was validated by the identification of experimentally characterized fusion/fission cases. However, around 25% of the identified fusion/fission families lack functional annotations for both composite and split states, showing the need for experimental characterization in Archaea.IMPORTANCEGenome-wide fusion screening has never been performed in Archaea on a broad taxonomic scale. The overlay of multiple computational techniques allows the detection of a fine-grained set of predicted fusion/fission families, instead of rough estimations based on conserved domain annotations only. The exhaustive mapping of fused proteins to bacterial organisms allows us to capture fusion/fission families that are specific to archaeal biology, as well as to identify links between bacterial and archaeal lineages based on cooccurrence of taxonomically restricted proteins and their sequence features. Furthermore, the identification of poorly characterized lineage-specific fusion proteins opens up possibilities for future experimental and computational investigations. This approach enhances our understanding of Archaea in general and provides potential candidates for in-depth studies in the future.