Physical and Functional Compartmentalization of Archaeal Chromosomes.

The three-dimensional organization of chromosomes can have a profound impact on their replication and expression. The chromosomes of higher eukaryotes possess discrete compartments that are characterized by differing transcriptional activities. Contrastingly, most bacterial chromosomes have simpler organization with local domains, the boundaries of which are influenced by gene expression. Numerous studies have revealed that the higher-order architectures of bacterial and eukaryotic chromosomes are dependent on the actions of structural maintenance of chromosomes (SMC) superfamily protein complexes, in particular, the near-universal condensin complex. Intriguingly, however, many archaea, including members of the genus Sulfolobus do not encode canonical condensin. We describe chromosome conformation capture experiments on Sulfolobus species. These reveal the presence of distinct domains along Sulfolobus chromosomes that undergo discrete and specific higher-order interactions, thus defining two compartment types. We observe causal linkages between compartment identity, gene expression, and binding of a hitherto uncharacterized SMC superfamily protein that we term "coalescin."

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 kb 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.

Chromosome segregation in Archaea: SegA- and SegB-DNA complex structures provide insights into segrosome assembly.

Genome segregation is a vital process in all organisms. Chromosome partitioning remains obscure in Archaea, the third domain of life. Here, we investigated the SegAB system from Sulfolobus solfataricus. SegA is a ParA Walker-type ATPase and SegB is a site-specific DNA-binding protein. We determined the structures of both proteins and those of SegA-DNA and SegB-DNA complexes. The SegA structure revealed an atypical, novel non-sandwich dimer that binds DNA either in the presence or in the absence of ATP. The SegB structure disclosed a ribbon-helix-helix motif through which the protein binds DNA site specifically. The association of multiple interacting SegB dimers with the DNA results in a higher order chromatin-like structure. The unstructured SegB N-terminus plays an essential catalytic role in stimulating SegA ATPase activity and an architectural regulatory role in segrosome (SegA-SegB-DNA) formation. Electron microscopy results also provide a compact ring-like segrosome structure related to chromosome organization. These findings contribute a novel mechanistic perspective on archaeal chromosome segregation.

CRISPRi as an efficient tool for gene repression in archaea.

In the years following its discovery and characterization, the CRISPR-Cas system has been modified and converted into a multitude of applications for eukaryotes and bacteria, such as genome editing and gene regulation. Since no such method has been available for archaea, we developed a tool for gene repression in the haloarchaeon Haloferax volcanii by repurposing its endogenous type I-B CRISPR-Cas system. Here, we present the two possible approaches for gene repression as well as our workflow to achieve and assess gene knockdown, offer recommendations on protospacer selection and give some examples of genes we have successfully silenced.

Role of free DNA ends and protospacer adjacent motifs for CRISPR DNA uptake in Pyrococcus furiosus.

To acquire CRISPR-Cas immunity against invasive mobile genetic elements, prokaryotes must first integrate fragments of foreign DNA into their genomic CRISPR arrays for use in future invader silencing. Here, we found that the hyperthermophilic archaeaon, Pyrococcus furiosus, actively incorporates DNA fragments (spacers) from both plasmid (foreign) and host genome (self) sequences into its seven CRISPR loci. The majority of new spacers were derived from DNA immediately downstream from a 5'-CCN-3' protospacer adjacent motif (PAM) that is critical for invader targeting. Interestingly, spacers were preferentially acquired from genome or plasmid regions corresponding to active transposons, CRISPR loci, ribosomal RNA genes, rolling circle origins of replication, and areas where plasmids recombined with the host chromosome. A common feature of the highly sampled spacers is that they arise from DNA regions expected to undergo DNA nicking and/or double-strand breaks. Taken together with recent results from bacterial systems, our findings indicate that free DNA termini and PAMs are conserved features important for CRISPR spacer uptake in diverse prokaryotes and CRISPR-Cas systems. Moreover, lethal self-targeting by CRISPR systems may contribute to host genome stability by eliminating cells undergoing active transposon mobility or chromosomal uptake of autonomously replicating foreign mobile genetic elements.

