Analysis of an Erwinia chrysanthemi gene cluster involved in pectin degradation

Genome-Wide Identification of Genes Important for Growth of Dickeya dadantii and Dickeya dianthicola in Potato (Solanum tuberosum) Tubers

Dickeya species are causal agents of soft rot diseases in many economically important crops, including soft rot disease of potato (Solanum tuberosum). Using random barcode transposon-site sequencing (RB-TnSeq), we generated genome-wide mutant fitness profiles of Dickeya dadantii 3937, Dickeya dianthicola ME23, and Dickeya dianthicola 67-19 isolates collected after passage through several in vitro and in vivo conditions. Though all three strains are pathogenic on potato, D. dadantii 3937 is a well-characterized model while D. dianthicola strains ME23 and 67-19 are recent isolates. Strain ME23 specifically was identified as a representative strain from a 2014 outbreak on potato. This study generated comparable gene fitness measurements across ecologically relevant conditions for both model and non-model strains. Tubers from the potato cultivars “Atlantic,” “Dark Red Norland,” and “Upstate Abundance” provided highly similar conditions for bacterial growth. Using the homolog detection software PyParanoid, we matched fitness values for orthologous genes in the three bacterial strains. Direct comparison of fitness among the strains highlighted shared and variable traits important for growth. Bacterial growth in minimal medium required many metabolic traits that were also essential for competitive growth in planta, such as amino acid, carbohydrate, and nucleotide biosynthesis. Growth in tubers specifically required the pectin degradation gene kduD. Disruption in three putative DNA-binding proteins had strain-specific effects on competitive fitness in tubers. Though the Soft Rot Pectobacteriaceae can cause disease with little host specificity, it remains to be seen the extent to which strain-level variation impacts virulence.

Separation and quantification of 2-keto-3-deoxy-gluconate (KDG) a major metabolite in pectin and alginate degradation pathways

A rapid and sensitive High Performance Liquid Chromatography (HPLC) method with photometric and fluorescence detection is developed for routine analysis of 2-Keto-3-deoxy-gluconate (KDG), a catabolite product of pectin and alginate. These polysaccharides are primary-based compounds for biofuel production and for generation of high-value-added products. HPLC is performed, after derivatization of the 2-oxo-acid groups of the metabolite with o-phenylenediamine (oPD), using a linear gradient of trifluoroacetic acid and acetonitrile. Quantification is accomplished with an internal standard method. The gradient is optimized to distinguish KDG from its close structural analogues such as 5-keto-4-deoxyuronate (DKI) and 2,5-diketo-3-deoxygluconate (DKII). The proposed method is simple, highly sensitive and accurate for time course analysis of pectin or alginate degradation.

Dickeya poaceiphila sp. nov., a plant-pathogenic bacterium isolated from sugar cane (Saccharum officinarum)

The genus Dickeya is an important group of plant pathogens that currently comprises 10 recognized species. Although most Dickeya isolates originated from infected cultivated plants, they are also isolated from water. The genomic sequence of the Australian strain NCPPB 569^T clearly established its separation from the previously characterized Dickeya species. The average nucleotide identity and digital DNA-DNA hybridization values obtained by comparing strain NCPPB 569^T with strains of characterized Dickeya species were lower than 87 and 32 %, respectively, supporting the delineation of a new species. The name Dickeya poaceiphila sp. nov. is proposed for this taxon with the type strain NCPPB 569^T (=CFBP 8731^T). Two other strains isolated in Australia, CFBP 1537 and CFBP 2040, also belong to this species. Phenotypic and genomic comparisons enabled the identification of traits distinguishing D. poaceiphila isolates from strains of other Dickeya species.

A novel, cupin-type phosphoglucose isomerase in Escherichia coli

BACKGROUND

Escherichia coli cells contain a homolog of presumed 5-keto-4-deoxyuronate isomerase (KduI) from pectin-degrading soil bacteria, but the catalytic activity of the E. coli protein (o-KduI) was never demonstrated.

METHODS

The known three-dimensional structure of E. coli o-KduI was compared with the available structures of sugar-converting enzymes. Based on the results of this analysis, sugar isomerization activity of recombinant o-KduI was tested against a panel of D-sugars and their derivatives.

RESULTS

The three-dimensional structure of o-KduI exhibits a close similarity with Pyrococcus furiosus cupin-type phosphoglucose isomerase. In accordance with this similarity, o-KduI was found to catalyze interconversion of glucose-6-phosphate and fructose-6-phosphate and, less efficiently, conversion of glucuronate to fructuronate. o-KduI was hexameric in crystals but represented a mixture of inactive hexamers and active dimers in solution and contained a tightly bound Zn²⁺ ion. Dilution, substrate binding and Zn²⁺ removal shifted the hexamer ⇆ dimer equilibrium to the dimers.

CONCLUSIONS

Our findings identify o-KduI as a novel phosphosugar isomerase in E. coli, whose activity may be regulated by changes in oligomeric structure.

GENERAL SIGNIFICANCE

More than 5700 protein sequences are annotated as KduI, but their enzymatic activity has not been directly demonstrated. E. coli o-KduI is the first characterized member of this group, and its enzymatic activity was found to be different from the predicted activity.

Microbial Metabolism Modulates Antibiotic Susceptibility within the Murine Gut Microbiome

Although antibiotics disturb the structure of the gut microbiota, factors that modulate these perturbations are poorly understood. Bacterial metabolism is an important regulator of susceptibility in vitro and likely plays a large role within the host. We applied a metagenomic and metatranscriptomic approach to link antibiotic-induced taxonomic and transcriptional responses within the murine microbiome. We found that antibiotics significantly alter the expression of key metabolic pathways at the whole-community and single-species levels. Notably, Bacteroides thetaiotaomicron, which blooms in response to amoxicillin, upregulated polysaccharide utilization. In vitro, we found that the sensitivity of this bacterium to amoxicillin was elevated by glucose and reduced by polysaccharides. Accordingly, we observed that dietary composition affected the abundance and expansion of B. thetaiotaomicron, as well as the extent of microbiome disruption with amoxicillin. Our work indicates that the metabolic environment of the microbiome plays a role in the response of this community to antibiotics.

