Genome-wide analysis of growth phase-dependent translational and transcriptional regulation in halophilic archaea

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

Selective Translation of Low Abundance and Upregulated Transcripts in Halobacterium salinarum

Our findings demonstrate conclusively that low abundance and upregulated transcripts are preferentially translated, potentially by environment-specific translation systems with distinct ribosomal protein composition. We show that a complex interplay of transcriptional and posttranscriptional regulation underlies the conditional and modular regulatory programs that generate ribosomes of distinct protein composition. The modular regulation of ribosomal proteins with other transcription, translation, and metabolic genes is generalizable to bacterial and eukaryotic microbes. These findings are relevant to how microorganisms adapt to unfavorable environments when they transition from active growth to quiescence by generating proteins from upregulated transcripts that are in considerably lower abundance relative to transcripts associated with the previous physiological state. Selective translation of transcripts by distinct ribosomes could form the basis for adaptive evolution to new environments through a modular regulation of the translational systems.

When organisms encounter an unfavorable environment, they transition to a physiologically distinct, quiescent state wherein abundant transcripts from the previous active growth state continue to persist, albeit their active transcription is downregulated. In order to generate proteins for the new quiescent physiological state, we hypothesized that the translation machinery must selectively translate upregulated transcripts in an intracellular milieu crowded with considerably higher abundance transcripts from the previous active growth state. Here, we have analyzed genome-wide changes in the transcriptome (RNA sequencing [RNA-seq]), changes in translational regulation and efficiency by ribosome profiling across all transcripts (ribosome profiling [Ribo-seq]), and protein level changes in assembled ribosomal proteins (sequential window acquisition of all theoretical mass spectra [SWATH-MS]) to investigate the interplay of transcriptional and translational regulation in Halobacterium salinarum as it transitions from active growth to quiescence. We have discovered that interplay of regulatory processes at different levels of information processing generates condition-specific ribosomal complexes to translate preferentially pools of low abundance and upregulated transcripts. Through analysis of the gene regulatory network architecture of H. salinarum, Escherichia coli, and Saccharomyces cerevisiae, we demonstrate that this conditional, modular organization of regulatory programs governing translational systems is a generalized feature across all domains of life.

IMPORTANCE Our findings demonstrate conclusively that low abundance and upregulated transcripts are preferentially translated, potentially by environment-specific translation systems with distinct ribosomal protein composition. We show that a complex interplay of transcriptional and posttranscriptional regulation underlies the conditional and modular regulatory programs that generate ribosomes of distinct protein composition. The modular regulation of ribosomal proteins with other transcription, translation, and metabolic genes is generalizable to bacterial and eukaryotic microbes. These findings are relevant to how microorganisms adapt to unfavorable environments when they transition from active growth to quiescence by generating proteins from upregulated transcripts that are in considerably lower abundance relative to transcripts associated with the previous physiological state. Selective translation of transcripts by distinct ribosomes could form the basis for adaptive evolution to new environments through a modular regulation of the translational systems.

Characterization of the transcriptome of Haloferax volcanii, grown under four different conditions, with mixed RNA-Seq

Haloferax volcanii is a well-established model species for haloarchaea. Small scale RNomics and bioinformatics predictions were used to identify small non-coding RNAs (sRNAs), and deletion mutants revealed that sRNAs have important regulatory functions. A recent dRNA-Seq study was used to characterize the primary transcriptome. Unexpectedly, it was revealed that, under optimal conditions, H. volcanii contains more non-coding sRNAs than protein-encoding mRNAs. However, the dRNA-Seq approach did not contain any length information. Therefore, a mixed RNA-Seq approach was used to determine transcript length and to identify additional transcripts, which are not present under optimal conditions. In total, 50 million paired end reads of 150 nt length were obtained. 1861 protein-coding RNAs (cdRNAs) were detected, which encoded 3092 proteins. This nearly doubled the coverage of cdRNAs, compared to the previous dRNA-Seq study. About 2/3 of the cdRNAs were monocistronic, and 1/3 covered more than one gene. In addition, 1635 non-coding sRNAs were identified. The highest fraction of non-coding RNAs were cis antisense RNAs (asRNAs). Analysis of the length distribution revealed that sRNAs have a median length of about 150 nt. Based on the RNA-Seq and dRNA-Seq results, genes were chosen to exemplify characteristics of the H. volcanii transcriptome by Northern blot analyses, e.g. 1) the transcript patterns of gene clusters can be straightforward, but also very complex, 2) many transcripts differ in expression level under the four analyzed conditions, 3) some genes are transcribed into RNA isoforms of different length, which can be differentially regulated, 4) transcripts with very long 5'-UTRs and with very long 3'-UTRs exist, and 5) about 30% of all cdRNAs have overlapping 3'-ends, which indicates, together with the asRNAs, that H. volcanii makes ample use of sense-antisense interactions. Taken together, this RNA-Seq study, together with a previous dRNA-Seq study, enabled an unprecedented view on the H. volcanii transcriptome.

Multiplexing polysome profiling experiments to study translation in Escherichia coli

