Reconstruction of the core and extended regulons of global transcription factors

Surface Plasmon Resonance as a Tool to Elucidate the Molecular Determinants of Key Transcriptional Regulators Controlling Rhizobial Lifestyles

Bacteria must be provided with a battery of tools integrated into regulatory networks, in order to respond and, consequently, adapt their physiology to changing environments. Within these networks, transcription factors finely orchestrate the expression of genes in response to a variety of signals, by recognizing specific DNA sequences at their promoter regions. Rhizobia are host-interacting soil bacteria that face severe changes to adapt their physiology from free-living conditions to the nitrogen-fixing endosymbiotic state inside root nodules associated with leguminous plants. One of these cues is the low partial pressure of oxygen within root nodules.Surface plasmon resonance (SPR) constitutes a technique that allows to measure molecular interactions dynamics at real time by detecting changes in the refractive index of a surface. Here, we implemented the SPR methodology to analyze the discriminatory determinants of transcription factors for specific interaction with their target genes. We focused on FixK2, a CRP/FNR-type protein with a central role in the complex oxygen-responsive regulatory network in the soybean endosymbiont Bradyrhizobium diazoefficiens. Our study unveiled relevant residues for protein-DNA interaction as well as allowed us to monitor kinetics and stability protein-DNA complex. We believe that this approach can be employed for the characterization of other relevant transcription factors which can assist to the better understanding of the adaptation of bacteria with agronomic or human interest to their different modes of life.

Data on the solution and processing time reached when constructing a phylogenetic tree using a quantum-inspired computer

Phylogenetic trees provide insight into the evolutionary trajectories of species and molecules. However, because (2n-5)! Phylogenetic trees can be constructed from a dataset containing n sequences, but this method of phylogenetic tree construction is not ideal from the viewpoint of a combinatorial explosion to determine the optimal tree using brute force. Therefore, we developed a method for constructing a phylogenetic tree using a Fujitsu Digital Annealer, a quantum-inspired computer that solves combinatorial optimization problems at a high speed. Specifically, phylogenetic trees are generated by repeating the process of partitioning a set of sequences into two parts (i.e., the graph-cut problem). Here, the optimality of the solution (normalized cut value) obtained by the proposed method was compared with the existing methods using simulated and real data. The simulation dataset contained 32-3200 sequences, and the average branch length according to a normal distribution or the Yule model ranged from 0.125 to 0.750, covering a wide range of sequence diversity. In addition, the statistical information of the dataset is described in terms of two indices: transitivity and average p-distance. As phylogenetic tree construction methods are expected to continue to improve, we believe that this dataset can be used as a reference for comparison and confirmation of the validity of the results. Further interpretation of these analyses is explained in W. Onodera, N. Hara, S. Aoki, T. Asahi, N. Sawamura, Phylogenetic tree reconstruction via graph cut presented using a quantum-inspired computer, Mol. Phylogenet. Evol. 178 (2023) 107636.

Phylogenetic tree reconstruction via graph cut presented using a quantum-inspired computer

Phylogenetic trees are essential tools in evolutionary biology that present information on evolutionary events among organisms and molecules. From a dataset of n sequences, a phylogenetic tree of (2n-5)!! possible topologies exists, and determining the optimum topology using brute force is infeasible. Recently, a recursive graph cut on a graph-represented-similarity matrix has proven accurate in reconstructing a phylogenetic tree containing distantly related sequences. However, identifying the optimum graph cut is challenging, and approximate solutions are currently utilized. Here, a phylogenetic tree was reconstructed with an improved graph cut using a quantum-inspired computer, the Fujitsu Digital Annealer (DA), and the algorithm was named the "Normalized-Minimum cut by Digital Annealer (NMcutDA) method". First, a criterion for the graph cut, the normalized cut value, was compared with existing clustering methods. Based on the cut, we verified that the simulated phylogenetic tree could be reconstructed with the highest accuracy when sequences were diverged. Moreover, for some actual data from the structure-based protein classification database, only NMcutDA could cluster sequences into correct superfamilies. Conclusively, NMcutDA reconstructed better phylogenetic trees than those using other methods by optimizing the graph cut. We anticipate that when the diversity of sequences is sufficiently high, NMcutDA can be utilized with high efficiency.

Fine-Tuning Modulation of Oxidation-Mediated Posttranslational Control of Bradyrhizobium diazoefficiens FixK2 Transcription Factor

FixK2 is a CRP/FNR-type transcription factor that plays a central role in a sophisticated regulatory network for the anoxic, microoxic and symbiotic lifestyles of the soybean endosymbiont Bradyrhizobium diazoefficiens. Aside from the balanced expression of the fixK2 gene under microoxic conditions (induced by the two-component regulatory system FixLJ and negatively auto-repressed), FixK2 activity is posttranslationally controlled by proteolysis, and by the oxidation of a singular cysteine residue (C183) near its DNA-binding domain. To simulate the permanent oxidation of FixK2, we replaced C183 for aspartic acid. Purified C183D FixK2 protein showed both low DNA binding and in vitro transcriptional activation from the promoter of the fixNOQP operon, required for respiration under symbiosis. However, in a B. diazoefficiens strain coding for C183D FixK2, expression of a fixNOQP'-'lacZ fusion was similar to that in the wild type, when both strains were grown microoxically. The C183D FixK2 encoding strain also showed a wild-type phenotype in symbiosis with soybeans, and increased fixK2 gene expression levels and FixK2 protein abundance in cells. These two latter observations, together with the global transcriptional profile of the microoxically cultured C183D FixK2 encoding strain, suggest the existence of a finely tuned regulatory strategy to counterbalance the oxidation-mediated inactivation of FixK2 in vivo.

Exploring the Meta-regulon of the CRP/FNR Family of Global Transcriptional Regulators in a Partial-Nitritation Anammox Microbiome

Microorganisms must respond to environmental changes to survive, often by controlling transcription initiation. Intermittent aeration during wastewater treatment presents a cyclically changing environment to which microorganisms must react. We used an intermittently aerated bioreactor performing partial nitritation and anammox (PNA) to investigate how the microbiome responds to recurring change. Meta-transcriptomic analysis revealed a dramatic disconnect between the relative DNA abundance and gene expression within the metagenome-assembled genomes (MAGs) of community members, suggesting the importance of transcriptional regulation in this microbiome. To explore how community members responded to cyclic aeration via transcriptional regulation, we searched for homologs of the catabolite repressor protein/fumarate and nitrate reductase regulatory protein (CRP/FNR) family of transcription factors (TFs) within the MAGs. Using phylogenetic analyses, evaluation of sequence conservation in important amino acid residues, and prediction of genes regulated by TFs in the MAGs, we identified homologs of the oxygen-sensing FNR in Nitrosomonas and Rhodocyclaceae, nitrogen-sensing dissimilative nitrate respiration regulator that responds to nitrogen species (DNR) in Rhodocyclaceae, and nitrogen-sensing nitrite and nitric oxide reductase regulator that responds to nitrogen species (NnrR) in Nitrospira MAGs. Our data also predict that CRP/FNR homologs in Ignavibacteria, Flavobacteriales, and Saprospiraceae MAGs sense carbon availability. In addition, a CRP/FNR homolog in a Brocadia MAG was most closely related to CRP TFs known to sense carbon sources in well-studied organisms. However, we predict that in autotrophic Brocadia, this TF most likely regulates a diverse set of functions, including a response to stress during the cyclic aerobic/anoxic conditions. Overall, this analysis allowed us to define a meta-regulon of the PNA microbiome that explains functions and interactions of the most active community members.

IMPORTANCE Microbiomes are important contributors to many ecosystems, including ones where nutrient cycling is stimulated by aeration control. Optimizing cyclic aeration helps reduce energy needs and maximize microbiome performance during wastewater treatment; however, little is known about how most microbial community members respond to these alternating conditions. We defined the meta-regulon of a PNA microbiome by combining existing knowledge of how the CRP/FNR family of bacterial TFs respond to stimuli, with metatranscriptomic analyses to characterize gene expression changes during aeration cycles. Our results indicated that, for some members of the community, prior knowledge is sufficient for high-confidence assignments of TF function, whereas other community members have CRP/FNR TFs for which inferences of function are limited by lack of prior knowledge. This study provides a framework to begin elucidating meta-regulons in microbiomes, where pure cultures are not available for traditional transcriptional regulation studies. Defining the meta-regulon can help in optimizing microbiome performance.

Dissection of FixK2 protein-DNA interaction unveils new insights into Bradyrhizobium diazoefficiens lifestyles control

