Deletion of the archaeal histone in Methanosarcina mazei Gö1 results in reduced growth and genomic transcription

DNA-bridging by an archaeal histone variant via a unique tetramerisation interface

In eukaryotes, histone paralogues form obligate heterodimers such as H3/H4 and H2A/H2B that assemble into octameric nucleosome particles. Archaeal histones are dimeric and assemble on DNA into 'hypernucleosome' particles of varying sizes with each dimer wrapping 30 bp of DNA. These are composed of canonical and variant histone paralogues, but the function of these variants is poorly understood. Here, we characterise the structure and function of the histone paralogue MJ1647 from Methanocaldococcus jannaschii that has a unique C-terminal extension enabling homotetramerisation. The 1.9 Å X-ray structure of a dimeric MJ1647 species, structural modelling of the tetramer, and site-directed mutagenesis reveal that the C-terminal tetramerization module consists of two alpha helices in a handshake arrangement. Unlike canonical histones, MJ1647 tetramers can bridge two DNA molecules in vitro. Using single-molecule tethered particle motion and DNA binding assays, we show that MJ1647 tetramers bind ~60 bp DNA and compact DNA in a highly cooperative manner. We furthermore show that MJ1647 effectively competes with the transcription machinery to block access to the promoter in vitro. To the best of our knowledge, MJ1647 is the first histone shown to have DNA bridging properties, which has important implications for genome structure and gene expression in archaea.

Histone variants in archaea - An undiscovered country

Exchanging core histones in the nucleosome for paralogous variants can have important functional ramifications. Many of these variants, and their physiological roles, have been characterized in exquisite detail in model eukaryotes, including humans. In comparison, our knowledge of histone biology in archaea remains rudimentary. This is true in particular for our knowledge of histone variants. Many archaea encode several histone genes that differ in sequence, but do these paralogs make distinct, adaptive contributions to genome organization and regulation in a manner comparable to eukaryotes? Below, we review what we know about histone variants in archaea at the level of structure, regulation, and evolution. In all areas, our knowledge pales when compared to the wealth of insight that has been gathered for eukaryotes. Recent findings, however, provide tantalizing glimpses into a rich and largely undiscovered country that is at times familiar and eukaryote-like and at times strange and uniquely archaeal. We sketch a preliminary roadmap for further exploration of this country; an undertaking that may ultimately shed light not only on chromatin biology in archaea but also on the origin of histone-based chromatin in eukaryotes.

The Hypersaline Archaeal Histones HpyA and HstA Are DNA Binding Proteins That Defy Categorization According to Commonly Used Functional Criteria

Histone proteins are found across diverse lineages of Archaea, many of which package DNA and form chromatin. However, previous research has led to the hypothesis that the histone-like proteins of high-salt-adapted archaea, or halophiles, function differently. The sole histone protein encoded by the model halophilic species Halobacterium salinarum, HpyA, is nonessential and expressed at levels too low to enable genome-wide DNA packaging. Instead, HpyA mediates the transcriptional response to salt stress. Here we compare the features of genome-wide binding of HpyA to those of HstA, the sole histone of another model halophile, Haloferax volcanii. hstA, like hpyA, is a nonessential gene. To better understand HpyA and HstA functions, protein-DNA binding data (chromatin immunoprecipitation sequencing [ChIP-seq]) of these halophilic histones are compared to publicly available ChIP-seq data from DNA binding proteins across all domains of life, including transcription factors (TFs), nucleoid-associated proteins (NAPs), and histones. These analyses demonstrate that HpyA and HstA bind the genome infrequently in discrete regions, which is similar to TFs but unlike NAPs, which bind a much larger genomic fraction. However, unlike TFs that typically bind in intergenic regions, HpyA and HstA binding sites are located in both coding and intergenic regions. The genome-wide dinucleotide periodicity known to facilitate histone binding was undetectable in the genomes of both species. Instead, TF-like and histone-like binding sequence preferences were detected for HstA and HpyA, respectively. Taken together, these data suggest that halophilic archaeal histones are unlikely to facilitate genome-wide chromatin formation and that their function defies categorization as a TF, NAP, or histone. IMPORTANCE Most cells in eukaryotic species-from yeast to humans-possess histone proteins that pack and unpack DNA in response to environmental cues. These essential proteins regulate genes necessary for important cellular processes, including development and stress protection. Although the histone fold domain originated in the domain of life Archaea, the function of archaeal histone-like proteins is not well understood relative to those of eukaryotes. We recently discovered that, unlike histones of eukaryotes, histones in hypersaline-adapted archaeal species do not package DNA and can act as transcription factors (TFs) to regulate stress response gene expression. However, the function of histones across species of hypersaline-adapted archaea still remains unclear. Here, we compare hypersaline histone function to a variety of DNA binding proteins across the tree of life, revealing histone-like behavior in some respects and specific transcriptional regulatory function in others.

Growth temperature and chromatinization in archaea

DNA in cells is associated with proteins that constrain its structure and affect DNA-templated processes including transcription and replication. HU and histones are the main constituents of chromatin in bacteria and eukaryotes, respectively, with few exceptions. Archaea, in contrast, have diverse repertoires of nucleoid-associated proteins (NAPs). To analyse the evolutionary and ecological drivers of this diversity, we combined a phylogenomic survey of known and predicted NAPs with quantitative proteomic data. We identify the Diaforarchaea as a hotbed of NAP gain and loss, and experimentally validate candidate NAPs in two members of this clade, Thermoplasma volcanium and Methanomassiliicoccus luminyensis. Proteomic analysis across a diverse sample of 19 archaea revealed that NAP investment varies from <0.03% to >5% of total protein. This variation is predicted by growth temperature. We propose that high levels of chromatinization have evolved as a mechanism to prevent uncontrolled helix denaturation at higher temperatures, with implications for the origin of chromatin in both archaea and eukaryotes.

