Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Host-Microbe Biology

Metabolic Model of the Phytophthora infestans-Tomato Interaction Reveals Metabolic Switches during Host Colonization

Sander Y. A. Rodenburg, Michael F. Seidl, Howard S. Judelson, Andrea L. Vu, Francine Govers, Dick de Ridder
John McDowell, Invited Editor, B. Gillian Turgeon, Editor
Sander Y. A. Rodenburg
aLaboratory of Phytopathology, Wageningen University, Wageningen, the Netherlands
bBioinformatics Group, Wageningen University, Wageningen, the Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sander Y. A. Rodenburg
Michael F. Seidl
aLaboratory of Phytopathology, Wageningen University, Wageningen, the Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Michael F. Seidl
Howard S. Judelson
cDepartment of Microbiology and Plant Pathology, University of California Riverside, Riverside, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Howard S. Judelson
Andrea L. Vu
cDepartment of Microbiology and Plant Pathology, University of California Riverside, Riverside, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Francine Govers
aLaboratory of Phytopathology, Wageningen University, Wageningen, the Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Francine Govers
Dick de Ridder
bBioinformatics Group, Wageningen University, Wageningen, the Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Dick de Ridder
John McDowell
Virginia Tech
Roles: Invited Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
B. Gillian Turgeon
Cornell University
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.00454-19
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

The oomycete pathogen Phytophthora infestans causes potato and tomato late blight, a disease that is a serious threat to agriculture. P. infestans is a hemibiotrophic pathogen, and during infection, it scavenges nutrients from living host cells for its own proliferation. To date, the nutrient flux from host to pathogen during infection has hardly been studied, and the interlinked metabolisms of the pathogen and host remain poorly understood. Here, we reconstructed an integrated metabolic model of P. infestans and tomato (Solanum lycopersicum) by integrating two previously published models for both species. We used this integrated model to simulate metabolic fluxes from host to pathogen and explored the topology of the model to study the dependencies of the metabolism of P. infestans on that of tomato. This showed, for example, that P. infestans, a thiamine auxotroph, depends on certain metabolic reactions of the tomato thiamine biosynthesis. We also exploited dual-transcriptome data of a time course of a full late blight infection cycle on tomato leaves and integrated the expression of metabolic enzymes in the model. This revealed profound changes in pathogen-host metabolism during infection. As infection progresses, P. infestans performs less de novo synthesis of metabolites and scavenges more metabolites from tomato. This integrated metabolic model for the P. infestans-tomato interaction provides a framework to integrate data and generate hypotheses about in planta nutrition of P. infestans throughout its infection cycle.

IMPORTANCE Late blight disease caused by the oomycete pathogen Phytophthora infestans leads to extensive yield losses in tomato and potato cultivation worldwide. To effectively control this pathogen, a thorough understanding of the mechanisms shaping the interaction with its hosts is paramount. While considerable work has focused on exploring host defense mechanisms and identifying P. infestans proteins contributing to virulence and pathogenicity, the nutritional strategies of the pathogen are mostly unresolved. Genome-scale metabolic models (GEMs) can be used to simulate metabolic fluxes and help in unravelling the complex nature of metabolism. We integrated a GEM of tomato with a GEM of P. infestans to simulate the metabolic fluxes that occur during infection. This yields insights into the nutrients that P. infestans obtains during different phases of the infection cycle and helps in generating hypotheses about nutrition in planta.

INTRODUCTION

Plants and pathogens maintain a complex relationship that generally involves the secretion of effector proteins by the pathogen to manipulate plant cell processes and the scavenging of nutrients from the host by the pathogen to support its growth and proliferation (1). While increasing knowledge has been gained on secreted effector proteins that facilitate host colonization (2), the nature of pathogen nutrition remains underexplored. A class of organisms comprising important plant and animal pathogens is the oomycetes (3). These share morphological characteristics with fungi, yet belong to the Stramenopiles, a eukaryotic lineage that besides oomycetes also includes diatoms and brown algae (4). The most well-known oomycete is Phytophthora infestans, the causal agent of late blight disease on tomato and potato, which leads to significant yield losses worldwide (5). P. infestans is challenging to control. The pathogen has a highly flexible genome facilitating rapid adaptation to control strategies, be it resistant cultivars or chemical agents (6), and its profuse sporulation causes P. infestans to spread extremely fast (7). Hence, there is a continuous quest for novel more durable control strategies.

P. infestans sporangia are aerially dispersed, and after landing on a plant surface, these can release flagellate zoospores (8). These encyst and germinate to form a germ tube with an appressorium to penetrate epidermal cells of the plant. From there, P. infestans colonizes the mesophyll; hyphae grow in the apoplast while forming intracellular feeding structures called haustoria. Both the apoplastic hyphae and haustoria provide close contact with the plant, facilitating the exchange of effectors and nutrients (9). Oomycetes are considered osmotrophs that extracellularly catabolize complex host polymers, such as proteins, sugars, and fatty acids, using an arsenal of secreted enzymes, followed by the import of degraded nutrients into the cell (10). Moreover, P. infestans is a hemibiotrophic pathogen that requires viable host cells during the initial stages of the infection cycle, the so-called biotrophic phase with minimal symptoms. Typically, after 3 to 6 days, this biotrophic phase is followed by a necrotrophic phase during which the lesion becomes necrotic and new sporangia emerge (11, 12). This implies that the physiology of the host tissue changes throughout the infection cycle, and consequently, also the nutrients available for the pathogen (13). Conceivably, P. infestans fine-tunes its metabolism to available nutrients, for example, by regulating the expression of enzyme-encoding genes through catabolite repression and/or substrate induction (14, 15).

The P. infestans metabolism is remarkably dynamic. Transcriptome-based studies revealed significant differences in the transcript abundances of enzyme-encoding genes throughout the asexual and sexual lifecycles and during plant infection (16–19). While these studies provide detailed insight into the potential dynamics of metabolic enzymes, the insight gained into the overall characteristics of the cell metabolism is very limited. Importantly, cell metabolism is not a static framework of reactions and pathways but adapts to different environments to allow for the uptake of metabolites or the production of required compounds. These dynamics are facilitated by regulation of the rates of individual reactions within the pathways. However, as many intrinsic (e.g., enzyme activity) or extrinsic (e.g., pH or temperature) factors influence these reaction rates, studying cell metabolism remains challenging.

Genome-scale metabolic models (GEMs) can be used to help understand cell metabolism. GEMs represent cell metabolism as a cell-scale network of biochemical reactions, typically distributed over several cellular compartments, that connects the uptake of nutrients to the production of biomass precursors (20, 21). Assuming metabolism to be in steady state allows for the derivation of possible rates (called metabolic fluxes) for each reaction in the network (see File S1 in the supplemental material). A GEM can be used to predict putative nutrients and essential enzymes/reactions in the cell, thereby generating hypotheses about the responses of the cell metabolism to perturbations (22). Therefore, metabolic models have great potential to aid in the development of novel control strategies against pathogens (23, 24). The metabolism of a pathogen is tightly interconnected with that of its host and can be regarded as a single system (25). Moreover, the relationship of a pathogen with its host shapes its metabolism through evolution, leading to enzyme gene loss, which makes the pathogen dependent on its host (26). An integrated host-pathogen metabolic model can yield insights into the system-wide metabolic fluxes that shape an infection. Thus far, only a few models exist that describe the joint metabolism of beneficial or pathogenic microbes and their hosts (27–33).

FILE S1

Principles of metabolic modeling and explanation of flux coupling analysis. Download TEXT S1, PDF file, 0.2 MB.
Copyright © 2019 Rodenburg et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Recently, we published the first GEM for P. infestans (34). We used this model to simulate the growth of P. infestans on minimal culture medium and predicted the essential genes and corresponding reactions that are required to convert its nutrients into biomass components (34). Here, we integrated the P. infestans GEM with a tomato GEM (35). Using these models, we developed a host-pathogen interaction metabolic model that allowed us to identify and study hallmarks of the P. infestans-tomato interaction. Exploiting gene expression data of late blight infections on tomato enabled us to dissect the changes in pathogen and host metabolism, providing novel insights into this host-pathogen interaction.

RESULTS and DISCUSSION

