Updated clusters of orthologous genes for Archaea: a complex ancestor of the Archaea and the byways of horizontal gene transfer

Long range segmentation of prokaryotic genomes by gene age and functionality

Bacterial and archaeal genomes encompass numerous operons that typically consist of two to five genes. On larger scales, however, gene order is poorly conserved through the evolution of prokaryotes. Nevertheless, non-random localization of different classes of genes on prokaryotic chromosomes could reflect important functional and evolutionary constraints. We explored the patterns of genomic localization of evolutionarily conserved (ancient) and variable (young) genes across the diversity of bacteria and archaea. Nearly all bacterial and archaeal chromosomes were found to encompass large segments of 100-300 kilobases that were significantly enriched in either ancient or young genes. Similar clustering of genes with lethal knockout phenotype (essential genes) was observed as well. Mathematical modeling of genome evolution suggests that this long-range gene clustering in prokaryotic chromosomes reflects perpetual genome rearrangement driven by a combination of selective and neutral processes rather than evolutionary conservation.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Computational analysis of genes with lethal knockout phenotype and prediction of essential genes in archaea

The identification of microbial genes essential for survival as those with lethal knockout phenotype (LKP) is a common strategy for functional interrogation of genomes. However, interpretation of the LKP is complicated because a substantial fraction of the genes with this phenotype remains poorly functionally characterized. Furthermore, many genes can exhibit LKP not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes (conditionally essential genes). We analyzed the sets of LKP genes for two archaea, Methanococcus maripaludis and Sulfolobus islandicus, using a variety of computational approaches aiming to differentiate between essential and conditionally essential genes and to predict at least a general function for as many of the proteins encoded by these genes as possible. This analysis allowed us to predict the functions of several LKP genes including previously uncharacterized subunit of the GINS protein complex with an essential function in genome replication and of the KEOPS complex that is responsible for an essential tRNA modification as well as GRP protease implicated in protein quality control. Additionally, several novel antitoxins (conditionally essential genes) were predicted, and this prediction was experimentally validated by showing that the deletion of these genes together with the adjacent genes apparently encoding the cognate toxins caused no growth defect. We applied principal component analysis based on sequence and comparative genomic features showing that this approach can separate essential genes from conditionally essential ones and used it to predict essential genes in other archaeal genomes.IMPORTANCEOnly a relatively small fraction of the genes in any bacterium or archaeon is essential for survival as demonstrated by the lethal effect of their disruption. The identification of essential genes and their functions is crucial for understanding fundamental cell biology. However, many of the genes with a lethal knockout phenotype remain poorly functionally characterized, and furthermore, many genes can exhibit this phenotype not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes. We applied state-of-the-art computational methods to predict the functions of a number of uncharacterized genes with the lethal knockout phenotype in two archaeal species and developed a computational approach to predict genes involved in essential functions. These findings advance the current understanding of key functionalities of archaeal cells.

TroR is the primary regulator of the iron homeostasis transcription network in the halophilic archaeon Haloferax volcanii

Maintaining the intracellular iron concentration within the homeostatic range is vital to meet cellular metabolic needs and reduce oxidative stress. Previous research revealed that the haloarchaeon Halobacterium salinarum encodes four diphtheria toxin repressor (DtxR) family transcription factors (TFs) that together regulate the iron response through an interconnected transcriptional regulatory network (TRN). However, the conservation of the TRN and the metal specificity of DtxR TFs remained poorly understood. Here we identified and characterized the TRN of Haloferax volcanii for comparison. Genetic analysis demonstrated that Hfx. volcanii relies on three DtxR transcriptional regulators (Idr, SirR, and TroR), with TroR as the primary regulator of iron homeostasis. Bioinformatics and molecular approaches revealed that TroR binds a conserved cis-regulatory motif located ∼100 nt upstream of the start codon of iron-related target genes. Transcriptomics analysis demonstrated that, under conditions of iron sufficiency, TroR repressed iron uptake and induced iron storage mechanisms. TroR repressed the expression of one other DtxR TF, Idr. This reduced DtxR TRN complexity relative to that of Hbt. salinarum appeared correlated with natural variations in iron availability. Based on these data, we hypothesize that variable environmental conditions such as iron availability appear to select for increasing TRN complexity.

The Raphidiopsis (= Cylindrospermopsis) raciborskii pangenome updated: Two new metagenome-assembled genomes from the South American clade

Two Raphidiopsis (=Cylindrospermopsis) raciborskii metagenome-assembled genomes (MAGs) were recovered from two freshwater metagenomic datasets sampled in 2011 and 2012 in Pampulha Lake, a hypereutrophic, artificial, shallow reservoir, located in the city of Belo Horizonte (MG), Brazil. Since the late 1970s, the lake has undergone increasing eutrophication pressure, due to wastewater input, leading to the occurrence of frequent cyanobacterial blooms. The major difference observed between PAMP2011 and PAMP2012 MAGs was the lack of the saxitoxin gene cluster in PAMP2012, which also presented a smaller genome, while PAMP2011 presented the complete sxt cluster and all essential proteins and clusters. The pangenome analysis was performed with all Raphidiopsis/Cylindrospermopsis genomes available at NCBI to date, with the addition of PAMP2011 and PAMP2012 MAGs (All33 subset), but also without the South American strains (noSA subset), and only among the South American strains (SA10 and SA8 subsets). We observed a substantial increase in the core genome size for the 'noSA' subset, in comparison to 'All33' subset, and since the core genome reflects the closeness among the pangenome members, the results strongly suggest that the conservation level of the essential gene repertoire seems to be affected by the geographic origin of the strains being analyzed, supporting the existence of a distinct SA clade. The Raphidiopsis pangenome comprised a total of 7943 orthologous protein clusters, and the two new MAGs increased the pangenome size by 11%. The pangenome based phylogenetic relationships among the 33 analyzed genomes showed that the SA genomes clustered together with 99% bootstrap support, reinforcing the metabolic particularity of the Raphidiopsis South American clade, related to its saxitoxin producing unique ability, while also indicating a different evolutionary history due to its geographic isolation.

Multi-environment ecogenomics analysis of the cosmopolitan phylum Gemmatimonadota

Gemmatimonadota is a diverse bacterial phylum commonly found in environments such as soils, rhizospheres, fresh waters, and sediments. So far, the phylum contains just six cultured species (five of them sequenced), which limits our understanding of their diversity and metabolism. Therefore, we analyzed over 400 metagenome-assembled genomes (MAGs) and 5 culture-derived genomes representing Gemmatimonadota from various aquatic environments, hydrothermal vents, sediments, soils, and host-associated (with marine sponges and coral) species. The principal coordinate analysis based on the presence/absence of genes in Gemmatimonadota genomes and phylogenomic analysis documented that marine and host-associated Gemmatimonadota were the most distant from freshwater and wastewater species. A smaller genome size and coding sequences (CDS) number reduction were observed in marine MAGs, pointing to an oligotrophic environmental adaptation. Several metabolic pathways are restricted to specific environments. For example, genes for anoxygenic phototrophy were found only in freshwater, wastewater, and soda lake sediment genomes. There were several genomes from soda lake sediments and wastewater containing type IC/ID ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO). Various genomes from wastewater harbored bacterial type II RuBisCO, whereas RuBisCO-like protein was found in genomes from fresh waters, soil, host-associated, and marine sediments. Gemmatimonadota does not contain nitrogen fixation genes; however, the nosZ gene, involved in the reduction of N2O, was present in genomes from most environments, missing only in marine water and host-associated Gemmatimonadota. The presented data suggest that Gemmatimonadota evolved as an organotrophic species relying on aerobic respiration and then remodeled its genome inventory when adapting to particular environments. IMPORTANCE Gemmatimonadota is a rarely studied bacterial phylum consisting of a handful of cultured species. Recent culture-independent studies documented that these organisms are distributed in many environments, including soil, marine, fresh, and waste waters. However, due to the lack of cultured species, information about their metabolic potential and environmental role is scarce. Therefore, we collected Gemmatimonadota metagenome-assembled genomes (MAGs) from different habitats and performed a systematic analysis of their genomic characteristics and metabolic potential. Our results show how Gemmatimonadota have adapted their genomes to different environments.

Metabolic versatility of Caldarchaeales from geothermal features of Hawai'i and Chile as revealed by five metagenome-assembled genomes

