Parallel evolution of transcriptome architecture during genome reorganization

A comprehensive spectral assay library to quantify the Halobacterium salinarum NRC-1 proteome by DIA/SWATH-MS

Data-Independent Acquisition (DIA) is a mass spectrometry-based method to reliably identify and reproducibly quantify large fractions of a target proteome. The peptide-centric data analysis strategy employed in DIA requires a priori generated spectral assay libraries. Such assay libraries allow to extract quantitative data in a targeted approach and have been generated for human, mouse, zebrafish, E. coli and few other organisms. However, a spectral assay library for the extreme halophilic archaeon Halobacterium salinarum NRC-1, a model organism that contributed to several notable discoveries, is not publicly available yet. Here, we report a comprehensive spectral assay library to measure 2,563 of 2,646 annotated H. salinarum NRC-1 proteins. We demonstrate the utility of this library by measuring global protein abundances over time under standard growth conditions. The H. salinarum NRC-1 library includes 21,074 distinct peptides representing 97% of the predicted proteome and provides a new, valuable resource to confidently measure and quantify any protein of this archaeon. Data and spectral assay libraries are available via ProteomeXchange (PXD042770, PXD042774) and SWATHAtlas (SAL00312-SAL00319).

Tuning Gene Expression by Phosphate in the Methanogenic Archaeon Methanococcus maripaludis

Methanococcus maripaludis is a rapidly growing, hydrogenotrophic, and genetically tractable methanogen with unique capabilities to convert formate and CO₂ to CH₄. The existence of genome-scale metabolic models and an established, robust system for both large-scale and continuous cultivation make it amenable for industrial applications. However, the lack of molecular tools for differential gene expression has hindered its application as a microbial cell factory to produce biocatalysts and biochemicals. In this study, a library of differentially regulated promoters was designed and characterized based on the pst promoter, which responds to the inorganic phosphate concentration in the growth medium. Gene expression increases by 4- to 6-fold when the medium phosphate drops to growth-limiting concentrations. Hence, this regulated system decouples growth from heterologous gene expression without the need for adding an inducer. The minimal pst promoter is identified and contains a conserved AT-rich region, a factor B recognition element, and a TATA box for phosphate-dependent regulation. Rational changes to the factor B recognition element and start codon had no significant impact on expression; however, changes to the transcription start site and the 5' untranslated region resulted in the differential protein production with regulation remaining intact. Compared to a previous expression system based upon the histone promoter, this regulated expression system resulted in significant improvements in the expression of a key methanogenic enzyme complex, methyl-coenzyme M reductase, and the potentially toxic arginine methyltransferase MmpX.

Proteomic perspectives on thermotolerant microbes: an updated review

INTRODUCTION

Thermotolerant microbes are a group of microorganisms that survive in elevated temperatures. The thermotolerant microbes, which are found in geothermal heat zones, grow at temperatures of or above 45°C. The proteins present in such microbes are optimally active at these elevated temperatures. Hence, therefore, serves as an advantage in various biotechnological applications. In the last few years, scientists have tried to understand the molecular mechanisms behind the maintenance of the structural integrity of the cell and to study the stability of various thermotolerant proteins at extreme temperatures. Proteomic analysis is the solution for this search. Applying novel proteomic tools determines the proteins involved in the thermostability of microbes at elevated temperatures.

METHODS

Advanced proteomic techniques like Mass spectrometry, nano-LC-MS, protein microarray, ICAT, iTRAQ, and SILAC could enable the screening and identification of novel thermostable proteins.

RESULTS

This review provides up-to-date details on the protein signature of various thermotolerant microbes analyzed through advanced proteomic tools concerning relevant research articles. The protein complex composition from various thermotolerant microbes cultured at different temperatures, their structural arrangement, and functional efficiency of the protein was reviewed and reported.

CONCLUSION

This review provides an overview of thermotolerant microbes, their enzymes, and the proteomic tools implemented to characterize them. This article also reviewed a comprehensive view of the current proteomic approaches for protein profiling in thermotolerant microbes.

The Exploration of Novel Regulatory Relationships Drives Haloarchaeal Operon-Like Structural Dynamics over Short Evolutionary Distances

Operons are a dominant feature of bacterial and archaeal genome organization. Numerous investigations have related aspects of operon structure to operon function, making operons exemplars for studies aimed at deciphering Nature's design principles for genomic organization at a local scale. We consider this understanding to be both fundamentally important and ultimately useful in the de novo design of increasingly complex synthetic circuits. Here we analyze the evolution of the genomic context of operon-like structures in a set of 76 sequenced and annotated species of halophilic archaea. The phylogenetic depth and breadth of this dataset allows insight into changes in operon-like structures over shorter evolutionary time scales than have been studied in previous cross-species analysis of operon evolution. Our analysis, implemented in the updated software package JContextExplorer finds that operon-like context as measured by changes in structure frequently differs from a sequence divergence model of whole-species phylogeny and that changes seem to be dominated by the exploration of novel regulatory relationships.

Whole-genome comparison between the type strain of Halobacterium salinarum (DSM 3754T ) and the laboratory strains R1 and NRC-1

Halobacterium salinarum is an extremely halophilic archaeon that is widely distributed in hypersaline environments and was originally isolated as a spoilage organism of salted fish and hides. The type strain 91-R6 (DSM 3754T ) has seldom been studied and its genome sequence has only recently been determined by our group. The exact relationship between the type strain and two widely used model strains, NRC-1 and R1, has not been described before. The genome of Hbt. salinarum strain 91-R6 consists of a chromosome (2.17 Mb) and two large plasmids (148 and 102 kb, with 39,230 bp being duplicated). Cytosine residues are methylated (m4 C) within CTAG motifs. The genomes of type and laboratory strains are closely related, their chromosomes sharing average nucleotide identity (ANIb) values of 98% and in silico DNA-DNA hybridization (DDH) values of 95%. The chromosomes are completely colinear, do not show genome rearrangement, and matching segments show <1% sequence difference. Among the strain-specific sequences are three large chromosomal replacement regions (>10 kb). The well-studied AT-rich island (61 kb) of the laboratory strains is replaced by a distinct AT-rich sequence (47 kb) in 91-R6. Another large replacement (91-R6: 78 kb, R1: 44 kb) codes for distinct homologs of proteins involved in motility and N-glycosylation. Most (107 kb) of plasmid pHSAL1 (91-R6) is very closely related to part of plasmid pHS3 (R1) and codes for essential genes (e.g. arginine-tRNA ligase and the pyrimidine biosynthesis enzyme aspartate carbamoyltransferase). Part of pHS3 (42.5 kb total) is closely related to the largest strain-specific sequence (164 kb) in the type strain chromosome. Genome sequencing unraveled the close relationship between the Hbt. salinarum type strain and two well-studied laboratory strains at the DNA and protein levels. Although an independent isolate, the type strain shows a remarkably low evolutionary difference to the laboratory strains.

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.

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.