A metabolic model for the tomato-P. infestans pathosystem.To reconstruct an integrated tomato-P. infestans metabolic model, we exploited a recently constructed tomato GEM (35). This tomato model comprised 2,143 reactions and 3,410 genes, whereas the previously constructed P. infestans GEM comprised 2,394 reactions and 1,301 genes, suggesting that P. infestans has less genetic redundancy for metabolic enzymes than does tomato (Table 1). Our original P. infestans GEM was based on a framework reconstructed by the annotation of enzyme orthologs in the P. infestans genome and linking those to the biochemical reactions in KEGG (34, 36). Based on the literature, we added a minimal growth medium, a set of known biomass precursors, and cellular compartments, and simulated in vitro growth (34). The tomato GEM was originally used to simulate photorespiratory fluxes under different conditions and was reconstructed from LycoCyc (35, 37). We manually improved the P. infestans and the tomato GEMs (see Materials and Methods). For instance, according to recent insights into the mitochondrial localization of glycolytic enzymes (17), we curated the gene-reaction association of P. infestans in this pathway. In the tomato GEM, we curated the thiamine biosynthesis, since P. infestans is a thiamine auxotroph (8). The two GEMs were connected by so-called transport reactions, representing the nutrient flux from tomato to P. infestans. Although the transporter repertoire of P. infestans has been predicted (16), we could not attribute specific substrates to transporters. Most functional annotations lack specificity, and even between oomycete species, transporter substrate specificity can vary (38). Therefore, we chose an unbiased approach, connecting the two GEMs by the addition of a hypothetical unidirectional transport reaction for each of the 520 metabolites that were found to be shared between the cytosol compartments, each one representing the flux of a single metabolite from the tomato cytosol to the P. infestans cytosol (Fig. 1). We realized that for a pathogen to maintain homeostasis during infection, fluxes may occur in both directions (i.e., from the host to the pathogen and vice versa). However, the lack of knowledge of this process only allowed us to model host-to-pathogen nutrient fluxes. The 520 added transport reactions predominantly concern metabolites taking part in primary metabolic subsystems (from KEGG). These include amino acid biosynthesis (n = 75), carbon metabolism (n = 61), purine and pyrimidine biosynthesis (n = 51 and 37, respectively), and ABC transporters (n = 42), which indicates that at least 42 metabolites are known ABC transporter substrates. Of the 520 introduced transport reactions, 179 could carry a flux, meaning that the respective metabolites can both be produced by the tomato plant and assimilated by P. infestans. The integration of the two GEMs resulted in an integrated metabolic model of P. infestans-tomato interaction, comprising a total of 4,695 reactions, 4,303 metabolites, and 4,578 genes (Table 1).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Integrated P. infestans-tomato model. (A) Schematic illustration. Dots are metabolites, arrows are reactions, and dotted lines represent the host-pathogen transport reactions. (B) Numbers of unique and shared cytosol metabolites.

View this table:
  • View inline
  • View popup
TABLE 1

Properties of the genome-scale metabolic models of P. infestans and tomato and of the integrated metabolic model

Flux simulations pinpoint nutrients utilized by P. infestans.Given a metabolic model and a metabolic objective for the cell, optimal values for the fluxes through all reactions can be calculated to simulate the production of P. infestans biomass precursors through the import of nutrients from the tomato cytosol (39). An objective could be, for instance, to maximize the production of biomass or to produce biomass at minimal enzyme expense. To study the versatility of the model, we composed four scenarios based on different cellular objectives (Table S2). In scenario I, we calculated the minimal set of nutrients P. infestans needs to import from tomato to form all its biomass precursors. This was already possible by importing three compounds, thiamine, l-gamma-glutamyl phosphate, and cysteine. While this scenario is unlikely to occur in planta, it does illustrate that P. infestans has a comprehensive enzyme repertoire, enabling it to form its biomass precursors based on a very small pool of nutrients. In scenario II, we calculated P. infestans biomass production using the minimum number of P. infestans reactions possible, to simulate growth with minimal enzyme cost for P. infestans. In this scenario, 48 nutrients were used, including ammonium, a known growth substrate (40). However, ammonium is an unfavored nitrogen source compared to amino acids (16). Among other imported nutrients in this scenario are the tricarboxylic acid (TCA) cycle intermediates malic acid and succinic acid, for which P. infestans has an experimentally verified transporter (38). Organic acids, such as succinic and malic acids, strongly promote P. infestans growth in vitro, especially in combination with ammonium as a source of nitrogen (41). Notably, 29 of the nutrients in this scenario are direct biomass precursors for P. infestans, including 19 amino acids (aspartate being the sole exception) (Fig. S1). In scenario III, we maximized the usage of P. infestans reactions to simulate a scenario in which the pathogen maximally exploits its metabolism. This simulation predicted a pool of 29 nutrients, including the inorganic compounds nitrite and hydrogen sulfide. Interestingly, Phytophthora spp. retained the complete assimilation pathways for these compounds (16), in contrast to multiple obligate biotrophic oomycetes which lost multiple genes throughout evolution (8). In scenario IV, we combined scenarios I and II, in which we calculated the minimal nutrient uptake combined with minimal usage of P. infestans reactions to produce biomass. Since both nutrient import and assimilation into biomass require energy, these processes likely require a trade-off; hence, we anticipate that this is a more realistic scenario. Here, 38 nutrients were used, 16 of which were amino acids.

FIG S1

Overlap in numbers of imported nutrients in three flux scenarios composed based on different objective functions and the set of P. infestans biomass precursors, as follows: Biomasspinf, P. infestans biomass precursors; Min vpinf, scenario II, minimization of P. infestans flux; Max vpinf, scenario III, maximization of P. infestans flux; and Min vpinf ∪ vtransport, scenario IV, minimize both P. infestans and host-pathogen transport flux. The biomass precursors and imported nutrients in each scenario are listed in Table S2. Download FIG S1, TIF file, 0.2 MB.
Copyright © 2019 Rodenburg et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The four scenarios described above have opposite objectives (minimize/maximize the fluxes of P. infestans) and simulate extreme circumstances not likely found in nature. The true pool of imported nutrients will likely be a combination of the pools predicted in each individual scenario. A comparison of the nutrients imported in all four scenarios revealed sets of nutrients in common between scenarios. There are 14 nutrients found in three out of four scenarios, i.e., II, III, and IV (Fig. S1), including several amino acids, nucleotide precursors, and glycerol 3-phosphate as a lipid precursor. Surprisingly, aspartate is not imported in any of the scenarios, even though it is a biomass precursor of P. infestans and serves as a precursor to a variety of other amino acids (42, 43). An oomycete-specific form of aspartate aminotransferase (EC 2.6.1.1) was found that seems to play a key role in pathogenicity in Phytophthora sojae, possibly by balancing nitrogen and carbon metabolism and facilitating the interconversion between amino acids (44). The presence of this enzyme in our model explains why aspartic acid import may not be necessary. Interestingly, thiamine is only imported in scenario III (maximizing the P. infestans fluxes), while thiamine pyrophosphate (TPP) is imported in scenarios II and IV (minimizing fluxes). It is assumed that many oomycetes import thiamine to form TPP (8), an essential cofactor in carbohydrate metabolism (45) and in TPP-responsive riboswitches (46). However, related Stramenopiles were found to also grow on thiamine alternatives (47), and possibly, P. infestans can import up- or downstream compounds instead. Taken together, these simulations suggest a range of nutrients from tomato that can be effectively assimilated by P. infestans. Nutrients in common between the simulated scenarios suggest that these are likely more versatile than others and hence most useful for P. infestans to import.

Network analysis identifies dependencies of P. infestans on tomato metabolism.To identify dependencies of P. infestans on tomato, we investigated the topology of the integrated model by looking for essential reactions (i.e., reactions that are indispensable in the model to form all biomass precursors) and (inter)dependencies of reactions (i.e., the flux of a particular reaction relies on the flux of another), also known as coupled reactions. Coupled reactions within a metabolic model can be identified using a method called flux coupling analysis (FCA; File S1) (48). FCA of the model identified 77 coupled P. infestans-tomato reaction pairs (Fig. 2) involving 53 unique P. infestans reactions and 49 unique tomato reactions (Table S3). One of these P. infestans-tomato reaction couplings comprises the biomass reaction of TPP in P. infestans, which is coupled to 17 tomato reactions, illustrating that P. infestans TPP is dependent on the thiamine biosynthesis pathway in tomatoes. In addition, of the 4,695 reactions in the model, 112 P. infestans reactions and 35 tomato reactions were found to be essential for P. infestans (Fig. 2A). Some reactions are essential and coupled to many other reactions at the same time (appearing as hubs in Fig. 2A), implying that these play a central role in the model, with potentially large biological implications. For instance, there is a tomato transport reaction of aspartate into the plastid compartment, supplying aspartate as an amino group donor for the synthesis of a thiamine precursor (aminoimidazole ribotide) (49, 50), thus indirectly making this tomato reaction essential for P. infestans. Similarly, ATP and phosphate transport between P. infestans cytosol and mitochondria is essential and seems to play an important role in the model. Phosphate and sulfate uptake by the tomato plant are part of the defined growth medium, hence coupled to reactions in both species. FCA also revealed several clusters of tightly interconnected coupled essential reactions (Fig. 2A). One of these represents P. infestans fatty acid biosynthesis. In the model, this comprises a mostly linear pathway, and consequently, most reactions are coupled (interdependent). Fatty acid biosynthesis in oomycetes is associated with fungicide resistance (51) and energy storage for sporangia (8).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Flux coupling between reactions in the P. infestans-tomato model. (A) Graph showing the coupled reactions in the model. Nodes represent reactions in tomato (red) or P. infestans (blue) and host-pathogen transport (green), and edges represent coupling between those reactions. Node size reflects essentiality for P. infestans biomass production. Stars represent transport reactions listed in panel B. Highly connected nodes (1 to 5) and clusters (I to V) are indicated and listed in the boxes on the right (see also main text). (B) Nutrients associated with the 12 most frequently coupled host-pathogen transport reactions. The bars are stacked and indicate the number of coupled reactions per species.

