Evolution of context dependent regulation by expansion of feast/famine regulatory proteins

Lrp as a potential transcriptional regulator involved in stress response in Haloferax mediterranei

The Archaea domain consists of a heterogeneous group of microorganisms with unique physiological properties that occupy a wide variety of niches in nature. Haloferax mediterranei is an extremely halophilic archaeon classified in the Phylum Euryarchaeota, which requires a high concentration of inorganic salts for optimal growth. In haloarchaea, transcription factors play a fundamental role in an adequate adaptation to environmental and nutritional changes, preserving the survival and integrity of the organism. To deepen knowledge of the Lrp/AsnC transcriptional regulator family, a lrp gene (HFX_RS01210) from this family has been studied. Site-directed mutagenesis has allowed us to identify the TATA-box and two potential sites of the transcriptional factor (TF) to its own promoter and autoregulate itself. Several approaches were carried out to elucidate whether this transcriptional regulator is involved in stresses due to heavy metals and limited nitrogen conditions. Characterization of the lrp deletion mutant and the Lrp overexpressed strain, suggests that the level of lrp expression depends on the nitrogen source and the presence of cobalt. The most striking results were obtained in the presence of nitrate as a nitrogen source due to the inability of the deletion mutant to grow. All these results confirm that Lrp is a powerful candidate for a regulatory role in the stress response, particularly under N-limiting conditions and the presence of cobalt.

Genomic analysis of heavy metal-resistant Halobacterium salinarum isolated from Sfax solar saltern sediments

The draft genome sequences of five archaeal strains, isolated from Sfax solar saltern sediments and affiliated with Halobacterium salinarum, were analyzed in order to reveal their adaptive strategies to live in hypersaline environments polluted with heavy metals. The genomes of the strains (named AS1, AS2, AS8, AS11, and AS19) are found to contain 2,060,688; 2,467,461; 2,236,624; 2,432,692; and 2,428,727 bp respectively, with a G + C content of 65.5, 66.0, 67.0, and 66.2%. The majority of these genes (43.69–55.65%) are annotated as hypothetical proteins. Growth under osmotic stress is possible by genes coding for potassium uptake, sodium efflux, and kinases, as well as stress proteins, DNA repair systems, and proteasomal components. These strains harbor many genes responsible for metal transport/resistance, such as: copper-translocating P-type ATPases, ABC transporter, and cobalt-zinc-cadmium resistance protein. In addition, detoxification enzymes and secondary metabolites are also identified. The results show strain AS1, as compared to the other strains, is more adapted to heavy metals and may be used in the bioremediation of multi-metal contaminated environments. This study highlights the presence of several commercially valuable bioproducts (carotenoids, retinal proteins, exopolysaccharide, stress proteins, squalene, and siderophores) and enzymes (protease, sulfatase, phosphatase, phosphoesterase, and chitinase) that can be used in many industrial applications.

Analysis of Haloferax mediterranei Lrp Transcriptional Regulator

Haloferax mediterranei is an extremely halophilic archaeon, able to live in hypersaline environments with versatile nutritional requirements, whose study represents an excellent basis in the field of biotechnology. The transcriptional machinery in Archaea combines the eukaryotic basal apparatus and the bacterial regulation mechanisms. However, little is known about molecular mechanisms of gene expression regulation compared with Bacteria, particularly in Haloarchaea. The genome of Hfx. mediterranei contains a gene, lrp (HFX_RS01210), which encodes a transcriptional factor belonging to Lrp/AsnC family. It is located downstream of the glutamine synthetase gene (HFX_RS01205), an enzyme involved in ammonium assimilation and amino acid metabolism. To study this transcriptional factor more deeply, the lrp gene has been homologously overexpressed and purified under native conditions by two chromatographic steps, namely nickel affinity and gel filtration chromatography, showing that Lrp behaves asa tetrameric protein of approximately 67 kDa. Its promoter region has been characterized under different growth conditions using bgaH as a reporter gene. The amount of Lrp protein was also analyzed by Western blotting in different nitrogen sources and under various stress conditions. To sum up, regarding its involvement in the nitrogen cycle, it has been shown that its expression profile does not change in response to the nitrogen sources tested. Differences in its expression pattern have been observed under different stress conditions, such as in the presence of hydrogen peroxide or heavy metals. According to these results, the Lrp seems to be involved in a general response against stress factors, acting as a first-line transcriptional regulator.

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

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

