Transcriptional regulation of transport and utilization systems for hexuronides, hexuronates and hexonates in gamma purple bacteria

Coordination of host and endosymbiont gene expression governs endosymbiont growth and elimination in the cereal weevil Sitophilus spp

BACKGROUND

Insects living in nutritionally poor environments often establish long-term relationships with intracellular bacteria that supplement their diets and improve their adaptive and invasive powers. Even though these symbiotic associations have been extensively studied on physiological, ecological, and evolutionary levels, few studies have focused on the molecular dialogue between host and endosymbionts to identify genes and pathways involved in endosymbiosis control and dynamics throughout host development.

RESULTS

We simultaneously analyzed host and endosymbiont gene expression during the life cycle of the cereal weevil Sitophilus oryzae, from larval stages to adults, with a particular emphasis on emerging adults where the endosymbiont Sodalis pierantonius experiences a contrasted growth-climax-elimination dynamics. We unraveled a constant arms race in which different biological functions are intertwined and coregulated across both partners. These include immunity, metabolism, metal control, apoptosis, and bacterial stress response.

CONCLUSIONS

The study of these tightly regulated functions, which are at the center of symbiotic regulations, provides evidence on how hosts and bacteria finely tune their gene expression and respond to different physiological challenges constrained by insect development in a nutritionally limited ecological niche. Video Abstract.

PcaR, a GntR/FadR Family Transcriptional Repressor Controls the Transcription of Phenazine-1-Carboxylic Acid 1,2-Dioxygenase Gene Cluster in Sphingomonas histidinilytica DS-9

In our previous study, the phenazine-1-carboxylic acid (PCA) 1,2-dioxygenase gene cluster (pcaA1A2A3A4 cluster) in Sphingomonas histidinilytica DS-9 was identified to be responsible for the conversion of PCA to 1,2-dihydroxyphenazine (Ren Y, Zhang M, Gao S, Zhu Q, et al. 2022. Appl Environ Microbiol 88:e00543-22). However, the regulatory mechanism of the pcaA1A2A3A4 cluster has not been elucidated yet. In this study, the pcaA1A2A3A4 cluster was found to be transcribed as two divergent operons: pcaA3-ORF5205 (named A3-5205 operon) and pcaA1A2-ORF5208-pcaA4-ORF5210 (named A1-5210 operon). The promoter regions of the two operons were overlapped. PcaR acts as a transcriptional repressor of the pcaA1A2A3A4 cluster, and it belongs to GntR/FadR family transcriptional regulator. Gene disruption of pcaR can shorten the lag phase of PCA degradation. The results of electrophoretic mobility shift assay and DNase I footprinting showed that PcaR binds to a 25-bp motif in the ORF5205-pcaA1 intergenic promoter region to regulate the expression of two operons. The 25-bp motif covers the -10 region of the promoter of A3-5205 operon and the -35 region and -10 region of the promoter of A1-5210 operon. The TNGT/ANCNA box within the motif was essential for PcaR binding to the two promoters. PCA acted as an effector of PcaR, preventing it from binding to the promoter region and repressing the transcription of the pcaA1A2A3A4 cluster. In addition, PcaR represses its own transcription, and this repression can be relieved by PCA. This study reveals the regulatory mechanism of PCA degradation in strain DS-9, and the identification of PcaR increases the variety of regulatory model of the GntR/FadR-type regulator. IMPORTANCE Sphingomonas histidinilytica DS-9 is a phenazine-1-carboxylic acid (PCA)-degrading strain. The 1,2-dioxygenase gene cluster (pcaA1A2A3A4 cluster, encoding dioxygenase PcaA1A2, reductase PcaA3, and ferredoxin PcaA4) is responsible for the initial degradation step of PCA and widely distributed in Sphingomonads, but its regulatory mechanism has not been investigated yet. In this study, a GntR/FadR-type transcriptional regulator PcaR repressing the transcription of pcaA1A2A3A4 cluster and pcaR gene was identified and characterized. The binding site of PcaR in ORF5205-pcaA1 intergenic promoter region contains a TNGT/ANCNA box, which is important for the binding. These findings enhance our understanding of the molecular mechanism of PCA degradation.

Differential Impact of Hexuronate Regulators ExuR and UxuR on the Escherichia coli Proteome

ExuR and UxuR are paralogous proteins belonging to the GntR family of transcriptional regulators. Both are known to control hexuronic acid metabolism in a variety of Gammaproteobacteria but the relative impact of each of them is still unclear. Here, we apply 2D difference electrophoresis followed by mass-spectrometry to characterise the changes in the Escherichia coli proteome in response to a uxuR or exuR deletion. Our data clearly show that the effects are different: deletion of uxuR resulted in strongly enhanced expression of D-mannonate dehydratase UxuA and flagellar protein FliC, and in a reduced amount of outer membrane porin OmpF, while the absence of ExuR did not significantly alter the spectrum of detected proteins. Consequently, the physiological roles of proteins predicted as homologs seem to be far from identical. Effects of uxuR deletion were largely dependent on the cultivation conditions: during growth with glucose, UxuA and FliC were dramatically altered, while during growth with glucuronate, activation of both was not so prominent. During the growth with glucose, maximal activation was detected for FliC. This was further confirmed by expression analysis and physiological tests, thus suggesting the involvement of UxuR in the regulation of bacterial motility and biofilm formation.

Use of Transposon Directed Insertion-Site Sequencing to Probe the Antibacterial Mechanism of a Model Honey on E. coli K-12

Antimicrobial resistance is an ever-growing health concern worldwide that has created renewed interest in the use of traditional anti-microbial treatments, including honey. However, understanding the underlying mechanism of the anti-microbial action of honey has been hampered due to the complexity of its composition. High throughput genetic tools could assist in understanding this mechanism. In this study, the anti-bacterial mechanism of a model honey, made of sugars, hydrogen peroxide, and gluconic acid, was investigated using genome-wide transposon mutagenesis combined with high-throughput sequencing (TraDIS), with the strain Escherichia coli K-12 MG1655 as the target organism. We identified a number of genes which when mutated caused a severe loss of fitness when cells were exposed to the model honey. These genes encode membrane proteins including those involved in uptake of essential molecules, and components of the electron transport chain. They are enriched for pathways involved in intracellular homeostasis and redox activity. Genes involved in assembly and activity of formate dehydrogenase O (FDH-O) were of particular note. The phenotypes of mutants in a subset of the genes identified were confirmed by phenotypic screening of deletion strains. We also found some genes which when mutated led to enhanced resistance to treatment with the model honey. This study identifies potential synergies between the main honey stressors and provides insights into the global antibacterial mechanism of this natural product.

Transcription factor allosteric regulation through substrate coordination to zinc

