Identification of a bacterial regulatory system for ribonucleotide reductases by phylogenetic profiling

The global population stru cture of Lacticaseibacillus rhamnosus and its application to an investigation of a rare case of infective endocarditis

BACKGROUND

Lacticaseibacillus (formerly Lactobacillus) rhamnosus is widely used in probiotics or food supplements to promote microbiome health and may also be part of the normal microbiota of the human gastrointestinal tract. However, it rarely also causes invasive or severe infections in patients. It has been postulated that these infections may originate from probiotics or from endogenous commensal reservoirs. In this report, we examine the population structure of Lacticaseibacillus rhamnosus and investigate the utility of using bacterial genomics to identify the source of invasive Lacticaseibacillus infections.

METHODS

Core genome phylogenetic analysis was performed on 602 L. rhamnosus genome sequences from the National Center for Biotechnology public database. This information was then used along with newly generated sequences of L. rhamnosus isolates from yogurt to investigate a fatal case of L. rhamnosus endocarditis.

RESULTS

Phylogenetic analysis demonstrated substantial genetic overlap of L. rhamnosus isolates cultured from food, probiotics, infected patients, and colonized individuals. This was applied to a patient who had both consumed yogurt and developed L. rhamnosus endocarditis to attempt to identify the source of his infection. The sequence of the isolate from the patient's bloodstream differed at only one nucleotide position from one of the yogurt isolates. Both isolates belonged to a clade, identified here as clade YC, composed of mostly gastrointestinal isolates from healthy individuals, some of which also differed by only a single nucleotide change from the patient's isolate.

CONCLUSIONS

As illustrated by this case, whole genome sequencing may be insufficient to reliably determine the source of invasive infections caused by L. rhamnosus.

Recombinant AcnB, NrdR and RibD of Acinetobacter baumannii and their potential interaction with DNA adenine methyltransferase AamA

In the last decades Acinetobacter baumannii developed into an increasingly challenging nosocomial pathogen. A. baumannii ATCC 17978 harbors a DNA-(adenine N6)-methyltransferase termed AamA. Previous studies revealed a low specific activity of AamA in vitro despite proven folding, which led us to speculate about possible interaction partners assisting AamA in targeting methylation sites. Here, applying a pulldown assay with subsequent mass spectrometry we identified aconitate hydratase 2 (AcnB) as possible interaction partner. In addition, we considered the putative transcriptional regulator gene nrdR (A1S_0220) and the pyrimidine deaminase/reductase gene ribD (A1S_0221) of A. baumannii strain ATCC 17978 to encode additional potential interaction partners due to their vicinity to the aamA gene (A1S_0222). Proteins were recombinantly produced in the milligram scale, purified to near homogeneity, and interactions with AamA were studied applying blue native gel electrophoreses, electrophoretic mobility shift assay, chemical cross-linking and co-immunoprecipitation. These analyses did not provide evidence of interaction between AamA and purified proteins. Solution structures of RibD, NrdR and AcnB were studied by small-angle X-ray scattering (SAXS) alone and in combination with AamA. While in the case of RibD and AcnB no evidence of an interaction with AamA was produced, addition of AamA to NrdR resulted in dissociation of long and rod-shaped polymeric NrdR structures, implying a specific but transient interaction. Moreover, we identified a molecular crowding effect possibly impeding the DNA methyltransferase activity in vivo and a sequence-independent DNA binding activity of AamA calling for continued efforts to identify the interaction network of AamA.

A nucleotide-sensing oligomerization mechanism that controls NrdR-dependent transcription of ribonucleotide reductases

Ribonucleotide reductase (RNR) is an essential enzyme that catalyzes the synthesis of DNA building blocks in virtually all living cells. NrdR, an RNR-specific repressor, controls the transcription of RNR genes and, often, its own, in most bacteria and some archaea. NrdR senses the concentration of nucleotides through its ATP-cone, an evolutionarily mobile domain that also regulates the enzymatic activity of many RNRs, while a Zn-ribbon domain mediates binding to NrdR boxes upstream of and overlapping the transcription start site of RNR genes. Here, we combine biochemical and cryo-EM studies of NrdR from Streptomyces coelicolor to show, at atomic resolution, how NrdR binds to DNA. The suggested mechanism involves an initial dodecamer loaded with two ATP molecules that cannot bind to DNA. When dATP concentrations increase, an octamer forms that is loaded with one molecule each of dATP and ATP per monomer. A tetramer derived from this octamer then binds to DNA and represses transcription of RNR. In many bacteria - including well-known pathogens such as Mycobacterium tuberculosis - NrdR simultaneously controls multiple RNRs and hence DNA synthesis, making it an excellent target for novel antibiotics development.

Potential Probiotic Pediococcus pentosaceus M41 Modulates Its Proteome Differentially for Tolerances Against Heat, Cold, Acid, and Bile Stresses

Probiotics containing functional food confer health benefits in addition to their nutritional properties. In this study, we have evaluated the differential proteomic responses of a potential novel probiotic Pediococcus pentosaceus M41 under heat, cold, acid, and bile stress conditions. We identified stress response proteins that could provide tolerances against these stresses and could be used as probiotic markers for evaluating stress tolerance. Pediococcus pentosaceus M41 was exposed for 2 h to each condition: 50°C (heat stress), 4°C (cold stress), pH 3.0 (acid stress) and 0.05% bile (bile stress). Proteomic analysis was carried out using 2D-IEF SDS PAGE and LC-MS/MS. Out of 60 identified proteins, 14 upregulated and 6 downregulated proteins were common among all the stress conditions. These proteins were involved in different biological functions such as translation-related proteins, carbohydrate metabolism (phosphoenolpyruvate phosphotransferase), histidine biosynthesis (imidazole glycerol phosphate synthase) and cell wall synthesis (tyrosine-protein kinase CapB). Proteins such as polysaccharide deacetylase, lactate oxidase, transcription repressor NrdR, dihydroxyacetone kinase were upregulated under three out of the four stress conditions. The differential expression of these proteins might be responsible for tolerance and protection of P. pentosaceus M41 against different stress conditions.

A second riboswitch class for the enzyme cofactor NAD. RNA (NEW YORK, N.Y.) 2021;27:99-105. [PMID: 33087526 PMCID: PMC7749635 DOI: 10.1261/rna.077891.120]

A bacterial noncoding RNA motif almost exclusively associated with pnuC genes was uncovered using comparative sequence analysis. Some PnuC proteins are known to transport nicotinamide riboside (NR), which is a component of the ubiquitous and abundant enzyme cofactor nicotinamide adenine dinucleotide (NAD⁺). Thus, we speculated that the newly found "pnuC motif" RNAs might function as aptamers for a novel class of NAD⁺-sensing riboswitches. RNA constructs that encompass the conserved nucleotides and secondary structure features that define the motif indeed selectively bind NAD⁺, nicotinamide mononucleotide (NMN), and NR. Mutations that disrupt strictly conserved nucleotides of the aptamer also disrupt ligand binding. These bioinformatic and biochemical findings indicate that pnuC motif RNAs are likely members of a second riboswitch class that regulates gene expression in response to NAD⁺ binding.

Comparative genomics of transport proteins in seven Bacteroides species