A phylogenetic and proteomic reconstruction of eukaryotic chromatin evolution

Histones and associated chromatin proteins have essential functions in eukaryotic genome organization and regulation. Despite this fundamental role in eukaryotic cell biology, we lack a phylogenetically-comprehensive understanding of chromatin evolution. Here, we combine comparative proteomics and genomics analysis of chromatin in eukaryotes and archaea. Proteomics uncovers the existence of histone post-translational modifications in Archaea. However, archaeal histone modifications are scarce, in contrast with the highly conserved and abundant marks we identify across eukaryotes. Phylogenetic analysis reveals that chromatin-associated catalytic functions (e.g., methyltransferases) have pre-eukaryotic origins, whereas histone mark readers and chaperones are eukaryotic innovations. We show that further chromatin evolution is characterized by expansion of readers, including capture by transposable elements and viruses. Overall, our study infers detailed evolutionary history of eukaryotic chromatin: from its archaeal roots, through the emergence of nucleosome-based regulation in the eukaryotic ancestor, to the diversification of chromatin regulators and their hijacking by genomic parasites.

An archaeal histone-like protein regulates gene expression in response to salt stress

Histones, ubiquitous in eukaryotes as DNA-packing proteins, find their evolutionary origins in archaea. Unlike the characterized histone proteins of a number of methanogenic and themophilic archaea, previous research indicated that HpyA, the sole histone encoded in the model halophile Halobacterium salinarum, is not involved in DNA packaging. Instead, it was found to have widespread but subtle effects on gene expression and to maintain wild type cell morphology. However, the precise function of halophilic histone-like proteins remain unclear. Here we use quantitative phenotyping, genetics, and functional genomics to investigate HpyA function. These experiments revealed that HpyA is important for growth and rod-shaped morphology in reduced salinity. HpyA preferentially binds DNA at discrete genomic sites under low salt to regulate expression of ion uptake, particularly iron. HpyA also globally but indirectly activates other ion uptake and nucleotide biosynthesis pathways in a salt-dependent manner. Taken together, these results demonstrate an alternative function for an archaeal histone-like protein as a transcriptional regulator, with its function tuned to the physiological stressors of the hypersaline environment.

The CARF Protein MM_0565 Affects Transcription of the Casposon-Encoded cas1-solo Gene in Methanosarcina mazei Gö1

Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) loci are found in bacterial and archaeal genomes where they provide the molecular machinery for acquisition of immunity against foreign DNA. In addition to the cas genes fundamentally required for CRISPR activity, a second class of genes is associated with the CRISPR loci, of which many have no reported function in CRISPR-mediated immunity. Here, we characterize MM_0565 associated to the type I-B CRISPR-locus of Methanosarcina mazei Gö1. We show that purified MM_0565 composed of a CRISPR-Cas Associated Rossmann Fold (CARF) and a winged helix-turn-helix domain forms a dimer in solution; in vivo, the dimeric MM_0565 is strongly stabilized under high salt stress. While direct effects on CRISPR-Cas transcription were not detected by genetic approaches, specific binding of MM_0565 to the leader region of both CRISPR-Cas systems was observed by microscale thermophoresis and electromobility shift assays. Moreover, overexpression of MM_0565 strongly induced transcription of the cas1-solo gene located in the recently reported casposon, the gene product of which shows high similarity to classical Cas1 proteins. Based on our findings, and taking the absence of the expressed CRISPR locus-encoded Cas1 protein into account, we hypothesize that MM_0565 might modulate the activity of the CRISPR systems on different levels.

Different Proteins Mediate Step-Wise Chromosome Architectures in Thermoplasma acidophilum and Pyrobaculum calidifontis

Archaeal species encode a variety of distinct lineage-specific chromosomal proteins. We have previously shown that in Thermococcus kodakarensis, histone, Alba, and TrmBL2 play distinct roles in chromosome organization. Although our understanding of individual archaeal chromosomal proteins has been advancing, how archaeal chromosomes are folded into higher-order structures and how they are regulated are largely unknown. Here, we investigated the primary and higher-order structures of archaeal chromosomes from different archaeal lineages. Atomic force microscopy of chromosome spreads out of Thermoplasma acidophilum and Pyrobaculum calidifontis cells revealed 10-nm fibers and 30–40-nm globular structures, suggesting the occurrence of higher-order chromosomal folding. Our results also indicated that chromosome compaction occurs toward the stationary phase. Micrococcal nuclease digestion indicated that fundamental structural units of the chromosome exist in T. acidophilum and T. kodakarensis but not in P. calidifontis or Sulfolobus solfataricus. In vitro reconstitution showed that, in T. acidophilum, the bacterial HU protein homolog HTa formed a 6-nm fiber by wrapping DNA, and that Alba was responsible for the formation of the 10-nm fiber by binding along the DNA without wrapping. Remarkably, Alba could form different higher-order complexes with histone or HTa on DNA in vitro. Mass spectrometry detected HTa and Rad50 in the T. acidophilum chromosome but not in other species. A putative transcriptional regulator of the AsnC/Lrp family (Pcal_1183) was detected on the P. calidifontis chromosome, but not on that of other species studied. Putative membrane-associated proteins were detected in the chromosomes of the three archaeal species studied, including T. acidophilum, P. calidifontis, and T. kodakarensis. Collectively, our data show that Archaea use different combinations of proteins to achieve chromosomal architecture and functional regulation.