Members of the archaeal order Caldarchaeales (previously the phylum Aigarchaeota) are poorly sampled and are represented in public databases by relatively few genomes. Additional representative genomes will help resolve their placement among all known members of Archaea and provide insights into their roles in the environment. In this study, we analyzed 16S rRNA gene amplicons belonging to the Caldarchaeales that are available in public databases, which demonstrated that archaea of the order Caldarchaeales are diverse, widespread, and most abundant in geothermal habitats. We also constructed five metagenome-assembled genomes (MAGs) of Caldarchaeales from two geothermal features to investigate their metabolic potential and phylogenomic position in the domain Archaea. Two of the MAGs were assembled from microbial community DNA extracted from fumarolic lava rocks from Mauna Ulu, Hawai'i, and three were assembled from DNA obtained from hot spring sinters from the El Tatio geothermal field in Chile. MAGs from Hawai'i are high quality bins with completeness >95% and contamination <1%, and one likely belongs to a novel species in a new genus recently discovered at a submarine volcano off New Zealand. MAGs from Chile have lower completeness levels ranging from 27 to 70%. Gene content of the MAGs revealed that these members of Caldarchaeales are likely metabolically versatile and exhibit the potential for both chemoorganotrophic and chemolithotrophic lifestyles. The wide array of metabolic capabilities exhibited by these members of Caldarchaeales might help them thrive under diverse harsh environmental conditions. All the MAGs except one from Chile harbor putative prophage regions encoding several auxiliary metabolic genes (AMGs) that may confer a fitness advantage on their Caldarchaeales hosts by increasing their metabolic potential and make them better adapted to new environmental conditions. Phylogenomic analysis of the five MAGs and over 3,000 representative archaeal genomes showed the order Caldarchaeales forms a monophyletic group that is sister to the clade comprising the orders Geothermarchaeales (previously Candidatus Geothermarchaeota), Conexivisphaerales and Nitrososphaerales (formerly known as Thaumarchaeota), supporting the status of Caldarchaeales members as a clade distinct from the Thaumarchaeota.

In-depth genome and pan-genome analysis of a metal-resistant bacterium Pseudomonas parafulva OS-1

A metal-resistant bacterium Pseudomonas parafulva OS-1 was isolated from waste-contaminated soil in Ranchi City, India. The isolated strain OS-1 showed its growth at 25-45°C, pH 5.0-9.0, and in the presence of ZnSO₄ (upto 5 mM). Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain OS-1 belonged to the genus Pseudomonas and was most closely related to parafulva species. To unravel the genomic features, we sequenced the complete genome of P. parafulva OS-1 using Illumina HiSeq 4,000 sequencing platform. The results of average nucleotide identity (ANI) analysis indicated the closest similarity of OS-1 to P. parafulva PRS09-11288 and P. parafulva DTSP2. The metabolic potential of P. parafulva OS-1 based on Clusters of Othologous Genes (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) indicated a high number of genes related to stress protection, metal resistance, and multiple drug-efflux, etc., which is relatively rare in P. parafulva strains. Compared with other parafulva strains, P. parafulva OS-1 was found to have the unique β-lactam resistance and type VI secretion system (T6SS) gene. Additionally, its genomes encode various CAZymes such as glycoside hydrolases and other genes associated with lignocellulose breakdown, suggesting that strain OS-1 have strong biomass degradation potential. The presence of genomic complexity in the OS-1 genome indicates that horizontal gene transfer (HGT) might happen during evolution. Therefore, genomic and comparative genome analysis of parafulva strains is valuable for further understanding the mechanism of resistance to metal stress and opens a perspective to exploit a newly isolated bacterium for biotechnological applications.

Nutrition or nature: using elementary flux modes to disentangle the complex forces shaping prokaryote pan-genomes

Background

Microbial pan-genomes are shaped by a complex combination of stochastic and deterministic forces. Even closely related genomes exhibit extensive variation in their gene content. Understanding what drives this variation requires exploring the interactions of gene products with each other and with the organism’s external environment. However, to date, conceptual models of pan-genome dynamics often represent genes as independent units and provide limited information about their mechanistic interactions.

Results

We simulated the stochastic process of gene-loss using the pooled genome-scale metabolic reaction networks of 46 taxonomically diverse bacterial and archaeal families as proxies for their pan-genomes. The frequency by which reactions are retained in functional networks when stochastic gene loss is simulated in diverse environments allowed us to disentangle the metabolic reactions whose presence depends on the metabolite composition of the external environment (constrained by “nutrition”) from those that are independent of the environment (constrained by “nature”). By comparing the frequency of reactions from the first group with their observed frequencies in bacterial and archaeal families, we predicted the metabolic niches that shaped the genomic composition of these lineages. Moreover, we found that the lineages that were shaped by a more diverse metabolic niche also occur in more diverse biomes as assessed by global environmental sequencing datasets.

Conclusion

We introduce a computational framework for analyzing and interpreting pan-reactomes that provides novel insights into the ecological and evolutionary drivers of pan-genome dynamics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02052-3.

Elucidating the picocyanobacteria salinity divide through ecogenomics of new freshwater isolates

BACKGROUND

Cyanobacteria are the major prokaryotic primary producers occupying a range of aquatic habitats worldwide that differ in levels of salinity, making them a group of interest to study one of the major unresolved conundrums in aquatic microbiology which is what distinguishes a marine microbe from a freshwater one? We address this question using ecogenomics of a group of picocyanobacteria (cluster 5) that have recently evolved to inhabit geographically disparate salinity niches. Our analysis is made possible by the sequencing of 58 new genomes from freshwater representatives of this group that are presented here, representing a 6-fold increase in the available genomic data.

RESULTS

Overall, freshwater strains had larger genomes (≈2.9 Mb) and %GC content (≈64%) compared to brackish (2.69 Mb and 64%) and marine (2.5 Mb and 58.5%) isolates. Genomic novelties/differences across the salinity divide highlighted acidic proteomes and specific salt adaptation pathways in marine isolates (e.g., osmolytes/compatible solutes - glycine betaine/ggp/gpg/gmg clusters and glycerolipids glpK/glpA), while freshwater strains possessed distinct ion/potassium channels, permeases (aquaporin Z), fatty acid desaturases, and more neutral/basic proteomes. Sulfur, nitrogen, phosphorus, carbon (photosynthesis), or stress tolerance metabolism while showing distinct genomic footprints between habitats, e.g., different types of transporters, did not obviously translate into major functionality differences between environments. Brackish microbes show a mixture of marine (salt adaptation pathways) and freshwater features, highlighting their transitional nature.

CONCLUSIONS

The plethora of freshwater isolates provided here, in terms of trophic status preference and genetic diversity, exemplifies their ability to colonize ecologically diverse waters across the globe. Moreover, a trend towards larger and more flexible/adaptive genomes in freshwater picocyanobacteria may hint at a wider number of ecological niches in this environment compared to the relatively homogeneous marine system.

TrmB Family Transcription Factor as a Thiol-Based Regulator of Oxidative Stress Response

Oxidative stress causes cellular damage, including DNA mutations, protein dysfunction, and loss of membrane integrity. Here, we discovered that a TrmB (transcription regulator of mal operon) family protein (Pfam PF01978) composed of a single winged-helix DNA binding domain (InterPro IPR002831) can function as thiol-based transcriptional regulator of oxidative stress response. Using the archaeon Haloferax volcanii as a model system, we demonstrate that the TrmB-like OxsR is important for recovery of cells from hypochlorite stress. OxsR is shown to bind specific regions of genomic DNA, particularly during hypochlorite stress. OxsR-bound intergenic regions were found proximal to oxidative stress operons, including genes associated with thiol relay and low molecular weight thiol biosynthesis. Further analysis of a subset of these sites revealed OxsR to function during hypochlorite stress as a transcriptional activator and repressor. OxsR was shown to require a conserved cysteine (C24) for function and to use a CG-rich motif upstream of conserved BRE/TATA box promoter elements for transcriptional activation. Protein modeling suggested the C24 is located at a homodimer interface formed by antiparallel α helices, and that oxidation of this cysteine would result in the formation of an intersubunit disulfide bond. This covalent linkage may promote stabilization of an OxsR homodimer with the enhanced DNA binding properties observed in the presence of hypochlorite stress. The phylogenetic distribution TrmB family proteins, like OxsR, that have a single winged-helix DNA binding domain and conserved cysteine residue suggests this type of redox signaling mechanism is widespread in Archaea.

PCycDB: a comprehensive and accurate database for fast analysis of phosphorus cycling genes