Identification by Tn-seq of Dickeya dadantii genes required for survival in chicory plants

The identification of the virulence factors of plant-pathogenic bacteria has relied on the testing of individual mutants on plants, a time-consuming process. Transposon sequencing (Tn-seq) is a very powerful method for the identification of the genes required for bacterial growth in their host. We used this method in a soft-rot pathogenic bacterium to identify the genes required for the multiplication of Dickeya dadantii in chicory. About 100 genes were identified showing decreased or increased fitness in the plant. Most had no previously attributed role in plant-bacterium interactions. Following our screening, in planta competition assays confirmed that the uridine monophosphate biosynthesis pathway and the purine biosynthesis pathway were essential to the survival of D. dadantii in the plant, as the mutants ∆carA, ∆purF, ∆purL, ∆guaB and ∆pyrE were unable to survive in the plant in contrast with the wild-type (WT) bacterium. This study also demonstrated that the biosynthetic pathways of leucine, cysteine and lysine were essential for bacterial survival in the plant and that RsmC and GcpA were important in the regulation of the infection process, as the mutants ∆rsmC and ∆gcpA were hypervirulent. Finally, our study showed that D. dadantii flagellin was glycosylated and that this modification conferred fitness to the bacterium during plant infection. Assay by this method of the large collections of environmental pathogenic strains now available will allow an easy and rapid identification of new virulence factors.

Microbial hexuronate catabolism in biotechnology

The most abundant hexuronate in plant biomass is D-galacturonate. D-Galacturonate is the main constituent of pectin. Pectin-rich biomass is abundantly available as sugar beet pulp or citrus processing waste and is currently mainly used as cattle feed. Other naturally occurring hexuronates are D-glucuronate, L-guluronate, D-mannuronate and L-iduronate. D-Glucuronate is a constituent of the plant cell wall polysaccharide glucuronoxylan and of the algal polysaccharide ulvan. Ulvan also contains L-iduronate. L-Guluronate and D-mannuronate are the monomers of alginate. These raw materials have the potential to be used as raw material in biotechnology-based production of fuels or chemicals. In this communication, we will review the microbial pathways related to these hexuronates and their potential use in biotechnology.

Biochemical Reconstruction of a Metabolic Pathway from a Marine Bacterium Reveals Its Mechanism of Pectin Depolymerization

Pectin is a complex uronic acid-containing polysaccharide typically found in plant cell walls, though forms of pectin are also found in marine diatoms and seagrasses. Genetic loci that target pectin have recently been identified in two phyla of marine bacteria. These loci appear to encode a pectin saccharification pathway that is distinct from the canonical pathway typically associated with phytopathogenic terrestrial bacteria. However, very few components of the marine pectin metabolism pathway have been experimentally validated. Here, we biochemically reconstructed the pectin saccharification pathway from a marine Pseudoalteromonas sp. in vitro and show that it results in the production of galacturonate and the key metabolic intermediate 5-keto-4-deoxyuronate (DKI). We demonstrate the sequential de-esterification and depolymerization of pectin into oligosaccharides and the synergistic action of glycoside hydrolases (GHs) to fully degrade these oligosaccharides into monosaccharides. Furthermore, we show that this pathway relies on enzymes belonging to GH family 105 to carry out the equivalent chemistry afforded by an exolytic polysaccharide lyase (PL) and KdgF in the canonical pectin pathway. Finally, we synthesize our findings into a model of marine pectin degradation and compare it with the canonical pathway. Our results underline the shifting view of pectin as a solely terrestrial polysaccharide and highlight the importance of marine pectin as a carbon source for suitably adapted marine heterotrophs. This alternate pathway has the potential to be exploited in the growing field of biofuel production from plant waste.IMPORTANCE Marine polysaccharides, found in the cell walls of seaweeds and other marine macrophytes, represent a vast sink of photosynthetically fixed carbon. As such, their breakdown by marine microbes contributes significantly to global carbon cycling. Pectin is an abundant polysaccharide found in the cell walls of terrestrial plants, but it has recently been reported that some marine bacteria possess the genetic capacity to degrade it. In this study, we biochemically characterized seven key enzymes from a marine bacterium that, together, fully degrade the backbone of pectin into its constituent monosaccharides. Our findings highlight the importance of pectin as a marine carbon source available to bacteria that possess this pathway. The characterized enzymes also have the potential to be utilized in the production of biofuels from plant waste.

Probiotics in human gut microbiota can degrade host glycosaminoglycans

Glycosaminoglycans (GAGs) (e.g. heparin, chondroitin sulfate, and hyaluronan) show various significant physiological functions as a major component of extracellular matrix in animals. Some bacteria target GAGs for adhesion and/or infection to host cells, although no probiotics have been known to degrade GAGs. Here, we show GAG degradation by probiotics from human gut microbiota and their adhesion to human intestinal cells through a GAG. GAG-degrading bacteria were isolated from human faeces and identified as Enterococcus faecium, and some typical probiotics such as Lactobacillus casei, Lactobacillus rhamnosus and Enterococcus faecalis were also found to degrade heparin. GAG-degrading lactobacilli and enterococci including the isolated E. faecium possessed a genetic cluster encoding GAG-degrading/metabolising enzymes in the bacterial genome. KduI and KduD enzymes encoded in the GAG cluster of L. rhamnosus functioned as 4-deoxy-l-threo-5-hexosulose-uronate ketol-isomerase and 2-keto-3-deoxy-d-gluconate dehydrogenase, respectively, both of which were crucial for GAG metabolism. GAG-degrading L. rhamnosus and E. faecium attached to human intestinal Caco-2 cells via heparin. Some species of Bacteroides, considered to be the next generation probiotics, degraded chondroitin sulfate C and hyaluronan, and genes coding for the Bacteroides GAG-degrading enzyme were frequently detected from human gut microbiota. This is the first report on GAG-degrading probiotics in human gut microbiota.