Couplings to the host-pathogen transport reactions in particular can provide information about the importance of transported metabolites; the more P. infestans reactions that are coupled to a transport reaction, the more likely that the associated substrate is important for P. infestans. We selected the most frequently coupled transport reactions from the model and assessed their couplings to other tomato and P. infestans reactions (Fig. 2B). The substrates of these transport reactions were associated with a diverse range of metabolic processes in P. infestans, i.e., pantothenate/coenzyme A (CoA) biosynthesis (2-dehydropantoate), de novo pyrimidine biosynthesis (S-dihydroorotate) (52), riboflavin metabolism [2,5-diamino-4-hydroxy-6-(5-phosphoribosylamino)pyrimidine], and inositol phosphate metabolism (myo-inositol 3-phosphate). Notably, three transporter substrates were associated with glycerophospholipid metabolism (phosphodimethylethanolamine N−methylethanolamine phosphate and choline). It is conceivable that P. infestans can take up glycerophospholipid precursors to facilitate the formation of membrane lipids (53). Seven P. infestans reactions were coupled to the reaction importing fructose, an efficient growth substrate (40). Among the most frequently coupled host-pathogen transport reactions are also those transporting dTTP, histidine, and leucine. Since these are mainly coupled to tomato reactions, the import of these nutrients seems to depend on their biosynthesis pathways in tomato. FCA can identify coupled reactions in metabolic models that are not necessarily directly connected in the network and can thus derive functionally related modules of (pathogen-host) metabolism that may otherwise be overlooked.

Transcriptome-based submodels provide insight into the dynamics of P. infestans-tomato metabolism.To obtain insight into the metabolism at subsequent developmental stages of the P. infestans-tomato interaction, we integrated the model with dual-transcriptome data obtained from RNA, isolated from tomato leaves inoculated with P. infestans strain 1306 (2 to 6 days postinoculation, sampled every 4 h). The infections showed profuse sporulation at around 4 days postinoculation (dpi). In previous studies, it was found that in infections by this P. infestans strain, marker genes for biotrophic growth peak at 2 to 3 dpi, and those for necrotrophic growth peak at 4 to 5 dpi (16), while marker genes for sporulation are increasingly expressed from 3 dpi on and peak at 4 dpi. Expectedly, as the proportion of reads of each sample mapping to the P. infestans genome steadily increased over time (from 7% to 84%), the amount of reads mapping to the tomato genome decreased accordingly (from 87% to 8%) (Fig. S2). The transcriptome data were used to generate time point-specific submodels, which are subsets of the full model according to the expression of the genes in the model. Since gene expression and metabolic activity are not directly related, we used the Integrative Network Inference for Tissues (INIT) algorithm (54) that calculates a submodel with maximal agreement to the expression of genes in the model such that biomass can be produced (see Materials and Methods). This resulted in 25 submodels containing 32 to 44% of the reactions of the full model (Fig. 3), indicating that only a subset of reactions is required to form all defined biomass precursors. Over time, the total number of reactions per submodel decreases (Fig. 3), mostly due to a reduction in the number of P. infestans reactions. The number of tomato reactions is relatively stable across the submodels, while the number of transport reactions increases slightly over time. These results suggest that as the infection progresses, P. infestans relies less on its own metabolism but increasingly imports metabolites from the necrotic tomato lesion. Conceivably, once P. infestans switches to its necrotrophic lifestyle, nutrients can be more easily obtained from the decaying leaf tissue (15). To obtain an impression of the metabolic changes in the infection over time, we calculated all pairwise distances of the reaction content of the submodels (Fig. 3). This revealed a gradient of similarity scores, clearly reflecting metabolic changes over time, possibly as a response to changing nutrient availability in the tomato leaf tissue. A relatively large distance can be observed between groups of submodels, suggesting two more profound switches in metabolism at roughly 3 days and 8 h (3d/8h) and 4 days and 12 h (4d/12h), possibly reflecting a transition from the pathogens’ biotrophic growth to necrotrophic growth and sporulation (11, 12, 16, 55).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Submodels based on dual-transcriptome data from a time course covering a full infection cycle of P. infestans on tomato leaf. Submodels are representative for the sampling time postinoculation shown in days (d) and hours on the y axis (middle) and x axis (left). Left, stacked bar graph indicating the number of reactions per species and number of host-pathogen transport reactions (not part of either species). Right, Jaccard similarity (intersection divided by union) between the reaction content of the submodels.

FIG S2

Statistics of the RNA sequencing (RNA-seq) read mapping. (A) Per-base quality values for the Illumina reads (all samples combined) derived from the FastQC software. (B) Mapping efficiency per sample for P. infestans and tomato. Download FIG S2, TIF file, 1.3 MB.
Copyright © 2019 Rodenburg et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

To obtain more insight into the metabolic processes of each submodel, we performed Fisher’s exact tests on the KEGG pathways in the submodels. A variety of metabolic pathways were overrepresented in submodels compared to the full model (Fig. 4). Consistent with the previously observed gradient of similarity scores (Fig. 3), groups of overrepresented pathways can be distinguished before and after 3d/8h. Notably, in the early submodels prior to the first transition, P. infestans reactions are enriched for several primary metabolic pathways, such as amino acid biosynthesis, glycolysis, and the TCA cycle. This could be a response to sugars in the apoplast that are still produced by photosynthesizing tomato leaf cells during the biotrophic phase of infection (56, 57). In contrast, late submodels, and in particular after the second transition (4d/12h), show that in tomato, several amino acid biosynthesis pathways are enriched, suggesting that during early infection, amino acids are mostly synthesized de novo by P. infestans but are later scavenged from tomato. This corroborates with a transcriptome study of P. infestans infecting potatoes, in which a diverse but generally higher expression of amino acid transporter genes at late infection time points was observed (16). In addition, early submodels are enriched in folate-related pathways on both sides, suggesting that folate plays an important role in early infection. This is in line with the finding that P. infestans has an unusually large repertoire of putative folate-biopterin transporter genes (18), although a divergent expression pattern was observed for these genes throughout infection (16). Arachidonic acid metabolism is enriched in midinfection submodels, which is a characteristic fatty acid in P. infestans known to elicit plant defense and to act as a signaling molecule in plants (58, 59). In summary, our analyses reveal profound changes in metabolic processes during infection of tomato by P. infestans and suggest that P. infestans reduces its metabolism at the expense of its host.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Enriched KEGG pathways in transcriptome-based submodels from a time course covering a full infection cycle of P. infestans on tomato. Submodels are representative for the sampling time postinoculation shown on the x axis in days (d) and hours. Enrichment is calculated based on the reaction content of each of the submodels compared to the full model. Colors scale to the adjusted P value. metab, metabolism; biosyn, biosynthesis; syn, synthesis; deg, degradation.

Import of specific nutrients becomes increasingly essential to P. infestans.The ability of a model to maintain its functionality (biomass production) under perturbations (e.g., simulated reaction deletions) is often referred to as robustness (60). We can express the robustness of each submodel as the fraction of reaction deletions that do not disturb P. infestans biomass production (60). This revealed that early infection submodels have higher robustness to reaction deletions than late infection submodels (Fig. 5). This suggests that as the infection progresses, P. infestans largely shuts down accessory/alternative pathways to synthesize its biomass precursors, rendering the remaining reactions essential. The number of essential genes for P. infestans remained stable over time, with an average of 84 ± 5 essential genes (Table S1).

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Robustness of transcriptome-based submodels from a time course covering a full infection cycle of P. infestans on tomato. The y axis shows the robustness as the fraction of single reaction deletions that do not disable the biomass flux for P. infestans. Submodels are representative of the sampling time postinoculation shown on the x axis in days and hours.

TABLE S1

Reactions in the P. infestans-tomato model iSR4578, the associated pathways, formulas, and genes, and their occurrence in the 25 submodels. Reactions that were added during the reconstruction process are listed, as well as essential genes per submodel. Download Table S1, XLSX file, 1.0 MB.
Copyright © 2019 Rodenburg et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S2

Nutrients imported per scenario (objective function). Download Table S2, XLSX file, 0.01 MB.
Copyright © 2019 Rodenburg et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S3

Flux coupling analysis results. The cross-species couplings and the flux couplings to transport reactions are shown. Download Table S3, XLSX file, 0.2 MB.
Copyright © 2019 Rodenburg et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

To evaluate the importance of nutrient transport for P. infestans while colonizing tomato, we assessed the essentiality of the host-pathogen transport reactions in each submodel (61). A host-pathogen transport reaction is essential when its deletion disables biomass production and is partially essential when deletion in combination with deletions of other reactions disables biomass production. We found no essential nutrients for P. infestans when considering the full model, yet the submodels display various patterns of nutrient essentiality (Fig. S3). Early infection submodels have just a few essential nutrients, while mid- and late-infection submodels, after the first transition point, show various essential amino acid transport reactions. This corroborates our previous observations (Fig. 4) where the first transition point marks the switch from de novo synthesis of amino acids in P. infestans to an increased uptake from tomato.

FIG S3