BACKGROUND

Phosphorus (P) is one of the most essential macronutrients on the planet, and microorganisms (including bacteria and archaea) play a key role in P cycling in all living things and ecosystems. However, our comprehensive understanding of key P cycling genes (PCGs) and microorganisms (PCMs) as well as their ecological functions remains elusive even with the rapid advancement of metagenome sequencing technologies. One of major challenges is a lack of a comprehensive and accurately annotated P cycling functional gene database.

RESULTS

In this study, we constructed a well-curated P cycling database (PCycDB) covering 139 gene families and 10 P metabolic processes, including several previously ignored PCGs such as pafA encoding phosphate-insensitive phosphatase, ptxABCD (phosphite-related genes), and novel aepXVWPS genes for 2-aminoethylphosphonate transporters. We achieved an annotation accuracy, positive predictive value (PPV), sensitivity, specificity, and negative predictive value (NPV) of 99.8%, 96.1%, 99.9%, 99.8%, and 99.9%, respectively, for simulated gene datasets. Compared to other orthology databases, PCycDB is more accurate, more comprehensive, and faster to profile the PCGs. We used PCycDB to analyze P cycling microbial communities from representative natural and engineered environments and showed that PCycDB could apply to different environments.

CONCLUSIONS

We demonstrate that PCycDB is a powerful tool for advancing our understanding of microbially driven P cycling in the environment with high coverage, high accuracy, and rapid analysis of metagenome sequencing data. The PCycDB is available at https://github.com/ZengJiaxiong/Phosphorus-cycling-database . Video Abstract.

Microbial diversity in extreme environments

A wide array of microorganisms, including many novel, phylogenetically deeply rooted taxa, survive and thrive in extreme environments. These unique and reduced-complexity ecosystems offer a tremendous opportunity for studying the structure, function and evolution of natural microbial communities. Marker gene surveys have resolved patterns and ecological drivers of these extremophile assemblages, revealing a vast uncultured microbial diversity and the often predominance of archaea in the most extreme conditions. New omics studies have uncovered linkages between community function and environmental variables, and have enabled discovery and genomic characterization of major new lineages that substantially expand microbial diversity and change the structure of the tree of life. These efforts have significantly advanced our understanding of the diversity, ecology and evolution of microorganisms populating Earth's extreme environments, and have facilitated the exploration of microbiota and processes in more complex ecosystems.

Influence of Non-canonical DNA Bases on the Genomic Diversity of Tevenvirinae

The Tevenvirinae viruses are some of the most common viruses on Earth. Representatives of this subfamily have long been used in the molecular biology studies as model organisms – since the emergence of the discipline. Tevenvirinae are promising agents for phage therapy in animals and humans, since their representatives have only lytic life cycle and many of their host bacteria are pathogens. As confirmed experimentally, some Tevenvirinae have non-canonical DNA bases. Non-canonical bases can play an essential role in the diversification of closely related viruses. The article performs a comparative and evolutionary analysis of Tevenvirinae genomes and components of Tevenvirinae genomes. A comparative analysis of these genomes and the genes associated with the synthesis of non-canonical bases allows us to conclude that non-canonical bases have a major influence on the divergence of Tevenvirinae viruses within the same habitats. Supposedly, Tevenvirinae developed a strategy for changing HGT frequency in individual populations, which was based on the accumulation of proteins for the synthesis of non-canonical bases and proteins that used those bases as substrates. Owing to this strategy, ancestors of Tevenvirinae with the highest frequency of HGT acquired genes that allowed them to exist in a certain niche, and ancestors with the lowest HGT frequency preserved the most adaptive of those genes. Given the origin and characteristics of genes associated with the synthesis of non-canonical bases in Tevenvirinae, one can assume that other phages may have similar strategies. The article demonstrates the dependence of genomic diversity of closely related Tevenvirinae on non-canonical bases.

Insights into synthesis and function of KsgA/Dim1-dependent rRNA modifications in archaea

Ribosomes are intricate molecular machines ensuring proper protein synthesis in every cell. Ribosome biogenesis is a complex process which has been intensively analyzed in bacteria and eukaryotes. In contrast, our understanding of the in vivo archaeal ribosome biogenesis pathway remains less characterized. Here, we have analyzed the in vivo role of the almost universally conserved ribosomal RNA dimethyltransferase KsgA/Dim1 homolog in archaea. Our study reveals that KsgA/Dim1-dependent 16S rRNA dimethylation is dispensable for the cellular growth of phylogenetically distant archaea. However, proteomics and functional analyses suggest that archaeal KsgA/Dim1 and its rRNA modification activity (i) influence the expression of a subset of proteins and (ii) contribute to archaeal cellular fitness and adaptation. In addition, our study reveals an unexpected KsgA/Dim1-dependent variability of rRNA modifications within the archaeal phylum. Combining structure-based functional studies across evolutionary divergent organisms, we provide evidence on how rRNA structure sequence variability (re-)shapes the KsgA/Dim1-dependent rRNA modification status. Finally, our results suggest an uncoupling between the KsgA/Dim1-dependent rRNA modification completion and its release from the nascent small ribosomal subunit. Collectively, our study provides additional understandings into principles of molecular functional adaptation, and further evolutionary and mechanistic insights into an almost universally conserved step of ribosome synthesis.

Assessment of assumptions underlying models of prokaryotic pangenome evolution

BACKGROUND

The genomes of bacteria and archaea evolve by extensive loss and gain of genes which, for any group of related prokaryotic genomes, result in the formation of a pangenome with the universal, asymmetrical U-shaped distribution of gene commonality. However, the evolutionary factors that define the specific shape of this distribution are not thoroughly understood.

RESULTS

We investigate the fit of simple models of genome evolution to the empirically observed gene commonality distributions and genome intersections for 33 groups of closely related bacterial genomes. A model with an infinite external gene pool available for gene acquisition and constant genome size (IGP-CGS model), and two gene turnover rates, one for slow- and the other one for fast-evolving genes, allows two approaches to estimate the parameters for gene content dynamics. One is by fitting the model prediction to the distribution of the number of genes shared by precisely k genomes (gene commonality distribution) and another by analyzing the distribution of the number of genes common for k genome sets (k-cores). Both approaches produce a comparable overall quality of fit, although the former significantly overestimates the number of the universally conserved genes, while the latter overestimates the number of singletons. We further explore the effect of dropping each of the assumptions of the IGP-CGS model on the fit to the gene commonality distributions and show that models with either a finite gene pool or unequal rates of gene loss and gain (greater gene loss rate) eliminate the overestimate of the number of singletons or the core genome size.

CONCLUSIONS

We examine the assumptions that are usually adopted for modeling the evolution of the U-shaped gene commonality distributions in prokaryote genomes, namely, those of infinitely many genes and constant genome size. The combined analysis of genome intersections and gene commonality suggests that at least one of these assumptions is invalid. The violation of both these assumptions reflects the limited ability of prokaryotes to gain new genes. This limitation seems to stem, at least partly, from the horizontal gene transfer barrier, i.e., the cost of accommodation of foreign genes by prokaryotes. Further development of models taking into account the complexity of microbial evolution is necessary for an improved understanding of the evolution of prokaryotes.

COG database update: focus on microbial diversity, model organisms, and widespread pathogens

The Clusters of Orthologous Genes (COG) database, also referred to as the Clusters of Orthologous Groups of proteins, was created in 1997 and went through several rounds of updates, most recently, in 2014. The current update, available at https://www.ncbi.nlm.nih.gov/research/COG, substantially expands the scope of the database to include complete genomes of 1187 bacteria and 122 archaea, typically, with a single genome per genus. In addition, the current version of the COGs includes the following new features: (i) the recently deprecated NCBI's gene index (gi) numbers for the encoded proteins are replaced with stable RefSeq or GenBank\ENA\DDBJ coding sequence (CDS) accession numbers; (ii) COG annotations are updated for >200 newly characterized protein families with corresponding references and PDB links, where available; (iii) lists of COGs grouped by pathways and functional systems are added; (iv) 266 new COGs for proteins involved in CRISPR-Cas immunity, sporulation in Firmicutes and photosynthesis in cyanobacteria are included; and (v) the database is made available as a web page, in addition to FTP. The current release includes 4877 COGs. Future plans include further expansion of the COG collection by adding archaeal COGs (arCOGs), splitting the COGs containing multiple paralogs, and continued refinement of COG annotations.