The FixK₂ protein plays a pivotal role in a complex regulatory network, which controls genes for microoxic, denitrifying, and symbiotic nitrogen-fixing lifestyles in Bradyrhizobium diazoefficiens. Among the microoxic-responsive FixK₂ -activated genes are the fixNOQP operon, indispensable for respiration in symbiosis, and the nnrR regulatory gene needed for the nitric-oxide dependent induction of the norCBQD genes encoding the denitrifying nitric oxide reductase. FixK₂ is a CRP/FNR-type transcription factor, which recognizes a 14 bp-palindrome (FixK₂ box) at the regulated promoters through three residues (L195, E196, and R200) within a C-terminal helix-turn-helix motif. Here, we mapped the determinants for discriminatory FixK₂ -mediated regulation. While R200 was essential for DNA binding and activity of FixK₂ , L195 was involved in protein-DNA complex stability. Mutation at positions 1, 3, or 11 in the genuine FixK₂ box at the fixNOQP promoter impaired transcription activation by FixK₂ , which was residual when a second mutation affecting the box palindromy was introduced. The substitution of nucleotide 11 within the NnrR box at the norCBQD promoter allowed FixK₂ -mediated activation in response to microoxia. Thus, position 11 within the FixK₂ /NnrR boxes constitutes a key element that changes FixK₂ targets specificity, and consequently, it might modulate B. diazoefficiens lifestyle as nitrogen fixer or as denitrifier.

Involvement of an FNR-like oxygen sensor in Komagataeibacter medellinensis for survival under oxygen depletion

During acetic acid fermentation, acetic acid bacteria face oxygen depletion stress caused by the vigorous oxidation of ethanol to acetic acid. However, the molecular mechanisms underlying the response to oxygen depletion stress remain largely unknown. Here, we focused on an oxygen-sensing FNR homolog, FnrG, in Komagataeibacter medellinensis. Comparative transcriptomic analysis between the wild-type and fnrG-disrupted strains revealed that FnrG upregulated eight genes (fold change > 3). Recombinant FnrG bound to a specific DNA sequence only when FnrG was reconstituted anaerobically. An operon consisting of acetate kinase and xylulose-5-phosphate/fructose-6-phosphate phosphoketolase genes was found to be an FnrG regulon involved in cell survival under oxygen-limiting conditions. Moreover, a strain that overexpressed these two genes accumulated more acetic acid than the wild-type strain harboring an empty vector. Thus, these two genes could be new targets for the molecular breeding of acetic acid bacteria with high acetic acid productivity.

RegAB Homolog of Burkholderia pseudomallei is the Master Regulator of Redox Control and involved in Virulence

Burkholderia pseudomallei, the etiological agent of melioidosis in humans and animals, often occupies environmental niches and infection sites characterized by limited concentrations of oxygen. Versatile genomic features enable this pathogen to maintain its physiology and virulence under hypoxia, but the crucial regulatory networks employed to switch from oxygen dependent respiration to alternative terminal electron acceptors (TEA) like nitrate, remains poorly understood. Here, we combined a Tn5 transposon mutagenesis screen and an anaerobic growth screen to identify a two-component signal transduction system with homology to RegAB. We show that RegAB is not only essential for anaerobic growth, but also for full virulence in cell lines and a mouse infection model. Further investigations of the RegAB regulon, using a global transcriptomic approach, identified 20 additional regulators under transcriptional control of RegAB, indicating a superordinate role of RegAB in the B. pseudomallei anaerobiosis regulatory network. Of the 20 identified regulators, NarX/L and a FNR homolog were selected for further analyses and a role in adaptation to anaerobic conditions was demonstrated. Growth experiments identified nitrate and intermediates of the denitrification process as the likely signal activateing RegAB, NarX/L, and probably of the downstream regulators Dnr or NsrR homologs. While deletions of individual genes involved in the denitrification process demonstrated their important role in anaerobic fitness, they showed no effect on virulence. This further highlights the central role of RegAB as the master regulator of anaerobic metabolism in B. pseudomallei and that the complete RegAB-mediated response is required to achieve full virulence. In summary, our analysis of the RegAB-dependent modulon and its interconnected regulons revealed a key role for RegAB of B. pseudomallei in the coordination of the response to hypoxic conditions and virulence, in the environment and the host.

Delila-PY, a Pipeline for Utilizing the Delila Suite of Software to Identify Potential DNA Binding Motifs

Predicting potential DNA binding motifs is a critical part of understanding gene expression across all domains of life. Here, we report the development of Delila-PY, an easy-to-use pipeline to utilize the Delila suite to identify DNA binding motifs.

Predicting potential DNA binding motifs is a critical part of understanding gene expression across all domains of life. Here, we report the development of Delila-PY, an easy-to-use pipeline for utilizing the Delila suite of software to identify DNA binding motifs.

Multiple sensors provide spatiotemporal oxygen regulation of gene expression in a Rhizobium-legume symbiosis

Regulation by oxygen (O₂) in rhizobia is essential for their symbioses with plants and involves multiple O₂ sensing proteins. Three sensors exist in the pea microsymbiont Rhizobium leguminosarum Rlv3841: hFixL, FnrN and NifA. At low O₂ concentrations (1%) hFixL signals via FxkR to induce expression of the FixK transcription factor, which activates transcription of downstream genes. These include fixNOQP, encoding the high-affinity cbb₃-type terminal oxidase used in symbiosis. In free-living Rlv3841, the hFixL-FxkR-FixK pathway was active at 1% O₂, and confocal microscopy showed hFixL-FxkR-FixK activity in the earliest stages of Rlv3841 differentiation in nodules (zones I and II). Work on Rlv3841 inside and outside nodules showed that the hFixL-FxkR-FixK pathway also induces transcription of fnrN at 1% O₂ and in the earliest stages of Rlv3841 differentiation in nodules. We confirmed past findings suggesting a role for FnrN in fixNOQP expression. However, unlike hFixL-FxkR-FixK, Rlv3841 FnrN was only active in the near-anaerobic zones III and IV of pea nodules. Quantification of fixNOQP expression in nodules showed this was driven primarily by FnrN, with minimal direct hFixL-FxkR-FixK induction. Thus, FnrN is key for full symbiotic expression of fixNOQP. Without FnrN, nitrogen fixation was reduced by 85% in Rlv3841, while eliminating hFixL only reduced fixation by 25%. The hFixL-FxkR-FixK pathway effectively primes the O₂ response by increasing fnrN expression in early differentiation (zones I-II). In zone III of mature nodules, near-anaerobic conditions activate FnrN, which induces fixNOQP transcription to the level required for wild-type nitrogen fixation activity. Modelling and transcriptional analysis indicates that the different O₂ sensitivities of hFixL and FnrN lead to a nuanced spatiotemporal pattern of gene regulation in different nodule zones in response to changing O₂ concentration. Multi-sensor O₂ regulation is prevalent in rhizobia, suggesting the fine-tuned control this enables is common and maximizes the effectiveness of the symbioses.

Rhizobia are soil bacteria that form a symbiosis with legume plants. In exchange for shelter from the plant, rhizobia provide nitrogen fertilizer, produced by nitrogen fixation. Fixation is catalysed by the nitrogenase enzyme, which is inactivated by oxygen. To prevent this, plants house rhizobia in root nodules, which create a low oxygen environment. However, rhizobia need oxygen, and must adapt to survive the low oxygen concentration in the nodule. Key to this is regulating their genes based on oxygen concentration. We studied one Rhizobium species which uses three different protein sensors of oxygen, each turning on at a different oxygen concentration. As the bacteria get deeper inside the plant nodule and the oxygen concentration drops, each sensor switches on in turn. Our results also show that the first sensor to turn on, hFixL, primes the second sensor, FnrN. This prepares the rhizobia for the core region of the nodule where oxygen concentration is lowest and most nitrogen fixation takes place. If both sensors are removed, the bacteria cannot fix nitrogen. Many rhizobia have several oxygen sensing proteins, so using multiple sensors is likely a common strategy enabling rhizobia to adapt to low oxygen precisely and in stages during symbiosis.

Genome-Wide Identification of Transcription Start Sites in Two Alphaproteobacteria, Rhodobacter sphaeroides 2.4.1 and Novosphingobium aromaticivorans DSM 12444

Here, we report the genome-wide identification of transcription start sites (TSSs) from two Alphaproteobacteria grown under conditions that result in significant changes in gene expression. TSSs that were identified as present in one condition or both will be an important resource for future studies of these, and possibly other, Alphaproteobacteria.

Here, we report the genome-wide identification of transcription start sites (TSSs) from two Alphaproteobacteria grown under conditions that result in significant changes in gene expression. TSSs that were identified as present in one condition or both will be an important resource for future studies of these, and possibly other, Alphaproteobacteria.

Graph Splitting: A Graph-Based Approach for Superfamily-Scale Phylogenetic Tree Reconstruction

