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 degrees 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

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

Facing rapid fluctuations in their natural environment, extremophiles, like the hyperthermophilic archaeon Pyrococcus furiosus, exhibit remarkable adaptability to extreme conditions. However, our understanding of their dynamic cellular responses remains limited. This study integrates RNA-sequencing and mass spectrometry data, thereby elucidating transcriptomic and proteomic responses to heat and cold shock stress in P. furiosus. Our results reveal rapid and dynamic changes in gene and protein expression following these stress responses. Heat shock triggers extensive transcriptome reprogramming, orchestrated by the transcriptional regulator Phr, targeting a broader gene repertoire than previously demonstrated. For heat shock signature genes, RNA levels swiftly return to baseline upon recovery, while protein levels remain persistently upregulated, reflecting a rapid but sustained response. Intriguingly, cold shock at 4°C elicits distinct short- and long-term responses at both RNA and protein levels. Cluster analysis identified gene sets with either congruent or contrasting trends in RNA and protein changes, representing well-separated arCOG groups tailored to their individual cellular responses. Particularly, upregulation of ribosomal proteins and significant enrichment of 5′-leadered sequences in cold-shock responsive genes suggest that translation regulation is important during cold shock adaption. Further investigating transcriptomic features, we reveal that thermal stress genes are equipped with basal sequence elements, such as strong promoter and poly(U)-terminators, facilitating a regulated response of the respective transcription units. Our study provides a comprehensive overview of the cellular response to temperature stress, advancing our understanding of stress response mechanisms in hyperthermophilic archaea and providing valuable insights into the molecular adaptations that facilitate life in extreme environments

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

Genetic engineering of Pyrococcus furiosus to use chitin as a carbon source

BACKGROUND: Bioinformatic analysis of the genes coding for the chitinase in Pyrococcus furiosus and Thermococcus kodakarensis revealed that most likely a one nucleotide insertion in Pyrococcus caused a frame shift in the chitinase gene. This splits the enzyme into two separate genes, PF1233 and PF1234, in comparison to Thermococcus kodakarensis. Furthermore, our attempts to grow the wild type strain of Pyrococcus furiosus on chitin were negative. From these data we assume that Pyrococcus furiosus is most likely unable to use chitin as a carbon source. The aim of this study was to analyze in vivo if the one nucleotide insertion is responsible for the inability to grow on chitin, using a recently described genetic system for Pyrococcus furiosus. RESULTS: A marker-less genetic system for Pyrococcus furiosus was developed using simvastatin for positive selection and 6-methylpurine for negative selection. Resistance against simvastatin was achieved by overexpression of the hydroxymethylglutaryl coenzyme A reductase gene. For the resistance to 6-methylpurine the hypoxanthine-guanine phosphoribosyltransferase gene was deleted. This system was used to delete the additional nucleotide at position 1006 in PF1234. The resulting chitinase in the mutant strain was a single subunit enzyme and aligns perfectly to the enzyme from Thermococcus kodakarensis. A detailed analysis of the wild type and the mutant using counted cell numbers as well as ATP and acetate production as growth indicators revealed that only the mutant is able to use chitin as a carbon source. An additional mutant strain containing a reduced chitinase version containing just one catalytic and one chitin-binding domain showed diminished growth on chitin in comparison to the mutant containing the single large enzyme. CONCLUSIONS: Wild type Pyrococcus furiosus is most likely unable to grow on chitin in the natural biotope due to a nucleotide insertion which separates the chitinase gene into two ORFs, whereas a genetically engineered strain with the deleted nucleotide is able to grow on chitin. The overall high sequence identity of the two chitinases between P. furiosus and T. kodakarensis indicates that this mutation occurred very recently or there is still some kind of selection pressure for a functional enzyme using programmed +/-1 frameshifting

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