The development of new synthetic biology circuits for biotechnology and medicine requires deeper mechanistic insight into allosteric transcription factors (aTFs). Here we studied the aTF UxuR, a homodimer of two domains connected by a highly flexible linker region. To explore how ligand binding to UxuR affects protein dynamics we performed molecular dynamics simulations in the free protein, the aTF bound to the inducer D-fructuronate or the structural isomer D-glucuronate. We then validated our results by constructing a sensor plasmid for D-fructuronate in Escherichia coli and performed site-directed mutagenesis. Our results show that zinc coordination is necessary for UxuR function since mutation to alanines prevents expression de-repression by D-fructuronate. Analyzing the different complexes, we found that the disordered linker regions allow the N-terminal domains to display fast and large movements. When the inducer is bound, UxuR can sample an open conformation with a more pronounced negative charge at the surface of the N-terminal DNA binding domains. In opposition, in the free and D-glucuronate bond forms the protein samples closed conformations, with a more positive character at the surface of the DNA binding regions. These molecular insights provide a new basis to harness these systems for biological systems engineering.

Regulatory Role of PlaR (YiaJ) for Plant Utilization in Escherichia coli K-12

Outside a warm-blooded animal host, the enterobacterium Escherichia coli K-12 is also able to grow and survive in stressful nature. The major organic substance in nature is plant, but the genetic system of E. coli how to utilize plant-derived materials as nutrients is poorly understood. Here we describe the set of regulatory targets for uncharacterized IclR-family transcription factor YiaJ on the E. coli genome, using gSELEX screening system. Among a total of 18 high-affinity binding targets of YiaJ, the major regulatory target was identified to be the yiaLMNOPQRS operon for utilization of ascorbate from fruits and galacturonate from plant pectin. The targets of YiaJ also include the genes involved in the utilization for other plant-derived materials as nutrients such as fructose, sorbitol, glycerol and fructoselysine. Detailed in vitro and in vivo analyses suggest that L-ascorbate and α-D-galacturonate are the effector ligands for regulation of YiaJ function. These findings altogether indicate that YiaJ plays a major regulatory role in expression of a set of the genes for the utilization of plant-derived materials as nutrients for survival. PlaR was also suggested to play protecting roles of E. coli under stressful environments in nature, including the formation of biofilm. We then propose renaming YiaJ to PlaR (regulator of plant utilization). The natural hosts of enterobacterium Escherichia coli are warm-blooded animals, but even outside hosts, E. coli can survive even under stressful environments. On earth, the most common organic materials to be used as nutrients by E. coli are plant-derived components, but up to the present time, the genetic system of E. coli for plant utilization is poorly understand. In the course of gSELEX screening of the regulatory targets for hitherto uncharacterized TFs, we identified in this study the involvement of the IclR-family YiaJ in the regulation of about 20 genes or operons, of which the majority are related to the catabolism of plant-derived materials such as ascorbate, galacturonate, sorbitol, fructose and fructoselysine. Therefore, we propose to rename YiaJ to PlaR (regulator of plant utilization).

Diet-derived galacturonic acid regulates virulence and intestinal colonization in enterohaemorrhagic Escherichia coli and Citrobacter rodentium

Enteric pathogens sense the complex chemistry within the gastrointestinal (GI) tract to efficiently compete with the resident microbiota and establish a colonization niche. Here we show that enterohemorrhagic E. coli (EHEC), and its surrogate murine infection model Citrobacter rodentium, sense galacturonic-acid to initiate a multi-layered program towards successful mammalian infection. Galacturonic-acid utilization as a carbon source aids the initial pathogen expansion. The main source of galacturonic-acid is dietary pectin, which is broken into galacturonic-acid by the prominent member of the microbiota, Bacteroides thetaiotamicron (Bt). This regulation occurs through the ExuR transcription factor. However, galacturonic-acid is also sensed as a signal through ExuR to modulate the expression of the genes encoding a molecular syringe known as a type three secretion system (T3SS) leading to infectious colitis and inflammation. Galacturonic-acid moonlights as a nutrient and a signal directing the exquisite microbiota-pathogen relationships within the GI tract. Importantly, this work highlights that differential dietary sugar availability impacts the relationship between the microbiota and enteric pathogens, as well as disease outcomes.

Adaptation to host in Vibrio vulnificus, a zoonotic pathogen that causes septicemia in fish and humans

Vibrio vulnificus is a siderophilic pathogen spreading due to global warming. The zoonotic strains constitute a clonal-complex related to fish farms that are distributed worldwide. In this study, we applied a transcriptomic and single gene approach and discover that the zoonotic strains bypassed the iron requirement of the species thanks to the acquisition of two iron-regulated outer membrane proteins (IROMPs) involved in resistance to fish innate immunity. Both proteins have been acquired by horizontal gene transfer and are contributing to the successful spreading of this clonal-complex. We have also discovered that the zoonotic strains express a virulent phenotype in the blood of its main susceptible hosts (iron-overloaded humans and healthy eels) by combining a host-specific protective envelope with the common expression of two toxins (VvhA and RtxA1), one of which (RtxA1) is directly involved in sepsis. Finally, we found that both IROMPs are also present in other fish pathogenic species and have recently been transmitted to the phylogenetic lineage involved in human primary sepsis after raw seafood ingestion. Together our results highlight the potential hazard that the aquaculture industry poses to public health, which is of particular relevance in the context of a warming world.

Overproduction and purification of the Escherichia coli transcription factors "toxic" to a bacterial cell

Transcription factors play a crucial role in control of life of a bacterial cell, working as switchers to a different life style or pathogenicity. To reconstruct the network of regulatory events taking place in changing growth conditions, we need to know regulons of as many transcription factors as possible, and motifs recognized by them. Experimentally this can be attained via ChIP-seq in vivo, SELEX and DNAse I footprinting in vitro. All these approaches require large amounts of purified proteins. However, overproduction of transcription factors leading to their extensive binding to the regulatory elements on the DNA make them toxic to a bacterial cell thus significantly complicating production of a soluble protein. Here, on the example of three regulators from Escherichia coli, UxuR, ExuR, and LeuO, we show that stable production of toxic transcription factors in a soluble fraction can be significantly enhanced by holding the expression of a recombinant protein back at the early stages of bacterial growth. This can be achieved by cloning genes together with their regulatory regions containing repressor sites, with subsequent growth in a very rich media where activity of excessive regulators is not crucial, followed by induction with a very low concentration of an inducer. Schemes of further purification of these proteins were developed, and functional activity was confirmed.

Molecular and Functional Insights into the Regulation of d-Galactonate Metabolism by the Transcriptional Regulator DgoR in Escherichia coli