A protein superfamily contains distantly related proteins that have acquired diverse biological functions through a long evolutionary history. Phylogenetic analysis of the early evolution of protein superfamilies is a key challenge because existing phylogenetic methods show poor performance when protein sequences are too diverged to construct an informative multiple sequence alignment (MSA). Here, we propose the Graph Splitting (GS) method, which rapidly reconstructs a protein superfamily-scale phylogenetic tree using a graph-based approach. Evolutionary simulation showed that the GS method can accurately reconstruct phylogenetic trees and be robust to major problems in phylogenetic estimation, such as biased taxon sampling, heterogeneous evolutionary rates, and long-branch attraction when sequences are substantially diverge. Its application to an empirical data set of the triosephosphate isomerase (TIM)-barrel superfamily suggests rapid evolution of protein-mediated pyrimidine biosynthesis, likely taking place after the RNA world. Furthermore, the GS method can also substantially improve performance of widely used MSA methods by providing accurate guide trees.

Shedding light on a Group IV (ECF11) alternative σ factor

This year marks the 50th anniversary of the discovery of σ70 as a protein factor that was needed for bacterial RNA polymerase to accurately transcribe a promoter in vitro. It was 25 years later that the Group IV alternative σs were described as a distinct family of proteins related to σ70 . In the intervening time, there has been an ever-growing list of Group IV σs, numbers of genes they transcribe, insight into the diverse suite of processes they control, and appreciation for their impact on bacterial lifestyles. This work summarizes knowledge of the Rhodobacter sphaeroides σE -ChrR pair, a member of the ECF11 subfamily of Group IV alternative σs, in protecting cells from the reactive oxygen species, singlet oxygen. It describes lessons learned from analyzing ChrR, a zinc-dependent anti-σ factor, that are generally applicable to Group IV σs and relevant to the response to single oxygen. This MicroReview also illustrates insights into stress responses in this and other bacteria that have been acquired by analyzing or modeling the activity of the σE -ChrR across the bacterial phylogeny.

Tracing the phylogenetic history of the Crl regulon through the Bacteria and Archaea genomes

Background

Crl, identified for curli production, is a small transcription factor that stimulates the association of the σ^S factor (RpoS) with the RNA polymerase core through direct and specific interactions, increasing the transcription rate of genes during the transition from exponential to stationary phase at low temperatures, using indole as an effector molecule. The lack of a comprehensive collection of information on the Crl regulon makes it difficult to identify a dominant function of Crl and to generate any hypotheses concerning its taxonomical distribution in archaeal and bacterial organisms.

Results

In this work, based on a systematic literature review, we identified the first comprehensive dataset of 86 genes under the control of Crl in the bacterium Escherichia coli K-12; those genes correspond to 40% of the σ^S regulon in this bacterium. Based on an analysis of orthologs in 18 archaeal and 69 bacterial taxonomical divisions and using E. coli K-12 as a framework, we suggest three main events that resulted in this regulon’s actual form: (i) in a first step, rpoS, a gene widely distributed in bacteria and archaea cellular domains, was recruited to regulate genes involved in ancient metabolic processes, such as those associated with glycolysis and the tricarboxylic acid cycle; (ii) in a second step, the regulon recruited those genes involved in metabolic processes, which are mainly taxonomically constrained to Proteobacteria, with some secondary losses, such as those genes involved in responses to stress or starvation and cell adhesion, among others; and (iii) in a posterior step, Crl might have been recruited in Enterobacteriaceae; because its taxonomical pattern constrained to this bacterial order, however further analysis are necessary.

Conclusions

Therefore, we suggest that the regulon Crl is highly flexible for phenotypic adaptation, probably as consequence of the diverse growth environments associated with all organisms in which members of this regulatory network are present.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5619-z) contains supplementary material, which is available to authorized users.

In silico simulations of occurrence of transcription factor binding sites in bacterial genomes

Background

Interactions between transcription factors and their specific binding sites are a key component of regulation of gene expression. Until recently, it was generally assumed that most bacterial transcription factor binding sites are located at or near promoters. However, several recent works utilizing high-throughput technology to detect transcription factor binding sites in bacterial genomes found a large number of binding sites in unexpected locations, particularly inside genes, as opposed to known or expected promoter regions. While some of these intragenic binding sites likely have regulatory functions, an alternative scenario is that many of these binding sites arise by chance in the absence of selective constraints. The latter possibility was supported by in silico simulations for σ⁵⁴ binding sites in Salmonella.

Results

In this work, we extend these simulations to more than forty transcription factors from E. coli and other bacteria. The results suggest that binding sites for all analyzed transcription factors are likely to arise throughout the genome by random genetic drift and many transcription factor binding sites found in genomes may not have specific regulatory functions. In addition, when comparing observed and expected patterns of occurrence of binding sites in genomes, we observed distinct differences among different transcription factors.

Conclusions

We speculate that transcription factor binding sites randomly occurring throughout the genome could be beneficial in promoting emergence of new regulatory interactions and thus facilitating evolution of gene regulatory networks.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1381-8) contains supplementary material, which is available to authorized users.

Reassessing the Structure and Function Relationship of the O2 Sensing Transcription Factor FNR

SIGNIFICANCE

The Escherichia coli regulatory protein fumarate nitrate reduction (FNR) mediates a global transcriptional response upon O₂ deprivation. Spanning nearly 40 years of research investigations, our understanding of how FNR senses and responds to O₂ has considerably progressed despite a lack of structural information for most of that period. This knowledge has established the paradigm for how facultative anaerobic bacteria sense changes in O₂ tension. Recent Advances: Recently, the X-ray crystal structure of Aliivibrio fischeri FNR with its [4Fe-4S] cluster cofactor was solved and has provided valuable new insight into FNR structure and function. These findings have alluded to the conformational changes that may occur to alter FNR activity in response to O₂.

CRITICAL ISSUES

Here, we review the major features of this structure in context of previously acquired data. In doing so, we discuss additional mechanistic aspects of FNR function that warrant further investigation.

FUTURE DIRECTIONS

To complement the [4Fe-4S]-FNR structure, the structures of apo-FNR and FNR bound to DNA or RNA polymerase are needed. Together, these structures would elevate our understanding of how ligation of its [4Fe-4S] cluster allows FNR to regulate transcription according to the level of environmental O₂.

Tuning the modular Paracoccus denitrificans respirome to adapt from aerobic respiration to anaerobic denitrification

Bacterial denitrification is a respiratory process that is a major source and sink of the potent greenhouse gas nitrous oxide. Many denitrifying bacteria can adjust to life in both oxic and anoxic environments through differential expression of their respiromes in response to environmental signals such as oxygen, nitrate and nitric oxide. We used steady-state oxic and anoxic chemostat cultures to demonstrate that the switch from aerobic to anaerobic metabolism is brought about by changes in the levels of expression of relatively few genes, but this is sufficient to adjust the configuration of the respirome to allow the organism to efficiently respire nitrate without the significant release of intermediates, such as nitrous oxide. The regulation of the denitrification respirome in strains deficient in the transcription factors FnrP, Nnr and NarR was explored and revealed that these have both inducer and repressor activities, possibly due to competitive binding at similar DNA binding sites. This may contribute to the fine tuning of expression of the denitrification respirome and so adds to the understanding of the regulation of nitrous oxide emission by denitrifying bacteria in response to different environmental signals.

Evolutionary Analysis of HIV-1 Pol Proteins Reveals Representative Residues for Viral Subtype Differentiation

RNA viruses have been used as model systems to understand the patterns and processes of molecular evolution because they have high mutation rates and are genetically diverse. Human immunodeficiency virus 1 (HIV-1), the etiological agent of acquired immune deficiency syndrome, is highly genetically diverse, and is classified into several groups and subtypes. However, it has been difficult to use its diverse sequences to establish the overall phylogenetic relationships of different strains or the trends in sequence conservation with the construction of phylogenetic trees. Our aims were to systematically characterize HIV-1 subtype evolution and to identify the regions responsible for HIV-1 subtype differentiation at the amino acid level in the Pol protein, which is often used to classify the HIV-1 subtypes. In this study, we systematically characterized the mutation sites in 2,052 Pol proteins from HIV-1 group M (144 subtype A; 1,528 subtype B; 380 subtype C), using sequence similarity networks. We also used spectral clustering to group the sequences based on the network graph structures. A stepwise analysis of the cluster hierarchies allowed us to estimate a possible evolutionary pathway for the Pol proteins. The subtype A sequences also clustered according to when and where the viruses were isolated, whereas both the subtype B and C sequences remained as single clusters. Because the Pol protein has several functional domains, we identified the regions that are discriminative by comparing the structures of the domain-based networks. Our results suggest that sequence changes in the RNase H domain and the reverse transcriptase (RT) connection domain are responsible for the subtype classification. By analyzing the different amino acid compositions at each site in both domain sequences, we found that a few specific amino acid residues (i.e., M357 in the RT connection domain and Q480, Y483, and L491 in the RNase H domain) represent the differences among the subtypes. These residues were located on the surface of the RT structure and in the vicinity of the amino acid sites responsible for RT enzymatic activity or function.

