Comparative genomics reveals conserved positioning of essential genomic clusters in highly rearranged Thermococcales chromosomes

Cooperation between two modes for DNA replication initiation in the archaeon Thermococcus barophilus

The mechanisms underpinning the replication of genomic DNA have recently been challenged in Archaea. Indeed, the lack of origin of replication has no deleterious effect on growth, suggesting that replication initiation relies on homologous recombination. Recombination-dependent replication (RDR) appears to be based on the recombinase RadA, which is of absolute requirement when no initiation origins are detected. The origin of this flexibility in the initiation of replication and the extent to which it is used in nature are yet to be understood. Here, we followed the process of DNA replication throughout the growth stages of Thermococcus barophilus. We combined deep sequencing and genetics to elucidate the dynamics of oriC utilization according to growth phases. We discovered that in T. barophilus, the use of oriC diminishes from the lag to the middle of the log phase, and subsequently increases gradually upon entering the stationary phase. Although oriC demonstrates no indispensability, RadA does exhibit essentiality. Notably, a knockdown mutant strain provides confirmation of the pivotal role of RadA in RDR for the first time. Thus, we demonstrate the existence of a tight combination between oriC utilization and homologous recombination to initiate DNA replication along the growth phases. Overall, this study demonstrates how diverse physiological states can influence the initiation of DNA replication, offering insights into how environmental sensing might impact this fundamental mechanism of life.

IMPORTANCE

Replication of DNA is highly important in all organisms. It initiates at a specific locus called ori, which serves as the binding site for scaffold proteins-either Cdc6 or DnaA-depending on the domain of life. However, recent studies have shown that the Archaea, Haloferax volcanii and Thermococcus kodakarensis could subsist without ori. Recombination-dependent replication (RDR), via the recombinase RadA, is the mechanism that uses homologous recombination to initiate DNA replication. The extent to which ori's use is necessary in natural growth remains to be characterized. In this study, using Thermococcus barophilus, we demonstrated that DNA replication initiation relies on both oriC and RDR throughout its physiological growth, each to varying degrees depending on the phase. Notably, a knockdown RadA mutant confirmed the prominent use of RDR during the log phase. Moreover, the study of ploidy in oriC and radA mutant strains showed that the number of chromosomes per cell is a critical proxy for ensuring proper growth and cell survival.

Membrane adaptation in the hyperthermophilic archaeon Pyrococcus furiosus relies upon a novel strategy involving glycerol monoalkyl glycerol tetraether lipids

Microbes preserve membrane functionality under fluctuating environmental conditions by modulating their membrane lipid composition. Although several studies have documented membrane adaptations in Archaea, the influence of most biotic and abiotic factors on archaeal lipid compositions remains underexplored. Here, we studied the influence of temperature, pH, salinity, the presence/absence of elemental sulfur, the carbon source, and the genetic background on the lipid core composition of the hyperthermophilic neutrophilic marine archaeon Pyrococcus furiosus. Every growth parameter tested affected the lipid core composition to some extent, the carbon source and the genetic background having the greatest influence. Surprisingly, P. furiosus appeared to only marginally rely on the two major responses implemented by Archaea, i.e., the regulation of the ratio of diether to tetraether lipids and that of the number of cyclopentane rings in tetraethers. Instead, this species increased the ratio of glycerol monoalkyl glycerol tetraethers (GMGT, aka. H-shaped tetraethers) to glycerol dialkyl glycerol tetrathers (GDGT) in response to decreasing temperature and pH and increasing salinity, thus providing for the first time evidence of adaptive functions for GMGT. Besides P. furiosus, numerous other species synthesize significant proportions of GMGT, which suggests that this unprecedented adaptive strategy might be common in Archaea. This article is protected by copyright. All rights reserved.

Archaeal tyrosine recombinases