Solution structure of an archaeal DUF61 family protein SSO0941 encoded by a gene in the operon of box C/D RNA protein complexes

Domain of unknown function 61 (DUF61) family proteins widely exist in archaea and the genes of DUF61 proteins in crenarchaea are in an operon containing two genes of box C/D RNA protein complexes. Here we report the solution NMR structure of DUF61 family member protein SSO0941, from the hyperthermophilic archaeon Sulfolobus solfataricus. SSO0941 has a rigid core structure and flexible N- and C-terminal regions as well as a negatively-charged independent C-terminal helix. The core structure consists of N- and C-terminal subdomains, in which the C-terminal subdomain shows significant structural similarity with several nucleic acid binding proteins. The structure of SSO0941 is the first representative structure of DUF61 family proteins.

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

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

A novel bifunctional transcriptional regulator of riboflavin metabolism in Archaea

Riboflavin (vitamin B2) is the precursor of flavin mononucleotide (FMN) and flavin adenine dinucleotide, which are essential coenzymes in all free-living organisms. Riboflavin biosynthesis in many Bacteria but not in Archaea is controlled by FMN-responsive riboswitches. We identified a novel bifunctional riboflavin kinase/regulator (RbkR), which controls riboflavin biosynthesis and transport genes in major lineages of Crenarchaeota, Euryarchaeota and Thaumarchaeota. RbkR proteins are composed of the riboflavin kinase domain and a DNA-binding winged helix-turn-helix-like domain. Using comparative genomics, we predicted RbkR operator sites and reconstructed RbkR regulons in 94 archaeal genomes. While the identified RbkR operators showed significant variability between archaeal lineages, the conserved core of RbkR regulons includes riboflavin biosynthesis genes, known/predicted vitamin uptake transporters and the rbkR gene. The DNA motifs and CTP-dependent riboflavin kinase activity of two RbkR proteins were experimentally validated in vitro. The DNA binding activity of RbkR was stimulated by CTP and suppressed by FMN, a product of riboflavin kinase. The crystallographic structure of RbkR from Thermoplasma acidophilum was determined in complex with CTP and its DNA operator revealing key residues for operator and ligand recognition. Overall, this study contributes to our understanding of metabolic and regulatory networks for vitamin homeostasis in Archaea.

Genomic and transcriptomic landscape of Escherichia coli BL21(DE3)

Escherichia coli BL21(DE3) has long served as a model organism for scientific research, as well as a workhorse for biotechnology. Here we present the most current genome annotation of E. coli BL21(DE3) based on the transcriptome structure of the strain that was determined for the first time. The genome was annotated using multiple automated pipelines and compared to the current genome annotation of the closely related strain, E. coli K-12. High-resolution tiling array data of E. coli BL21(DE3) from several different stages of cell growth in rich and minimal media were analyzed to characterize the transcriptome structure and to provide supporting evidence for open reading frames. This new integrated analysis of the genomic and transcriptomic structure of E. coli BL21(DE3) has led to the correction of translation initiation sites for 88 coding DNA sequences and provided updated information for most genes. Additionally, 37 putative genes and 66 putative non-coding RNAs were also identified. The panoramic landscape of the genome and transcriptome of E. coli BL21(DE3) revealed here will allow us to better understand the fundamental biology of the strain and also advance biotechnological applications in industry.

Ethanol production by the hyperthermophilic archaeon Pyrococcus furiosus by expression of bacterial bifunctional alcohol dehydrogenases

Ethanol is an important target for the renewable production of liquid transportation fuels. It can be produced biologically from pyruvate, via pyruvate decarboxylase, or from acetyl‐CoA, by alcohol dehydrogenase E (AdhE). Thermophilic bacteria utilize AdhE, which is a bifunctional enzyme that contains both acetaldehyde dehydrogenase and alcohol dehydrogenase activities. Many of these organisms also contain a separate alcohol dehydrogenase (AdhA) that generates ethanol from acetaldehyde, although the role of AdhA in ethanol production is typically not clear. As acetyl‐CoA is a key central metabolite that can be generated from a wide range of substrates, AdhE can serve as a single gene fuel module to produce ethanol through primary metabolic pathways. The focus here is on the hyperthermophilic archaeon Pyrococcus furiosus, which grows by fermenting sugar to acetate, CO₂ and H₂. Previously, by the heterologous expression of adhA from a thermophilic bacterium, P. furiosus was shown to produce ethanol by a novel mechanism from acetate, mediated by AdhA and the native enzyme aldehyde oxidoreductase (AOR). In this study, the AOR gene was deleted from P. furiosus to evaluate ethanol production directly from acetyl‐CoA by heterologous expression of the adhE gene from eight thermophilic bacteria. Only AdhEs from two Thermoanaerobacter strains showed significant activity in cell‐free extracts of recombinant P. furiosus and supported ethanol production in vivo. In the AOR deletion background, the highest amount of ethanol (estimated 61% theoretical yield) was produced when adhE and adhA from Thermoanaerobacter were co‐expressed.

Detecting differential growth of microbial populations with Gaussian process regression

Microbial growth curves are used to study differential effects of media, genetics, and stress on microbial population growth. Consequently, many modeling frameworks exist to capture microbial population growth measurements. However, current models are designed to quantify growth under conditions for which growth has a specific functional form. Extensions to these models are required to quantify the effects of perturbations, which often exhibit nonstandard growth curves. Rather than assume specific functional forms for experimental perturbations, we developed a general and robust model of microbial population growth curves using Gaussian process (GP) regression. GP regression modeling of high-resolution time-series growth data enables accurate quantification of population growth and allows explicit control of effects from other covariates such as genetic background. This framework substantially outperforms commonly used microbial population growth models, particularly when modeling growth data from environmentally stressed populations. We apply the GP growth model and develop statistical tests to quantify the differential effects of environmental perturbations on microbial growth across a large compendium of genotypes in archaea and yeast. This method accurately identifies known transcriptional regulators and implicates novel regulators of growth under standard and stress conditions in the model archaeal organism Halobacterium salinarum For yeast, our method correctly identifies known phenotypes for a diversity of genetic backgrounds under cyclohexamide stress and also detects previously unidentified oxidative stress sensitivity across a subset of strains. Together, these results demonstrate that the GP models are interpretable, recapitulating biological knowledge of growth response while providing new insights into the relevant parameters affecting microbial population growth.