X-Ray Solution Scattering Study of Four Escherichia coli Enzymes Involved in Stationary-Phase Metabolism

The structural analyses of four metabolic enzymes that maintain and regulate the stationary growth phase of Escherichia coli have been performed primarily drawing on the results obtained from solution small angle X-ray scattering (SAXS) and other structural techniques. The proteins are (i) class I fructose-1,6-bisphosphate aldolase (FbaB); (ii) inorganic pyrophosphatase (PPase); (iii) 5-keto-4-deoxyuronate isomerase (KduI); and (iv) glutamate decarboxylase (GadA). The enzyme FbaB, that until now had an unknown structure, is predicted to fold into a TIM-barrel motif that form globular protomers which SAXS experiments show associate into decameric assemblies. In agreement with previously reported crystal structures, PPase forms hexamers in solution that are similar to the previously reported X-ray crystal structure. Both KduI and GadA that are responsible for carbohydrate (pectin) metabolism and acid stress responses, respectively, form polydisperse mixtures consisting of different oligomeric states. Overall the SAXS experiments yield additional insights into shape and organization of these metabolic enzymes and further demonstrate the utility of hybrid methods, i.e., solution SAXS combined with X-ray crystallography, bioinformatics and predictive 3D-structural modeling, as tools to enrich structural studies. The results highlight the structural complexity that the protein components of metabolic networks may adopt which cannot be fully captured using individual structural biology techniques.

KdgF, the missing link in the microbial metabolism of uronate sugars from pectin and alginate

Uronates are charged sugars that form the basis of two abundant sources of biomass-pectin and alginate-found in the cell walls of terrestrial plants and marine algae, respectively. These polysaccharides represent an important source of carbon to those organisms with the machinery to degrade them. The microbial pathways of pectin and alginate metabolism are well studied and essentially parallel; in both cases, unsaturated monouronates are produced and processed into the key metabolite 2-keto-3-deoxygluconate (KDG). The enzymes required to catalyze each step have been identified within pectinolytic and alginolytic microbes; yet the function of a small ORF, kdgF, which cooccurs with the genes for these enzymes, is unknown. Here we show that KdgF catalyzes the conversion of pectin- and alginate-derived 4,5-unsaturated monouronates to linear ketonized forms, a step in uronate metabolism that was previously thought to occur spontaneously. Using enzyme assays, NMR, mutagenesis, and deletion of kdgF, we show that KdgF proteins from both pectinolytic and alginolytic bacteria catalyze the ketonization of unsaturated monouronates and contribute to efficient production of KDG. We also report the X-ray crystal structures of two KdgF proteins and propose a mechanism for catalysis. The discovery of the function of KdgF fills a 50-y-old gap in the knowledge of uronate metabolism. Our findings have implications not only for the understanding of an important metabolic pathway, but also the role of pectinolysis in plant-pathogen virulence and the growing interest in the use of pectin and alginate as feedstocks for biofuel production.

Structural determinants in bacterial 2-keto-3-deoxy-D-gluconate dehydrogenase KduD for dual-coenzyme specificity

