Molecular and structural considerations of TF-DNA binding for the generation of biologically meaningful and accurate phylogenetic footprinting analysis: the LysR-type transcriptional regulator family as a study model.

BACKGROUND: The goal of most programs developed to find transcription factor binding sites (TFBSs) is the identification of discrete sequence motifs that are significantly over-represented in a given set of sequences where a transcription factor (TF) is expected to bind. These programs assume that the nucleotide conservation of a specific motif is indicative of a selective pressure required for the recognition of a TF for its corresponding TFBS. Despite their extensive use, the accuracies reached with these programs remain low. In many cases, true TFBSs are excluded from the identification process, especially when they correspond to low-affinity but important binding sites of regulatory systems. RESULTS: We developed a computational protocol based on molecular and structural criteria to perform biologically meaningful and accurate phylogenetic footprinting analyses. Our protocol considers fundamental aspects of the TF-DNA binding process, such as: i) the active homodimeric conformations of TFs that impose symmetric structures on the TFBSs, ii) the cooperative binding of TFs, iii) the effects of the presence or absence of co-inducers, iv) the proximity between two TFBSs or one TFBS and a promoter that leads to very long spurious motifs, v) the presence of AT-rich sequences not recognized by the TF but that are required for DNA flexibility, and vi) the dynamic order in which the different binding events take place to determine a regulatory response (i.e., activation or repression). In our protocol, the abovementioned criteria were used to analyze a profile of consensus motifs generated from canonical Phylogenetic Footprinting Analyses using a set of analysis windows of incremental sizes. To evaluate the performance of our protocol, we analyzed six members of the LysR-type TF family in Gammaproteobacteria. CONCLUSIONS: The identification of TFBSs based exclusively on the significance of the over-representation of motifs in a set of sequences might lead to inaccurate results. The consideration of different molecular and structural properties of the regulatory systems benefits the identification of TFBSs and enables the development of elaborate, biologically meaningful and precise regulatory models that offer a more integrated view of the dynamics of the regulatory process of transcription.

RNA-Seq analysis of the multipartite genome of Rhizobium etli CE3 shows different replicon contributions under heat and saline shock.

BACKGROUND: Regulation of transcription is essential for any organism and Rhizobium etli (a multi-replicon, nitrogen-fixing symbiotic bacterium) is no exception. This bacterium is commonly found in the rhizosphere (free-living) or inside of root-nodules of the common bean (Phaseolus vulgaris) in a symbiotic relationship. Abiotic stresses, such as high soil temperatures and salinity, compromise the genetic stability of R. etli and therefore its symbiotic interaction with P. vulgaris. However, it is still unclear which genes are up- or down-regulated to cope with these stress conditions. The aim of this study was to identify the genes and non-coding RNAs (ncRNAs) that are differentially expressed under heat and saline shock, as well as the promoter regions of the up-regulated loci. RESULTS: Analysing the heat and saline shock responses of R. etli CE3 through RNA-Seq, we identified 756 and 392 differentially expressed genes, respectively, and 106 were up-regulated under both conditions. Notably, the set of genes over-expressed under either condition was preferentially encoded on plasmids, although this observation was more significant for the heat shock response. In contrast, during either saline shock or heat shock, the down-regulated genes were principally chromosomally encoded. Our functional analysis shows that genes encoding chaperone proteins were up-regulated during the heat shock response, whereas genes involved in the metabolism of compatible solutes were up-regulated following saline shock. Furthermore, we identified thirteen and nine ncRNAs that were differentially expressed under heat and saline shock, respectively, as well as eleven ncRNAs that had not been previously identified. Finally, using an in silico analysis, we studied the promoter motifs in all of the non-coding regions associated with the genes and ncRNAs up-regulated under both conditions. CONCLUSIONS: Our data suggest that the replicon contribution is different for different stress responses and that the heat shock response is more complex than the saline shock response. In general, this work exemplifies how strategies that not only consider differentially regulated genes but also regulatory elements of the stress response provide a more comprehensive view of bacterial gene regulation.

Analysis of Spo0M function in Bacillus subtilis.