Essentiality score (ESS) for each host-pathogen transport reaction per submodel. The x axis represents transcriptome-based submodels, the y axis indicates the nutrient that is transported. The transport reactions of nutrients indicated in black are fully essential, and those indicated in various shades of gray are partially essential. Download FIG S3, TIF file, 0.5 MB.
Copyright © 2019 Rodenburg et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Metabolomics can be used to refine transcriptome-based submodels.To assess to what extent the transcriptome-based submodels are coherent with metabolome data, we utilized untargeted metabolome data of tomato leaves colonized by P. infestans at 2 days and 12 h (2d/12h) and 5 days and 12 h (5d/12h) postinoculation. To relate the detected metabolites to metabolic fluxes, we hypothesized that metabolites strongly decreasing in abundance (log2 fold change <−2) between these two time points are produced in the submodel of 2d/12h, and metabolites that increase (log2 fold change >2) are produced in the submodel of 5d/12h (Fig. S4). Subsequently, we generated two additional submodels using INIT (54) based on the transcriptome data at 2d/12h and 5d/12h postinoculation, while enforcing the presence of the detected metabolites in the respective submodels. The submodels of the transcriptome-only (T) and the metabolome-guided (T+M) submodels differed in 4% (2d/12h) and 11% (5d/12h) of the reactions and in 3% and 7% of the metabolites, respectively (Fig. 6). Overall, the addition of metabolomics data leads to a net increase of reactions and metabolites in particular at the late submodels, suggesting that the T submodels are slightly too conservative; there may be a more active metabolism than suggested by mRNA levels of enzyme-encoding genes. It should be emphasized that both GEMs and metabolomics data only provide a partial description of the cell metabolism at the different infection time points. Discrepancies may arise since transcriptome-based submodels do not take into account any posttranscriptional regulation or posttranslational modification of enzymes, and misannotations and missing/blocked reactions potentially increase uncertainty in submodels. Moreover, since detected metabolites are not necessarily continuously metabolized, steady-state fluxes are not directly linked to measured metabolite abundances (62). Vice versa, metabolic fluxes do not necessarily yield detectable metabolites, for example, when reaction products are immediately consumed in downstream reactions. To predict metabolic fluxes more reliably, higher-resolution metabolome data and sufficient data points over time are prerequisites which could provide sufficient data points in time to calculate metabolite coefficients (63). Nonetheless, our initial results suggest that the integration of high-resolution metabolome data of tomato infection has the potential to refine stage-specific patterns that are embedded in the joint metabolism of the P. infestans-tomato interaction.

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Overlap of reaction and metabolite content of submodels, based on transcriptome data only (T) or transcriptome data and metabolome data (T+M). Data used to generate these submodels were obtained from P. infestans-infected tomato leaves at 2d/12h and 5d/12h postinoculation.

FIG S4

Integration of metabolome data in the model. (A) Volcano plot of the metabolome data, highlighting metabolites that show a significant log2 fold change of <−2 or >2 between two time points postinoculation, 2d/12h and 5d/12h. The x axis represents the log2 fold change, and the y axis represents the significance with the −log10 P value. (B) Adding metabolome data leads to an increase or decrease in reactions in various pathways. The analysis was done on two submodels, 2d/12h (left) and 5d/12h (right). Negative bars indicate reactions that were removed in the metabolome-guided (T + M) submodels compared to the transcriptome-only (T) submodels, and positive bars indicate reactions that were added. Download FIG S4, TIF file, 0.8 MB.
Copyright © 2019 Rodenburg et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Conclusions.Pathogens scavenge metabolites from their hosts to support their growth and proliferation. Here, we took a systems biology approach and modeled the joint metabolism of tomato and P. infestans to generate hypotheses about their relationship. Our metabolic model of the P. infestans-tomato interaction represents one of the few integrated pathogen-host metabolic models published to date (23). Regarding pathogen and host as one entity can yield hypotheses about the combined metabolism at a single metabolic equilibrium. The modeling allowed us to infer a conceivable pool of nutrients, mainly consisting of amino acids, lipid precursors, and a TPP precursor that is likely exploited by P. infestans while infecting the tomato. The model also helped us further characterize host-pathogen dependencies, such as the long-known thiamine dependency of Phytophthora spp. (64). Modification of the thiamine biosynthesis pathway in tomato might be an interesting strategy for controlling late blight. Similarly, P. infestans lipid metabolism might be an interesting pathway for pathogen intervention, as our model predicts that P. infestans takes up lipids as a membrane precursor, coherent with previous experiments (53). Alternatively, fatty acid biosynthesis likely has an important role in synthesizing fuel reserves for spores (8). Our analyses showed that fatty acid biosynthesis is a largely linear pathway with high interdependency of the participating reactions (34).

Clearly, many parameters that can further improve this model are still unknown, for example, the biomass composition of P. infestans growing in planta and different metabolic processes in different zones of a lesion and during the infection cycle (65). As such, the model would benefit from more extensive and in-depth metabolomics analyses (66). There is also a need to extend our knowledge on the potential role of transporters in this system, as these play a key role in infection and are potential control targets (67, 68). This could be done, for example, by 13C-flux spectral analyses to monitor nutrient fluxes and to validate transporter functions in the pathogen (38, 69). Altogether, this model provides insights into P. infestans-tomato metabolism and serves as a stepping stone for the design of novel control strategies for this devastating pathogen.

MATERIALS AND METHODS

Model reconstruction.To improve the predictive capabilities of the previously constructed P. infestans iSR1301 GEM (34), we performed several literature- and protocol-based curations (70). As a starting point, we used our published genome-scale metabolic model (GEM) of P. infestans iSR1301 (34). We removed the sink- and gap-filling reactions that were previously added (34). The sink reactions of the generic metabolites fatty acid, holo-acyl-carrier protein, and apoprotein have no mass or formula but are needed to calculate fatty acid biosynthesis fluxes and were therefore retained. Recently, it was confirmed that some metabolic enzymes of glycolysis and serine biosynthesis in P. infestans reside in mitochondria (17). We manually corrected the gene-reaction associations accordingly to assign these enzymes to the appropriate compartment. The KEGG reaction identifiers of the P. infestans GEM were matched to MetaNetX identifiers (71), and the associated reaction formula was retrieved. MetaNetX indicates for each reaction whether it is mass balanced or not; all reactions that were not mass balanced were removed from the model. The conversion of arachidonic acid to eicosapentaenoic acid (EPA) was manually added according to the literature, whereafter EPA was added to the biomass precursors (59, 72). Thiamine diphosphate was added to the set of P. infestans biomass precursors to represent the thiamine auxotrophy of P. infestans. Water and proton metabolites were removed from the model, as they do not fulfill a meaningful function in the P. infestans GEM but do have a function in the tomato GEM, where they simulate proton pumps (35). The thermodynamic constraints (directionality) were by default inferred from KEGG maps and manually adjusted according to standard model reconstruction protocols (70), such that ATP/GTP-consuming reactions were unidirectional, except for ATP synthase (EC 3.6.3.14), nucleotide diphosphate kinase (EC 2.7.4.6), and succinate coenzyme A synthase (EC 6.2.1.4).

To identify possible missing enzymes for P. infestans in the model, we annotated KEGG enzyme orthologs (KOs) in 19 oomycete proteomes which were downloaded from FungiDB (2 June 2018) (73) (Table S1), as previously described (34). We selected enzymes present in at least two oomycetes and not yet in the model. The encoding protein sequences were aligned to the P. infestans T30-4 genome sequence (74), masked for its annotated open reading frames, using tblastn (v2.2.31+). Sequences with an E value of ≤1e−50 and query coverage of >90% were submitted to the LocTree3 Web server (75) to predict their subcellular localization, and the associated reactions were added to the model. Additionally, P. infestans gene models were manually corrected and reannotated for enzyme orthologs (34). To avoid spurious addition of reactions, any of the added reactions described above were removed if they were eventually found to be unable to carry flux in the model.

The tomato GEM was inferred from supplemental files of the associated publication (35). Since the tomato GEM iHY3510 was built by a different lab using a different reconstruction methodology, we tried to keep adaptations to a minimal level. As was done for the P. infestans model identifiers, the reaction identifiers of the tomato model were matched to MetaNetX reaction identifiers for which the associated formula was retrieved. According to knowledge of thiamine biosynthesis in plants (50), the thiamine biosynthesis was manually curated, such that thiamine is formed in the tomato cytosol from a plastidial thiazole and pyrimidine precursor (76). The P. infestans and tomato GEMs were connected by unidirectional transport reactions for all metabolites shared between the tomato cytosol and the P. infestans cytosol compartments (Table S1).

Constraint-based modeling and optimization.Modeling was performed in MATLAB (R2017B) using the RAVEN (v2) (77) and COBRA (v3) toolboxes (78). Flux coupling analysis (FCA) was performed using F2C2 (v0.91) (79). Optimization was performed using Gurobi (v8.0) for flux balance analysis (FBA) and the Integrative Network Inference for Tissues (INIT) (54). GLPK (v2.8) was used for FCA. Essential genes/reactions and essentiality scores of the transport reactions were calculated using ESS (61), which implements Fast-SL (80), constraining the total flux of P. infestans biomass production to >5%. The global robustness statistic was calculated according to Peyraud et al. (60), defined as the fraction of reaction deletions that did not render the biomass flux of P. infestans <5%.

RNA sequencing and mass spectrometry.RNA and metabolites were isolated from tomato leaflets (cv. New Yorker) inoculated with P. infestans isolate 1306. For the plant infection assays, a zoospore suspension was prepared as described previously (18) and adjusted to a concentration of 5 × 104 per ml, applied to detached leaves placed on 1.5% water agar. The zoospores were sprayed on the leaves with a hand sprayer until run-off. The leaves were then incubated at 18°C under high humidity in plastic bags containing wet paper towels with a 12-h light/dark cycle. Between 2 and 6 days postinoculation, leaves were harvested every 4 h and flash-frozen until further use. After library construction, single-end 75-nucleotide (nt) sequence reads were obtained using an Illumina NextSeq 500 platform, and quality was assessed using FastQC (v0.11.8) (Fig. S2) (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). The reads were independently mapped to the P. infestans T30-4 (74) and tomato ITAG2.3 (81) genomes using HISAT2 (v2.1.0) (82), and mapping efficiencies were retrieved using SAMtools flagstat (v0.1.19) (83). Transcript abundance was quantified and normalized using cuffnorm (v2.2.1) (84) that implements the DESeq2 normalization procedure (85), which divides the read counts by a factor calculated from the median of geometric means across samples, accounting for differences in sequencing depth and RNA composition.