FixK2 Is the Main Transcriptional Activator of Bradyrhizobium diazoefficiens nosRZDYFLX Genes in Response to Low Oxygen

The powerful greenhouse gas, nitrous oxide (N₂O) has a strong potential to drive climate change. Soils are the major source of N₂O and microbial nitrification and denitrification the main processes involved. The soybean endosymbiont Bradyrhizobium diazoefficiens is considered a model to study rhizobial denitrification, which depends on the napEDABC, nirK, norCBQD, and nosRZDYFLX genes. In this bacterium, the role of the regulatory cascade FixLJ-FixK₂-NnrR in the expression of napEDABC, nirK, and norCBQD genes involved in N₂O synthesis has been previously unraveled. However, much remains to be discovered regarding the regulation of the respiratory N₂O reductase (N₂OR), the key enzyme that mitigates N₂O emissions. In this work, we have demonstrated that nosRZDYFLX genes constitute an operon which is transcribed from a major promoter located upstream of the nosR gene. Low oxygen was shown to be the main inducer of expression of nosRZDYFLX genes and N₂OR activity, FixK₂ being the regulatory protein involved in such control. Further, by using an in vitro transcription assay with purified FixK₂ protein and B. diazoefficiens RNA polymerase we were able to show that the nosRZDYFLX genes are direct targets of FixK₂.

FnrL and Three Dnr Regulators Are Used for the Metabolic Adaptation to Low Oxygen Tension in Dinoroseobacter shibae

The heterotrophic marine bacterium Dinoroseobacter shibae utilizes aerobic respiration and anaerobic denitrification supplemented with aerobic anoxygenic photosynthesis for energy generation. The aerobic to anaerobic transition is controlled by four Fnr/Crp family regulators in a unique cascade-type regulatory network. FnrL is utilizing an oxygen-sensitive Fe-S cluster for oxygen sensing. Active FnrL is inducing most operons encoding the denitrification machinery and the corresponding heme biosynthesis. Activation of gene expression of the high oxygen affinity cbb3-type and repression of the low affinity aa3-type cytochrome c oxidase is mediated by FnrL. Five regulator genes including dnrE and dnrF are directly controlled by FnrL. Multiple genes of the universal stress protein (USP) and cold shock response are further FnrL targets. DnrD, most likely sensing NO via a heme cofactor, co-induces genes of denitrification, heme biosynthesis, and the regulator genes dnrE and dnrF. DnrE is controlling genes for a putative Na+/H+ antiporter, indicating a potential role of a Na+ gradient under anaerobic conditions. The formation of the electron donating primary dehydrogenases is coordinated by FnrL and DnrE. Many plasmid encoded genes were DnrE regulated. DnrF is controlling directly two regulator genes including the Fe-S cluster biosynthesis regulator iscR, genes of the electron transport chain and the glutathione metabolism. The genes for nitrate reductase and CO dehydrogenase are repressed by DnrD and DnrF. Both regulators in concert with FnrL are inducing the photosynthesis genes. One of the major denitrification operon control regions, the intergenic region between nirS and nosR2, contains one Fnr/Dnr binding site. Using regulator gene mutant strains, lacZ-reporter gene fusions in combination with promoter mutagenesis, the function of the single Fnr/Dnr binding site for FnrL-, DnrD-, and partly DnrF-dependent nirS and nosR2 transcriptional activation was shown. Overall, the unique regulatory network of the marine bacterium D. shibae for the transition from aerobic to anaerobic growth composed of four Crp/Fnr family regulators was elucidated.

Assessment of transfer methods for comparative genomics of regulatory networks in bacteria

Background

Comparative genomics can leverage the vast amount of available genomic sequences to reconstruct and analyze transcriptional regulatory networks in Bacteria, but the efficacy of this approach hinges on the ability to transfer regulatory network information from reference species to the genomes under analysis. Several methods have been proposed to transfer regulatory information between bacterial species, but the paucity and distributed nature of experimental information on bacterial transcriptional networks have prevented their systematic evaluation.

Results

We report the compilation of a large catalog of transcription factor-binding sites across Bacteria and its use to systematically benchmark proposed transfer methods across pairs of bacterial species. We evaluate motif- and accuracy-based metrics to assess the results of regulatory network transfer and we identify the precision-recall area-under-the-curve as the best metric for this purpose due to the large class-imbalanced nature of the problem. Methods assuming conservation of the transcription factor-binding motif (motif-based) are shown to substantially outperform those assuming conservation of regulon composition (network-based), even though their efficiency can decrease sharply with increasing phylogenetic distance. Variations of the basic motif-based transfer method do not yield significant improvements in transfer accuracy. Our results indicate that detection of a large enough number of regulated orthologs is critical for network-based transfer methods, but that relaxing orthology requirements does not improve results. Using the transcriptional regulators LexA and Fur as case examples, we also show how DNA-binding domain sequence similarity can yield confounding results as an indicator of transfer efficiency for motif-based methods.

Conclusions

Counter to standard practice, our evaluation of metrics to assess the efficiency of methods for regulatory network information transfer reveals that the area under precision-recall (PR) curves is a more precise and informative metric than that of receiver-operating-characteristic (ROC) curves, confirming similar findings in other class-imbalanced settings. Our systematic assessment of transfer methods reveals that simple approaches to both motif- and network-based transfer of regulatory information provide equal or better results than more elaborate methods. We also show that there are not effective predictors of transfer efficacy, substantiating the long-standing practice of manual curation in comparative genomics analyses.

Electronic supplementary material

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

Nitrogen Fixation and Molecular Oxygen: Comparative Genomic Reconstruction of Transcription Regulation in Alphaproteobacteria

Biological nitrogen fixation plays a crucial role in the nitrogen cycle. An ability to fix atmospheric nitrogen, reducing it to ammonium, was described for multiple species of Bacteria and Archaea. The transcriptional regulatory network for nitrogen fixation was extensively studied in several representatives of the class Alphaproteobacteria. This regulatory network includes the activator of nitrogen fixation NifA, working in tandem with the alternative sigma-factor RpoN as well as oxygen-responsive regulatory systems, one-component regulators FnrN/FixK and two-component system FixLJ. Here we used a comparative genomics approach for in silico study of the transcriptional regulatory network in 50 genomes of Alphaproteobacteria. We extended the known regulons and proposed the scenario for the evolution of the nitrogen fixation transcriptional network. The reconstructed network substantially expands the existing knowledge of transcriptional regulation in nitrogen-fixing microorganisms and can be used for genetic experiments, metabolic reconstruction, and evolutionary analysis.

Transcriptional Regulation of the Outer Membrane Porin Gene ompW Reveals its Physiological Role during the Transition from the Aerobic to the Anaerobic Lifestyle of Escherichia coli

Understanding bacterial physiology relies on elucidating the regulatory mechanisms and cellular functions of those differentially expressed genes in response to environmental changes. A widespread Gram-negative bacterial outer membrane protein OmpW has been implicated in the adaptation to stresses in various species. It is recently found to be present in the regulon of the global anaerobic transcription factor FNR and ArcA in Escherichia coli. However, little is known about the physiological implications of this regulatory disposition. In this study, we demonstrate that transcription of ompW is indeed mediated by a series of global regulators involved in the anaerobiosis of E. coli. We show that FNR can both activate and repress the expression of ompW through its direct binding to two distinctive sites, -81.5 and -126.5 bp respectively, on ompW promoter. ArcA also participates in repression of ompW under anaerobic condition, but in an FNR dependent manner. Additionally, ompW is also subject to the regulation by CRP and NarL which senses the availability and types of carbon sources and respiration electron acceptors in the environment respectively, implying a role of OmpW in the carbon and energy metabolism of E. coli during its anaerobic adaptation. Molecular docking reveals that OmpW can bind fumarate, an alternative electron acceptor in anaerobic respiration, with sufficient affinity. Moreover, supplement of fumarate or succinate which belongs to the C₄-dicarboxylates family of metabolite, to E. coli culture rescues OmpW-mediated colicin S4 killing. Taken together, we propose that OmpW is involved in anaerobic carbon and energy metabolism to mediate the transition from aerobic to anaerobic lifestyle in E. coli.

Different Functions of Phylogenetically Distinct Bacterial Complex I Isozymes