Short-chain dehydrogenase/reductase (SDR) is distributed in many organisms, from bacteria to humans, and has significant roles in metabolism of carbohydrates, lipids, amino acids, and other biomolecules. An important intermediate in acidic polysaccharide metabolism is 2-keto-3-deoxy-d-gluconate (KDG). Recently, two short and long loops in Sphingomonas KDG-producing SDR enzymes (NADPH-dependent A1-R and NADH-dependent A1-R') involved in alginate metabolism were shown to be crucial for NADPH or NADH coenzyme specificity. Two SDR family enzymes-KduD from Pectobacterium carotovorum (PcaKduD) and DhuD from Streptococcus pyogenes (SpyDhuD)-prefer NADH as coenzyme, although only PcaKduD can utilize both NADPH and NADH. Both enzymes reduce 2,5-diketo-3-deoxy-d-gluconate to produce KDG. Tertiary and quaternary structures of SpyDhuD and PcaKduD and its complex with NADH were determined at high resolution (approximately 1.6 Å) by X-ray crystallography. Both PcaKduD and SpyDhuD consist of a three-layered structure, α/β/α, with a coenzyme-binding site in the Rossmann fold; similar to enzymes A1-R and A1-R', both arrange the two short and long loops close to the coenzyme-binding site. The primary structures of the two loops in PcaKduD and SpyDhuD were similar to those in A1-R' but not A1-R. Charge neutrality and moderate space at the binding site of the nucleoside ribose 2' coenzyme region were determined to be structurally crucial for dual-coenzyme specificity in PcaKduD by structural comparison of the NADH- and NADPH-specific SDR enzymes. The corresponding site in SpyDhuD was negatively charged and spatially shallow. This is the first reported study on structural determinants in SDR family KduD related to dual-coenzyme specificity. Proteins 2016; 84:934-947. © 2016 Wiley Periodicals, Inc.

Metabolic fate of unsaturated glucuronic/iduronic acids from glycosaminoglycans: molecular identification and structure determination of streptococcal isomerase and dehydrogenase

Glycosaminoglycans in mammalian extracellular matrices are degraded to their constituents, unsaturated uronic (glucuronic/iduronic) acids and amino sugars, through successive reactions of bacterial polysaccharide lyase and unsaturated glucuronyl hydrolase. Genes coding for glycosaminoglycan-acting lyase, unsaturated glucuronyl hydrolase, and the phosphotransferase system are assembled into a cluster in the genome of pathogenic bacteria, such as streptococci and clostridia. Here, we studied the streptococcal metabolic pathway of unsaturated uronic acids and the structure/function relationship of its relevant isomerase and dehydrogenase. Two proteins (gbs1892 and gbs1891) of Streptococcus agalactiae strain NEM316 were overexpressed in Escherichia coli, purified, and characterized. 4-Deoxy-l-threo-5-hexosulose-uronate (Dhu) nonenzymatically generated from unsaturated uronic acids was converted to 2-keto-3-deoxy-d-gluconate via 3-deoxy-d-glycero-2,5-hexodiulosonate through successive reactions of gbs1892 isomerase (DhuI) and gbs1891 NADH-dependent reductase/dehydrogenase (DhuD). DhuI and DhuD enzymatically corresponded to 4-deoxy-l-threo-5-hexosulose-uronate ketol-isomerase (KduI) and 2-keto-3-deoxy-d-gluconate dehydrogenase (KduD), respectively, involved in pectin metabolism, although no or low sequence identity was observed between DhuI and KduI or between DhuD and KduD, respectively. Genes for DhuI and DhuD were found to be included in the streptococcal genetic cluster, whereas KduI and KduD are encoded in clostridia. Tertiary and quaternary structures of DhuI and DhuD were determined by x-ray crystallography. Distinct from KduI β-barrels, DhuI adopts an α/β/α-barrel structure as a basic scaffold similar to that of ribose 5-phosphate isomerase. The structure of DhuD is unable to accommodate the substrate/cofactor, suggesting that conformational changes are essential to trigger enzyme catalysis. This is the first report on the bacterial metabolism of glycosaminoglycan-derived unsaturated uronic acids by isomerase and dehydrogenase.

Discovery of function in the enolase superfamily: D-mannonate and d-gluconate dehydratases in the D-mannonate dehydratase subgroup

The continued increase in the size of the protein sequence databases as a result of advances in genome sequencing technology is overwhelming the ability to perform experimental characterization of function. Consequently, functions are assigned to the vast majority of proteins via automated, homology-based methods, with the result that as many as 50% are incorrectly annotated or unannotated (Schnoes et al. PLoS Comput. Biol.2009, 5 (12), e100060520011109). This manuscript describes a study of the d-mannonate dehydratase (ManD) subgroup of the enolase superfamily (ENS) to investigate how function diverges as sequence diverges. Previously, one member of the subgroup had been experimentally characterized as ManD [dehydration of d-mannonate to 2-keto-3-deoxy-d-mannonate (equivalently, 2-keto-3-deoxy-d-gluconate)]. In this study, 42 additional members were characterized to sample sequence–function space in the ManD subgroup. These were found to differ in both catalytic efficiency and substrate specificity: (1) high efficiency (k_cat/K_M = 10³ to 10⁴ M^–1 s^–1) for dehydration of d-mannonate, (2) low efficiency (k_cat/K_M = 10¹ to 10² M^–1 s^–1) for dehydration of d-mannonate and/or D-gluconate, and 3) no-activity with either d-mannonate or d-gluconate (or any other acid sugar tested). Thus, the ManD subgroup is not isofunctional and includes d-gluconate dehydratases (GlcDs) that are divergent from the GlcDs that have been characterized in the mandelate racemase subgroup of the ENS (Lamble et al. FEBS Lett.2004, 576, 133–13615474024) (Ahmed et al. Biochem. J.2005, 390, 529–54015869466). These observations signal caution for functional assignment based on sequence homology and lay the foundation for the studies of the physiological functions of the GlcDs and the promiscuous ManDs/GlcDs.

Galactaro δ-lactone isomerase: lactone isomerization by a member of the amidohydrolase superfamily

Agrobacterium tumefaciens strain C58 can utilize d-galacturonate as a sole source of carbon via a pathway in which the first step is oxidation of d-galacturonate to d-galactaro-1,5-lactone. We have identified a novel enzyme, d-galactarolactone isomerase (GLI), that catalyzes the isomerizaton of d-galactaro-1,5-lactone to d-galactaro-1,4-lactone. GLI, a member of the functionally diverse amidohydrolase superfamily, is a homologue of LigI that catalyzes the hydrolysis of 2-pyrone-4,6-dicarboxylate in lignin degradation. The ability of GLI to catalyze lactone isomerization instead of hydrolysis can be explained by the absence of the general basic catalysis used by 2-pyrone-4,6-dicarboxylate lactonase.

Evolution of the gut microbiota and the influence of diet

Diet is a major force that shapes the composition and activity of the gut microbiota. This is evident from alterations in gut microbiota composition after weaning or drastic dietary changes. Owing to the complexity of the microbiota, interactions of intestinal bacteria with the host are difficult to study. Gnotobiotic animal models offer the opportunity to reduce the complexity and the interindividual variability of the intestinal microbiota. Germ-free animals were associated with a simplified microbial community consisting of eight bacterial species, that are found in the human gut. These microbes were selected because their genome sequences are available, and they mimic to some extent the metabolic activity of the human gut microbiota. The microbiota responded to dietary modifications by changes in the relative proportions of the community members. This model offers the chance to better define the role of intestinal bacteria in obesity development, but little is known on how diet affects intestinal bacteria at the cellular level. Mice monoassociated with Escherichia coli were used as a simplified model to investigate the influence of dietary factors on bacterial protein expression in the intestine. The mice were fed three different diets: a carbohydrate (lactose)-rich diet, a protein-rich diet and a diet rich in starch. The lactose-rich diet led to an induction of proteins involved in E. coli's oxidative stress response (Fur, AhpF, Dps). The corresponding genes are under control of the OxyR transcriptional regulator which is activated by oxidative stress. Further experiments demonstrated that osmotic stress exerted by various carbohydrates leads to an upregulation of proteins belonging to the oxyR regulon. The data suggest that the upregulated proteins enable intestinal E. coli to better cope with diet-induced osmotic stress. These examples demonstrate that gnotobiotic animal models are a valuable tool for studying diet-induced changes at the community and the cell level.

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 evolution of metabolic networks of E. coli

Background

Despite the availability of numerous complete genome sequences from E. coli strains, published genome-scale metabolic models exist only for two commensal E. coli strains. These models have proven useful for many applications, such as engineering strains for desired product formation, and we sought to explore how constructing and evaluating additional metabolic models for E. coli strains could enhance these efforts.

Results