Polysome profiling is a widely used method to monitor the translation status of mRNAs. Although it is theoretically a simple technique, it is labor intensive. Repetitive polysome fractionation rapidly generates a large number of samples to be handled in the downstream processes of protein elimination, RNA extraction and quantification. Here, we propose a multiplex polysome profiling experiment in which distinct cellular extracts are pooled before loading on the sucrose gradient for fractionation. We used the multiplexing method to study translation in E. coli. Multiplexing polysome profiling experiments provided similar mRNA translation status to that obtained with the non-multiplex method with comparable distribution of mRNA copies between the polysome profiling fractions, similar ribosome occupancy and ribosome density. The multiplexing method was used for parallel characterization of gene translational responses to changing mRNA levels. When the mRNA level of two native genes, cysZ and lacZ was increased by transcription induction, their global translational response was similar, with a higher ribosome load leading to increased ribosome occupancy and ribosome densities. However the pattern and the magnitude of the translational response were gene specific. By reducing the number of polysome profiling experiments, the multiplexing method saved time and effort and reduced cost and technical bias. This method would be useful to study the translational effect of mRNA sequence-dependent parameters that often require testing multiple samples and conditions in parallel.

Proteomic Study of the Exponential-Stationary Growth Phase Transition in the Haloarchaea Natrialba magadii and Haloferax volcanii

The dynamic changes that take place along the phases of microbial growth (lag, exponential, stationary, and death) have been widely studied in bacteria at the molecular and cellular levels, but little is known for archaea. In this study, a high-throughput approach was used to analyze and compare the proteomes of two haloarchaea during exponential and stationary growth: the neutrophilic Haloferax volcanii and the alkaliphilic Natrialba magadii. Almost 2000 proteins were identified in each species (≈50% of the predicted proteome). Among them, 532 and 432 were found to be differential between growth phases in H. volcanii and N. magadii, respectively. Changes upon entrance into stationary phase included an overall increase in proteins involved in the transport of small molecules and ions, stress response, and fatty acid catabolism. Proteins related to genetic processes and cell division showed a notorious decrease in amount. The data reported in this study not only contributes to our understanding of the exponential-stationary growth phase transition in extremophilic archaea but also provides the first comprehensive analysis of the proteome composition of N. magadii. The MS proteomics data have been deposited in the ProteomeXchange Consortium with the dataset identifier JPST000395.

Differential Gene Expression in Response to Salinity and Temperature in a Haloarcula Strain from Great Salt Lake, Utah

Haloarchaea that inhabit Great Salt Lake (GSL), a thalassohaline terminal lake, must respond to the fluctuating climate conditions of the elevated desert of Utah. We investigated how shifting environmental factors, specifically salinity and temperature, affected gene expression in the GSL haloarchaea, NA6-27, which we isolated from the hypersaline north arm of the lake. Combined data from cultivation, microscopy, lipid analysis, antibiotic sensitivity, and 16S rRNA gene alignment, suggest that NA6-27 is a member of the Haloarcula genus. Our prior study demonstrated that archaea in the Haloarcula genus were stable in the GSL microbial community over seasons and years. In this study, RNA arbitrarily primed PCR (RAP-PCR) was used to determine the transcriptional responses of NA6-27 grown under suboptimal salinity and temperature conditions. We observed alteration of the expression of genes related to general stress responses, such as transcription, translation, replication, signal transduction, and energy metabolism. Of the ten genes that were expressed differentially under stress, eight of these genes responded in both conditions, highlighting this general response. We also noted gene regulation specific to salinity and temperature conditions, such as osmoregulation and transport. Taken together, these data indicate that the GSL Haloarcula strain, NA6-27, demonstrates both general and specific responses to salinity and/or temperature stress, and suggest a mechanistic model for homeostasis that may explain the stable presence of this genus in the community as environmental conditions shift.

Response of extreme haloarchaeon Haloarcula argentinensis RR10 to simulated microgravity in clinorotation

Gravity is the fundamental force that may have operated during the evolution of life on Earth. It is thus important to understand as to what the effects of gravity are on cellular life. The studies related to effect of microgravity on cells may provide greater insights in understanding of how the physical force of gravity shaped life on Earth. The present study focuses on a unique group of organisms called the Haloarchaea, which are known for their extreme resistance to survive in stress-induced environments. The aim of the present investigation was to study the effect of simulated microgravity on physiological response of extremely halophilic archaeon, Haloarcula argentinensis RR10, under slow clinorotation. The growth kinetics of the archaeon in microgravity was studied using the Baryani model and the viable and apoptotic cells were assessed using propidium iodide fluorescent microscopic studies. The physiological mechanism of adaptation was activation of ‘salt-in’ strategy by intracellular sequestration of sodium ions as detected by EDAX. The organism upregulated the production of ribosomal proteins in simulated microgravity as evidenced by Matrix-assisted laser desorption ionization Time of flight–Mass Spectrophotometry. Simulated microgravity altered the antibiotic susceptibility of the haloarchaeon and it developed resistance to Augmentin, Norfloxacin, Tobramycin and Cefoperazone, rendering it a multidrug resistant strain. The presence of antibiotic efflux pump was detected in the haloarchaeon and it also enhanced production of protective carotenoid pigment in simulated microgravity. The present study is presumably the first report of physiological response of H. argentinensis RR10 in microgravity simulated under slow clinorotation.

Effects of Kasugamycin on the Translatome of Escherichia coli

