Genome-wide binding analysis of the transcriptional regulator TrmBL1 in Pyrococcus furiosus

A conserved transcription factor controls gluconeogenesis via distinct targets in hypersaline-adapted archaea with diverse metabolic capabilities

Timely regulation of carbon metabolic pathways is essential for cellular processes and to prevent futile cycling of intracellular metabolites. In Halobacterium salinarum, a hypersaline adapted archaeon, a sugar-sensing TrmB family protein controls gluconeogenesis and other biosynthetic pathways. Notably, Hbt. salinarum does not utilize carbohydrates for energy, uncommon among Haloarchaea. We characterized a TrmB-family transcriptional regulator in a saccharolytic generalist, Haloarcula hispanica, to investigate whether the targets and function of TrmB, or its regulon, is conserved in related species with distinct metabolic capabilities. In Har. hispanica, TrmB binds to 15 sites in the genome and induces the expression of genes primarily involved in gluconeogenesis and tryptophan biosynthesis. An important regulatory control point in Hbt. salinarum, activation of ppsA and repression of pykA, is absent in Har. hispanica. Contrary to its role in Hbt. salinarum and saccharolytic hyperthermophiles, TrmB does not act as a global regulator: it does not directly repress the expression of glycolytic enzymes, peripheral pathways such as cofactor biosynthesis, or catabolism of other carbon sources in Har. hispanica. Cumulatively, these findings suggest rewiring of the TrmB regulon alongside metabolic network evolution in Haloarchaea.

Uncovering the temporal dynamics and regulatory networks of thermal stress response in a hyperthermophile using transcriptomics and proteomics

IMPORTANCE

Extreme environments provide unique challenges for life, and the study of extremophiles can shed light on the mechanisms of adaptation to such conditions. Pyrococcus furiosus, a hyperthermophilic archaeon, is a model organism for studying thermal stress response mechanisms. In this study, we used an integrated analysis of RNA-sequencing and mass spectrometry data to investigate the transcriptomic and proteomic responses of P. furiosus to heat and cold shock stress and recovery. Our results reveal the rapid and dynamic changes in gene and protein expression patterns associated with these stress responses, as well as the coordinated regulation of different gene sets in response to different stressors. These findings provide valuable insights into the molecular adaptations that facilitate life in extreme environments and advance our understanding of stress response mechanisms in hyperthermophilic archaea.

The archaeal Lsm protein from Pyrococcus furiosus binds co-transcriptionally to poly(U)-rich target RNAs

Posttranscriptional processes in Bacteria include the association of small regulatory RNAs (sRNA) with a target mRNA. The sRNA/mRNA annealing process is often mediated by an RNA chaperone called Hfq. The functional role of bacterial and eukaryotic Lsm proteins is partially understood, whereas knowledge about archaeal Lsm proteins is scarce. Here, we used the genetically tractable archaeal hyperthermophile Pyrococcus furiosus to identify the protein interaction partners of the archaeal Sm-like proteins (PfuSmAP1) using mass spectrometry and performed a transcriptome-wide binding site analysis of PfuSmAP1. Most of the protein interaction partners we found are part of the RNA homoeostasis network in Archaea including ribosomal proteins, the exosome, RNA-modifying enzymes, but also RNA polymerase subunits, and transcription factors. We show that PfuSmAP1 preferentially binds messenger RNAs and antisense RNAs recognizing a gapped poly(U) sequence with high affinity. Furthermore, we found that SmAP1 co-transcriptionally associates with target RNAs. Our study reveals that in contrast to bacterial Hfq, PfuSmAP1 does not affect the transcriptional activity or the pausing behaviour of archaeal RNA polymerases. We propose that PfuSmAP1 recruits antisense RNAs to target mRNAs and thereby executes its putative regulatory function on the posttranscriptional level.

Optimizing Strategies for Bio-Based Ethanol Production Using Genome-Scale Metabolic Modeling of the Hyperthermophilic Archaeon, Pyrococcus furiosus