We used the genomic information from 16 E. coli strains to generate an E. coli pangenome metabolic network by evaluating their collective 76,990 ORFs. Each of these ORFs was assigned to one of 17,647 ortholog groups including ORFs associated with reactions in the most recent metabolic model for E. coli K-12. For orthologous groups that contain an ORF already represented in the MG1655 model, the gene to protein to reaction associations represented in this model could then be easily propagated to other E. coli strain models. All remaining orthologous groups were evaluated to see if new metabolic reactions could be added to generate a pangenome-scale metabolic model (iEco1712_pan). The pangenome model included reactions from a metabolic model update for E. coli K-12 MG1655 (iEco1339_MG1655) and enabled development of five additional strain-specific genome-scale metabolic models. These additional models include a second K-12 strain (iEco1335_W3110) and four pathogenic strains (two enterohemorrhagic E. coli O157:H7 and two uropathogens). When compared to the E. coli K-12 models, the metabolic models for the enterohemorrhagic (iEco1344_EDL933 and iEco1345_Sakai) and uropathogenic strains (iEco1288_CFT073 and iEco1301_UTI89) contained numerous lineage-specific gene and reaction differences. All six E. coli models were evaluated by comparing model predictions to carbon source utilization measurements under aerobic and anaerobic conditions, and to batch growth profiles in minimal media with 0.2% (w/v) glucose. An ancestral genome-scale metabolic model based on conserved ortholog groups in all 16 E. coli genomes was also constructed, reflecting the conserved ancestral core of E. coli metabolism (iEco1053_core). Comparative analysis of all six strain-specific E. coli models revealed that some of the pathogenic E. coli strains possess reactions in their metabolic networks enabling higher biomass yields on glucose. Finally the lineage-specific metabolic traits were compared to the ancestral core model predictions to derive new insight into the evolution of metabolism within this species.

Conclusion

Our findings demonstrate that a pangenome-scale metabolic model can be used to rapidly construct additional E. coli strain-specific models, and that quantitative models of different strains of E. coli can accurately predict strain-specific phenotypes. Such pangenome and strain-specific models can be further used to engineer metabolic phenotypes of interest, such as designing new industrial E. coli strains.

The oligogalacturonate-specific porin KdgM of Erwinia chrysanthemi belongs to a new porin family

The phytopathogenic Gram-negative bacteria Erwinia chrysanthemi secretes pectinases, which are able to degrade the pectic polymers of plant cell walls, and uses the degradation products as a carbon source for growth. We characterized a major outer membrane protein, KdgM, whose synthesis is strongly induced in the presence of pectic derivatives. The corresponding gene was characterized. Analysis of transcriptional fusions showed that the kdgM expression is controlled by the general repressor of pectinolytic genes, KdgR, by the repressor of hexuronate catabolism genes, ExuR, by the pectinase gene repressor, PecS, and by catabolite repression via the cyclic AMP receptor protein (CRP) transcriptional activator. A kdgM mutant is unable to grow on oligogalacturonides longer than trimers, and its virulence is affected. Electrophysiological experiments with planar lipid bilayers showed that KdgM behaves like a voltage-dependent porin that is slightly selective for anions and that exhibits fast block in the presence of trigalacturonate. In contrast to most porins, KdgM seems to be monomeric. KdgM has no homology with currently known porins, but proteins similar to KdgM are present in several bacteria. Therefore, these proteins might constitute a new family of porin channels.

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.

The exopolygalacturonate lyase PelW and the oligogalacturonate lyase Ogl, two cytoplasmic enzymes of pectin catabolism in Erwinia chrysanthemi 3937

Erwinia chrysanthemi 3937 secretes into the external medium several pectinolytic enzymes, among which are eight isoenzymes of the endo-cleaving pectate lyases: PelA, PelB, PelC, PelD, and PelE (family 1); PelI (family 4); PelL (family 3); and PelZ (family 5). In addition, one exo-cleaving pectate lyase, PelX (family 3), has been found in the periplasm of E. chrysanthemi. The E. chrysanthemi 3937 gene kdgC has been shown to exhibit a high degree of similarity to the genes pelY of Yersinia pseudotuberculosis and pelB of Erwinia carotovora, which encode family 2 pectate lyases. However, no pectinolytic activity has been assigned to the KdgC protein. After verification of the corresponding nucleotide sequence, we cloned a longer DNA fragment and showed that this gene encodes a 553-amino-acid protein exhibiting an exo-cleaving pectate lyase activity. Thus, the kdgC gene was renamed pelW. PelW catalyzes the formation of unsaturated digalacturonates from polygalacturonate or short oligogalacturonates. PelW is located in the bacterial cytoplasm. In this compartment, PelW action could complete the degradation of pectic oligomers that was initiated by the extracellular or periplasmic pectinases and precede the action of the cytoplasmic oligogalacturonate lyase, Ogl. Both cytoplasmic pectinases, PelW and Ogl, seem to act in sequence during oligogalacturonate depolymerization, since oligomers longer than dimers are very poor substrates for Ogl but are good substrates for PelW. The estimated number of binding subsites for PelW is three, extending from subsite -2 to +1, while it is probably two for Ogl, extending from subsite -1 to +1. The activities of the two cytoplasmic lyases, PelW and Ogl, are dependent on the presence of divalent cations, since both enzymes are inhibited by EDTA. In contrast to the extracellular pectate lyases, Ca2+ is unable to restore the activity of PelW or Ogl, while several other cations, including Co2+, Mn2+, and Ni2+, can activate both cytoplasmic lyases.

Processing of the pectate lyase PelI by extracellular proteases of Erwinia chrysanthemi 3937

Erwinia chrysanthemi causes soft rot on various plants. The maceration of plant tissues is mainly due to the action of endopectate lyases. The E. chrysanthemi strain 3937 produces eight endopectate lyases (PelA, PelB, PelC, PelD, PelE, PelI, PelL and PelZ) that are secreted by the Out pathway. The necrotic response elicited by the wild-type E. chrysanthemi strain on tobacco leaves is due to an extracellular protein secreted by the Out machinery. Purification of the active factor revealed that it corresponds to a pectate lyase presenting immunological cross-reaction with PelI. Analysis of pelI and out mutants indicated that the necrosis-inducing pectate lyase results from a post-translational modification of PelI occurring extracellularly both in culture media and in planta. This modification consists of the cleavage of 97 N-terminal amino acids by the extracellular proteases of E. chrysanthemi. The enzymatic properties of the maturated form, PelI-3, are not, or only weakly, modified. However, this maturation gives rise to a small size and basic form that is active as a defence elicitor in plants.