It is long known that Kasugamycin inhibits translation of canonical transcripts containing a 5’-UTR with a Shine Dalgarno (SD) motif, but not that of leaderless transcripts. To gain a global overview of the influence of Kasugamycin on translation efficiencies, the changes of the translatome of Escherichia coli induced by a 10 minutes Kasugamycin treatment were quantified. The effect of Kasugamycin differed widely, 102 transcripts were at least twofold more sensitive to Kasugamycin than average, and 137 transcripts were at least twofold more resistant, and there was a more than 100-fold difference between the most resistant and the most sensitive transcript. The 5’-ends of 19 transcripts were determined from treated and untreated cultures, but Kasugamycin resistance did neither correlate with the presence or absence of a SD motif, nor with differences in 5’-UTR lengths or GC content. RNA Structure Logos were generated for the 102 Kasugamycin-sensitive and for the 137 resistant transcripts. For both groups a short Shine Dalgarno (SD) motif was retrieved, but no specific motifs associated with resistance or sensitivity could be found. Notably, this was also true for the region -3 to -1 upstream of the start codon and the presence of an extended SD motif, which had been proposed to result in Kasugamycin resistance. Comparison of the translatome results with the database RegulonDB showed that the transcript with the highest resistance was leaderless, but no further leaderless transcripts were among the resistant transcripts. Unexpectedly, it was found that translational coupling might be a novel feature that is associated with Kasugamycin resistance. Taken together, Kasugamycin has a profound effect on translational efficiencies of E. coli transcripts, but the mechanism of action is different than previously described.

Genome-wide gene expression and RNA half-life measurements allow predictions of regulation and metabolic behavior in Methanosarcina acetivorans

BACKGROUND

While a few studies on the variations in mRNA expression and half-lives measured under different growth conditions have been used to predict patterns of regulation in bacterial organisms, the extent to which this information can also play a role in defining metabolic phenotypes has yet to be examined systematically. Here we present the first comprehensive study for a model methanogen.

RESULTS

We use expression and half-life data for the methanogen Methanosarcina acetivorans growing on fast- and slow-growth substrates to examine the regulation of its genes. Unlike Escherichia coli where only small shifts in half-lives were observed, we found that most mRNA have significantly longer half-lives for slow growth on acetate compared to fast growth on methanol or trimethylamine. Interestingly, half-life shifts are not uniform across functional classes of enzymes, suggesting the existence of a selective stabilization mechanism for mRNAs. Using the transcriptomics data we determined whether transcription or degradation rate controls the change in transcript abundance. Degradation was found to control abundance for about half of the metabolic genes underscoring its role in regulating metabolism. Genes involved in half of the metabolic reactions were found to be differentially expressed among the substrates suggesting the existence of drastically different metabolic phenotypes that extend beyond just the methanogenesis pathways. By integrating expression data with an updated metabolic model of the organism (iST807) significant differences in pathway flux and production of metabolites were predicted for the three growth substrates.

CONCLUSIONS

This study provides the first global picture of differential expression and half-lives for a class II methanogen, as well as provides the first evidence in a single organism that drastic genome-wide shifts in RNA half-lives can be modulated by growth substrate. We determined which genes in each metabolic pathway control the flux and classified them as regulated by transcription (e.g. transcription factor) or degradation (e.g. post-transcriptional modification). We found that more than half of genes in metabolism were controlled by degradation. Our results suggest that M. acetivorans employs extensive post-transcriptional regulation to optimize key metabolic steps, and more generally that degradation could play a much greater role in optimizing an organism's metabolism than previously thought.

Post-transcriptional RNA regulons affecting cell cycle and proliferation

The cellular growth cycle is initiated and maintained by punctual, yet agile, regulatory events involving modifications of cell cycle proteins as well as coordinated gene expression to support cyclic checkpoint decisions. Recent evidence indicates that post-transcriptional partitioning of messenger RNA subsets by RNA-binding proteins help physically localize, temporally coordinate, and efficiently translate cell cycle proteins. This dynamic organization of mRNAs encoding cell cycle components contributes to the overall economy of the cell cycle consistent with the post-transcriptional RNA regulon model of gene expression. This review examines several recent studies demonstrating the coordination of mRNA subsets encoding cell cycle proteins during nuclear export and subsequent coupling to protein synthesis, and discusses evidence for mRNA coordination of p53 targets and the DNA damage response pathway. We consider how these observations may connect to upstream and downstream post-transcriptional coordination and coupling of splicing, export, localization, and translation. Published examples from yeast, nematode, insect, and mammalian systems are discussed, and we consider genetic evidence supporting the conclusion that dysregulation of RNA regulons may promote pathogenic states of growth such as carcinogenesis.

An evolutionarily conserved RNase-based mechanism for repression of transcriptional positive autoregulation

It is known that environmental context influences the degree of regulation at the transcriptional and post-transcriptional levels. However, the principles governing the differential usage and interplay of regulation at these two levels are not clear. Here, we show that the integration of transcriptional and post-transcriptional regulatory mechanisms in a characteristic network motif drives efficient environment-dependent state transitions. Through phenotypic screening, systems analysis, and rigorous experimental validation, we discovered an RNase (VNG2099C) in Halobacterium salinarum that is transcriptionally co-regulated with genes of the aerobic physiologic state but acts on transcripts of the anaerobic state. Through modelling and experimentation we show that this arrangement generates an efficient state-transition switch, within which RNase-repression of a transcriptional positive autoregulation (RPAR) loop is critical for shutting down ATP-consuming active potassium uptake to conserve energy required for salinity adaptation under aerobic, high potassium, or dark conditions. Subsequently, we discovered that many Escherichia coli operons with energy-associated functions are also putatively controlled by RPAR indicating that this network motif may have evolved independently in phylogenetically distant organisms. Thus, our data suggest that interplay of transcriptional and post-transcriptional regulation in the RPAR motif is a generalized principle for efficient environment-dependent state transitions across prokaryotes.

