Virulence of Pseudomonas syringae pv. tomato DC3000 Is Influenced by the Catabolite Repression Control Protein Crc

Chemosensory systems interact to shape relevant traits for bacterial plant pathogenesis

Chemosensory systems allow bacteria to respond and adapt to environmental conditions. Many bacteria contain more than one chemosensory system, but knowledge of their specific roles in regulating different functions remains scarce. Here, we address this issue by analyzing the function of the F6, F8, and alternative (non-motility) cellular functions (ACF) chemosensory systems of the model plant pathogen Pseudomonas syringae pv. tomato. In this work, we assign PsPto chemoreceptors to each chemosensory system, and we visualize for the first time the F6 and F8 chemosensory systems of PsPto using cryo-electron tomography. We confirm that chemotaxis and swimming motility are controlled by the F6 system, and we demonstrate how different components from the F8 and ACF systems also modulate swimming motility. We also determine how the kinase and response regulators from the F6 and F8 chemosensory systems do not work together in the regulation of biofilm, whereas both components from the ACF system contribute together to regulate these traits. Furthermore, we show how the F6, F8, and ACF kinases interact with the ACF response regulator WspR, supporting crosstalk among chemosensory systems. Finally, we reveal how all chemosensory systems play a role in regulating virulence.

IMPORTANCE

Chemoperception through chemosensory systems is an essential feature for bacterial survival, as it allows bacterial interaction with its surrounding environment. In the case of plant pathogens, it is especially relevant to enter the host and achieve full virulence. Multiple chemosensory systems allow bacteria to display a wider plasticity in their response to external signals. Here, we perform a deep characterization of the F6, F8, and alternative (non-motility) cellular functions chemosensory systems in the model plant pathogen Pseudomonas syringae pv. tomato DC3000. These chemosensory systems regulate key virulence-related traits, like motility and biofilm formation. Furthermore, we unveil an unexpected crosstalk among these chemosensory systems at the level of the interaction between kinases and response regulators. This work shows novel results that contribute to the knowledge of chemosensory systems and their role in functions alternative to chemotaxis.

bacLIFE: a user-friendly computational workflow for genome analysis and prediction of lifestyle-associated genes in bacteria