The integration of mobile genetic elements into their host chromosome influences the immediate fate of cellular organisms and gradually shapes their evolution. Site-specific recombinases catalyzing this integration have been extensively characterized both in bacteria and eukarya. More recently, a number of reports provided the in-depth characterization of archaeal tyrosine recombinases and highlighted new particular features not observed in the other two domains. In addition to being active in extreme environments, archaeal integrases catalyze reactions beyond site-specific recombination. Some of these integrases can catalyze low-sequence specificity recombination reactions with the same outcome as homologous recombination events generating deep rearrangements of their host genome. A large proportion of archaeal integrases are termed suicidal due to the presence of a specific recombination target within their own gene. The paradoxical maintenance of integrases that disrupt their gene upon integration implies novel mechanisms for their evolution. In this review, we assess the diversity of the archaeal tyrosine recombinases using a phylogenomic analysis based on an exhaustive similarity network. We outline the biochemical, ecological and evolutionary properties of these enzymes in the context of the families we identified and emphasize similarities and differences between archaeal recombinases and their bacterial and eukaryal counterparts.

Extended Archaeal Histone-Based Chromatin Structure Regulates Global Gene Expression in Thermococcus kodakarensis

Histone proteins compact and organize DNA resulting in a dynamic chromatin architecture impacting DNA accessibility and ultimately gene expression. Eukaryotic chromatin landscapes are structured through histone protein variants, epigenetic marks, the activities of chromatin-remodeling complexes, and post-translational modification of histone proteins. In most Archaea, histone-based chromatin structure is dominated by the helical polymerization of histone proteins wrapping DNA into a repetitive and closely gyred configuration. The formation of the archaeal-histone chromatin-superhelix is a regulatory force of adaptive gene expression and is likely critical for regulation of gene expression in all histone-encoding Archaea. Single amino acid substitutions in archaeal histones that block formation of tightly packed chromatin structures have profound effects on cellular fitness, but the underlying gene expression changes resultant from an altered chromatin landscape have not been resolved. Using the model organism Thermococcus kodakarensis, we genetically alter the chromatin landscape and quantify the resultant changes in gene expression, including unanticipated and significant impacts on provirus transcription. Global transcriptome changes resultant from varying chromatin landscapes reveal the regulatory importance of higher-order histone-based chromatin architectures in regulating archaeal gene expression.

Metagenomics survey unravels diversity of biogas microbiomes with potential to enhance productivity in Kenya

The obstacle to optimal utilization of biogas technology is poor understanding of biogas microbiomes diversities over a wide geographical coverage. We performed random shotgun sequencing on twelve environmental samples. Randomized complete block design was utilized to assign the twelve treatments to four blocks, within eastern and central regions of Kenya. We obtained 42 million paired-end reads that were annotated against sixteen reference databases using two ENVO ontologies, prior to β-diversity studies. We identified 37 phyla, 65 classes and 132 orders. Bacteria dominated and comprised 28 phyla, 42 classes and 92 orders, conveying substrate's versatility in the treatments. Though, Fungi and Archaea comprised 5 phyla, the Fungi were richer; suggesting the importance of hydrolysis and fermentation in biogas production. High β-diversity within the taxa was largely linked to communities' metabolic capabilities. Clostridiales and Bacteroidales, the most prevalent guilds, metabolize organic macromolecules. The identified Cytophagales, Alteromonadales, Flavobacteriales, Fusobacteriales, Deferribacterales, Elusimicrobiales, Chlamydiales, Synergistales to mention but few, also catabolize macromolecules into smaller substrates to conserve energy. Furthermore, δ-Proteobacteria, Gloeobacteria and Clostridia affiliates syntrophically regulate PH2 and reduce metal to provide reducing equivalents. Methanomicrobiales and other Methanomicrobia species were the most prevalence Archaea, converting formate, CO2(g), acetate and methylated substrates into CH4(g). Thermococci, Thermoplasmata and Thermoprotei were among the sulfur and other metal reducing Archaea that contributed to redox balancing and other metabolism within treatments. Eukaryotes, mainly fungi were the least abundant guild, comprising largely Ascomycota and Basidiomycota species. Chytridiomycetes, Blastocladiomycetes and Mortierellomycetes were among the rare species, suggesting their metabolic and substrates limitations. Generally, we observed that environmental and treatment perturbations influenced communities' abundance, β-diversity and reactor performance largely through stochastic effect. Understanding diversity of biogas microbiomes over wide environmental variables and its' productivity provided insights into better management strategies that ameliorate biochemical limitations to effective biogas production.

Pan-genome study of Thermococcales reveals extensive genetic diversity and genetic evidence of thermophilic adaption