A Novel Transcriptional Regulator Related to Thiamine Phosphate Synthase Controls Thiamine Metabolism Genes in Archaea

Thiamine (vitamin B₁) is a precursor of thiamine pyrophosphate (TPP), an essential coenzyme in the central metabolism of all living organisms. Bacterial thiamine biosynthesis and salvage genes are controlled at the RNA level by TPP-responsive riboswitches. In Archaea, TPP riboswitches are restricted to the Thermoplasmatales order. Mechanisms of transcriptional control of thiamine genes in other archaeal lineages remain unknown. Using the comparative genomics approach, we identified a novel family of transcriptional regulators (named ThiR) controlling thiamine biosynthesis and transport genes in diverse lineages in the Crenarchaeota phylum as well as in the Halobacteria and Thermococci classes of the Euryarchaeota ThiR regulators are composed of an N-terminal DNA-binding domain and a C-terminal ligand-binding domain, which is similar to the archaeal thiamine phosphate synthase ThiN. By using comparative genomics, we predicted ThiR-binding DNA motifs and reconstructed ThiR regulons in 67 genomes representing all above-mentioned lineages. The predicted ThiR-binding motifs are characterized by palindromic symmetry with several distinct lineage-specific consensus sequences. In addition to thiamine biosynthesis genes, the reconstructed ThiR regulons include various transporters for thiamine and its precursors. Bioinformatics predictions were experimentally validated by in vitro DNA-binding assays with the recombinant ThiR protein from the hyperthermophilic archaeon Metallosphaera yellowstonensis MK1. Thiamine phosphate and, to some extent, TPP and hydroxyethylthiazole phosphate were required for the binding of ThiR to its DNA targets, suggesting that ThiR is derepressed by limitation of thiamine phosphates. The thiamine phosphate-binding residues previously identified in ThiN are highly conserved in ThiR regulators, suggesting a conserved mechanism for effector recognition.

IMPORTANCE

Thiamine pyrophosphate is a cofactor for many essential enzymes for glucose and energy metabolism. Thiamine or vitamin B₁ biosynthesis and its transcriptional regulation in Archaea are poorly understood. We applied the comparative genomics approach to identify a novel family of regulators for the transcriptional control of thiamine metabolism genes in Archaea and reconstructed the respective regulons. The predicted ThiR regulons in archaeal genomes control the majority of thiamine biosynthesis genes. The reconstructed regulon content suggests that numerous uptake transporters for thiamine and/or its precursors are encoded in archaeal genomes. The ThiR regulon was experimentally validated by DNA-binding assays with Metallosphaera spp. These discoveries contribute to our understanding of metabolic and regulatory networks involved in vitamin homeostasis in diverse lineages of Archaea.

Identification of the first transcriptional activator of an archaellum operon in a euryarchaeon

The archaellum is the swimming organelle of the third domain, the Archaea. In the euryarchaeon Methanococcus maripaludis, genes involved in archaella formation, including the three archaellins flaB1, flaB2 and flaB3, are mainly located in the fla operon. Previous studies have shown that transcription of fla genes and expression of Fla proteins are regulated under different growth conditions. In this study, we identify MMP1718 as the first transcriptional activator that directly regulates the fla operon in M. maripaludis. Mutants carrying an in-frame deletion in mmp1718 did not express FlaB2 detected by western blotting. Quantitative reverse transcription PCR analysis of purified RNA from the Δmmp1718 mutant showed that transcription of flaB2 was negligible compared to wildtype cells. In addition, no archaella were observed on the cell surface of the Δmmp1718 mutant. FlaB2 expression and archaellation were restored when the Δmmp1718 mutant was complemented with mmp1718 in trans. Electrophoretic motility shift assay and isothermal titration calorimetry results demonstrated the specific binding of purified MMP1718 to DNA fragments upstream of the fla promoter. Four 6 bp consensus sequences were found immediately upstream of the fla promoter and are considered the putative MMP1718-binding sites. Herein, we designate MMP1718 as EarA, the first euryarchaeal archaellum regulator.

Ancillary contributions of heterologous biotin protein ligase and carbonic anhydrase for CO2 incorporation into 3-hydroxypropionate by metabolically engineered Pyrococcus furiosus

Acetyl-Coenzyme A carboxylase (ACC), malonyl-CoA reductase (MCR), and malonic semialdehyde reductase (MRS) convert HCO₃^- and acetyl-CoA into 3-hydroxypropionate (3HP) in the 3-hydroxypropionate/4-hydroxybutyrate carbon fixation cycle resident in the extremely thermoacidophilic archaeon Metallosphaera sedula. These three enzymes, when introduced into the hyperthermophilic archaeon Pyrococcus furiosus, enable production of 3HP from maltose and CO₂ . Sub-optimal function of ACC was hypothesized to be limiting for production of 3HP, so accessory enzymes carbonic anhydrase (CA) and biotin protein ligase (BPL) from M. sedula were produced recombinantly in Escherichia coli to assess their function. P. furiosus lacks a native, functional CA, while the M. sedula CA (Msed_0390) has a specific activity comparable to other microbial versions of this enzyme. M. sedula BPL (Msed_2010) was shown to biotinylate the β-subunit (biotin carboxyl carrier protein) of the ACC in vitro. Since the native BPLs in E. coli and P. furiosus may not adequately biotinylate the M. sedula ACC, the carboxylase was produced in P. furiosus by co-expression with the M. sedula BPL. The baseline production strain, containing only the ACC, MCR, and MSR, grown in a CO₂ -sparged bioreactor reached titers of approximately 40 mg/L 3HP. Strains in which either the CA or BPL accessory enzyme from M. sedula was added to the pathway resulted in improved titers, 120 or 370 mg/L, respectively. The addition of both M. sedula CA and BPL, however, yielded intermediate titers of 3HP (240 mg/L), indicating that the effects of CA and BPL on the engineered 3HP pathway were not additive, possible reasons for which are discussed. While further efforts to improve 3HP production by regulating gene dosage, improving carbon flux and optimizing bioreactor operation are needed, these results illustrate the ancillary benefits of accessory enzymes for incorporating CO₂ into 3HP production in metabolically engineered P. furiosus, and hint at the important role that CA and BPL likely play in the native 3HP/4HB pathway in M. sedula. Biotechnol. Bioeng. 2016;113: 2652-2660. © 2016 Wiley Periodicals, Inc.