Bacteria have an extensive adaptive ability to live in close association with eukaryotic hosts, exhibiting detrimental, neutral or beneficial effects on host growth and health. However, the genes involved in niche adaptation are mostly unknown and their functions poorly characterized. Here, we present bacLIFE ( https://github.com/Carrion-lab/bacLIFE ) a streamlined computational workflow for genome annotation, large-scale comparative genomics, and prediction of lifestyle-associated genes (LAGs). As a proof of concept, we analyzed 16,846 genomes from the Burkholderia/Paraburkholderia and Pseudomonas genera, which led to the identification of hundreds of genes potentially associated with a plant pathogenic lifestyle. Site-directed mutagenesis of 14 of these predicted LAGs of unknown function, followed by plant bioassays, showed that 6 predicted LAGs are indeed involved in the phytopathogenic lifestyle of Burkholderia plantarii and Pseudomonas syringae pv. phaseolicola. These 6 LAGs encompassed a glycosyltransferase, extracellular binding proteins, homoserine dehydrogenases and hypothetical proteins. Collectively, our results highlight bacLIFE as an effective computational tool for prediction of LAGs and the generation of hypotheses for a better understanding of bacteria-host interactions.

Plant myo-inositol transport influences bacterial colonization phenotypes

Plant microbiomes are assembled and modified through a complex milieu of biotic and abiotic factors. Despite dynamic and fluctuating contributing variables, specific host metabolites are consistently identified as important mediators of microbial interactions. We combine information from a large-scale metatranscriptomic dataset from natural poplar trees and experimental genetic manipulation assays in seedlings of the model plant Arabidopsis thaliana to converge on a conserved role for transport of the plant metabolite myo-inositol in mediating host-microbe interactions. While microbial catabolism of this compound has been linked to increased host colonization, we identify bacterial phenotypes that occur in both catabolism-dependent and -independent manners, suggesting that myo-inositol may additionally serve as a eukaryotic-derived signaling molecule to modulate microbial activities. Our data suggest host control of this compound and resulting microbial behavior are important mechanisms at play surrounding the host metabolite myo-inositol.

Regulation of hierarchical carbon substrate utilization, nitrogen fixation, and root colonization by the Hfq/Crc/CrcZY genes in Pseudomonas stutzeri

Bacteria of the genus Pseudomonas consume preferred carbon substrates in nearly reverse order to that of enterobacteria, and this process is controlled by RNA-binding translational repressors and regulatory ncRNA antagonists. However, their roles in microbe-plant interactions and the underlying mechanisms remain uncertain. Here we show that root-associated diazotrophic Pseudomonas stutzeri A1501 preferentially catabolizes succinate, followed by the less favorable substrate citrate, and ultimately glucose. Furthermore, the Hfq/Crc/CrcZY regulatory system orchestrates this preference and contributes to optimal nitrogenase activity and efficient root colonization. Hfq has a central role in this regulatory network through different mechanisms of action, including repressing the translation of substrate-specific catabolic genes, activating the nitrogenase gene nifH posttranscriptionally, and exerting a positive effect on the transcription of an exopolysaccharide gene cluster. Our results illustrate an Hfq-mediated mechanism linking carbon metabolism to nitrogen fixation and root colonization, which may confer rhizobacteria competitive advantages in rhizosphere environments.

Implications of carbon catabolite repression for plant-microbe interactions

Carbon catabolite repression (CCR) plays a key role in many physiological and adaptive responses in a broad range of microorganisms that are commonly associated with eukaryotic hosts. When a mixture of different carbon sources is available, CCR, a global regulatory mechanism, inhibits the expression and activity of cellular processes associated with utilization of secondary carbon sources in the presence of the preferred carbon source. CCR is known to be executed by completely different mechanisms in different bacteria, yeast, and fungi. In addition to regulating catabolic genes, CCR also appears to play a key role in the expression of genes involved in plant-microbe interactions. Here, we present a detailed overview of CCR mechanisms in various bacteria. We highlight the role of CCR in beneficial as well as deleterious plant-microbe interactions based on the available literature. In addition, we explore the global distribution of known regulatory mechanisms within bacterial genomes retrieved from public repositories and within metatranscriptomes obtained from different plant rhizospheres. By integrating the available literature and performing targeted meta-analyses, we argue that CCR-regulated substrate use preferences of microorganisms should be considered an important trait involved in prevailing plant-microbe interactions.

Influence of glucose on swarming and quorum sensing of Dickeya solani

Dickeya solani is a pathogen most frequently responsible for infecting potato plants in Europe. As in the case of most plant pathogens, its ability to colonize and invade the host depends on chemotaxis and motility. The coordinated movement of Dickeya over solid surfaces is governed by a quorum sensing mechanism. In D. solani motility is regulated by ExpI-ExpR proteins, homologous to luxI-luxR system from Vibrio fisheri, in which N-acyl-homoserine lactones (AHLs) serve as signaling molecules. Moreover, in many Gram-negative bacteria motility is coupled with central metabolism via carbon catabolite repression. This enables them to reach more nutrient-efficient niches. The aim of this study was to analyze the swarming motility of D. solani depending on the volume of the medium in the cultivation plate and glucose content. We show that the ability of this bacterium to move is strictly dependent on both these factors. Moreover, we analyze the production of AHLs and show that the quorum sensing mechanism in D. solani is also influenced by the availability of glucose in the medium and that the distribution of these signaling molecules are different depending on the volume of the medium in the plate.

The Regulatory Hierarchy Following Signal Integration by the CbrAB Two-Component System: Diversity of Responses and Functions

CbrAB is a two-component system, unique to bacteria of the family Pseudomonaceae, capable of integrating signals and involved in a multitude of physiological processes that allow bacterial adaptation to a wide variety of varying environmental conditions. This regulatory system provides a great metabolic versatility that results in excellent adaptability and metabolic optimization. The two-component system (TCS) CbrA-CbrB is on top of a hierarchical regulatory cascade and interacts with other regulatory systems at different levels, resulting in a robust output. Among the regulatory systems found at the same or lower levels of CbrAB are the NtrBC nitrogen availability adaptation system, the Crc/Hfq carbon catabolite repression cascade in Pseudomonas, or interactions with the GacSA TCS or alternative sigma ECF factor, such as SigX. The interplay between regulatory mechanisms controls a number of physiological processes that intervene in important aspects of bacterial adaptation and survival. These include the hierarchy in the use of carbon sources, virulence or resistance to antibiotics, stress response or definition of the bacterial lifestyle. The multiple actions of the CbrAB TCS result in an important competitive advantage.

Identification of IAA-regulated genes in Pseudomonas syringae pv. tomato strain DC3000

The auxin indole-3-acetic acid (IAA) is a plant hormone that not only regulates plant growth and development but also plays important roles in plant-microbe interactions. We previously reported that IAA alters expression of several virulence-related genes in the plant pathogen Pseudomonas syringae pv. tomato strain DC3000 (PtoDC3000). To learn more about the impact of IAA on regulation of PtoDC3000 gene expression we performed a global transcriptomic analysis of bacteria grown in culture, in the presence or absence of exogenous IAA. We observed that IAA repressed expression of genes involved in the Type III secretion (T3S) system and motility and promoted expression of several known and putative transcriptional regulators. Several of these regulators are orthologs of factors known to regulate stress responses and accordingly expression of several stress response-related genes was also upregulated by IAA. Similar trends in expression for several genes were also observed by RT-qPCR. Using an Arabidopsis thaliana auxin receptor mutant that accumulates elevated auxin, we found that many of the P. syringae genes regulated by IAA in vitro were also regulated by auxin in planta. Collectively the data indicate that IAA modulates many aspects of PtoDC3000 biology, presumably to promote both virulence and survival under stressful conditions, including those encountered in or on plant leaves. IMPORTANCE Indole-3-acetic acid (IAA), a form of the plant hormone auxin, is used by many plant-associated bacteria as a cue to sense the plant environment. Previously, we showed that IAA can promote disease in interactions between the plant pathogen Pseudomonas syringae strain PtoDC000 and one of its hosts, Arabidopsis thaliana. However, the mechanisms by which IAA impacts the biology of PtoDC3000 and promotes disease are not well understood. Here we demonstrate that IAA is a signal molecule that regulates gene expression in PtoDC3000. The presence of exogenous IAA affects expression of over 700 genes in the bacteria, including genes involved in Type III secretion and genes involved in stress response. This work offers insight into the roles of auxin promoting pathogenesis.

Bacterial-induced pH shifts link individual cell physiology to macroscale collective behavior

Bacteria have evolved a diverse array of signaling pathways that enable them to quickly respond to environmental changes. Understanding how these pathways reflect environmental conditions and produce an orchestrated response is an ongoing challenge. Herein, we present a role for collective modifications of environmental pH carried out by microbial colonies living on a surface. We show that by collectively adjusting the local pH value, Paenibacillus spp., specifically, regulate their swarming motility. Moreover, we show that such pH-dependent regulation can converge with the carbon repression pathway to down-regulate flagellin expression and inhibit swarming in the presence of glucose. Interestingly, our results demonstrate that the observed glucose-dependent swarming repression is not mediated by the glucose molecule per se, as commonly thought to occur in carbon repression pathways, but rather is governed by a decrease in pH due to glucose metabolism. In fact, modification of the environmental pH by neighboring bacterial species could override this glucose-dependent repression and induce swarming of Paenibacillus spp. away from a glucose-rich area. Our results suggest that bacteria can use local pH modulations to reflect nutrient availability and link individual bacterial physiology to macroscale collective behavior.

The application of the CRISPR–Cas9 system in Pseudomonas syringae pv. actinidiae

Introduction.Pseudomonas syringaepv. actinidiae (Psa) has emerged as a major bacterial pathogen of kiwifruit cultivation throughout the world.

Aim.We aim to introduce a CRISPR–Cas9 system, a commonly used genome editing tool, into Psa. The protocols may also be useful in otherPseudomonasspecies.

Methodology.Using standard molecular biology techniques, we modified plasmid pCas9, which carries the CRISPR–Cas9 sequences fromStreptococcus pyogenes,for use in Psa. The final plasmid, pJH1, was produced in a series of steps and is maintained with selection in bothEscherichia coliand Psa.

Results.We have constructed plasmids carrying a CRISPR–Cas9 system based on that ofS. pyogenes, which can be maintained, under selection, in Psa. We have shown that the gene targeting capacity of the CRISPR–Cas9 system is active and that the Cas9 protein is able to cleave the targeted sites. The Cas9 was directed to several different sites in theP. syringaegenome. Using Cas9 we have generated Psa transformants that no longer carry the native plasmid present in Psa, and other transformants that lack the integrative, conjugative element, Pac_ICE1. Targeting of a specific gene, a chromosomal non-ribosomal peptide synthetase, led to gene knockouts with the transformants having deletions encompassing the target site.

Conclusion.We have constructed shuttle plasmids carrying a CRISPR–Cas9 system that are maintained in bothE. coliandP. syringaepv. actinidiae. We have used this gene editing system to eliminate features of the accessory genome (plasmids or ICEs) from Psa and to target a single chromosomal gene.

Full Transcriptomic Response of Pseudomonas aeruginosa to an Inulin-Derived Fructooligosaccharide

Pseudomonas aeruginosa is an ubiquitous gram-negative opportunistic human pathogen which is not considered part of the human commensal gut microbiota. However, depletion of the intestinal microbiota (Dysbiosis) following antibiotic treatment facilitates the colonization of the intestinal tract by Multidrug-Resistant P. aeruginosa. One possible strategy is based on the use of functional foods with prebiotic activity. The bifidogenic effect of the prebiotic inulin and its hydrolyzed form (fructooligosaccharide: FOS) is well established since they promote the growth of specific beneficial (probiotic) gut bacteria such as bifidobacteria. Previous studies of the opportunistic nosocomial pathogen Pseudomonas aeruginosa PAO1 have shown that inulin and to a greater extent FOS reduce growth and biofilm formation, which was found to be due to a decrease in motility and exotoxin secretion. However, the transcriptional basis for these phenotypic alterations remains unclear. To address this question we conducted RNA-sequence analysis. Changes in the transcript level induced by inulin and FOS were similar, but a set of transcript levels were increased in response to inulin and reduced in the presence of FOS. In the presence of inulin or FOS, 260 and 217 transcript levels, respectively, were altered compared to the control to which no polysaccharide was added. Importantly, changes in transcript levels of 57 and 83 genes were found to be specific for either inulin or FOS, respectively, indicating that both compounds trigger different changes. Gene pathway analyses of differentially expressed genes (DEG) revealed a specific FOS-mediated reduction in transcript levels of genes that participate in several canonical pathways involved in metabolism and growth, motility, biofilm formation, β-lactamase resistance, and in the modulation of type III and VI secretion systems; results that have been partially verified by real time quantitative PCR measurements. Moreover, we have identified a genomic island formed by a cluster of 15 genes, encoding uncharacterized proteins, which were repressed in the presence of FOS. The analysis of isogenic mutants has shown that genes of this genomic island encode proteins involved in growth, biofilm formation and motility. These results indicate that FOS selectively modulates bacterial pathogenicity by interfering with different signaling pathways.

Pseudomonad reverse carbon catabolite repression, interspecies metabolite exchange, and consortial division of labor

Microorganisms acquire energy and nutrients from dynamic environments, where substrates vary in both type and abundance. The regulatory system responsible for prioritizing preferred substrates is known as carbon catabolite repression (CCR). Two broad classes of CCR have been documented in the literature. The best described CCR strategy, referred to here as classic CCR (cCCR), has been experimentally and theoretically studied using model organisms such as Escherichia coli. cCCR phenotypes are often used to generalize universal strategies for fitness, sometimes incorrectly. For instance, extremely competitive microorganisms, such as Pseudomonads, which arguably have broader global distributions than E. coli, have achieved their success using metabolic strategies that are nearly opposite of cCCR. These organisms utilize a CCR strategy termed 'reverse CCR' (rCCR), because the order of preferred substrates is nearly reverse that of cCCR. rCCR phenotypes prefer organic acids over glucose, may or may not select preferred substrates to optimize growth rates, and do not allocate intracellular resources in a manner that produces an overflow metabolism. cCCR and rCCR have traditionally been interpreted from the perspective of monocultures, even though most microorganisms live in consortia. Here, we review the basic tenets of the two CCR strategies and consider these phenotypes from the perspective of resource acquisition in consortia, a scenario that surely influenced the evolution of cCCR and rCCR. For instance, cCCR and rCCR metabolism are near mirror images of each other; when considered from a consortium basis, the complementary properties of the two strategies can mitigate direct competition for energy and nutrients and instead establish cooperative division of labor.

Chemoperception of Specific Amino Acids Controls Phytopathogenicity in Pseudomonas syringae pv. tomato

There is substantive evidence that chemotaxis is a key requisite for efficient pathogenesis in plant pathogens. However, information regarding particular bacterial chemoreceptors and the specific plant signal that they sense is scarce. Our work shows that the phytopathogenic bacterium Pseudomonas syringae pv. tomato mediates not only chemotaxis but also the control of pathogenicity through the perception of the plant abundant amino acids Asp and Glu. We describe the specificity of the perception of l- and d-Asp and l-Glu by the PsPto-PscA chemoreceptor and the involvement of this perception in the regulation of pathogenicity-related traits. Moreover, a saturating concentration of d-Asp reduces bacterial virulence, and we therefore propose that ligand-mediated interference of key chemoreceptors may be an alternative strategy to control virulence.

Chemotaxis has been associated with the pathogenicity of bacteria in plants and was found to facilitate bacterial entry through stomata and wounds. However, knowledge regarding the plant signals involved in this process is scarce. We have addressed this issue using Pseudomonas syringae pv. tomato, which is a foliar pathogen that causes bacterial speck in tomato. We show that the chemoreceptor P. syringae pv. tomato PscA (PsPto-PscA) recognizes specifically and with high affinity l-Asp, l-Glu, and d-Asp. The mutation of the chemoreceptor gene largely reduced chemotaxis to these ligands but also altered cyclic di-GMP (c-di-GMP) levels, biofilm formation, and motility, pointing to cross talk between different chemosensory pathways. Furthermore, the PsPto-PscA mutant strain showed reduced virulence in tomato. Asp and Glu are the most abundant amino acids in plants and in particular in tomato apoplasts, and we hypothesize that this receptor may have evolved to specifically recognize these compounds to facilitate bacterial entry into the plant. Infection assays with the wild-type strain showed that the presence of saturating concentrations of d-Asp also reduced bacterial virulence.

Pseudomonas donghuensis HYS virulence towards Caenorhabditis elegans is regulated by the Cbr/Crc system

Pseudomonas donghuensis HYS is the type strain of a recently identified species, P. donghuensis, which has pathogenic potential with an unclear virulence mechanism. In this study, we used Caenorhabditis elegans as a host to explore the virulence mechanism of P. donghuensis HYS. Based on a correlation between P. donghuensis HYS virulence and its repellence property, we identified 68 potential virulence-related genes, among them the Cbr/Crc system, which regulates the virulence of prokaryotic microorganisms. Slow-killing assays indicated that cbrA, cbrB, or specific sRNA-encoding genes all affected P. donghuensis virulence positively, whereas crc affected it negatively. Transcriptome analyses demonstrated that the Cbr/Crc system played an important role in the pathogenesis of P. donghuensis. In addition, experiments using the worm mutant KU25 pmk-1(km25) showed a correlation between P. donghuensis HYS virulence and the PMK-1/p38 MAPK pathway in C. elegans. In conclusion, our data show that Crc plays a novel role in the Cbr/Crc system, and the P. donghuensis virulence phenotype therefore differs from that of P. aeruginosa. This process also involves C. elegans innate immunity. These findings significantly increase the available information about Cbr/Crc-based virulence mechanisms in the genus Pseudomonas.

Prevention of Surface-Associated Calcium Phosphate by the Pseudomonas syringae Two-Component System CvsSR

CvsSR is a Ca²⁺-induced two-component system (TCS) in the plant pathogen Pseudomonas syringae pv. tomato DC3000. Here, we discovered that CvsSR is induced by Fe³⁺, Zn²⁺, and Cd²⁺ However, only supplementation of Ca²⁺ to medium resulted in rugose, opaque colonies in ΔcvsS and ΔcvsR strains. This phenotype corresponded to formation of calcium phosphate precipitation on the surface of ΔcvsS and ΔcvsR colonies. CvsSR regulated swarming motility in P. syringae pv. tomato in a Ca²⁺-dependent manner, but swarming behavior was not influenced by Fe³⁺, Zn²⁺, or Cd²⁺ We hypothesized that reduced swarming displayed by ΔcvsS and ΔcvsR strains was due to precipitation of calcium phosphate on the surface of ΔcvsS and ΔcvsR cells grown on agar medium supplemented with Ca²⁺ By reducing the initial pH or adding glucose to the medium, calcium precipitation was inhibited, and swarming was restored to ΔcvsS and ΔcvsR strains, suggesting that calcium precipitation influences swarming ability. Constitutive expression of a CvsSR-regulated carbonic anhydrase and a CvsSR-regulated putative sulfate major facilitator superfamily transporter in ΔcvsS and ΔcvsR strains inhibited formation of calcium precipitates and restored the ability of ΔcvsS and ΔcvsR bacteria to swarm. Lastly, we found that glucose inhibited Ca²⁺-based induction of CvsSR. Hence, CvsSR is a key regulator that controls calcium precipitation on the surface of bacterial cells.IMPORTANCE Bacteria are capable of precipitating and dissolving minerals. We previously reported the characterization of the two-component system CvsSR in the plant-pathogenic bacterium Pseudomonas syringae CvsSR responds to the presence of calcium and is important for causing disease. Here, we show that CvsSR controls the ability of the bacterium to prevent calcium phosphate precipitation on the surface of cells. We also identified a carbonic anhydrase and transporter that modulate formation of surface-associated calcium precipitates. Furthermore, our results demonstrate that the ability of the bacterium to swarm is controlled by the formation and dissolution of calcium precipitates on the surface of cells. Our study describes new mechanisms for microbially induced mineralization and provides insights into the role of mineral deposits on bacterial physiology. The discoveries may lead to new technological and environmental applications.

Application of Data-Independent Acquisition Approach to Study the Proteome Change from Early to Later Phases of Tomato Pathogenesis Responses

Plants and pathogens are entangled in a continual arms race. Plants have evolved dynamic defence and immune mechanisms to resist infection and enhance immunity for second wave attacks from the same or different types of pathogenic species. In addition to evolutionarily and physiological changes, plant-pathogen interaction is also highly dynamic at the molecular level. Recently, an emerging quantitative mass spectrometry-based proteomics approach named data-independent acquisition (DIA), has been developed for the analysis of the proteome in a high-throughput fashion. In this study, the DIA approach was applied to quantitatively trace the change in the plant proteome from the early to the later stage of pathogenesis progression. This study revealed that at the early stage of the pathogenesis response, proteins directly related to the chaperon were regulated for the defence proteins. At the later stage, not only the defence proteins but also a set of the pathogen-associated molecular pattern-triggered immunity (PTI) and effector triggered immunity (ETI)-related proteins were highly induced. Our findings show the dynamics of the plant regulation of pathogenesis at the protein level and demonstrate the potential of using the DIA approach for tracing the dynamics of the plant proteome during pathogenesis responses.

Physiological and genetic characterization of calcium phosphate precipitation by Pseudomonas species

Microbial biomineralization is a widespread phenomenon. The ability to induce calcium precipitation around bacterial cells has been reported in several Pseudomonas species but has not been thoroughly tested. We assayed 14 Pseudomonas strains representing five different species for the ability to precipitate calcium. Calcium phosphate precipitated adjacent to the colonies of all the Pseudomonas strains tested and also precipitated on the surface of colonies for several of the Pseudomonas strains assayed. The precipitate was commonly precipitated as amorphous calcium phosphate, however seven of the 14 Pseudomonas strains tested precipitated amorphous apatite in agar adjacent to the colonies. Out of the seven Pseudomonas strains that precipitated amorphous apatite, six are plant pathogenic. The formation of amorphous apatite was commonly observed in the area of the agar where amorphous calcium phosphate had previously formed. A transposon mutagenesis screen in Pseudomonas syringae pv. tomato DC3000 revealed genes involved in general metabolism, lipopolysaccharide and cell wall biogenesis, and in regulation of virulence play a role in calcium precipitation. These results shed light on the common ability of Pseudomonas species to perform calcium precipitation and the underlying genetic regulation involved in biomineralization.

Control of primary metabolism by a virulence regulatory network promotes robustness in a plant pathogen

Robustness is a key system-level property of living organisms to maintain their functions while tolerating perturbations. We investigate here how a regulatory network controlling multiple virulence factors impacts phenotypic robustness of a bacterial plant pathogen. We reconstruct a cell-scale model of Ralstonia solanacearum connecting a genome-scale metabolic network, a virulence macromolecule network, and a virulence regulatory network, which includes 63 regulatory components. We develop in silico methods to quantify phenotypic robustness under a broad set of conditions in high-throughput simulation analyses. This approach reveals that the virulence regulatory network exerts a control of the primary metabolism to promote robustness upon infection. The virulence regulatory network plugs into the primary metabolism mainly through the control of genes likely acquired via horizontal gene transfer, which results in a functional overlay with ancestral genes. These results support the view that robustness may be a selected trait that promotes pathogenic fitness upon infection.