Small regulatory RNAs in Archaea

Small regulatory RNAs (sRNAs) are universally distributed in all three domains of life, Archaea, Bacteria, and Eukaryotes. In bacteria, sRNAs typically function by binding near the translation start site of their target mRNAs and thereby inhibit or activate translation. In eukaryotes, miRNAs and siRNAs typically bind to the 3′-untranslated region (3′-UTR) of their target mRNAs and influence translation efficiency and/or mRNA stability. In archaea, sRNAs have been identified in all species investigated using bioinformatic approaches, RNomics, and RNA-Seq. Their size can vary significantly between less than 50 to more than 500 nucleotides. Differential expression of sRNA genes has been studied using northern blot analysis, microarrays, and RNA-Seq. In addition, biological functions have been unraveled by genetic approaches, i.e., by characterization of designed mutants. As in bacteria, it was revealed that archaeal sRNAs are involved in many biological processes, including metabolic regulation, adaptation to extreme conditions, stress responses, and even in regulation of morphology and cellular behavior. Recently, the first target mRNAs were identified in archaea, including one sRNA that binds to the 5′-region of two mRNAs in Methanosarcina mazei Gö1 and a few sRNAs that bind to 3′-UTRs in Sulfolobus solfataricus, three Pyrobaculum species, and Haloferax volcanii, indicating that archaeal sRNAs appear to be able to target both the 5′-UTR or the 3′-UTRs of their respective target mRNAs. In addition, archaea contain tRNA-derived fragments (tRFs), and one tRF has been identified as a major ribosome-binding sRNA in H. volcanii, which downregulates translation in response to stress. Besides regulatory sRNAs, archaea contain further classes of sRNAs, e.g., CRISPR RNAs (crRNAs) and snoRNAs.

A comprehensive analysis of the importance of translation initiation factors for Haloferax volcanii applying deletion and conditional depletion mutants

Translation is an important step in gene expression. The initiation of translation is phylogenetically diverse, since currently five different initiation mechanisms are known. For bacteria the three initiation factors IF1 – IF3 are described in contrast to archaea and eukaryotes, which contain a considerably higher number of initiation factor genes. As eukaryotes and archaea use a non-overlapping set of initiation mechanisms, orthologous proteins of both domains do not necessarily fulfill the same function. The genome of Haloferax volcanii contains 14 annotated genes that encode (subunits of) initiation factors. To gain a comprehensive overview of the importance of these genes, it was attempted to construct single gene deletion mutants of all genes. In 9 cases single deletion mutants were successfully constructed, showing that the respective genes are not essential. In contrast, the genes encoding initiation factors aIF1, aIF2γ, aIF5A, aIF5B, and aIF6 were found to be essential. Factors aIF1A and aIF2β are encoded by two orthologous genes in H. volcanii. Attempts to generate double mutants failed in both cases, indicating that also these factors are essential. A translatome analysis of one of the single aIF2β deletion mutants revealed that the translational efficiency of the second ortholog was enhanced tenfold and thus the two proteins can replace one another. The phenotypes of the single deletion mutants also revealed that the two aIF1As and aIF2βs have redundant but not identical functions. Remarkably, the gene encoding aIF2α, a subunit of aIF2 involved in initiator tRNA binding, could be deleted. However, the mutant had a severe growth defect under all tested conditions. Conditional depletion mutants were generated for the five essential genes. The phenotypes of deletion mutants and conditional depletion mutants were compared to that of the wild-type under various conditions, and growth characteristics are discussed.

Integrative "omics"-approach discovers dynamic and regulatory features of bacterial stress responses

Bacteria constantly face stress conditions and therefore mount specific responses to ensure adaptation and survival. Stress responses were believed to be predominantly regulated at the transcriptional level. In the phototrophic bacterium Rhodobacter sphaeroides the response to singlet oxygen is initiated by alternative sigma factors. Further adaptive mechanisms include post-transcriptional and post-translational events, which have to be considered to gain a deeper understanding of how sophisticated regulation networks operate. To address this issue, we integrated three layers of regulation: (1) total mRNA levels at different time-points revealed dynamics of the transcriptome, (2) mRNAs in polysome fractions reported on translational regulation (translatome), and (3) SILAC-based mass spectrometry was used to quantify protein abundances (proteome). The singlet oxygen stress response exhibited highly dynamic features regarding short-term effects and late adaptation, which could in part be assigned to the sigma factors RpoE and RpoH2 generating distinct expression kinetics of corresponding regulons. The occurrence of polar expression patterns of genes within stress-inducible operons pointed to an alternative of dynamic fine-tuning upon stress. In addition to transcriptional activation, we observed significant induction of genes at the post-transcriptional level (translatome), which identified new putative regulators and assigned genes of quorum sensing to the singlet oxygen stress response. Intriguingly, the SILAC approach explored the stress-dependent decline of photosynthetic proteins, but also identified 19 new open reading frames, which were partly validated by RNA-seq. We propose that comparative approaches as presented here will help to create multi-layered expression maps on the system level ("expressome"). Finally, intense mass spectrometry combined with RNA-seq might be the future tool of choice to re-annotate genomes in various organisms and will help to understand how they adapt to alternating conditions.