Thermococcales has a strong adaptability to extreme environments, which is of profound interest in explaining how complex life forms emerge on earth. However, their gene composition, thermal stability and evolution in hyperthermal environments are still little known. Here, we characterized the pan-genome architecture of 30 Thermococcales species to gain insight into their genetic properties, evolutionary patterns and specific metabolisms adapted to niches. We revealed an open pan-genome of Thermococcales comprising 6070 gene families that tend to increase with the availability of additional genomes. The genome contents of Thermococcales were flexible, with a series of genes experienced gene duplication, progressive divergence, or gene gain and loss events exhibiting distinct functional features. These archaea had concise types of heat shock proteins, such as HSP20, HSP60 and prefoldin, which were constrained by strong purifying selection that governed their conservative evolution. Furthermore, purifying selection forced genes involved in enzyme, motility, secretion system, defence system and chaperones to differ in functional constraints and their disparity in the rate of evolution may be related to adaptation to specific niche. These results deepened our understanding of genetic diversity and adaptation patterns of Thermococcales, and provided valuable research models for studying the metabolic traits of early life forms.

G-Quadruplexes in the Archaea Domain

The importance of unusual DNA structures in the regulation of basic cellular processes is an emerging field of research. Amongst local non-B DNA structures, G-quadruplexes (G4s) have gained in popularity during the last decade, and their presence and functional relevance at the DNA and RNA level has been demonstrated in a number of viral, bacterial, and eukaryotic genomes, including humans. Here, we performed the first systematic search of G4-forming sequences in all archaeal genomes available in the NCBI database. In this article, we investigate the presence and locations of G-quadruplex forming sequences using the G4Hunter algorithm. G-quadruplex-prone sequences were identified in all archaeal species, with highly significant differences in frequency, from 0.037 to 15.31 potential quadruplex sequences per kb. While G4 forming sequences were extremely abundant in Hadesarchaea archeon (strikingly, more than 50% of the Hadesarchaea archaeon isolate WYZ-LMO6 genome is a potential part of a G4-motif), they were very rare in the Parvarchaeota phylum. The presence of G-quadruplex forming sequences does not follow a random distribution with an over-representation in non-coding RNA, suggesting possible roles for ncRNA regulation. These data illustrate the unique and non-random localization of G-quadruplexes in Archaea.

Genome Replication in Thermococcus kodakarensis Independent of Cdc6 and an Origin of Replication

The initiation of DNA replication is typically tightly regulated by proteins that form initiation complexes at specific sequences known as replication origins. In Archaea and Eukaryotes, Cdc6, a near-universally conserved protein binds and facilitates the origin-dependent assembly of the replicative apparatus. TK1901 encodes Cdc6 in Thermococcus kodakarensis but, as we report here, TK1901 and the presumed origin of replication can be deleted from the genome of this hyperthermophilic Archaeon without any detectable effects on growth, genetic competence or the ability to support autonomous plasmid replication. All regions of the genome were equally represented in the sequences generated by whole genome sequencing of DNA isolated from T. kodakarensis strains with or without TK1901, inconsistent with DNA initiation occurring at one or few origins, and instead suggestive of replication initiating at many sites distributed throughout the genome. We were unable to generate strains lacking the recombination factors, RadA or RadB, consistent with T. kodakarensis cells, that are oligoploid (7–19 genomes per cell), employing a recombination-based mechanism of DNA replication. Deletion of the previously presumed origin region reduced the long-term viability of cultures supporting the possibility that retaining an origin-based mechanism of DNA initiation provides a survival mechanism for stationary phase cells with only one genome.

Flipping chromosomes in deep-sea archaea

One of the major mechanisms driving the evolution of all organisms is genomic rearrangement. In hyperthermophilic Archaea of the order Thermococcales, large chromosomal inversions occur so frequently that even closely related genomes are difficult to align. Clearly not resulting from the native homologous recombination machinery, the causative agent of these inversions has remained elusive. We present a model in which genomic inversions are catalyzed by the integrase enzyme encoded by a family of mobile genetic elements. We characterized the integrase from Thermococcus nautili plasmid pTN3 and showed that besides canonical site-specific reactions, it catalyzes low sequence specificity recombination reactions with the same outcome as homologous recombination events on DNA segments as short as 104bp both in vitro and in vivo, in contrast to other known tyrosine recombinases. Through serial culturing, we showed that the integrase-mediated divergence of T. nautili strains occurs at an astonishing rate, with at least four large-scale genomic inversions appearing within 60 generations. Our results and the ubiquitous distribution of pTN3-like integrated elements suggest that a major mechanism of evolution of an entire order of Archaea results from the activity of a selfish mobile genetic element.