NADH:quinone oxidoreductase (complex I) is a bioenergetic enzyme that transfers electrons from NADH to quinone, conserving the energy of this reaction by contributing to the proton motive force. While the importance of NADH oxidation to mitochondrial aerobic respiration is well documented, the contribution of complex I to bacterial electron transport chains has been tested in only a few species. Here, we analyze the function of two phylogenetically distinct complex I isozymes in Rhodobacter sphaeroides, an alphaproteobacterium that contains well-characterized electron transport chains. We found that R. sphaeroides complex I activity is important for aerobic respiration and required for anaerobic dimethyl sulfoxide (DMSO) respiration (in the absence of light), photoautotrophic growth, and photoheterotrophic growth (in the absence of an external electron acceptor). Our data also provide insight into the functions of the phylogenetically distinct R. sphaeroidescomplex I enzymes (complex IA and complex IE) in maintaining a cellular redox state during photoheterotrophic growth. We propose that the function of each isozyme during photoheterotrophic growth is either NADH synthesis (complex IA) or NADH oxidation (complex IE). The canonical alphaproteobacterial complex I isozyme (complex IA) was also shown to be important for routing electrons to nitrogenase-mediated H2 production, while the horizontally acquired enzyme (complex IE) was dispensable in this process. Unlike the singular role of complex I in mitochondria, we predict that the phylogenetically distinct complex I enzymes found across bacterial species have evolved to enhance the functions of their respective electron transport chains.

IMPORTANCE

Cells use a proton motive force (PMF), NADH, and ATP to support numerous processes. In mitochondria, complex I uses NADH oxidation to generate a PMF, which can drive ATP synthesis. This study analyzed the function of complex I in bacteria, which contain more-diverse and more-flexible electron transport chains than mitochondria. We tested complex I function in Rhodobacter sphaeroides, a bacterium predicted to encode two phylogenetically distinct complex I isozymes. R. sphaeroides cells lacking both isozymes had growth defects during all tested modes of growth, illustrating the important function of this enzyme under diverse conditions. We conclude that the two isozymes are not functionally redundant and predict that phylogenetically distinct complex I enzymes have evolved to support the diverse lifestyles of bacteria.

Global transcriptional analysis of Escherichia coli expressing IrrE, a regulator from Deinococcus radiodurans, in response to NaCl shock

Improving the microbial tolerance to stresses is very important for bioprocesses. Our previous study showed that IrrE, a global regulator from the extremely radioresistant bacterium Deinococcus radiodurans, dramatically enhanced the multi-stress tolerance of Escherichia coli when expressed exogenously. However, the function of IrrE is still unclear. In this study, we used whole-genome microarray assays to profile the global gene expression of the IrrE-expressing E. coli strain MGE and the control strain MGT with or without salt shock. The analysis showed that IrrE expression led to many differentially expressed genes in E. coli, which were responsible for the transport and metabolism of trehalose and glycerol, nucleotide biosynthesis, carbon source utilization, amino acid utilization, and acid resistance, including many RpoS-dependent genes, e.g., the trehalose biosynthesis genes otsAB, the acid-resistance genes gadABC and uspB, the osmotic and oxidative stress response genes katE (response to DNA damage stimulus and stress) and osmBC (response to stress), and gadWX (which controls the transcription of pH-inducible genes). The intracellular content of trehalose and glycerol increased significantly in the IrrE-expressing strain after NaCl treatment for 0 and 60 min as determined by HPLC. These results indicated the possibility that IrrE regulates the global regulator RpoS. Interestingly, we found that although IrrE did not affect the level of the rpoS transcript, it enhanced the accumulation of the RpoS protein by increasing the expression of the antiadaptors, AppY, IraM and IraD, which inhibit RpoS degradation, suggesting that the accumulation of RpoS due to IrrE regulation is an important way to improve tolerance to salt and other stresses in E. coli.

Oxygen-dependent regulation of bacterial lipid production

Understanding the mechanisms of lipid accumulation in microorganisms is important for several reasons. In addition to providing insight into assembly of biological membranes, lipid accumulation has important applications in the production of renewable fuels and chemicals. The photosynthetic bacterium Rhodobacter sphaeroides is an attractive organism to study lipid accumulation, as it has the ability to increase membrane production at low O2 tensions. Under these conditions, R. sphaeroides develops invaginations of the cytoplasmic membrane to increase its membrane surface area for housing of the membrane-bound components of its photosynthetic apparatus. Here we use fatty acid levels as a reporter of membrane lipid content. We show that, under low-O2 and anaerobic conditions, the total fatty acid content per cell increases 3-fold. We also find that the increases in the amount of fatty acid and photosynthetic pigment per cell are correlated as O2 tensions or light intensity are changed. To ask if lipid and pigment accumulation were genetically separable, we analyzed strains with mutations in known photosynthetic regulatory pathways. While a strain lacking AppA failed to induce photosynthetic pigment-protein complex accumulation, it increased fatty acid content under low-O2 conditions. We also found that an intact PrrBA pathway is required for low-O2-induced fatty acid accumulation. Our findings suggest a previously unknown role of R. sphaeroides transcriptional regulators in increasing fatty acid and phospholipid accumulation in response to decreased O2 tension.

IMPORTANCE

Lipids serve important functions in living systems, either as structural components of membranes or as a form of carbon storage. Understanding the mechanisms of lipid accumulation in microorganisms is important for providing insight into the assembly of biological membranes and additionally has important applications in the production of renewable fuels and chemicals. In this study, we investigate the ability of Rhodobacter sphaeroides to increase membrane production at low O2 tensions in order to house its photosynthetic apparatus. We demonstrate that this bacterium has a mechanism to increase lipid content in response to decreased O2 tension and identify a transcription factor necessary for this response. This is significant because it identifies a transcriptional regulatory pathway that can increase microbial lipid content.

An integrated approach to reconstructing genome-scale transcriptional regulatory networks

Transcriptional regulatory networks (TRNs) program cells to dynamically alter their gene expression in response to changing internal or environmental conditions. In this study, we develop a novel workflow for generating large-scale TRN models that integrates comparative genomics data, global gene expression analyses, and intrinsic properties of transcription factors (TFs). An assessment of this workflow using benchmark datasets for the well-studied γ-proteobacterium Escherichia coli showed that it outperforms expression-based inference approaches, having a significantly larger area under the precision-recall curve. Further analysis indicated that this integrated workflow captures different aspects of the E. coli TRN than expression-based approaches, potentially making them highly complementary. We leveraged this new workflow and observations to build a large-scale TRN model for the α-Proteobacterium Rhodobacter sphaeroides that comprises 120 gene clusters, 1211 genes (including 93 TFs), 1858 predicted protein-DNA interactions and 76 DNA binding motifs. We found that ~67% of the predicted gene clusters in this TRN are enriched for functions ranging from photosynthesis or central carbon metabolism to environmental stress responses. We also found that members of many of the predicted gene clusters were consistent with prior knowledge in R. sphaeroides and/or other bacteria. Experimental validation of predictions from this R. sphaeroides TRN model showed that high precision and recall was also obtained for TFs involved in photosynthesis (PpsR), carbon metabolism (RSP_0489) and iron homeostasis (RSP_3341). In addition, this integrative approach enabled generation of TRNs with increased information content relative to R. sphaeroides TRN models built via other approaches. We also show how this approach can be used to simultaneously produce TRN models for each related organism used in the comparative genomics analysis. Our results highlight the advantages of integrating comparative genomics of closely related organisms with gene expression data to assemble large-scale TRN models with high-quality predictions.

The ever growing amount of genomic data enables the assembly of large-scale network models that can provide important new insights into living systems. However, assembly and validation of such large-scale models can be challenging, since we often lack sufficient information to make accurate predictions. This work describes a new approach for constructing large-scale transcriptional regulatory networks of individual cells. We show that the reconstructed network captures a significantly larger fraction of cellular regulatory processes than networks generated by other existing approaches. We predict this approach, with appropriate refinements, will allow reconstruction of large-scale transcriptional network models for a variety of other organisms. As we work towards modeling the function of cells or complex ecosystems, individually reconstructed network models of signaling, information transfer and metabolism, can be integrated to provide high information predictions and insights not otherwise obtainable.

Global analysis of photosynthesis transcriptional regulatory networks