Influence of Origin Recognition Complex Proteins on the Copy Numbers of Three Chromosomes in Haloferax volcanii.

Replication initiation in archaea involves a protein named ORC, Cdc6, or ORC1/Cdc6, which is homologous to the eukaryotic origin recognition complex (ORC) proteins and to the eukaryotic Cdc6. Archaeal replication origins are comprised of origin repeat regions and adjacent orc genes. Some archaea contain a single replication origin and a single orc gene, while others have more than one of each. Haloferax volcanii is exceptional because it contains, in total, six replication origins on three chromosomes and 16 orc genes. Phylogenetic trees were constructed that showed that orc gene duplications occurred at very different times in evolution. To unravel the influence of the ORC proteins on chromosome copy number and cellular fitness, it was attempted to generate deletion mutants of all 16 genes. A total of 12 single-gene deletion mutants could be generated, and only three orc gene turned out to be essential. For one gene, the deletion analysis failed. Growth analyses revealed that no deletion mutant had a growth defect, but some had a slight growth advantage compared to the wild type. Quantification of the chromosome copy numbers in the deletion mutants showed that all 12 ORC proteins influenced the copy numbers of one, two, or all three chromosomes. The lack of an ORC led to an increase or decrease of chromosome copy number. Therefore, chromosome copy numbers in Hfxvolcanii are regulated by an intricate network of ORC proteins. This is in contrast to other archaea, in which ORC proteins typically bind specifically to the adjacent origin.IMPORTANCE The core origins of archaea are comprised of a repeat region and an adjacent gene for an origin recognition complex (ORC) protein, which is homologous to eukaryotic ORC proteins. Haloferax volcanii is exceptional because it contains six replication origins on three chromosomes and an additional 10 orc genes that are not adjacent to an origin. This unique ORC protein repertoire was used to unravel the importance of core origin orc genes and of origin-remote orc genes. Remarkably, all ORC proteins influenced the copy number of at least one chromosome. Some of them influenced those of all three chromosomes, showing that cross-regulation in trans exists in Hfx. volcanii Furthermore, the evolution of the archaeal ORC protein family was analyzed.

Transcriptional Repressor TrmBL2 from Thermococcus kodakarensis Forms Filamentous Nucleoprotein Structures and Competes with Histones for DNA Binding in a Salt- and DNA Supercoiling-dependent Manner.

Architectural DNA proteins play important roles in the chromosomal DNA organization and global gene regulation in living cells. However, physiological functions of some DNA-binding proteins from archaea remain unclear. Recently, several abundant DNA-architectural proteins including histones, Alba, and TrmBL2 have been identified in model euryarchaeon Thermococcus kodakarensis. Although histones and Alba proteins have been previously characterized, the DNA binding properties of TrmBL2 and its interplay with the other major architectural proteins in the chromosomal DNA organization and gene transcription regulation remain largely unexplored. Here, we report single-DNA studies showing that at low ionic strength (<300 mM KCl), TrmBL2 binds to DNA largely in non-sequence-specific manner with positive cooperativity, resulting in formation of stiff nucleoprotein filamentous patches, whereas at high ionic strength (>300 mM KCl) TrmBL2 switches to more sequence-specific interaction, suggesting the presence of high affinity TrmBL2-filament nucleation sites. Furthermore, in vitro assays indicate the existence of DNA binding competition between TrmBL2 and archaeal histones B from T. kodakarensis, which can be strongly modulated by DNA supercoiling and ionic strength of surrounding solution. Overall, these results advance our understanding of TrmBL2 DNA binding properties and provide important insights into potential functions of architectural proteins in nucleoid organization and gene regulation in T. kodakarensis.

The chromosome copy number of the hyperthermophilic archaeon Thermococcus kodakarensis KOD1.