The DNA-binding protein HTa from Thermoplasma acidophilum is an archaeal histone analog

Histones are a principal constituent of chromatin in eukaryotes and fundamental to our understanding of eukaryotic gene regulation. In archaea, histones are widespread but not universal: several lineages have lost histone genes. What prompted or facilitated these losses and how archaea without histones organize their chromatin remains largely unknown. Here, we elucidate primary chromatin architecture in an archaeon without histones, Thermoplasma acidophilum, which harbors a HU family protein (HTa) that protects part of the genome from micrococcal nuclease digestion. Charting HTa-based chromatin architecture in vitro, in vivo and in an HTa-expressing E. coli strain, we present evidence that HTa is an archaeal histone analog. HTa preferentially binds to GC-rich sequences, exhibits invariant positioning throughout the growth cycle, and shows archaeal histone-like oligomerization behavior. Our results suggest that HTa, a DNA-binding protein of bacterial origin, has converged onto an architectural role filled by histones in other archaea.

Chromatinization of Escherichia coli with archaeal histones

Nucleosomes restrict DNA accessibility throughout eukaryotic genomes, with repercussions for replication, transcription, and other DNA-templated processes. How this globally restrictive organization emerged during evolution remains poorly understood. Here, to better understand the challenges associated with establishing globally restrictive chromatin, we express histones in a naive system that has not evolved to deal with nucleosomal structures: Escherichia coli. We find that histone proteins from the archaeon Methanothermus fervidus assemble on the E. coli chromosome in vivo and protect DNA from micrococcal nuclease digestion, allowing us to map binding footprints genome-wide. We show that higher nucleosome occupancy at promoters is associated with lower transcript levels, consistent with local repressive effects. Surprisingly, however, this sudden enforced chromatinization has only mild repercussions for growth unless cells experience topological stress. Our results suggest that histones can become established as ubiquitous chromatin proteins without interfering critically with key DNA-templated processes.

Global Transcriptional Programs in Archaea Share Features with the Eukaryotic Environmental Stress Response

The environmental stress response (ESR), a global transcriptional program originally identified in yeast, is characterized by a rapid and transient transcriptional response composed of large, oppositely regulated gene clusters. Genes induced during the ESR encode core components of stress tolerance, macromolecular repair, and maintenance of homeostasis. In this review, we investigate the possibility for conservation of the ESR across the eukaryotic and archaeal domains of life. We first re-analyze existing transcriptomics data sets to illustrate that a similar transcriptional response is identifiable in Halobacterium salinarum, an archaeal model organism. To substantiate the archaeal ESR, we calculated gene-by-gene correlations, gene function enrichment, and comparison of temporal dynamics. We note reported examples of variation in the ESR across fungi, then synthesize high-level trends present in expression data of other archaeal species. In particular, we emphasize the need for additional high-throughput time series expression data to further characterize stress-responsive transcriptional programs in the Archaea. Together, this review explores an open question regarding features of global transcriptional stress response programs shared across domains of life.

The Role of Archaeal Chromatin in Transcription

Genomic organization impacts accessibility and movement of information processing systems along DNA. DNA-bound proteins dynamically dictate gene expression and provide regulatory potential to tune transcription rates to match ever-changing environmental conditions. Archaeal genomes are typically small, circular, gene dense, and organized either by histone proteins that are homologous to their eukaryotic counterparts, or small basic proteins that function analogously to bacterial nucleoid proteins. We review here how archaeal genomes are organized and how such organization impacts archaeal gene expression, focusing on conserved DNA-binding proteins within the clade and the factors that are known to impact transcription initiation and elongation within protein-bound genomes.

Small RNAs Involved in Regulation of Nitrogen Metabolism

Global (metabolic) regulatory networks allow microorganisms to survive periods of nitrogen starvation or general nutrient stress. Uptake and utilization of various nitrogen sources are thus commonly tightly regulated in Prokarya (Bacteria and Archaea) in response to available nitrogen sources. Those well-studied regulations occur mainly at the transcriptional and posttranslational level. Surprisingly, and in contrast to their involvement in most other stress responses, small RNAs (sRNAs) involved in the response to environmental nitrogen fluctuations are only rarely reported. In addition to sRNAs indirectly affecting nitrogen metabolism, only recently it was demonstrated that three sRNAs were directly involved in regulation of nitrogen metabolism in response to changes in available nitrogen sources. All three trans-acting sRNAs are under direct transcriptional control of global nitrogen regulators and affect expression of components of nitrogen metabolism (glutamine synthetase, nitrogenase, and PII-like proteins) by either masking the ribosome binding site and thus inhibiting translation initiation or stabilizing the respective target mRNAs. Most likely, there are many more sRNAs and other types of noncoding RNAs, e.g., riboswitches, involved in the regulation of nitrogen metabolism in Prokarya that remain to be uncovered. The present review summarizes the current knowledge on sRNAs involved in nitrogen metabolism and their biological functions and targets.

Cross-cleavage activity of Cas6b in crRNA processing of two different CRISPR-Cas systems in Methanosarcina mazei Gö1