Pectate lyase PelI of Erwinia chrysanthemi 3937 belongs to a new family

Erwinia chrysanthemi 3937 secretes five major isoenzymes of pectate lyases encoded by the pel4, pelB, pelC, pelD, and pelE genes and a set of secondary pectate lyases, two of which, pelL and pelZ, have been already identified. We cloned the pelI gene, encoding a ninth pectate lyase of E. chrysanthemi 3937. The pelI reading frame is 1,035 bases long, corresponding to a protein of 344 amino acids including a typical amino-terminal signal sequence of 19 amino acids. The purified mature PelI protein has an isoelectric point of about 9 and an apparent molecular mass of 34 kDa. PelI has a preference for partially methyl esterified pectin and presents an endo-cleaving activity with an alkaline pH optimum and an absolute requirement for Ca2+ ions. PelI is an extracellular protein secreted by the Out secretory pathway of E. chrysanthemi. The PelI protein is very active in the maceration of plant tissues. A pelI mutant displayed reduced pathogenicity on chicory leaves, but its virulence did not appear to be affected on potato tubers or Saintpaulia ionantha plants. The pelI gene constitutes an independent transcriptional unit. As shown for the other pel genes, the transcription of pelI is dependent on various environmental conditions. It is induced by pectic catabolic products and affected by growth phase, oxygen limitation, temperature, nitrogen starvation, and catabolite repression. Regulation of pelI expression appeared to be dependent on the three repressors of pectinase synthesis, KdgR, PecS, and PecT, and on the global activator of sugar catabolism, cyclic AMP receptor protein. A functional KdgR binding site was identified close to the putative pelI promoter. Analysis of the amino acid sequence of PelI revealed high homology with a pectate lyase from Erwinia carotovora subsp. carotovora (65% identity) and low homology with pectate lyases of the phytopathogenic fungus Nectria haematococca (Fusarium solani). This finding indicates that PelI belongs to pectate lyase class III. Using immunoblotting experiments, we detected PelI homologs in various strains of E. chrysanthemi and E. carotovora subsp. carotovora but not in E. carotovora subsp. atroseptica.

Characterization of the pecT control region from Erwinia chrysanthemi 3937

Erwinia chrysanthemi synthesizes and secretes pectate lyases that attack components of the plant cell wall and, therefore, play a major role in the pathogenesis of soft rot disease. We isolated a new mutant (designated pec-1), by Tn5 mutagenesis, that displays weak pectate lyase production and decreased motility and mucoidicity. Maceration and pathogenicity tests done on different plant organs showed that the pec-1 strain displays a reduced virulence compared to that of the parental strain. The Tn5 insertion was localized between the pelL and the out loci and defines a new regulatory region. Sequencing of the pec-1::Tn5 insertion revealed that pec-1 is tightly linked to the pecT regulatory gene that also controls pectate lyase synthesis. Moreover, the pecT mutation is dominant over the pec-1 mutation, suggesting that these two loci are involved in the same regulatory network. We demonstrated, by Northern blot analysis, that the pec-1::Tn5 insertion provokes derepression of pecT transcription and defines a cis-acting element. Introduction of the pecT gene in trans of a pecT::uidA fusion induced a decrease of pecT::uidA transcription, indicating a negative autoregulation. Band shift experiments confirmed that the PecT repressor specifically interacts with the pecT regulatory region. We also demonstrated that the PecT protein interacts with the regulatory region of the pelD gene encoding a pectate lyase. Therefore, the abolition of the pecT autoregulation in the pec-1 mutant provokes an overproduction of the PecT repressor that is responsible for the decrease of pectate lyase synthesis. Mutagenesis of the pecT regulatory region revealed the presence of two sites in which insertions reproduced the pec-1 phenotype. This result suggests that pecT autoregulation requires the presence of two functional operator sites. From this study, we propose that the PecT repressor binds to these two sites, generating a loop that blocks pecT transcription.

The cyclic AMP receptor protein is the main activator of pectinolysis genes in Erwinia chrysanthemi

The main virulence factors of the phytopathogenic bacterium Erwinia chrysanthemi are pectinases that cleave pectin, a major constituent of the plant cell wall. Although physiological studies suggested that pectinase production in Erwinia species is subjected to catabolite repression, the direct implication of the cyclic AMP receptor protein (CRP) in this regulation has never been demonstrated. To investigate the role of CRP in pectin catabolism, we cloned the E. chrysanthemi crp gene by complementation of an Escherichia coli crp mutation and then constructed E. chrysanthemi crp mutants by reverse genetics. The carbohydrate fermentation phenotype of the E. chrysanthemi crp mutants is similar to that of an E. coli crp mutant. Furthermore, these mutants are unable to grow on pectin or polygalacturonate as the sole carbon source. Analysis of the nucleotide sequence of the E. chrysanthemi crp gene revealed the presence of a 630-bp open reading frame (ORF) that codes for a protein highly similar to the CRP of E. coli. Using a crp::uidA transcriptional fusion, we demonstrated that the E. chrysanthemi CRP represses its own expression, probably via a mechanism similar to that described for the E. coli crp gene. Moreover, in the E. chrysanthemi crp mutants, expression of pectinase genes (pemA, pelB, pelC, pelD, and pelE) and of genes of the intracellular part of the pectin degradation pathway (ogl, kduI, and kdgT), which are important for inducer formation and transport, is dramatically reduced in induced conditions. In contrast, expression of pelA, which encodes a pectate lyase important for E. chrysanthemi pathogenicity, seems to be negatively regulated by CRP. The E. chrysanthemi crp mutants have greatly decreased maceration capacity in potato tubers, chicory leaves, and celery petioles as well as highly diminished virulence on saintpaulia plants. These findings demonstrate that CRP plays a crucial role in expression of the pectinolysis genes and in the pathogenicity of E. chrysanthemi.