Contribution of transcriptomics to systems-level understanding of methanogenic Archaea

Methane-producing Archaea are of interest due to their contribution to atmospheric change and for their roles in technological applications including waste treatment and biofuel production. Although restricted to anaerobic environments, methanogens are found in a wide variety of habitats, where they commonly live in syntrophic relationships with bacterial partners. Owing to tight thermodynamic constraints of methanogenesis alone or in syntrophic metabolism, methanogens must carefully regulate their catabolic pathways including the regulation of RNA transcripts. The transcriptome is a dynamic and important control point in microbial systems. This paper assesses the impact of mRNA (transcriptome) studies on the understanding of methanogenesis with special consideration given to how methanogenesis is regulated to cope with nutrient limitation, environmental variability, and interactions with syntrophic partners. In comparison with traditional microarray-based transcriptome analyses, next-generation high-throughput RNA sequencing is greatly advantageous in assessing transcription start sites, the extent of 5′ untranslated regions, operonic structure, and the presence of small RNAs. We are still in the early stages of understanding RNA regulation but it is already clear that determinants beyond transcript abundance are highly relevant to the lifestyles of methanogens, requiring further study.

Regulatory RNAs in archaea: first target identification in Methanoarchaea

sRNAs (small non-coding RNAs) representing important players in many cellular and regulatory processes have been identified in all three domains of life. In Eukarya and Bacteria, functions have been assigned for many sRNAs, whereas the sRNA populations in Archaea are considerably less well characterized. Recent analyses on a genome-wide scale particularly using high-throughput sequencing techniques demonstrated the presence of high numbers of sRNA candidates in several archaea. However, elucidation of the molecular mechanism of sRNA action, as well as understanding their physiological roles, is in general still challenging, particularly in Archaea, since efficient genetic tools are missing. The identification of cellular targets of identified archaeal sRNAs by experimental approaches or computational prediction programs has begun only recently. At present, targets have been identified for one archaeal sRNA, sRNA162 in Methanosarcina mazei, which interacts with the 5′ region of its targets, a cis-encoded and a trans-encoded target, blurring the paradigm of a border between cis- and trans-encoded sRNAs. Besides, the first experimental implications have been obtained in Haloarchaea and Pyrobaculum that archaeal sRNAs also target 3′ regions of mRNAs. The present review summarizes our current knowledge on archaeal sRNAs and their biological functions and targets.

Bacterial translational regulations: high diversity between all mRNAs and major role in gene expression

BACKGROUND

In bacteria, the weak correlations at the genome scale between mRNA and protein levels suggest that not all mRNAs are translated with the same efficiency. To experimentally explore mRNA translational level regulation at the systemic level, the detailed translational status (translatome) of all mRNAs was measured in the model bacterium Lactococcus lactis in exponential phase growth.

RESULTS

Results demonstrated that only part of the entire population of each mRNA species was engaged in translation. For transcripts involved in translation, the polysome size reached a maximum of 18 ribosomes. The fraction of mRNA engaged in translation (ribosome occupancy) and ribosome density were not constant for all genes. This high degree of variability was analyzed by bioinformatics and statistical modeling in order to identify general rules of translational regulation. For most of the genes, the ribosome density was lower than the maximum value revealing major control of translation by initiation. Gene function was a major translational regulatory determinant. Both ribosome occupancy and ribosome density were particularly high for transcriptional regulators, demonstrating the positive role of translational regulation in the coordination of transcriptional networks. mRNA stability was a negative regulatory factor of ribosome occupancy and ribosome density, suggesting antagonistic regulation of translation and mRNA stability. Furthermore, ribosome occupancy was identified as a key component of intracellular protein levels underlining the importance of translational regulation.

CONCLUSIONS

We have determined, for the first time in a bacterium, the detailed translational status for all mRNAs present in the cell. We have demonstrated experimentally the high diversity of translational states allowing individual gene differentiation and the importance of translation-level regulation in the complex process linking gene expression to protein synthesis.

High throughput sequencing reveals a plethora of small RNAs including tRNA derived fragments in Haloferax volcanii

To define the complete sRNA population of the halophilic archaeon Haloferax volcanii, we employed high throughput sequencing. cDNAs were generated from RNA ranging in size from 17 to 500 nucleotides isolated from cells grown at three different conditions to exponential and stationary phase, respectively. Altogether, 145 intergenic and 45 antisense sRNAs were identified. Comparison of the expression profile showed different numbers of reads at the six different conditions for the majority of sRNAs. A striking difference in the number of sRNA reads was observed between cells grown under standard vs. low salt conditions. Furthermore, the six highest numbers of reads were found for low salt conditions. In contrast, only slight differences between sRNA reads at different growth temperatures were detected. Attempts to delete four sRNA genes revealed that one sRNA gene is essential. The three viable sRNA gene deletion mutants possessed distinct phenotypes. According to microarray analyses, the removal of the sRNA gene resulted in a profound change of the transcriptome when compared with the wild type. High throughput sequencing also showed the presence of high concentrations of tRNA derived fragments in H. volcanii. These tRF molecules were shown to have different amounts of reads at the six conditions analyzed. Northern analysis was used to confirm the presence of the tRNA-derived fragments.