The euryarchaeon Thermococcus kodakarensis is a well-characterized anaerobic hyperthermophilic heterotroph and due to the availability of genetic engineering systems it has become one of the model organisms for studying Archaea. Despite this prominent role among the Euryarchaeota, no data about the ploidy level of this species is available. While polyploidy has been shown to exist in various Euryarchaeota, especially Halobacteria, the chromosome copy number of species belonging to one of the major orders within that phylum, i.e., the Thermococcales (including Thermococcus spp. and Pyrococcus spp.), has never been determined. This prompted us to investigate the chromosome copy number of T. kodakarensis. In this study, we demonstrate that T. kodakarensis is polyploid with a chromosome copy number that varies between 7 and 19 copies, depending on the growth phase. An apparent correlation between the presence of histones and polyploidy in Archaea is observed.

Four chromosome replication origins in the archaeon Pyrobaculum calidifontis.

Replication origins were mapped in hyperthermophilic crenarchaea, using high-throughput sequencing-based marker frequency analysis. We confirm previous origin mapping in Sulfolobus acidocaldarius, and demonstrate that the single chromosome of Pyrobaculum calidifontis contains four replication origins, the highest number detected in a prokaryotic organism. The relative positions of the origins in both organisms coincided with regions enriched in highly conserved (core) archaeal genes. We show that core gene distribution provides a useful tool for origin identification in archaea, and predict multiple replication origins in a range of species. One of the P. calidifontis origins was mapped in detail, and electrophoretic mobility shift assays demonstrated binding of the Cdc6/Orc1 replication initiator protein to a repeated sequence element, denoted Orb-1, within the origin. The high-throughput sequencing approach also allowed for an annotation update of both genomes, resulting in the restoration of open reading frames encoding proteins involved in, e.g., sugar, nitrate and energy metabolism, as well as in glycosylation and DNA repair.

Homologous recombination in the archaeon Sulfolobus acidocaldarius: effects of DNA substrates and mechanistic implications.

Although homologous recombination (HR) is known to influence the structure, stability, and evolution of microbial genomes, few of its functional properties have been measured in cells of hyperthermophilic archaea. The present study manipulated various properties of the parental DNAs in high-resolution assays of Sulfolobus acidocaldarius transformation, and measured the impact on the efficiency and pattern of marker transfer to the recipient chromosome. The relative orientation of homologous sequences, the type and position of chromosomal mutation being replaced, and the length of DNA flanking the marked region all affected the efficiency, linkage, tract continuity, and other parameters of marker transfer. Effects predicted specifically by the classical reciprocal-exchange model of HR were not observed. One analysis observed only 90 % linkage between markers defined by adjacent bases; in another series of experiments, sequence divergence up to 4 % had no detectable impact on overall efficiency of HR or on the co-transfer of a distal non-selected marker. The effects of introducing DNA via conjugation, rather than transformation, were more difficult to assess, but appeared to increase co-transfer (i.e. linkage) of relatively distant non-selected markers. The results indicate that HR events between gene-sized duplex DNAs and the S. acidocaldarius chromosome typically involve neither crossing over nor interference from a mismatch-activated anti-recombination system. Instead, the donor DNA may anneal to a transient chromosomal gap, as in the mechanism proposed for oligonucleotide-mediated transformation of Sulfolobus and other micro-organisms.

Evidence for a Xer/dif system for chromosome resolution in archaea.