Photosynthesis is a crucial biological process that depends on the interplay of many components. This work analyzed the gene targets for 4 transcription factors: FnrL, PrrA, CrpK and MppG (RSP_2888), which are known or predicted to control photosynthesis in Rhodobacter sphaeroides. Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) identified 52 operons under direct control of FnrL, illustrating its regulatory role in photosynthesis, iron homeostasis, nitrogen metabolism and regulation of sRNA synthesis. Using global gene expression analysis combined with ChIP-seq, we mapped the regulons of PrrA, CrpK and MppG. PrrA regulates ∼34 operons encoding mainly photosynthesis and electron transport functions, while CrpK, a previously uncharacterized Crp-family protein, regulates genes involved in photosynthesis and maintenance of iron homeostasis. Furthermore, CrpK and FnrL share similar DNA binding determinants, possibly explaining our observation of the ability of CrpK to partially compensate for the growth defects of a ΔFnrL mutant. We show that the Rrf2 family protein, MppG, plays an important role in photopigment biosynthesis, as part of an incoherent feed-forward loop with PrrA. Our results reveal a previously unrealized, high degree of combinatorial regulation of photosynthetic genes and significant cross-talk between their transcriptional regulators, while illustrating previously unidentified links between photosynthesis and the maintenance of iron homeostasis.

Photosynthetic organisms are among the most abundant life forms on earth. Their unique ability to harvest solar energy and use it to fix atmospheric carbon dioxide is at the foundation of the global food chain. This paper reports the first comprehensive analysis of networks that control expression of photosynthesis genes using Rhodobacter sphaeroides, a microbe that has been studied for decades as a model of solar energy capture and other aspects of the photosynthetic lifestyle. We find a previously unappreciated complexity in the level of control of photosynthetic genes, while identifying new links between photosynthesis and central processes like iron availability. This organism is an ancestor of modern day plants, so our data can inform studies in other photosynthetic organisms and improve our ability to harness solar energy for food and industrial processes.

Fe-S proteins that regulate gene expression

Iron-sulfur (Fe-S) cluster containing proteins that regulate gene expression are present in most organisms. The innate chemistry of their Fe-S cofactors makes these regulatory proteins ideal for sensing environmental signals, such as gases (e.g. O2 and NO), levels of Fe and Fe-S clusters, reactive oxygen species, and redox cycling compounds, to subsequently mediate an adaptive response. Here we review the recent findings that have provided invaluable insight into the mechanism and function of these highly significant Fe-S regulatory proteins. This article is part of a Special Issue entitled: Fe/S proteins: Analysis, structure, function, biogenesis and diseases.

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.

Construction and verification of the transcriptional regulatory response network of Streptococcus mutans upon treatment with the biofilm inhibitor carolacton

Background

Carolacton is a newly identified secondary metabolite causing altered cell morphology and death of Streptococcus mutans biofilm cells. To unravel key regulators mediating these effects, the transcriptional regulatory response network of S. mutans biofilms upon carolacton treatment was constructed and analyzed. A systems biological approach integrating time-resolved transcriptomic data, reverse engineering, transcription factor binding sites, and experimental validation was carried out.

Results

The co-expression response network constructed from transcriptomic data using the reverse engineering algorithm called the Trend Correlation method consisted of 8284 gene pairs. The regulatory response network inferred by superimposing transcription factor binding site information into the co-expression network comprised 329 putative transcriptional regulatory interactions and could be classified into 27 sub-networks each co-regulated by a transcription factor. These sub-networks were significantly enriched with genes sharing common functions. The regulatory response network displayed global hierarchy and network motifs as observed in model organisms. The sub-networks modulated by the pyrimidine biosynthesis regulator PyrR, the glutamine synthetase repressor GlnR, the cysteine metabolism regulator CysR, global regulators CcpA and CodY and the two component system response regulators VicR and MbrC among others could putatively be related to the physiological effect of carolacton. The predicted interactions from the regulatory network between MbrC, known to be involved in cell envelope stress response, and the murMN-SMU_718c genes encoding peptidoglycan biosynthetic enzymes were experimentally confirmed using Electro Mobility Shift Assays. Furthermore, gene deletion mutants of five predicted key regulators from the response networks were constructed and their sensitivities towards carolacton were investigated. Deletion of cysR, the node having the highest connectivity among the regulators chosen from the regulatory network, resulted in a mutant which was insensitive to carolacton thus demonstrating not only the essentiality of cysR for the response of S. mutans biofilms to carolacton but also the relevance of the predicted network.

Conclusion

The network approach used in this study revealed important regulators and interactions as part of the response mechanisms of S. mutans biofilm cells to carolacton. It also opens a door for further studies into novel drug targets against streptococci.

Electronic supplementary material

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

A Rhodobacter sphaeroides protein mechanistically similar to Escherichia coli DksA regulates photosynthetic growth

ABSTRACT DksA is a global regulatory protein that, together with the alarmone ppGpp, is required for the "stringent response" to nutrient starvation in the gammaproteobacterium Escherichia coli and for more moderate shifts between growth conditions. DksA modulates the expression of hundreds of genes, directly or indirectly. Mutants lacking a DksA homolog exhibit pleiotropic phenotypes in other gammaproteobacteria as well. Here we analyzed the DksA homolog RSP2654 in the more distantly related Rhodobacter sphaeroides, an alphaproteobacterium. RSP2654 is 42% identical and similar in length to E. coli DksA but lacks the Zn finger motif of the E. coli DksA globular domain. Deletion of the RSP2654 gene results in defects in photosynthetic growth, impaired utilization of amino acids, and an increase in fatty acid content. RSP2654 complements the growth and regulatory defects of an E. coli strain lacking the dksA gene and modulates transcription in vitro with E. coli RNA polymerase (RNAP) similarly to E. coli DksA. RSP2654 reduces RNAP-promoter complex stability in vitro with RNAPs from E. coli or R. sphaeroides, alone and synergistically with ppGpp, suggesting that even though it has limited sequence identity to E. coli DksA (DksAEc), it functions in a mechanistically similar manner. We therefore designate the RSP2654 protein DksARsp. Our work suggests that DksARsp has distinct and important physiological roles in alphaproteobacteria and will be useful for understanding structure-function relationships in DksA and the mechanism of synergy between DksA and ppGpp. IMPORTANCE The role of DksA has been analyzed primarily in the gammaproteobacteria, in which it is best understood for its role in control of the synthesis of the translation apparatus and amino acid biosynthesis. Our work suggests that DksA plays distinct and important physiological roles in alphaproteobacteria, including the control of photosynthesis in Rhodobacter sphaeroides. The study of DksARsp, should be useful for understanding structure-function relationships in the protein, including those that play a role in the little-understood synergy between DksA and ppGpp.

Comprehensive computational analysis of bacterial CRP/FNR superfamily and its target motifs reveals stepwise evolution of transcriptional networks

The cAMP receptor protein (CRP)/fumarate and nitrate reduction regulatory protein (FNR)-type transcription factors (TFs) are members of a well-characterized global TF family in bacteria and have two conserved domains: the N-terminal ligand-binding domain for small molecules (e.g., cAMP, NO, or O₂) and the C-terminal DNA-binding domain. Although the CRP/FNR-type TFs recognize very similar consensus DNA target sequences, they can regulate different sets of genes in response to environmental signals. To clarify the evolution of the CRP/FNR-type TFs throughout the bacterial kingdom, we undertook a comprehensive computational analysis of a large number of annotated CRP/FNR-type TFs and the corresponding bacterial genomes. Based on the amino acid sequence similarities among 1,455 annotated CRP/FNR-type TFs, spectral clustering classified the TFs into 12 representative groups, and stepwise clustering allowed us to propose a possible process of protein evolution. Although each cluster mainly consists of functionally distinct members (e.g., CRP, NTC, FNR-like protein, and FixK), FNR-related TFs are found in several groups and are distributed in a wide range of bacterial phyla in the sequence similarity network. This result suggests that the CRP/FNR-type TFs originated from an ancestral FNR protein, involved in nitrogen fixation. Furthermore, a phylogenetic profiling analysis showed that combinations of TFs and their target genes have fluctuated dynamically during bacterial evolution. A genome-wide analysis of TF-binding sites also suggested that the diversity of the transcriptional regulatory system was derived by the stepwise adaptation of TF-binding sites to the evolution of TFs.

Piecing Together Cell-like Systems

Several laboratories are pursuing the synthesis of cellular systems from different directions, including those that begin with simple chemicals to those that exploit existing cells. The methods that begin with nonliving components tend to focus on mimicking specific features of life, such as genomic replication, protein synthesis, sensory systems, and compartment formation, growth, and division. Conversely, the more prevalent synthetic biology approaches begin with something that is already alive and seek to impart new behavior on existing cells. Here we discuss advances in building cell-like systems that mimic key features of life with defined components.

How life changes itself: the Read-Write (RW) genome

The genome has traditionally been treated as a Read-Only Memory (ROM) subject to change by copying errors and accidents. In this review, I propose that we need to change that perspective and understand the genome as an intricately formatted Read-Write (RW) data storage system constantly subject to cellular modifications and inscriptions. Cells operate under changing conditions and are continually modifying themselves by genome inscriptions. These inscriptions occur over three distinct time-scales (cell reproduction, multicellular development and evolutionary change) and involve a variety of different processes at each time scale (forming nucleoprotein complexes, epigenetic formatting and changes in DNA sequence structure). Research dating back to the 1930s has shown that genetic change is the result of cell-mediated processes, not simply accidents or damage to the DNA. This cell-active view of genome change applies to all scales of DNA sequence variation, from point mutations to large-scale genome rearrangements and whole genome duplications (WGDs). This conceptual change to active cell inscriptions controlling RW genome functions has profound implications for all areas of the life sciences.