A genome-scale metabolic model, encompassing a total of 623 genes, 727 reactions, and 865 metabolites, was developed for Pyrococcus furiosus, an archaeon that grows optimally at 100°C by carbohydrate and peptide fermentation. The model uses subsystem-based genome annotation, along with extensive manual curation of 237 gene-reaction associations including those involved in central carbon metabolism, amino acid metabolism, and energy metabolism. The redox and energy balance of P. furiosus was investigated through random sampling of flux distributions in the model during growth on disaccharides. The core energy balance of the model was shown to depend on high acetate production and the coupling of a sodium-dependent ATP synthase and membrane-bound hydrogenase, which generates a sodium gradient in a ferredoxin-dependent manner, aligning with existing understanding of P. furiosus metabolism. The model was utilized to inform genetic engineering designs that favor the production of ethanol over acetate by implementing an NADPH and CO-dependent energy economy. The P. furiosus model is a powerful tool for understanding the relationship between generation of end products and redox/energy balance at a systems-level that will aid in the design of optimal engineering strategies for production of bio-based chemicals and fuels. IMPORTANCE The bio-based production of organic chemicals provides a sustainable alternative to fossil-based production in the face of today's climate challenges. In this work, we present a genome-scale metabolic reconstruction of Pyrococcus furiosus, a well-established platform organism that has been engineered to produce a variety of chemicals and fuels. The metabolic model was used to design optimal engineering strategies to produce ethanol. The redox and energy balance of P. furiosus was examined in detail, which provided useful insights that will guide future engineering designs.

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

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

ChIP-Seq Occupancy Mapping of the Archaeal Transcription Machinery

Genome-wide occupancy studies for RNA polymerases and their basal transcription factors deliver information about transcription dynamics and the recruitment of transcription elongation and termination factors in eukaryotes and prokaryotes. The primary method to determine genome-wide occupancies is chromatin immunoprecipitation combined with deep sequencing (ChIP-seq). Archaea possess a transcription machinery that is evolutionarily closer related to its eukaryotic counterpart but it operates in a prokaryotic cellular context. Studies on archaeal transcription brought insight into the evolution of transcription machineries and the universality of transcription mechanisms. Because of the limited resolution of ChIP-seq, the close spacing of promoters and transcription units found in archaeal genomes pose a challenge for ChIP-seq and the ensuing data analysis. The extreme growth temperature of many established archaeal model organisms necessitates further adaptations. This chapter describes a version of ChIP-seq adapted for the basal transcription machinery of thermophilic archaea and some modifications to the data analysis.

CopR, a Global Regulator of Transcription to Maintain Copper Homeostasis in Pyrococcus furiosus

Although copper is in many cases an essential micronutrient for cellular life, higher concentrations are toxic. Therefore, all living cells have developed strategies to maintain copper homeostasis. In this manuscript, we have analyzed the transcriptome-wide response of Pyrococcus furiosus to increased copper concentrations and described the essential role of the putative copper-sensing metalloregulator CopR in the detoxification process. To this end, we employed biochemical and biophysical methods to characterize the role of CopR. Additionally, a copR knockout strain revealed an amplified sensitivity in comparison to the parental strain towards increased copper levels, which designates an essential role of CopR for copper homeostasis. To learn more about the CopR-regulated gene network, we performed differential gene expression and ChIP-seq analysis under normal and 20 μM copper-shock conditions. By integrating the transcriptome and genome-wide binding data, we found that CopR binds to the upstream regions of many copper-induced genes. Negative-stain transmission electron microscopy and 2D class averaging revealed an octameric assembly formed from a tetramer of dimers for CopR, similar to published crystal structures from the Lrp family. In conclusion, we propose a model for CopR-regulated transcription and highlight the regulatory network that enables Pyrococcus to respond to increased copper concentrations.

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.