The clustered regularly interspaced short palindromic repeat (CRISPR) system is a prokaryotic adaptive defense system against foreign nucleic acids. In the methanoarchaeon Methanosarcina mazei Gö1, two types of CRISPR-Cas systems are present (type I-B and type III-C). Both loci encode a Cas6 endonuclease, Cas6b-IB and Cas6b-IIIC, typically responsible for maturation of functional short CRISPR RNAs (crRNAs). To evaluate potential cross cleavage activity, we biochemically characterized both Cas6b proteins regarding their crRNA binding behavior and their ability to process pre-crRNA from the respective CRISPR array in vivo. Maturation of crRNA was studied in the respective single deletion mutants by northern blot and RNA-Seq analysis demonstrating that in vivo primarily Cas6b-IB is responsible for crRNA processing of both CRISPR arrays. Tentative protein level evidence for the translation of both Cas6b proteins under standard growth conditions was detected, arguing for different activities or a potential non-redundant role of Cas6b-IIIC within the cell. Conservation of both Cas6 endonucleases was observed in several other M. mazei isolates, though a wide variety was displayed. In general, repeat and leader sequence conservation revealed a close correlation in the M. mazei strains. The repeat sequences from both CRISPR arrays from M. mazei Gö1 contain the same sequence motif with differences only in two nucleotides. These data stand in contrast to all other analyzed M. mazei isolates, which have at least one additional CRISPR array with repeats belonging to another sequence motif. This conforms to the finding that Cas6b-IB is the crucial and functional endonuclease in M. mazei Gö1. Abbreviations: sRNA: small RNA; crRNA: CRISPR RNA; pre-crRNAs: Precursor CRISPR RNA; CRISPR: clustered regularly interspaced short palindromic repeats; Cas: CRISPR associated; nt: nucleotide; RNP: ribonucleoprotein; RBS: ribosome binding site.

Transcriptional Regulation in Archaea: From Individual Genes to Global Regulatory Networks

Archaea are major contributors to biogeochemical cycles, possess unique metabolic capabilities, and resist extreme stress. To regulate the expression of genes encoding these unique programs, archaeal cells use gene regulatory networks (GRNs) composed of transcription factor proteins and their target genes. Recent developments in genetics, genomics, and computational methods used with archaeal model organisms have enabled the mapping and prediction of global GRN structures. Experimental tests of these predictions have revealed the dynamical function of GRNs in response to environmental variation. Here, we review recent progress made in this area, from investigating the mechanisms of transcriptional regulation of individual genes to small-scale subnetworks and genome-wide global networks. At each level, archaeal GRNs consist of a hybrid of bacterial, eukaryotic, and uniquely archaeal mechanisms. We discuss this theme from the perspective of the role of individual transcription factors in genome-wide regulation, how these proteins interact to compile GRN topological structures, and how these topologies lead to emergent, high-level GRN functions. We conclude by discussing how systems biology approaches are a fruitful avenue for addressing remaining challenges, such as discovering gene function and the evolution of GRNs.

Methanosarcina Spherical Virus, a Novel Archaeal Lytic Virus Targeting Methanosarcina Strains

A novel archaeal lytic virus targeting species of the genus Methanosarcina was isolated using Methanosarcina mazei strain Gö1 as the host. Due to its spherical morphology, the virus was designated Methanosarcina spherical virus (MetSV). Molecular analysis demonstrated that MetSV contains double-stranded linear DNA with a genome size of 10,567 bp containing 22 open reading frames (ORFs), all oriented in the same direction. Functions were predicted for some of these ORFs, i.e., such as DNA polymerase, ATPase, and DNA-binding protein as well as envelope (structural) protein. MetSV-derived spacers in CRISPR loci were detected in several published Methanosarcina draft genomes using bioinformatic tools, revealing a potential protospacer-adjacent motif (PAM) motif (TTA/T). Transcription and expression of several predicted viral ORFs were validated by reverse transcription-PCR (RT-PCR), PAGE analysis, and liquid chromatography-mass spectrometry (LC-MS)-based proteomics. Analysis of core lipids by atmospheric pressure chemical ionization (APCI) mass spectrometry showed that MetSV and Methanosarcina mazei both contain archaeol and glycerol dialkyl glycerol tetraether without a cyclopentane moiety (GDGT-0). The MetSV host range is limited to Methanosarcina strains growing as single cells (M. mazei, Methanosarcina barkeri and Methanosarcina soligelidi). In contrast, strains growing as sarcina-like aggregates were apparently protected from infection. Heterogeneity related to morphology phases in M. mazei cultures allowed acquisition of resistance to MetSV after challenge by growing cultures as sarcina-like aggregates. CRISPR/Cas-mediated resistance was excluded since neither of the two CRISPR arrays showed MetSV-derived spacer acquisition. Based on these findings, we propose that changing the morphology from single cells to sarcina-like aggregates upon rearrangement of the envelope structure prevents infection and subsequent lysis by MetSV.IMPORTANCE Methanoarchaea are among the most abundant organisms on the planet since they are present in high numbers in major anaerobic environments. They convert various carbon sources, e.g., acetate, methylamines, or methanol, to methane and carbon dioxide; thus, they have a significant impact on the emission of major greenhouse gases. Today, very little is known about viruses specifically infecting methanoarchaea that most probably impact the abundance of methanoarchaea in microbial consortia. Here, we characterize the first identified Methanosarcina-infecting virus (MetSV) and show a mechanism for acquiring resistance against MetSV. Based on our results, we propose that growth as sarcina-like aggregates prevents infection and subsequent lysis. These findings allow new insights into the virus-host relationship in methanogenic community structures, their dynamics, and their phase heterogeneity. Moreover, the availability of a specific virus provides new possibilities to deepen our knowledge of the defense mechanisms of potential hosts and offers tools for genetic manipulation.