The structure of Bradyrhizobium japonicum transcription factor FixK2 unveils sites of DNA binding and oxidation

FixK2 is a regulatory protein that activates a large number of genes for the anoxic and microoxic, endosymbiotic, and nitrogen-fixing life styles of the α-proteobacterium Bradyrhizobium japonicum. FixK2 belongs to the cAMP receptor protein (CRP) superfamily. Although most CRP family members are coregulated by effector molecules, the activity of FixK2 is negatively controlled by oxidation of its single cysteine (Cys-183) located next to the DNA-binding domain and possibly also by proteolysis. Here, we report the three-dimensional x-ray structure of FixK2, a representative of the FixK subgroup of the CRP superfamily. Crystallization succeeded only when (i) an oxidation- and protease-insensitive protein variant (FixK2(C183S)-His6) was used in which Cys-183 was replaced with serine and the C terminus was fused with a hexahistidine tag and (ii) this protein was allowed to form a complex with a 30-mer double-stranded target DNA. The structure of the FixK2-DNA complex was solved at a resolution of 1.77 Å, at which the protein formed a homodimer. The precise protein-DNA contacts were identified, which led to an affirmation of the canonical target sequence, the so-called FixK2 box. The C terminus is surface-exposed, which might explain its sensitivity to specific cleavage and degradation. The oxidation-sensitive Cys-183 is also surface-exposed and in close proximity to DNA. Therefore, we propose a mechanism whereby the oxo acids generated after oxidation of the cysteine thiol cause an electrostatic repulsion, thus preventing specific DNA binding.

Signal correlations in ecological niches can shape the organization and evolution of bacterial gene regulatory networks

Transcriptional regulation plays a significant role in the biological response of bacteria to changing environmental conditions. Therefore, mapping transcriptional regulatory networks is an important step not only in understanding how bacteria sense and interpret their environment but also to identify the functions involved in biological responses to specific conditions. Recent experimental and computational developments have facilitated the characterization of regulatory networks on a genome-wide scale in model organisms. In addition, the multiplication of complete genome sequences has encouraged comparative analyses to detect conserved regulatory elements and infer regulatory networks in other less well-studied organisms. However, transcription regulation appears to evolve rapidly, thus, creating challenges for the transfer of knowledge to nonmodel organisms. Nevertheless, the mechanisms and constraints driving the evolution of regulatory networks have been the subjects of numerous analyses, and several models have been proposed. Overall, the contributions of mutations, recombination, and horizontal gene transfer are complex. Finally, the rapid evolution of regulatory networks plays a significant role in the remarkable capacity of bacteria to adapt to new or changing environments. Conversely, the characteristics of environmental niches determine the selective pressures and can shape the structure of regulatory network accordingly.

Proteins needed to activate a transcriptional response to the reactive oxygen species singlet oxygen

Singlet oxygen ((1)O(2)) is a reactive oxygen species generated by energy transfer from one or more excited donors to molecular oxygen. Many biomolecules are prone to oxidation by (1)O(2), and cells have evolved systems to protect themselves from damage caused by this compound. One way that the photosynthetic bacterium Rhodobacter sphaeroides protects itself from (1)O(2) is by inducing a transcriptional response controlled by ChrR, an anti-σ factor which releases an alternative sigma factor, σ(E), in the presence of (1)O(2). Here we report that induction of σ(E)-dependent gene transcription is decreased in the presence of (1)O(2) when two conserved genes in the σ(E) regulon are deleted, including one encoding a cyclopropane fatty acid synthase homologue (RSP2144) or one encoding a protein of unknown function (RSP1091). Thus, we conclude that RSP2144 and RSP1091 are each necessary to increase σ(E) activity in the presence of (1)O(2). In addition, we found that unlike in wild-type cells, where ChrR is rapidly degraded when (1)O(2) is generated, turnover of this anti-σ factor is slowed when cells lacking RSP2144, RSP1091, or both of these proteins are exposed to (1)O(2). Further, we demonstrate that the organic hydroperoxide tert-butyl hydroperoxide promotes ChrR turnover in both wild-type cells and mutants lacking RSP2144 or RSP1091, suggesting differences in the ways different types of oxidants increase σ(E) activity.

IMPORTANCE

Oxygen serves many crucial functions on Earth; it is produced during photosynthesis and needed for other pathways. While oxygen is relatively inert, it can be converted to reactive oxygen species (ROS) that destroy biomolecules, cause disease, or kill cells. When energy is transferred to oxygen, the ROS singlet oxygen is generated. To understand how singlet oxygen impacts cells, we study the stress response to this ROS in Rhodobacter sphaeroides, a bacterium that, like plants, generates this compound as a consequence of photosynthesis. This paper identifies proteins that activate a stress response to singlet oxygen and shows that they act in a specific response to this ROS. The identified proteins are found in many free-living, symbiotic, or pathogenic bacteria that can encounter singlet oxygen in nature. Thus, our findings provide new information about a stress response to a ROS of broad biological, agricultural, and biomedical importance.

Unusual heme binding properties of the dissimilative nitrate respiration regulator, a bacterial nitric oxide sensor

AIMS

In the opportunistic pathogen Pseudomonas aeruginosa, nitric oxide (NO) triggers the respiration of nitrate (denitrification), thus allowing survival in chronic infection sites as a microaerobic-anaerobic biofilm. The NO-dependent induction of denitrification is mediated by the dissimilative nitrate respiration regulator (DNR), a transcription factor forming a stable complex with heme, which is required to sense the physiological messenger (i.e., NO). The molecular details of NO sensing in DNR and, more in general, in this class of sensors are largely unknown, and a study aimed at integrating microbiology and biochemistry is needed.

RESULTS

Here we present a comprehensive study, including in vivo results and spectroscopy, kinetics, and protein engineering, that demonstrates the direct involvement of a histidine residue in heme iron coordination. Moreover, a peculiar phenomenon of ligand switching around heme iron, which hampers the identification of the second heme axial ligand, is also suggested. These results indicate that DNR is characterized by a remarkable flexibility in solution, as observed for other cAMP receptor protein/fumarate and nitrate reductase regulators (CRP-FNR) to which DNR belongs.

INNOVATION

The present work represents one of the few studies focused on the biochemistry of NO sensing by bacterial transcriptional regulators. The data presented demonstrate that structural plasticity of DNR is crucial for the sensing activity and confers to the protein unusual heme binding properties.

CONCLUSIONS

Protein flexibility and dynamics is a key structural feature essential to explain the evolutionary success and adaptability of CRP-FNR, and may represent a common strategy employed by heme-based redox sensors, which presents features deeply different from those of canonical hemeproteins.

Transcription factor family-based reconstruction of singleton regulons and study of the Crp/Fnr, ArsR, and GntR families in Desulfovibrionales genomes