The Role of Stress Proteins in Haloarchaea and Their Adaptive Response to Environmental Shifts

Over the years, in order to survive in their natural environment, microbial communities have acquired adaptations to nonoptimal growth conditions. These shifts are usually related to stress conditions such as low/high solar radiation, extreme temperatures, oxidative stress, pH variations, changes in salinity, or a high concentration of heavy metals. In addition, climate change is resulting in these stress conditions becoming more significant due to the frequency and intensity of extreme weather events. The most relevant damaging effect of these stressors is protein denaturation. To cope with this effect, organisms have developed different mechanisms, wherein the stress genes play an important role in deciding which of them survive. Each organism has different responses that involve the activation of many genes and molecules as well as downregulation of other genes and pathways. Focused on salinity stress, the archaeal domain encompasses the most significant extremophiles living in high-salinity environments. To have the capacity to withstand this high salinity without losing protein structure and function, the microorganisms have distinct adaptations. The haloarchaeal stress response protects cells against abiotic stressors through the synthesis of stress proteins. This includes other heat shock stress proteins (Hsp), thermoprotectants, survival proteins, universal stress proteins, and multicellular structures. Gene and family stress proteins are highly conserved among members of the halophilic archaea and their study should continue in order to develop means to improve for biotechnological purposes. In this review, all the mechanisms to cope with stress response by haloarchaea are discussed from a global perspective, specifically focusing on the role played by universal stress proteins.

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.

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.

Conserved principles of transcriptional networks controlling metabolic flexibility in archaea

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

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

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

Systematic Discovery of Archaeal Transcription Factor Functions in Regulatory Networks through Quantitative Phenotyping Analysis

To ensure survival in the face of stress, microorganisms employ inducible damage repair pathways regulated by extensive and complex gene networks. Many archaea, microorganisms of the third domain of life, persist under extremes of temperature, salinity, and pH and under other conditions. In order to understand the cause-effect relationships between the dynamic function of the stress network and ultimate physiological consequences, this study characterized the physiological role of nearly one-third of all regulatory proteins known as transcription factors (TFs) in an archaeal organism. Using a unique quantitative phenotyping approach, we discovered functions for many novel TFs and revealed important secondary functions for known TFs. Surprisingly, many TFs are required for resisting multiple stressors, suggesting cross-regulation of stress responses. Through extensive validation experiments, we map the physiological roles of these novel TFs in stress response back to their position in the regulatory network wiring. This study advances understanding of the mechanisms underlying how microorganisms resist extreme stress. Given the generality of the methods employed, we expect that this study will enable future studies on how regulatory networks adjust cellular physiology in a diversity of organisms.

Gene regulatory networks (GRNs) are critical for dynamic transcriptional responses to environmental stress. However, the mechanisms by which GRN regulation adjusts physiology to enable stress survival remain unclear. Here we investigate the functions of transcription factors (TFs) within the global GRN of the stress-tolerant archaeal microorganism Halobacterium salinarum. We measured growth phenotypes of a panel of TF deletion mutants in high temporal resolution under heat shock, oxidative stress, and low-salinity conditions. To quantitate the noncanonical functional forms of the growth trajectories observed for these mutants, we developed a novel modeling framework based on Gaussian process regression and functional analysis of variance (FANOVA). We employ unique statistical tests to determine the significance of differential growth relative to the growth of the control strain. This analysis recapitulated known TF functions, revealed novel functions, and identified surprising secondary functions for characterized TFs. Strikingly, we observed that the majority of the TFs studied were required for growth under multiple stress conditions, pinpointing regulatory connections between the conditions tested. Correlations between quantitative phenotype trajectories of mutants are predictive of TF-TF connections within the GRN. These phenotypes are strongly concordant with predictions from statistical GRN models inferred from gene expression data alone. With genome-wide and targeted data sets, we provide detailed functional validation of novel TFs required for extreme oxidative stress and heat shock survival. Together, results presented in this study suggest that many TFs function under multiple conditions, thereby revealing high interconnectivity within the GRN and identifying the specific TFs required for communication between networks responding to disparate stressors.