The communities of beneficial bacteria that live in our intestines, the gut microbiome, are important for the development and function of the immune system. Bacteroides species make up a significant fraction of the human gut microbiome, and can be probiotic and pathogenic, depending upon various genetic and environmental factors. These can cause disease conditions such as intra-abdominal sepsis, appendicitis, bacteremia, endocarditis, pericarditis, skin infections, brain abscesses and meningitis. In this study, we identify the transport systems and predict their substrates within seven Bacteroides species, all shown to be probiotic; however, four of them (B. thetaiotaomicron, B. vulgatus, B. ovatus, B. fragilis) can be pathogenic (probiotic and pathogenic; PAP), while B. cellulosilyticus, B. salanitronis and B. dorei are believed to play only probiotic roles (only probiotic; OP). The transport system characteristics of the four PAP and three OP strains were identified and tabulated, and results were compared among the seven strains, and with E. coli and Salmonella strains. The Bacteroides strains studied contain similarities and differences in the numbers and types of transport proteins tabulated, but both OP and PAP strains contain similar outer membrane carbohydrate receptors, pore-forming toxins and protein secretion systems, the similarities were noteworthy, but these Bacteroides strains showed striking differences with probiotic and pathogenic enteric bacteria, particularly with respect to their high affinity outer membrane receptors and auxiliary proteins involved in complex carbohydrate utilization. The results reveal striking similarities between the PAP and OP species of Bacteroides, and suggest that OP species may possess currently unrecognized pathogenic potential.

Transcriptomic analysis of Bifidobacterium longum subsp. longum BBMN68 in response to oxidative shock

Bifidobacterium longum strain BBMN68 is sensitive to low concentrations of oxygen. A transcriptomic study was performed to identify candidate genes for B. longum BBMN68's response to oxygen treatment (3%, v/v). Expression of genes and pathways of B. longum BBMN68 involved in nucleotide metabolism, amino acid transport, protein turnover and chaperones increased, and that of carbohydrate metabolism, translation and biogenesis decreased to adapt to the oxidative stress. Notably, expression of two classes of ribonucleotide reductase (RNR), which are important for deoxyribonucleotide biosynthesis, was rapidly and persistently induced. First, the class Ib RNR NrdHIEF was immediately upregulated after 5 min oxygen exposure, followed by the class III RNR NrdDG, which was upregulated after 20 min of exposure. The upregulated expression of branched-chain amino acids and tetrahydrofolate biosynthesis-related genes occurred in bifidobacteria in response to oxidative stress. These change toward to compensate for DNA and protein damaged by reactive oxygen species (ROS). In addition, oxidative stress resulted in improved B. longum BBMN68 cell hydrophobicity and autoaggregation. These results provide a rich resource for our understanding of the response mechanisms to oxidative stress in bifidobacteria.

LC-MS/MS proteomic analysis of starved Bacillus subtilis cells overexpressing ribonucleotide reductase (nrdEF): implications in stress-associated mutagenesis

The non-appropriate conditions faced by nutritionally stressed bacteria propitiate error-prone repair events underlying stationary-phase- or stress-associated mutagenesis (SPM). The genetic and molecular mechanisms involved in SPM have been deeply studied but the biochemical aspects of this process have so far been less explored. Previous evidence showed that under conditions of nutritional stress, non-dividing cells of strain B. subtilis YB955 overexpressing ribonucleotide reductase (RNR) exhibited a strong propensity to generate true reversions in the hisC952 (amber), metB5 (ochre) and leuC425 (missense) mutant alleles. To further advance our knowledge on the metabolic conditions underlying this hypermutagenic phenotype, a high-throughput LC-MS/MS proteomic analysis was performed in non-dividing cells of an amino acid-starved strain, deficient for NrdR, the RNR repressor. Compared with the parental strain, the level of 57 proteins was found to increase and of 80 decreases in the NrdR-deficient strain. The proteomic analysis revealed an altered content in proteins associated with the stringent response, nucleotide metabolism, DNA repair, and cell signaling in amino acid-starved cells of the ∆nrdR strain. Overall, our results revealed that amino acid-starved cells of strain B. subtilis ∆nrdR that escape from growth-limiting conditions exhibit a complex proteomic pattern reminiscent of a disturbed metabolism. Future experiments aimed to understand the consequences of disrupting the cell signaling pathways unveiled in this study, will advance our knowledge on the genetic adaptations deployed by bacteria to escape from growth-limiting environments.

Role of Ribonucleotide Reductase in Bacillus subtilis Stress-Associated Mutagenesis

The Gram-positive microorganism Bacillus subtilis relies on a single class Ib ribonucleotide reductase (RNR) to generate 2'-deoxyribonucleotides (dNDPs) for DNA replication and repair. In this work, we investigated the influence of RNR levels on B. subtilis stationary-phase-associated mutagenesis (SPM). Since RNR is essential in this bacterium, we engineered a conditional mutant of strain B. subtilis YB955 (hisC952 metB5 leu427) in which expression of the nrdEF operon was modulated by isopropyl-β-d-thiogalactopyranoside (IPTG). Moreover, genetic inactivation of ytcG, predicted to encode a repressor (NrdR) of nrdEF in this strain, dramatically increased the expression levels of a transcriptional nrdE-lacZ fusion. The frequencies of mutations conferring amino acid prototrophy in three genes were measured in cultures under conditions that repressed or induced RNR-encoding genes. The results revealed that RNR was necessary for SPM and overexpression of nrdEF promoted growth-dependent mutagenesis and SPM. We also found that nrdEF expression was induced by H₂O₂ and such induction was dependent on the master regulator PerR. These observations strongly suggest that the metabolic conditions operating in starved B. subtilis cells increase the levels of RNR, which have a direct impact on SPM.

IMPORTANCE

Results presented in this study support the concept that the adverse metabolic conditions prevailing in nutritionally stressed bacteria activate an oxidative stress response that disturbs ribonucleotide reductase (RNR) levels. Such an alteration of RNR levels promotes mutagenic events that allow Bacillus subtilis to escape from growth-limited conditions.

NrdR Transcription Regulation: Global Proteome Analysis and Its Role in Escherichia coli Viability and Virulence

Bacterial ribonucleotide reductases (RNRs) play an important role in the synthesis of dNTPs and their expression is regulated by the transcription factors, NrdR and Fur. Recent transcriptomic studies using deletion mutants have indicated a role for NrdR in bacterial chemotaxis and in the maintenance of topoisomerase levels. However, NrdR deletion alone has no effect on bacterial growth or virulence in infected flies or in human blood cells. Furthermore, transcriptomic studies are limited to the deletion strain alone, and so are inadequate for drawing biological implications when the NrdR repressor is active or abundant. Therefore, further examination is warranted of changes in the cellular proteome in response to both NrdR overexpression, as well as deletion, to better understand its functional relevance as a bacterial transcription repressor. Here, we profile bacterial fate under conditions of overexpression and deletion of NrdR in E. coli. Biochemical assays show auxiliary zinc enhances the DNA binding activity of NrdR. We also demonstrate at the physiological level that increased nrdR expression causes a significant reduction in bacterial growth and fitness even at normal temperatures, and causes lethality at elevated temperatures. Corroborating these direct effects, global proteome analysis following NrdR overexpression showed a significant decrease in global protein expression. In parallel, studies on complementary expression of downregulated essential genes polA, eno and thiL showed partial rescue of the fitness defect caused by NrdR overexpression. Deletion of downregulated non-essential genes ygfK and trxA upon NrdR overexpression resulted in diminished bacterial growth and fitness suggesting an additional role for NrdR in regulating other genes. Moreover, in comparison with NrdR deletion, E. coli cells overexpressing NrdR showed significantly diminished adherence to human epithelial cells, reflecting decreased bacterial virulence. These results suggest that elevated expression of NrdR could be a suitable means to retard bacterial growth and virulence, as its elevated expression reduces bacterial fitness and impairs host cell adhesion.