Homologous recombination events between circular chromosomes, occurring during or after replication, can generate dimers that need to be converted to monomers prior to their segregation at cell division. In Escherichia coli, chromosome dimers are converted to monomers by two paralogous site-specific tyrosine recombinases of the Xer family (XerC/D). The Xer recombinases act at a specific dif site located in the replication termination region, assisted by the cell division protein FtsK. This chromosome resolution system has been predicted in most Bacteria and further characterized for some species. Archaea have circular chromosomes and an active homologous recombination system and should therefore resolve chromosome dimers. Most archaea harbour a single homologue of bacterial XerC/D proteins (XerA), but not of FtsK. Therefore, the role of XerA in chromosome resolution was unclear. Here, we have identified dif-like sites in archaeal genomes by using a combination of modeling and comparative genomics approaches. These sites are systematically located in replication termination regions. We validated our in silico prediction by showing that the XerA protein of Pyrococcus abyssi specifically recombines plasmids containing the predicted dif site in vitro. In contrast to the bacterial system, XerA can recombine dif sites in the absence of protein partners. Whereas Archaea and Bacteria use a completely different set of proteins for chromosome replication, our data strongly suggest that XerA is most likely used for chromosome resolution in Archaea.

Activity and transcriptional regulation of bacterial protein-like glycerol-3-phosphate dehydrogenase of the haloarchaea in Haloferax volcanii.

Glycerol is a primary energy source for heterotrophic haloarchaea and a major component of "salty" biodiesel waste. Glycerol is catabolized solely by glycerol kinase (encoded by glpK) to glycerol-3-phosphate (G3P) in Haloferax volcanii. Here we characterized the next critical step of this metabolic pathway: the conversion of G3P to dihydroxyacetone phosphate by G3P dehydrogenase (G3PDH). H. volcanii harbors two putative G3PDH operons: (i) glpA1B1C1, located on the chromosome within the neighborhood of glpK, and (ii) glpA2B2C2, on megaplasmid pHV4. Analysis of knockout strains revealed that glpA1(and not glpA2) is required for growth on glycerol. However, both glpA1 and glpA2 could complement a glpA1 knockout strain (when expressed from a strong promoter in trans) and were required for the total G3PDH activity of cell lysates. The glpA1B1C1, glpK, glpF(encoding a putative glycerol facilitator), and ptsH2(encoding a homolog of the bacterial phosphotransferase system protein Hpr) genes were transcriptionally linked and appeared to be under the control of a strong, G3P-inducible promoter upstream of glpA1. Overall, this study provides fundamental insights into glycerol metabolism in H. volcanii and enhances our understanding of central metabolic pathways of haloarchaea.

Generation of gene-level resolution chromosome contact maps in bacteria and archaea.

Chromosome conformation capture (Hi-C) has become a routine method for probing the 3D organization of genomes. However, when applied to bacteria and archaea, current protocols are expensive and limited in their resolution. By dissecting the different steps of published eukaryotic and prokaryotic Hi-C protocols, we have developed a cost- and time-effective approach to generate high-resolution (down to 500 bp - 1 kb) contact matrices of both bacteria and archaea genomes. For complete details on the use and execution of this protocol, please refer to Cockram et al. (2020).

Characterization of NADH oxidase/NADPH polysulfide oxidoreductase and its unexpected participation in oxygen sensitivity in an anaerobic hyperthermophilic archaeon.

Many genomes of anaerobic hyperthermophiles encode multiple homologs of NAD(P)H oxidase that are thought to function in response to oxidative stress. We investigated one of the seven NAD(P)H oxidase homologs (TK1481) in the sulfur-reducing hyperthermophilic archaeon Thermococcus kodakarensis, focusing on the catalytic properties and roles in oxidative-stress defense and sulfur-dependent energy conservation. The recombinant form of TK1481 exhibited both NAD(P)H oxidase and NAD(P)H:polysulfide oxidoreductase activities. The enzyme also possessed low NAD(P)H peroxidase and NAD(P)H:elemental sulfur oxidoreductase activities under anaerobic conditions. A mutant form of the enzyme, in which the putative redox-active residue Cys43 was replaced by Ala, still showed NADH-dependent flavin adenine dinucleotide (FAD) reduction activity. Although it also retained successive oxidase and anaerobic peroxidase activities, the ability to reduce polysulfide and sulfur was completely lost, suggesting the specific reactivity of the Cys43 residue for sulfur. To evaluate the physiological function of TK1481, we constructed a gene deletant, ΔTK1481, and mutant KUTK1481C43A, into which two base mutations altering Cys43 of TK1481 to Ala were introduced. ΔTK1481 exhibited growth properties nearly identical to those of the parent strain, KU216, in sulfur-containing media. Interestingly, in the absence of elemental sulfur, the growth of ΔTK1481 was not affected by dissolved oxygen, whereas the growth of KU216 and KUTK1481C43A was significantly impaired. These results indicate that although TK1481 does not play a critical role in either sulfur reduction or the response to oxidative stress, the NAD(P)H oxidase activity of TK1481 unexpectedly participates in the oxygen sensitivity of the hyperthermophilic archaeon T. kodakarensis in the absence of sulfur.