Newly discovered Asgard archaea Hermodarchaeota potentially degrade alkanes and aromatics via alkyl/benzyl-succinate synthase and benzoyl-CoA pathway

Asgard archaea are widely distributed in anaerobic environments. Previous studies revealed the potential capability of Asgard archaea to utilize various organic substrates including proteins, carbohydrates, fatty acids, amino acids and hydrocarbons, suggesting that Asgard archaea play an important role in sediment carbon cycling. Here, we describe a previously unrecognized archaeal phylum, Hermodarchaeota, affiliated with the Asgard superphylum. The genomes of these archaea were recovered from metagenomes generated from mangrove sediments, and were found to encode alkyl/benzyl-succinate synthases and their activating enzymes that are similar to those identified in alkane-degrading sulfate-reducing bacteria. Hermodarchaeota also encode enzymes potentially involved in alkyl-coenzyme A and benzoyl-coenzyme A oxidation, the Wood–Ljungdahl pathway and nitrate reduction. These results indicate that members of this phylum have the potential to strictly anaerobically degrade alkanes and aromatic compounds, coupling the reduction of nitrate. By screening Sequence Read Archive, additional genes encoding 16S rRNA and alkyl/benzyl-succinate synthases analogous to those in Hermodarchaeota were identified in metagenomic datasets from a wide range of marine and freshwater sediments. These findings suggest that Asgard archaea capable of degrading alkanes and aromatics via formation of alkyl/benzyl-substituted succinates are ubiquitous in sediments.

A Novel Family of Acinetobacter Mega-Plasmids Are Disseminating Multi-Drug Resistance Across the Globe While Acquiring Location-Specific Accessory Genes

Acinetobacter species are emerging as major nosocomial pathogens, aided by their ability to acquire resistance to all classes of antibiotics. A key factor leading to their multi-drug resistance phenotypes is the acquisition of a wide variety of mobile genetic elements, particularly large conjugative plasmids. Here, we characterize a family of 21 multi-drug resistance mega-plasmids in 11 different Acinetobacter species isolated from various locations across the globe. The plasmid family exhibits a highly dynamic and diverse accessory genome, including 221 antibiotic resistance genes (ARGs) that confer resistance to 13 classes of antibiotics. We show that plasmids isolated within the same geographic region are often evolutionarily divergent members of this family based on their core-genome, yet they exhibit a more similar accessory genome. Individual plasmids, therefore, can disseminate to different locations around the globe, where they then appear to acquire diverse sets of accessory genes from their local surroundings. Further, we show that plasmids from several geographic regions were enriched with location-specific functional traits. Together, our findings show that these mega-plasmids can transmit across species boundaries, have the capacity for global dissemination, can accumulate a diverse suite of location-specific accessory genes, and can confer multi-drug resistance phenotypes of significant concern for human health. We therefore highlight this previously undescribed plasmid family as a serious threat to healthcare systems worldwide. These findings also add to the growing concern that mega-plasmids are key disseminators of antibiotic resistance and require global surveillance.

Biochemical and Taxonomic Characterization of Novel Haloarchaeal Strains and Purification of the Recombinant Halotolerant α-Amylase Discovered in the Isolate

Haloarchaea are salt-loving archaea and potential source of industrially relevant halotolerant enzymes. In the present study, three reddish-pink, extremely halophilic archaeal strains, namely wsp1 (wsp-water sample Pondicherry), wsp3, and wsp4, were isolated from the Indian Solar saltern. The phylogenetic analysis based on 16S rRNA gene sequences suggests that both wsp3 and wsp4 strains belong to Halogeometricum borinquense while wsp1 is closely related to Haloferax volcanii species. The comparative genomics revealed an open pangenome for both genera investigated here. Whole-genome sequence analysis revealed that these isolates have multiple copies of industrially/biotechnologically important unique genes and enzymes. Among these unique enzymes, for recombinant expression and purification, we selected four putative α-amylases identified in these three isolates. We successfully purified functional halotolerant recombinant Amy2, from wsp1 using pelB signal sequence-based secretion strategy using Escherichia coli as an expression host. This method may prove useful to produce functional haloarchaeal secretory recombinant proteins suitable for commercial or research applications. Biochemical analysis of Amy2 suggests the halotolerant nature of the enzyme having maximum enzymatic activity observed at 1 M NaCl. We also report the isolation and characterization of carotenoids purified from these isolates. This study highlights the presence of several industrially important enzymes in the haloarchaeal strains which may potentially have improved features like stability and salt tolerance suitable for industrial applications.

Microbial genome analysis: the COG approach

For the past 20 years, the Clusters of Orthologous Genes (COG) database had been a popular tool for microbial genome annotation and comparative genomics. Initially created for the purpose of evolutionary classification of protein families, the COG have been used, apart from straightforward functional annotation of sequenced genomes, for such tasks as (i) unification of genome annotation in groups of related organisms; (ii) identification of missing and/or undetected genes in complete microbial genomes; (iii) analysis of genomic neighborhoods, in many cases allowing prediction of novel functional systems; (iv) analysis of metabolic pathways and prediction of alternative forms of enzymes; (v) comparison of organisms by COG functional categories; and (vi) prioritization of targets for structural and functional characterization. Here we review the principles of the COG approach and discuss its key advantages and drawbacks in microbial genome analysis.

FttA is a CPSF73 homologue that terminates transcription in Archaea

Regulated gene expression is largely achieved by controlling the activities of essential, multisubunit RNA polymerase transcription elongation complexes (TECs). The extreme stability required of TECs to processively transcribe large genomic regions necessitates robust mechanisms to terminate transcription. Efficient transcription termination is particularly critical for gene-dense bacterial and archaeal genomes1-3 in which continued transcription would necessarily transcribe immediately adjacent genes and result in conflicts between the transcription and replication apparatuses4-6; the coupling of transcription and translation7,8 would permit the loading of ribosomes onto aberrant transcripts. Only select sequences or transcription termination factors can disrupt the otherwise extremely stable TEC and we demonstrate that one of the last universally conserved archaeal proteins with unknown biological function is the Factor that terminates transcription in Archaea (FttA). FttA resolves the dichotomy of a prokaryotic gene structure (operons and polarity) and eukaryotic molecular homology (general transcription apparatus) that is observed in Archaea. This missing link between prokaryotic and eukaryotic transcription regulation provides the most parsimonious link to the evolution of the processing activities involved in RNA 3'-end formation in Eukarya.

Comparative Genomics Reveals Ecological and Evolutionary Insights into Sponge-Associated Thaumarchaeota

Thaumarchaeota are frequently reported to associate with marine sponges (phylum Porifera); however, little is known about the features that distinguish them from their free-living thaumarchaeal counterparts. In this study, thaumarchaeal metagenome-assembled genomes (MAGs) were reconstructed from metagenomic data sets derived from the marine sponges Hexadella detritifera, Hexadella cf. detritifera, and Stylissa flabelliformis Phylogenetic and taxonomic analyses revealed that the three thaumarchaeal MAGs represent two new species within the genus Nitrosopumilus and one novel genus, for which we propose the names "Candidatus UNitrosopumilus hexadellus," "Candidatus UNitrosopumilus detritiferus," and "Candidatus UCenporiarchaeum stylissum" (the U superscript indicates that the taxon is uncultured). Comparison of these genomes to data from the Sponge Earth Microbiome Project revealed that "Ca UCenporiarchaeum stylissum" has been exclusively detected in sponges and can hence be classified as a specialist, while "Ca UNitrosopumilus detritiferus" and "Ca UNitrosopumilus hexadellus" are also detected outside the sponge holobiont and likely lead a generalist lifestyle. Comparison of the sponge-associated MAGs to genomes of free-living Thaumarchaeota revealed signatures that indicate functional features of a sponge-associated lifestyle, and these features were related to nutrient transport and metabolism, restriction-modification, defense mechanisms, and host interactions. Each species exhibited distinct functional traits, suggesting that they have reached different stages of evolutionary adaptation and/or occupy distinct ecological niches within their sponge hosts. Our study therefore offers new evolutionary and ecological insights into the symbiosis between sponges and their thaumarchaeal symbionts.IMPORTANCE Sponges represent ecologically important models to understand the evolution of symbiotic interactions of metazoans with microbial symbionts. Thaumarchaeota are commonly found in sponges, but their potential adaptations to a host-associated lifestyle are largely unknown. Here, we present three novel sponge-associated thaumarchaeal species and compare their genomic and predicted functional features with those of closely related free-living counterparts. We found different degrees of specialization of these thaumarchaeal species to the sponge environment that is reflected in their host distribution and their predicted molecular and metabolic properties. Our results indicate that Thaumarchaeota may have reached different stages of evolutionary adaptation in their symbiosis with sponges.