Heterologous Production of an Energy-Conserving Carbon Monoxide Dehydrogenase Complex in the Hyperthermophile Pyrococcus furiosus

Carbon monoxide (CO) is an important intermediate in anaerobic carbon fixation pathways in acetogenesis and methanogenesis. In addition, some anaerobes can utilize CO as an energy source. In the hyperthermophilic archaeon Thermococcus onnurineus, which grows optimally at 80°C, CO oxidation and energy conservation is accomplished by a respiratory complex encoded by a 16-gene cluster containing a CO dehydrogenase, a membrane-bound [NiFe]-hydrogenase and a Na⁺/H⁺ antiporter module. This complex oxidizes CO, evolves CO₂ and H₂, and generates a Na⁺ motive force that is used to conserve energy by a Na⁺-dependent ATP synthase. Herein we used a bacterial artificial chromosome to insert the 13.2 kb gene cluster encoding the CO-oxidizing respiratory complex of T. onnurineus into the genome of the heterotrophic archaeon, Pyrococcus furiosus, which grows optimally at 100°C. P. furiosus is normally unable to utilize CO, however, the recombinant strain readily oxidized CO and generated H₂ at 80°C. Moreover, CO also served as an energy source and allowed the P. furiosus strain to grow with a limiting concentration of sugar or with peptides as the carbon source. Moreover, CO oxidation by P. furiosus was also coupled to the re-utilization, presumably for biosynthesis, of acetate generated by fermentation. The functional transfer of CO utilization between Thermococcus and Pyrococcus species demonstrated herein is representative of the horizontal gene transfer of an environmentally relevant metabolic capability. The transfer of CO utilizing, hydrogen-producing genetic modules also has applications for biohydrogen production and a CO-based industrial platform for various thermophilic organisms.

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.

Electronic supplementary material

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

The survival mechanisms of thermophiles at high temperatures: an angle of omics

Thermophiles are referred to as microorganisms with optimal growth temperatures of >60 °C. Over the past few years, a number of studies have been conducted regarding thermophiles, especially using the omics strategies. This review provides a systematic view of the survival physiology of thermophiles from an "omics" perspective, which suggests that the adaptive ability of thermophiles is based on a cooperative mode with multi-dimensional regulations integrating genomics, transcriptomics, and proteomics.

Recombinant expression library of Pyrococcus furiosus constructed by high-throughput cloning: a useful tool for functional and structural genomics

Hyperthermophile Pyrococcus furiosus grows optimally near 100°C and is an important resource of many industrial and molecular biological enzymes. To study the structure and function of P. furiosus proteins at whole genome level, we constructed expression plasmids of each P. furiosus gene using a ligase-independent cloning method, which was based on amplifying target gene and vector by PCR using phosphorothioate-modified primers and digesting PCR products by λ exonuclease. Our cloning method had a positive clone percentage of ≥ 80% in 96-well plate cloning format. Small-scale expression experiment showed that 55 out of 80 genes were efficiently expressed in Escherichia coli Strain Rosetta 2(DE3)pLysS. In summary, this recombinant expression library of P. furiosus provides a platform for functional and structural studies, as well as developing novel industrial enzymes. Our cloning scheme is adaptable to constructing recombinant expression library of other sequenced organisms.

Bioprocessing analysis of Pyrococcus furiosus strains engineered for CO₂-based 3-hydroxypropionate production

Metabolically engineered strains of the hyperthermophile Pyrococcus furiosus (T(opt) 95-100°C), designed to produce 3-hydroxypropionate (3HP) from maltose and CO2 using enzymes from the Metallosphaera sedula (T(opt) 73°C) carbon fixation cycle, were examined with respect to the impact of heterologous gene expression on metabolic activity, fitness at optimal and sub-optimal temperatures, gas-liquid mass transfer in gas-intensive bioreactors, and potential bottlenecks arising from product formation. Transcriptomic comparisons of wild-type P. furiosus, a genetically-tractable, naturally-competent mutant (COM1), and COM1-based strains engineered for 3HP production revealed numerous differences after being shifted from 95°C to 72°C, where product formation catalyzed by the heterologously-produced M. sedula enzymes occurred. At 72°C, significantly higher levels of metabolic activity and a stress response were evident in 3HP-forming strains compared to the non-producing parent strain (COM1). Gas-liquid mass transfer limitations were apparent, given that 3HP titers and volumetric productivity in stirred bioreactors could be increased over 10-fold by increased agitation and higher CO2 sparging rates, from 18 mg/L to 276 mg/L and from 0.7 mg/L/h to 11 mg/L/h, respectively. 3HP formation triggered transcription of genes for protein stabilization and turnover, RNA degradation, and reactive oxygen species detoxification. The results here support the prospects of using thermally diverse sources of pathways and enzymes in metabolically engineered strains designed for product formation at sub-optimal growth temperatures.

A method to analyze, sort, and retain viability of obligate anaerobic microorganisms from complex microbial communities

A high speed flow cytometric cell sorter was modified to maintain a controlled anaerobic environment. This technology enabled coupling of the precise high-throughput analytical and cell separation capabilities of flow cytometry to the assessment of cell viability of evolved lineages of obligate anaerobic organisms from cocultures.

Systems biology approaches to defining transcription regulatory networks in halophilic archaea