IMPORTANCE To ensure survival in the face of stress, microorganisms employ inducible damage repair pathways regulated by extensive and complex gene networks. Many archaea, microorganisms of the third domain of life, persist under extremes of temperature, salinity, and pH and under other conditions. In order to understand the cause-effect relationships between the dynamic function of the stress network and ultimate physiological consequences, this study characterized the physiological role of nearly one-third of all regulatory proteins known as transcription factors (TFs) in an archaeal organism. Using a unique quantitative phenotyping approach, we discovered functions for many novel TFs and revealed important secondary functions for known TFs. Surprisingly, many TFs are required for resisting multiple stressors, suggesting cross-regulation of stress responses. Through extensive validation experiments, we map the physiological roles of these novel TFs in stress response back to their position in the regulatory network wiring. This study advances understanding of the mechanisms underlying how microorganisms resist extreme stress. Given the generality of the methods employed, we expect that this study will enable future studies on how regulatory networks adjust cellular physiology in a diversity of organisms.

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.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Evolution of a global regulator: Lrp in four orders of γ-Proteobacteria

Background

Bacterial global regulators each regulate the expression of several hundred genes. In Escherichia coli, the top seven global regulators together control over half of all genes. Leucine-responsive regulatory protein (Lrp) is one of these top seven global regulators. Lrp orthologs are very widely distributed, among both Bacteria and Archaea. Surprisingly, even within the phylum γ-Proteobacteria (which includes E. coli), Lrp is a global regulator in some orders and a local regulator in others. This raises questions about the evolution of Lrp and, more broadly, of global regulators.

Results

We examined Lrp sequences from four bacterial orders of the γ-Proteobacteria using phylogenetic and Logo analyses. The orders studied were Enterobacteriales and Vibrionales, in which Lrp plays a global role in tested species; Pasteurellales, in which Lrp is a local regulator in the tested species; and Alteromonadales, an order closely related to the other three but in which Lrp has not yet been studied. For comparison, we analyzed the Lrp paralog AsnC, which in all tested cases is a local regulator. The Lrp and AsnC phylogenetic clusters each divided, as expected, into subclusters representing the Enterobacteriales, Vibrionales, and Pasteuralles. However the Alteromonadales did not yield coherent clusters for either Lrp or AsnC. Logo analysis revealed signatures associated with globally- vs. locally- acting Lrp orthologs, providing testable hypotheses for which portions of Lrp are responsible for a global vs. local role. These candidate regions include both ends of the Lrp polypeptide but not, interestingly, the highly-conserved helix-turn-helix motif responsible for DNA sequence specificity.

Conclusions

Lrp and AsnC have conserved sequence signatures that allow their unambiguous annotation, at least in γ-Proteobacteria. Among Lrp orthologs, specific residues correlated with global vs. local regulatory roles, and can now be tested to determine which are functionally relevant and which simply reflect divergence. In the Alteromonadales, it appears that there are different subgroups of Lrp orthologs, one of which may act globally while the other may act locally. These results suggest experiments to improve our understanding of the evolution of bacterial global regulators.

Electronic supplementary material

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

Inference of expanded Lrp-like feast/famine transcription factor targets in a non-model organism using protein structure-based prediction

Widespread microbial genome sequencing presents an opportunity to understand the gene regulatory networks of non-model organisms. This requires knowledge of the binding sites for transcription factors whose DNA-binding properties are unknown or difficult to infer. We adapted a protein structure-based method to predict the specificities and putative regulons of homologous transcription factors across diverse species. As a proof-of-concept we predicted the specificities and transcriptional target genes of divergent archaeal feast/famine regulatory proteins, several of which are encoded in the genome of Halobacterium salinarum. This was validated by comparison to experimentally determined specificities for transcription factors in distantly related extremophiles, chromatin immunoprecipitation experiments, and cis-regulatory sequence conservation across eighteen related species of halobacteria. Through this analysis we were able to infer that Halobacterium salinarum employs a divergent local trans-regulatory strategy to regulate genes (carA and carB) involved in arginine and pyrimidine metabolism, whereas Escherichia coli employs an operon. The prediction of gene regulatory binding sites using structure-based methods is useful for the inference of gene regulatory relationships in new species that are otherwise difficult to infer.