For the metabolome analyses, the flash-frozen leaves were lyophilized, weighed, and provided to Metabolon, Inc. for further handling. After tissue grinding, proteins were removed by methanol precipitation. After eliminating the solvent, the samples were analyzed by reverse-phase (RP) ultraperformance liquid chromatography (UPLC)-tandem mass spectrometry (MS/MS) using positive-ion mode electrospray ionization (ESI), RP-UPLC-MS/MS with negative-ion mode ESI, and hydrophilic interaction-UPLC-MS/MS with negative-ion mode ESI. The area under the curve method was used to quantify peaks. For each time point, four biological replicates were analyzed. Differential abundance of metabolites was determined using t tests, selecting metabolites with a log2 fold change of >2 or <−2 with a Benjamini-Hochberg-adjusted P value of <0.05 (86). We imputed missing values among replicates using the minimal value across other replicates.

Submodel generation.The transcriptome-based submodels were generated using the INIT algorithm (54), which poses a mixed-integer linear optimization problem (MILP) that aims to optimize a global score based on reaction weights in the model. Reaction weights were calculated according to an adapted version of the formula described by Agren et al. (54):

wij=5log⁡(Eijx¯i + 1), where wij is a weight for each reaction i and condition j in the full model, Eij is the transcript abundance value summarized per reaction (i.e., the maximum expression value of the genes associated with each reaction), and x̄i is the mean expression of that reaction across samples. An inflation factor (+1) was added in the denominator to prevent inflated weights for extremely low expression values. According to Agren et al. (54), we maintained the same minimum and maximum cutoffs for w of −5 ≤ wij ≤ 10. Reactions with an expression lower than the average (thus negative weight) will be less likely to be included in the submodels, and reactions expressed above average (thus a positive weight) will be more likely to be included. Transport reactions were assigned a weight of −0.1, based on the hypothesis that these should not be “free” to include (i.e., weight 0) but should only have a minimal influence on the submodel solution. Reactions without associated genes were given a weight of 0. For the MILP optimization, the flux of all biomass precursors of both tomato and P. infestans was constrained to be >1.

Data availability.Transcriptome data were deposited in the NCBI Sequence Read Archive under BioProject number PRJNA516028.

TABLE S4

Gene expression values summarized per reaction and the quantitative metabolome data of the tomato leaf infection. Download Table S4, XLSX file, 0.7 MB.
Copyright © 2019 Rodenburg et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

ACKNOWLEDGMENTS

This work was funded by the Food-for-Thought campaign from Wageningen University Fund, by a VENI grant (to M.F.S.) from the Netherlands Organization for Scientific Research, and by awards to H.S.J. from the National Institute of Food and Agriculture of the U.S. Department of Agriculture and the National Science Foundation.