Regulation of pelZ, a gene of the pelB-pelC cluster encoding a new pectate lyase of Erwinia chrysanthemi 3937

The phytopathogenic enterobacterium Erwinia chrysanthemi 3937 produces five major and several secondary endo-pectate lyases encoded by the pel genes. Most of these genes are arranged in clusters on the bacterial chromosome. The genomic region surrounding the pelB-pelC cluster was supposed to be involved in the regulation of PelB and PelC synthesis. We demonstrated that the variation of pelB expression resulted from the titration of a regulatory protein by the gene adjacent to pelC. This gene was renamed pelZ since it encodes a protein of 420 amino acids with an endo-pectate lyase activity. Regulation of pelZ expression was investigated by using transcriptional fusions and a study of mRNA synthesis. Its transcription depends on different environmental conditions. It is induced in planta and in the presence of pectic catabolite products. This induction seems to be partially mediated by the KdgR protein but does not result from a direct interaction of KdgR with the pelZ 5' region. The transcription of pelZ leads to the synthesis of a monocistronic mRNA. However, the synthesis of a polycistronic mRNA from the pelC promoter, regulated by KdgR, is responsible for increased production of PelZ under inducing conditions. pelZ transcription is also controlled by pecT, which regulates some other pel genes, but it is independent of the pecS regulatory locus. The pelZ gene appears to be widespread in different strains of E. chrysanthemi. Moreover, a gene homologous to pelZ exists in Erwinia carotovora subsp. atroseptica adjacent to the cluster containing the pectate lyase-encoding genes pel1, pel2, and pel3. This conservation could reflect a significant role of PelZ in the pectinolytic system of Erwiniae. We showed pelZ is not a predominant virulence factor of E. chrysanthemi but is involved in host specificity.

Regulation of pectinolysis in Erwinia chrysanthemi

Erwinia chrysanthemi is an enterobacterium that causes various plant diseases. Its pathogenicity results from the secretion of pectinolytic enzymes responsible for the disorganization of the plant cell wall. The E. chrysanthemi strain 3937 produces two pectin methylesterases, at least seven pectate lyases, a polygalacturonase, and a pectin lyase. The extracellular degradation of the pectin leads to the formation of oligogalacturonides that are catabolized through an intracellular pathway. The pectinase genes are expressed from independent cistrons, and their transcription is favored by environmental conditions such as presence of pectin and plant extracts, stationary growth phase, low temperature, oxygen or iron limitation, and so on. Moreover, transcription of the pectin lyase gene responds to DNA-damaging agents. The differential expressions of individual pectinase genes presumably reflect their role during plant infection. The regulation of pel genes requires several regulatory systems, including the KdgR repressor, which mediates the induction of all the pectinolysis genes in the presence of pectin catabolites. KdgR also controls the genes necessary for pectinase secretion and other pectin-inducible genes not yet characterized. PecS, a cytoplasmic protein homologous to other transcriptional regulators, can bind in vitro to the regulatory regions of pectinase and cellulase genes. The PecT protein, a member of the LysR family of transcriptional regulators, represses the expression of some pectinase genes and also affects other metabolic pathways of the bacteria. Other proteins involved in global regulations, such as CRP or HNS, can bind to the regulatory regions of the pectinase genes and affect their transcription.

The Erwinia chrysanthemi pecT gene regulates pectinase gene expression

A new type of Erwinia chrysanthemi mutant displaying a derepressed synthesis of pectate lyase was isolated. The gene mutated in these strains, pecT, encodes a 316-amino-acid protein with a size of 34,761 Da that belongs to the LysR family of transcriptional activators and presents 61% identity with the E. coli protein LrhA. PecT represses the expression of pectate lyase genes pelC, pelD, pelE, pelL, and kdgC, activates pelB, and has no effect on the expression of pelA or the pectin methylesterase genes pemA and pemB. PecT activiates its own expression. The mechanism by which PecT regulates pectate lyase synthesis is independent of that of the two characterized regulators of pectate lyase genes, KdgR and PecS. In contrast to most of the members of the LysR family, pecT is not transcribed in a direction opposite that of a gene that it regulates. pecT mutants are mucoid when grown on minimal medium plates and flocculate when grown in liquid minimal medium, unless leucine or alanine is added to the medium. Thus, pecT may regulate other functions in the bacterium.

Molecular characterization of the Erwinia chrysanthemi kdgK gene involved in pectin degradation

The pathways of pectin and galacturonate catabolism in Erwinia chrysanthemi converge to form a common intermediate, 2-keto-3-deoxygluconate (KDG), which is phosphorylated by KDG kinase encoded by the kdgK gene. We cloned the kdgK gene of E. chrysanthemi 3937 by complementing an Escherichia coli kdgK mutation, using an RP4-derivative plasmid. One of the kdgK R-prime plasmids harbored a DNA insert of about 80 kb and carried the uxuA and uxuB genes involved in glucuronate catabolism and the celY gene coding for an E. chrysanthemi cellulase. The kdgK and celY genes were precisely located on this plasmid, and their respective transcriptional directions were determined. The nucleotide sequence of the kdgK region indicated that the kdgK reading frame is 981 bases long, corresponding to a protein of 329 amino acids with a molecular mass of 36,377 Da. Analysis of the deduced primary amino acid sequence showed that this enzyme is a new member of the PfkB family of carbohydrate kinases. Expression of kdgK is controlled by a negative regulatory gene, kdgR, which represses all the steps of pectin degradation. Near the putative promoter of the kdgK gene, we identified a putative KdgR-binding site and demonstrated that the KdgR protein specifically binds in vitro to this DNA region. The KdgR-KDG couple directly mediates the phenomenon of repression or induction. The KDG kinase, by limiting the intracellular inducer concentration, appears to be a key enzyme in induction of the whole catabolic pathway.