Next Generation DNA-Seq and Differential RNA-Seq Allow Re-annotation of the Pyrococcus furiosus DSM 3638 Genome and Provide Insights Into Archaeal Antisense Transcription

Pyrococcus furiosus DSM 3638 is a model organism for hyperthermophilic archaea with an optimal growth temperature near 100°C. The genome was sequenced about 18 years ago. However, some publications suggest that in contrast to other Pyrococcus species, the genome of P. furiosus DSM 3638 is prone to genomic rearrangements. Therefore, we re-sequenced the genome using third generation sequencing techniques. The new de novo assembled genome is 1,889,914 bp in size and exhibits high sequence identity to the published sequence. However, two major deviations were detected: (1) The genome is 18,342 bp smaller than the NCBI reference genome due to a recently described deletion. (2) The region between PF0349 and PF0388 is inverted most likely due an assembly problem for the original sequence. In addition, numerous minor variations, ranging from single nucleotide exchanges, deletions or insertions were identified. The total number of insertion sequence (IS) elements is also reduced from 30 to 24 in the new sequence. Re-sequencing of a 2-year-old “lab culture” using Nanopore sequencing confirmed the overall stability of the P. furiosus DSM 3638 genome even under normal lab conditions without taking any special care. To improve genome annotation, the updated DNA sequence was combined with an RNA sequencing approach. Here, RNAs from eight different growth conditions were pooled to increase the number of detected transcripts. Furthermore, a differential RNA-Seq approach was employed for the identification of transcription start sites (TSSs). In total, 2515 TSSs were detected and classified into 834 primary (pTSS), 797 antisense (aTSS), 739 internal and 145 secondary TSSs. Our analysis of the upstream regions revealed a well conserved archaeal promoter structure. Interrogation of the distances between pTSSs and aTSSs revealed a significant number of antisense transcripts, which are a result of bidirectional transcription from the same TATA box. This mechanism of antisense transcript production could be further confirmed by in vitro transcription experiments. We assume that bidirectional transcription gives rise to non-functional antisense RNAs and that this is a widespread phenomenon in archaea due to the architecture of the TATA element and the symmetric structure of the TATA-binding protein.

Conserved principles of transcriptional networks controlling metabolic flexibility in archaea

Gene regulation is intimately connected with metabolism, enabling the appropriate timing and tuning of biochemical pathways to substrate availability. In microorganisms, such as archaea and bacteria, transcription factors (TFs) often directly sense external cues such as nutrient substrates, metabolic intermediates, or redox status to regulate gene expression. Intense recent interest has characterized the functions of a large number of such regulatory TFs in archaea, which regulate a diverse array of unique archaeal metabolic capabilities. However, it remains unclear how the co-ordinated activity of the interconnected metabolic and transcription networks produces the dynamic flexibility so frequently observed in archaeal cells as they respond to energy limitation and intermittent substrate availability. In this review, we communicate the current state of the art regarding these archaeal networks and their dynamic properties. We compare the topology of these archaeal networks to those known for bacteria to highlight conserved and unique aspects. We present a new computational model for an exemplar archaeal network, aiming to lay the groundwork toward understanding general principles that unify the dynamic function of integrated metabolic-transcription networks across archaea and bacteria.

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.

Chromatin Immunoprecipitation Assay in the Hyperthermoacidophilic Crenarchaeon, Sulfolobus acidocaldarius

Chromatin immunoprecipitation (ChIP) is a powerful method used for identifying genome-wide DNA-protein interactions in vivo. A large number of essential intracellular processes such as DNA replication, transcription regulation, chromatin stability, and others are all dependent on protein interactions with DNA. The DNA fragments enriched from the ChIP assay are analyzed by downstream applications, for example, microarray hybridization (ChIP-chip), quantitative PCR (ChIP-qPCR), or deep sequencing (ChIP-seq). This chapter presents a stepwise protocol for ChIP performed in hyperthermophilic archaea that we have successfully used in the hyperthermoacidophilic crenarchaeon Sulfolobus acidocaldarius.