Functional genomic and advanced genetic studies reveal novel insights into the metabolism, regulation, and biology of Haloferax volcanii

The genome sequence of Haloferax volcanii is available and several comparative genomic in silico studies were performed that yielded novel insight for example into protein export, RNA modifications, small non-coding RNAs, and ubiquitin-like Small Archaeal Modifier Proteins. The full range of functional genomic methods has been established and results from transcriptomic, proteomic and metabolomic studies are discussed. Notably, Hfx. volcanii is together with Halobacterium salinarum the only prokaryotic species for which a translatome analysis has been performed. The results revealed that the fraction of translationally-regulated genes in haloarchaea is as high as in eukaryotes. A highly efficient genetic system has been established that enables the application of libraries as well as the parallel generation of genomic deletion mutants. Facile mutant generation is complemented by the possibility to culture Hfx. volcanii in microtiter plates, allowing the phenotyping of mutant collections. Genetic approaches are currently used to study diverse biological questions–from replication to posttranslational modification—and selected results are discussed. Taken together, the wealth of functional genomic and genetic tools make Hfx. volcanii a bona fide archaeal model species, which has enabled the generation of important results in recent years and will most likely generate further breakthroughs in the future.

New insights on gene regulation in archaea

Archaea represent an important and vast domain of life. This cellular domain includes a large diversity of organisms characterized as prokaryotes with basal transcriptional machinery similar to eukarya. In this work we explore the most recent findings concerning the transcriptional regulatory organization in archaeal genomes since the perspective of the DNA-binding transcription factors (TFs), such as the high proportion of archaeal TFs homologous to bacteria, the apparent deficit of TFs, only comparable to the proportion of TFs in parasites or intracellular pathogenic bacteria, suggesting a deficit in this class of proteins. We discuss an appealing hypothesis to explain the apparent deficit of TFs in archaea, based on their characteristics, such as their small length sizes. The hypothesis suggests that a large fraction of these small-sized TFs could supply the deficit of TFs in archaea, by forming different combinations of monomers similar to that observed in eukaryotic transcriptional machinery, where a wide diversity of protein-protein interactions could act as mediators of regulatory feedback, indicating a chimera of bacterial and eukaryotic TFs' functionality. Finally, we discuss how global experiments can help to understand in a global context the role of TFs in these organisms.

Growth phase-dependent gene regulation in vivo in Sulfolobus solfataricus

Ribosomal genes are strongly regulated dependent on growth phase in all organisms, but this regulation is poorly understood in Archaea. Moreover, very little is known about growth phase-dependent gene regulation in Archaea. SSV1-based lacS reporter gene constructs containing the Sulfolobus 16S/23S rRNA gene core promoter, the TF55α core promoter, or the native lacS promoter were tested in Sulfolobus solfataricus cells lacking the lacS gene. The 42-bp 16S/23S rRNA gene and 39-bp TF55α core promoters are sufficient for gene expression in S. solfataricus. However, only gene expression driven by the 16S/23S rRNA gene core promoter is dependent on the culture growth phase. This is the smallest known regulated promoter in Sulfolobus. To our knowledge, this is the first study to show growth phase-dependent rRNA gene regulation in Archaea.

Towards a systems approach in the genetic analysis of archaea: Accelerating mutant construction and phenotypic analysis in Haloferax volcanii

With the availability of a genome sequence and increasingly sophisticated genetic tools, Haloferax volcanii is becoming a model for both Archaea and halophiles. In order for H. volcanii to reach a status equivalent to Escherichia coli, Bacillus subtilis, or Saccharomyces cerevisiae, a gene knockout collection needs to be constructed in order to identify the archaeal essential gene set and enable systematic phenotype screens. A streamlined gene-deletion protocol adapted for potential automation was implemented and used to generate 22 H. volcanii deletion strains and identify several potentially essential genes. These gene deletion mutants, generated in this and previous studies, were then analyzed in a high-throughput fashion to measure growth rates in different media and temperature conditions. We conclude that these high-throughput methods are suitable for a rapid investigation of an H. volcanii mutant library and suggest that they should form the basis of a larger genome-wide experiment.

Large scale physiological readjustment during growth enables rapid, comprehensive and inexpensive systems analysis

Background

Rapidly characterizing the operational interrelationships among all genes in a given organism is a critical bottleneck to significantly advancing our understanding of thousands of newly sequenced microbial and eukaryotic species. While evolving technologies for global profiling of transcripts, proteins, and metabolites are making it possible to comprehensively survey cellular physiology in newly sequenced organisms, these experimental techniques have not kept pace with sequencing efforts. Compounding these technological challenges is the fact that individual experiments typically only stimulate relatively small-scale cellular responses, thus requiring numerous expensive experiments to survey the operational relationships among nearly all genetic elements. Therefore, a relatively quick and inexpensive strategy for observing changes in large fractions of the genetic elements is highly desirable.

Results