To survive complex and changing environmental conditions, microorganisms use gene regulatory networks (GRNs) composed of interacting regulatory transcription factors (TFs) to control the timing and magnitude of gene expression. Genome-wide datasets; such as transcriptomics and protein-DNA interactions; and experiments such as high throughput growth curves; facilitate the construction of GRNs and provide insight into TF interactions occurring under stress. Systems biology approaches integrate these datasets into models of GRN architecture as well as statistical and/or dynamical models to understand the function of networks occurring in cells. Previously, these types of studies have focused on traditional model organisms (e.g. Escherichia coli, yeast). However, recent advances in archaeal genetics and other tools have enabled a systems approach to understanding GRNs in these relatively less studied archaeal model organisms. In this report, we outline a systems biology workflow for generating and integrating data focusing on the TF regulator. We discuss experimental design, outline the process of data collection, and provide the tools required to produce high confidence regulons for the TFs of interest. We provide a case study as an example of this workflow, describing the construction of a GRN centered on multi-TF coordinate control of gene expression governing the oxidative stress response in the hypersaline-adapted archaeon Halobacterium salinarum.

Overlapping genes: a new strategy of thermophilic stress tolerance in prokaryotes

Overlapping genes (OGs) draw the focus of recent day's research. However, the significance of OGs in prokaryotic genomes remained unexplored. As an adaptation to high temperature, thermophiles were shown to eliminate their intergenic regions. Therefore, it could be possible that prokaryotes would increase their OG content to adapt to high temperature. To test this hypothesis, we carried out a comparative study on OG frequency of 256 prokaryotic genomes comprising both thermophiles and non-thermophiles. It was found that thermophiles exhibit higher frequency of overlapping genes than non-thermophiles. Moreover, overlap frequency was found to correlate with optimal growth temperature (OGT) in prokaryotes. Long overlap frequency was found to hold a positive correlation with OGT resulting in an abundance of long overlaps in thermophiles compared to non-thermophiles. On the other hand, short overlap (1-4 nucleotides) frequency (SOF) did not yield any direct correlation with OGT. However, the correlation of SOF with CAIavg (extent of variation of codon usage bias measured as the mean of codon adaptation index of all genes in a given genome) and IG% (proportion of intergenic regions) indicate that they might upregulate the aforementioned factors (CAIavg and IG%) which are already known to be vital forces for thermophilic adaptation. From these evidences, we propose that the OG content bears a strong link to thermophily. Long overlaps are important for their genome compaction and short overlaps are important to uphold high CAIavg. Our findings will surely help in better understanding of the significance of overlapping gene content in prokaryotic genomes.

PAIDB v2.0: exploration and analysis of pathogenicity and resistance islands

Pathogenicity is a complex multifactorial process confounded by the concerted activity of genetic regions associated with virulence and/or resistance determinants. Pathogenicity islands (PAIs) and resistance islands (REIs) are key to the evolution of pathogens and appear to play complimentary roles in the process of bacterial infection. While PAIs promote disease development, REIs give a fitness advantage to the host against multiple antimicrobial agents. The Pathogenicity Island Database (PAIDB, http://www.paidb.re.kr) has been the only database dedicated to providing comprehensive information on all reported PAIs and candidate PAIs in prokaryotic genomes. In this study, we present PAIDB v2.0, whose functionality is extended to incorporate REIs. PAIDB v2.0 contains 223 types of PAIs with 1331 accessions, and 88 types of REIs with 108 accessions. With an improved detection scheme, 2673 prokaryotic genomes were analyzed to locate candidate PAIs and REIs. With additional quantitative and qualitative advancements in database content and detection accuracy, PAIDB will continue to facilitate pathogenomic studies of both pathogenic and non-pathogenic organisms.

Transcription start site associated RNAs (TSSaRNAs) are ubiquitous in all domains of life

A plethora of non-coding RNAs has been discovered using high-resolution transcriptomics tools, indicating that transcriptional and post-transcriptional regulation is much more complex than previously appreciated. Small RNAs associated with transcription start sites of annotated coding regions (TSSaRNAs) are pervasive in both eukaryotes and bacteria. Here, we provide evidence for existence of TSSaRNAs in several archaeal transcriptomes including: Halobacterium salinarum, Pyrococcus furiosus, Methanococcus maripaludis, and Sulfolobus solfataricus. We validated TSSaRNAs from the model archaeon Halobacterium salinarum NRC-1 by deep sequencing two independent small-RNA enriched (RNA-seq) and a primary-transcript enriched (dRNA-seq) strand-specific libraries. We identified 652 transcripts, of which 179 were shown to be primary transcripts (∼7% of the annotated genome). Distinct growth-associated expression patterns between TSSaRNAs and their cognate genes were observed, indicating a possible role in environmental responses that may result from RNA polymerase with varying pausing rhythms. This work shows that TSSaRNAs are ubiquitous across all domains of life.

Primary transcriptome map of the hyperthermophilic archaeon Thermococcus kodakarensis

Background

Prokaryotes have relatively small genomes, densely-packed with protein-encoding sequences. RNA sequencing has, however, revealed surprisingly complex transcriptomes and here we report the transcripts present in the model hyperthermophilic Archaeon, Thermococcus kodakarensis, under different physiological conditions.

Results

Sequencing cDNA libraries, generated from RNA isolated from cells under different growth and metabolic conditions has identified >2,700 sites of transcription initiation, established a genome-wide map of transcripts, and consensus sequences for transcription initiation and post-transcription regulatory elements. The primary transcription start sites (TSS) upstream of 1,254 annotated genes, plus 644 primary TSS and their promoters within genes, are identified. Most mRNAs have a 5'-untranslated region (5'-UTR) 10 to 50 nt long (median = 16 nt), but ~20% have 5'-UTRs from 50 to 300 nt long and ~14% are leaderless. Approximately 50% of mRNAs contain a consensus ribosome binding sequence. The results identify TSS for 1,018 antisense transcripts, most with sequences complementary to either the 5'- or 3'-region of a sense mRNA, and confirm the presence of transcripts from all three CRISPR loci, the RNase P and 7S RNAs, all tRNAs and rRNAs and 69 predicted snoRNAs. Two putative riboswitch RNAs were present in growing but not in stationary phase cells. The procedure used is designed to identify TSS but, assuming that the number of cDNA reads correlates with transcript abundance, the results also provide a semi-quantitative documentation of the differences in T. kodakarensis genome expression under different growth conditions and confirm previous observations of substrate-dependent specific gene expression. Many previously unanticipated small RNAs have been identified, some with relative low GC contents (≤50%) and sequences that do not fold readily into base-paired secondary structures, contrary to the classical expectations for non-coding RNAs in a hyperthermophile.

Conclusion