Structure, function, evolution, and application of bacterial Pnu-type vitamin transporters

Many bacteria can take up vitamins from the environment via specific transport machineries. Uptake is essential for organisms that lack complete vitamin biosynthesis pathways, but even in the presence of biosynthesis routes uptake is likely preferred, because it is energetically less costly. Pnu transporters represent a class of membrane transporters for a diverse set of B-type vitamins. They were identified 30 years ago and catalyze transport by the mechanism of facilitated diffusion, without direct coupling to ATP hydrolysis or transport of coupling ions. Instead, directionality is achieved by metabolic trapping, in which the vitamin substrate is converted into a derivative that cannot be transported, for instance by phosphorylation. The recent crystal structure of the nicotinamide riboside transporter PnuC has provided the first insights in substrate recognition and selectivity. Here, we will summarize the current knowledge about the function, structure, and evolution of Pnu transporters. Additionally, we will highlight their role for potential biotechnological and pharmaceutical applications.

The Novel Gene CRNDE Encodes a Nuclear Peptide (CRNDEP) Which Is Overexpressed in Highly Proliferating Tissues

CRNDE, recently described as the lncRNA-coding gene, is overexpressed at RNA level in human malignancies. Its role in gametogenesis, cellular differentiation and pluripotency has been suggested as well. Herein, we aimed to verify our hypothesis that the CRNDE gene may encode a protein product, CRNDEP. By using bioinformatics methods, we identified the 84-amino acid ORF encoded by one of two CRNDE transcripts, previously described by our research team. This ORF was cloned into two expression vectors, subsequently utilized in localization studies in HeLa cells. We also developed a polyclonal antibody against CRNDEP. Its specificity was confirmed in immunohistochemical, cellular localization, Western blot and immunoprecipitation experiments, as well as by showing a statistically significant decrease of endogenous CRNDEP expression in the cells with transient shRNA-mediated knockdown of CRNDE. Endogenous CRNDEP localizes predominantly to the nucleus and its expression seems to be elevated in highly proliferating tissues, like the parabasal layer of the squamous epithelium, intestinal crypts or spermatocytes. After its artificial overexpression in HeLa cells, in a fusion with either the EGFP or DsRed Monomer fluorescent tag, CRNDEP seems to stimulate the formation of stress granules and localize to them. Although the exact role of CRNDEP is unknown, our preliminary results suggest that it may be involved in the regulation of the cell proliferation. Possibly, CRNDEP also participates in oxygen metabolism, considering our in silico results, and the correlation between its enforced overexpression and the formation of stress granules. This is the first report showing the existence of a peptide encoded by the CRNDE gene.

Function of the Pseudomonas aeruginosa NrdR Transcription Factor: Global Transcriptomic Analysis and Its Role on Ribonucleotide Reductase Gene Expression

Ribonucleotide reductases (RNRs) are a family of sophisticated enzymes responsible for the synthesis of the deoxyribonucleotides (dNTPs), the building blocks for DNA synthesis and repair. Although any living cell must contain one RNR activity to continue living, bacteria have the capacity to encode different RNR classes in the same genome, allowing them to adapt to different environments and growing conditions. Pseudomonas aeruginosa is well known for its adaptability and surprisingly encodes all three known RNR classes (Ia, II and III). There must be a complex transcriptional regulation network behind this RNR activity, dictating which RNR class will be expressed according to specific growing conditions. In this work, we aim to uncover the role of the transcriptional regulator NrdR in P. aeruginosa. We demonstrate that NrdR regulates all three RNR classes, being involved in differential control depending on whether the growth conditions are aerobic or anaerobic. Moreover, we also identify for the first time that NrdR is not only involved in controlling RNR expression but also regulates topoisomerase I (topA) transcription. Finally, to obtain the entire picture of NrdR regulon, we performed a global transcriptomic analysis comparing the transcription profile of wild-type and nrdR mutant strains. The results provide many new data about the regulatory network that controls P. aeruginosa RNR transcription, bringing us a step closer to the understanding of this complex system.

An ancient riboswitch class in bacteria regulates purine biosynthesis and one-carbon metabolism