Phylogenomics suggests oxygen availability as a driving force in Thaumarchaeota evolution

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread in marine and terrestrial habitats, playing a major role in the global nitrogen cycle. However, their evolutionary history remains unexplored, which limits our understanding of their adaptation mechanisms. Here, our comprehensive phylogenomic tree of Thaumarchaeota supports three sequential events: origin of AOA from terrestrial non-AOA ancestors, colonization of the shallow ocean, and expansion to the deep ocean. Careful molecular dating suggests that these events coincided with the Great Oxygenation Event around 2300 million years ago (Mya), and oxygenation of the shallow and deep ocean around 800 and 635-560 Mya, respectively. The first transition was likely enabled by the gain of an aerobic pathway for energy production by ammonia oxidation and biosynthetic pathways for cobalamin and biotin that act as cofactors in aerobic metabolism. The first transition was also accompanied by the loss of dissimilatory nitrate and sulfate reduction, loss of oxygen-sensitive pyruvate oxidoreductase, which reduces pyruvate to acetyl-CoA, and loss of the Wood-Ljungdahl pathway for anaerobic carbon fixation. The second transition involved gain of a K⁺ transporter and of the biosynthetic pathway for ectoine, which may function as an osmoprotectant. The third transition was accompanied by the loss of the uvr system for repairing ultraviolet light-induced DNA lesions. We conclude that oxygen availability drove the terrestrial origin of AOA and their expansion to the photic and dark oceans, and that the stressors encountered during these events were partially overcome by gene acquisitions from Euryarchaeota and Bacteria, among other sources.

Towards functional characterization of archaeal genomic dark matter

A substantial fraction of archaeal genes, from ∼30% to as much as 80%, encode ‘hypothetical' proteins or genomic ‘dark matter'. Archaeal genomes typically contain a higher fraction of dark matter compared with bacterial genomes, primarily, because isolation and cultivation of most archaea in the laboratory, and accordingly, experimental characterization of archaeal genes, are difficult. In the present study, we present quantitative characteristics of the archaeal genomic dark matter and discuss comparative genomic approaches for functional prediction for ‘hypothetical' proteins. We propose a list of top priority candidates for experimental characterization with a broad distribution among archaea and those that are characteristic of poorly studied major archaeal groups such as Thaumarchaea, DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota and Nanohaloarchaeota) and Asgard.

Oxygen-mediated growth enhancement of an obligate anaerobic archaeon Thermococcus onnurineus NA1

Thermococcus onnurineus NA1, an obligate anaerobic hyperthermophilic archaeon, showed variable oxygen (O₂) sensitivity depending on the types of substrate employed as an energy source. Unexpectedly, the culture with yeast extract as a sole energy source showed enhanced growth by 2-fold in the presence of O₂. Genome-wide transcriptome analysis revealed the upregulation of several antioxidant-related genes encoding thioredoxin peroxidase (TON_0862), rubrerythrin (TON_0864), rubrerythrin-related protein (TON_0873), NAD(P)H rubredoxin oxidoreductase (TON_0865), or thioredoxin reductase (TON_1603), which can couple the detoxification of reactive oxygen species with the regeneration of NAD(P)⁺ from NAD(P)H. We present a plausible mechanism by which O₂ serves to maintain the intracellular redox balance. This study demonstrates an unusual strategy of an obligate anaerobe underlying O₂-mediated growth enhancement despite not having heme-based or cytochrome-type proteins.

Early Response of Sulfolobus acidocaldarius to Nutrient Limitation

In natural environments microorganisms encounter extreme changes in temperature, pH, osmolarities and nutrient availability. The stress response of many bacterial species has been described in detail, however, knowledge in Archaea is limited. Here, we describe the cellular response triggered by nutrient limitation in the thermoacidophilic crenarchaeon Sulfolobus acidocaldarius. We measured changes in gene transcription and protein abundance upon nutrient depletion up to 4 h after initiation of nutrient depletion. Transcript levels of 1118 of 2223 protein coding genes and abundance of approximately 500 proteins with functions in almost all cellular processes were affected by nutrient depletion. Our study reveals a significant rerouting of the metabolism with respect to degradation of internal as well as extracellular-bound organic carbon and degradation of proteins. Moreover, changes in membrane lipid composition were observed in order to access alternative sources of energy and to maintain pH homeostasis. At transcript level, the cellular response to nutrient depletion in S. acidocaldarius seems to be controlled by the general transcription factors TFB2 and TFEβ. In addition, ribosome biogenesis is reduced, while an increased protein degradation is accompanied with a loss of protein quality control. This study provides first insights into the early cellular response of Sulfolobus to organic carbon and organic nitrogen depletion.

Genome rearrangements and selection in multi-chromosome bacteria Burkholderia spp

BACKGROUND

The genus Burkholderia consists of species that occupy remarkably diverse ecological niches. Its best known members are important pathogens, B. mallei and B. pseudomallei, which cause glanders and melioidosis, respectively. Burkholderia genomes are unusual due to their multichromosomal organization, generally comprised of 2-3 chromosomes.

RESULTS

We performed integrated genomic analysis of 127 Burkholderia strains. The pan-genome is open with the saturation to be reached between 86,000 and 88,000 genes. The reconstructed rearrangements indicate a strong avoidance of intra-replichore inversions that is likely caused by selection against the transfer of large groups of genes between the leading and the lagging strands. Translocated genes also tend to retain their position in the leading or the lagging strand, and this selection is stronger for large syntenies. Integrated reconstruction of chromosome rearrangements in the context of strains phylogeny reveals parallel rearrangements that may indicate inversion-based phase variation and integration of new genomic islands. In particular, we detected parallel inversions in the second chromosomes of B. pseudomallei with breakpoints formed by genes encoding membrane components of multidrug resistance complex, that may be linked to a phase variation mechanism. Two genomic islands, spreading horizontally between chromosomes, were detected in the B. cepacia group.

CONCLUSIONS

This study demonstrates the power of integrated analysis of pan-genomes, chromosome rearrangements, and selection regimes. Non-random inversion patterns indicate selective pressure, inversions are particularly frequent in a recent pathogen B. mallei, and, together with periods of positive selection at other branches, may indicate adaptation to new niches. One such adaptation could be a possible phase variation mechanism in B. pseudomallei.

Species Widely Distributed in Halophilic Archaea Exhibit Opsin-Mediated Inhibition of Bacterioruberin Biosynthesis

Halophilic Archaea are a distinctive pink color due to a carotenoid pigment called bacterioruberin. To sense or utilize light, many halophilic Archaea also produce rhodopsins, complexes of opsin proteins with a retinal prosthetic group. Both bacterioruberin and retinal are synthesized from isoprenoid precursors, with lycopene as the last shared intermediate. We previously described a regulatory mechanism by which Halobacterium salinarum bacterioopsin and Haloarcula vallismortis cruxopsin inhibit bacterioruberin synthesis catalyzed by lycopene elongase. In this work, we found that opsins in all three major Halobacteria clades inhibit bacterioruberin synthesis, suggesting that this regulatory mechanism existed in the common Halobacteria ancestor. Halophilic Archaea, which are generally heterotrophic and aerobic, likely evolved from an autotrophic, anaerobic methanogenic ancestor by acquiring many genes from Bacteria via lateral gene transfer. These bacterial "imports" include genes encoding opsins and lycopene elongases. To determine if opsins from Bacteria inhibit bacterioruberin synthesis, we tested bacterial opsins and found that an opsin from Curtobacterium, in the Actinobacteria phylum, inhibits bacterioruberin synthesis catalyzed by its own lycopene elongase, as well as that catalyzed by several archaeal enzymes. We also determined that the lycopene elongase from Halococcus salifodinae, a species from a family of Halobacteria lacking opsin homologs, retained the capacity to be inhibited by opsins. Together, our results indicate that opsin-mediated inhibition of bacterioruberin biosynthesis is a widely distributed mechanism found in both Archaea and Bacteria, possibly predating the divergence of the two domains. Further analysis may provide insight into the acquisition and evolution of the genes and their host species.IMPORTANCE All organisms use a variety of mechanisms to allocate limited resources to match their needs in their current environment. Here, we explore how halophilic microbes use a novel mechanism to allow efficient production of rhodopsin, a complex of an opsin protein and a retinal prosthetic group. We previously demonstrated that Halobacterium salinarum bacterioopsin directs available resources toward retinal by inhibiting synthesis of bacterioruberin, a molecule that shares precursors with retinal. In this work, we show that this mechanism can be carried out by proteins from halophilic Archaea that are not closely related to H. salinarum and those in at least one species of Bacteria Therefore, opsin-mediated inhibition of bacterioruberin synthesis may be a highly conserved, ancient regulatory mechanism.