The results identify >2,700 TSS, including almost all of the primary sites of transcription initiation upstream of annotated genes, plus many secondary sites, sites within genes and sites resulting in antisense transcripts. The T. kodakarensis genome is small (~2.1 Mbp) and tightly packed with protein-encoding genes, but the transcriptomes established also contain many non-coding RNAs and predict extensive RNA-based regulation in this model Archaeon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-684) contains supplementary material, which is available to authorized users.

Conservation of transcription start sites within genes across a bacterial genus

Transcription start sites (TSSs) lying inside annotated genes, on the same or opposite strand, have been observed in diverse bacteria, but the function of these unexpected transcripts is unclear. Here, we use the metal-reducing bacterium Shewanella oneidensis MR-1 and its relatives to study the evolutionary conservation of unexpected TSSs. Using high-resolution tiling microarrays and 5'-end RNA sequencing, we identified 2,531 TSSs in S. oneidensis MR-1, of which 18% were located inside coding sequences (CDSs). Comparative transcriptome analysis with seven additional Shewanella species revealed that the majority (76%) of the TSSs within the upstream regions of annotated genes (gTSSs) were conserved. Thirty percent of the TSSs that were inside genes and on the sense strand (iTSSs) were also conserved. Sequence analysis around these iTSSs showed conserved promoter motifs, suggesting that many iTSS are under purifying selection. Furthermore, conserved iTSSs are enriched for regulatory motifs, suggesting that they are regulated, and they tend to eliminate polar effects, which confirms that they are functional. In contrast, the transcription of antisense TSSs located inside CDSs (aTSSs) was significantly less likely to be conserved (22%). However, aTSSs whose transcription was conserved often have conserved promoter motifs and drive the expression of nearby genes. Overall, our findings demonstrate that some internal TSSs are conserved and drive protein expression despite their unusual locations, but the majority are not conserved and may reflect noisy initiation of transcription rather than a biological function. Importance: The first step of gene expression is the initiation of transcription from promoters, which have been traditionally thought to be located upstream of genes. Recently, studies showed that in diverse bacteria, promoters are often located inside genes. It has not been clear if these unexpected promoters are important to the organism or if they result from transcriptional noise. Here, we identify and examine promoters in eight related bacterial species. Promoters that lie within genes on the sense strand are often conserved as locations and in their sequences. Furthermore, these promoters often affect the bacterium's growth. Thus, many of these unexpected promoters are likely functional. Fewer promoters that lie within genes on the antisense strand are conserved, but the conserved ones seem to drive the expression of nearby genes.

Engineering hydrogen gas production from formate in a hyperthermophile by heterologous production of an 18-subunit membrane-bound complex

Biohydrogen gas has enormous potential as a source of reductant for the microbial production of biofuels, but its low solubility and poor gas mass transfer rates are limiting factors. These limitations could be circumvented by engineering biofuel production in microorganisms that are also capable of generating H2 from highly soluble chemicals such as formate, which can function as an electron donor. Herein, the model hyperthermophile, Pyrococcus furiosus, which grows optimally near 100 °C by fermenting sugars to produce H2, has been engineered to also efficiently convert formate to H2. Using a bacterial artificial chromosome vector, the 16.9-kb 18-gene cluster encoding the membrane-bound, respiratory formate hydrogen lyase complex of Thermococcus onnurineus was inserted into the P. furiosus chromosome and expressed as a functional unit. This enabled P. furiosus to utilize formate as well as sugars as an H2 source and to do so at both 80° and 95 °C, near the optimum growth temperature of the donor (T. onnurineus) and engineered host (P. furiosus), respectively. This accomplishment also demonstrates the versatility of P. furiosus for metabolic engineering applications.

A systems level predictive model for global gene regulation of methanogenesis in a hydrogenotrophic methanogen

Methanogens catalyze the critical methane-producing step (called methanogenesis) in the anaerobic decomposition of organic matter. Here, we present the first predictive model of global gene regulation of methanogenesis in a hydrogenotrophic methanogen, Methanococcus maripaludis. We generated a comprehensive list of genes (protein-coding and noncoding) for M. maripaludis through integrated analysis of the transcriptome structure and a newly constructed Peptide Atlas. The environment and gene-regulatory influence network (EGRIN) model of the strain was constructed from a compendium of transcriptome data that was collected over 58 different steady-state and time-course experiments that were performed in chemostats or batch cultures under a spectrum of environmental perturbations that modulated methanogenesis. Analyses of the EGRIN model have revealed novel components of methanogenesis that included at least three additional protein-coding genes of previously unknown function as well as one noncoding RNA. We discovered that at least five regulatory mechanisms act in a combinatorial scheme to intercoordinate key steps of methanogenesis with different processes such as motility, ATP biosynthesis, and carbon assimilation. Through a combination of genetic and environmental perturbation experiments we have validated the EGRIN-predicted role of two novel transcription factors in the regulation of phosphate-dependent repression of formate dehydrogenase—a key enzyme in the methanogenesis pathway. The EGRIN model demonstrates regulatory affiliations within methanogenesis as well as between methanogenesis and other cellular functions.

Charting microbial phenotypes in multiplex nanoliter batch bioreactors

High-throughput growth phenotyping is receiving great attention for establishing the genotype-phenotype map of sequenced organisms owing to the ready availability of complete genome sequences. To date, microbial growth phenotypes have been investigated mostly by the conventional method of batch cultivation using test tubes, Erlenmeyer flasks, or the recently available microwell plates. However, the current batch cultivation methods are time- and labor-intensive and often fail to consider sophisticated environmental changes. The implementation of batch cultures at the nanoliter scale has been difficult because of the quick evaporation of the culture medium inside the reactors. Here, we report a microfluidic system that allows independent cell cultures in evaporation-free multiplex nanoliter reactors under different culture conditions to assess the behavior of cells. The design allows three experimental replicates for each of eight culture environments in a single run. We demonstrate the versatility of the device by performing growth curve experiments with Escherichia coli and microbiological assays of antibiotics against the opportunistic pathogen Pseudomonas aeruginosa. Our study highlights that the microfluidic system can effectively replace the traditional batch culture methods with nanoliter volumes of bacterial cultivations, and it may be therefore promising for high-throughput growth phenotyping as well as for single-cell analyses.

Quantifying homologous replacement of loci between haloarchaeal species