Replication-biased genome organisation in the crenarchaeon Sulfolobus.

BACKGROUND: Species of the crenarchaeon Sulfolobus harbour three replication origins in their single circular chromosome that are synchronously initiated during replication. RESULTS: We demonstrate that global gene expression in two Sulfolobus species is highly biased, such that early replicating genome regions are more highly expressed at all three origins. The bias by far exceeds what would be anticipated by gene dosage effects alone. In addition, early replicating regions are denser in archaeal core genes (enriched in essential functions), display lower intergenic distances, and are devoid of mobile genetic elements. CONCLUSION: The strong replication-biased structuring of the Sulfolobus chromosome implies that the multiple replication origins serve purposes other than simply shortening the time required for replication. The higher-level chromosomal organisation could be of importance for minimizing the impact of DNA damage, and may also be linked to transcriptional regulation.

Genetic and physical mapping of DNA replication origins in Haloferax volcanii.

The halophilic archaeon Haloferax volcanii has a multireplicon genome, consisting of a main chromosome, three secondary chromosomes, and a plasmid. Genes for the initiator protein Cdc6/Orc1, which are commonly located adjacent to archaeal origins of DNA replication, are found on all replicons except plasmid pHV2. However, prediction of DNA replication origins in H. volcanii is complicated by the fact that this species has no less than 14 cdc6/orc1 genes. We have used a combination of genetic, biochemical, and bioinformatic approaches to map DNA replication origins in H. volcanii. Five autonomously replicating sequences were found adjacent to cdc6/orc1 genes and replication initiation point mapping was used to confirm that these sequences function as bidirectional DNA replication origins in vivo. Pulsed field gel analyses revealed that cdc6/orc1-associated replication origins are distributed not only on the main chromosome (2.9 Mb) but also on pHV1 (86 kb), pHV3 (442 kb), and pHV4 (690 kb) replicons. Gene inactivation studies indicate that linkage of the initiator gene to the origin is not required for replication initiation, and genetic tests with autonomously replicating plasmids suggest that the origin located on pHV1 and pHV4 may be dominant to the principal chromosomal origin. The replication origins we have identified appear to show a functional hierarchy or differential usage, which might reflect the different replication requirements of their respective chromosomes. We propose that duplication of H. volcanii replication origins was a prerequisite for the multireplicon structure of this genome, and that this might provide a means for chromosome-specific replication control under certain growth conditions. Our observations also suggest that H. volcanii is an ideal organism for studying how replication of four replicons is regulated in the context of the archaeal cell cycle.

Quantitative analysis of replication-related mutation and selection pressures in bacterial chromosomes and plasmids using generalised GC skew index.