Ecological and genomic features of two widespread freshwater picocyanobacteria

We present two genomes of widespread freshwater picocyanobacteria isolated by extinction dilution from a Spanish oligotrophic reservoir. Based on microscopy and genomic properties, both picocyanobacteria were tentatively designated Synechococcus lacustris Tous, formerly described as a metagenome assembled genome (MAG) from the same habitat, and Cyanobium usitatum Tous, described here for the first time. Both strains were purified in unicyanobacterial cultures, and their genomes were sequenced. They are broadly distributed in freshwater systems; the first seems to be a specialist on temperate reservoirs (Tous, Amadorio, Dexter, Lake Lanier, Sparkling), and the second appears to also be abundant in cold environments including ice-covered lakes such as Lake Baikal, Lake Erie or the brackish Baltic Sea. Having complete genomes provided access to the flexible genome that does not assemble in MAGs. We found several genomic islands in both genomes, within which there were genes for nitrogen acquisition, transporters for a wide set of compounds and biosynthesis of phycobilisomes in both strains. Some of these regions of low coverage in metagenomes also included antimicrobial compounds, transposases and phage defence systems, including a novel type III CRISPR-Cas phage defence system that was only detected in Synechococcus lacustris Tous.

Dissecting the Repertoire of DNA-Binding Transcription Factors of the Archaeon Pyrococcus furiosus DSM 3638

In recent years, there has been a large increase in the amount of experimental evidence for diverse archaeal organisms, and these findings allow for a comprehensive analysis of archaeal genetic organization. However, studies about regulatory mechanisms in this cellular domain are still limited. In this context, we identified a repertoire of 86 DNA-binding transcription factors (TFs) in the archaeon Pyrococcus furiosus DSM 3638, that are clustered into 32 evolutionary families. In structural terms, 45% of these proteins are composed of one structural domain, 41% have two domains, and 14% have three structural domains. The most abundant DNA-binding domain corresponds to the winged helix-turn-helix domain; with few alternative DNA-binding domains. We also identified seven regulons, which represent 13.5% (279 genes) of the total genes in this archaeon. These analyses increase our knowledge about gene regulation in P. furiosus DSM 3638 and provide additional clues for comprehensive modeling of transcriptional regulatory networks in the Archaea cellular domain.

Copy number variation is associated with gene expression change in archaea

Genomic instability, although frequently deleterious, is also an important mechanism for microbial adaptation to environmental change. Although widely studied in bacteria, in archaea the effect of genomic instability on organism phenotypes and fitness remains unclear. Here we use DNA segmentation methods to detect and quantify genome-wide copy number variation (CNV) in large compendia of high-throughput datasets in a model archaeal species, Halobacterium salinarum. CNV hotspots were identified throughout the genome. Some hotspots were strongly associated with changes in gene expression, suggesting a mechanism for phenotypic innovation. In contrast, CNV hotspots in other genomic loci left expression unchanged, suggesting buffering of certain phenotypes. The correspondence of CNVs with gene expression was validated with strain- and condition-matched transcriptomics and DNA quantification experiments at specific loci. Significant correlation of CNV hotspot locations with the positions of known insertion sequence (IS) elements suggested a mechanism for generating genomic instability. Given the efficient recombination capabilities in H. salinarum despite stability at the single nucleotide level, these results suggest that genomic plasticity mediated by IS element activity can provide a source of phenotypic innovation in extreme environments.

GlpR Is a Direct Transcriptional Repressor of Fructose Metabolic Genes in Haloferax volcanii

DeoR-type helix-turn-helix (HTH) domain proteins are transcriptional regulators of sugar and nucleoside metabolism in diverse bacteria and also occur in select archaea. In the model archaeon Haloferax volcanii, previous work implicated GlpR, a DeoR-type transcriptional regulator, in the transcriptional repression of glpR and the gene encoding the fructose-specific phosphofructokinase (pfkB) during growth on glycerol. However, the global regulon governed by GlpR remained unclear. Here, we compared transcriptomes of wild-type and ΔglpR mutant strains grown on glycerol and glucose to detect significant transcript level differences for nearly 50 new genes regulated by GlpR. By coupling computational prediction of GlpR binding sequences with in vivo and in vitro DNA binding experiments, we determined that GlpR directly controls genes encoding enzymes involved in fructose degradation, including fructose bisphosphate aldolase, a central control point in glycolysis. GlpR also directly controls other transcription factors. In contrast, other metabolic pathways appear to be under the indirect influence of GlpR. In vitro experiments demonstrated that GlpR purifies to function as a tetramer that binds the effector molecule fructose-1-phosphate (F1P). These results suggest that H. volcanii GlpR functions as a direct negative regulator of fructose degradation during growth on carbon sources other than fructose, such as glucose and glycerol, and that GlpR bears striking functional similarity to bacterial DeoR-type regulators.IMPORTANCE Many archaea are extremophiles, able to thrive in habitats of extreme salinity, pH and temperature. These biological properties are ideal for applications in biotechnology. However, limited knowledge of archaeal metabolism is a bottleneck that prevents the broad use of archaea as microbial factories for industrial products. Here, we characterize how sugar uptake and use are regulated in a species that lives in high salinity. We demonstrate that a key sugar regulatory protein in this archaeal species functions using molecular mechanisms conserved with distantly related bacterial species.

Evolution of Genome Architecture in Archaea: Spontaneous Generation of a New Chromosome in Haloferax volcanii

The common ancestry of archaea and eukaryotes is evident in their genome architecture. All eukaryotic and several archaeal genomes consist of multiple chromosomes, each replicated from multiple origins. Three scenarios have been proposed for the evolution of this genome architecture: 1) mutational diversification of a multi-copy chromosome; 2) capture of a new chromosome by horizontal transfer; 3) acquisition of new origins and splitting into two replication-competent chromosomes. We report an example of the third scenario: the multi-origin chromosome of the archaeon Haloferax volcanii has split into two elements via homologous recombination. The newly generated elements are bona fide chromosomes, because each bears "chromosomal" replication origins, rRNA loci, and essential genes. The new chromosomes were stable during routine growth but additional genetic manipulation, which involves selective bottlenecks, provoked further rearrangements. To the best of our knowledge, rearrangement of a naturally evolved prokaryotic genome to generate two new chromosomes has not been described previously.

Phyletic Distribution and Lineage-Specific Domain Architectures of Archaeal Two-Component Signal Transduction Systems

The two-component signal transduction (TCS) machinery is a key mechanism of sensing environmental changes in the prokaryotic world. TCS systems have been characterized thoroughly in bacteria but to a much lesser extent in archaea. Here, we provide an updated census of more than 2,000 histidine kinases and response regulators encoded in 218 complete archaeal genomes, as well as unfinished genomes available from metagenomic data. We describe the domain architectures of the archaeal TCS components, including several novel output domains, and discuss the evolution of the archaeal TCS machinery. The distribution of TCS systems in archaea is strongly biased, with high levels of abundance in haloarchaea and thaumarchaea but none detected in the sequenced genomes from the phyla Crenarchaeota, Nanoarchaeota, and Korarchaeota The archaeal sensor histidine kinases are generally similar to their well-studied bacterial counterparts but are often located in the cytoplasm and carry multiple PAS and/or GAF domains. In contrast, archaeal response regulators differ dramatically from the bacterial ones. Most archaeal genomes do not encode any of the major classes of bacterial response regulators, such as the DNA-binding transcriptional regulators of the OmpR/PhoB, NarL/FixJ, NtrC, AgrA/LytR, and ActR/PrrA families and the response regulators with GGDEF and/or EAL output domains. Instead, archaea encode multiple copies of response regulators containing either the stand-alone receiver (REC) domain or combinations of REC with PAS and/or GAF domains. Therefore, the prevailing mechanism of archaeal TCS signaling appears to be via a variety of protein-protein interactions, rather than direct transcriptional regulation.IMPORTANCE Although the Archaea represent a separate domain of life, their signaling systems have been assumed to be closely similar to the bacterial ones. A study of the domain architectures of the archaeal two-component signal transduction (TCS) machinery revealed an overall similarity of archaeal and bacterial sensory modules but substantial differences in the signal output modules. The prevailing mechanism of archaeal TCS signaling appears to involve various protein-protein interactions rather than direct transcription regulation. The complete list of histidine kinases and response regulators encoded in the analyzed archaeal genomes is available online at http://www.ncbi.nlm.nih.gov/Complete_Genomes/TCSarchaea.html.