In vitro studies of the haloarchaeal genus Haloferax have demonstrated their ability to frequently exchange DNA between species, whereas rates of homologous recombination estimated from natural populations in the genus Halorubrum are high enough to maintain random association of alleles between five loci. To quantify the effects of gene transfer and recombination of commonly held (relaxed core) genes during the evolution of the class Halobacteria (haloarchaea), we reconstructed the history of 21 genomes representing all major groups. Using a novel algorithm and a concatenated ribosomal protein phylogeny as a reference, we created a directed horizontal genetic transfer (HGT) network of contemporary and ancestral genomes. Gene order analysis revealed that 90% of testable HGTs were by direct homologous replacement, rather than nonhomologous integration followed by a loss. Network analysis revealed an inverse log-linear relationship between HGT frequency and ribosomal protein evolutionary distance that is maintained across the deepest divergences in Halobacteria. We use this mathematical relationship to estimate the total transfers and amino acid substitutions delivered by HGTs in each genome, providing a measure of chimerism. For the relaxed core genes of each genome, we conservatively estimate that 11–20% of their evolution occurred in other haloarchaea. Our findings are unexpected, because the transfer and homologous recombination of relaxed core genes between members of the class Halobacteria disrupts the coevolution of genes; however, the generation of new combinations of divergent but functionally related genes may lead to adaptive phenotypes not available through cumulative mutations and recombination within a single population.

The genome organization of Thermotoga maritima reflects its lifestyle

The generation of genome-scale data is becoming more routine, yet the subsequent analysis of omics data remains a significant challenge. Here, an approach that integrates multiple omics datasets with bioinformatics tools was developed that produces a detailed annotation of several microbial genomic features. This methodology was used to characterize the genome of Thermotoga maritima—a phylogenetically deep-branching, hyperthermophilic bacterium. Experimental data were generated for whole-genome resequencing, transcription start site (TSS) determination, transcriptome profiling, and proteome profiling. These datasets, analyzed in combination with bioinformatics tools, served as a basis for the improvement of gene annotation, the elucidation of transcription units (TUs), the identification of putative non-coding RNAs (ncRNAs), and the determination of promoters and ribosome binding sites. This revealed many distinctive properties of the T. maritima genome organization relative to other bacteria. This genome has a high number of genes per TU (3.3), a paucity of putative ncRNAs (12), and few TUs with multiple TSSs (3.7%). Quantitative analysis of promoters and ribosome binding sites showed increased sequence conservation relative to other bacteria. The 5′UTRs follow an atypical bimodal length distribution comprised of “Short” 5′UTRs (11–17 nt) and “Common” 5′UTRs (26–32 nt). Transcriptional regulation is limited by a lack of intergenic space for the majority of TUs. Lastly, a high fraction of annotated genes are expressed independent of growth state and a linear correlation of mRNA/protein is observed (Pearson r = 0.63, p<2.2×10⁻¹⁶ t-test). These distinctive properties are hypothesized to be a reflection of this organism's hyperthermophilic lifestyle and could yield novel insights into the evolutionary trajectory of microbial life on earth.

Genomic studies have greatly benefited from the advent of high-throughput technologies and bioinformatics tools. Here, a methodology integrating genome-scale data and bioinformatics tools is developed to characterize the genome organization of the hyperthermophilic, phylogenetically deep-branching bacterium Thermotoga maritima. This approach elucidates several features of the genome organization and enables comparative analysis of these features across diverse taxa. Our results suggest that the genome of T. maritima is reflective of its hyperthermophilic lifestyle. Ultimately, constraints imposed on the genome have negative impacts on regulatory complexity and phenotypic diversity. Investigating the genome organization of Thermotogae species will help resolve various causal factors contributing to the genome organization such as phylogeny and environment. Applying a similar analysis of the genome organization to numerous taxa will likely provide insights into microbial evolution.

Exploiting microbial hyperthermophilicity to produce an industrial chemical, using hydrogen and carbon dioxide

Microorganisms can be engineered to produce useful products, including chemicals and fuels from sugars derived from renewable feedstocks, such as plant biomass. An alternative method is to use low potential reducing power from nonbiomass sources, such as hydrogen gas or electricity, to reduce carbon dioxide directly into products. This approach circumvents the overall low efficiency of photosynthesis and the production of sugar intermediates. Although significant advances have been made in manipulating microorganisms to produce useful products from organic substrates, engineering them to use carbon dioxide and hydrogen gas has not been reported. Herein, we describe a unique temperature-dependent approach that confers on a microorganism (the archaeon Pyrococcus furiosus, which grows optimally on carbohydrates at 100°C) the capacity to use carbon dioxide, a reaction that it does not accomplish naturally. This was achieved by the heterologous expression of five genes of the carbon fixation cycle of the archaeon Metallosphaera sedula, which grows autotrophically at 73°C. The engineered P. furiosus strain is able to use hydrogen gas and incorporate carbon dioxide into 3-hydroxypropionic acid, one of the top 12 industrial chemical building blocks. The reaction can be accomplished by cell-free extracts and by whole cells of the recombinant P. furiosus strain. Moreover, it is carried out some 30°C below the optimal growth temperature of the organism in conditions that support only minimal growth but maintain sufficient metabolic activity to sustain the production of 3-hydroxypropionate. The approach described here can be expanded to produce important organic chemicals, all through biological activation of carbon dioxide.

Genome-scale analysis of gene function in the hydrogenotrophic methanogenic archaeon Methanococcus maripaludis

A comprehensive whole-genome analysis of gene function by transposon mutagenesis and deep sequencing methodology has been implemented successfully in a representative of the Archaea domain. Libraries of transposon mutants were generated for the hydrogenotrophic, methanogenic archaeon Methanococcus maripaludis S2 using a derivative of the Tn5 transposon. About 89,000 unique insertions were mapped to the genome, which allowed for the classification of 526 genes or about 30% of the genome as possibly essential or strongly advantageous for growth in rich medium. Many of these genes were homologous to eukaryotic genes that encode fundamental processes in replication, transcription, and translation, providing direct evidence for their importance in Archaea. Some genes classified as possibly essential were unique to the archaeal or methanococcal lineages, such as that encoding DNA polymerase PolD. In contrast, the archaeal homolog to the gene encoding DNA polymerase B was not essential for growth, a conclusion confirmed by construction of an independent deletion mutation. Thus PolD, and not PolB, likely plays a fundamental role in DNA replication in methanococci. Similarly, 121 hypothetical ORFs were classified as possibly essential and likely play fundamental roles in methanococcal information processing or metabolism that are not established outside this group of prokaryotes.