The Transcriptional Regulator TFB-RF1 Activates Transcription of a Putative ABC Transporter in Pyrococcus furiosus

Transcription factor B recruiting factor 1 (TFB-RF1; PF1088) is a transcription regulator which activates transcription on archaeal promoters containing weak TFB recognition elements (BRE) by recruiting TFB to the promoter. The mechanism of activation is described in detail, but nothing is known about the biological function of this protein in Pyrococcus furiosus. The protein is located in an operon structure together with the hypothetical gene pf1089 and western blot as well as end-point RT-PCR experiments revealed an extremely low expression rate of both proteins. Furthermore, conditions to induce the expression of the operon are not known. By introducing an additional copy of tfb-RF1 using a Pyrococcus shuttle vector we could circumvent the lacking expression of both proteins under standard growth conditions as indicated by western blot as well as end-point RT-PCR experiments. A ChIP-seq experiment revealed an additional binding site of TFB-RF1 in the upstream region of the pf1011/1012 operon, beside the expected target of the pf1089/tfb-RF1 region. This operon codes for a putative ABC transporter which is most-related to a multidrug export system and in vitro analysis using gel shift assays, DNase I footprinting and in vitro transcription confirmed the activator function of TFB-RF1 on the corresponding promoter. These findings are also in agreement with in vivo data, as RT-qPCR experiments also indicate transcriptional activation of both operons. Taken together, the overexpression strategy of tfb-RF1 enabled the identification of an additional operon of the TFB-RF1 regulon which indicates a transport-related function and provides a promising starting position to decipher the physiological function of the TFB-RF1 gene regulatory network in P. furiosus.

Structure and in situ organisation of the Pyrococcus furiosus archaellum machinery

The archaellum is the macromolecular machinery that Archaea use for propulsion or surface adhesion, enabling them to proliferate and invade new territories. The molecular composition of the archaellum and of the motor that drives it appears to be entirely distinct from that of the functionally equivalent bacterial flagellum and flagellar motor. Yet, the structure of the archaellum machinery is scarcely known. Using combined modes of electron cryo-microscopy (cryoEM), we have solved the structure of the Pyrococcus furiosus archaellum filament at 4.2 Å resolution and visualise the architecture and organisation of its motor complex in situ. This allows us to build a structural model combining the archaellum and its motor complex, paving the way to a molecular understanding of archaeal swimming motion.

A global analysis of transcription reveals two modes of Spt4/5 recruitment to archaeal RNA polymerase

The archaeal transcription apparatus is closely related to the eukaryotic RNA polymerase (RNAP) II system, while archaeal genomes are more similar to bacteria with densely packed genes organized in operons. This makes understanding transcription in archaea vital, both in terms of molecular mechanisms and evolution. Very little is known about how archaeal cells orchestrate transcription on a systems level. We have characterized the genome-wide occupancy of the Methanocaldococcus jannaschii transcription machinery and its transcriptome. Our data reveal how the TATA and BRE promoter elements facilitate recruitment of the essential initiation factors TATA-binding protein and transcription factor B, respectively, which in turn are responsible for the loading of RNAP into the transcription units. The occupancies of RNAP and Spt4/5 strongly correlate with each other and with RNA levels. Our results show that Spt4/5 is a general elongation factor in archaea as its presence on all genes matches RNAP. Spt4/5 is recruited proximal to the transcription start site on the majority of transcription units, while on a subset of genes, including rRNA and CRISPR loci, Spt4/5 is recruited to the transcription elongation complex during early elongation within 500 base pairs of the transcription start site and akin to its bacterial homologue NusG.

Global transcriptional regulator TrmB family members in prokaryotes