Mobile elements (MEs) such as viruses, plasmids and transposons infect most living organisms and often encode recombinases promoting their insertion into cellular genomes. These insertions alter the genome of their host according to two main mechanisms. First, MEs provide new functions to the cell by integrating their own genetic information into the DNA of the host, at one or more locations. Secondly, cellular homologous recombination will act upon multiple integrated copies and produce a variety of large-scale chromosomal rearrangements. If such modifications are advantageous, they will spread into the population by natural selection. Typically, enzymes involved in cellular homologous recombination and the integration of MEs are distinct. We describe here a novel plasmid-encoded archaeal integrase which in addition to site-specific recombination can catalyze low sequence specificity recombination reactions akin to homologous recombination.

In vitro site-specific recombination mediated by the tyrosine recombinase XerA of Thermoplasma acidophilum

In organisms with circular chromosomes, such as bacteria and archaea, an odd number of homologous recombination events can generate a chromosome dimer. Such chromosome dimers cannot be segregated unless they are converted to monomers before cell division. In Escherichia coli, dimer-to-monomer conversion is mediated by the paralogous XerC and XerD recombinases at a specific dif site in the replication termination region. Dimer resolution requires the highly conserved cell division protein/chromosome translocase FtsK, and this site-specific chromosome resolution system is present or predicted in most bacteria. However, most archaea have only XerA, a homologue of the bacterial XerC/D proteins, but no homologues of FtsK. In addition, the molecular mechanism of XerA-mediated chromosome resolution in archaea has been less thoroughly elucidated than those of the corresponding bacterial systems. In this study, we identified two XerA-binding sites (dif1 and dif2) in the Thermoplasma acidophilum chromosome. In vitro site-specific recombination assays showed that dif2, but not dif1, serves as a target site for XerA-mediated chromosome resolution. Mutational analysis indicated that not only the core consensus sequence of dif2, but also its flanking regions play important roles in the recognition and recombination reactions mediated by XerA.

Elevated Rate of Genome Rearrangements in Radiation-Resistant Bacteria

A number of bacterial, archaeal, and eukaryotic species are known for their resistance to ionizing radiation. One of the challenges these species face is a potent environmental source of DNA double-strand breaks, potential drivers of genome structure evolution. Efficient and accurate DNA double-strand break repair systems have been demonstrated in several unrelated radiation-resistant species and are putative adaptations to the DNA damaging environment. Such adaptations are expected to compensate for the genome-destabilizing effect of environmental DNA damage and may be expected to result in a more conserved gene order in radiation-resistant species. However, here we show that rates of genome rearrangements, measured as loss of gene order conservation with time, are higher in radiation-resistant species in multiple, phylogenetically independent groups of bacteria. Comparison of indicators of selection for genome organization between radiation-resistant and phylogenetically matched, nonresistant species argues against tolerance to disruption of genome structure as a strategy for radiation resistance. Interestingly, an important mechanism affecting genome rearrangements in prokaryotes, the symmetrical inversions around the origin of DNA replication, shapes genome structure of both radiation-resistant and nonresistant species. In conclusion, the opposing effects of environmental DNA damage and DNA repair result in elevated rates of genome rearrangements in radiation-resistant bacteria.

Complete Genome Sequence of Hyperthermophilic Piezophilic Archaeon Palaeococcus pacificus DY20341T, Isolated from Deep-Sea Hydrothermal Sediments

We report the genome sequence of Palaeococcus pacificus DY20341^T, isolated from a sediment sample collected from eastern Pacific Ocean hydrothermal fields, which is the first report of a complete genome for a Palaeococcus species. The genome sequence will help to better understand differentiation phylogenetic relationships and evolution of several Thermococcales species.