sRNA154 a newly identified regulator of nitrogen fixation in Methanosarcina mazei strain Gö1

Trans-encoded sRNA₁₅₄ is exclusively expressed under nitrogen (N)-deficiency in Methanosarcina mazei strain Gö1. The sRNA₁₅₄ deletion strain showed a significant decrease in growth under N-limitation, pointing toward a regulatory role of sRNA₁₅₄ in N-metabolism. Aiming to elucidate its regulatory function we characterized sRNA₁₅₄ by means of biochemical and genetic approaches. 24 homologs of sRNA₁₅₄ were identified in recently reported draft genomes of Methanosarcina strains, demonstrating high conservation in sequence and predicted secondary structure with two highly conserved single stranded loops. Transcriptome studies of sRNA₁₅₄ deletion mutants by an RNA-seq approach uncovered nifH- and nrpA-mRNA, encoding the α-subunit of nitrogenase and the transcriptional activator of the nitrogen fixation (nif)-operon, as potential targets besides other components of the N-metabolism. Furthermore, results obtained from stability, complementation and western blot analysis, as well as in silico target predictions combined with electrophoretic mobility shift-assays_, argue for a stabilizing effect of sRNA₁₅₄ on the polycistronic nif-mRNA and nrpA-mRNA by binding with both loops. Further identified N-related targets were studied, which demonstrates that translation initiation of glnA₂-mRNA, encoding glutamine synthetase2, appears to be affected by sRNA₁₅₄ masking the ribosome binding site, whereas glnA₁-mRNA appears to be stabilized by sRNA₁₅₄. Overall, we propose that sRNA₁₅₄ has a crucial regulatory role in N-metabolism in M. mazei by stabilizing the polycistronic mRNA encoding nitrogenase and glnA₁-mRNA, as well as allowing a feed forward regulation of nif-gene expression by stabilizing nrpA-mRNA. Consequently, sRNA₁₅₄ represents the first archaeal sRNA, for which a positive posttranscriptional regulation is demonstrated as well as inhibition of translation initiation.

Transcription Regulation in Archaea

The known diversity of metabolic strategies and physiological adaptations of archaeal species to extreme environments is extraordinary. Accurate and responsive mechanisms to ensure that gene expression patterns match the needs of the cell necessitate regulatory strategies that control the activities and output of the archaeal transcription apparatus. Archaea are reliant on a single RNA polymerase for all transcription, and many of the known regulatory mechanisms employed for archaeal transcription mimic strategies also employed for eukaryotic and bacterial species. Novel mechanisms of transcription regulation have become apparent by increasingly sophisticated in vivo and in vitro investigations of archaeal species. This review emphasizes recent progress in understanding archaeal transcription regulatory mechanisms and highlights insights gained from studies of the influence of archaeal chromatin on transcription.

Growth-Phase-Specific Modulation of Cell Morphology and Gene Expression by an Archaeal Histone Protein

In all three domains of life, organisms use nonspecific DNA-binding proteins to compact and organize the genome as well as to regulate transcription on a global scale. Histone is the primary eukaryotic nucleoprotein, and its evolutionary roots can be traced to the archaea. However, not all archaea use this protein as the primary DNA-packaging component, raising questions regarding the role of histones in archaeal chromatin function. Here, quantitative phenotyping, transcriptomic, and proteomic assays were performed on deletion and overexpression mutants of the sole histone protein of the hypersaline-adapted haloarchaeal model organism Halobacterium salinarum. This protein is highly conserved among all sequenced haloarchaeal species and maintains hallmark residues required for eukaryotic histone functions. Surprisingly, despite this conservation at the sequence level, unlike in other archaea or eukaryotes, H. salinarum histone is required to regulate cell shape but is not necessary for survival. Genome-wide expression changes in histone deletion strains were global, significant but subtle in terms of fold change, bidirectional, and growth phase dependent. Mass spectrometric proteomic identification of proteins from chromatin enrichments yielded levels of histone and putative nucleoid-associated proteins similar to those of transcription factors, consistent with an open and transcriptionally active genome. Taken together, these data suggest that histone in H. salinarum plays a minor role in DNA compaction but important roles in growth-phase-dependent gene expression and regulation of cell shape. Histone function in haloarchaea more closely resembles a regulator of gene expression than a chromatin-organizing protein like canonical eukaryotic histone.

Histones comprise the major protein component of eukaryotic chromatin and are required for both genome packaging and global regulation of expression. The current paradigm maintains that archaea whose genes encode histone also use these proteins to package DNA. In contrast, here we demonstrate that the sole histone encoded in the genome of the salt-adapted archaeon Halobacterium salinarum is both unessential and unlikely to be involved in DNA compaction despite conservation of residues important for eukaryotic histones. Rather, H. salinarum histone is required for global regulation of gene expression and cell shape. These data are consistent with the hypothesis that H. salinarum histone, strongly conserved across all other known salt-adapted archaea, serves a novel role in gene regulation and cell shape maintenance. Given that archaea possess the ancestral form of eukaryotic histone, this study has important implications for understanding the evolution of histone function.

Horizontal gene transfer and genome evolution in Methanosarcina

Background