Spo0M has been previously reported as a regulator of sporulation in Bacillus subtilis; however, little is known about the mechanisms through which it participates in sporulation, and there is no information to date that relates this protein to other processes in the bacterium. In this work we present evidence from proteomic, protein-protein interaction, morphological, subcellular localization microscopy and bioinformatics studies which indicate that Spo0M function is not necessarily restricted to sporulation, and point towards its involvement in other stages of the vegetative life cycle. In the current study, we provide evidence that Spo0M interacts with cytoskeletal proteins involved in cell division, which suggest a function additional to that previously described in sporulation. Spo0M expression is not restricted to the transition phase or sporulation; rather, its expression begins during the early stages of growth and Spo0M localization in B. subtilis depends on the bacterial life cycle and could be related to an additional proposed function. This is supported by our discovery of homologs in a broad distribution of bacterial genera, even in non-sporulating species. Our work paves the way for re-evaluation of the role of Spo0M in bacterial cell.

The cation diffusion facilitator protein EmfA of Rhizobium etli belongs to a novel subfamily of Mn(2+)/Fe(2+) transporters conserved in α-proteobacteria.

Manganese (Mn(2+)) plays a key role in important cellular functions such as oxidative stress response and bacterial virulence. The mechanisms of Mn(2+) homeostasis are not fully understood, there are few data regarding the functional and taxonomic diversity of Mn(2+) exporters. Our recent phylogeny of the cation diffusion facilitator (CDF) family of transporters classified the bacterial Mn(2+)-CDF transporters characterized to date, Streptococcus pneumoniae MntE and Deinococcus radiodurans DR1236, into two monophyletic groups. DR1236 was shown to belong to the highly-diverse metal specificity clade VI, together with TtCzrB, a Zn(2+)/Cd(2+) transporter from Thermus thermophilus, the Fe(2+) transporter Sll1263 from Synechocystis sp and eight uncharacterized homologs whose potential Mn(2+)/Zn(2+)/Cd(2+)/Fe(2+) specificities could not be accurately inferred because only eleven proteins were grouped in this clade. A new phylogeny inferred from the alignment of 197 clade VI homologs revealed three novel subfamilies of uncharacterized proteins. Remarkably, one of them contained 91 uncharacterized α-proteobacteria transporters (46% of the protein data set) grouped into a single subfamily. The Mn(2+)/Fe(2+) specificity of this subfamily was proposed through the functional characterization of the Rhizobium etli RHE_CH03072 gene. This gene was upregulated by Mn(2+), Zn(2+), Cd(2+) and Fe(2+) but conferred only Mn(2+) resistance to R. etli. The expression of the RHE_CH03072 gene in an E. coli mntP/zitB/zntA mutant did not relieve either Zn(2+) or Mn(2+) stress but slightly increased its Fe(2+) resistance. These results indicate that the RHE_CH03072 gene, now designated as emfA, encodes for a bacterial Mn(2+)/Fe(2+) resistance CDF protein, having orthologs in more than 60 α-proteobacterial species.

Phylogenomic analysis of Cation Diffusion Facilitator proteins uncovers Ni2+/Co2+ transporters.

The ubiquitous Cation Diffusion Facilitator proteins (CDF) play a key role in maintaining the cellular homeostasis of essential metal ions. Previous neighbor-joining phylogenetic analysis classified CDF proteins into three substrate-defined groups: Zn(2+), Fe(2+)/Zn(2+) and Mn(2+). These studies were unable to discern substrate-defined clades for Ni(2+), Co(2+), Cd(2+) and Cu(2+) transporters, despite their existence in this family. In this study we improved the accuracy of this previous functional classification using a phylogenomic approach based on a thorough maximum-likelihood phylogeny and the inclusion of recently characterized CDF transporters. The inference of CDF protein function predicted novel clades for Zn(2+), Fe(2+), Cd(2+) and Mn(2+). The Ni(2+)/Co(2+) and Co(2+) substrate specificities of two clades containing uncharacterized proteins were defined through the functional characterization of nepA and cepA metal inducible genes which independently conferred Ni(2+) and Co(2+) resistances to Rhizobium etli CFN42 and increased, respectively, Ni(2+)/Co(2+) and Co(2+) resistances to Escherichia coli. Neither NepA nor CepA confer Zn(2+), Fe(2+) and Mn(2+) resistances. The ability of NepA to confer Ni(2+)/Co(2+) resistance is dependent on clade-specific residues Asn(88) and Arg(197) whose mutations produce a non-functional protein.