Background: Several in vitro studies document the function of the transcriptional regulator TrmBL1 of Pyrococcus furiosus. These data indicate that the protein can act as repressor or activator and is mainly involved in transcriptional control of sugar uptake and in the switch between glycolysis and gluconeogenesis. The aim of this study was to complement the in vitro data with an in vivo analysis using ChIP-seq to explore the genome-wide binding profile of TrmBL1 under glycolytic and gluconeogenic growth conditions. Results: The ChIP-seq analysis revealed under gluconeogenic growth conditions 28 TrmBL1 binding sites where the TGM is located upstream of coding regions and no binding sites under glycolytic conditions. The experimental confirmation of the binding sites using qPCR, EMSA, DNase I footprinting and in vitro transcription experiments validated the in vivo identified TrmBL1 binding sites. Furthermore, this study provides evidence that TrmBL1 is also involved in transcriptional regulation of additional cellular processes e. g. amino acid metabolism, transcriptional control or metabolic pathways. In the initial setup we were interested to include the binding analysis of TrmB, an additional member of the TrmB family, but western blot experiments and the ChIP-seq data indicated that the corresponding gene is deleted in our Pyrococcus strain. A detailed analysis of a new type strain demonstrated that a 16 kb fragment containing the trmb gene is almost completely deleted after the first re-cultivation. Conclusions: The identified binding sites in the P. furiosus genome classified TrmBL1 as a more global regulator as hitherto known. Furthermore, the high resolution of the mapped binding positions enabled reliable predictions, if TrmBL1 activates (binding site upstream of the promoter) or represses transcription (binding site downstream) of the corresponding genes

Structural Insights into Nonspecific Binding of DNA by TrmBL2, an Archaeal Chromatin Protein

The crystal structure of TrmBL2 from the archaeon Pyrococcus furiosus shows an association of two pseudosymmetric dimers. The dimers follow the prototypical design of known bacterial repressors with two helix–turn–helix (HTH) domains binding to successive major grooves of the DNA. However, in TrmBL2, the two dimers are arranged at a mutual displacement of approximately 2 bp so that they associate with the DNA along the double-helical axis at an angle of approximately 80°.While the deoxyribose phosphate groups of the double-stranded DNA (dsDNA) used for co-crystallization are clearly seen in the electron density map, most of the nucleobases are averaged out. Refinement required to assume a superposition of at least three mutually displaced dsDNAs. The HTH domains interact primarily with the deoxyribose phosphate groups and polar interactions with the nucleobases are almost absent.This hitherto unseen mode of DNA binding by TrmBL2 seems to arise from nonoptimal protein–DNA contacts made by its four HTH domains resulting in a low-affinity, nonspecific binding to DNA

Shuttle Vector-Based Transformation System for Pyrococcus furiosus▿

Pyrococcus furiosus is a model organism for analyses of molecular biology and biochemistry of archaea, but so far no useful genetic tools for this species have been described. We report here a genetic transformation system for P. furiosus based on the shuttle vector system pYS2 from Pyrococcus abyssi. In the redesigned vector, the pyrE gene from Sulfolobus was replaced as a selectable marker by the 3-hydroxy-3-methylglutaryl coenzyme A reductase gene (HMG-CoA) conferring resistance of transformants to the antibiotic simvastatin. Use of this modified plasmid resulted in the overexpression of the HMG-CoA reductase in P. furiosus, allowing the selection of strains by growth in the presence of simvastatin. The modified shuttle vector replicated in P. furiosus, but the copy number was only one to two per chromosome. This system was used for overexpression of His6-tagged subunit D of the RNA polymerase (RNAP) in Pyrococcus cells. Functional RNAP was purified from transformed cells in two steps by Ni-NTA and gel filtration chromatography. Our data provide evidence that expression of transformed genes can be controlled from a regulated gluconeogenetic promoter

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

Additional file 6: of Genome-wide binding analysis of the transcriptional regulator TrmBL1 in Pyrococcus furiosus

Gene-annotation enrichment analysis. (XLSX 13 kb