Members of the TrmB family act as global transcriptional regulators for the activation or repression of sugar ABC transporters and central sugar metabolic pathways, including glycolytic, gluconeogenic, and other metabolic pathways, and also as chromosomal stabilizers in archaea. As a relatively newly classified transcriptional regulator family, there is limited experimental evidence for their role in Thermococcales, halophilic archaeon Halobacterium salinarum NRC1, and crenarchaea Sulfolobus strains, despite being one of the extending protein families in archaea. Recently, the protein structures of Pyrococcus furiosus TrmB and TrmBL2 were solved, and the transcriptomic data uncovered by microarray and ChIP-Seq were published. In the present review, recent evidence of the functional roles of TrmB family members in archaea is explained and extended to bacteria.

The genome-scale DNA-binding profile of BarR, a β-alanine responsive transcription factor in the archaeon Sulfolobus acidocaldarius

Background

The Leucine-responsive Regulatory Protein (Lrp) family is a widespread family of regulatory transcription factors in prokaryotes. BarR is an Lrp-like transcription factor in the model archaeon Sulfolobus acidocaldarius that activates the expression of a β-alanine aminotransferase gene, which is involved in β-alanine degradation. In contrast to classical Lrp-like transcription factors, BarR is not responsive to any of the α-amino acids but interacts specifically with β-alanine. Besides the juxtaposed β-alanine aminotransferase gene, other regulatory targets of BarR have not yet been identified although β-alanine is the precursor of coenzyme A and thus an important central metabolite. The aim of this study is to extend the knowledge of the DNA-binding characteristics of BarR and of its corresponding regulon from a local to a genome-wide perspective.

Results

We characterized the genome-wide binding profile of BarR using chromatin immunoprecipation combined with high-throughput sequencing (ChIP-seq). This revealed 21 genomic binding loci. High-enrichment binding regions were validated to interact with purified BarR protein in vitro using electrophoretic mobility shift assays and almost all targets were also shown to harbour a conserved semi-palindromic binding motif. Only a small subset of enriched genomic sites are located in intergenic regions at a relative short distance to a promoter, and qRT-PCR analysis demonstrated that only one additional operon is under activation of BarR, namely the glutamine synthase operon. The latter is also a target of other Lrp-like transcription factors. Detailed inspection of the BarR ChIP-seq profile at the β-alanine aminotransferase promoter region in combination with binding motif predictions indicate that the operator structure is more complicated than previously anticipated, consisting of multiple (major and auxiliary) operators.

Conclusions

BarR has a limited regulon, and includes also glutamine synthase genes besides the previously characterized β-alanine aminotransferase. Regulation of glutamine synthase is suggestive of a link between β-alanine and α-amino acid metabolism in S. acidocaldarius. Furthermore, this work reveals that the BarR regulon overlaps with that of other Lrp-like regulators.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2890-0) contains supplementary material, which is available to authorized users.

TrmBL2 from Pyrococcus furiosus Interacts Both with Double-Stranded and Single-Stranded DNA

In many hyperthermophilic archaea the DNA binding protein TrmBL2 or one of its homologues is abundantly expressed. TrmBL2 is thought to play a significant role in modulating the chromatin architecture in combination with the archaeal histone proteins and Alba. However, its precise physiological role is poorly understood. It has been previously shown that upon binding TrmBL2 covers double-stranded DNA, which leads to the formation of a thick and fibrous filament. Here we investigated the filament formation process as well as the stabilization of DNA by TrmBL2 from Pyroccocus furiosus in detail. We used magnetic tweezers that allow to monitor changes of the DNA mechanical properties upon TrmBL2 binding on the single-molecule level. Extended filaments formed in a cooperative manner and were considerably stiffer than bare double-stranded DNA. Unlike Alba, TrmBL2 did not form DNA cross-bridges. The protein was found to bind double- and single-stranded DNA with similar affinities. In mechanical disruption experiments of DNA hairpins this led to stabilization of both, the double- (before disruption) and the single-stranded (after disruption) DNA forms. Combined, these findings suggest that the biological function of TrmBL2 is not limited to modulating genome architecture and acting as a global repressor but that the protein acts additionally as a stabilizer of DNA secondary structure.