FOOTNOTES

    • Received 18 February 2019
    • Accepted 3 June 2019
    • Published 9 July 2019
  • Copyright © 2019 Rodenburg et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Rovenich H,
    2. Boshoven JC,
    3. Thomma B
    . 2014. Filamentous pathogen effector functions: of pathogens, hosts and microbiomes. Curr Opin Plant Biol 20:96–103. doi:10.1016/j.pbi.2014.05.001.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Whisson SC,
    2. Boevink PC,
    3. Wang S,
    4. Birch P
    . 2016. The cell biology of late blight disease. Curr Opin Microbiol 34:127–135. doi:10.1016/j.mib.2016.09.002.
    OpenUrlCrossRef
  3. 3.↵
    1. Derevnina L,
    2. Petre B,
    3. Kellner R,
    4. Dagdas YF,
    5. Sarowar MN,
    6. Giannakopoulou A,
    7. de la Concepcion JC,
    8. Chaparro-Garcia A,
    9. Pennington HG,
    10. van West P,
    11. Kamoun S
    . 2016. Emerging oomycete threats to plants and animals. Philos Trans R Soc B Biol Sci 371:20150459. doi:10.1098/rstb.2015.0459.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. McCarthy CGP,
    2. Fitzpatrick DA
    . 2017. Phylogenomic reconstruction of the oomycete phylogeny derived from 37 genomes. mSphere 2:e00095-17. doi:10.1128/mSphere.00095-17.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Kamoun S,
    2. Furzer O,
    3. Jones JDG,
    4. Judelson HS,
    5. Ali GS,
    6. Dalio RJD,
    7. Roy SG,
    8. Schena L,
    9. Zambounis A,
    10. Panabières F,
    11. Cahill D,
    12. Ruocco M,
    13. Figueiredo A,
    14. Chen XR,
    15. Hulvey J,
    16. Stam R,
    17. Lamour K,
    18. Gijzen M,
    19. Tyler BM,
    20. Grünwald NJ,
    21. Mukhtar MS,
    22. Tomé DFA,
    23. Tör M,
    24. Van Den Ackerveken G,
    25. Mcdowell J,
    26. Daayf F,
    27. Fry WE,
    28. Lindqvist-Kreuze H,
    29. Meijer HJG,
    30. Petre B,
    31. Ristaino J,
    32. Yoshida K,
    33. Birch PRJ,
    34. Govers F
    . 2015. The top 10 oomycete pathogens in molecular plant pathology. Mol Plant Pathol 16:413–434. doi:10.1111/mpp.12190.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Fry WE
    . 2016. Phytophthora infestans: new tools (and old ones) lead to new understanding and precision management. Annu Rev Phytopathol 54:529–547. doi:10.1146/annurev-phyto-080615-095951.
    OpenUrlCrossRef
  7. 7.↵
    1. Leesutthiphonchai W,
    2. Vu AL,
    3. Ah-Fong AM,
    4. Judelson HS
    . 2018. How does Phytophthora infestans evade control efforts? Modern insight into the late blight disease. Phytopathology 108:916–924. doi:10.1094/PHYTO-04-18-0130-IA.
    OpenUrlCrossRef
  8. 8.↵
    1. Judelson HS
    . 2017. Metabolic diversity and novelties in the oomycetes. Annu Rev Microbiol 71:21–39. doi:10.1146/annurev-micro-090816-093609.
    OpenUrlCrossRef
  9. 9.↵
    1. Judelson HS,
    2. Ah-Fong AM
    . 2018. Exchanges at the plant-oomycete interface that influence disease. Plant Physiol 179:1198–1211. doi:10.1104/pp.18.00979.
    OpenUrlCrossRef
  10. 10.↵
    1. Richards TA,
    2. Talbot NJ
    . 2013. Horizontal gene transfer in osmotrophs: playing with public goods. Nat Rev Microbiol 11:720–727. doi:10.1038/nrmicro3108.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Zuluaga AP,
    2. Vega-Arreguin JC,
    3. Fei Z,
    4. Matas AJ,
    5. Patev S,
    6. Fry WE,
    7. Rose J
    . 2016. Analysis of the tomato leaf transcriptome during successive hemibiotrophic stages of a compatible interaction with the oomycete pathogen Phytophthora infestans. Mol Plant Pathol 17:42–54. doi:10.1111/mpp.12260.
    OpenUrlCrossRef
  12. 12.↵
    1. van West P,
    2. De Jong AJ,
    3. Judelson HS,
    4. Emons AMC,
    5. Govers F
    . 1998. The ipiO gene of Phytophthora infestans is highly expressed in invading hyphae during infection. Fungal Genet Biol 23:126–138. doi:10.1006/fgbi.1998.1036.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    1. Ah-Fong AM,
    2. Shrivastava J,
    3. Judelson HS
    . 2017. Lifestyle, gene gain and loss, and transcriptional remodeling cause divergence in the transcriptomes of Phytophthora infestans and Pythium ultimum during potato tuber colonization. BMC Genomics 18:764. doi:10.1186/s12864-017-4151-2.
    OpenUrlCrossRef
  14. 14.↵
    1. Judelson HS,
    2. Tani S,
    3. Narayan RD
    . 2009. Metabolic adaptation of Phytophthora infestans during growth on leaves, tubers and artificial media. Mol Plant Pathol 10:843–855. doi:10.1111/j.1364-3703.2009.00570.x.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Divon HH,
    2. Fluhr R
    . 2007. Nutrition acquisition strategies during fungal infection of plants. FEMS Microbiol Lett 266:65–74. doi:10.1111/j.1574-6968.2006.00504.x.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Abrahamian M,
    2. Ah-Fong AM,
    3. Davis C,
    4. Andreeva K,
    5. Judelson HS
    . 2016. Gene expression and silencing studies in Phytophthora infestans reveal infection-specific nutrient transporters and a role for the nitrate reductase pathway in plant pathogenesis. PLoS Pathog 12:e1006097. doi:10.1371/journal.ppat.1006097.
    OpenUrlCrossRef
  17. 17.↵
    1. Abrahamian M,
    2. Kagda M,
    3. Ah-Fong AM,
    4. Judelson HS
    . 2017. Rethinking the evolution of eukaryotic metabolism: novel cellular partitioning of enzymes in stramenopiles links serine biosynthesis to glycolysis in mitochondria. BMC Evol Biol 17:241. doi:10.1186/s12862-017-1087-8.
    OpenUrlCrossRef
  18. 18.↵
    1. Ah-Fong AM,
    2. Kim KS,
    3. Judelson HS
    . 2017. RNA-seq of life stages of the oomycete Phytophthora infestans reveals dynamic changes in metabolic, signal transduction, and pathogenesis genes and a major role for calcium signaling in development. BMC Genomics 18:198. doi:10.1186/s12864-017-3585-x.
    OpenUrlCrossRef
  19. 19.↵
    1. Niu X,
    2. Ah-Fong AM,
    3. Lopez LA,
    4. Judelson HS
    . 2018. Transcriptomic and proteomic analysis reveals wall-associated and glucan-degrading proteins with potential roles in Phytophthora infestans sexual spore development. PLoS One 13:e0198186. doi:10.1371/journal.pone.0198186.
    OpenUrlCrossRef
  20. 20.↵
    1. Nielsen J
    . 2017. Systems biology of metabolism. Annu Rev Biochem 86:245–275. doi:10.1146/annurev-biochem-061516-044757.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. O’Brien EJ,
    2. Monk JM,
    3. Palsson BO
    . 2015. Using genome-scale models to predict biological capabilities. Cell 161:971–987. doi:10.1016/j.cell.2015.05.019.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Chavali AK,
    2. D’Auria KM,
    3. Hewlett EL,
    4. Pearson RD,
    5. Papin JA
    . 2012. A metabolic network approach for the identification and prioritization of antimicrobial drug targets. Trends Microbiol 20:113–123. doi:10.1016/j.tim.2011.12.004.
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.↵
    1. Cesur MF,
    2. Abdik E,
    3. Güven-Gülhan Ü,
    4. Durmus S,
    5. Çakir T
    . 2018. Computational systems biology of metabolism in infection, p 235–282. In Silvestre R, Torrado E (ed), Metabolic interaction in infection. Springer, Cham, Switzerland.
  24. 24.↵
    1. Peyraud R,
    2. Dubiella U,
    3. Barbacci A,
    4. Genin S,
    5. Raffaele S,
    6. Roby D
    . 2017. Advances on plant-pathogen interactions from molecular toward systems biology perspectives. Plant J 90:720–737. doi:10.1111/tpj.13429.
    OpenUrlCrossRef
  25. 25.↵
    1. Olive AJ,
    2. Sassetti CM
    . 2016. Metabolic crosstalk between host and pathogen: sensing, adapting and competing. Nat Rev Microbiol 14:221–234. doi:10.1038/nrmicro.2016.12.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Spanu PD
    . 2012. The genomics of obligate (and nonobligate) biotrophs. Annu Rev Phytopathol 50:91–109. doi:10.1146/annurev-phyto-081211-173024.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.↵
    1. Huthmacher C,
    2. Hoppe A,
    3. Bulik S,
    4. Holzhütter HG
    . 2010. Antimalarial drug targets in Plasmodium falciparum predicted by stage-specific metabolic network analysis. BMC Syst Biol 4:120. doi:10.1186/1752-0509-4-120.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Magnúsdóttir S,
    2. Heinken A,
    3. Kutt L,
    4. Ravcheev DA,
    5. Bauer E,
    6. Noronha A,
    7. Greenhalgh K,
    8. Jäger C,
    9. Baginska J,
    10. Wilmes P,
    11. Fleming RMT,
    12. Thiele I
    . 2017. Generation of genome-scale metabolic reconstructions for 773 members of the human gut microbiota. Nat Biotechnol 35:81–89. doi:10.1038/nbt.3703.
    OpenUrlCrossRef
  29. 29.↵
    1. Raghunathan A,
    2. Reed J,
    3. Shin S,
    4. Palsson BO,
    5. Daefler S
    . 2009. Constraint-based analysis of metabolic capacity of Salmonella typhimurium during host-pathogen interaction. BMC Syst Biol 3:38. doi:10.1186/1752-0509-3-38.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Raghunathan A,
    2. Shin S,
    3. Daefler S
    . 2010. Systems approach to investigating host-pathogen interactions in infections with the biothreat agent Francisella. Constraints-based model of Francisella tularensis. BMC Syst Biol 4:118. doi:10.1186/1752/0509-4-118.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Bordbar A,
    2. Lewis NE,
    3. Schellenberger J,
    4. Palsson BO,
    5. Jamshidi N
    . 2010. Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions. Mol Syst Biol 6:422. doi:10.1038/msb.2010.68.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Aller S,
    2. Scott A,
    3. Sarkar-Tyson M,
    4. Soyer OS
    . 2018. Integrated human-virus metabolic stoichiometric modelling predicts host-based antiviral targets against Chikungunya, Dengue and Zika viruses. J R Soc Interface 15:20180125. doi:10.1098/rsif.2018.0125.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Ankrah NYD,
    2. Chouaia B,
    3. Douglas AE
    . 2018. The cost of metabolic interactions in symbioses between insects and bacteria with reduced genomes. mBio 9:e01433-18. doi:10.1128/mBio.01433-18.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Rodenburg SYA,
    2. Seidl MF,
    3. de Ridder D,
    4. Govers F
    . 2018. Genome-wide characterization of Phytophthora infestans metabolism: a systems biology approach. Mol Plant Pathol 19:1403–1413. doi:10.1111/mpp.12623.
    OpenUrlCrossRef
  35. 35.↵
    1. Yuan H,
    2. Cheung CYM,
    3. Poolman MG,
    4. Hilbers PAJ,
    5. Van Riel N
    . 2016. A genome-scale metabolic network reconstruction of tomato (Solanum lycopersicum L.) and its application to photorespiratory metabolism. Plant J 85:289–304. doi:10.1111/tpj.13075.
    OpenUrlCrossRef
  36. 36.↵
    1. Kanehisa M,
    2. Sato Y,
    3. Kawashima M,
    4. Furumichi M,
    5. Tanabe M
    . 2015. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res 44:D457–D462. doi:10.1093/nar/gkv1070.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Caspi R,
    2. Altman T,
    3. Billington R,
    4. Dreher K,
    5. Foerster H,
    6. Fulcher CA,
    7. Holland TA,
    8. Keseler IM,
    9. Kothari A,
    10. Kubo A,
    11. Krummenacker M,
    12. Latendresse M,
    13. Mueller LA,
    14. Ong Q,
    15. Paley S,
    16. Subhraveti P,
    17. Weaver DS,
    18. Weerasinghe D,
    19. Zhang P,
    20. Karp PD
    . 2014. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nucleic Acids Res 42:D459–D471. doi:10.1093/nar/gkt1103.
    OpenUrlCrossRefPubMedWeb of Science
  38. 38.↵
    1. Savory FR,
    2. Milner DS,
    3. Miles DC,
    4. Richards TA
    . 2018. Ancestral function and diversification of a horizontally acquired oomycete carboxylic acid transporter. Mol Biol Evol 35:1887–1900. doi:10.1093/molbev/msy082.
    OpenUrlCrossRef
  39. 39.↵
    1. Orth JD,
    2. Thiele I,
    3. Palsson BO
    . 2010. What is flux balance analysis? Nat Biotechnol 28:245–248. doi:10.1038/nbt.1614.
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.↵
    1. Hodgson WA
    . 1958. Growth of four races of Phytophthora infestans (Mont.) De Bary in synthetic media. Can J Plant Sci 38:145–154. doi:10.4141/cjps58-026.
    OpenUrlCrossRef
  41. 41.↵
    1. Hohl HR
    . 1991. Nutrition, p 53–83. In Ingram DS, Williams PH (ed), Advances in plant pathology, vol 7. Academic Press, London, United Kingdom.
    OpenUrl
  42. 42.↵
    1. Azevedo RA,
    2. Lancien M,
    3. Lea PJ
    . 2006. The aspartic acid metabolic pathway, an exciting and essential pathway in plants. Amino Acids 30:143–162. doi:10.1007/s00726-005-0245-2.
    OpenUrlCrossRefPubMedWeb of Science
  43. 43.↵
    1. Jander G,
    2. Joshi V
    . 2009. Aspartate-derived amino acid biosynthesis in Arabidopsis thaliana. Arabidopsis Book 7:e0121. doi:10.1199/tab.0121.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Wang R,
    2. Zhang M,
    3. Liu H,
    4. Xu J,
    5. Yu J,
    6. He F,
    7. Zhang X,
    8. Dong S,
    9. Dou D
    . 2016. PsAAT3, an oomycete-specific aspartate aminotransferase, is required for full pathogenicity of the oomycete pathogen Phytophthora sojae. Fungal Biol 120:620–630. doi:10.1016/j.funbio.2016.01.005. (Corrigendum, 122:1228–1229, 2018, doi:10.1016/j.funbio.2018.09.004.)
    OpenUrlCrossRef
  45. 45.↵
    1. Chan XWA,
    2. Wrenger C,
    3. Stahl K,
    4. Bergmann B,
    5. Winterberg M,
    6. Müller IB,
    7. Saliba KJ
    . 2013. Chemical and genetic validation of thiamine utilization as an antimalarial drug target. Nat Commun 4:2060. doi:10.1038/ncomms3060.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Mukherjee S,
    2. Retwitzer MD,
    3. Barash D,
    4. Sengupta S
    . 2018. Phylogenomic and comparative analysis of the distribution and regulatory patterns of TPP riboswitches in fungi. Sci Rep 8:5563. doi:10.1038/s41598-018-23900-7.
    OpenUrlCrossRef
  47. 47.↵
    1. McRose D,
    2. Guo J,
    3. Monier A,
    4. Sudek S,
    5. Wilken S,
    6. Yan S,
    7. Mock T,
    8. Archibald JM,
    9. Begley TP,
    10. Reyes-Prieto A,
    11. Worden AZ
    . 2014. Alternatives to vitamin B1 uptake revealed with discovery of riboswitches in multiple marine eukaryotic lineages. ISME J 8:2517–2529. doi:10.1038/ismej.2014.146.
    OpenUrlCrossRef
  48. 48.↵
    1. Burgard AP,
    2. Nikolaev EV,
    3. Schilling CH,
    4. Maranas CD
    . 2004. Flux coupling analysis of genome-scale metabolic network reconstructions. Genome Res 14:301–312. doi:10.1101/gr.1926504.
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    1. Nelson DL,
    2. Cox MM
    . 2017. Lehninger principles of biochemistry. W. H. Freeman, New York, NY.
  50. 50.↵
    1. Goyer A
    . 2010. Thiamine in plants: aspects of its metabolism and functions. Phytochemistry 71:1615–1624. doi:10.1016/j.phytochem.2010.06.022.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Maridueña-Zavala MG,
    2. Freire-Peñaherrera A,
    3. Cevallos-Cevallos JM,
    4. Peralta EL
    . 2017. GC-MS metabolite profiling of Phytophthora infestans resistant to metalaxyl. Eur J Plant Pathol 149:563–574. doi:10.1007/s10658-017-1204-y.
    OpenUrlCrossRef
  52. 52.↵
    1. García-Bayona L,
    2. Garavito MF,
    3. Lozano GL,
    4. Vasquez JJ,
    5. Myers KL,
    6. Fry WE,
    7. Bernal A,
    8. Zimmermann BH,
    9. Restrepo S
    . 2014. De novo pyrimidine biosynthesis in the oomycete plant pathogen Phytophthora infestans. Gene 537:312–321. doi:10.1016/j.gene.2013.12.009.
    OpenUrlCrossRef
  53. 53.↵
    1. Griffiths RG,
    2. Dancer J,
    3. O'Neill E,
    4. Harwood JL
    . 2003. A mandelamide pesticide alters lipid metabolism in Phytophthora infestans. New Phytol 158:345–353. doi:10.1046/j.1469-8137.2003.00739.x.
    OpenUrlCrossRef
  54. 54.↵
    1. Agren R,
    2. Bordel S,
    3. Mardinoglu A,
    4. Pornputtapong N,
    5. Nookaew I,
    6. Nielsen J
    . 2012. Reconstruction of genome-scale active metabolic networks for 69 human cell types and 16 cancer types using INIT. PLoS Comput Biol 8:e1002518. doi:10.1371/journal.pcbi.1002518.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Judelson HS,
    2. Narayan RD,
    3. Ah-Fong AM,
    4. Kim KS
    . 2009. Gene expression changes during asexual sporulation by the late blight agent Phytophthora infestans occur in discrete temporal stages. Mol Genet Genomics 281:193–206. doi:10.1007/s00438-008-0407-5.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Fatima U,
    2. Senthil-Kumar M
    . 2015. Plant and pathogen nutrient acquisition strategies. Front Plant Sci 6:750. doi:10.3389/fpls.2015.00750.
    OpenUrlCrossRef
  57. 57.↵
    1. Chen LQ
    . 2014. SWEET sugar transporters for phloem transport and pathogen nutrition. New Phytol 201:1150–1155. doi:10.1111/nph.12445.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Savchenko T,
    2. Walley JW,
    3. Chehab EW,
    4. Xiao Y,
    5. Kaspi R,
    6. Pye MF,
    7. Mohamed ME,
    8. Lazarus CM,
    9. Bostock RM,
    10. Dehesh K
    . 2010. Arachidonic acid: an evolutionarily conserved signaling molecule modulates plant stress signaling networks. Plant Cell 22:3193–3205. doi:10.1105/tpc.110.073858.
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    1. Robinson SM,
    2. Bostock RM
    . 2015. β-Glucans and eicosapolyenoic acids as MAMPs in plant-oomycete interactions: past and present. Front Plant Sci 5:797. doi:10.3389/fpls.2014.00797.
    OpenUrlCrossRef
  60. 60.↵
    1. Peyraud R,
    2. Cottret L,
    3. Marmiesse L,
    4. Genin S
    . 2018. Control of primary metabolism by a virulence regulatory network promotes robustness in a plant pathogen. Nat Commun 9:418. doi:10.1038/s41467-017-02660-4.
    OpenUrlCrossRef
  61. 61.↵
    1. Zhang C,
    2. Bidkhori G,
    3. Benfeitas R,
    4. Lee S,
    5. Arif M,
    6. Uhlén M,
    7. Mardinoglu A
    . 2018. ESS: a tool for genome-scale quantification of essentiality score for reaction/genes in constraint-based modeling. Front Physiol 9:1355. doi:10.3389/fphys.2018.01355.
    OpenUrlCrossRef
  62. 62.↵
    1. Jang C,
    2. Chen L,
    3. Rabinowitz JD
    . 2018. Metabolomics and isotope tracing. Cell 173:822–837. doi:10.1016/j.cell.2018.03.055.
    OpenUrlCrossRefPubMed
  63. 63.↵
    1. Kleessen S,
    2. Irgang S,
    3. Klie S,
    4. Giavalisco P,
    5. Nikoloski Z
    . 2015. Integration of transcriptomics and metabolomics data specifies the metabolic response of Chlamydomonas to rapamycin treatment. Plant J 81:822–835. doi:10.1111/tpj.12763.
    OpenUrlCrossRefPubMed
  64. 64.↵
    1. Robbins WJ
    . 1938. Thiamin and growth of species of Phytophthora. Bull Torrey Bot Club 65:267. doi:10.2307/2481231.
    OpenUrlCrossRef
  65. 65.↵
    1. Peyraud R,
    2. Mbengue M,
    3. Barbacci A,
    4. Raffaele S
    . 2019. Intercellular cooperation in a fungal plant pathogen facilitates host colonization. Proc Natl Acad Sci U S A 116:3193–3201. doi:10.1073/pnas.1811267116.
    OpenUrlAbstract/FREE Full Text
  66. 66.↵
    1. Galeano Garcia P,
    2. Neves dos Santos F,
    3. Zanotta S,
    4. Eberlin M,
    5. Carazzone C
    . 2018. Metabolomics of Solanum lycopersicum infected with Phytophthora infestans leads to early detection of late blight in asymptomatic plants. Molecules 23:3330. doi:10.3390/molecules23123330.
    OpenUrlCrossRef
  67. 67.↵
    1. Blume M,
    2. Seeber F
    . 2018. Metabolic interactions between Toxoplasma gondii and its host. F1000Res 7:1719. doi:10.12688/f1000research.16021.1.
    OpenUrlCrossRef
  68. 68.↵
    1. Meier A,
    2. Erler H,
    3. Beitz E
    . 2018. Targeting channels and transporters in protozoan parasite infections. Front Chem 6:88. doi:10.3389/fchem.2018.00088.
    OpenUrlCrossRef
  69. 69.↵
    1. Beste DV,
    2. Nöh K,
    3. Niedenführ S,
    4. Mendum TA,
    5. Hawkins ND,
    6. Ward JL,
    7. Beale MH,
    8. Wiechert W,
    9. McFadden J
    . 2013. 13C-flux spectral analysis of host-pathogen metabolism reveals a mixed diet for intracellular Mycobacterium tuberculosis. Chem Biol 20:1012–1021. doi:10.1016/j.chembiol.2013.06.012.
    OpenUrlCrossRefPubMed
  70. 70.↵
    1. Thiele I,
    2. Palsson BO
    . 2010. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc 5:93–121. doi:10.1038/nprot.2009.203.
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    1. Moretti S,
    2. Martin O,
    3. Van Du Tran T,
    4. Bridge A,
    5. Morgat A,
    6. Pagni M
    . 2016. MetaNetX/MNXref–reconciliation of metabolites and biochemical reactions to bring together genome-scale metabolic networks. Nucleic Acids Res 44:D523–D526. doi:10.1093/nar/gkv1117.
    OpenUrlCrossRefPubMed
  72. 72.↵
    1. Yilmaz JL,
    2. Lim ZL,
    3. Beganovic M,
    4. Breazeale S,
    5. Andre C,
    6. Stymne S,
    7. Vrinten P,
    8. Senger T
    . 2017. Determination of substrate preferences for desaturases and elongases for production of docosahexaenoic acid from oleic acid in engineered canola. Lipids 52:207–222. doi:10.1007/s11745-017-4235-4.
    OpenUrlCrossRef
  73. 73.↵
    1. Basenko E,
    2. Pulman J,
    3. Shanmugasundram A,
    4. Harb O,
    5. Crouch K,
    6. Starns D,
    7. Warrenfeltz S,
    8. Aurrecoechea C,
    9. Stoeckert C, Jr.,
    10. Kissinger J,
    11. Roos DS,
    12. Hertz-Fowler C
    . 2018. FungiDB: an integrated bioinformatic resource for fungi and oomycetes. J Fungi 4:39. doi:10.3390/jof4010039.
    OpenUrlCrossRef
  74. 74.↵
    1. Haas BJ,
    2. Kamoun S,
    3. Zody MC,
    4. Jiang RHY,
    5. Handsaker RE,
    6. Cano LM,
    7. Grabherr M,
    8. Kodira CD,
    9. Raffaele S,
    10. Torto-Alalibo TA,
    11. Bozkurt TO,
    12. Ah-Fong AM,
    13. Alvarado L,
    14. Anderson VL,
    15. Armstrong MR,
    16. Avrova AO,
    17. Baxter L,
    18. Beynon JL,
    19. Boevink PC,
    20. Bollmann SR,
    21. Bos JIB,
    22. Bulone V,
    23. Cai G,
    24. Cakir C,
    25. Carrington JC,
    26. Chawner M,
    27. Conti L,
    28. Costanzo S,
    29. Ewan R,
    30. Fahlgren N,
    31. Fischbach MA,
    32. Fugelstad J,
    33. Gilroy EM,
    34. Gnerre S,
    35. Green PJ,
    36. Grenville-Briggs LJ,
    37. Griffith J,
    38. Grunwald NJ,
    39. Horn K,
    40. Horner NR,
    41. Hu CH,
    42. Huitema E,
    43. Jeong DH,
    44. Jones AME,
    45. Jones JDG,
    46. Jones RW,
    47. Karlsson EK,
    48. Kunjeti SG,
    49. Lamour KH,
    50. Liu Z, et al
    . 2009. Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans. Nature 461:393–398. doi:10.1038/nature08358.
    OpenUrlCrossRefPubMedWeb of Science
  75. 75.↵
    1. Goldberg T,
    2. Hecht M,
    3. Hamp T,
    4. Karl T,
    5. Yachdav G,
    6. Ahmed N,
    7. Altermann U,
    8. Angerer P,
    9. Ansorge S,
    10. Balasz K,
    11. Bernhofer M,
    12. Betz A,
    13. Cizmadija L,
    14. Do KT,
    15. Gerke J,
    16. Greil R,
    17. Joerdens V,
    18. Hastreiter M,
    19. Hembach K,
    20. Herzog M,
    21. Kalemanov M,
    22. Kluge M,
    23. Meier A,
    24. Nasir H,
    25. Neumaier U,
    26. Prade V,
    27. Reeb J,
    28. Sorokoumov A,
    29. Troshani I,
    30. Vorberg S,
    31. Waldraff S,
    32. Zierer J,
    33. Nielsen H,
    34. Rost B
    . 2014. LocTree3 prediction of localization. Nucleic Acids Res 42:W350–W355. doi:10.1093/nar/gku396.
    OpenUrlCrossRefPubMedWeb of Science
  76. 76.↵
    1. Guan J-C,
    2. Hasnain G,
    3. Garrett TJ,
    4. Chase CD,
    5. Gregory J,
    6. Hanson AD,
    7. McCarty DR
    . 2014. Divisions of labor in the thiamin biosynthetic pathway among organs of maize. Front Plant Sci 5:370. doi:10.3389/fpls.2014.00370.
    OpenUrlCrossRefPubMed
  77. 77.↵
    1. Wang H,
    2. Marcišauskas S,
    3. Sánchez BJ,
    4. Domenzain I,
    5. Hermansson D,
    6. Agren R,
    7. Nielsen J,
    8. Kerkhoven EJ
    . 2018. RAVEN 2.0: a versatile toolbox for metabolic network reconstruction and a case study on Streptomyces coelicolor. PLoS Comput Biol 14:e1006541. doi:10.1371/journal.pcbi.1006541.
    OpenUrlCrossRefPubMed
  78. 78.↵
    1. Heirendt L,
    2. Arreckx S,
    3. Pfau T,
    4. Mendoza SN,
    5. Richelle A,
    6. Heinken A,
    7. Haraldsdóttir HS,
    8. Wachowiak J,
    9. Keating SM,
    10. Vlasov V,
    11. Magnúsdóttir S,
    12. Ng CY,
    13. Preciat G,
    14. Žagare A,
    15. Chan SHJ,
    16. Aurich MK,
    17. Clancy CM,
    18. Modamio J,
    19. Sauls JT,
    20. Noronha A,
    21. Bordbar A,
    22. Cousins B,
    23. El Assal DC,
    24. Valcarcel LV,
    25. Apaolaza I,
    26. Ghaderi S,
    27. Ahookhosh M,
    28. Ben Guebila M,
    29. Kostromins A,
    30. Sompairac N,
    31. Le HM,
    32. Ma D,
    33. Sun Y,
    34. Wang L,
    35. Yurkovich JT,
    36. Oliveira MAP,
    37. Vuong PT,
    38. El Assal LP,
    39. Kuperstein I,
    40. Zinovyev A,
    41. Hinton HS,
    42. Bryant WA,
    43. Artacho FJA,
    44. Planes FJ,
    45. Stalidzans E,
    46. Maass A,
    47. Vempala S,
    48. Hucka M,
    49. Saunders MA,
    50. Maranas CD, et al
    . 2017. Creation and analysis of biochemical constraint-based models: the COBRA Toolbox v3.0. arXiv arXiv:1710.04038.
    OpenUrl
  79. 79.↵
    1. Larhlimi A,
    2. David L,
    3. Selbig J,
    4. Bockmayr A
    . 2012. F2C2: a fast tool for the computation of flux coupling in genome-scale metabolic networks. BMC Bioinformatics 13:57. doi:10.1186/1471-2105-13-57.
    OpenUrlCrossRefPubMed
  80. 80.↵
    1. Pratapa A,
    2. Balachandran S,
    3. Raman K
    . 2015. Fast-SL: an efficient algorithm to identify synthetic lethal sets in metabolic networks. Bioinformatics 31:3299–3305. doi:10.1093/bioinformatics/btv352.
    OpenUrlCrossRefPubMed
  81. 81.↵
    1. Sato S,
    2. Tabata S,
    3. Hirakawa H,
    4. Asamizu E,
    5. Shirasawa K,
    6. Isobe S,
    7. Kaneko T,
    8. Nakamura Y,
    9. Shibata D,
    10. Aoki K,
    11. Egholm M,
    12. Knight J,
    13. Bogden R,
    14. Li C,
    15. Shuang Y,
    16. Xu X,
    17. Pan S,
    18. Cheng S,
    19. Liu X,
    20. Ren Y,
    21. Wang J,
    22. Albiero A,
    23. Dal Pero F,
    24. Todesco S,
    25. Van Eck J,
    26. Buels RM,
    27. Bombarely A,
    28. Gosselin JR,
    29. Huang M,
    30. Leto JA,
    31. Menda N,
    32. Strickler S,
    33. Mao L,
    34. Gao S,
    35. Tecle IY,
    36. York T,
    37. Zheng Y,
    38. Vrebalov JT,
    39. Lee J,
    40. Zhong S,
    41. Mueller LA,
    42. Stiekema WJ,
    43. Ribeca P,
    44. Alioto T,
    45. Yang W,
    46. Huang S,
    47. Du Y,
    48. Zhang Z,
    49. Gao J,
    50. Guo Y, et al
    . 2012. The tomato genome sequence provides insights into fleshy fruit evolution. Nature 485:635–641. doi:10.1038/nature11119.
    OpenUrlCrossRefPubMedWeb of Science
  82. 82.↵
    1. Kim D,
    2. Langmead B,
    3. Salzberg SL
    . 2015. HISAT: a fast spliced aligner with low memory requirements. Nat Methods 12:357–360. doi:10.1038/nmeth.3317.
    OpenUrlCrossRefPubMed
  83. 83.↵
    1. Li H,
    2. Handsaker B,
    3. Wysoker A,
    4. Fennell T,
    5. Ruan J,
    6. Homer N,
    7. Marth G,
    8. Abecasis G,
    9. Durbin R
    , 1000 Genome Project Data Processing Subgroup. 2009. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078–2079. doi:10.1093/bioinformatics/btp352.
    OpenUrlCrossRefPubMedWeb of Science
  84. 84.↵
    1. Trapnell C,
    2. Williams BA,
    3. Pertea G,
    4. Mortazavi A,
    5. Kwan G,
    6. Van Baren MJ,
    7. Salzberg SL,
    8. Wold BJ,
    9. Pachter L
    . 2010. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28:511–515. doi:10.1038/nbt.1621.
    OpenUrlCrossRefPubMedWeb of Science
  85. 85.↵
    1. Robinson MD,
    2. Oshlack A
    . 2010. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol 11:R25. doi:10.1186/gb-2010-11-3-r25.
    OpenUrlCrossRefPubMed
  86. 86.↵
    1. Hochberg B
    . 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol 57:289–300.
    OpenUrl
PreviousNext
Back to top
Download PDF
Citation Tools
Metabolic Model of the Phytophthora infestans-Tomato Interaction Reveals Metabolic Switches during Host Colonization
Sander Y. A. Rodenburg, Michael F. Seidl, Howard S. Judelson, Andrea L. Vu, Francine Govers, Dick de Ridder
mBio Jul 2019, 10 (4) e00454-19; DOI: 10.1128/mBio.00454-19

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Metabolic Model of the Phytophthora infestans-Tomato Interaction Reveals Metabolic Switches during Host Colonization
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Metabolic Model of the Phytophthora infestans-Tomato Interaction Reveals Metabolic Switches during Host Colonization
Sander Y. A. Rodenburg, Michael F. Seidl, Howard S. Judelson, Andrea L. Vu, Francine Govers, Dick de Ridder
mBio Jul 2019, 10 (4) e00454-19; DOI: 10.1128/mBio.00454-19
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS and DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Phytophthora infestans
metabolic modeling
metabolism
oomycetes
tomato

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511