Proteomic Analysis of Methanonatronarchaeum thermophilum AMET1, a Representative of a Putative New Class of Euryarchaeota, "Methanonatronarchaeia"

The recently discovered Methanonatronarchaeia are extremely halophilic and moderately thermophilic methyl-reducing methanogens representing a novel class-level lineage in the phylum Euryarchaeota related to the class Halobacteria. Here we present a detailed analysis of 1D-nano liquid chromatography-electrospray ionization tandem mass spectrometry data obtained for "Methanonatronarchaeum thermophilum" AMET1 grown in different physiological conditions, including variation of the growth temperature and substrates. Analysis of these data allows us to refine the current understanding of the key biosynthetic pathways of this triple extremophilic methanogenic euryarchaeon and identify proteins that are likely to be involved in its response to growth condition changes.

Genetic Competence Drives Genome Diversity in Bacillus subtilis

Prokaryote genomes are the result of a dynamic flux of genes, with increases achieved via horizontal gene transfer and reductions occurring through gene loss. The ecological and selective forces that drive this genomic flexibility vary across species. Bacillus subtilis is a naturally competent bacterium that occupies various environments, including plant-associated, soil, and marine niches, and the gut of both invertebrates and vertebrates. Here, we quantify the genomic diversity of B. subtilis and infer the genome dynamics that explain the high genetic and phenotypic diversity observed. Phylogenomic and comparative genomic analyses of 42 B. subtilis genomes uncover a remarkable genome diversity that translates into a core genome of 1,659 genes and an asymptotic pangenome growth rate of 57 new genes per new genome added. This diversity is due to a large proportion of low-frequency genes that are acquired from closely related species. We find no gene-loss bias among wild isolates, which explains why the cloud genome, 43% of the species pangenome, represents only a small proportion of each genome. We show that B. subtilis can acquire xenologous copies of core genes that propagate laterally among strains within a niche. While not excluding the contributions of other mechanisms, our results strongly suggest a process of gene acquisition that is largely driven by competence, where the long-term maintenance of acquired genes depends on local and global fitness effects. This competence-driven genomic diversity provides B. subtilis with its generalist character, enabling it to occupy a wide range of ecological niches and cycle through them.

Functional divergence and comparative in-silico study of Cas4 proteins of DUF83 class

Clustered Regularly Interspaced Short Palindromic Repeats-CRISPR associated (CRISPR-Cas) systems present in genomes of bacteria and archaea have been the focus of many research studies recently. The Cas4 proteins of these systems are thought to be responsible for the adaptation step in the CRISPR mechanism. Cas4 proteins exhibit low sequence similarity among themselves and are currently classified into 2 main classes: DUF83 and DUF911. The characteristic features of Cas4 proteins belonging to DUF83 class have been elucidated by determining the structures of Cas4 protein from Sulfolobus solfataricus and Pyrobaculum calidifontis. Although, both Cas4 proteins characterized structurally are of same DUF83 class, these 2 proteins do exhibit significant biochemical and functional differences. The aim of the present study was to explore the structural and evolutionary features responsible for these differences. Our study predicts residues which might be responsible for such differences. Functional divergence analysis was used to predict sites exhibiting type I divergence, where certain amino acids are conserved in 1 clade whereas the same site is highly variable in the other clade. Our intra-molecular interaction analysis reinforces the influence of such divergence sites on the other functionally important amino acids. In general, this study identifies some of the divergence hotspots that could be the focus of future experimental studies for better understanding of Cas4 enzymatic activity in CRISPR mechanism.

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

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

Effect of the environment on horizontal gene transfer between bacteria and archaea

Background

Horizontal gene transfer, the transfer and incorporation of genetic material between different species of organisms, has an important but poorly quantified role in the adaptation of microbes to their environment. Previous work has shown that genome size and the number of horizontally transferred genes are strongly correlated. Here we consider how genome size confuses the quantification of horizontal gene transfer because the number of genes an organism accumulates over time depends on its evolutionary history and ecological context (e.g., the nutrient regime for which it is adapted).

Results

We investigated horizontal gene transfer between archaea and bacteria by first counting reciprocal BLAST hits among 448 bacterial and 57 archaeal genomes to find shared genes. Then we used the DarkHorse algorithm, a probability-based, lineage-weighted method (Podell & Gaasterland, 2007), to identify potential horizontally transferred genes among these shared genes. By removing the effect of genome size in the bacteria, we have identified bacteria with unusually large numbers of shared genes with archaea for their genome size. Interestingly, archaea and bacteria that live in anaerobic and/or high temperature conditions are more likely to share unusually large numbers of genes. However, high salt was not found to significantly affect the numbers of shared genes. Numbers of shared (genome size-corrected, reciprocal BLAST hits) and transferred genes (identified by DarkHorse) were strongly correlated. Thus archaea and bacteria that live in anaerobic and/or high temperature conditions are more likely to share horizontally transferred genes. These horizontally transferred genes are over-represented by genes involved in energy conversion as well as the transport and metabolism of inorganic ions and amino acids.

Conclusions

Anaerobic and thermophilic bacteria share unusually large numbers of genes with archaea. This is mainly due to horizontal gene transfer of genes from the archaea to the bacteria. In general, these transfers are from archaea that live in similar oxygen and temperature conditions as the bacteria that receive the genes. Potential hotspots of horizontal gene transfer between archaea and bacteria include hot springs, marine sediments, and oil wells. Cold spots for horizontal transfer included dilute, aerobic, mesophilic environments such as marine and freshwater surface waters.

Systematic Discovery of Archaeal Transcription Factor Functions in Regulatory Networks through Quantitative Phenotyping Analysis

To ensure survival in the face of stress, microorganisms employ inducible damage repair pathways regulated by extensive and complex gene networks. Many archaea, microorganisms of the third domain of life, persist under extremes of temperature, salinity, and pH and under other conditions. In order to understand the cause-effect relationships between the dynamic function of the stress network and ultimate physiological consequences, this study characterized the physiological role of nearly one-third of all regulatory proteins known as transcription factors (TFs) in an archaeal organism. Using a unique quantitative phenotyping approach, we discovered functions for many novel TFs and revealed important secondary functions for known TFs. Surprisingly, many TFs are required for resisting multiple stressors, suggesting cross-regulation of stress responses. Through extensive validation experiments, we map the physiological roles of these novel TFs in stress response back to their position in the regulatory network wiring. This study advances understanding of the mechanisms underlying how microorganisms resist extreme stress. Given the generality of the methods employed, we expect that this study will enable future studies on how regulatory networks adjust cellular physiology in a diversity of organisms.

Gene regulatory networks (GRNs) are critical for dynamic transcriptional responses to environmental stress. However, the mechanisms by which GRN regulation adjusts physiology to enable stress survival remain unclear. Here we investigate the functions of transcription factors (TFs) within the global GRN of the stress-tolerant archaeal microorganism Halobacterium salinarum. We measured growth phenotypes of a panel of TF deletion mutants in high temporal resolution under heat shock, oxidative stress, and low-salinity conditions. To quantitate the noncanonical functional forms of the growth trajectories observed for these mutants, we developed a novel modeling framework based on Gaussian process regression and functional analysis of variance (FANOVA). We employ unique statistical tests to determine the significance of differential growth relative to the growth of the control strain. This analysis recapitulated known TF functions, revealed novel functions, and identified surprising secondary functions for characterized TFs. Strikingly, we observed that the majority of the TFs studied were required for growth under multiple stress conditions, pinpointing regulatory connections between the conditions tested. Correlations between quantitative phenotype trajectories of mutants are predictive of TF-TF connections within the GRN. These phenotypes are strongly concordant with predictions from statistical GRN models inferred from gene expression data alone. With genome-wide and targeted data sets, we provide detailed functional validation of novel TFs required for extreme oxidative stress and heat shock survival. Together, results presented in this study suggest that many TFs function under multiple conditions, thereby revealing high interconnectivity within the GRN and identifying the specific TFs required for communication between networks responding to disparate stressors.