d-Galactonate, an aldonic sugar acid, is used as a carbon source by Escherichia coli, and the structural dgo genes involved in its metabolism have previously been investigated. Here, using genetic, biochemical and bioinformatics approaches, we present the first detailed molecular and functional insights into the regulation of d-galactonate metabolism in E. coli K-12 by the transcriptional regulator DgoR. We found that dgoR deletion accelerates the growth of E. coli in d-galactonate concomitant with the strong constitutive expression of dgo genes. In the dgo locus, sequence upstream of dgoR alone harbors the d-galactonate-inducible promoter that likely drives the expression of all dgo genes. DgoR exerts repression on the dgo operon by binding two inverted repeats overlapping the dgo promoter. Binding of d-galactonate induces a conformational change in DgoR to derepress the dgo operon. The findings from our work firmly place DgoR in the GntR family of transcriptional regulators: DgoR binds an operator sequence [5'-TTGTA(G/C)TACA(A/T)-3'] matching the signature of GntR family members that recognize inverted repeats [5'-(N) _y GT(N) _x AC(N) _y -3', where x and y indicate the number of nucleotides, which varies], and it shares critical protein-DNA contacts. We also identified features in DgoR that are otherwise less conserved in the GntR family. Recently, missense mutations in dgoR were recovered in a natural E. coli isolate adapted to the mammalian gut. Our results show these mutants to be DNA binding defective, emphasizing that mutations in the dgo-regulatory elements are selected in the host to allow simultaneous induction of dgo genes. The present study sets the basis to explore the regulation of dgo genes in additional enterobacterial strains where they have been implicated in host-bacterium interactions.IMPORTANCE d-Galactonate is a widely prevalent aldonic sugar acid. Despite the proposed significance of the d-galactonate metabolic pathway in the interaction of enteric bacteria with their hosts, there are no details on its regulation even in Escherichia coli, which has been known to utilize d-galactonate since the 1970s. Here, using multiple methodologies, we identified the promoter, operator, and effector of DgoR, the transcriptional repressor of d-galactonate metabolism in E. coli We establish DgoR as a GntR family transcriptional regulator. Recently, a human urinary tract isolate of E. coli introduced in the mouse gut was found to accumulate missense mutations in dgoR Our results show these mutants to be DNA binding defective, hence emphasizing the role of the d-galactonate metabolic pathway in bacterial colonization of the mammalian gut.

How Bioinformatic Tools Guide Experiments To Resolve the Chaos of Apparently Unlimited Metabolic Variation