BACKGROUND: Due to their bi-directional replication machinery starting from a single finite origin, bacterial genomes show characteristic nucleotide compositional bias between the two replichores, which can be visualised through GC skew or (C-G)/(C+G). Although this polarisation is used for computational prediction of replication origins in many bacterial genomes, the degree of GC skew visibility varies widely among different species, necessitating a quantitative measurement of GC skew strength in order to provide confidence measures for GC skew-based predictions of replication origins. RESULTS: Here we discuss a quantitative index for the measurement of GC skew strength, named the generalised GC skew index (gGCSI), which is applicable to genomes of any length, including bacterial chromosomes and plasmids. We demonstrate that gGCSI is independent of the window size and can thus be used to compare genomes with different sizes, such as bacterial chromosomes and plasmids. It can suggest the existence of different replication mechanisms in archaea and of rolling-circle replication in plasmids. Correlation of gGCSI values between plasmids and their corresponding host chromosomes suggests that within the same strain, these replicons have reproduced using the same replication machinery and thus exhibit similar strengths of replication strand skew. CONCLUSIONS: gGCSI can be applied to genomes of any length and thus allows comparative study of replication-related mutation and selection pressures in genomes of different lengths such as bacterial chromosomes and plasmids. Using gGCSI, we showed that replication-related mutation or selection pressure is similar for replicons with similar machinery.

Recombination of synthetic oligonucleotides with prokaryotic chromosomes: substrate requirements of the Escherichia coli/lambdaRed and Sulfolobus acidocaldarius recombination systems.

In order to reveal functional properties of recombination involving short ssDNAs in hyperthermophilic archaea, we evaluated oligonucleotide-mediated transformation (OMT) in Sulfolobus acidocaldarius and Escherichia coli as a function of the molecular properties of the ssDNA substrates. Unmodified ssDNAs as short as 20-22 nt yielded recombinants in both organisms, as did longer DNAs forming as few as 2-5 base pairs on one side of the genomic mutation. The two OMT systems showed similar responses to certain end modifications of the oligonucleotides, but E. coli was found to require a 5' phosphate on 5'-limited ssDNA whereas this requirement was not evident in S. acidocaldarius. The ability of both E. coli and S. acidocaldarius to incorporate short, mismatched ssDNAs into their genomes raises questions about the biological significance of this capability, including its phylogenetic distribution among microorganisms and its impact on genome stability. These questions seem particularly relevant for S. acidocaldarius, as this archaeon has natural competence for OMT, encodes no MutSL homologues and thrives under environmental conditions that accelerate DNA decomposition.

Expanding and understanding the genetic toolbox of the hyperthermophilic genus Sulfolobus.

Although Sulfolobus species are among the best studied archaeal micro-organisms, the development and availability of genetic tools has lagged behind. In the present paper, we discuss the latest progress in understanding recombination events of exogenous DNA into the chromosomes of Sulfolobus solfataricus and Sulfolobus acidocaldarius and their application in the construction of targeted-deletion mutant strains.

Ecophysiological significance of scale-dependent patterns in prokaryotic genomes unveiled by a combination of statistic and genometric analyses.

We combined genometric (DNA walks) and statistical (detrended fluctuation analysis) methods on 456 prokaryotic chromosomes from 309 different bacterial and archaeal species to look for specific patterns and long-range correlations along the genome and relate them to ecological lifestyles. The position of each nucleotide along the complete genome sequence was plotted on an orthogonal plane (DNA landscape), and fluctuation analysis applied to the DNA walk series showed a long-range correlation in contrast to the lack of correlation for artificially generated genomes. Different features in the DNA landscapes among genomes from different ecological and metabolic groups of prokaryotes appeared with the combined analysis. Transition from hyperthermophilic to psychrophilic environments could have been related to more complex structural adaptations in microbial genomes, whereas for other environmental factors such as pH and salinity this effect would have been smaller. Prokaryotes with domain-specific metabolisms, such as photoautotrophy in Bacteria and methanogenesis in Archaea, showed consistent differences in genome correlation structure. Overall, we show that, beyond the relative proportion of nucleotides, correlation properties derived from their sequential position within the genome hide relevant phylogenetic and ecological information. This can be studied by combining genometric and statistical physics methods, leading to a reduction of genome complexity to a few useful descriptors.