IMPORTANCE To ensure survival in the face of stress, microorganisms employ inducible damage repair pathways regulated by extensive and complex gene networks. Many archaea, microorganisms of the third domain of life, persist under extremes of temperature, salinity, and pH and under other conditions. In order to understand the cause-effect relationships between the dynamic function of the stress network and ultimate physiological consequences, this study characterized the physiological role of nearly one-third of all regulatory proteins known as transcription factors (TFs) in an archaeal organism. Using a unique quantitative phenotyping approach, we discovered functions for many novel TFs and revealed important secondary functions for known TFs. Surprisingly, many TFs are required for resisting multiple stressors, suggesting cross-regulation of stress responses. Through extensive validation experiments, we map the physiological roles of these novel TFs in stress response back to their position in the regulatory network wiring. This study advances understanding of the mechanisms underlying how microorganisms resist extreme stress. Given the generality of the methods employed, we expect that this study will enable future studies on how regulatory networks adjust cellular physiology in a diversity of organisms.

'ARMAN' archaea depend on association with euryarchaeal host in culture and in situ

Intriguing, yet uncultured ‘ARMAN’-like archaea are metabolically dependent on other members of the microbial community. It remains uncertain though which hosts they rely upon, and, because of the lack of complete genomes, to what extent. Here, we report the co-culturing of ARMAN-2-related organism, Mia14, with Cuniculiplasma divulgatum PM4 during the isolation of this strain from acidic streamer in Parys Mountain (Isle of Anglesey, UK). Mia14 is highly enriched in the binary culture (ca. 10% genomic reads) and its ungapped 0.95 Mbp genome points at severe voids in central metabolic pathways, indicating dependence on the host, C. divulgatum PM4. Analysis of C. divulgatum isolates from different sites and shotgun sequence data of Parys Mountain samples suggests an extensive genetic exchange between Mia14 and hosts in situ. Within the subset of organisms with high-quality genomic assemblies representing the ‘DPANN’ superphylum, the Mia14 lineage has had the largest gene flux, with dozens of genes gained that are implicated in the host interaction.

In the absence of complete genomes, the metabolic capabilities of uncultured ARMAN-like archaea have been uncertain. Here, Golyshina et al. apply an enrichment culture technique and find that the ungapped genome of the ARMAN-like archaeon Mia14 has lost key metabolic pathways, suggesting dependence on the host archaeon Cuniculiplasma divulgatum.

A Complex Endomembrane System in the Archaeon Ignicoccus hospitalis Tapped by Nanoarchaeum equitans

Based on serial sectioning, focused ion beam scanning electron microscopy (FIB/SEM), and electron tomography, we depict in detail the highly unusual anatomy of the marine hyperthermophilic crenarchaeon, Ignicoccus hospitalis. Our data support a complex and dynamic endomembrane system consisting of cytoplasmic protrusions, and with secretory function. Moreover, we reveal that the cytoplasm of the putative archaeal ectoparasite Nanoarchaeum equitans can get in direct contact with this endomembrane system, complementing and explaining recent proteomic, transcriptomic and metabolomic data on this inter-archaeal relationship. In addition, we identified a matrix of filamentous structures and/or tethers in the voluminous inter-membrane compartment (IMC) of I. hospitalis, which might be responsible for membrane dynamics. Overall, this unusual cellular compartmentalization, ultrastructure and dynamics in an archaeon that belongs to the recently proposed TACK superphylum prompts speculation that the eukaryotic endomembrane system might originate from Archaea.

Discovery of extremely halophilic, methyl-reducing euryarchaea provides insights into the evolutionary origin of methanogenesis

Methanogenic archaea are major players in the global carbon cycle and in the biotechnology of anaerobic digestion. The phylum Euryarchaeota includes diverse groups of methanogens that are interspersed with non-methanogenic lineages. So far, methanogens inhabiting hypersaline environments have been identified only within the order Methanosarcinales. We report the discovery of a deep phylogenetic lineage of extremophilic methanogens in hypersaline lakes and present analysis of two nearly complete genomes from this group. Within the phylum Euryarchaeota, these isolates form a separate, class-level lineage 'Methanonatronarchaeia' that is most closely related to the class Halobacteria. Similar to the Halobacteria, 'Methanonatronarchaeia' are extremely halophilic and do not accumulate organic osmoprotectants. The high intracellular concentration of potassium implies that 'Methanonatronarchaeia' employ the 'salt-in' osmoprotection strategy. These methanogens are heterotrophic methyl-reducers that use C₁-methylated compounds as electron acceptors and formate or hydrogen as electron donors. The genomes contain an incomplete and apparently inactivated set of genes encoding the upper branch of methyl group oxidation to CO₂ as well as membrane-bound heterodisulfide reductase and cytochromes. These features differentiate 'Methanonatronarchaeia' from all known methyl-reducing methanogens. The discovery of extremely halophilic, methyl-reducing methanogens related to haloarchaea provides insights into the origin of methanogenesis and shows that the strategies employed by methanogens to thrive in salt-saturating conditions are not limited to the classical methylotrophic pathway.

Reconstruction of the evolution of microbial defense systems

Background

Evolution of bacterial and archaeal genomes is a highly dynamic process that involves intensive loss of genes as well as gene gain via horizontal transfer, with a lesser contribution from gene duplication. The rates of these processes can be estimated by comparing genomes that are linked by an evolutionary tree. These estimated rates of genome dynamics events substantially differ for different functional classes of genes. The genes involved in defense against viruses and other invading DNA are among those that are gained and lost at the highest rates.

Results

We employed a stochastic birth-and-death model to obtain maximum likelihood estimates of the rates of gain and loss of defense genes in 35 groups of closely related bacterial genomes and one group of archaeal genomes. We find that on average, the defense genes experience 1.4 fold higher flux than the rest of microbial genes. This excessive flux of defense genes over the genomic mean is consistent across diverse microbial groups. The few exceptions include intracellular parasites with small, degraded genomes that possess few defense systems which are more stable than in other microbes. Generally, defense genes follow the previously established pattern of genome dynamics, with gene family loss being about 3 times more common than gain and an order of magnitude more common than expansion or contraction of gene families. Case by case analysis of the evolutionary dynamics of defense genes indicates frequent multiple events in the same locus and widespread involvement of mobile elements in the gain and loss of defense genes.

Conclusions

Evolution of microbial defense systems is highly dynamic but, notwithstanding the host-parasite arms race, generally follows the same trends that have been established for the rest of the genes. Apart from the paucity and the low flux of defense genes in parasitic bacteria with deteriorating genomes, there is no clear connection between the evolutionary regime of defense systems and microbial life style.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0942-y) contains supplementary material, which is available to authorized users.

Transcriptomic profiling and its implications for the H2 production of a non-methanogen deficient in the frhAGB-encoding hydrogenase

The F₄₂₀-reducing hydrogenase of methanogens functions in methanogenesis by providing reduced coenzyme F₄₂₀ (F₄₂₀H₂) as an electron donor. In non-methanogens, however, their physiological function has not been identified yet. In this study, we constructed an ΔfrhA mutant, whose frhA gene encoding the hydrogenase α subunit was deleted, in the non-methanogenic Thermococcus onnurineus NA1 as a model organism. There was no significant difference in the formate-dependent growth between the mutant and the wild-type strains. Interestingly, the mutation in the frhA gene affected the expression of genes involved in various cellular functions such as H₂ oxidation, chemotactic signal transduction, and carbon monoxide (CO) metabolism. Among these genes, the CO oxidation gene cluster, enabling CO-dependent growth and H₂ production, showed a 2.8- to 7.0-fold upregulation by microarray-based whole transcriptome expression profiling. The levels of proteins produced by this gene cluster were also significantly increased not only under the formate condition but also under the CO condition. In a controlled bioreactor, where 100% CO was continuously fed, the ΔfrhA mutant exhibited significant increases in cell growth (2.8-fold) and H₂ production (3.4-fold). These findings strongly imply that this hydrogenase is functional in non-methanogens and is related to various cellular metabolic processes through an unidentified mechanism. An understanding of the mechanism by which the frhA gene deletion affected the expression of other genes will provide insights that can be applied to the development of strategies for the enhancement of H₂ production using CO as a substrate.