Genomes of Methanosarcina spp. are among the largest archaeal genomes. One suggested reason for that is massive horizontal gene transfer (HGT) from bacteria. Genes of bacterial origin may be involved in the central metabolism and solute transport, in particular sugar synthesis, sulfur metabolism, phosphate metabolism, DNA repair, transport of small molecules etc. Horizontally transferred (HT) genes are considered to play the key role in the ability of Methanosarcina spp. to inhabit diverse environments. At the moment, genomes of three Methanosarcina spp. have been sequenced, and while these genomes vary in length and number of protein-coding genes, they all have been shown to accumulate HT genes. However, previous estimates had been made when fewer archaeal genomes were known. Moreover, several Methanosarcinaceae genomes from other genera have been sequenced recently. Here, we revise the census of genes of bacterial origin in Methanosarcinaceae.

Results

About 5 % of Methanosarcina genes have been shown to be horizontally transferred from various bacterial groups to the last common ancestor either of Methanosarcinaceae, or Methanosarcina, or later in the evolution. Simulation of the composition of the NCBI protein non-redundant database for different years demonstrates that the estimates of the HGT rate have decreased drastically since 2002, the year of publication of the first Methanosarcina genome.

The phylogenetic distribution of HT gene donors is non-uniform. Most HT genes were transferred from Firmicutes and Proteobacteria, while no HGT events from Actinobacteria to the common ancestor of Methanosarcinaceae were found. About 50 % of HT genes are involved in metabolism. Horizontal transfer of transcription factors is not common, while 46 % of horizontally transferred genes have demonstrated differential expression in a variety of conditions. HGT of complete operons is relatively infrequent and half of HT genes do not belong to operons.

Conclusions

While genes of bacterial origin are still more frequent in Methanosarcinaceae than in other Archaea, most HGT events described earlier as Methanosarcina-specific seem to have occurred before the divergence of Methanosarcinaceae. Genes horizontally transferred from bacteria to archaea neither tend to be transferred with their regulators, nor in long operons.

Electronic supplementary material

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

The transcriptional activator NrpA is crucial for inducing nitrogen fixation in Methanosarcina mazei Gö1 under nitrogen-limited conditions

With the aim of unraveling their potential involvement in the regulation of nitrogen metabolism in Methanosarcina mazei strain Gö1, we characterized five genes that are differentially transcribed in response to changing nitrogen availability and encoding putative transcriptional regulators. Study of the respective mutant strains under nitrogen-limited conditions revealed a growth delay for M. mazei MM0444::pac and MM1708::pac, and strongly reduced diazotrophic growth for MM0872::pac, whereas the absence of MM2441 or MM2525 did not affect growth behaviour. Transcriptome analyses further demonstrated that only MM1708 - encoding a CxxCG zinc finger protein - plays a regulatory role in nitrogen metabolism, most likely by specifically enhancing transcription of the N2 fixation (nif) operon under nitrogen-limited conditions. In agreement with this, a palindromic binding motif was predicted in silico in the nifH promoter region, nine nucleotides upstream of the BRE box, and confirmed to bind purified maltose-binding protein-MM1708 by electromobility shift assays. As MM1708 itself is under the control of the global nitrogen repressor NrpR, this adds a secondary level to the transcriptional regulation of the nif genes, and is most likely crucial for maximal nif induction under nitrogen-limited conditions. This is in accordance with the finding that protein expression of NifH is highly reduced in the absence of MM1708 under nitrogen-limited conditions. On the basis of our findings, we hypothesize that, in M. mazei, nitrogen fixation is controlled by a hierarchical network of two transcriptional regulators, the global nitrogen repressor NrpR, and the newly identified activator NrpA (MM1708), thereby providing tight control of N2 fixation.

Two CRISPR-Cas systems in Methanosarcina mazei strain Gö1 display common processing features despite belonging to different types I and III