Characterization and overexpression of the pem gene encoding pectin methylesterase of Erwinia chrysanthemi strain 3937

The pem gene encoding the pectin methylesterase (PME) of Erwinia chrysanthemi strain 3937 was subcloned and its nucleotide sequence determined. The gene contains an open reading frame of 1098 bp and codes for a protein of 366 amino acids (aa). The mature 37-kDa form of the protein is 342 aa long and has a calculated isoelectric point of 9.64. A plasmid was constructed to overproduce PME: a DNA fragment carrying pem was amplified by the polymerase chain reaction and cloned downstream from the pL promoter of the lambda phage, in a high-copy-number plasmid. In an Escherichia coli strain transformed with this plasmid, an increase in PME production of more than 60-fold was obtained, compared with the wild-type Er. chrysanthemi strain. PME represents about 5% of the total protein content of the cells. Comparison of this PME sequence with six PMEs from prokaryotic or eukaryotic organisms showed six highly conserved segments whose possible role in enzyme activity are discussed.

Some of the out genes involved in the secretion of pectate lyases in Erwinia chrysanthemi are regulated by kdgR

The out genes of Erwinia chrysanthemi are required for the translocation across the outer membrane of pectate lyases and cellulases. We present the characterization and the nucleotide sequence of five genes of the out cluster. The products of outS, B, C, D and E have significant homology with the PulS, B, C, D and E proteins necessary to the secretion of pullulanase in Klebsiella pneumoniae. An open reading frame, outT, located between outB and outC has no homology with the pul cluster but is involved in secretion. outC, outD and outE form an operon while outS, outB and outT constitute independent transcription units. outT and the outCDE operon are regulated by kdgR, the negative regulatory gene controlling pectinase production. outB and outS seem to be expressed constitutively.

Analysis of the regulation of the pelBC genes in Erwinia chrysanthemi 3937

Erwinia chrysanthemi secretes five major isoenzymes of pectate lyases encoded by the pelABCDE genes. The nucleotide sequence of the region surrounding the pelB gene of E. chrysanthemi 3937 was determined, including the regulatory regions involved in pelB and pelC expression. Analysis of the transcripts showed that transcription of pelB or pelC gave, in both cases, only one transcript. The transcription initiation sites of both pelB and pelC were precisely determined as well as the position of the transcription termination of pelB. The pelB and pelC promoters are very similar, showing a good homology with the -35 consensus region but low homology with the -10 consensus. In both cases a KdgR-box overlaps the -35 region. The pelC gene may have two KdgR operators. Moreover, the pelB and pelC genes are preceded by other sequences presenting the typical symmetry of operator sites that could be involved in more specific regulations. Comparison of E. chyrsanthemi pel regulatory regions revealed three classes of homology: pelA, pelB-pelC and pelD-pelE. The sole regulatory sequence conserved among the three classes corresponds to the KdgR-binding site. Moreover, all the pel regulatory regions are AT-rich in contrast to the coding regions which are GC-rich. Gel retardation experiments with fragments overlapping the pelB or pelC regulatory regions demonstrated that the KdgR protein specifically binds to these regions. Other proteins probably also interact with these DNA fragments. Transcription of pelB terminates in a region corresponding to a GC-rich inverted repeat followed by a run of T residues, typical of rho-independent transcription termination sites. Moreover, preliminary results imply that a region adjacent to pelC provoke, directly or indirectly, the repression of pelB and pelC expression.

Purification and functional characterization of the KdgR protein, a major repressor of pectinolysis genes of Erwinia chrysanthemi

The phytopathogenicity of the enterobacterium Erwinia chrysanthemi chiefly results from its capacity to degrade pectin, which is the major component of plant cell walls. This degradation requires the product of 12 genes which constitute independent transcriptional units. All these genes, including kdgT which encodes the 2-keto-3-deoxygluconate (KDG) transport system, are negatively regulated by the KdgR protein. The E. chrysanthemi kdgR gene was cloned into an expression vector and overexpressed in Escherichia coli. KdgR was then purified to homogeneity by two chromatographic steps as a dimer of approximately 62 kDa. Using gel retardation assays, we demonstrated that this purified repressor binds to the 25bp oligonucleotide (AAAAAAGAAACATTGTTTCATTTGT) present in the kdgT regulatory region. Dimethyl sulphate interference experiments revealed that the repressor interacts with four guanine bases and 10 adenine bases in the two strands of this KdgR box. KDG, an actual inducer of pectinolysis, releases the repressor from the operator complexes, whereas galacturonate, which is the precursor of the actual inducer, does not. These results suggest the existence of a specific interaction between KDG and KdgR protein. This study opens discussion on the relative affinity of the KdgR protein for the different operators of pectinolysis genes which are interpreted in terms of differential regulation.

Characterization of kdgR, a gene of Erwinia chrysanthemi that regulates pectin degradation

Erwinia chrysanthemi is a phytopathogenic enterobacterium able to degrade the pectic fraction of plant cell walls. The kdgR negative regulatory gene controls all the genes involved in pectin catabolism, including the pel genes encoding pectate lyases. The E. chrysanthemi kdgR regulatory gene was subcloned in Escherichia coli where it was shown to be functional, since it repressed the expression of a pelE::uidA fusion. The nucleotide sequence of kdgR contained an open reading frame of 918bp preceded by classical transcriptional initiation signals. KdgR shows similarity to two other regulatory proteins, namely GylR, encoding an activator protein of the glycerol operon in Streptomyces coelicolor, and IclR, encoding a repressor of the acetate operon in Salmonella typhimurium and in Escherichia coli. Previously, comparison of regulatory regions of several genes controlled by kdgR revealed the existence of a conserved region which was proposed as a KdgR-binding site. The 25 bp oligonucleotide AAAAAAGAAACATTGTTTCATTTGT corresponding to this consensus was substituted to the lac operator, at the beginning of transcription of the lacZ gene. This construct functioned as an operator for binding of the KdgR protein in vivo.