Expanded target and cofactor repertoire for the transcriptional activator LysM from Sulfolobus

Previously, Lrp-like transcriptional regulator LysM from the hyperthermoacidophilic crenarchaeon Sulfolobus solfataricus was proposed to have a single target, the lysWXJK operon of lysine biosynthesis, and a single effector molecule, l-lysine. Here we identify ∼70 novel binding sites for LysM in the S. solfataricus genome with a LysM-specific nanobody-based chromatin immunoprecipitation assay coupled to microarray hybridization (ChIP-chip) and in silico target site prediction using an energy-based position weight matrix, and validate these findings with in vitro binding. LysM binds to intergenic and coding regions, including promoters of various amino acid biosynthesis and transport genes. We confirm that l-lysine is the most potent effector molecule that reduces, but does not completely abolish, LysM binding, and show that several other amino acids and derivatives, including d-lysine, l-arginine, l-homoarginine, l-glutamine and l-methionine and branched-chain amino acids l-leucine, l-isoleucine and l-valine, significantly affect DNA-binding properties of LysM. Therefore, it appears from this study that LysM is a much more versatile regulator than previously thought, and that it uses a variety of amino acids to sense nutritional quality of the environment and to modulate expression of the metabolic machinery of Sulfolobus accordingly.

Effects of H2 and formate on growth yield and regulation of methanogenesis in Methanococcus maripaludis

Hydrogenotrophic methanogenic Archaea are defined by an H2 requirement for growth. Despite this requirement, many hydrogenotrophs are also capable of growth with formate as an electron donor for methanogenesis. While certain responses of these organisms to hydrogen availability have been characterized, responses to formate starvation have not been reported. Here we report that during continuous culture of Methanococcus maripaludis under defined nutrient conditions, growth yields relative to methane production decreased markedly with either H2 excess or formate excess. Analysis of the growth yields of several mutants suggests that this phenomenon occurs independently of the storage of intracellular carbon or a transcriptional response to methanogenesis. Using microarray analysis, we found that the expression of genes encoding coenzyme F420-dependent steps of methanogenesis, including one of two formate dehydrogenases, increased with H2 starvation but with formate occurred at high levels regardless of limitation or excess. One gene, encoding H2-dependent methylene-tetrahydromethanopterin dehydrogenase, decreased in expression with either H2 limitation or formate limitation. Expression of genes for the second formate dehydrogenase, molybdenum-dependent formylmethanofuran dehydrogenase, and molybdenum transport increased specifically with formate limitation. Of the two formate dehydrogenases, only the first could support growth on formate in batch culture where formate was in excess.

Architecture for interoperable software in biology

Understanding biological complexity demands a combination of high-throughput data and interdisciplinary skills. One way to bring to bear the necessary combination of data types and expertise is by encapsulating domain knowledge in software and composing that software to create a customized data analysis environment. To this end, simple flexible strategies are needed for interconnecting heterogeneous software tools and enabling data exchange between them. Drawing on our own work and that of others, we present several strategies for interoperability and their consequences, in particular, a set of simple data structures—list, matrix, network, table and tuple—that have proven sufficient to achieve a high degree of interoperability. We provide a few guidelines for the development of future software that will function as part of an interoperable community of software tools for biological data analysis and visualization.

Genome sequencing of a genetically tractable Pyrococcus furiosus strain reveals a highly dynamic genome

The model archaeon Pyrococcus furiosus grows optimally near 100°C on carbohydrates and peptides. Its genome sequence (NCBI) was determined 12 years ago. A genetically tractable strain, COM1, was very recently reported, and here we describe its genome sequence. Of 1,909,827 bp in size, it is 1,571 bp longer (0.1%) than the reference NCBI sequence. The COM1 genome contains numerous chromosomal rearrangements, deletions, and single base changes. COM1 also has 45 full or partial insertion sequences (ISs) compared to 35 in the reference NCBI strain, and these have resulted in the direct deletion or insertional inactivation of 13 genes. Another seven genes were affected by chromosomal deletions and are predicted to be nonfunctional. In addition, the amino acid sequences of another 102 of the 2,134 predicted gene products are different in COM1. These changes potentially impact various cellular functions, including carbohydrate, peptide, and nucleotide metabolism; DNA repair; CRISPR-associated defense; transcriptional regulation; membrane transport; and growth at 72°C. For example, the IS-mediated inactivation of riboflavin synthase in COM1 resulted in a riboflavin requirement for growth. Nevertheless, COM1 grew on cellobiose, malto-oligosaccharides, and peptides in complex and minimal media at 98 and 72°C to the same extent as did both its parent strain and a new culture collection strain (DSMZ 3638). This was in spite of COM1 lacking several metabolic enzymes, including nonphosphorylating glyceraldehyde-3-phosphate dehydrogenase and beta-glucosidase. The P. furiosus genome is therefore of high plasticity, and the availability of the COM1 sequence will be critical for the future studies of this model hyperthermophile.

A workflow for genome-wide mapping of archaeal transcription factors with ChIP-seq

Deciphering the structure of gene regulatory networks across the tree of life remains one of the major challenges in postgenomic biology. We present a novel ChIP-seq workflow for the archaea using the model organism Halobacterium salinarum sp. NRC-1 and demonstrate its application for mapping the genome-wide binding sites of natively expressed transcription factors. This end-to-end pipeline is the first protocol for ChIP-seq in archaea, with methods and tools for each stage from gene tagging to data analysis and biological discovery. Genome-wide binding sites for transcription factors with many binding sites (TfbD) are identified with sensitivity, while retaining specificity in the identification the smaller regulons (bacteriorhodopsin-activator protein). Chromosomal tagging of target proteins with a compact epitope facilitates a standardized and cost-effective workflow that is compatible with high-throughput immunoprecipitation of natively expressed transcription factors. The Pique package, an open-source bioinformatics method, is presented for identification of binding events. Relative to ChIP-Chip and qPCR, this workflow offers a robust catalog of protein–DNA binding events with improved spatial resolution and significantly decreased cost. While this study focuses on the application of ChIP-seq in H. salinarum sp. NRC-1, our workflow can also be adapted for use in other archaea and bacteria with basic genetic tools.