The clustered regularly interspaced short palindromic repeats (CRISPR) system represents a highly adaptive and heritable defense system against foreign nucleic acids in bacteria and archaea. We analyzed the two CRISPR-Cas systems in Methanosarcina mazei strain Gö1. Although belonging to different subtypes (I-B and III-B), the leaders and repeats of both loci are nearly identical. Also, despite many point mutations in each array, a common hairpin motif was identified in the repeats by a bioinformatics analysis and in vitro structural probing. The expression and maturation of CRISPR-derived RNAs (crRNAs) were studied in vitro and in vivo. Both respective potential Cas6b-type endonucleases were purified and their activity tested in vitro. Each protein showed significant activity and could cleave both repeats at the same processing site. Cas6b of subtype III-B, however, was significantly more efficient in its cleavage activity compared with Cas6b of subtype I-B. Northern blot and differential RNAseq analyses were performed to investigate in vivo transcription and maturation of crRNAs, revealing generally very low expression of both systems, whereas significant induction at high NaCl concentrations was observed. crRNAs derived proximal to the leader were generally more abundant than distal ones and in vivo processing sites were clarified for both loci, confirming the previously well-established 8 nt 5' repeat tags. The 3'-ends were more diverse, but generally ended in a prefix of the following repeat sequence (3'-tag). The analysis further revealed a 5'-hydroxy and 3'-phosphate termini architecture of small crRNAs specific for cleavage products of Cas6 endonucleases from type I-E and I-F and type III-B.

An archaeal histone is required for transformation of Thermococcus kodakarensis

Archaeal histones wrap DNA into complexes, designated archaeal nucleosomes, that resemble the tetrasome core of a eukaryotic nucleosome. Therefore, all DNA interactions in vivo in Thermococcus kodakarensis, the most genetically versatile model species for archaeal research, must occur in the context of a histone-bound genome. Here we report the construction and properties of T. kodakarensis strains that have TK1413 or TK2289 deleted, the genes that encode HTkA and HTkB, respectively, the two archaeal histones present in this archaeon. All attempts to generate a strain with both TK1413 and TK2289 deleted were unsuccessful, arguing that a histone-mediated event(s) in T. kodakarensis is essential. The HTkA and HTkB amino acid sequences are 84% identical (56 of 67 residues) and 94% similar (63 of 67 residues), but despite this homology and their apparent redundancy in terms of supporting viability, the absence of HTkA and HTkB resulted in differences in growth and in quantitative and qualitative differences in genome transcription. A most surprising result was that the deletion of TK1413 (ΔhtkA) resulted in a T. kodakarensis strain that was no longer amenable to transformation, whereas the deletion of TK2289 (ΔhtkB) had no detrimental effects on transformation. Potential roles for the archaeal histones in regulating gene expression and for HTkA in DNA uptake and recombination are discussed.

Overview of the genetic tools in the Archaea

This section provides an overview of the genetic systems developed in the Archaea. Genetic manipulation is possible in many members of the halophiles, methanogens, Sulfolobus, and Thermococcales. We describe the selection/counterselection principles utilized in each of these groups, which consist of antibiotics and their resistance markers, and auxotrophic host strains and complementary markers. The latter strategy utilizes techniques similar to those developed in yeast. However, Archaea are resistant to many of the antibiotics routinely used for selection in the Bacteria, and a number of strategies specific to the Archaea have been developed. In addition, examples utilizing the genetic systems developed for each group will be briefly described.

An archaeal sRNA targeting cis- and trans-encoded mRNAs via two distinct domains

We report on the characterization and target analysis of the small (s)RNA₁₆₂ in the methanoarchaeon Methanosarcina mazei. Using a combination of genetic approaches, transcriptome analysis and computational predictions, the bicistronic MM2441-MM2440 mRNA encoding the transcription factor MM2441 and a protein of unknown function was identified as a potential target of this sRNA, which due to processing accumulates as three stabile 5′ fragments in late exponential growth. Mobility shift assays using various mutants verified that the non-structured single-stranded linker region of sRNA₁₆₂ (SLR) base-pairs with the MM2440-MM2441 mRNA internally, thereby masking the predicted ribosome binding site of MM2441. This most likely leads to translational repression of the second cistron resulting in dis-coordinated operon expression. Analysis of mutant RNAs in vivo confirmed that the SLR of sRNA₁₆₂ is crucial for target interactions. Furthermore, our results indicate that sRNA₁₆₂-controlled MM2441 is involved in regulating the metabolic switch between the carbon sources methanol and methylamine. Moreover, biochemical studies demonstrated that the 5′ end of sRNA₁₆₂ targets the 5′-untranslated region of the cis-encoded MM2442 mRNA. Overall, this first study of archaeal sRNA/mRNA-target interactions unraveled that sRNA₁₆₂ acts as an antisense (as)RNA on cis- and trans-encoded mRNAs via two distinct domains, indicating that cis-encoded asRNAs can have larger target regulons than previously anticipated.

Archaeal chromatin proteins

Archaea, along with Bacteria and Eukarya, are the three domains of life. In all living cells, chromatin proteins serve a crucial role in maintaining the integrity of the structure and function of the genome. An array of small, abundant and basic DNA-binding proteins, considered candidates for chromatin proteins, has been isolated from the Euryarchaeota and the Crenarchaeota, the two major phyla in Archaea. While most euryarchaea encode proteins resembling eukaryotic histones, crenarchaea appear to synthesize a number of unique DNA-binding proteins likely involved in chromosomal organization. Several of these proteins (e.g., archaeal histones, Sac10b homologs, Sul7d, Cren7, CC1, etc.) have been extensively studied. However, whether they are chromatin proteins and how they function in vivo remain to be fully understood. Future investigation of archaeal chromatin proteins will lead to a better understanding of chromosomal organization and gene expression in Archaea and provide valuable information on the evolution of DNA packaging in cellular life.

The Histone Database: an integrated resource for histones and histone fold-containing proteins

Eukaryotic chromatin is composed of DNA and protein components-core histones-that act to compactly pack the DNA into nucleosomes, the fundamental building blocks of chromatin. These nucleosomes are connected to adjacent nucleosomes by linker histones. Nucleosomes are highly dynamic and, through various core histone post-translational modifications and incorporation of diverse histone variants, can serve as epigenetic marks to control processes such as gene expression and recombination. The Histone Sequence Database is a curated collection of sequences and structures of histones and non-histone proteins containing histone folds, assembled from major public databases. Here, we report a substantial increase in the number of sequences and taxonomic coverage for histone and histone fold-containing proteins available in the database. Additionally, the database now contains an expanded dataset that includes archaeal histone sequences. The database also provides comprehensive multiple sequence alignments for each of the four core histones (H2A, H2B, H3 and H4), the linker histones (H1/H5) and the archaeal histones. The database also includes current information on solved histone fold-containing structures. The Histone Sequence Database is an inclusive resource for the analysis of chromatin structure and function focused on histones and histone fold-containing proteins.

Establishing a markerless genetic exchange system for Methanosarcina mazei strain Gö1 for constructing chromosomal mutants of small RNA genes

A markerless genetic exchange system was successfully established in Methanosarcina mazei strain Gö1 using the hpt gene coding for hypoxanthine phosphoribosyltransferase. First, a chromosomal deletion mutant of the hpt gene was generated conferring resistance to the purine analog 8-aza-2,6-diaminopurine (8-ADP). The nonreplicating allelic exchange vector (pRS345) carrying the pac-resistance cassette for direct selection of chromosomal integration, and the hpt gene for counterselection was introduced into this strain. By a pop-in and ultimately pop-out event of the plasmid from the chromosome, allelic exchange is enabled. Using this system, we successfully generated a M. mazei deletion mutant of the gene encoding the regulatory non-coding RNA sRNA₁₅₄. Characterizing M. mazeiΔsRNA₁₅₄ under nitrogen limiting conditions demonstrated differential expression of at least three cytoplasmic proteins and reduced growth strongly arguing for a prominent role of sRNA₁₅₄ in regulation of nitrogen fixation by posttranscriptional regulation.

Refined solution structure and backbone dynamics of the archaeal MC1 protein

The 3D structure of methanogen chromosomal protein 1 (MC1), determined with heteronuclear NMR methods, agrees with its function in terms of the shape and nature of the binding surface, whereas the 3D structure determined with homonuclear NMR does not. The structure features five loops, which show a large distribution in the ensemble of 3D structures. Evidence for the fact that this distribution signifies internal mobility on the nanosecond time scale was provided by using (15)N-relaxation and molecular dynamics simulations. Structural variations of the arm (11 residues) induced large shape anisotropy variations on the nanosecond time scale that ruled out the use of the model-free formalism to analyze the relaxation data. The backbone dynamics analysis of MC1 was achieved by comparison with 20 ns molecular dynamics trajectories. Two β-bulges showed that hydrogen bond formation correlated with ϕ and ψ dihedral angle transitions. These jumps were observed on the nanosecond time scale, in agreement with a large decrease in (15)N-NOE for Gly17 and Ile89. One water molecule bridging NH(Glu87) and CO(Val57) through hydrogen bonding contributed to these dynamics. Nanosecond slow motions observed in loops LP3 (35-42) and LP5 (67-77) reflected the lack of stable hydrogen bonds, whereas the other loops, LP1 (10-14), LP2 (22-24), and LP4 (50-53), were stabilized by several hydrogen bonds. Dynamics are often directly related to function. Our data strongly suggest that residues belonging to the flexible regions of MC1 could be involved in the interaction with DNA.

In vivo effects of NbSiR silencing on chloroplast development in Nicotiana benthamiana

Sulfite reductase (SiR) performs dual functions, acting as a sulfur assimilation enzyme and as a chloroplast (cp-) nucleoid binding protein. In this study, we examined the in vivo effects of SiR deficiency on chloroplast development in Nicotiana benthamiana. Virus-induced gene silencing of NbSiR resulted in leaf yellowing and growth retardation phenotypes, which were not rescued by cysteine supplementation. NbSiR:GFP fusion protein was targeted to chloroplasts and colocalized with cp-nucleoids. Recombinant full-length NbSiR protein and the C-terminal half of NbSiR possessed cp-DNA compaction activities in vitro, and expression of full-length NbSiR in E. coli caused condensation of genomic DNA. NbSiR silencing differentially affected expression of plastid-encoded genes, inhibiting expression of several genes more severely than others. In the later stages, depletion of NbSiR resulted in chloroplast ablation. In NbSiR-silenced plants, enlarged cp-nucleoids containing an increased amount of cp-DNA were observed in the middle of the abnormal chloroplasts, and the cp-DNAs were predominantly of subgenomic sizes based on pulse field gel electrophoresis. The abnormal chloroplasts developed prolamellar body-like cubic lipid structures in the light without accumulating NADPH:protochlorophyllide oxidoreductase proteins. Our results suggest that NbSiR plays a role in cp-nucleoid metabolism, plastid gene expression, and thylakoid membrane development.

Archaeal chromatin proteins histone HMtB and Alba have lost DNA-binding ability in laboratory strains of Methanothermobacter thermautotrophicus

Alignments of the sequences of the all members of the archaeal histone and Alba1 families of chromatin proteins identified isoleucine residues, I19 in HMtB and I39 in MtAlba, in Methanothermobacter thermautotrophicus, at locations predicted to be directly involved in DNA binding. In all other HMfB family members, residue 19 is an arginine (R19), and either arginine or lysine is present in almost all other Alba1 family members at the structural site equivalent to I39 in MtAlba. Electrophoretic mobility shift assays revealed that recombinant HMtB and MtAlba do not bind DNA, but variants constructed with R19 and R39, respectively, bound DNA; and whereas MtAlba(I19) did not bind RNA, MtAlba(R19) bound both single stranded RNA and tRNA. Amplification and sequencing of MT0254 (encodes HMtB) and MT1483 (encodes MtAlba) from several Methanothermobacter thermautotrophicus lineages has revealed that HMtB and MtAlba had arginine residues at positions 19 and 39, respectively, in the original isolate and that spontaneous mutations must have occurred, and been fixed, in some laboratory lineages that now have HMtB(I19) and MtAlba(I39). The retention of these variants suggests some continuing functions and fusion of the HMtB(I19) sequence to HMtA2 resulted in a protein that folds to form a histone fold heterodimer that binds and compacts DNA. The loss of DNA binding by HMtB(I19) does not therefore prevent HMtB from participating in DNA interactions as one partner of an archaeal histone heterodimer.