Over 30 years ago, ZTP (5-aminoimidazole-4-carboxamide riboside 5'-triphosphate), a modified purine biosynthetic intermediate, was proposed to signal 10-formyl-tetrahydrofolate (10f-THF) deficiency in bacteria. However, the mechanisms by which this putative alarmone or its precursor ZMP (5-aminoimidazole-4-carboxamide ribonucleotide, also known as AICAR) brings about any metabolic changes remain unexplained. Herein, we report the existence of a widespread riboswitch class that is most commonly associated with genes related to de novo purine biosynthesis and one-carbon metabolism. Biochemical data confirm that members of this riboswitch class selectively bind ZMP and ZTP with nanomolar affinity while strongly rejecting numerous natural analogs. Indeed, increases in the ZMP/ZTP pool, caused by folate stress in bacterial cells, trigger changes in the expression of a reporter gene fused to representative ZTP riboswitches in vivo. The wide distribution of this riboswitch class suggests that ZMP/ZTP signaling is important for species in numerous bacterial lineages.

Phylogenetic analysis of vitamin B12-related metabolism in Mycobacterium tuberculosis

Comparison of genome sequences from clinical isolates of Mycobacterium tuberculosis with phylogenetically-related pathogens Mycobacterium marinum, Mycobacterium kansasii, and Mycobacterium leprae reveals diversity amongst genes associated with vitamin B12-related metabolism. Diversity is generated by gene deletion events, differential acquisition of genes by horizontal transfer, and single nucleotide polymorphisms (SNPs) with predicted impact on protein function and transcriptional regulation. Differences in the B12 synthesis pathway, methionine biosynthesis, fatty acid catabolism, and DNA repair and replication are consistent with adaptations to different environmental niches and pathogenic lifestyles. While there is no evidence of further gene acquisition during expansion of the M. tuberculosis complex, the emergence of other forms of genetic diversity provides insights into continuing host-pathogen co-evolution and has the potential to identify novel targets for disease intervention.

Phylogenetic profiling: how much input data is enough?

Phylogenetic profiling is a well-established approach for predicting gene function based on patterns of gene presence and absence across species. Much of the recent developments have focused on methodological improvements, but relatively little is known about the effect of input data size on the quality of predictions. In this work, we ask: how many genomes and functional annotations need to be considered for phylogenetic profiling to be effective? Phylogenetic profiling generally benefits from an increased amount of input data. However, by decomposing this improvement in predictive accuracy in terms of the contribution of additional genomes and of additional annotations, we observed diminishing returns in adding more than ∼100 genomes, whereas increasing the number of annotations remained strongly beneficial throughout. We also observed that maximising phylogenetic diversity within a clade of interest improves predictive accuracy, but the effect is small compared to changes in the number of genomes under comparison. Finally, we show that these findings are supported in light of the Open World Assumption, which posits that functional annotation databases are inherently incomplete. All the tools and data used in this work are available for reuse from http://lab.dessimoz.org/14_phylprof. Scripts used to analyse the data are available on request from the authors.

Analysis of the roles of NrdR and DnaB from Streptococcus pyogenes in response to host defense

Toxic shock syndrome caused by Streptococcus pyogenes (S. pyogenes) is a re-emerging infectious disease. Many virulence-associated proteins play important roles in its pathogenesis and the production of these proteins is controlled by many regulatory factors. CovS is one of the most important two-component sensor proteins in S. pyogenes, and it has been analyzed extensively. Our recent analyses revealed the existence of a transposon between covS and nrdR in several strains, and we speculated that this insertion has some importance. Hence, we examined the significances of the NrdR stand-alone regulator and DnaB, which is encoded by the gene located immediately downstream of nrdR in S. pyogenes infection. We established an nrdR-only knockout strain, and both nrdR and partial dnaB knockout strain. These established knockout strains exhibited a deteriorated response to H2 O2 exposure. nrdR and partial dnaB knockout strain was more easily killed by human polynuclear blood cells, but the nrdR-only knockout strain had no significant difference compared to wild type in contrast to the combined knockout strain. In addition, the mouse infection model experiment illustrated that nrdR and partial dnaB knockout strain, but not the nrdR-only knockout strain, was less virulent compared with the parental strain. These results suggest that DnaB is involved in response to host defense.

Genome-scale co-expression network comparison across Escherichia coli and Salmonella enterica serovar Typhimurium reveals significant conservation at the regulon level of local regulators despite their dissimilar lifestyles

Availability of genome-wide gene expression datasets provides the opportunity to study gene expression across different organisms under a plethora of experimental conditions. In our previous work, we developed an algorithm called COMODO (COnserved MODules across Organisms) that identifies conserved expression modules between two species. In the present study, we expanded COMODO to detect the co-expression conservation across three organisms by adapting the statistics behind it. We applied COMODO to study expression conservation/divergence between Escherichia coli, Salmonella enterica, and Bacillus subtilis. We observed that some parts of the regulatory interaction networks were conserved between E. coli and S. enterica especially in the regulon of local regulators. However, such conservation was not observed between the regulatory interaction networks of B. subtilis and the two other species. We found co-expression conservation on a number of genes involved in quorum sensing, but almost no conservation for genes involved in pathogenicity across E. coli and S. enterica which could partially explain their different lifestyles. We concluded that despite their different lifestyles, no significant rewiring have occurred at the level of local regulons involved for instance, and notable conservation can be detected in signaling pathways and stress sensing in the phylogenetically close species S. enterica and E. coli. Moreover, conservation of local regulons seems to depend on the evolutionary time of divergence across species disappearing at larger distances as shown by the comparison with B. subtilis. Global regulons follow a different trend and show major rewiring even at the limited evolutionary distance that separates E. coli and S. enterica.

Comparative genomics and evolution of regulons of the LacI-family transcription factors

DNA-binding transcription factors (TFs) are essential components of transcriptional regulatory networks in bacteria. LacI-family TFs (LacI-TFs) are broadly distributed among certain lineages of bacteria. The majority of characterized LacI-TFs sense sugar effectors and regulate carbohydrate utilization genes. The comparative genomics approaches enable in silico identification of TF-binding sites and regulon reconstruction. To study the function and evolution of LacI-TFs, we performed genomics-based reconstruction and comparative analysis of their regulons. For over 1300 LacI-TFs from over 270 bacterial genomes, we predicted their cognate DNA-binding motifs and identified target genes. Using the genome context and metabolic subsystem analyses of reconstructed regulons, we tentatively assigned functional roles and predicted candidate effectors for 78 and 67% of the analyzed LacI-TFs, respectively. Nearly 90% of the studied LacI-TFs are local regulators of sugar utilization pathways, whereas the remaining 125 global regulators control large and diverse sets of metabolic genes. The global LacI-TFs include the previously known regulators CcpA in Firmicutes, FruR in Enterobacteria, and PurR in Gammaproteobacteria, as well as the three novel regulators—GluR, GapR, and PckR—that are predicted to control the central carbohydrate metabolism in three lineages of Alphaproteobacteria. Phylogenetic analysis of regulators combined with the reconstructed regulons provides a model of evolutionary diversification of the LacI protein family. The obtained genomic collection of in silico reconstructed LacI-TF regulons in bacteria is available in the RegPrecise database (http://regprecise.lbl.gov). It provides a framework for future structural and functional classification of the LacI protein family and identification of molecular determinants of the DNA and ligand specificity. The inferred regulons can be also used for functional gene annotation and reconstruction of sugar catabolic networks in diverse bacterial lineages.

Ribonucleotide reductases: essential enzymes for bacterial life

Ribonucleotide reductase (RNR) is a key enzyme that mediates the synthesis of deoxyribonucleotides, the DNA precursors, for DNA synthesis in every living cell. This enzyme converts ribonucleotides to deoxyribonucleotides, the building blocks for DNA replication, and repair. Clearly, RNR enzymes have contributed to the appearance of genetic material that exists today, being essential for the evolution of all organisms on Earth. The strict control of RNR activity and dNTP pool sizes is important, as pool imbalances increase mutation rates, replication anomalies, and genome instability. Thus, RNR activity should be finely regulated allosterically and at the transcriptional level. In this review we examine the distribution, the evolution, and the genetic regulation of bacterial RNRs. Moreover, this enzyme can be considered an ideal target for anti-proliferative compounds designed to inhibit cell replication in eukaryotic cells (cancer cells), parasites, viruses, and bacteria.

RegPrecise 3.0--a resource for genome-scale exploration of transcriptional regulation in bacteria

Background

Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in prokaryotes is one of the critical tasks of modern genomics. Bacteria from different taxonomic groups, whose lifestyles and natural environments are substantially different, possess highly diverged transcriptional regulatory networks. The comparative genomics approaches are useful for in silico reconstruction of bacterial regulons and networks operated by both transcription factors (TFs) and RNA regulatory elements (riboswitches).

Description

RegPrecise (http://regprecise.lbl.gov) is a web resource for collection, visualization and analysis of transcriptional regulons reconstructed by comparative genomics. We significantly expanded a reference collection of manually curated regulons we introduced earlier. RegPrecise 3.0 provides access to inferred regulatory interactions organized by phylogenetic, structural and functional properties. Taxonomy-specific collections include 781 TF regulogs inferred in more than 160 genomes representing 14 taxonomic groups of Bacteria. TF-specific collections include regulogs for a selected subset of 40 TFs reconstructed across more than 30 taxonomic lineages. Novel collections of regulons operated by RNA regulatory elements (riboswitches) include near 400 regulogs inferred in 24 bacterial lineages. RegPrecise 3.0 provides four classifications of the reference regulons implemented as controlled vocabularies: 55 TF protein families; 43 RNA motif families; ~150 biological processes or metabolic pathways; and ~200 effectors or environmental signals. Genome-wide visualization of regulatory networks and metabolic pathways covered by the reference regulons are available for all studied genomes. A separate section of RegPrecise 3.0 contains draft regulatory networks in 640 genomes obtained by an conservative propagation of the reference regulons to closely related genomes.

Conclusions

RegPrecise 3.0 gives access to the transcriptional regulons reconstructed in bacterial genomes. Analytical capabilities include exploration of: regulon content, structure and function; TF binding site motifs; conservation and variations in genome-wide regulatory networks across all taxonomic groups of Bacteria. RegPrecise 3.0 was selected as a core resource on transcriptional regulation of the Department of Energy Systems Biology Knowledgebase, an emerging software and data environment designed to enable researchers to collaboratively generate, test and share new hypotheses about gene and protein functions, perform large-scale analyses, and model interactions in microbes, plants, and their communities.

Cooperative and allosterically controlled nucleotide binding regulates the DNA binding activity of NrdR

Ribonucleotide reductases (RNRs) are required for the synthesis of deoxyribonucleoside triphosphates (dNTPs) from ribonucleotides. In Escherichia coli, regulation of RNR expression is co-ordinated with the cell cycle, and involves several regulatory proteins. One of these, NrdR, has recently been shown to regulate all three nrd operons that encode RNR isoenzymes. Repression by NrdR is believed to be stimulated by elevated dNTPs, although there is no direct evidence for this model. Here, we sought to elucidate the mechanism by which NrdR regulates nrd expression according to the abundance of (d)NTPs. We determined that ATP and dATP bind to NrdR in a negatively cooperative fashion, such that neither can fully occupy the protein. Both nucleotides also appear to act as positive heterotropic effectors, since the binding of one stimulates binding of the other. Nucleotide binding stimulates self-association of NrdR, with tri- and diphosphates stimulating oligomerization more effectively than monophosphates. As-prepared NrdR contains (deoxy)nucleoside monophosphates, diphosphates and triphosphates, and its DNA binding activity is inhibited by triphosphates and diphosphates but not by monophosphates. We propose a model in which NrdR selectively binds (deoxy)nucleoside triphosphates, which are hydrolysed to their monophosphate counterparts in order to regulate DNA binding.

MGcV: the microbial genomic context viewer for comparative genome analysis

Background

Conserved gene context is used in many types of comparative genome analyses. It is used to provide leads on gene function, to guide the discovery of regulatory sequences, but also to aid in the reconstruction of metabolic networks. We present the Microbial Genomic context Viewer (MGcV), an interactive, web-based application tailored to strengthen the practice of manual comparative genome context analysis for bacteria.

Results

MGcV is a versatile, easy-to-use tool that renders a visualization of the genomic context of any set of selected genes, genes within a phylogenetic tree, genomic segments, or regulatory elements. It is tailored to facilitate laborious tasks such as the interactive annotation of gene function, the discovery of regulatory elements, or the sequence-based reconstruction of gene regulatory networks. We illustrate that MGcV can be used in gene function annotation by visually integrating information on prokaryotic genes, like their annotation as available from NCBI with other annotation data such as Pfam domains, sub-cellular location predictions and gene-sequence characteristics such as GC content. We also illustrate the usefulness of the interactive features that allow the graphical selection of genes to facilitate data gathering (e.g. upstream regions, ID’s or annotation), in the analysis and reconstruction of transcription regulation. Moreover, putative regulatory elements and their corresponding scores or data from RNA-seq and microarray experiments can be uploaded, visualized and interpreted in (ranked-) comparative context maps. The ranked maps allow the interpretation of predicted regulatory elements and experimental data in light of each other.

Conclusion

MGcV advances the manual comparative analysis of genes and regulatory elements by providing fast and flexible integration of gene related data combined with straightforward data retrieval. MGcV is available at http://mgcv.cmbi.ru.nl.

ProPhylo: partial phylogenetic profiling to guide protein family construction and assignment of biological process

Background

Phylogenetic profiling is a technique of scoring co-occurrence between a protein family and some other trait, usually another protein family, across a set of taxonomic groups. In spite of several refinements in recent years, the technique still invites significant improvement. To be its most effective, a phylogenetic profiling algorithm must be able to examine co-occurrences among protein families whose boundaries are uncertain within large homologous protein superfamilies.

Results

Partial Phylogenetic Profiling (PPP) is an iterative algorithm that scores a given taxonomic profile against the taxonomic distribution of families for all proteins in a genome. The method works through optimizing the boundary of each protein family, rather than by relying on prebuilt protein families or fixed sequence similarity thresholds. Double Partial Phylogenetic Profiling (DPPP) is a related procedure that begins with a single sequence and searches for optimal granularities for its surrounding protein family in order to generate the best query profiles for PPP. We present ProPhylo, a high-performance software package for phylogenetic profiling studies through creating individually optimized protein family boundaries. ProPhylo provides precomputed databases for immediate use and tools for manipulating the taxonomic profiles used as queries.

Conclusion

ProPhylo results show universal markers of methanogenesis, a new DNA phosphorothioation-dependent restriction enzyme, and efficacy in guiding protein family construction. The software and the associated databases are freely available under the open source Perl Artistic License from ftp://ftp.jcvi.org/pub/data/ppp/.

Class I ribonucleotide reductases: metallocofactor assembly and repair in vitro and in vivo

Incorporation of metallocofactors essential for the activity of many enyzmes is a major mechanism of posttranslational modification. The cellular machinery required for these processes in the case of mono- and dinuclear nonheme iron and manganese cofactors has remained largely elusive. In addition, many metallocofactors can be converted to inactive forms, and pathways for their repair have recently come to light. The class I ribonucleotide reductases (RNRs) catalyze the conversion of nucleotides to deoxynucleotides and require dinuclear metal clusters for activity: an Fe(III)Fe(III)-tyrosyl radical (Y•) cofactor (class Ia), a Mn(III)Mn(III)-Y• cofactor (class Ib), and a Mn(IV)Fe(III) cofactor (class Ic). The class Ia, Ib, and Ic RNRs are structurally homologous and contain almost identical metal coordination sites. Recent progress in our understanding of the mechanisms by which the cofactor of each of these RNRs is generated in vitro and in vivo and by which the damaged cofactors are repaired is providing insight into how nature prevents mismetallation and orchestrates active cluster formation in high yields.

CT406 encodes a chlamydial ortholog of NrdR, a repressor of ribonucleotide reductase

Chlamydia trachomatis is an obligate intracellular bacterium that is dependent on its host cell for nucleotides. Chlamydia imports ribonucleotide triphosphates (NTPs) but not deoxyribonucleotide triphosphates (dNTPs) and instead uses ribonucleotide reductase to convert imported ribonucleotides into deoxyribonucleotides for DNA synthesis. The genes encoding ribonucleotide reductase have been recently shown to be negatively controlled by a conserved regulator called NrdR. In this study, we provide direct evidence that Escherichia coli NrdR is a transcriptional repressor and that C. trachomatis CT406 encodes its chlamydial ortholog. We showed that CT406 binds specifically to two NrdR boxes upstream of the nrdAB operon in C. trachomatis. Using an in vitro transcription assay, we confirmed that these NrdR boxes function as an operator since they were necessary and sufficient for CT406-mediated repression. We validated our in vitro findings with reporter studies in E. coli showing that both E. coli NrdR and CT406 repressed transcription from the E. coli nrdH and C. trachomatis nrdAB promoters in vivo. This in vivo repression was reversed by hydroxyurea treatment. Since hydroxyurea inhibits ribonucleotide reductase and reduces intracellular deoxyribonucleotide levels, these results suggest that NrdR activity is modulated by a deoxyribonucleotide corepressor.

Comparative genomic reconstruction of transcriptional networks controlling central metabolism in the Shewanella genus

Background

Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in bacteria is one of the critical tasks of modern genomics. The Shewanella genus is comprised of metabolically versatile gamma-proteobacteria, whose lifestyles and natural environments are substantially different from Escherichia coli and other model bacterial species. The comparative genomics approaches and computational identification of regulatory sites are useful for the in silico reconstruction of transcriptional regulatory networks in bacteria.

Results

To explore conservation and variations in the Shewanella transcriptional networks we analyzed the repertoire of transcription factors and performed genomics-based reconstruction and comparative analysis of regulons in 16 Shewanella genomes. The inferred regulatory network includes 82 transcription factors and their DNA binding sites, 8 riboswitches and 6 translational attenuators. Forty five regulons were newly inferred from the genome context analysis, whereas others were propagated from previously characterized regulons in the Enterobacteria and Pseudomonas spp.. Multiple variations in regulatory strategies between the Shewanella spp. and E. coli include regulon contraction and expansion (as in the case of PdhR, HexR, FadR), numerous cases of recruiting non-orthologous regulators to control equivalent pathways (e.g. PsrA for fatty acid degradation) and, conversely, orthologous regulators to control distinct pathways (e.g. TyrR, ArgR, Crp).

Conclusions

We tentatively defined the first reference collection of ~100 transcriptional regulons in 16 Shewanella genomes. The resulting regulatory network contains ~600 regulated genes per genome that are mostly involved in metabolism of carbohydrates, amino acids, fatty acids, vitamins, metals, and stress responses. Several reconstructed regulons including NagR for N-acetylglucosamine catabolism were experimentally validated in S. oneidensis MR-1. Analysis of correlations in gene expression patterns helps to interpret the reconstructed regulatory network. The inferred regulatory interactions will provide an additional regulatory constrains for an integrated model of metabolism and regulation in S. oneidensis MR-1.

Ribonucleotide reductases of Salmonella typhimurium: transcriptional regulation and differential role in pathogenesis

Ribonucleotide reductases (RNRs) are essential enzymes that carry out the de novo synthesis of deoxyribonucleotides by reducing ribonucleotides. There are three different classes of RNRs (I, II and III), all having different oxygen dependency and biochemical characteristics. Salmonella enterica serovar Typhimurium (S. Typhimurium) harbors class Ia, class Ib and class III RNRs in its genome. We have studied the transcriptional regulation of these three RNR classes in S. Typhimurium as well as their differential function during infection of macrophage and epithelial cells. Deletion of both NrdR and Fur, two main transcriptional regulators, indicates that Fur specifically represses the class Ib enzyme and that NrdR acts as a global repressor of all three classes. A Fur recognition sequence within the nrdHIEF promoter has also been described and confirmed by electrophoretic mobility shift assays (EMSA). In order to elucidate the role of each RNR class during infection, S. Typhimurium single and double RNR mutants (as well as Fur and NrdR mutants) were used in infection assays with macrophage and epithelial cell lines. Our results indicate class Ia to be mainly responsible for deoxyribonucleotide production during invasion and proliferation inside macrophages and epithelial cells. Neither class Ib nor class III seem to be essential for growth under these conditions. However, class Ib is able to maintain certain growth in an nrdAB mutant during the first hours of macrophage infection. Our results suggest that, during the early stages of macrophage infection, class Ib may contribute to deoxyribonucleotide synthesis by means of both an NrdR and a Fur-dependent derepression of nrdHIEF due to hydrogen peroxide production and DNA damage associated with the oxidative burst, thus helping to overcome the host defenses.

RegPredict: an integrated system for regulon inference in prokaryotes by comparative genomics approach

RegPredict web server is designed to provide comparative genomics tools for reconstruction and analysis of microbial regulons using comparative genomics approach. The server allows the user to rapidly generate reference sets of regulons and regulatory motif profiles in a group of prokaryotic genomes. The new concept of a cluster of co-regulated orthologous operons allows the user to distribute the analysis of large regulons and to perform the comparative analysis of multiple clusters independently. Two major workflows currently implemented in RegPredict are: (i) regulon reconstruction for a known regulatory motif and (ii) ab initio inference of a novel regulon using several scenarios for the generation of starting gene sets. RegPredict provides a comprehensive collection of manually curated positional weight matrices of regulatory motifs. It is based on genomic sequences, ortholog and operon predictions from the MicrobesOnline. An interactive web interface of RegPredict integrates and presents diverse genomic and functional information about the candidate regulon members from several web resources. RegPredict is freely accessible at http://regpredict.lbl.gov.

DnaA-ATP acts as a molecular switch to control levels of ribonucleotide reductase expression in Escherichia coli

Ribonucleotide reductase (RNR) is the bottleneck enzyme in the synthesis of dNTPs required for DNA replication. In order to avoid the mutagenic effects of imbalances in dNTPs the amount and activity of RNR enzyme in the cell is tightly regulated. RNR expression from the nrdAB operon is thus coupled to coincide with the initiation of DNA replication. However, the mechanism for the co-ordination of gene transcription and DNA replication remains to be elucidated. The timing and synchrony of DNA replication initiation in Escherichia coli is controlled in part by the binding of the DnaA protein to the origin of replication. DnaA is also a transcription factor of the nrdAB operon and could thus be the link between these two processes. Here we show that RNA polymerase can form a stable transcription initiation complex at the nrdAB promoter by direct interaction with the far upstream sites required for the timing of expression as a function of DNA replication. In addition, we show that the binding of DnaA on the promoter can either activate or repress transcription as a function of its concentration and its nucleotide-bound state. However, transcription regulation by DnaA does not significantly affect the timing of expression of RNR from the nrdAB operon.

RNRdb, a curated database of the universal enzyme family ribonucleotide reductase, reveals a high level of misannotation in sequences deposited to Genbank

Background

Ribonucleotide reductases (RNRs) catalyse the only known de novo pathway for deoxyribonucleotide synthesis, and are therefore essential to DNA-based life. While ribonucleotide reduction has a single evolutionary origin, significant differences between RNRs nevertheless exist, notably in cofactor requirements, subunit composition and allosteric regulation. These differences result in distinct operational constraints (anaerobicity, iron/oxygen dependence and cobalamin dependence), and form the basis for the classification of RNRs into three classes.

Description

In RNRdb (Ribonucleotide Reductase database), we have collated and curated all known RNR protein sequences with the aim of providing a resource for exploration of RNR diversity and distribution. By comparing expert manual annotations with annotations stored in Genbank, we find that significant inaccuracies exist in larger databases. To our surprise, only 23% of protein sequences included in RNRdb are correctly annotated across the key attributes of class, role and function, with 17% being incorrectly annotated across all three categories. This illustrates the utility of specialist databases for applications where a high degree of annotation accuracy may be important. The database houses information on annotation, distribution and diversity of RNRs, and links to solved RNR structures, and can be searched through a BLAST interface. RNRdb is accessible through a public web interface at http://rnrdb.molbio.su.se.

Conclusion

RNRdb is a specialist database that provides a reliable annotation and classification resource for RNR proteins, as well as a tool to explore distribution patterns of RNR classes. The recent expansion in available genome sequence data have provided us with a picture of RNR distribution that is more complex than believed only a few years ago; our database indicates that RNRs of all three classes are found across all three cellular domains. Moreover, we find a number of organisms that encode all three classes.

RegPrecise: a database of curated genomic inferences of transcriptional regulatory interactions in prokaryotes

The RegPrecise database (http://regprecise.lbl.gov) was developed for capturing, visualization and analysis of predicted transcription factor regulons in prokaryotes that were reconstructed and manually curated by utilizing the comparative genomic approach. A significant number of high-quality inferences of transcriptional regulatory interactions have been already accumulated for diverse taxonomic groups of bacteria. The reconstructed regulons include transcription factors, their cognate DNA motifs and regulated genes/operons linked to the candidate transcription factor binding sites. The RegPrecise allows for browsing the regulon collections for: (i) conservation of DNA binding sites and regulated genes for a particular regulon across diverse taxonomic lineages; (ii) sets of regulons for a family of transcription factors; (iii) repertoire of regulons in a particular taxonomic group of species; (iv) regulons associated with a metabolic pathway or a biological process in various genomes. The initial release of the database includes ∼11 500 candidate binding sites for ∼400 orthologous groups of transcription factors from over 350 prokaryotic genomes. Majority of these data are represented by genome-wide regulon reconstructions in Shewanella and Streptococcus genera and a large-scale prediction of regulons for the LacI family of transcription factors. Another section in the database represents the results of accurate regulon propagation to the closely related genomes.

Genetic makeup of the Corynebacterium glutamicum LexA regulon deduced from comparative transcriptomics and in vitro DNA band shift assays

The lexA gene of Corynebacterium glutamicum ATCC 13032 was deleted to create the mutant strain C. glutamicum NJ2114, which has an elongated cell morphology and an increased doubling time. To characterize the SOS regulon in C. glutamicum, the transcriptomes of NJ2114 and a DNA-damage-induced wild-type strain were compared with that of a wild-type control using DNA microarray hybridization. The expression data were combined with bioinformatic pattern searches for LexA binding sites, leading to the detection of 46 potential SOS boxes located upstream of differentially expressed transcription units. Binding of a hexahistidyl-tagged LexA protein to 40 double-stranded oligonucleotides containing the potential SOS boxes was demonstrated in vitro by DNA band shift assays. It turned out that LexA binds not only to SOS boxes in the promoter-operator region of upregulated genes, but also to SOS boxes detected upstream of downregulated genes. These results demonstrated that LexA controls directly the expression of at least 48 SOS genes organized in 36 transcription units. The deduced genes encode a variety of physiological functions, many of them involved in DNA repair and survival after DNA damage, but nearly half of them have hitherto unknown functions. Alignment of the LexA binding sites allowed the corynebacterial SOS box consensus sequence TcGAA(a/c)AnnTGTtCGA to be deduced. Furthermore, the common intergenic region of lexA and the differentially expressed divS-nrdR operon, encoding a cell division suppressor and a regulator of deoxyribonucleotide biosynthesis, was characterized in detail. Promoter mapping revealed differences in divS-nrdR expression during SOS response and normal growth conditions. One of the four LexA binding sites detected in the intergenic region is involved in regulating divS-nrdR transcription, whereas the other sites are apparently used for negative autoregulation of lexA expression.

Regulation by transcription factors in bacteria: beyond description

Transcription is an essential step in gene expression and its understanding has been one of the major interests in molecular and cellular biology. By precisely tuning gene expression, transcriptional regulation determines the molecular machinery for developmental plasticity, homeostasis and adaptation. In this review, we transmit the main ideas or concepts behind regulation by transcription factors and give just enough examples to sustain these main ideas, thus avoiding a classical ennumeration of facts. We review recent concepts and developments: cis elements and trans regulatory factors, chromosome organization and structure, transcriptional regulatory networks (TRNs) and transcriptomics. We also summarize new important discoveries that will probably affect the direction of research in gene regulation: epigenetics and stochasticity in transcriptional regulation, synthetic circuits and plasticity and evolution of TRNs. Many of the new discoveries in gene regulation are not extensively tested with wetlab approaches. Consequently, we review this broad area in Inference of TRNs and Dynamical Models of TRNs. Finally, we have stepped backwards to trace the origins of these modern concepts, synthesizing their history in a timeline schema.

Functional analysis of the Streptomyces coelicolor NrdR ATP-cone domain: role in nucleotide binding, oligomerization, and DNA interactions

Ribonucleotide reductases (RNRs) are essential enzymes in all living cells, providing the only known de novo pathway for the biosynthesis of deoxyribonucleotides (dNTPs), the immediate precursors of DNA synthesis and repair. RNRs catalyze the controlled reduction of all four ribonucleotides to maintain a balanced pool of dNTPs during the cell cycle. Streptomyces species contain genes, nrdAB and nrdJ, coding for oxygen-dependent class I and oxygen-independent class II RNRs, either of which is sufficient for vegetative growth. Both sets of genes are transcriptionally repressed by NrdR. NrdR contains a zinc ribbon DNA-binding domain and an ATP-cone domain similar to that present in the allosteric activity site of many class I and class III RNRs. Purified NrdR contains up to 1 mol of tightly bound ATP or dATP per mol of protein and binds to tandem 16-bp sequences, termed NrdR-boxes, present in the upstream regulatory regions of bacterial RNR operons. Previously, we showed that the ATP-cone domain alone determines nucleotide binding and that an NrdR mutant defective in nucleotide binding was unable to bind to DNA probes containing NrdR-boxes. These observations led us to propose that when NrdR binds ATP/dATP it undergoes a conformational change that affects DNA binding and hence RNR gene expression. In this study, we analyzed a collection of ATP-cone mutant proteins containing changes in residues inferred to be implicated in nucleotide binding and show that they result in pleiotrophic effects on ATP/dATP binding, on protein oligomerization, and on DNA binding. A model is proposed to integrate these observations.

Function and regulation of class I ribonucleotide reductase-encoding genes in mycobacteria

Ribonucleotide reductases (RNRs) are crucial to all living cells, since they provide deoxyribonucleotides (dNTPs) for DNA synthesis and repair. In Mycobacterium tuberculosis, a class Ib RNR comprising nrdE- and nrdF2-encoded subunits is essential for growth in vitro. Interestingly, the genome of this obligate human pathogen also contains the nrdF1 (Rv1981c) and nrdB (Rv0233) genes, encoding an alternate class Ib RNR small (R2) subunit and a putative class Ic RNR R2 subunit, respectively. However, the role(s) of these subunits in dNTP provision during M. tuberculosis pathogenesis is unknown. In this study, we demonstrate that nrdF1 and nrdB are dispensable for the growth and survival of M. tuberculosis after exposure to various stresses in vitro and, further, that neither gene is required for growth and survival in mice. These observations argue against a specialist role for the alternate R2 subunits under the conditions tested. Through the construction of nrdR-deficient mutants of M. tuberculosis and Mycobacterium smegmatis, we establish that the genes encoding the essential class Ib RNR subunits are specifically regulated by an NrdR-type repressor. Moreover, a strain of M. smegmatis mc(2)155 lacking the 56-kb chromosomal region, which includes duplicates of nrdHIE and nrdF2, and a mutant retaining only one copy of nrdF2 are shown to be hypersensitive to the class I RNR inhibitor hydroxyurea as a result of depleted levels of the target. Together, our observations identify a potential vulnerability in dNTP provision in mycobacteria and thereby offer a compelling rationale for pursuing the class Ib RNR as a target for drug discovery.

Modes of overinitiation, dnaA gene expression, and inhibition of cell division in a novel cold-sensitive hda mutant of Escherichia coli

The chromosomal replication cycle is strictly coordinated with cell cycle progression in Escherichia coli. ATP-DnaA initiates replication, leading to loading of the DNA polymerase III holoenzyme. The DNA-loaded form of the beta clamp subunit of the polymerase binds the Hda protein, which promotes ATP-DnaA hydrolysis, yielding inactive ADP-DnaA. This regulation is required to repress overinitiation. In this study, we have isolated a novel cold-sensitive hda mutant, the hda-185 mutant. The hda-185 mutant caused overinitiation of chromosomal replication at 25 degrees C, which most likely led to blockage of replication fork progress. Consistently, the inhibition of colony formation at 25 degrees C was suppressed by disruption of the diaA gene, an initiation stimulator. Disruption of the seqA gene, an initiation inhibitor, showed synthetic lethality with hda-185 even at 42 degrees C. The cellular ATP-DnaA level was increased in an hda-185-dependent manner. The cellular concentrations of DnaA protein and dnaA mRNA were comparable at 25 degrees C to those in a wild-type hda strain. We also found that multiple copies of the ribonucleotide reductase genes (nrdAB or nrdEF) or dnaB gene repressed overinitiation. The cellular levels of dATP and dCTP were elevated in cells bearing multiple copies of nrdAB. The catalytic site within NrdA was required for multicopy suppression, suggesting the importance of an active form of NrdA or elevated levels of deoxyribonucleotides in inhibition of overinitiation in the hda-185 cells. Cell division in the hda-185 mutant was inhibited at 25 degrees C in a LexA regulon-independent manner, suggesting that overinitiation in the hda-185 mutant induced a unique division inhibition pathway.

Proteomic analysis of stationary phase in the marine bacterium "Candidatus Pelagibacter ubique"

"Candidatus Pelagibacter ubique," an abundant marine alphaproteobacterium, subsists in nature at low ambient nutrient concentrations and may often be exposed to nutrient limitation, but its genome reveals no evidence of global regulatory mechanisms for adaptation to stationary phase. High-resolution capillary liquid chromatography coupled online to an LTQ mass spectrometer was used to build an accurate mass and time (AMT) tag library that enabled quantitative examination of proteomic differences between exponential- and stationary-phase "Ca. Pelagibacter ubique" cells cultivated in a seawater medium. The AMT tag library represented 65% of the predicted protein-encoding genes. "Ca. Pelagibacter ubique" appears to respond adaptively to stationary phase by increasing the abundance of a suite of proteins that contribute to homeostasis rather than undergoing a major remodeling of its proteome. Stationary-phase abundances increased significantly for OsmC and thioredoxin reductase, which may mitigate oxidative damage in "Ca. Pelagibacter," as well as for molecular chaperones, enzymes involved in methionine and cysteine biosynthesis, proteins involved in rho-dependent transcription termination, and the signal transduction enzyme CheY-FisH. We speculate that this limited response may enable "Ca. Pelagibacter ubique" to cope with ambient conditions that deprive it of nutrients for short periods and, furthermore, that the ability to resume growth overrides the need for a more comprehensive global stationary-phase response to create a capacity for long-term survival.

Transcriptional regulation of NAD metabolism in bacteria: genomic reconstruction of NiaR (YrxA) regulon

A comparative genomic approach was used to reconstruct transcriptional regulation of NAD biosynthesis in bacteria containing orthologs of Bacillus subtilis gene yrxA, a previously identified niacin-responsive repressor of NAD de novo synthesis. Members of YrxA family (re-named here NiaR) are broadly conserved in the Bacillus/Clostridium group and in the deeply branching Fusobacteria and Thermotogales lineages. We analyzed upstream regions of genes associated with NAD biosynthesis to identify candidate NiaR-binding DNA motifs and assess the NiaR regulon content in these species. Representatives of the two distinct types of candidate NiaR-binding sites, characteristic of the Firmicutes and Thermotogales, were verified by an electrophoretic mobility shift assay. In addition to transcriptional control of the nadABC genes, the NiaR regulon in some species extends to niacin salvage (the pncAB genes) and includes uncharacterized membrane proteins possibly involved in niacin transport. The involvement in niacin uptake proposed for one of these proteins (re-named NiaP), encoded by the B. subtilis gene yceI, was experimentally verified. In addition to bacteria, members of the NiaP family are conserved in multicellular eukaryotes, including human, pointing to possible NaiP involvement in niacin utilization in these organisms. Overall, the analysis of the NiaR and NrtR regulons (described in the accompanying paper) revealed mechanisms of transcriptional regulation of NAD metabolism in nearly a hundred diverse bacteria.

Transcriptional regulation of NAD metabolism in bacteria: NrtR family of Nudix-related regulators

A novel family of transcription factors responsible for regulation of various aspects of NAD synthesis in a broad range of bacteria was identified by comparative genomics approach. Regulators of this family (here termed NrtR for Nudix-related transcriptional regulators), currently annotated as ADP-ribose pyrophosphatases from the Nudix family, are composed of an N-terminal Nudix-like effector domain and a C-terminal DNA-binding HTH-like domain. NrtR regulons were reconstructed in diverse bacterial genomes by identification and comparative analysis of NrtR-binding sites upstream of genes involved in NAD biosynthetic pathways. The candidate NrtR-binding DNA motifs showed significant variability between microbial lineages, although the common consensus sequence could be traced for most of them. Bioinformatics predictions were experimentally validated by gel mobility shift assays for two NrtR family representatives. ADP-ribose, the product of glycohydrolytic cleavage of NAD, was found to suppress the in vitro binding of NrtR proteins to their DNA target sites. In addition to a major role in the direct regulation of NAD homeostasis, some members of NrtR family appear to have been recruited for the regulation of other metabolic pathways, including sugar pentoses utilization and biogenesis of phosphoribosyl pyrophosphate. This work and the accompanying study of NiaR regulon demonstrate significant variability of regulatory strategies for control of NAD metabolic pathway in bacteria.

Practical and theoretical advances in predicting the function of a protein by its phylogenetic distribution

The gap between the amount of genome information released by genome sequencing projects and our knowledge about the proteins' functions is rapidly increasing. To fill this gap, various 'genomic-context' methods have been proposed that exploit sequenced genomes to predict the functions of the encoded proteins. One class of methods, phylogenetic profiling, predicts protein function by correlating the phylogenetic distribution of genes with that of other genes or phenotypic characteristics. The functions of a number of proteins, including ones of medical relevance, have thus been predicted and subsequently confirmed experimentally. Additionally, various approaches to measure the similarity of phylogenetic profiles and to account for the phylogenetic bias in the data have been proposed. We review the successful applications of phylogenetic profiling and analyse the performance of various profile similarity measures with a set of one microsporidial and 25 fungal genomes. In the fungi, phylogenetic profiling yields high-confidence predictions for the highest and only the highest scoring gene pairs illustrating both the power and the limitations of the approach. Both practical examples and theoretical considerations suggest that in order to get a reliable and specific picture of a protein's function, results from phylogenetic profiling have to be combined with other sources of evidence.

A model of evolution with constant selective pressure for regulatory DNA sites

BACKGROUND

Molecular evolution is usually described assuming a neutral or weakly non-neutral substitution model. Recently, new data have become available on evolution of sequence regions under a selective pressure, e.g. transcription factor binding sites. To reconstruct the evolutionary history of such sequences, one needs evolutionary models that take into account a substantial constant selective pressure.

RESULTS

We present a simple evolutionary model with a single preferred (consensus) nucleotide and the neutral substitution model adopted for all other nucleotides. This evolutionary model has a rate matrix in which all substitutions that do not involve the consensus nucleotide occur with the same rate. The model has two time scales for achieving a stationary distribution; in the general case only one of the two rate parameters can be evaluated from the stationary distribution. In the middle-time zone, a counterintuitive behavior was observed for some parameter values, with a probability of conservation for a non-consensus nucleotide greater than that for the consensus nucleotide. Such an effect can be observed only in the case of weak preference for the consensus nucleotide, when the probability to observe the consensus nucleotide in the stationary distribution is less than 1/2. If the substitution rate is represented as a product of mutation and fixation, only the fixation can be calculated from the stationary distribution. The exhibited conservation of non-consensus nucleotides does not take place if the elements of mutation matrix are identical, and can be related to the reduced mutation rate between the non-consensus nucleotides. This bias can have no effect on the stationary distribution of nucleotide frequencies calculated over the ensemble of multiple alignments, e.g. transcription factor binding sites upstream of different sets of co-regulated orthologous genes.

CONCLUSION

The derived model can be used as a null model when analyzing the evolution of orthologous transcription factor binding sites. In particular, our findings show that a nucleotide preferred at some position of a multiple alignment of binding sites for some transcription factor in the same genome is not necessarily the most conserved nucleotide in an alignment of orthologous sites from different species. However, this effect can take place only in the case of a mutation matrix whose elements are not identical.