We have discovered in the model organism Halobacterium salinarum NRC-1 that batch culturing in complex medium stimulates meaningful changes in the expression of approximately two thirds of all genes. While the majority of these changes occur during transition from rapid exponential growth to the stationary phase, several transient physiological states were detected beyond what has been previously observed. In sum, integrated analysis of transcript and metabolite changes has helped uncover growth phase-associated physiologies, operational interrelationships among two thirds of all genes, specialized functions for gene family members, waves of transcription factor activities, and growth phase associated cell morphology control.

Conclusions

Simple laboratory culturing in complex medium can be enormously informative regarding the activities of and interrelationships among a large fraction of all genes in an organism. This also yields important baseline physiological context for designing specific perturbation experiments at different phases of growth. The integration of such growth and perturbation studies with measurements of associated environmental factor changes is a practical and economical route for the elucidation of comprehensive systems-level models of biological systems.

A new hybrid protein for production of recombinant bacteriorhodopsin in Escherichia coli

Unique properties of bacteriorhodopsin, namely, photochromism and high thermal stability, make this protein an attractive target for physico-chemical studies, as well as for various biotechnological applications. Using Mistic as a suitable carrier for insertion of recombinant membrane proteins into cytoplasmic membrane of Escherichia coli, we developed a system for overexpression of bacteriorhodopsin and worked out an efficient procedure for its purification and renaturation with the final yield of 120 mg/l of refolded protein, which is the highest value reported to date for bacteriorhodopsin produced in E. coli. Functional activity of recombinant bacteriorhodopsin was confirmed by spectroscopic and electrochemical assays.

Identification and genomic analysis of transcription factors in archaeal genomes exemplifies their functional architecture and evolutionary origin

Archaea, which represent a large fraction of the phylogenetic diversity of organisms, are prokaryotes with eukaryote-like basal transcriptional machinery. This organization makes the study of their DNA-binding transcription factors (TFs) and their transcriptional regulatory networks particularly interesting. In addition, there are limited experimental data regarding their TFs. In this work, 3,918 TFs were identified and exhaustively analyzed in 52 archaeal genomes. TFs represented less than 5% of the gene products in all the studied species comparable with the number of TFs identified in parasites or intracellular pathogenic bacteria, suggesting a deficit in this class of proteins. A total of 75 families were identified, of which HTH_3, AsnC, TrmB, and ArsR families were universally and abundantly identified in all the archaeal genomes. We found that archaeal TFs are significantly small compared with other protein-coding genes in archaea as well as bacterial TFs, suggesting that a large fraction of these small-sized TFs could supply the probable deficit of TFs in archaea, by possibly forming different combinations of monomers similar to that observed in eukaryotic transcriptional machinery. Our results show that although the DNA-binding domains of archaeal TFs are similar to bacteria, there is an underrepresentation of ligand-binding domains in smaller TFs, which suggests that protein–protein interactions may act as mediators of regulatory feedback, indicating a chimera of bacterial and eukaryotic TFs’ functionality. The analysis presented here contributes to the understanding of the details of transcriptional apparatus in archaea and provides a framework for the analysis of regulatory networks in these organisms.

Hot Transcriptomics

DNA microarray technology allows for a quick and easy comparison of complete transcriptomes, resulting in improved molecular insight in fluctuations of gene expression. After emergence of the microarray technology about a decade ago, the technique has now matured and has become routine in many molecular biology laboratories. Numerous studies have been performed that have provided global transcription patterns of many organisms under a wide range of conditions. Initially, implementation of this high-throughput technology has lead to high expectations for ground breaking discoveries. Here an evaluation is performed of the insight that transcriptome analysis has brought about in the field of hyperthermophilic archaea. The examples that will be discussed have been selected on the basis of their impact, in terms of either biological insight or technological progress.

Transcriptome and proteome exploration to model translation efficiency and protein stability in Lactococcus lactis

This genome-scale study analysed the various parameters influencing protein levels in cells. To achieve this goal, the model bacterium Lactococcus lactis was grown at steady state in continuous cultures at different growth rates, and proteomic and transcriptomic data were thoroughly compared. Ratios of mRNA to protein were highly variable among proteins but also, for a given gene, between the different growth conditions. The modeling of cellular processes combined with a data fitting modeling approach allowed both translation efficiencies and degradation rates to be estimated for each protein in each growth condition. Estimated translational efficiencies and degradation rates strongly differed between proteins and were tested for their biological significance through statistical correlations with relevant parameters such as codon or amino acid bias. These efficiencies and degradation rates were not constant in all growth conditions and were inversely proportional to the growth rate, indicating a more efficient translation at low growth rate but an antagonistic higher rate of protein degradation. Estimated protein median half-lives ranged from 23 to 224 min, underlying the importance of protein degradation notably at low growth rates. The regulation of intracellular protein level was analysed through regulatory coefficient calculations, revealing a complex control depending on protein and growth conditions. The modeling approach enabled translational efficiencies and protein degradation rates to be estimated, two biological parameters extremely difficult to determine experimentally and generally lacking in bacteria. This method is generic and can now be extended to other environments and/or other micro-organisms.