Hexuronic acids, oxidation products of common sugars, are widespread in eukaryotic cells. Galacturonic acid is the main carbohydrate component of pectin found in plant cell walls and glucuronic acid is a component of proteoglycans in animals. However, despite their importance as carbohydrate substrates, metabolism of hexuronic acids has long remained a poorly studied corner of the bacterial metabolic map. In the current issue of Journal of Bacteriology, Bouvier and coworkers present a detailed analysis of genes involved in hexuronate utilization in various proteobacteria and report the verification of their bioinformatics predictions by carefully designed experiments (J. T. Bouvier et al., J Bacteriol 201:e00431-18, 2019, https://doi.org/10.1128/JB.00431-18). This study provides a solid basis for understanding hexuronate metabolism and its regulation in other bacterial phyla.

The Role of OmpR in the Expression of Genes of the KdgR Regulon Involved in the Uptake and Depolymerization of Oligogalacturonides in Yersinia enterocolitica

Oligogalacturonide (OGA)-specific porins of the KdgM family have previously been identified and characterized in enterobacterial plant pathogens. We found that deletion of the gene encoding response regulator OmpR causes the porin KdgM2 to become one of the most abundant proteins in the outer membrane of the human enteropathogen Yersinia enterocolitica. Reporter gene fusion and real-time PCR analysis confirmed that the expression of kdgM2 is repressed by OmpR. We also found that kdgM2 expression is subject to negative regulation by KdgR, a specific repressor of genes involved in the uptake and metabolism of pectin derivatives in plant pathogens. The additive effect of kdgR and ompR mutations suggested that KdgR and OmpR regulate kdgM2 expression independently. We confirmed that kdgM2 occurs in an operon with the pelP gene, encoding the periplasmic pectate lyase PelP. A pectinolytic assay showed strong upregulation of PelP production/activity in a Y. enterocolitica strain lacking OmpR and KdgR, which corroborates the repression exerted by these regulators on kdgM2. In addition, our data showed that OmpR is responsible for up regulation of the kdgM1 gene encoding the second specific oligogalacturonide porin KdgM1. This indicates the involvement of OmpR in the reciprocal regulation of both KdgM1 and KdgM2. Moreover, we demonstrated the negative impact of OmpR on kdgR transcription, which might positively affect the expression of genes of the KdgR regulon. Binding of OmpR to the promoter regions of the kdgM2-pelP-sghX operon, and kdgM1 and kdgR genes was confirmed using the electrophoretic mobility shift assay, suggesting that OmpR can directly regulate their transcription. We also found that the overexpression of porin KdgM2 increases outer membrane permeability. Thus, OmpR-mediated regulation of the KdgM porins may contribute to the fitness of Y. enterocolitica in particular local environments.

Comparative genomics and transcriptomics of Pichia pastoris

Background

Pichia pastoris has emerged as an important alternative host for producing recombinant biopharmaceuticals, owing to its high cultivation density, low host cell protein burden, and the development of strains with humanized glycosylation. Despite its demonstrated utility, relatively little strain engineering has been performed to improve Pichia, due in part to the limited number and inconsistent frameworks of reported genomes and transcriptomes. Furthermore, the co-mingling of genomic, transcriptomic and fermentation data collected about Komagataella pastoris and Komagataella phaffii, the two strains co-branded as Pichia, has generated confusion about host performance for these genetically distinct species. Generation of comparative high-quality genomes and transcriptomes will enable meaningful comparisons between the organisms, and potentially inform distinct biotechnological utilies for each species.

Results

Here, we present a comprehensive and standardized comparative analysis of the genomic features of the three most commonly used strains comprising the tradename Pichia: K. pastoris wild-type, K. phaffii wild-type, and K. phaffii GS115. We used a combination of long-read (PacBio) and short-read (Illumina) sequencing technologies to achieve over 1000X coverage of each genome. Construction of individual genomes was then performed using as few as seven individual contigs to create gap-free assemblies. We found substantial syntenic rearrangements between the species and characterized a linear plasmid present in K. phaffii. Comparative analyses between K. phaffii genomes enabled the characterization of the mutational landscape of the GS115 strain. We identified and examined 35 non-synonomous coding mutations present in GS115, many of which are likely to impact strain performance. Additionally, we investigated transcriptomic profiles of gene expression for both species during cultivation on various carbon sources. We observed that the most highly transcribed genes in both organisms were consistently highly expressed in all three carbon sources examined. We also observed selective expression of certain genes in each carbon source, including many sequences not previously reported as promoters for expression of heterologous proteins in yeasts.

Conclusions

Our studies establish a foundation for understanding critical relationships between genome structure, cultivation conditions and gene expression. The resources we report here will inform and facilitate rational, organism-wide strain engineering for improved utility as a host for protein production.

Electronic supplementary material

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

Catabolism of Hexuronides, Hexuronates, Aldonates, and Aldarates

Following elucidation of the regulation of the lactose operon in Escherichia coli, studies on the metabolism of many sugars were initiated in the early 1960s. The catabolic pathways of D-gluconate and of the two hexuronates, D-glucuronate and D-galacturonate, were investigated. The post genomic era has renewed interest in the study of these sugar acids and allowed the complete characterization of the D-gluconate pathway and the discovery of the catabolic pathways for L-idonate, D-glucarate, galactarate, and ketogluconates. Among the various sugar acids that are utilized as sole carbon and energy sources to support growth of E. coli, galacturonate, glucuronate, and gluconate were shown to play an important role in the colonization of the mammalian large intestine. In the case of sugar acid degradation, the regulators often mediate negative control and are inactivated by interaction with a specific inducer, which is either the substrate or an intermediate of the catabolism. These regulators coordinate the synthesis of all the proteins involved in the same pathway and, in some cases, exert crosspathway control between related catabolic pathways. This is particularly well illustrated in the case of hexuronide and hexuronate catabolism. The structural genes encoding the different steps of hexuronate catabolism were identified by analysis of numerous mutants affected for growth with galacturonate or glucuronate. E. coli is able to use the diacid sugars D-glucarate and galactarate (an achiral compound) as sole carbon source for growth. Pyruvate and 2-phosphoglycerate are the final products of the D-glucarate/galactarate catabolism.

Signature-tagged mutagenesis and co-infection studies demonstrate the importance of P fimbriae in a murine model of urinary tract infection

Escherichia coli is the leading cause of urinary tract infections (UTIs), one of the most common infections in humans. P fimbria was arguably the first proposed virulence factor for uropathogenic E. coli, based on the capacity of E. coli isolated from UTIs to adhere to exfoliated epithelial cells in higher numbers than fecal strains of E. coli. Overwhelming epidemiologic evidence has been presented for involvement of P fimbriae in colonization. It has been difficult, however, to demonstrate this requirement for uropathogenic strains in animal models of infections or in humans. In this study, a signature-tagged mutagenesis screen identified a P-fimbrial gene (papC) and 18 other genes as being among those required for full fitness of cystitis isolate E. coli F11. A P-fimbrial mutant was outcompeted by the wild-type strain in cochallenge in the murine model of ascending UTI, and this colonization defect could be complemented with the cloned pap operon. To our knowledge, this study is the first to fulfill molecular Koch's postulates in which a pathogenic strain was attenuated by mutation of pap genes and then complemented to restore fitness, confirming P fimbria as a virulence factor in a pathogenic clinical isolate.

NAD captureSeq indicates NAD as a bacterial cap for a subset of regulatory RNAs

A distinctive feature of prokaryotic gene expression is the absence of 5'-capped RNA. In eukaryotes, 5',5'-triphosphate-linked 7-methylguanosine protects messenger RNA from degradation and modulates maturation, localization and translation. Recently, the cofactor nicotinamide adenine dinucleotide (NAD) was reported as a covalent modification of bacterial RNA. Given the central role of NAD in redox biochemistry, posttranslational protein modification and signalling, its attachment to RNA indicates that there are unknown functions of RNA in these processes and undiscovered pathways in RNA metabolism and regulation. The unknown identity of NAD-modified RNAs has so far precluded functional analyses. Here we identify NAD-linked RNAs from bacteria by chemo-enzymatic capture and next-generation sequencing (NAD captureSeq). Among those identified, specific regulatory small RNAs (sRNAs) and sRNA-like 5'-terminal fragments of certain mRNAs are particularly abundant. Analogous to a eukaryotic cap, 5'-NAD modification is shown in vitro to stabilize RNA against 5'-processing by the RNA-pyrophosphohydrolase RppH and against endonucleolytic cleavage by ribonuclease (RNase) E. The nudix phosphohydrolase NudC decaps NAD-RNA and thereby triggers RNase-E-mediated RNA decay, while being inactive against triphosphate-RNA. In vivo, ∼13% of the abundant sRNA RNAI is NAD-capped in the presence, and ∼26% in the absence, of functional NudC. To our knowledge, this is the first description of a cap-like structure and a decapping machinery in bacteria.

Identification of host-microbe interaction factors in the genomes of soft rot-associated pathogens Dickeya dadantii 3937 and Pectobacterium carotovorum WPP14 with supervised machine learning

Background

A wealth of genome sequences has provided thousands of genes of unknown function, but identification of functions for the large numbers of hypothetical genes in phytopathogens remains a challenge that impacts all research on plant-microbe interactions. Decades of research on the molecular basis of pathogenesis focused on a limited number of factors associated with long-known host-microbe interaction systems, providing limited direction into this challenge. Computational approaches to identify virulence genes often rely on two strategies: searching for sequence similarity to known host-microbe interaction factors from other organisms, and identifying islands of genes that discriminate between pathogens of one type and closely related non-pathogens or pathogens of a different type. The former is limited to known genes, excluding vast collections of genes of unknown function found in every genome. The latter lacks specificity, since many genes in genomic islands have little to do with host-interaction.

Result

In this study, we developed a supervised machine learning approach that was designed to recognize patterns from large and disparate data types, in order to identify candidate host-microbe interaction factors. The soft rot Enterobacteriaceae strains Dickeya dadantii 3937 and Pectobacterium carotovorum WPP14 were used for development of this tool, because these pathogens are important on multiple high value crops in agriculture worldwide and more genomic and functional data is available for the Enterobacteriaceae than any other microbial family. Our approach achieved greater than 90% precision and a recall rate over 80% in 10-fold cross validation tests.

Conclusion

Application of the learning scheme to the complete genome of these two organisms generated a list of roughly 200 candidates, many of which were previously not implicated in plant-microbe interaction and many of which are of completely unknown function. These lists provide new targets for experimental validation and further characterization, and our approach presents a promising pattern-learning scheme that can be generalized to create a resource to study host-microbe interactions in other bacterial phytopathogens.

Electronic supplementary material

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

Polysaccharides utilization in human gut bacterium Bacteroides thetaiotaomicron: comparative genomics reconstruction of metabolic and regulatory networks

Background

Bacteroides thetaiotaomicron, a predominant member of the human gut microbiota, is characterized by its ability to utilize a wide variety of polysaccharides using the extensive saccharolytic machinery that is controlled by an expanded repertoire of transcription factors (TFs). The availability of genomic sequences for multiple Bacteroides species opens an opportunity for their comparative analysis to enable characterization of their metabolic and regulatory networks.

Results

A comparative genomics approach was applied for the reconstruction and functional annotation of the carbohydrate utilization regulatory networks in 11 Bacteroides genomes. Bioinformatics analysis of promoter regions revealed putative DNA-binding motifs and regulons for 31 orthologous TFs in the Bacteroides. Among the analyzed TFs there are 4 SusR-like regulators, 16 AraC-like hybrid two-component systems (HTCSs), and 11 regulators from other families. Novel DNA motifs of HTCSs and SusR-like regulators in the Bacteroides have the common structure of direct repeats with a long spacer between two conserved sites.

Conclusions

The inferred regulatory network in B. thetaiotaomicron contains 308 genes encoding polysaccharide and sugar catabolic enzymes, carbohydrate-binding and transport systems, and TFs. The analyzed TFs control pathways for utilization of host and dietary glycans to monosaccharides and their further interconversions to intermediates of the central metabolism. The reconstructed regulatory network allowed us to suggest and refine specific functional assignments for sugar catabolic enzymes and transporters, providing a substantial improvement to the existing metabolic models for B. thetaiotaomicron. The obtained collection of reconstructed TF regulons is available in the RegPrecise database (http://regprecise.lbl.gov).

Novel insights into E. coli's hexuronate metabolism: KduI facilitates the conversion of galacturonate and glucuronate under osmotic stress conditions

Using a gnotobiotic mouse model, we previously observed the upregulation of 2-deoxy-D-gluconate 3-dehydrogenase (KduD) in intestinal E. coli of mice fed a lactose-rich diet and the downregulation of this enzyme and of 5-keto 4-deoxyuronate isomerase (KduI) on a casein-rich diet. The present study aimed to define the role of the so far poorly characterized E. coli proteins KduD and KduI in vitro. Galacturonate and glucuronate induced kduD and kduI gene expression 3-fold and 7 to 11-fold, respectively, under aerobic conditions as well as 9 to 20-fold and 19 to 54-fold, respectively, under anaerobic conditions. KduI facilitated the breakdown of these hexuronates. In E. coli, galacturonate and glucuronate are normally degraded by UxaABC and UxuAB. However, osmotic stress represses the expression of the corresponding genes in an OxyR-dependent manner. When grown in the presence of galacturonate or glucuronate, kduID-deficient E. coli had a 30% to 80% lower maximal cell density and 1.5 to 2-fold longer doubling times under osmotic stress conditions than wild type E. coli. Growth on lactose promoted the intracellular formation of hexuronates, which possibly explain the induction of KduD on a lactose-rich diet. These results indicate a novel function of KduI and KduD in E. coli and demonstrate the crucial influence of osmotic stress on the gene expression of hexuronate degrading enzymes.

The Entner-Doudoroff pathway is obligatory for gluconate utilization and contributes to the pathogenicity of Vibrio cholerae

The Entner-Doudoroff (ED) pathway has recently been shown to play an important role in sugar catabolism for many organisms although very little information is available on the functionality of this pathway in Vibrio cholerae, the causative agent of cholera. In this study, activation of the genes edd and eda, encoding 6-phosphogluconate dehydratase and 2-keto-3-deoxy-6-phosphogluconate aldolase, was used as a marker of a functional ED pathway in V. cholerae. Transcriptional activation analyses and gene silencing experiments with cells grown in sugar-supplemented M9 medium demonstrated that the ED pathway is functional in V. cholerae and is obligatory for gluconate catabolism. Importantly, selective activation of the ED pathway led to concurrent elevation of transcripts of prime virulence genes (ctxA and tcpA) and their regulator (toxT). Further, lowering of these transcript levels and cholera toxin production in vitro by an ED pathway-defective mutant (strain N16961 with a Δedd mutation [Δedd(N16961) strain]) suggested the importance of this pathway in regulating V. cholerae virulence. The in vivo relevance of these data was established as the mutant failed to colonize in suckling mice intestine or to induce fluid accumulation in ligated rabbit ileal loops. Activation of the ED pathway in V. cholerae was shown to inhibit biofilm formation in vitro that could be reversed in the mutant. As further support for these results, comparative transcriptome analysis with cells grown in the presence of glucose or gluconate revealed that a functional ED pathway led to activation of a subset of previously reported in vivo expressed genes. All of these results suggest the importance of the ED pathway in V. cholerae pathogenesis.

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.

Comparative genomic analysis of the hexuronate metabolism genes and their regulation in gammaproteobacteria

The hexuronate metabolism in Escherichia coli is regulated by two related transcription factors from the FadR subfamily of the GntR family, UxuR and ExuR. UxuR controls the d-glucuronate metabolism, while ExuR represses genes involved in the metabolism of all hexuronates. We use a comparative genomics approach to reconstruct the hexuronate metabolic pathways and transcriptional regulons in gammaproteobacteria. We demonstrate differences in the binding motifs of UxuR and ExuR, identify new candidate members of the UxuR/ExuR regulons, and describe the links between the UxuR/ExuR regulons and the adjacent regulons UidR, KdgR, and YjjM. We provide experimental evidence that two predicted members of the UxuR regulon, yjjM and yjjN, are the subject of complex regulation by this transcription factor in E. coli.

Untangling metabolic and communication networks: interactions of enterics with phytobacteria and their implications in produce safety

Recent outbreaks of vegetable-borne gastrointestinal illnesses across the globe demonstrate that human enteric pathogens can contaminate produce at any stage of production. Interactions of enterics with native plant-associated microbiota influence the microbiological safety of produce by affecting the attachment, persistence and proliferation of human pathogens on plants. Supermarket surveys have revealed that bacteria, but not fungi or mechanical damage, promote the growth of Salmonella enterica on produce. Field and laboratory studies have indicated that some plant pathogenic bacteria and fungi facilitate the entry and internalization of human pathogens in plants. Conversely, some phytobacteria, including those involved in biocontrol of plant diseases, significantly inhibit attachment and plant colonization by non-typhoidal Salmonella and enterovirulent Escherichia coli by producing antibiotics or competing for nutrients in the phyllosphere. In this review, we attempt to elucidate the mechanisms of interactions between human enteric pathogens and plant-associated microbiota, and describe how these interactions affect produce safety.

Comparative proteomic and genetic analyses reveal unidentified mutations in Escherichia coli XL1-Blue and DH5α

Escherichia coli has been used widely in laboratory and the biotech industry. However, the genetic and metabolic characteristics remain inadequately studied, particularly for those strains with extensive genetic manipulations that might have resulted in unknown mutations. Here, we demonstrate a comparative proteomics and genetics approach to identify unknown mutations in E. coli K-12 derivatives. The comparative proteomic and genetic analyses revealed an IS5 disruption of the kdgR gene in two commonly used derivative strains of E. coli K-12, XL1-Blue and DH5α, compared with K-12 wild-type strain W3110. In addition, a controversial deoR mutation was clarified as a wild type in E. coli DH5α using the same approach. This approach should be useful in characterizing the unknown mutations in various mutant strains developed. At the same time, comparative proteomic analysis also revealed the distinct metabolic characteristic of the two derivatives: higher biosynthetic flux to purine nucleotides. This is potentially beneficial for the synthesis of plasmid DNA.

The active site of oligogalacturonate lyase provides unique insights into cytoplasmic oligogalacturonate beta-elimination

Oligogalacturonate lyases (OGLs; now also classified as pectate lyase family 22) are cytoplasmic enzymes found in pectinolytic members of Enterobacteriaceae, such as the enteropathogen Yersinia enterocolitica. OGLs utilize a β-elimination mechanism to preferentially catalyze the conversion of saturated and unsaturated digalacturonate into monogalacturonate and the 4,5-unsaturated monogalacturonate-like molecule, 5-keto-4-deoxyuronate. To provide mechanistic insights into the specificity of this enzyme activity, we have characterized the OGL from Y. enterocolitica, YeOGL, on oligogalacturonides and determined its three-dimensional x-ray structure to 1.65 Å. The model contains a Mn(2+) atom in the active site, which is coordinated by three histidines, one glutamine, and an acetate ion. The acetate mimics the binding of the uronate group of galactourono-configured substrates. These findings, in combination with enzyme kinetics and metal supplementation assays, provide a framework for modeling the active site architecture of OGL. This enzyme appears to contain a histidine for the abstraction of the α-proton in the -1 subsite, a residue that is highly conserved throughout the OGL family and represents a unique catalytic base among pectic active lyases. In addition, we present a hypothesis for an emerging relationship observed between the cellular distribution of pectate lyase folding and the distinct metal coordination chemistries of pectate lyases.

Engineering transcription factors with novel DNA-binding specificity using comparative genomics

The transcriptional program for a gene consists of the promoter necessary for recruiting RNA polymerase along with neighboring operator sites that bind different activators and repressors. From a synthetic biology perspective, if the DNA-binding specificity of these proteins can be changed, then they can be used to reprogram gene expression in cells. While many experimental methods exist for generating such specificity-altering mutations, few computational approaches are available, particularly in the case of bacterial transcription factors. In a previously published computational study of nitrogen oxide metabolism in bacteria, a small number of amino-acid residues were found to determine the specificity within the CRP (cAMP receptor protein)/FNR (fumarate and nitrate reductase regulatory protein) family of transcription factors. By analyzing how these amino acids vary in different regulators, a simple relationship between the identity of these residues and their target DNA-binding sequence was constructed. In this article, we experimentally tested whether this relationship could be used to engineer novel DNA–protein interactions. Using Escherichia coli CRP as a template, we tested eight designs based on this relationship and found that four worked as predicted. Collectively, these results in this work demonstrate that comparative genomics can inform the design of bacterial transcription factors.

Structural biology of pectin degradation by Enterobacteriaceae

SUMMARY

Pectin is a structural polysaccharide that is integral for the stability of plant cell walls. During soft rot infection, secreted virulence factors from pectinolytic bacteria such as Erwinia spp. degrade pectin, resulting in characteristic plant cell necrosis and tissue maceration. Catabolism of pectin and its breakdown products by pectinolytic bacteria occurs within distinct cellular environments. This process initiates outside the cell, continues within the periplasmic space, and culminates in the cytoplasm. Although pectin utilization is well understood at the genetic and biochemical levels, an inclusive structural description of pectinases and pectin binding proteins by both extracellular and periplasmic enzymes has been lacking, especially following the recent characterization of several periplasmic components and protein-oligogalacturonide complexes. Here we provide a comprehensive analysis of the protein folds and mechanisms of pectate lyases, polygalacturonases, and carbohydrate esterases and the binding specificities of two periplasmic pectic binding proteins from Enterobacteriaceae. This review provides a structural understanding of the molecular determinants of pectin utilization and the mechanisms driving catabolite selectivity and flow through the pathway.

Regulation of the kduID operon of Bacillus subtilis by the KdgR repressor and the ccpA gene: identification of two KdgR-binding sites within the kdgR-kduI intergenic region

Transcription of the Bacillus subtilis kdgRKAT operon, which comprises genes involved in the late stage of galacturonate utilization, is known to be negatively regulated by the KdgR repressor. In this study, Northern analysis was carried out to demonstrate that the kdgR gene also negatively regulates the kduID operon, encoding ketodeoxyuronate isomerase and ketodeoxygluconate reductase. It has also been demonstrated that expression of the kduID operon can be induced by galacturonate and is subject to catabolite repression by glucose. The ccpA gene was found to be involved in this catabolite repression. Primer extension analysis identified a sigma(A)-like promoter sequence preceding kduI. Gel mobility shift assays and DNase I footprinting analyses indicated that KdgR is capable of binding specifically to two sites within the kdgR-kduI intergenic region in vitro. Reporter gene analysis revealed that these two KdgR-binding sites function in vivo. One site is centred 33.5 bp upstream of the translational start site of kdgR and can serve as an operator for controlling expression of the kdgRKAT operon. The other is centred 57.5 bp upstream of the translational start site of kduI and can serve as an operator for controlling expression of the kduID operon. Possible physiological significance of this regulation is discussed.

PhyloScan: identification of transcription factor binding sites using cross-species evidence

Background

When transcription factor binding sites are known for a particular transcription factor, it is possible to construct a motif model that can be used to scan sequences for additional sites. However, few statistically significant sites are revealed when a transcription factor binding site motif model is used to scan a genome-scale database.

Methods

We have developed a scanning algorithm, PhyloScan, which combines evidence from matching sites found in orthologous data from several related species with evidence from multiple sites within an intergenic region, to better detect regulons. The orthologous sequence data may be multiply aligned, unaligned, or a combination of aligned and unaligned. In aligned data, PhyloScan statistically accounts for the phylogenetic dependence of the species contributing data to the alignment and, in unaligned data, the evidence for sites is combined assuming phylogenetic independence of the species. The statistical significance of the gene predictions is calculated directly, without employing training sets.

Results

In a test of our methodology on synthetic data modeled on seven Enterobacteriales, four Vibrionales, and three Pasteurellales species, PhyloScan produces better sensitivity and specificity than MONKEY, an advanced scanning approach that also searches a genome for transcription factor binding sites using phylogenetic information. The application of the algorithm to real sequence data from seven Enterobacteriales species identifies novel Crp and PurR transcription factor binding sites, thus providing several new potential sites for these transcription factors. These sites enable targeted experimental validation and thus further delineation of the Crp and PurR regulons in E. coli.

Conclusion

Better sensitivity and specificity can be achieved through a combination of (1) using mixed alignable and non-alignable sequence data and (2) combining evidence from multiple sites within an intergenic region.

Evolution of the NadR regulon in Enterobacteriaceae

The NAD biosynthetic pathway and NAD transformations in E. coli and S. typhi are well characterized. Using comparative genomics methods we describe the NadR regulon in other Enterobacteriaceae, identity new candidate regulon members and demonstrate that even a very simple regulon covering an essential methabolic pathway could be different in closely related genomes.

Glycolysis and Flux Control

Central metabolism of carbohydrates uses the Embden-Meyerhof-Parnas (EMP), pentose phosphate (PP), and Entner-Doudoroff (ED) pathways. This review reviews the biological roles of the enzymes and genes of these three pathways of E. coli. Glucose, pentoses, and gluconate are primarily discussed as the initial substrates of the three pathways, respectively. The genetic and allosteric regulatory mechanisms of glycolysis and the factors that affect metabolic flux through the pathways are considered here. Despite the fact that a lot of information on each of the reaction steps has been accumulated over the years for E. coli, surprisingly little quantitative information has been integrated to analyze glycolysis as a system. Therefore, the review presents a detailed description of each of the catalytic steps by a systemic approach. It considers both structural and kinetic aspects. Models that include kinetic information of the reaction steps will always contain the reaction stoichiometry and therefore follow the structural constraints, but in addition to these also kinetic rate laws must be fulfilled. The kinetic information obtained on isolated enzymes can be integrated using computer models to simulate behavior of the reaction network formed by these enzymes. Successful examples of such approaches are the modeling of glycolysis in S. cerevisiae, the parasite Trypanosoma brucei, and the red blood cell. With the rapid developments in the field of Systems Biology many new methods have been and will be developed, for experimental and theoretical approaches, and the authors expect that these will be applied to E. coli glycolysis in the near future.

Comparative genomics of the KdgR regulon in Erwinia chrysanthemi 3937 and other gamma-proteobacteria

In the plant-pathogenic enterobacterium Erwinia chrysanthemi, almost all known genes involved in pectin catabolism are controlled by the transcriptional regulator KdgR. In this study, the comparative genomics approach was used to analyse the KdgR regulon in completely sequenced genomes of eight enterobacteria, including Erw. chrysanthemi, and two Vibrio species. Application of a signal recognition procedure complemented by operon structure and protein sequence analysis allowed identification of new candidate genes of the KdgR regulon. Most of these genes were found to be controlled by the cAMP-receptor protein, a global regulator of catabolic genes. At the next step, regulation of these genes in Erw. chrysanthemi was experimentally verified using in vivo transcriptional fusions and an attempt was made to clarify the functional role of the predicted genes in pectin catabolism. Interestingly, it was found that the KdgR protein, previously known as a repressor, positively regulates expression of two new members of the regulon, phosphoenolpyruvate synthase gene ppsA and an adjacent gene, ydiA, of unknown function. Other predicted regulon members, namely chmX, dhfX, gntB, pykF, spiX, sotA, tpfX, yeeO and yjgK, were found to be subject to classical negative regulation by KdgR. Possible roles of newly identified members of the Erw. chrysanthemi KdgR regulon, chmX, dhfX, gntDBMNAC, spiX, tpfX, ydiA, yeeO, ygjV and yjgK, in pectin catabolism are discussed. Finally, complete reconstruction of the KdgR regulons in various gamma-proteobacteria yielded a metabolic map reflecting a globally conserved pathway for the catabolism of pectin and its derivatives with variability in transport and enzymic capabilities among species. In particular, possible non-orthologous substitutes of isomerase KduI and a new oligogalacturonide transporter in the Vibrio species were detected.

Multiple regulators control expression of the Entner-Doudoroff aldolase (Eda) of Escherichia coli

The Escherichia coli eda gene, which encodes the Entner-Doudoroff aldolase, is central to the catabolism of several sugar acids. Here, we show that Eda synthesis is induced by growth on gluconate, glucuronate, or methyl-beta-D-glucuronide; phosphate limitation; and carbon starvation. Transcription of eda initiates from three promoters, designated P1, P2, and P4, each of which is responsible for induction under different growth conditions. P1 controls eda induction on gluconate and is regulated by GntR. P2 controls eda induction on glucuronate and galacturonate and is regulated by KdgR. P4 is active under conditions of phosphate starvation and is directly controlled by PhoB. In addition, CsrA activates Eda synthesis, apparently by an indirect mechanism that may be involved in the modest changes in expression level that are associated with carbon starvation. The complex regulation of eda is discussed with respect to its several physiological roles, which apparently accommodate not only sugar acid catabolism but also detoxification of metabolites that could accumulate during starvation-induced stress.

Evolution of transcription factor DNA binding sites

In bioinformatics, binding of transcription regulatory factors to the cognate binding sites is usually described by sequence-specific binding energy, which is estimated from a training sample of sites. This model implies that all binding sites with binding energy above some threshold are functional and site sequence variations should be considered neutral until they do not reduce this energy below the threshold. To quantify this energy, the binding profile (positional weight matrix, PWM) model or consensus-based model is usually applied. Here we show that in many cases available data are not sufficient to construct a relevant PWM, and modified consensus-based model could be more effective to describe binding properties. Further, using the data about binding sites of several transcription factors, we demonstrate that some non-consensus nucleotides in "orthologous sites" (that is, binding sites of the same factor upstream of orthologous genes), which have been believed to be irrelevant or even hindering the regulation, are evolutionary very stable and specific for the regulated gene. For each two considered genomes, the number of substitutions between non-consensus nucleotides is far less than the expected number of neutral substitutions. Moreover, in several positions of binding sites regulating different genes, there are non-consensus nucleotides conserved in distant genomes. It means that there exists a selection pressure, which results in the stability of non-consensus nucleotides.

GntP is the Escherichia coli Fructuronic acid transporter and belongs to the UxuR regulon

Escherichia coli has four gluconate transporters, GntP, GntU, GntT, and IdnT, which are members of the major facilitator superfamily. The physiological function of GntP was previously unknown and is the subject of this study. GntP is not induced by gluconate, and despite being located adjacent to genes involved in glucuronate catabolism, gntP does not encode a glucuronate transporter. Here we identify gntP as the gene which encodes the fructuronate transporter. We show that gntP is induced by fructuronate and is a new member of the UxuR regulon: gntP is derepressed in an uxuR strain, UxuR binds in vitro specifically to an operator site that overlaps the gntP promoter, and UxuR binding is eliminated by fructuronate. Transcription of gntP requires activation by cyclic AMP (cAMP)-cAMP receptor protein. A gntP mutant cannot grow on fructuronate but grows normally on glucuronate and gluconate. Thus, the UxuR regulon is a module of sugar acid catabolism whose physiological role is for growth on fructuronate. Glucuronate, because it proceeds through a fructuronate intermediate, must induce the UxuR regulon and must also induce the ExuR regulon, which encodes the glucuronate transporter, ExuT, and the first step in its catabolism, UxaC. Thus, hexuronate catabolism in E. coli requires both the ExuR and UxuR regulons, while fructuronate catabolism requires only the UxuR regulon.

Reconstruction of regulatory and metabolic pathways in metal-reducing delta-proteobacteria

A study of the genetic and regulatory factors in several biosynthesis, metal ion homeostasis, stress response, and energy metabolism pathways suggests that phylogenetically diverse δ-proteobacteria have homologous regulatory components.

Background

Relatively little is known about the genetic basis for the unique physiology of metal-reducing genera in the delta subgroup of the proteobacteria. The recent availability of complete finished or draft-quality genome sequences for seven representatives allowed us to investigate the genetic and regulatory factors in a number of key pathways involved in the biosynthesis of building blocks and cofactors, metal-ion homeostasis, stress response, and energy metabolism using a combination of regulatory sequence detection and analysis of genomic context.

Results

In the genomes of δ-proteobacteria, we identified candidate binding sites for four regulators of known specificity (BirA, CooA, HrcA, sigma-32), four types of metabolite-binding riboswitches (RFN-, THI-, B12-elements and S-box), and new binding sites for the FUR, ModE, NikR, PerR, and ZUR transcription factors, as well as for the previously uncharacterized factors HcpR and LysX. After reconstruction of the corresponding metabolic pathways and regulatory interactions, we identified possible functions for a large number of previously uncharacterized genes covering a wide range of cellular functions.

Conclusions

Phylogenetically diverse δ-proteobacteria appear to have homologous regulatory components. This study for the first time demonstrates the adaptability of the comparative genomic approach to de novo reconstruction of a regulatory network in a poorly studied taxonomic group of bacteria. Recent efforts in large-scale functional genomic characterization of Desulfovibrio species will provide a unique opportunity to test and expand our predictions.

Factors influencing the identification of transcription factor binding sites by cross-species comparison

As the number of sequenced genomes has grown, the questions of which species are most useful and how many genomes are sufficient for comparison have become increasingly important for comparative genomics studies. We have systematically addressed these questions with respect to phylogenetic footprinting of transcription factor (TF) binding sites in the gamma-proteobacteria, and have evaluated the statistical significance of our motif predictions. We used a study set of 166 Escherichia coli genes that have experimentally identified TF binding sites upstream of the gene, with orthologous data from nine additional gamma-proteobacteria for phylogenetic footprinting. Just three species were sufficient for approximately 74.0% of the motif predictions to correspond to the experimentally reported E. coli sites, and important characteristics to consider when choosing species were phylogenetic distance, genome size, and natural habitat. We also performed simulations using randomized data to determine the critical maximum a posteriori probability (MAP) values for statistical significance of our motif predictions (P = 0.05). Approximately 60% of motif predictions containing sites from just three species had average MAP values above these critical MAP values. The inclusion of a species very closely related to E. coli increased the number of statistically significant motif predictions, despite substantially increasing the critical MAP value.

Subdivision of the helix-turn-helix GntR family of bacterial regulators in the FadR, HutC, MocR, and YtrA subfamilies

Haydon and Guest (Haydon, D. J, and Guest, J. R. (1991) FEMS Microbiol. Lett. 63, 291-295) first described the helix-turn-helix GntR family of bacterial regulators. They presented them as transcription factors sharing a similar N-terminal DNA-binding (d-b) domain, but they observed near-maximal divergence in the C-terminal effector-binding and oligomerization (E-b/O) domain. To elucidate this C-terminal heterogeneity, structural, phylogenetic, and functional analyses were performed on a family that now comprises about 270 members. Our comparative study first focused on the C-terminal E-b/O domains and next on DNA-binding domains and palindromic operator sequences, has classified the GntR members into four subfamilies that we called FadR, HutC, MocR, and YtrA. Among these subfamilies a degree of similarity of about 55% was observed throughout the entire sequence. Structure/function associations were highlighted although they were not absolutely stringent. The consensus sequences deduced for the DNA-binding domain were slightly different for each subfamily, suggesting that fusion between the D-b and E-b/O domains have occurred separately, with each subfamily having its own D-b domain ancestor. Moreover, the compilation of the known or predicted palindromic cis-acting elements has highlighted different operator sequences according to our subfamily subdivision. The observed C-terminal E-b/O domain heterogeneity was therefore reflected on the DNA-binding domain and on the cis-acting elements, suggesting the existence of a tight link between the three regions involved in the regulating process.

Transcriptional regulation of pentose utilisation systems in the Bacillus/Clostridium group of bacteria

In Bacillus subtilis, utilisation of xylose, arabinose and ribose is controlled by the transcriptional factors XylR, AraR and RbsR, respectively. Here we apply the comparative approach to the analysis of these regulons in the Bacillus/Clostridium group. Evolutionary variability of operon structures is demonstrated and operator sites for the main transcription factors are predicted. The consensus sequences for the XylR and RbsR binding sites vary in different subgroups of genomes. The functional coupling of gene clusters and the conservation of regulatory sites allow for detection of non-orthologous gene displacement of ribulose kinase in Enterococcus faecium and Clostridium acetobutylicum. Moreover, candidate catabolite responsive elements found upstream of most pentose-utilising genes suggest CcpA-mediated catabolite repression.

Comparative analysis of FUR regulons in gamma-proteobacteria

Iron is an essential element for the survival and pathogenesis of bacteria. The strict control of iron homeostasis is mediated by the FUR repressor, which is highly conserved among various bacterial species. Here we apply the comparative genomics approach to analyze candidate Fur-binding sites in the genomes of Escherichia coli (K12 and O157:H7), Salmonella typhi, Yersinia pestis and Vibrio cholerae. We describe a number of new loci encoding siderophore biosynthesis and transport proteins. A new regulator of iron-acquisition systems was found in S.typhi. We predict FUR regulation for several virulence systems. We also predict FUR regulation for the chemotaxis system of V.cholerae that is probably involved in the process of pathogenesis.

BphS, a key transcriptional regulator of bph genes involved in polychlorinated biphenyl/biphenyl degradation in Pseudomonas sp. KKS102

The bph genes in Pseudomonas sp. KKS102, which are involved in the degradation of polychlorinated biphenyl/biphenyl, are induced in the presence of biphenyl. In this study our goal was to understand the regulatory mechanisms involved in the inducible expression. The bph genes (bphEGF(orf4)A1A2A3BCD(orf1)A4R) constitute an operon, and its expression is strongly dependent on the pE promoter located upstream of the bphE gene. A bphS gene, whose deduced amino acid sequence showed homology with the GntR family transcriptional repressors, was identified at the upstream region of the bphE gene. Disruption of the bphS gene resulted in constitutive expression of bph genes, suggesting that the bphS gene product negatively regulated the pE promoter. The gel retardation and DNase footprinting analyses demonstrated specific binding of BphS to the pE promoter region and identified four BphS binding sites that were located within and immediately downstream of the -10 box of the pE promoter. The four binding sites were functional in repression because their respective elimination resulted in derepression of the pE promoter. The binding of BphS was abolished in the presence of 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid, an intermediate compound in the biphenyl degradation pathway. We concluded that the negative regulator BphS plays a central role in the regulation of bph gene expression through its action at the pE promoter.

Identification of TogMNAB, an ABC transporter which mediates the uptake of pectic oligomers in Erwinia chrysanthemi 3937

The bacterium Erwinia chrysanthemi, which causes soft rot disease on various plants, is able to use pectin as a carbon source for growth. Knowledge of the critical step in pectin catabolism which allows the entry of pectic oligomers into the cells is scarce. We report here the first example of a transport system involved in the uptake of pectic oligomers. The TogMNAB transporter of E. chrysanthemi is a member of the ATP-binding cassette (ABC) superfamily. TogM and TogN are homologous to the inner membrane components, TogA exhibits the signature of ABC ATPases and TogB shows similarity with periplasmic ligand-binding proteins. The TogMNAB transporter is a new member of the carbohydrate uptake transporter-1 family (CUT1, TC no. 3.1.1), which is specialized in the transport of complex sugars. The four genes, togM, togN, togA and togB, are apparently co-transcribed in a large operon which also includes the pectate lyase gene pelW. The transcription of the tog operon is induced in the presence of pectic derivatives and is affected by catabolite repression. It is controlled by the KdgR repressor and the CRP activator. The TogMNAB system is able to provide Escherichia coli with the ability to transport oligogalacturonides. In E. chrysanthemi, the TogMNAB system seems to play a major role in switching on the induction of pectin catabolism. TogB also acts as a specific receptor for chemotaxis towards oligogalacturonides. The decreased capacity of maceration of a togM mutant indicates the importance of transport and/or attraction of oligogalacturonides for E. chrysanthemi pathogenicity.