Accurate detection of transcriptional regulatory elements is essential for high-quality genome annotation, metabolic reconstruction, and modeling of regulatory networks. We developed a computational approach for reconstruction of regulons operated by transcription factors (TFs) from large protein families and applied this novel approach to three TF families in 10 Desulfovibrionales genomes. Phylogenetic analyses of 125 regulators from the ArsR, Crp/Fnr, and GntR families revealed that 65% of these regulators (termed reference TFs) are well conserved in Desulfovibrionales, while the remaining 35% of regulators (termed singleton TFs) are species specific and show a mosaic distribution. For regulon reconstruction in the group of singleton TFs, the standard orthology-based approach was inefficient, and thus, we developed a novel approach based on the simultaneous study of all homologous TFs from the same family in a group of genomes. As a result, we identified binding for 21 singleton TFs and for all reference TFs in all three analyzed families. Within each TF family we observed structural similarities between DNA-binding motifs of different reference and singleton TFs. The collection of reconstructed regulons is available at the RegPrecise database (http://regprecise.lbl.gov/RegPrecise/Desulfovibrionales.jsp).

Convergence of the transcriptional responses to heat shock and singlet oxygen stresses

Cells often mount transcriptional responses and activate specific sets of genes in response to stress-inducing signals such as heat or reactive oxygen species. Transcription factors in the RpoH family of bacterial alternative σ factors usually control gene expression during a heat shock response. Interestingly, several α-proteobacteria possess two or more paralogs of RpoH, suggesting some functional distinction. We investigated the target promoters of Rhodobacter sphaeroides RpoH_I and RpoH_II using genome-scale data derived from gene expression profiling and the direct interactions of each protein with DNA in vivo. We found that the RpoH_I and RpoH_II regulons have both distinct and overlapping gene sets. We predicted DNA sequence elements that dictate promoter recognition specificity by each RpoH paralog. We found that several bases in the highly conserved TTG in the −35 element are important for activity with both RpoH homologs; that the T-9 position, which is over-represented in the RpoH_I promoter sequence logo, is critical for RpoH_I–dependent transcription; and that several bases in the predicted −10 element were important for activity with either RpoH_II or both RpoH homologs. Genes that are transcribed by both RpoH_I and RpoH_II are predicted to encode for functions involved in general cell maintenance. The functions specific to the RpoH_I regulon are associated with a classic heat shock response, while those specific to RpoH_II are associated with the response to the reactive oxygen species, singlet oxygen. We propose that a gene duplication event followed by changes in promoter recognition by RpoH_I and RpoH_II allowed convergence of the transcriptional responses to heat and singlet oxygen stress in R. sphaeroides and possibly other bacteria.

An important property of living systems is their ability to survive under conditions of stress such as increased temperature or the presence of reactive oxygen species. Central to the function of these stress responses are transcription factors that activate specific sets of genes needed for this response. Despite the central role of stress responses across all forms of life, the processes driving their organization and evolution across organisms are poorly understood. This paper uses genomic, computational, and mutational analyses to dissect stress responses controlled by two proteins that are each members of the RpoH family of alternative σ factors. RpoH family members usually control gene expression during a heat shock response. However, the photosynthetic bacterium Rhodobacter sphaeroides and several other α-proteobacteria possess two or more paralogs of RpoH, suggesting some functional distinction. Our findings predict that a gene duplication event followed by changes in DNA recognition by RpoH_I and RpoH_II allowed convergence of the transcriptional responses to heat and singlet oxygen stress in R. sphaeroides and possibly other bacteria. Our approach and findings should interest those studying the evolution of transcription factors or the signal transduction pathways that control stress responses.

Extracytoplasmic function σ factors of the widely distributed group ECF41 contain a fused regulatory domain

Bacteria need signal transducing systems to respond to environmental changes. Next to one- and two-component systems, alternative σ factors of the extra-cytoplasmic function (ECF) protein family represent the third fundamental mechanism of bacterial signal transduction. A comprehensive classification of these proteins identified more than 40 phylogenetically distinct groups, most of which are not experimentally investigated. Here, we present the characterization of such a group with unique features, termed ECF41. Among analyzed bacterial genomes, ECF41 σ factors are widely distributed with about 400 proteins from 10 different phyla. They lack obvious anti-σ factors that typically control activity of other ECF σ factors, but their structural genes are often predicted to be cotranscribed with carboxymuconolactone decarboxylases, oxidoreductases, or epimerases based on genomic context conservation. We demonstrate for Bacillus licheniformis and Rhodobacter sphaeroides that the corresponding genes are preceded by a highly conserved promoter motif and are the only detectable targets of ECF41-dependent gene regulation. In contrast to other ECF σ factors, proteins of group ECF41 contain a large C-terminal extension, which is crucial for σ factor activity. Our data demonstrate that ECF41 σ factors are regulated by a novel mechanism based on the presence of a fused regulatory domain.

Conservation of thiol-oxidative stress responses regulated by SigR orthologues in actinomycetes

Numerous thiol-reactive compounds cause oxidative stress where cells counteract by activation of survival strategies regulated by thiol-based sensors. In Streptomyces coelicolor, a model actinomycete, a sigma/antisigma pair SigR/RsrA controls the response to thiol-oxidative stress. To unravel its full physiological functions, chromatin immuno-precipitation combined with sequence and transcript analyses were employed to identify 108 SigR target genes in S. coelicolor and to predict orthologous regulons across actinomycetes. In addition to reported genes for thiol homeostasis, protein degradation and ribosome modulation, 64 additional operons were identified suggesting new functions of this global regulator. We demonstrate that SigR maintains the level and activity of the housekeeping sigma factor HrdB during thiol-oxidative stress, a novel strategy for stress responses. We also found that SigR defends cells against UV and thiol-reactive damages, in which repair UvrA takes a part. Using a refined SigR-binding sequence model, SigR orthologues and their targets were predicted in 42 actinomycetes. This revealed a conserved core set of SigR targets to function for thiol homeostasis, protein quality control, possible modulation of transcription and translation, flavin-mediated redox reactions, and Fe-S delivery. The composition of the SigR regulon reveals a robust conserved physiological mechanism to deal with thiol-oxidative stress from bacteria to human.

Bacterial adaptation of respiration from oxic to microoxic and anoxic conditions: redox control

Under a shortage of oxygen, bacterial growth can be faced mainly by two ATP-generating mechanisms: (i) by synthesis of specific high-affinity terminal oxidases that allow bacteria to use traces of oxygen or (ii) by utilizing other substrates as final electron acceptors such as nitrate, which can be reduced to dinitrogen gas through denitrification or to ammonium. This bacterial respiratory shift from oxic to microoxic and anoxic conditions requires a regulatory strategy which ensures that cells can sense and respond to changes in oxygen tension and to the availability of other electron acceptors. Bacteria can sense oxygen by direct interaction of this molecule with a membrane protein receptor (e.g., FixL) or by interaction with a cytoplasmic transcriptional factor (e.g., Fnr). A third type of oxygen perception is based on sensing changes in redox state of molecules within the cell. Redox-responsive regulatory systems (e.g., ArcBA, RegBA/PrrBA, RoxSR, RegSR, ActSR, ResDE, and Rex) integrate the response to multiple signals (e.g., ubiquinone, menaquinone, redox active cysteine, electron transport to terminal oxidases, and NAD/NADH) and activate or repress target genes to coordinate the adaptation of bacterial respiration from oxic to anoxic conditions. Here, we provide a compilation of the current knowledge about proteins and regulatory networks involved in the redox control of the respiratory adaptation of different bacterial species to microxic and anoxic environments.

Iron sulfur cluster proteins and microbial regulation: implications for understanding tuberculosis

All pathogenic and nonpathogenic microbes are continuously exposed to environmental or endogenous reactive oxygen and nitrogen species, which can critically effect survival and disease. Iron-sulfur [Fe-S] cluster containing prosthetic groups provide the microbial cell with a unique capacity to sense and transcriptionally respond to diatomic gases (e.g. NO and O2) and redox-cycling agents. Recent advances in our understanding of the mechanisms for how the FNR and SoxR [Fe-S] cluster proteins respond to NO and O2 have provided new insights into the biochemical mechanism of action of the Mycobacterium tuberculosis (Mtb) family of WhiB [Fe-S] cluster proteins. These insights have provided the basis for establishing a unifying paradigm for the Mtb WhiB family of proteins. Mtb is the etiological agent for tuberculosis (TB), a disease that affects nearly one-third of the world's population.

Genome-wide mapping of the binding sites of proteins that interact with DNA

The coordinated regulation of the expression of a group of genes by a specific transcription factor frequently lies at the heart of the ability of a bacterium to respond to an environmental signal, or to progress through a developmental program. Thus, in many situations, it is of interest to identify all of the genes that are under the control of a particular regulatory protein. This chapter begins with a brief overview of some of the methods that have been used in attempts to identify some or all of the members of a regulon (i.e., those genes that are the targets for a transcriptional activator or repressor). Thereafter, the chapter will focus on one technique, chromatin immunoprecipitation and microarray analysis (ChIP-chip) and some of its variants. Design considerations and some protocols for ChIP-chip experiments are provided, along with some considerations related to downstream data analysis. ChIP-chip is a method for the genome-wide localization of protein-binding sites. In a typical ChIP-chip protocol, proteins are cross-linked nonspecifically to DNA in vivo. Chromatin is extracted and sheared, and specific protein-DNA complexes are immunoprecipitated with a suitable antibody. After purification, the DNA is hybridized to a microarray (after an amplification step in some protocols), together with a differentially labeled reference sample. Features on the microarray that show an elevated fluorescence ratio reveal DNA sequences that were enriched by immunoprecipitation. The corresponding genomic locations are those that were enriched, and are therefore close to sites of binding. The use of high-density tiled microarrays allows for binding site localization with quite high resolution. It is likely that ChIP-chip will soon be superseded by ChIP-seq, in which the immunoprecipitated DNA is analyzed directly by next-generation sequencing technologies. ChIP-chip and ChIP-seq applications are not confined to regulatory proteins, since they can be used with any protein that binds to DNA, either directly, or indirectly via an interaction with another protein. Thus, ChIP-chip has been used successfully to map binding sites for nucleoid proteins, and proteins involved in DNA replication.