This work is in the field of systems biology. Via an in-depth comparison of proteomic and transcriptomic data in various culture conditions, our objective was to better understand the regulation of protein levels. We have demonstrated that bacteria exert a tight control on intracellular protein levels, through a multi-level regulation involving translation but also dilution due to growth and protein degradation. We have estimated translational efficiencies and protein degradation rates by modeling. These two biological parameters are extremely difficult to measure experimentally and have not been previously determined in bacteria. We have found that they are growth rate dependent, indicating a fine control of translation and degradation processes. We have worked with the small genome bacterium Lactococcus lactis on a limited number of mRNA-protein couples but keeping in mind that this approach could be extended to other micro-organisms and biological phenomena. We have exhibited that mathematical modeling associated to experimental steady-states cultures is a powerful tool to understand microbial physiology.

Regulation of translation in haloarchaea: 5'- and 3'-UTRs are essential and have to functionally interact in vivo

Recently a first genome-wide analysis of translational regulation using prokaryotic species had been performed which revealed that regulation of translational efficiency plays an important role in haloarchaea. In fact, the fractions of genes under differential growth phase-dependent translational control in the two species Halobacterium salinarum and Haloferax volcanii were as high as in eukaryotes. However, nothing is known about the mechanisms of translational regulation in archaea. Therefore, two genes exhibiting opposing directions of regulation were selected to unravel the importance of untranslated regions (UTRs) for differential translational control in vivo.

Differential translational regulation in exponentially growing versus stationary phase cells was studied by comparing translational efficiencies using a reporter gene system. Translational regulation was not observed when 5′-UTRs or 3′-UTRs alone were fused to the reporter gene. However, their simultaneous presence was sufficient to transfer differential translational control from the native transcript to the reporter transcript. This was true for both directions of translational control. Translational regulation was completely abolished when stem loops in the 5′-UTR were changed by mutagenesis. An “UTR-swap” experiment demonstrated that the direction of translational regulation is encoded in the 3′-UTR, not in the 5′-UTR. While much is known about 5′-UTR-dependent translational control in bacteria, the reported findings provide the first examples that both 5′- and 3′-UTRs are essential and sufficient to drive differential translational regulation in a prokaryote and therefore have to functionally interact in vivo. The current results indicate that 3′-UTR-dependent translational control had already evolved before capping and polyadenylation of transcripts were invented, which are essential for circularization of transcripts in eukaryotes.

Characterization of a Haloferax volcanii member of the enolase superfamily: deletion mutant construction, expression analysis, and transcriptome comparison

The enolase superfamily (COG4948) contains proteins with very different biological functions including regulators like the Escherichia coli RspA and metabolic enzymes like enolase. To unravel the biological function of an archaeal family member, an in frame deletion mutant of a gene encoding a COG4948 protein of Haloferax volcanii was generated. The mutant had a lag phase of 3 days after transition from a richer to a poorer medium, in contrast to the wild-type, and the gene was therefore named "important for transition" (iftA). After inoculation of fresh casamino acids or complex medium with stationary phase wild-type cells, the transcript level of iftA was transiently induced at the onset of growth. In contrast, in minimal (or "poor") glucose medium, both transcript and protein were present throughout growth, even in late stationary phase. A comparison of the transcriptomes of deletion mutant and wild-type revealed that transcript levels of a very restricted set of genes were differentially regulated, including genes encoding proteins involved in phosphate metabolism, regulators and stress response proteins. Taken together, the results indicate that IftA might have a dual function, i.e., transiently after transition to fresh medium and permanently during growth in glucose medium.

Genomics and functional genomics with haloarchaea

The first haloarchaeal genome was published in 2000 and today five genome sequences are available. Transcriptome and proteome analyses have been established for two and three haloarchaeal species, respectively, and more than 20 studies using these functional genomic approaches have been published in the last two years. These studies gave global overviews of metabolic regulation (aerobic and anaerobic respiration, phototrophy, carbon source usage), stress response (UV, X-rays, transition metals, osmotic and temperature stress), cell cycle-dependent transcript level regulation, and transcript half-lives. The only translatome analysis available for any prokaryotic species revealed that 10 and 20% of all transcripts are translationally regulated in Haloferax volcanii and Halobacterium salinarum, respectively. Very effective methods for the construction of in frame deletion mutants have been established recently for haloarchaea and are intensively used to unravel the biological roles of genes in this group. Bioinformatic analyses include both cross-genome comparisons as well as integration of genomic data with experimental results. The first systems biology approaches have been performed that used experimental data to construct predictive models of gene expression and metabolism, respectively. In this contribution the current status of genomics, functional genomics, and molecular genetics of haloarchaea is summarized and selected examples are discussed.

Evolution in the laboratory: the genome of Halobacterium salinarum strain R1 compared to that of strain NRC-1

We report the sequence of the Halobacterium salinarum strain R1 chromosome and its four megaplasmids. Our set of protein-coding genes is supported by extensive proteomic and sequence homology data. The structures of the plasmids, which show three large-scale duplications (adding up to 100 kb), were unequivocally confirmed by cosmid analysis. The chromosome of strain R1 is completely colinear and virtually identical to that of strain NRC-1. Correlation of the plasmid sequences revealed 210 kb of sequence that occurs only in strain R1. The remaining 350 kb shows virtual sequence identity in the two strains. Nevertheless, the number and overall structure of the plasmids are largely incompatible. Also, 20% of the protein sequences differ despite the near identity at the DNA sequence level. Finally, we report genome-wide mobility data for insertion sequences from which we conclude that strains R1 and NRC-1 originate from the same natural isolate. This exemplifies evolution in the laboratory.