Skip to main content

Main menu

  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
  • About
    • Editorial Board and Staff
    • About the Journal
    • Terms & Privacy
  • More
    • Alerts
    • Contact Us
  • Submit a Manuscript
    • Instructions for Authors
    • Submit a Manuscript
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Teaching Tools in Plant Biology
    • ASPB
    • Plantae

User menu

  • My alerts
  • Log in

Search

  • Advanced search
Plant Cell
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Teaching Tools in Plant Biology
    • ASPB
    • Plantae
  • My alerts
  • Log in
Plant Cell

Advanced Search

  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
  • About
    • Editorial Board and Staff
    • About the Journal
    • Terms & Privacy
  • More
    • Alerts
    • Contact Us
  • Submit a Manuscript
    • Instructions for Authors
    • Submit a Manuscript
  • Follow PlantCell on Twitter
  • Visit PlantCell on Facebook
  • Visit Plantae
Research ArticleLARGE-SCALE BIOLOGY ARTICLES
Open Access

Plastic Transcriptomes Stabilize Immunity to Pathogen Diversity: The Jasmonic Acid and Salicylic Acid Networks within the Arabidopsis/Botrytis Pathosystem

Wei Zhang, Jason A. Corwin, Daniel Copeland, Julie Feusier, Robert Eshbaugh, Fang Chen, Susana Atwell, Daniel J. Kliebenstein
Wei Zhang
aDepartment of Plant Sciences, University of California, Davis, California 95616
bNational and Local Joint Engineering Laboratory for Energy Plant Bio-oil Production and Application, Key Laboratory of Bio-resource and Eco-environment, Ministry of Education, College of Life Sciences, Sichuan University, Chengdu 610064, P.R. China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Wei Zhang
Jason A. Corwin
cDepartment of Ecology and Evolution Biology, University of Colorado, Boulder, Colorado 80309-0334
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jason A. Corwin
Daniel Copeland
aDepartment of Plant Sciences, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Daniel Copeland
Julie Feusier
aDepartment of Plant Sciences, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Julie Feusier
Robert Eshbaugh
aDepartment of Plant Sciences, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fang Chen
bNational and Local Joint Engineering Laboratory for Energy Plant Bio-oil Production and Application, Key Laboratory of Bio-resource and Eco-environment, Ministry of Education, College of Life Sciences, Sichuan University, Chengdu 610064, P.R. China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Susana Atwell
aDepartment of Plant Sciences, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel J. Kliebenstein
aDepartment of Plant Sciences, University of California, Davis, California 95616
dDynaMo Center of Excellence, University of Copenhagen, DK-1871 Frederiksberg C, Denmark
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Daniel J. Kliebenstein
  • For correspondence: kliebenstein@ucdavis.edu

Published November 2017. DOI: https://doi.org/10.1105/tpc.17.00348

  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading
  • © 2017 American Society of Plant Biologists. All rights reserved.

Abstract

To respond to pathogen attack, selection and associated evolution has led to the creation of plant immune system that are a highly effective and inducible defense system. Central to this system are the plant defense hormones jasmonic acid (JA) and salicylic acid (SA) and crosstalk between the two, which may play an important role in defense responses to specific pathogens or even genotypes. Here, we used the Arabidopsis thaliana-Botrytis cinerea pathosystem to test how the host’s defense system functions against genetic variation in a pathogen. We measured defense-related phenotypes and transcriptomic responses in Arabidopsis wild-type Col-0 and JA- and SA-signaling mutants, coi1-1 and npr1-1, individually challenged with 96 diverse B. cinerea isolates. Those data showed genetic variation in the pathogen influences on all components within the plant defense system at the transcriptional level. We identified four gene coexpression networks and two vectors of defense variation triggered by genetic variation in B. cinerea. This showed that the JA and SA signaling pathways functioned to constrain/canalize the range of virulence in the pathogen population, but the underlying transcriptomic response was highly plastic. These data showed that plants utilize major defense hormone pathways to buffer disease resistance, but not the metabolic or transcriptional responses to genetic variation within a pathogen.

INTRODUCTION

The photosynthetic capacity of plants to generate much of their own energy and resources makes them attractive targets for exploitative species, such as microbial pathogens. These attacks can lead to strong selective forces on the plant genome in wild ecosystems and can cause significant financial and yield losses in agricultural crops (Bergelson et al., 2001; Van der Hoorn et al., 2002; Mauricio et al., 2003; Rose et al., 2007; Fisher et al., 2012; Hörger et al., 2012; Züst et al., 2012; Karasov et al., 2014). In turn, pathogens have evolved an array of virulence mechanisms based on their infection strategy that causes a prolonged selection pressure on plants to evolve complex and coherent innate immune systems. The resulting immune system must simultaneously detect and interpret a breadth of potential pathogen signals to mount an appropriate response tailored to the pathogen’s infection strategy. For instance, pathogens can be classified based on their feeding requirements (biotrophic versus necrotrophic pathogen) or by their host plant range (specialist versus generalist). Biotrophic pathogens, such as powdery mildews and rusts, infect the host and strive to circumvent triggering the host immune system through the use of effector proteins, which inhibit the detection of the pathogen and/or triggering of the plant immune system (Stergiopoulos and de Wit, 2009; Spanu et al., 2010; Okmen and Doehlemann, 2014). This allows the pathogen to absorb nutrients from the living host cell while avoiding host defenses. This strategy can drive specialization in the pathogen to focus on infecting one or a few species of host plants. By contrast, necrotrophic pathogens are often generalists that can have a broad host range throughout plants and are frequently characterized by a high level of intraspecific genetic diversity relative to specialist pathogens (Bolton et al., 2006; Williamson et al., 2007). These generalist pathogens, such as Botrytis cinerea and Sclerotinia sclerotiorum, actively kill host cells through a combination of toxic metabolites, lytic enzymes, and microRNAs (Lyu et al., 2016). The ubiquitous presence of microbes that differ in host range and infection strategy within an environment means that an individual plant likely interacts with multiple organisms that have a diverse array of infection strategies, and the plant must interpret these disparate signals to mount a coherent defense. To optimize the coordination of these defense responses, the plant utilizes a complex array of mechanisms to perceive the presence of the pathogen and then to transmit this signal via an interacting set of hormone signal transduction pathways to specifically modulate actual defense outputs. However, it is not presently known how these defense signaling pathways coordinately function in the presence of pathogens that have diverse virulence strategies.

The first line of active defense against specialist biotrophic pathogens is the plant eukaryotic innate immune system, which detects the presence of a pathogen and activates pattern-triggered immunity (PTI) (Jones and Dangl, 2006). PTI is a general form of qualitative resistance that functions against nonadapted, nonspecialist pathogens. To circumvent PTI responses, successful specialist pathogens often secrete or inject effector proteins that target diverse host cellular processes to evade or suppress PTI signaling. Selection then drives the formation of plant resistance genes (R genes), such as nucleotide binding site leucine-rich repeat proteins, that can either directly or indirectly detect the presence of the pathogen effectors and trigger effector-triggered immunity (ETI). This creates an evolutionary Red Queen selective cycle (Z-scheme), whereby more resistance genes are created in an ongoing arms race between plant defense and pathogen host immune suppression (Van Valen, 1973; Bergelson et al., 2001; Dangl and Jones, 2001; Jones and Dangl, 2006; Benton, 2009; Pritchard and Birch, 2014; Fei et al., 2016; Iakovidis et al., 2016; Keller et al., 2016). The foundation of this highly specific and mechanistic model of pathogen detection assumes an ecologically intimate specialist connection between host and pathogen, where the individual pathogen species or isolate exists primarily on one or a few plant species. By contrast, the vast ecological host ranges of even individual isolates of generalist pathogens like B. cinerea do not readily fit into the assumptions of this mechanistic model, as this pathogen can simply evade any new resistance mechanism by moving to another equally suitable plant host. While some molecules, such as chitin and mannitol, can trigger PTI to provide resistance against generalist fungi, it is not fully understood how generalist pathogens fit into the intricate host/pathogen coevolution assumption foundational to ETI/PTI (Voegele et al., 2005; Staal et al., 2008; Stergiopoulos and de Wit, 2009; Liu et al., 2012; Meena et al., 2015; Patel and Williamson, 2016).

After the initial pathogen detection by the PTI/ETI systems, these signals integrate via complex plant signaling networks that control the expression of plant defenses to shape downstream defense strategies. These defense transduction and coordination pathways are shaped by the defense hormones salicylic acid (SA) and jasmonic acid (JA) that function to shape downstream immune responses. Commonly, the literature links the SA pathway to defense responses against biotrophic pathogens and the JA pathway to necrotroph/herbivore responses (Vlot et al., 2009; He et al., 2012; Yang et al., 2015). The accumulation of SA and the change of cellular redox state release the defense regulatory protein NONEXPRESSOR OF PATHOGENESIS GENES1 (NPR1) to translocate to the nucleus and interact with members of the basic leucine zipper transcription factor family and induce defense responses (Clarke et al., 1998; Mou et al., 2003; Pieterse and Van Loon, 2004; Mukhtar et al., 2009; Moreau et al., 2012; Wu et al., 2012; Yan and Dong, 2014; Saleh et al., 2015). The JA response is mediated by the bioactive JA-Ile conjugate stimulating the interaction of the E3 ubiquitin ligase SCFCOI1 with the JA-ZIM-domain (JAZ) (Xie et al., 1998; Xu et al., 2002; Wasternack, 2007; Browse, 2009). This leads to the degradation of the JAZ proteins that repress transcription of JA-responsive genes by binding transcriptional activators (i.e., MYC2 and the APETALA2/ethylene-responsive element binding factor [AP2/ERF] transcription factor ORA59) (Thines et al., 2007; Pré et al., 2008; Chini et al., 2009; Sheard et al., 2010; Zarei et al., 2011; Zander et al., 2014; Zhang et al., 2015b). At the level of individual gene expression, the JA and SA signaling pathways appear to competitively intercommunicate, or crosstalk, which can result in antagonistic molecular interactions between these two pathways, potentially to optimize the defense to the specific attacker (Thaler et al., 2012; Caarls et al., 2015). Most studies of these signaling pathways have largely been conducted using specialist pathogens and/or individual genotypes of other pathogens with the assumption that these pathways function similarly across a diverse array of pathogen genotypes of the same or different taxonomic groups. However, given that generalists can have genetically diverse pathogen genotypes and contain significant variation in virulence mechanisms, it is possible that studying how this variation interacts with plant defense pathways could affect our view of how plant defense signaling functions. One possibility is that the defense signaling pathways help the plant to canalize its defense response across a broad array of pathogen genetic diversity. In this context, we define canalization as production by the plant host of a narrow range of phenotypes when confronted with a wide range of pathogen diversity (Waddington, 1942, 1959; Queitsch et al., 2002; Lachowiec et al., 2016). Another response to a diverse array of virulence mechanisms could be a plastic response where the range remains the same but the specific response is shifted, or finally, no change in canalization or plasticity. The degree to which each of these phenomena contribute to host defense responses remains an open question.

The Arabidopsis thaliana-B. cinerea pathosystem is an ideal system to investigate how the plant innate immune system coordinates decisions in the presence of natural intraspecific genetic diversity within the pathogen. B. cinerea is a necrotrophic pathogen that causes significant crop losses both pre- and postharvest, most dramatically in the developing world (Williamson et al., 2007). Additionally, B. cinerea is a true haploid ascomycete that infects a wide range of evolutionarily distinct plant hosts, from bryophytes to eudicots. B. cinerea has elevated natural genetic variation that results in multiple major-effect polymorphisms in known virulence mechanisms, including the production of phytotoxic metabolites (Colmenares et al., 2002; Dalmais et al., 2011), enzymes that detoxify plant defense metabolites (Ferrari et al., 2003; Pedras et al., 2005, 2007, 2008, 2009, 2011; Stefanato et al., 2009; Rowe et al., 2010), and the ability to degrade plant cell walls (Rowe and Kliebenstein, 2007; Schumacher et al., 2012, 2015; Kumari et al., 2014). Because wild B. cinerea isolates have recombination and random mating, a population of isolates is a random intermixed sample of the diverse virulence mechanisms (Rowe and Kliebenstein, 2007; Kretschmer et al., 2009; Rowe et al., 2010; Kumari et al., 2014; Atwell et al., 2015; Corwin et al., 2016a, 2016b; Zhang et al., 2016). This allows us to use a population of isolates to challenge a plant’s innate immune system with a diverse array of natural inputs from a single pathogen species and assess the flexibility of the signaling network (Finkers et al., 2007a, 2007b; Rowe and Kliebenstein, 2008; Davis et al., 2009; Corwin et al., 2016b; Zhang et al., 2016).

To test how plants respond to a high level of standing genetic variation in a generalist necrotrophic pathogen, we infected Arabidopsis with 96 genetically diverse natural B. cinerea isolates. This was done by infecting three Arabidopsis genotypes, including wild-type Columbia-0 (Col-0) and two immune-deficient mutants, coi1-1 (jasmonate insensitive) and npr1-1 (defective in some SA-mediated defenses), to allow us to measure the quantitative influence of the SA and JA defense pathways in the response to pathogen diversity and their role in canalization of defense responses. Defense response traits, such as lesion size, camalexin production, and the plant’s transcriptomic response were measured in a replicated and randomized complete block design, which allowed for all Arabidopsis × B. cinerea genotypic combinations to be tested. Using this experimental design, we partitioned the effects of genetic variation in the host and pathogen variation as well as their interactions. The variation for the vast majority of Arabidopsis transcripts was controlled by an equal interaction of hormone-perception variation in Arabidopsis and genetic variation in B. cinerea. By coupling principal component analysis (PCA) with gene coexpression networks, we identified four main host response networks that were the predominant targets of the genetic interaction between host and pathogen. These networks included genes focused on camalexin production, a suite of plastid-encoded genes and two other defense pathways. Furthermore, COI1 (CORONATINE INSENSITIVE1) and NPR1 appeared to canalize the response to the diverse pathogens when measuring resistance, but not when measuring transcriptomic responses. Instead, the SA and JA pathways appear to provide alternative molecular signaling routes to obtain similar phenotypic outcomes in response to the diverse pathogen isolates.

RESULTS

Natural Genetic Variation of B. cinerea Alters Lesion Area on Arabidopsis

To test how genetic variation within the species B. cinerea differentially probes the major plant defense pathways, we infected a collection of 96 diverse B. cinerea isolates onto the Col-0 accession of Arabidopsis (Supplemental Data Set 1) (Atwell et al., 2015; Corwin et al., 2016a). This collection of natural isolates was obtained mainly from California with a sampling of available international germplasm. They represent a wide range of genomic, phenotypic, and virulence variation (Rowe and Kliebenstein, 2007; Atwell et al., 2015; Corwin et al., 2016a). In conjunction with wild-type Col-0, we also included the immuno-compromised, single-gene knockout mutants coi1-1 and npr1-1, which abolish or diminish the major defense processing of the JA and SA pathways respectively (Xie et al., 1998; Xu et al., 2002; Thines et al., 2007; Chini et al., 2009; Sheard et al., 2010; Yan and Dong, 2014). Plants were infected with 6-fold replication across two independent experiments using a randomized complete block design for a total of 12 infected leaves per isolate/mutant pair. Two infected leaves per experiment were collected at 16 h postinfection (HPI) to measure plant transcript abundance in response to the diverse B. cinerea genotypes. The remaining four infected leaves within each experiment were imaged at 72 HPI to measure lesion area and collected for camalexin extraction. Each individual infection consists of a single 4 μL droplet of 40 spores diluted in 0.5× organic grape juice. While this provides a small energy boost to the germinating spores, B. cinerea also typically requires sugars as germination signals and this helps to ensure consistent germination (Cotoras et al., 2009). This detached leaf assay system has been successfully utilized to understand the Arabidopsis-B. cinerea pathosystem (ten Have et al., 1998; Thomma et al., 1999a; Ferrari et al., 2003; Denby et al., 2004; Kliebenstein et al., 2005; Han et al., 2010; Windram et al., 2012). As with any experimental system, changing the inoculum media may alter the specific patterns observed.

To determine whether natural genetic variation in the B. cinerea population has heritable effects on lesion area and camalexin abundance, we used a general linear model (GLM), which assumes a Gaussian distribution, to estimate broad sense heritability (H2) as well as model-corrected means (LSMeans) and se (Supplemental Data Set 2). The residuals from these models adhere to Gaussian distribution confirming the a priori assumptions. Heritability due to variation across the B. cinerea isolates measures the fraction of lesion area variation that can be contributed to natural genetic variation in this B. cinerea population. The Arabidopsis genotype-specific heritability is the phenotypic variation solely controlled by the differences among Arabidopsis genotypes (i.e., wild-type Col-0, JA-insensitive mutant coi1-1, and SA-insensitive mutant npr1-1). By contrast, the heritability attributed to the interaction between pathogen isolate and host genotype illustrates the degree to which plant genotypes respond differentially among the B. cinerea isolates in a nonadditive fashion. Natural genetic variation in the pathogen population and variation among the immuno-compromised Arabidopsis defense mutants significantly influenced lesion area to a similar extent, as they account for the same proportion of variance (H2Isolate = 0.164, H2HostGenotype = 0.151, respectively; Table 1, Figure 1A). There was also a significant interaction between pathogen and plant genotypes that contributed about one third of the heritability of the host or pathogen alone, showing that the pathogen’s genetic variation differentially interacts with the host defense pathways (H2Isolate x HostGenotype = 0.047; Table 1, Figure 1A). To explore how isolates differentially interact with the host defense pathways, we used model-corrected means and standard errors for each isolate in downstream analyses.

View this table:
  • View inline
  • View popup
Table 1. F-Tables of the GLM for Lesion Area and Induced Camalexin Accumulation
Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Distribution of Genetic Variation Controlling Defense Responses.

Violin plots illustrating the distribution of estimated broad-sense heritability (H2) values for transcripts responding to 96 B. cinerea isolates across Arabidopsis genotypes. Heritability is partitioned across the different sources, 96 pathogen genotypes (Isolate), the Col-0, coi1-1, and npr1-1 plant genotypes (Host), and the corresponding interaction. The red lines indicate heritability values for lesion area, while blue shows those for camalexin accumulation in the same experiments. Plant defense-related phenotypes of camalexin accumulation and lesion area were measured on leaves from three Arabidopsis genotypes at 72 HPI with 96 diverse B. cinerea isolates. The transcriptomic analysis was conducted by sequencing mRNA extracted from infected Arabidopsis leaves at 16 HPI.

(A) Heritability distributions for all 23,898 detected Arabidopsis transcripts in B. cinerea infected plant tissues.

(B) Heritability values for Arabidopsis transcripts in B. cinerea infected plant tissues estimated by individually modeling the Col-0, coi1-1, and npr1-1 genotypes. The Arabidopsis genotypes are presented on the horizontal axis.

Using the entire population of 96 B. cinerea isolates to compare the Arabidopsis genotypes’ effect on B. cinerea showed that both the coi1-1 and npr1-1 mutations have a significant effect on lesion area compared with wild-type Col-0 (Table 1, Figure 2; an interactive plot, see Supplemental Data Set 3). The dramatic increase in lesion area of coi1-1 agrees with previous work showing that JA signaling is important for Arabidopsis defense against necrotrophs (Figure 2B) (Thomma et al., 1998, 1999b; Denby et al., 2004; Ferrari et al., 2007; Rowe et al., 2010). Interestingly, the npr1-1 mutant in SA signaling also showed a significantly increased susceptibility to the genetically diverse B. cinerea population, albeit with a smaller lesion relative to coi1-1 (Figure 2B). This shows that the SA-mediated innate immune responses also play an important role in plant resistance to necrotroph pathogens. Additionally, the significant interaction of pathogen and host implied that the pathogen genotypes elicit different responses across the immune compromised mutants, which led us to search for alternative patterns using hierarchical clustering (Figure 2A; for an interactive plot, see Supplemental Data Set 3). One major pattern is represented by the common B05.10 model isolate, where virulence is similar on the wild-type Col-0 and npr1-1 genotypes, with the coi1-1 mutant having increased susceptibility (Figure 2F) (Rowe et al., 2010). A contrasting pattern is illustrated by isolate 1.03.22 that maintains the enhanced virulence on coi1-1, but also displayed enhanced virulence on npr1-1 (Figure 2C). A third pattern as shown by isolate Apple 404 supports the theoretical antagonism, where there is increased virulence on the JA signaling mutant, coi1-1, and decreased virulence on the SA signaling mutant, npr1-1 (Figure 2D). A final pattern as shown by 2.04.11 is one in which the coi1-1 and npr1-1 mutants did not have any effect on the interaction (Figure 2E). We were unable to identify broader patterns in the isolates as there were no significant associations with any metadata related to these isolates, such as geography or host origin. Thus, this shows that there is variation within this population in how B. cinerea interacts with the Arabidopsis wild-type Col-0, npr1-1, and coi1-1.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Variation in Arabidopsis Susceptibility across Host Genotypes Driven by Natural Genetic Variation in B. cinerea.

Model-corrected lesion area means were estimated using the linear model with all the data from three Arabidopsis genotypes at 72 HPI with 96 B. cinerea isolates. The three Arabidopsis genotypes are either labeled at the bottom of the figure in (A) or shown with different color bars in (B) to (F), with red for wild-type Col-0, blue for the JA-insensitive mutant coi1-1, and brown for the SA mutant npr1-1. Error bars represent se for eight total samples across two independent experiments (n = 8).

(A) Hierarchical clustering of lesion areas produced by 96 B. cinerea isolates across Arabidopsis genotypes. Color key indicates the lesion area value. Red boxes show the position of isolates highlighted in (C) to (E).

(B) The average lesion area produced by 96 B. cinerea isolates within each Arabidopsis genotype.

(C) Lesion area of isolate 1.03.22.

(D) Lesion area of isolate Apple 404.

(E) Lesion area of isolate 2.04.11.

(F) Lesion area of isolate B05.10. Different letters above each bar show significant differences as determined by post-hoc Tukey’s HSD test following linear modeling.

Variation in Camalexin Accumulation in Response to B. cinerea Genetic Variation

To compare lesion formation to plant biochemical defense responses to genetic variation in B. cinerea, we measured the induction of camalexin accumulation across the three host genotypes at 72 HPI. Camalexin is a known phytoalexin in Arabidopsis that provides defense against B. cinerea depending upon the ability of the fungus to detoxify camalexin (Zhou et al., 1999; Stefanato et al., 2009; Rowe et al., 2010). Genetic variation in the host and pathogen significantly affected camalexin accumulation (Table 1; Supplemental Figure 1 and for an interactive plot, see Supplemental Data Set 4). However, unlike lesion area, there was no significant evidence for an interaction between host and pathogen genetic variation. Compared with heritability for lesion area, camalexin displayed a lower pathogen-dependent heritability (H2Isolate = 0.043), but with a similar plant-dependent heritability (H2HostGenotype = 0.152) (Table 1, Figure 1A). This indicates that variation in camalexin accumulation across this experiment is more dependent on differences between the wild-type Col-0, npr1-1, and coi1-1 plants and differences between the isolates than on their nonadditive interaction (Figure 3A).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Correlation between Camalexin Accumulation and Lesion Area across Arabidopsis Genotypes.

The relationship of camalexin and lesion area measured on leaves from three Arabidopsis genotypes at 72 HPI with 96 B. cinerea isolates was compared using model-corrected means. The three Arabidopsis genotypes are wild-type Col-0 (red dot), JA-insensitive mutant coi1-1 (blue triangle), and SA mutant npr1-1 (brown diamond).

(A) Scatterplots are shown comparing camalexin accumulation and lesion area across the 96 isolates within each of the three Arabidopsis genotypes. The 90% confidence ellipse intervals are drawn for each Arabidopsis genotype for reference. The regressions with the highest correlation are shown, with the equations being as follows: coi1-1, y = −2.1x + 40.7, P = 0.013; Col-0, y = −0.06x2 + 0.88x + 14.4, P = 0.014; npr1-1, y = −0.16x2 + 2.2x + 14.6, P = 0.012.

(B) Scatterplots are shown comparing lesion area (top left) and camalexin accumulation (bottom right) across pairwise Arabidopsis genotypes.

Camalexin Has Differential Effectiveness across B. cinerea Genetic Variation

Previous studies using one or a few B. cinerea isolates to focus on genetic difference in the Arabidopsis host have shown a negative correlation between camalexin content and lesion area across host genotypes (Ferrari et al., 2003; Denby et al., 2004; Corwin et al., 2016a, 2016b). This larger population of B. cinerea isolates allows us to test how genetic diversity in the pathogen influences the effectiveness of plant molecular defenses, specifically camalexin. We used a linear model to test correlations between camalexin and lesion area across the isolates in the three Arabidopsis genotypes, which showed a significant negative correlation in coi1-1 (r = −0.28, P valuecoi1-1 = 0.01, n = 96), but no significant linear correlation in wild-type Col-0 and npr1-1 (P valueCol-0 = 0.63, n = 96; P valuenpr1-1 = 0.90, n = 96, respectively) (Figure 3A; Supplemental Table 1). A visual analysis of the camalexin × lesion area plot for Col-0 infected leaves suggested the potential for a curvilinear relationship, which was confirmed by fitting a significant second-degree polynomial model between the two phenotypes in wild-type Col-0 and npr1-1 (P value ≤ 0.01 for both) (Figure 3A; Supplemental Table 1). There was no significant second-degree polynomial correlation between camalexin and lesion area in coi1-1. The curvilinear relationship in Col-0 and npr1-1 is a positive correlation within lower levels of camalexin accumulation, suggesting that camalexin is responding to increasing pathogen virulence and is not a major contributor to resistance in Col-0 and npr1-1 for those isolates. As camalexin accumulates beyond the inflection point of the curve (∼8 ng/mm), the relationship turns negative suggesting that camalexin is beginning to contribute to host resistance (Figure 3A). Interestingly, these different shaped correlations between camalexin and lesion size across the host genotypes was not reflected in different correlations of lesion across the host genotypes (Figure 3B). However, it should be noted that the residuals are not evenly distributed in these correlations suggesting a higher-order structure to lesion size (Figure 3B). This contradiction in relationships between camalexin and lesion area between the genotypes indicates that while JA-dependent signaling is important for camalexin defense, there are other signals able to induce camalexin in the absence of JA. Thus, the pathogen’s genetic variation differentially affects lesion area and camalexin accumulation across the host genotypes.

Plasticity in the Arabidopsis Transcriptomic Response to B. cinerea Genetic Variation

The camalexin analysis suggested that Arabidopsis biochemical responses to B. cinerea natural genetic variation may highlight different aspects of the interaction than are observable with lesion area. To conduct a broader survey of the molecular basis of the plant defense response against the B. cinerea invasion, we performed transcriptomic analysis by sequencing mRNA extracted from infected Arabidopsis leaves at 16 HPI. This time point was chosen as there should be minimal lesion outgrowth from the initial infection site and represents a significant switch in the Arabidopsis transcriptome in response to B. cinerea infection (Windram et al., 2012). The transcriptome of each Arabidopsis genotype × B. cinerea isolate combination was measured in fourfold replication across two independent experiments, including uninfected control (mock-treated) samples, for each Arabidopsis genotype. Specifically, there are two independent samples for each Arabidopsis genotype × B. cinerea isolate in each of two independent experiments providing four samples total. Each sample was made into an independent library and sequenced as a pool of 96 providing on average 1,172,000 ± 26,000 mappable Arabidopsis reads per library, or ∼4,500,000 reads per Arabidopsis genotype × B. cinerea isolate interaction. Using a negative binomial GLM, we tested for significance of all factors on each Arabidopsis transcript and obtained the proportion of variation (η2) due to the pathogen genotypes and plant mutants, least square means (log2), and se for transcript abundance (Table 2; Supplemental Data Sets 5 to 7).

View this table:
  • View inline
  • View popup
Table 2. Hub Genes and Bottlenecks Identified from Four Gene Coexpression Networks across Arabidopsis Genotypes

Of the 23,898 Arabidopsis genes that were quantifiable across infected samples across all three host genotypes, 20,328 (85%) transcript abundances responded differentially among the 96 diverse B. cinerea genotypes (false discovery rate [FDR]-adjusted P value < 0.001; Supplemental Table 2), indicating that diverse B. cinerea isolates triggered large-scale genome-wide differences in host gene expression in a pathogen-genotype-dependent manner during infection. In a similar fashion, the host genetic variation between the three plant genotypes significantly altered the expression of 88% of Arabidopsis genes within this experiment. Importantly, 34% of the Arabidopsis host genes showed significantly altered expression due to the interaction of pathogen and host genotypes (FDR-adjusted P value < 0.001; Supplemental Table 2 and Supplemental Data Set 7), which further demonstrates the differential interaction of natural variation within B. cinerea to these defense hormone pathways. Thus, the Arabidopsis genotypes and the diversity among B. cinerea isolates leads to a dramatic reshaping of the plant’s transcriptomic defense response.

To quantify the transcriptomic impact of genetic variation in the host and pathogen, we estimated the proportion of variance (η2) of each Arabidopsis transcript due to variation among the host genotypes (wild-type Col-0, coi1-1, and npr1-1), genetic variation in the pathogen (96 diverse isolates), and the interaction of host x pathogen (Figure 1). The average Arabidopsis transcript was equally affected by genetic variation in the pathogen and the interaction of pathogen × host (H2Isolate x HostGenotype = 0.111, H2Isolate = 0.108) (Figure 1). By contrast, the effect of the host’s variation on the average transcript was lower (H2HostGenotype = 0.05). This contrasts with both lesion size and camalexin, where the interaction of host and pathogen had the smallest effect. Thus, the interaction between pathogen and plant genetic variation significantly influences the host’s transcriptional responses when being attacked by diverse B. cinerea isolates. This indicates that the plant transcriptome is highly sensitive to genetic variation in the pathogen and how natural variation in the pathogen influences the SA and JA pathways in the plant (Figure 1A).

Plasticity of Known Arabidopsis Defense Pathways and Genes

To visualize the response of well-characterized plant defense pathways across the diverse B. cinerea genotypes in wild-type Col-0, coi1-1, and npr1-1 Arabidopsis genotypes, we examined transcript levels of individual genes in these pathways. We initially focused on the transcriptional response of specific pathways associated with camalexin production and defense signaling (Supplemental Data Set 7). In wild-type Col-0 and coi1-1 genotypes, all B. cinerea isolates gave rise to increased camalexin related transcript abundance for these genes, while some isolates lead to significant decrease from the mock level in the npr1-1 genotype (Figure 4A; Supplemental Figure 2D and Supplemental Data Set 7) (Zhou et al., 1999; Schuhegger et al., 2006; Liu et al., 2010).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Plasticity in Arabidopsis Defense-Related Genes in Response to B. cinerea Genetic Variation.

The Arabidopsis genotypes, wild-type Col-0, JA-insensitive mutant coi1-1, and SA mutant npr1-1, are shown on the x axis. Violin plots show the distribution of transcript accumulation in response to the 96 B. cinerea isolates for specific transcripts. Red lines indicate the average transcript accumulation in mock-treated plant tissues for each Arabidopsis genotype. The transcripts shown are PAD3 (A), AOS (B), PDF1.2b (C), PR1(D), WRKY54 (E), ycf2-A (F), PNSB2 (G), and LHCA2 (H).

We further compared the individual transcriptional performance of genes involved in JA and SA biosynthesis, perception, and downstream responses. The transcript abundance for most tested genes associated with JA synthesis and signaling was induced in wild-type Col-0 in response to infection by all B. cinerea isolates and showed opposing responses in the coi1-1 and npr1-1 mutants (Figure 4; Supplemental Figure 2). However, JASMONATE RESISTANT1 (JAR1), showed both activated and repressed transcriptional responses depending upon the specific B. cinerea isolate (Supplemental Figure 2C) (Staswick et al., 2002). More critically, there were isolates that had the ability to induce JA response genes, such as PLANT DEFENSIN1.2 (PDF1.2), VEGETATIVE STORAGE PROTEIN1 (VSP1), and VSP2 in the JA signaling deficient coi1-1 background (Figure 4C) (Penninckx et al., 1996; Ellis and Turner, 2001).

Most SA-responsive genes showed transcriptional responses that differed across the diverse B. cinerea isolates and were opposite of the above JA responses. They showed the highest transcript accumulation in the coi1-1 mutant, moderate in Col-0, and downregulation in the npr1-1 mutant (Figures 4E and 5; Supplemental Figure 2E). A key exception to this pattern is ENHANCED DISEASE SUSCEPTIBILITY5 (EDS5), the multidrug and toxin extrusion transporter that exports SA from the plastid, which showed a pattern that was largely independent of plant genotype and dominated by pathogen genetic variation (Supplemental Figure 2G) (Nawrath et al., 2002; Ferrari et al., 2003). Importantly, there were isolates that could shift this pattern such as Noble Rot that can actually induce transcripts like PR1 in an npr1-1 mutant background (Figure 5). This suggests that these isolates contain signals that can bypass NPR1 dependency for PR1 and other SA-responsive gene regulation (Figure 5). It remains to be characterized what these other signals may involve. Thus, the SA and JA related transcripts show a similar pattern of expression that is a mix of effects from the host and pathogen genetic variation.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

PR1 Expression across Arabidopsis Genotypes.

Rank plot showing the relationship of PR1 expression across three host genotypes (x axis) in control and in response to 30 diverse B. cinerea isolates (right). The isolates shown on the plot were chosen to provide an image of the major patterns found in the full collection of 96 isolates. The isolate names are shown in the same order as the lines for the npr1-1 mutant. The three Arabidopsis genotypes are wild-type Col-0 (red dot), JA-insensitive mutant coi1-1 (blue triangle), and SA mutant npr1-1 (brown diamond). The model-corrected means for transcript of PR1 are utilized for this plot. The transcript expression levels of PR1 across three Arabidopsis genotypes in control or induced by the same isolate are represented with colored connecting lines. Blue line indicates the expression levels of PR1 in the control treated samples. Purple lines indicate isolates where the expression levels of PR1 are higher in coi1-1 and npr1-1 than in Col-0. Black lines indicate isolates where the PR1 expression levels are higher in coi1-1 than in Col-0 but lower in npr1-1. Orange lines indicate isolates where the highest PR1 expression levels are in Col-0. The green line indicates an isolate with slightly higher PR1 expression levels in npr1-1 and slightly lower in coi1-1.

We finally investigated transcriptomic performance of genes in the Arabidopsis basal defense and plant innate immune system (Supplemental Figure 3 and Supplemental Data Set 5). Most transcripts associated with PTI were positively, albeit differentially, regulated in response to infection by diverse B. cinerea isolates on the three Arabidopsis genotypes as found for SA- and JA-responsive genes (Supplemental Figure 3). A key outlier is ACCELERATED CELL DEATH11 (ACD11), which was strongly suppressed in response to all 96 B. cinerea isolates in coi1-1 while showing no response to B. cinerea infection on wild-type Col-0 and npr1-1 (Supplemental Figure 3H). Interestingly, ACD11 suppresses the hypersensitive response and cell death, suggesting that in the coi1-1 background, this suppression of ACD11 function could enhance the hypersensitive response and sensitivity to B. cinerea (Brodersen et al., 2005). Thus, the use of a population of B. cinerea and plant signaling mutants illustrates that there is a wide diversity of host transcriptional networks involved in the response to B. cinerea infection. Additionally, the npr1-1 and coi1-1 mutants do not completely abolish host defense responses to the pathogen, but instead reprogram this response and some isolates can bypass the need for COI1 in JA signaling (Rowe et al., 2010).

Individual Gene Association

To test if any individual transcripts are highly correlated with lesion area across all the Arabidopsis genotypes and B. cinerea isolates, we utilized a Spearman’s rank correlation with all of the data. This showed that in any individual Arabidopsis genotype, no transcript had a correlation above r = 0.49 and using all the genotypes led to the maximum r = 0.52 (Supplemental Data Set 8). Thus, it does not appear that there is any singular gene or subset of genes that are highly correlated with resistance and lesion area. This supports the perspective that host resistance to B. cinerea is extremely polygenic (Corwin et al., 2016a). Interestingly, three of the top 10 highest positively correlated genes encode a chitinase, a chitin-inducible gene (AT5G25260), and a PAMP-inducible kinase FRK1 (Passarinho et al., 2001; Millet et al., 2010; Yeh et al., 2016; Nie et al., 2017). A fourth gene in the top 10 is UGT76B1, which controls crosstalk between the SA and JA pathways (von Saint Paul et al., 2011). In support of these observations was the finding that the top 100 correlated genes are enriched in biotic response terms when conducting a Gene Ontology (GO) enrichment analysis. At the other end of the spectrum, the gene with the highest negative correlation to lesion area was MUSE7, which is a conserved kinase substrate that influences NLR protein regulation (Johnson et al., 2017). The only enriched GO terms in the genes with a strong negative correlation to lesion size relate to photosynthesis. Thus, while the correlation of gene expression for individual genes with lesion area is low, genes that are the most positively correlated often appear to have functions related to pathogen resistance.

Variation of Arabidopsis Gene Coexpression Networks

To identify collections of genes in Arabidopsis that have a common transcriptional response to the genotypic variation in B. cinerea, we used a gene coexpression network approach (Zhang and Horvath, 2005; Dong and Horvath, 2007; Kliebenstein, 2009). In contrast to previous studies that used a single genetic isolate of the pathogen and different time points as source of variation, this analysis is based on natural genetic variation in the pathogen, which allows for a new perspective (Windram et al., 2012; Howard et al., 2013; Lewis et al., 2015). We first calculated the nonparametric Spearman’s rank correlation coefficient of model-corrected transcript means for all gene pairs within each Arabidopsis genotype and filtered for highly correlated gene pairs (R > 0.95). Using correlation thresholds at 0.9, 0.8, and 0.7 produced comparable network structures, albeit with greatly increased gene membership. Analysis of the pairwise correlation networks identified a set of 131 genes that were commonly coexpressed among the Arabidopsis genotypes. These 131 genes coalesced into four coexpression networks that had five or more genes and were subsequently named Networks I, II, III, and IV, respectively (Figure 6; Supplemental Data Sets 9 to 11 and Supplemental Figure 4). Using these genes as kernels, we extended these four core response networks within each genotype to include all coexpressed genes within each Arabidopsis genotype (R > 0.95). This showed that while the four networks contain a similar core of genes, their network architecture changes between the Arabidopsis wild type and the immune-deficient mutants. For example, the total gene membership (nodes) in Network I was 123 in wild-type Col-0, 34 in coi1-1, and 403 in npr1-1 (Figure 6; Supplemental Data Sets 9 to 11 and Supplemental Figure 4). The dramatic change in network membership across the Arabidopsis mutant backgrounds demonstrates the ability of plant defense hormones to shape coordinated defense responses across the diverse B. cinerea genotypes.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Gene Coexpression Networks Associated with Variation in Arabidopsis Transcriptomic Responses to Natural Genetic Variation in B. cinerea.

The four core transcript networks as estimated both by PCA and correlation are shown with the suggested biological function determined by GO enrichment. These are as follows: Network I, JA and SA signaling processes and camalexin biosynthesis; Network III, defense and cell cycle; Networks II and IV, photosynthesis. Nodes within each network represent Arabidopsis transcripts. Purple nodes show transcripts encoded by the plastid genome. Red and yellow nodes represent biosynthetic genes for camalexin and tryptophan, respectively. Blue and orange nodes represent genes in the JA- and SA-signaling pathways, respectively. The interrelationship between the nodes within each network represents similarity between the expression levels of transcript. The similarity matrix is computed using Spearman’s rank correlation coefficient. The interrelationship between the genes within these networks is shown for the Arabidopsis genotypes wild-type Col-0 (A), JA-insensitive mutant coi1-1 (B), and SA-insensitive mutant npr1-1 (C).

We tested for two types of nodes within each network, hub genes (highly connected nodes) and bottlenecks (high centrality hubs), to identify the potential biological function of each network and the potential drivers of interactions. Hub genes, genes with a high degree of connectivity, in a biological network are nodes with comparatively high levels of connectivity and are often assumed to be vital for network function (Yu et al., 2007). Using both connectivity and centrality scores, we identified the top 10% of hub and bottleneck genes in each core network within each genotype (Table 2). This showed that the hub and bottleneck genes for Network I were predominantly biosynthetic enzymes, such as PHYTOALEXIN DEFICIENT3 (PAD3) and TRYPTOPHAN SYNTHASE ALPHA1 (TSA1), required for the production of camalexin and its precursor tryptophan (Supplemental Data Set 12). This agreed with GO analysis and KEGG analysis that Network I is dominated by JA and SA signaling processes and indole biosynthesis (Supplemental Data Set 13).

Networks II and IV were both associated with chloroplast function (Figure 6; Supplemental Data Sets 9 to 12). Network II consists entirely of genes encoded within the plastid genome, indicating its biological function might be associated with defense-associated roles for plastids (Supplemental Data Set 12). It needs to be noted that as the methodology utilized poly(A) purification, the plastid genes are dependent on poly(A) stretches within the gene and as such, the quantification of these genes may have a differential and unknown bias than the nuclear genes. The only hub genes in this network are the YCF2.1 and YCF2.2 loci, whose function is presently unclear (Table 2). Similarly, genes in Network IV are highly enriched in nuclear-encoded photosynthetic genes that are centered upon the photosystem reaction centers located in the chloroplast (Supplemental Data Set 13). Together, these networks suggest that B. cinerea genetic variation directly affects the functioning of the photosynthetic apparatus within the Arabidopsis plastid during the pathogenic interaction. While both networks influence plastid function, they are differently regulated in response to genetic diversity among the B. cinerea isolates, as Network II’s interconnectivity decreases in an npr1-1 background but increases in a coi1-1 background in comparison to the wild type. By contrast, Network IV’s connectivity increases in both npr1-1 and coi1-1 (Figure 6; Supplemental Data Sets 9 to 11). The expression of genes of Network II was repressed in infected wild-type Col-0, and this repression was enhanced in npr1-1 (Supplemental Data Set 7). This suggests that JA functions to repress these plastid-encoded genes, while SA activates their expression. The network and GO analysis showed that Network III was associated with defense and cell cycle genes previously unlinked to B. cinerea interactions (Supplemental Data Set 13). Thus, the four core host transcript coexpression networks that coordinately respond to genetic variation in the pathogen are focused on networks that function in metabolic host defense responses and chloroplast gene regulation.

Arabidopsis Wild Type Transcriptomic Variation Can Identify Altered Resistance

To quantify the network’s transcriptional response, we utilized PCA to obtain weighted eigengenes to measure the network responses of Arabidopsis transcripts to variation across the 96 B. cinerea isolates within the wild-type Col-0 genome (Alter et al., 2000). An eigengene is a vector that allows the direct measurement of a network of genes as a single value. It builds from a set of weightings for each gene in the network. The eigengene value/score is calculated by multiplying the gene’s weight by the gene’s transcript level and summing this for all the genes in the network. The weightings provide some indication of how the gene connects to the network. This showed that the first four principal components corresponded to Networks I to IV, respectively, i.e., PC1 is an eigenvector describing the variance in expression of the genes in Network I (Figure 7; Supplemental Figure 5). Comparing the PCA to the coexpression networks showed that genes within Network I largely contribute to PC1, while Network II largely contributes to PC2, suggesting that the principle components in PC-ordination space are being driven in a network dependent manner (Supplemental Figure 5 and Supplemental Data Set 12). We then utilized the loading vector for each eigengene to calculate the corresponding value for the control samples and the B. cinerea infected coi1-1 and npr1-1 transcriptomes (Figure 7). Network 1/PC1 associates with 63% of the transcriptome variation in Col-0 and the projected eigengenes for this vector link to 53% and 83% of the transcriptome variation in coi1-1 and npr1-1, respectively. Network 2/PC2 explained 5% of the Col-0 variation and 3% of the variation in both coi1-1 and npr1-1. Interestingly, the uninfected control samples are not extreme outliers within the ordination field among the first two components as they reside just barely outside the 90% confidence interval (Figure 7). Thus, the transcriptomes for the wild-type Col-0, coi1-1, and npr1-1 genotypes are responding to the pathogen’s genetic diversity within a transcriptomic space that is not being significantly expanded, or altered, by the loss of either JA or SA signaling. This further suggests that the coi1-1 and npr1-1 mutants do not abolish the ability to respond to B. cinerea, but instead help transform and shape the transcriptional response.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

PCA Comparison of Arabidopsis Transcriptomic Profiles Responding to Natural Genetic Variation in B. cinerea.

Scatterplots show the scores of the first two principal components estimated using the response of the top 2000 transcripts responding to B. cinerea within the total data set using the mean response in Col-0. Red dots show the scores for the response to the individual isolates in Arabidopsis wild-type Col-0. Blue triangles show the response of coi1-1, and brown diamonds show the response of the SA-insensitive mutant npr1-1. Red, blue, and brown crosses represent the imputed score from mock treated samples from Col-0, coi1-1, and npr1-1, respectively. The 90% confidence ellipse intervals are drawn on scores of the first two principal components from each Arabidopsis genotype for reference.

Using the eigenvector-derived weighted values for the different networks, we tested whether the transcriptomic responses at an early infection time point (16 HPI) were linked with later lesion development and camalexin induction (72 HPI) (Supplemental Table 3). Linear modeling with all three host genotypes showed that the early transcriptome response of Network I at 16 HPI is highly linked to resistance at 72 HPI. This relationship was significantly different between the host genotypes, with wild-type Col-0 having a negative correlation and npr1-1 having a positive correlation between PC1 of the 16 HPI transcriptional response and the 72 HPI lesion area (Figure 8; Supplemental Tables 1 and 3). Including the eigengene for Network II/PC2 in the model as an interaction term with Network I/PC1 improved the model’s ability to link the transcriptomic response to lesion size and camalexin accumulation, suggesting that resistance is connected to an interaction of these two networks as well as other networks (Supplemental Tables 1 and 3). Thus, by combining information across the transcripts using eigengenes, we can develop a statistical model that describes lesion size variation at 72 HPI using transcriptomic variation measured at 16 HPI.

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

Correlation between Lesion Area and Transcriptome Responses of Arabidopsis Genotypes.

Scatterplots are shown comparing lesion area as measured at 72 HPI to the transcriptomic response at 16 HPI as measured using the first two principal components. The three Arabidopsis genotypes are wild-type Col-0 (red dot), JA-insensitive mutant coi1-1 (blue triangle), and SA mutant npr1-1 (brown diamond). Each point represents the model-corrected mean for each of the 96 isolates. The 90% confidence ellipse intervals are drawn for reference lesion area versus transcriptome PC1 (A) and lesion area versus transcriptome PC2 (B).

Genetic Variation in B. cinerea Identifies Different Canalization and Plasticity Roles for COI1 and NPR1

The significant interaction of the Arabidopsis defense hormone perception with the B. cinerea isolates suggests that COI1 and NPR1 control the response of Arabidopsis when faced with a diversity of pathogen virulence mechanisms via one of two options (Weinig, 2000; Liefting et al., 2009; Palacio-López et al., 2015). First, COI1/NPR1 may canalize the plant’s response by channeling the response to a diversity of stimuli/attacks into a single focused, coherent outcome. Alternatively, COI1/NPR1 may control the plasticity of the plant’s response to this diverse mixture of attacks without canalizing the response. These alternatives can be visualized using a rank order plot. For example, if the isolates are plotted in their rank across the three Arabidopsis genotypes, NPR1/COI1-mediated canalization would generate an image where the rank of isolates is the same (i.e., no crossing lines), but with an increase or decrease in the range of the lesion area. By contrast, plasticity within a rank order plot would show the same general range of isolates between the Arabidopsis genotypes, but the rank of the isolates would dramatically change, leading to crossing lines. Present theory about SA and JA simply argues that the plant hormones control the response to pathogens, but it is unknown whether they control canalization and/or plasticity in response to diverse isolates (Figure 1). To assess this possibility, we visualized the rank order of diverse B. cinerea isolates, to see if this fluctuates among the three Arabidopsis genotypes depending on the defense phenotypes, from resistance to the transcriptional network (Figure 9). Rank plots of lesion area showed that the significant isolate × host genotype effect is a result of a lack of canalization/buffering in the npr1-1 and coi1-1 genotypes. The range of lesion sizes is larger in these two mutant host genotypes, but there is little change in the isolate’s relative rank in comparison to the wild-type Col-0 host, suggesting that COI1/NPR1 largely canalize host resistance (Figure 9A). Comparing camalexin accumulation to lesion area, the rank order of the isolates is largely fluctuating, especially in the npr1-1 background, suggesting that the defense hormones contribute to plasticity but not canalization in camalexin production (Figure 9B). Proceeding to the transcriptome-derived eigenvectors for Network I/PC1 and Network II/PC2 showed that there was no evidence of coi1-1 or npr1-1 canalizing the transcriptomic response, as the range of the response is similar between all three host genotypes. Instead, the coi1-1 and npr1-1 mutants illuminate a highly plastic transcriptome where there are dramatic rank order shifts in how the plant’s transcriptome responds to the genetic diversity of the pathogen (Figures 9C and 9D). Interestingly, npr1-1 seems to function to invert the rank order of the wild-type Col-0 responses. Thus, the transcriptional response to the diverse B. cinerea isolates across the Arabidopsis hosts shows greater isolate × host genotype interactions than the resulting lesion size. This suggests that the JA and SA pathways function to canalize the final resistance phenotypes (lesion area) while not similarly buffering the initial transcriptional responses.

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9.

Differential Interactions of Host × Pathogen Genotype across Host Defense Responses.

Rank plots showing the relationship of Arabidopsis defense responses across three host genotypes (x axis) in response to 96 diverse B. cinerea isolates. The three Arabidopsis genotypes are wild-type Col-0 (red dot), JA-insensitive mutant coi1-1 (blue triangle), and SA mutant npr1-1 (brown diamond). The model-corrected means for all traits are utilized for these plots. The connecting lines are colored in a gradient from green to purple based on their value in the Col-0 host genotype.

(A) Lesion area at 72 HPI.

(B) Camalexin accumulation at 72 HPI.

(C) Transcriptome PC1 at 16 HPI.

(D) Transcriptome PC2 at 16 HPI.

DISCUSSION

Using a large collection of natural B. cinerea isolates, we assessed how Arabidopsis responds to quantitative genetic diversity in a generalist necrotrophic pathogen. This showed extensive genetic variation among pathogen genotypes for lesion formation, camalexin accumulation, and host-transcriptomic responses that were dependent upon the interaction of the host innate immune system and pathogen genotype. The host and pathogen genotypes had similar effects, which was unexpected as we utilized two major effect Arabidopsis mutants that largely eliminate JA (coi1-1) and SA (npr1-1) downstream signaling, which should overwhelm the phenotypic variance among B. cinerea genotypes (Xie et al., 1998; Xu et al., 2002; Mou et al., 2003; Sheard et al., 2010; Wu et al., 2012). While most models assume that these mutants abolish responses, the varied response to the diverse isolates indicates that while these mutants may diminish their specific response (i.e., SA or JA dependent), they lead to new and increased responses in other networks. This indicates that coi1-1 does not abolish the plant’s ability to detect and respond to B. cinerea but simply alters the overall defense response. It should be noted that both coi1-1 and npr1-1 are pleiotropic mutations that have significant effects on all downstream processes, which includes crosstalk between the SA and JA pathways. As such, we are interpreting these data concerning the genetic effect of removing these key nodes, rather than specific mechanistic molecular insight.

In contrast to the common model where SA and JA function to induce or repress specific plant responses, the lesion area measurements show that a key function of the SA and JA pathways is to canalize the overall defense response to the diverse B. cinerea isolates. However, this canalization function did not similarly influence the transcriptional response to these same B. cinerea isolates. Instead, the SA and JA signaling mutants create dramatically plastic transcriptomic responses to the pathogen’s genetic variation that occur within the same range as the wild-type Col-0 (Figure 9). This suggests that the SA and JA pathways canalize lesion formation, not by creating a single coordinated transcriptomic response, but by choosing between different molecular responses to the diverse isolates. As such, the wild-type Col-0 genotype gives the appearance of canalized responses to the diverse B. cinerea isolates that is achieved by the downstream signaling pathways picking and choosing among a myriad of response modules, some with positive and some with negative effects on the interaction, that together sum up to a net lesion size. This summation of positive and negative effects agrees with the observation that the coi1-1 and npr1-1 mutants both lead to decreased resistance to the pathogen population even though they are widely believed to be antagonistic.

Natural Genetic Variation in B. cinerea Influences Plant Transcriptomic Response

The interaction of plants with pathogens involves the passage of molecular signals back and forth between the two organisms. With specialist pathogens, this is frequently highlighted by the presence of single large-effect qualitative virulence genes (Glazebrook, 2005). By contrast, the Arabidopsis response to the 96 diverse B. cinerea isolates is quantitative and showed no evidence of bimodality or multimodality (Table 1, Figures 2 and 3A). This suggests that this collection of B. cinerea genotypes does not contain large-effect qualitative virulence loci; instead, this generalist uses a highly polygenic virulence architecture (Denby et al., 2004; Finkers et al., 2007b; Rowe and Kliebenstein, 2008; Corwin et al., 2016b; Zhang et al., 2016). This agrees with previous B. cinerea mechanistic studies that identify a large collection of standing genetic variation that controls a wide range of specific virulence mechanisms (Choquer et al., 2007; Rowe and Kliebenstein, 2007; Noda et al., 2010; Rowe et al., 2010; Dalmais et al., 2011; Michielse et al., 2011; Shlezinger et al., 2011; Windram et al., 2012; Pearson and Bailey, 2013; Kumari et al., 2014; An et al., 2015; Atwell et al., 2015; Hevia et al., 2015; Plaza et al., 2015; Schumacher et al., 2015; Zhang et al., 2015a; Corwin et al., 2016a, 2016b; Lopez-Cruz et al., 2017; Zhang et al., 2016). Taken together, this supports the perspective that the Arabidopsis-B. cinerea interaction involves a large polygenic pool of genetic variation in the host and pathogen that interact to determine the outcome (Figures 2 and 9A) (van Kan, 2006).

Relative Importance of Camalexin Mediated Defense

The dominant network with altered responses to the genetic variation in the pathogen isolates was Network I, which contains the majority of biosynthetic genes for camalexin and its precursor, tryptophan, along with key transcription factors known to regulate camalexin production, including MYB122, WRKY40, and ANAC042 (Pandey et al., 2010). The regulatory hierarchy of this biosynthetic pathway is highly flexible with different pathogen species or stresses frequently identifying different key transcription factors or regulatory processes (Glazebrook and Ausubel, 1994; Thomma et al., 1999a; Glazebrook, 2001; Mengiste et al., 2003; Ren et al., 2008; Mao et al., 2011; Birkenbihl et al., 2012; Saga et al., 2012; Frerigmann et al., 2015). It is possible that this regulatory flexibility is why the camalexin network is a major network responding to pathogen genetic diversity.

While camalexin controls resistance to some B. cinerea isolates depending upon their genotype, little is known about how this phytoalexin may influence the interaction with a broad population of isolates (Kliebenstein et al., 2005; Bednarek et al., 2011; Rajniak et al., 2015; Frerigmann et al., 2016). Comparing the transcript accumulation for all the enzymes in the biochemical pathway and the final defense metabolite showed that while camalexin accumulation is highly diminished in the coi1-1 background, the transcriptional response of the genes encoding the biosynthetic enzyme is largely intact in this JA mutant (Figure 4A; Supplemental Figure 1). Interestingly, the npr1-1 mutant also showed a general decrease in camalexin levels that positively correlated with a reduced transcriptional response. Thus, JA and SA are both positive regulators of camalexin accumulation in response to the population of B. cinerea isolates via transcriptional and posttranscriptional processes. Interestingly, expression of camalexin genes is more sensitive to genetic variation among the pathogen isolates than to the major mutants in the SA/JA perception pathways in the plant (Supplemental Data Set 7). Further use of the pathogen’s natural variation in stimulating signaling and biosynthesis of camalexin in Arabidopsis may provide new insights into plant-pathogen interactions that have not previously been studied.

JA/SA Molecular Antagonism versus Additivity

JA and SA signaling pathways are classified as the major plant defense hormones against necrotrophs and biotrophs, respectively (Vlot et al., 2009; Mahesh et al., 2012; Mengiste, 2012; Mulema and Denby, 2012; Windram et al., 2012; Yang et al., 2015). In response to the population of B. cinerea isolates, marker genes in JA and SA responses behaved largely as expected with JA-responsive genes downregulated or off in the coi1-1 background and vice versa for SA/npr1-1 (Figures 4 and 5; Supplemental Figure 2). Interestingly, this contrasts with our observation that the resistance to the average B. cinerea isolate is significantly increased in both npr1-1 and coi1-1. Thus, SA and JA signaling both provide resistance to B. cinerea and are not in antagonism at this level. Further complicating our ability to interpret the SA/JA signaling pathways is the fact that there were isolates that could elicit JA and SA signaling within coi1-1 and npr1-1 mutant backgrounds and bypass JA/SA perception. For example, there were B. cinerea isolates that could induce the JA marker PDF1.2 in the coi1-1 mutant, while some isolates could stimulate SA responses, such as PR1, in the npr1-1 background (Figures 4 and 5). As both npr1-1 and coi1-1 block perception of SA and JA respectively, this is not caused by fungus exogenously producing the hormone to attenuate host defenses. Instead, these individual isolates suggest that these classical JA- and SA-specific transcripts can be stimulated by alternative pathways (Sato et al., 2010; Tsuda et al., 2013; Kim et al., 2014). Thus, these specific isolates can provide unique tools to query previously unrecognized defense signaling pathways.

The Plastid Transcriptome and B. cinerea Interactions

Two of the four major plant transcriptional networks that were influenced by natural genetic variation in B. cinerea were directly associated with plastid function. Network II (PC2) solely contained genes encoded in the plastid genome with no associated nuclear locus (Figure 6). Most genes in Network II are involved in photosystem I, with the uncharacterized transcript YFC2 acting as a hub (Table 2, Figure 6). The expression and network structure of the plastid-encoded Network II was influenced by the npr1-1 and coi1-1 mutants, suggesting that the SA-JA signal systems that function within the nucleus and cytosol can lead to rapid shifts in the expression of plastid-encoded genes. Network IV contained nuclear-encoded chloroplast-localized genes involved in photosynthesis and reactive oxygen species production (Figure 6; Supplemental Data Sets 9 to 12) or in PC2 (Supplemental Data Set 14). Furthermore, these results are consistent with downregulation of genes associated with chloroplast organization after infection (Windram et al., 2012). While the importance of the plastid in plant defense is beginning to be highlighted, the roles of these plastid systems in altering interactions with B. cinerea have barely been investigated to date (Serrano et al., 2016). One possibility comes from the observation that JA signaling induces chlorosis surrounding the developing B. cinerea lesion in a pathogen genotype-dependent manner (Rowe et al., 2010; Corwin et al., 2016a). This chlorosis surrounding the lesion may play a role either directly or indirectly in defense. Alternatively, the plastids are major sites of amino acid synthesis that lead to the production of precursors for a number of plant defense metabolites and signals. These results indicate that there is a broader need to understand how the plastid responses are coordinated into the broader pathogen response.

Summary

We generated phenotypic and transcriptomic data in Arabidopsis wild-type Col-0 and JA- and SA-insensitive mutants challenged by 96 diverse B. cinerea isolates and applied a variety of computational and statistical tools to obtain a detailed comprehension of varied defense responses in multiple layers of the defense system in plants infected with necrotrophic pathogens. We identified four gene coexpression networks as well as eigengene vectors in Arabidopsis responsible for the defense response triggered by natural genetic variation in necrotrophic pathogens. An exciting future challenge will be the identification of integrating transcriptional network models from both pathogen and plant, which may provide more insights on the strategies employed by two organisms for attack and defense. Further biochemical and evolutionary analysis of genes predicted by the network models will drive the discovery of novel pathogen targets in plants and facilitate resistance crop breeding in sustainable agriculture.

METHODS

Plants Material and Growth Conditions

For all experiments listed, we used the Arabidopsis thaliana accession Col-0 and the well-characterized mutants coi1-1 and npr1-1, two mutant plant genotypes that abolish the major defense perception pathways of JA and SA, respectively. We grew the three genotypes over two independent experiments as part of a larger experiment. The plants were planted within 18 pot flats with one plant per pot using a randomized complete block design. Each flat contained three independent plants of each genotype with three other genotypes filling the remaining pots. Using a total of 30 flats per experiment, we grew 90 plants per genotype per experiment and two independent experiments were conducted separately. The Arabidopsis seeds were vernalized in 0.1% phytoagar at 4°C for 4 d in the dark. Three seeds of each genotype were placed in the center of a cell filled with soil (Sunshine Mix #1; Sun Gro Horticulture) and the flat covered with a humidity dome to maintain high humidity during germination. The humidity dome was removed one week after germination and the plants were thinned to one plant per cell. Plants were watered twice a week and grown in a growth chamber at 22°C under a photoperiod of 10 h light/14 h dark. High output fluorescent bulbs were set for an output of 150 µE. Five weeks after sowing, we chose the first fully mature adult leaf and the next five to six leaves in the developmental cycle as these leaves are highly similar and placed them on 1% phytoagar in large plastic flats identical to those used under the pots (agar flats) prior to infection with Botrytis cinerea. Under our growth conditions, flowering does not occur until week 9-10 and no senescence was observed on any of the leaves used.

Fungal Isolate Growing Conditions and Inoculations

For growth and inoculation of a pathogen, we utilized a long-established detached leaf assay that allows for high-throughput analysis and is typically consistent with whole plant assays (Govrin and Levine, 2000; Ferrari et al., 2003; Mengiste et al., 2003; Denby et al., 2004; Sharma et al., 2005; Mulema and Denby, 2012; Boydom, 2013). A total of 96 B. cinerea isolates were selected for their phenotypic and genotypic differences (Supplemental Data Set 1). These isolates are all natural isolates that were largely isolated in California from a diversity of crops. There is also a set of international isolates obtained from other labs (Supplemental Data Set 1). Previous work with gene-specific sequencing and genome sequencing of a subset of these isolates has shown that there is no significant population structure when comparing California to the international isolates, suggesting that they represent a sample of the global intermating population (Atwell et al., 2015; Corwin et al., 2016a).

A summary of the pathogen inoculation and growth is as follows. Small canned peach slices (∼3 cm3) in Petri plates were inoculated from frozen glycerol stocks of isolated spores and were incubated on the bench top for one week. Spores were collected by submerging the sporulating peach slices in 5 mL filter-sterilized ultrapure water and agitated with a flame-sterilized glass rod. The solution was filtered through a syringe containing a small plug of sterilized glass wool to remove hyphal and plant particles followed by a mild centrifugation in a bucket rotor at room temperature at 1000 rpm for 15 min. The supernatant was discarded and the spore pellet was resuspended in sterilized 0.5× organic grape juice (Santa Cruz Organics). Spore concentration was determined using a hemacytometer and spore solutions were diluted to 10 spores/μL using the 0.5× organic grape juice. Leaves detached from 5-week-old Col-0, coi1-1, and npr1-1 were inoculated with 4 μL of the diluted spore solution and were incubated at room temperature on a flat containing ∼2 cm of 1% phytoagar covered with a humidity dome. B. cinerea isolates were inoculated in a randomized complete block design across the six trays. The inoculations all occurred within an hour of dawn for the leaves and the light period of the leaves was maintained. Two of the six blocks were harvested at 16 HPI for RNA extraction by transferring the infected leaf into a 2-mL Eppendorf tube containing one 4-mm and four 2.4-mm stainless steel ball bearings and immediately submerging the tube in liquid nitrogen. Frozen tubes were stored at −80°C until RNA extraction. The infected leaves were first photographed for lesion analysis and then placed in 400 μL of 90% methanol for camalexin extraction at 72 HPI.

RNA-Seq Library Construction, Sequencing, and Alignment

RNA-seq library construction was mainly based on previous methods (Kumar et al., 2012). Frozen infected leaves were homogenized by rapid agitation in a bead beater followed by direct mRNA isolation and purification from tissue lysate using the DynaBeads Oligo(dT)25 Kit (Invitrogen). The first-strand synthesis of mRNA was created using the SuperScript III kit (Invitrogen) with random primers and RNaseOut. The second strand synthesis was accomplished using 50 units of DNA Pol I (Fermentas) with an equal mix of 2.5 mM dNTPs and 1.6 units of RNase H. The resulting cDNA was purified using 0.8× of AMPure beads and fragmented using the Fragmentase enzyme (New England Biolabs) for 20 min at 37°C. Fragmented cDNA was repaired using the End Repair Kit (New England Biolabs) and an A-base was added using the exo-Klenow fragment (New England Biolabs) with 0.8× AMPure bead cleanup steps between both reactions. A barcoded sequencing adapter (NEXTflex DNA Indexed Adapters; Bioo Scientific) was ligated to sample fragments using the NEB Quick Ligase ligation kit (New England Biolabs). Adapter-ligated fragments were size selected using 0.8× Ampure beads and PCR-enriched using the Phusion Master Mix (New England Biolabs) and 10 μM PE sequencing primers. PCR conditions ran for 98°C for 30 min followed by 14 cycles of 98°C for 30 s, 62°C for 30 s, 72°C for 30 s, and final extension at 72°C for 5 min. A total of 2 μL of the PCR-enriched samples was run on a 1.5% agarose gel to test for successful amplification. Successfully amplified samples were size-selected using 0.7× AMPure beads. Amplified and size selected libraries were pooled in groups of 96 libraries according to the unique indexed barcodes within the adapter and submitted for single-read, 50-bp sequencing on a single lane of Illumina HiSeq 2500 at the U.C. Davis Genome Center-DNA Technologies Core. A pool of 96 libraries represented a single replicate of one plant genotype across all isolates and was sequenced as a unit. As there were four independent samples for each isolate on each genotype, there were four pools of libraries for Col-0, four for coi1-1, and four for npr1-1. This is a total of 1164 libraries accounting for all RNA samples. After sequencing, 61 libraries were discarded due to poor quality. This left at least three independent samples for all host genotype x pathogen isolate combinations.

Fastq files from individual HiSeq lanes were separated by adapter index into individual RNA-seq library samples. Individual libraries were qualitatively assessed for overall read quality and overrepresented sequences using FastQC software (version 0.11.3; www.bioinformatics.babraham.ac.uk/projects/). To align the reads to the Arabidopsis TAIR10.25 cDNA reference genome, we used Bowtie 1 V.1.1.2 with phred33 quality scores (http://sourceforge.net/projects/bowtie-bio/files/bowtie/1.1.2/) (Langmead et al., 2009). The first 10 bp of the reads were trimmed by the fastx toolkit due to low quality at the beginning of reads (http://hannonlab.cshl.edu/fastx_toolkit/commandline.html). Gene counts were pulled from the resulting sam file using a combination of SAMtools (Li et al., 2009) and custom R scripts and then summed across gene models to reduce overrepresentation of genes with multiple splice variants.

Picture Analysis and Lesion Measurement

Pictures of 72 HPI leaves were analyzed using a custom semiautomated image analysis pipeline implemented in the R V3.2.1 statistical environment (R Core Team, 2011) to measure lesion traits as previously described (Corwin et al., 2016b). Briefly, leaves were identified as highly saturated objects with a green hue in HSV colorspace using the EBImage and CRImage packages in R (Failmezger et al., 2010; Pau et al., 2010). Lesions were identified as low saturation, brown objects within the leaf area, and independent leaf and lesion masks were generated. Lesion masks were refined individually by a technician due to a high variance in lesion appearance. Lesion areas and perimeters were then calculated using EBImage to calculate the number of pixels in the lesion object and convert to area using a ruler reference within each image to determine the number pixels per centimeter.

Camalexin Detection

Camalexin was extracted from mock treated and B. cinerea-infected leaves and measured via HPLC analysis as previously described (Kliebenstein et al., 2005). Briefly, leaves were homogenized in 400 μL of 90% methanol followed by centrifugation to remove precipitates. The supernatant was transferred to a fresh 96-well autosampler plate and 50 μL of the extract was run on an Agilent 1100 series HPLC equipped with an Agilent Lichrocart 250-4 RP18e 5-µm column. Camalexin was detected using a fluorescence detector with an excitation at 318 nm and emission at 385 nm. Separation of camalexin from the complex mixture was accomplished using an aqueous to acetonitrile gradient (5 min gradient from 63 to 69% acetonitrile, 30-s gradient from 69 to 99% acetonitrile, 2 min at 99% acetonitrile, and a postrun equilibration of 3.5 min at 63% acetonitrile). A standard curve from purified camalexin had been previously determined for this method and used to quantify camalexin. Camalexin content was normalized to the perimeter of the lesion (ng/mm) (Kliebenstein et al., 2005).

Statistical Analysis Methods

All the analyses were conducted using the R V3.2.1 statistical environment (R Core Team, 2011). To test for the ability of natural variation within B. cinerea to impact disease-related phenotypes, we ran the following generalized linear model with a Gaussian link function for camalexin and lesion area, and with a negative binomial link function for all transcripts:

camalexin/lesion area:

Embedded Image

transcripts in a single host genotype:

Embedded Image

transcripts among infected samples in host genotypes:

Embedded Image

transcripts among infected and uninfected samples in host genotypes:

Embedded Image

where the main effects E, N, I, and H are denoted as experiment, infection status, isolate genotype, and plant host genotype, respectively. Nested effects of the growing flat (Gf) within the experimental replicates and agar flat (Af) nested within growing flat are also accounted. For camalexin and lesion area, mock controls were zero for all samples tested and therefore were not included in the model and did not necessitate accounting for variation in infection status (N). However, infection status was included among the transcripts to account for the large amount of variation between infected and uninfected leaves. For host transcripts, gene count data obtained from RNA-seq experiments were subjected to analysis using multiple statistical tools. Normalization on gene counts was first conducted using the TMM method in function calcNormFactors() from the “edgeR” package (Robinson and Smyth, 2008; Robinson et al., 2010; Robinson and Oshlack, 2010; Nikolayeva and Robinson, 2014), and normalized pseudo-counts were obtained for downstream analysis. The linear model was conducted on normalized gene counts using function glm.nb() from the “MASS” package (Venables and Ripley, 2002). Model-corrected means and standard errors for camalexin, lesion area, and each transcript within each isolate were determined using the “lsmeans” V2.19 package (Lenth, 2016). Raw P values for the F-test and χ2 test were determined by Type II sums of squares using the ANOVA function from the “car” package (Fox and Weisberg, 2011). Transcript P values were FDR corrected (P value < 0.05) for multiple tests of significance (Benjamini et al., 2001; Strimmer, 2008). Broad-sense heritability of each phenotype and transcripts as the proportion of variance attributed to B. cinerea genotype, Arabidopsis genotype, or their interaction effect to the total variance within the model were estimated. We then compared broad-sense heritability of transcripts involved in the camalexin, tryptophan, JA, and SA biosynthesis pathways in each Arabidopsis genotype. The violin plots were generated to depict the distribution of heritability using the “vioplot” package (Hintze and Nelson, 1998).

Gene coexpression networks for all Arabidopsis transcripts were generated using the model-corrected means for each isolate. Spearman’s rank correlation coefficient of these model-corrected means of transcript in plant tissues derived from each infected Arabidopsis genotype were calculated for all gene pairs. Gene pairs with a positive Spearman’s correlation greater than 0.95 in each data set were selected as coexpressed genes and overlapped genes from three selected gene data sets were used to construct gene coexpression networks. Gene coexpression networks were visualized using Cytoscape V3.2.1 (Java version 1.8.0_60) (Shannon et al., 2003). Genes identified from each Arabidopsis genotype represented the nodes in the networks and were named based on AGI locus IDs as published by TAIR10 (Lamesch et al., 2012). The B. cinerea induced transcript profiles of genes involved in camalexin, tryptophan, SA, and JA pathways, as well as in plant immune systems, were illustrated by violin plots (Hintze and Nelson, 1998). Tukey’s multiple comparison was performed using the HSD test function from the “agricolae” R package (Mendiburu, 2016).

To investigate the associations among Arabidopsis wild-type Col-0, coi1-1, and npr1-1 challenged by 96 diverse B. cinerea isolates, we conducted a Mantel test on correlation matrices generated by model-corrected transcript means in each Arabidopsis genotype using a custom R script based on the function getPermuteMatrix() from the “vegan” V2.3.0 package (Mantel, 1967; Oksanen et al., 2016). The Mantel test was inferred following 999 permutations using Spearman’s Rank coefficient method.

We selected the top 2000 genes whose expression was significantly regulated by B. cinerea infection. PCA was performed on model-corrected means of each transcript in Arabidopsis wild-type Col-0 using the princomp function in R to capture the unobserved data structure (R Core Team, 2011). The first two principal component score values for two mutant data sets, as well as transcripts in the mock control samples, were imputed using wild-type Col-0 as training model. We then conducted PCA on transcripts of each Arabidopsis mutant to obtain the data structure. We employed an ANOVA model to identify the statistical significance of influence of host genotype, each principal component, or a specific interaction between host genotype and a given principal component. The ANOVA model was as follows:

Embedded Image

Heat map visualization used for clustering analysis of lesion and camalexin was conducted in R using the heatmap2 function from the “gplot” R package (Warnes et al., 2016). Gene annotation was obtained using https://www.araport.org/. GO enrichment analysis was performed with a BiNGO plug-in in the Cytoscape environment using Fisher’s extract test with multiple testing correction of FDR at 0.05 (Benjamini and Hochberg, 1995; Maere et al., 2005).

Accession Numbers

Sequence data from this article can be found in the GenBank/EMBL data libraries under the following accession numbers: ACD11, AT2G34690; ANAC042, AT2G43000; AOS, AT5G42650; chitinase, AT2G43570; COI1, AT2G39940; EDS5, AT4G39030; FRK1, AT2G19190; JAR1, AT2G46370; LHCA2, AT3G61470; MUSE7, AT5G46020; MYB122, AT1G74080; NPR1, AT1G64280; UGT76B1, AT3G11340; PAD3, AT3G26830; PDF1.2, AT5G44420; PNSB2, AT1G64770; PR1, AT2G14610; VSP1, AT5G24780; VSP2, AT5G24770; WRKY40, AT1G80840; and WRKY54, AT2G40750.

Supplemental Data

  • Supplemental Figure 1. Variation in camalexin accumulation across Arabidopsis genotypes induced by natural genetic variation in B. cinerea.

  • Supplemental Figure 2. Plasticity in Arabidopsis jasmonate- and salicylic acid-responsive gene expression regulated by diversity in B. cinerea.

  • Supplemental Figure 3. Plasticity in Arabidopsis basal defense-related gene expression regulated by diversity in B. cinerea.

  • Supplemental Figure 4. Overlap of network associated host genes across three Arabidopsis genotypes.

  • Supplemental Figure 5. PCA on induced transcriptome in Arabidopsis wild-type Col-0.

  • Supplemental Table 1. Correlation analysis of lesion size and camalexin against the first two principal components of the transcriptome.

  • Supplemental Table 2. Summary of negative binomial general linear models for Arabidopsis transcripts.

  • Supplemental Table 3. ANOVA model tables of lesion area and camalexin induced by B. cinerea across Arabidopsis genotypes.

  • The following materials have been deposited in the DRYAD repository under accession number http://dx.doi.org/10.5061/dryad.7gd5q.

  • Supplemental Data Set 1. Description of B. cinerea isolates.

  • Supplemental Data Set 2. Model adjusted means for defense phenotypes and score values of PCA.

  • Supplemental Data Set 3. Variation in Arabidopsis susceptibility across host genotypes driven by natural genetic variation in B. cinerea.

  • Supplemental Data Set 4. Variation in Arabidopsis camalexin accumulation across host genotypes driven by natural genetic variation in B. cinerea.

  • Supplemental Data Set 5. Model adjusted means for all transcripts measured in all genotypes of this study.

  • Supplemental Data Set 6. se for all transcripts measured in all genotypes of this study.

  • Supplemental Data Set 7. Transcriptomic analysis of significance and heritability.

  • Supplemental Data Set 8. Spearman Rank correlation.

  • Supplemental Data Set 9. Gene coexpression architecture for responses in the Arabidopsis wild-type Col-0.

  • Supplemental Data Set 10. Gene coexpression architecture for responses in the Arabidopsis mutant coi1-1.

  • Supplemental Data Set 11. Gene coexpression architecture for responses in the Arabidopsis mutant npr1-1.

  • Supplemental Data Set 12. Gene membership of the four Arabidopsis coexpression networks.

  • Supplemental Data Set 13. GO enrichment of the four Arabidopsis coexpression networks.

  • Supplemental Data Set 14. Top 5% genes regulated by B. cinerea isolates in each Arabidopsis genotype, clustered by PCA.

Acknowledgments

This work was supported by NSF Division of Integrative Organismal Systems Award 1339125 to D.J.K., by the USDA National Institute of Food and Agriculture, Hatch project number CA-D-PLS-7033-H to D.J.K., by Danish National Research Foundation Grant DNRF99 to D.J.K., and by China Scholarship Council Grant 20130624 to W.Z.

AUTHOR CONTRIBUTIONS

J.A.C. and D.J.K. conceived and designed the experiments. W.Z., R.E., D.C., and J.F. performed the experiments. W.Z., J.A.C., F.C., S.A., and D.J.K. analyzed the data. W.Z., J.A.C., and D.J.K. wrote the article.

Footnotes

  • www.plantcell.org/cgi/doi/10.1105/tpc.17.00348

  • The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantcell.org) is: Daniel J. Kliebenstein (kliebenstein{at}ucdavis.edu).

  • ↵[OPEN] Articles can be viewed without a subscription.

  • Received May 9, 2017.
  • Revised September 22, 2017.
  • Accepted October 13, 2017.
  • Published October 17, 2017.

References

  1. ↵
    1. Alter, O.,
    2. Brown, P.O.,
    3. Botstein, D.
    (2000). Singular value decomposition for genome-wide expression data processing and modeling. Proc. Natl. Acad. Sci. USA 97: 10101–10106.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. An, B.,
    2. Li, B.,
    3. Qin, G.,
    4. Tian, S.
    (2015). Function of small GTPase Rho3 in regulating growth, conidiation and virulence of Botrytis cinerea. Fungal Genet. Biol. 75: 46–55.
  3. ↵
    1. Atwell, S.,
    2. Corwin, J.A.,
    3. Soltis, N.E.,
    4. Subedy, A.,
    5. Denby, K.J.,
    6. Kliebenstein, D.J.
    (2015). Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity. Front. Microbiol. 6: 996.
    OpenUrl
  4. ↵
    1. Bednarek, P.,
    2. Piślewska-Bednarek, M.,
    3. Ver Loren van Themaat, E.,
    4. Maddula, R.K.,
    5. Svatoš, A.,
    6. Schulze-Lefert, P.
    (2011). Conservation and clade-specific diversification of pathogen-inducible tryptophan and indole glucosinolate metabolism in Arabidopsis thaliana relatives. New Phytol. 192: 713–726.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Benjamini, Y.,
    2. Hochberg, Y.
    (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. A Stat. Soc. 57: 289–300.
    OpenUrl
  6. ↵
    1. Benjamini, Y.,
    2. Drai, D.,
    3. Elmer, G.,
    4. Kafkafi, N.,
    5. Golani, I.
    (2001). Controlling the false discovery rate in behavior genetics research. Behav. Brain Res. 125: 279–284.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Benton, M.J.
    (2009). The Red Queen and the Court Jester: species diversity and the role of biotic and abiotic factors through time. Science 323: 728–732.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Bergelson, J.,
    2. Kreitman, M.,
    3. Stahl, E.A.,
    4. Tian, D.
    (2001). Evolutionary dynamics of plant R-genes. Science 292: 2281–2285.
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Birkenbihl, R.P.,
    2. Diezel, C.,
    3. Somssich, I.E.
    (2012). Arabidopsis WRKY33 is a key transcriptional regulator of hormonal and metabolic responses toward Botrytis cinerea infection. Plant Physiol. 159: 266–285.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Bolton, M.D.,
    2. Thomma, B.P.,
    3. Nelson, B.D.
    (2006). Sclerotinia sclerotiorum (Lib.) de Bary: biology and molecular traits of a cosmopolitan pathogen. Mol. Plant Pathol. 7: 1–16.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Boydom, A.
    (2013). Evaluation of detached leaf assay for assessing leaf rust (Puccinia triticina Eriks.) resistance in wheat. J. Plant Pathol. Microbiol. 4: 176–181.
  12. ↵
    1. Brodersen, P.,
    2. Malinovsky, F.G.,
    3. Hématy, K.,
    4. Newman, M.A.,
    5. Mundy, J.
    (2005). The role of salicylic acid in the induction of cell death in Arabidopsis acd11. Plant Physiol. 138: 1037–1045.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Browse, J.
    (2009). Jasmonate passes muster: a receptor and targets for the defense hormone. Annu. Rev. Plant Biol. 60: 183–205.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Caarls, L.,
    2. Pieterse, C.M.,
    3. Van Wees, S.C.
    (2015). How salicylic acid takes transcriptional control over jasmonic acid signaling. Front. Plant Sci. 6: 170.
    OpenUrlPubMed
  15. ↵
    1. Chini, A.,
    2. Boter, M.,
    3. Solano, R.
    (2009). Plant oxylipins: COI1/JAZs/MYC2 as the core jasmonic acid-signalling module. FEBS J. 276: 4682–4692.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Choquer, M.,
    2. Fournier, E.,
    3. Kunz, C.,
    4. Levis, C.,
    5. Pradier, J.M.,
    6. Simon, A.,
    7. Viaud, M.
    (2007). Botrytis cinerea virulence factors: new insights into a necrotrophic and polyphageous pathogen. FEMS Microbiol. Lett. 277: 1–10.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Clarke, J.D.,
    2. Liu, Y.,
    3. Klessig, D.F.,
    4. Dong, X.
    (1998). Uncoupling PR gene expression from NPR1 and bacterial resistance: characterization of the dominant Arabidopsis cpr6-1 mutant. Plant Cell 10: 557–569.
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Colmenares, A.J.,
    2. Aleu, J.,
    3. Durán-Patrón, R.,
    4. Collado, I.G.,
    5. Hernández-Galán, R.
    (2002). The putative role of botrydial and related metabolites in the infection mechanism of Botrytis cinerea. J. Chem. Ecol. 28: 997–1005.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Corwin, J.A.,
    2. Subedy, A.,
    3. Eshbaugh, R.,
    4. Kliebenstein, D.J.
    (2016a). Expansive phenotypic landscape of Botrytis cinerea shows differential contribution of genetic diversity and plasticity. Mol. Plant Microbe Interact. 29: 287–298.
    OpenUrl
  20. ↵
    1. Corwin, J.A.,
    2. Copeland, D.,
    3. Feusier, J.,
    4. Subedy, A.,
    5. Eshbaugh, R.,
    6. Palmer, C.,
    7. Maloof, J.,
    8. Kliebenstein, D.J.
    (2016b). The quantitative basis of the Arabidopsis innate immune system to endemic pathogens depends on pathogen genetics. PLoS Genet. 12: e1005789.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Cotoras, M.,
    2. García, C.,
    3. Mendoza, L.
    (2009). Botrytis cinerea isolates collected from grapes present different requirements for conidia germination. Mycologia 101: 287–295.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Dalmais, B.,
    2. Schumacher, J.,
    3. Moraga, J.,
    4. LE Pêcheur, P.,
    5. Tudzynski, B.,
    6. Collado, I.G.,
    7. Viaud, M.
    (2011). The Botrytis cinerea phytotoxin botcinic acid requires two polyketide synthases for production and has a redundant role in virulence with botrydial. Mol. Plant Pathol. 12: 564–579.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Dangl, J.L.,
    2. Jones, J.D.G.
    (2001). Plant pathogens and integrated defence responses to infection. Nature 411: 826–833.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Davis, J.,
    2. Yu, D.,
    3. Evans, W.,
    4. Gokirmak, T.,
    5. Chetelat, R.T.,
    6. Stotz, H.U.
    (2009). Mapping of loci from Solanum lycopersicoides conferring resistance or susceptibility to Botrytis cinerea in tomato. Theor. Appl. Genet. 119: 305–314.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Denby, K.J.,
    2. Kumar, P.,
    3. Kliebenstein, D.J.
    (2004). Identification of Botrytis cinerea susceptibility loci in Arabidopsis thaliana. Plant J. 38: 473–486.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Dong, J.,
    2. Horvath, S.
    (2007). Understanding network concepts in modules. BMC Syst. Biol. 1: 24.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Ellis, C.,
    2. Turner, J.G.
    (2001). The Arabidopsis mutant cev1 has constitutively active jasmonate and ethylene signal pathways and enhanced resistance to pathogens. Plant Cell 13: 1025–1033.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Failmezger, H.,
    2. Yuan, Y.,
    3. Rueda, O.,
    4. Markowetz, F.
    (2010). CRImage: CRImage a packgae to classify cells and calculate tumour cellularity. R Package version 1.4.0, https://rdrr.io/bioc/CRImage/f/inst/doc/CRImage.pdf.
  29. ↵
    1. Fei, Q.,
    2. Zhang, Y.,
    3. Xia, R.,
    4. Meyers, B.C.
    (2016). Small RNAs add zing to the zig-zag-zig model of plant defenses. Mol. Plant Microbe Interact. 29: 165–169.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Ferrari, S.,
    2. Plotnikova, J.M.,
    3. De Lorenzo, G.,
    4. Ausubel, F.M.
    (2003). Arabidopsis local resistance to Botrytis cinerea involves salicylic acid and camalexin and requires EDS4 and PAD2, but not SID2, EDS5 or PAD4. Plant J. 35: 193–205.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Ferrari, S.,
    2. Galletti, R.,
    3. Denoux, C.,
    4. De Lorenzo, G.,
    5. Ausubel, F.M.,
    6. Dewdney, J.
    (2007). Resistance to Botrytis cinerea induced in Arabidopsis by elicitors is independent of salicylic acid, ethylene, or jasmonate signaling but requires PHYTOALEXIN DEFICIENT3. Plant Physiol. 144: 367–379.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Finkers, R.,
    2. van Heusden, A.W.,
    3. Meijer-Dekens, F.,
    4. van Kan, J.A.,
    5. Maris, P.,
    6. Lindhout, P.
    (2007a). The construction of a Solanum habrochaites LYC4 introgression line population and the identification of QTLs for resistance to Botrytis cinerea. Theor. Appl. Genet. 114: 1071–1080.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Finkers, R.,
    2. van den Berg, P.,
    3. van Berloo, R.,
    4. ten Have, A.,
    5. van Heusden, A.W.,
    6. van Kan, J.A.,
    7. Lindhout, P.
    (2007b). Three QTLs for Botrytis cinerea resistance in tomato. Theor. Appl. Genet. 114: 585–593.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Fisher, M.C.,
    2. Henk, D.A.,
    3. Briggs, C.J.,
    4. Brownstein, J.S.,
    5. Madoff, L.C.,
    6. McCraw, S.L.,
    7. Gurr, S.J.
    (2012). Emerging fungal threats to animal, plant and ecosystem health. Nature 484: 186–194.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Fox, J.,
    2. Weisberg, S.
    (2011). An R Companion to Applied Regression. (Thousand Oaks, CA: SAGE).
  36. ↵
    1. Frerigmann, H.,
    2. Glawischnig, E.,
    3. Gigolashvili, T.
    (2015). The role of MYB34, MYB51 and MYB122 in the regulation of camalexin biosynthesis in Arabidopsis thaliana. Front. Plant Sci. 6: 654.
    OpenUrl
  37. ↵
    1. Frerigmann, H.,
    2. Piślewska-Bednarek, M.,
    3. Sánchez-Vallet, A.,
    4. Molina, A.,
    5. Glawischnig, E.,
    6. Gigolashvili, T.,
    7. Bednarek, P.
    (2016). Regulation of pathogen-triggered tryptophan metabolism in Arabidopsis thaliana by MYB transcription factors and indole glucosinolate conversion products. Mol. Plant 9: 682–695.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Glazebrook, J.
    (2001). Genes controlling expression of defense responses in Arabidopsis--2001 status. Curr. Opin. Plant Biol. 4: 301–308.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Glazebrook, J.
    (2005). Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu. Rev. Phytopathol. 43: 205–227.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Glazebrook, J.,
    2. Ausubel, F.M.
    (1994). Isolation of phytoalexin-deficient mutants of Arabidopsis thaliana and characterization of their interactions with bacterial pathogens. Proc. Natl. Acad. Sci. USA 91: 8955–8959.
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Govrin, E.M.,
    2. Levine, A.
    (2000). The hypersensitive response facilitates plant infection by the necrotrophic pathogen Botrytis cinerea. Curr. Biol. 10: 751–757.
    OpenUrlCrossRefPubMed
  42. ↵
    1. Han, L.,
    2. Li, G.-J.,
    3. Yang, K.-Y.,
    4. Mao, G.,
    5. Wang, R.,
    6. Liu, Y.,
    7. Zhang, S.
    (2010). Mitogen-activated protein kinase 3 and 6 regulate Botrytis cinerea-induced ethylene production in Arabidopsis. Plant J. 64: 114–127.
    OpenUrlPubMed
  43. ↵
    1. He, Y.,
    2. Chung, E.H.,
    3. Hubert, D.A.,
    4. Tornero, P.,
    5. Dangl, J.L.
    (2012). Specific missense alleles of the arabidopsis jasmonic acid co-receptor COI1 regulate innate immune receptor accumulation and function. PLoS Genet. 8: e1003018.
    OpenUrlPubMed
  44. ↵
    1. Hevia, M.A.,
    2. Canessa, P.,
    3. Müller-Esparza, H.,
    4. Larrondo, L.F.
    (2015). A circadian oscillator in the fungus Botrytis cinerea regulates virulence when infecting Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 112: 8744–8749.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Hintze, J.L.,
    2. Nelson, R.D.
    (1998). Violin plots: a box plot-density trace synergism. Am. Stat. 52: 181–184.
    OpenUrlCrossRef
  46. ↵
    1. Hörger, A.C.,
    2. Ilyas, M.,
    3. Stephan, W.,
    4. Tellier, A.,
    5. van der Hoorn, R.A.,
    6. Rose, L.E.
    (2012). Balancing selection at the tomato RCR3 Guardee gene family maintains variation in strength of pathogen defense. PLoS Genet. 8: e1002813.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Howard, B.E.,
    2. Hu, Q.,
    3. Babaoglu, A.C.,
    4. Chandra, M.,
    5. Borghi, M.,
    6. Tan, X.,
    7. He, L.,
    8. Winter-Sederoff, H.,
    9. Gassmann, W.,
    10. Veronese, P.,
    11. Heber, S.
    (2013). High-throughput RNA sequencing of pseudomonas-infected Arabidopsis reveals hidden transcriptome complexity and novel splice variants. PLoS One 8: e74183.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Iakovidis, M.,
    2. Teixeira, P.J.,
    3. Exposito-Alonso, M.,
    4. Cowper, M.G.,
    5. Law, T.F.,
    6. Liu, Q.,
    7. Vu, M.C.,
    8. Dang, T.M.,
    9. Corwin, J.A.,
    10. Weigel, D.,
    11. Dangl, J.L.,
    12. Grant, S.R.
    (2016). Effector-triggered immune response in Arabidopsis thaliana is a quantitative trait. Genetics 204: 337–353.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Johnson, K.C.M.,
    2. Zhao, J.,
    3. Wu, Z.,
    4. Roth, C.,
    5. Lipka, V.,
    6. Wiermer, M.,
    7. Li, X.
    (2017). The putative kinase substrate MUSE7 negatively impacts the accumulation of NLR proteins. Plant J. 89: 1174–1183.
    OpenUrl
  50. ↵
    1. Jones, J.D.,
    2. Dangl, J.L.
    (2006). The plant immune system. Nature 444: 323–329.
    OpenUrlCrossRefPubMed
  51. ↵
    1. Karasov, T.L., et al
    . (2014). The long-term maintenance of a resistance polymorphism through diffuse interactions. Nature 512: 436–440.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Keller, H.,
    2. Boyer, L.,
    3. Abad, P.
    (2016). Disease susceptibility in the Zig-Zag model of host-microbe interactions: only a consequence of immune suppression? Mol. Plant Pathol. 17: 475–479.
    OpenUrl
  53. ↵
    1. Kim, Y.,
    2. Tsuda, K.,
    3. Igarashi, D.,
    4. Hillmer, R.A.,
    5. Sakakibara, H.,
    6. Myers, C.L.,
    7. Katagiri, F.
    (2014). Mechanisms underlying robustness and tunability in a plant immune signaling network. Cell Host Microbe 15: 84–94.
    OpenUrlCrossRefPubMed
  54. ↵
    1. Kliebenstein, D.J.
    (2009). Quantification of variation in expression networks. Methods Mol. Biol. 553: 227–245.
    OpenUrlPubMed
  55. ↵
    1. Kliebenstein, D.J.,
    2. Rowe, H.C.,
    3. Denby, K.J.
    (2005). Secondary metabolites influence Arabidopsis/Botrytis interactions: variation in host production and pathogen sensitivity. Plant J. 44: 25–36.
    OpenUrlCrossRefPubMed
  56. ↵
    1. Kretschmer, M.,
    2. Leroch, M.,
    3. Mosbach, A.,
    4. Walker, A.S.,
    5. Fillinger, S.,
    6. Mernke, D.,
    7. Schoonbeek, H.J.,
    8. Pradier, J.M.,
    9. Leroux, P.,
    10. De Waard, M.A.,
    11. Hahn, M.
    (2009). Fungicide-driven evolution and molecular basis of multidrug resistance in field populations of the grey mould fungus Botrytis cinerea. PLoS Pathog. 5: e1000696.
    OpenUrlCrossRefPubMed
  57. ↵
    1. Kumar, R.,
    2. Ichihashi, Y.,
    3. Kimura, S.,
    4. Chitwood, D.H.,
    5. Headland, L.R.,
    6. Peng, J.,
    7. Maloof, J.N.,
    8. Sinha, N.R.
    (2012). A high-throughput method for Illumina RNA-seq library preparation. Front. Plant Sci. 3: 202.
    OpenUrlPubMed
  58. ↵
    1. Kumari, S.,
    2. Tayal, P.,
    3. Sharma, E.,
    4. Kapoor, R.
    (2014). Analyses of genetic and pathogenic variability among Botrytis cinerea isolates. Microbiol. Res. 169: 862–872.
    OpenUrl
  59. ↵
    1. Lachowiec, J.,
    2. Queitsch, C.,
    3. Kliebenstein, D.J.
    (2016). Molecular mechanisms governing differential robustness of development and environmental responses in plants. Ann. Bot. 117: 795–809.
    OpenUrlCrossRefPubMed
  60. ↵
    1. Lamesch, P., et al
    . (2012). The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 40: D1202–D1210.
    OpenUrlCrossRefPubMed
  61. ↵
    1. Langmead, B.,
    2. Trapnell, C.,
    3. Pop, M.,
    4. Salzberg, S.L.
    (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10: R25.
    OpenUrlCrossRefPubMed
  62. ↵
    1. Lenth, R.V.
    (2016). Least-squares means: The R package lsmeans. J. Stat. Softw. 69: 1–33.
    OpenUrlCrossRef
  63. ↵
    1. Lewis, L.A., et al
    . (2015). Transcriptional dynamics driving MAMP-triggered immunity and pathogen effector-mediated immunosuppression in Arabidopsis leaves following infection with Pseudomonas syringae pv tomato DC3000. Plant Cell 27: 3038–3064.
    OpenUrlAbstract/FREE Full Text
  64. ↵
    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.
    OpenUrlCrossRefPubMed
  65. ↵
    1. Liefting, M.,
    2. Hoffmann, A.A.,
    3. Ellers, J.
    (2009). Plasticity versus environmental canalization: population differences in thermal responses along a latitudinal gradient in Drosophila serrata. Evolution 63: 1954–1963.
    OpenUrlCrossRefPubMed
  66. ↵
    1. Liu, F.,
    2. Jiang, H.,
    3. Ye, S.,
    4. Chen, W.P.,
    5. Liang, W.,
    6. Xu, Y.,
    7. Sun, B.,
    8. Sun, J.,
    9. Wang, Q.,
    10. Cohen, J.D.,
    11. Li, C.
    (2010). The Arabidopsis P450 protein CYP82C2 modulates jasmonate-induced root growth inhibition, defense gene expression and indole glucosinolate biosynthesis. Cell Res. 20: 539–552.
    OpenUrlCrossRefPubMed
  67. ↵
    1. Liu, T.,
    2. Liu, Z.,
    3. Song, C.,
    4. Hu, Y.,
    5. Han, Z.,
    6. She, J.,
    7. Fan, F.,
    8. Wang, J.,
    9. Jin, C.,
    10. Chang, J.,
    11. Zhou, J.M.,
    12. Chai, J.
    (2012). Chitin-induced dimerization activates a plant immune receptor. Science 336: 1160–1164.
    OpenUrlAbstract/FREE Full Text
  68. ↵
    1. Lopez-Cruz, J.,
    2. Crespo-Salvador, O.,
    3. Fernandez-Crespo, E.,
    4. Garcia-Agustin, P.,
    5. Gonzalez-Bosch, C.
    (2017). Absence of Cu-Zn-superoxide dismutase BCSOD1 reduces Botrytis cinerea virulence in Arabidopsis and in tomato plants, which reveals interplay among ROS, callose and signaling pathways. Mol. Plant Pathol. 18: 16–31.
    OpenUrl
  69. ↵
    1. Lyu, X.,
    2. Shen, C.,
    3. Fu, Y.,
    4. Xie, J.,
    5. Jiang, D.,
    6. Li, G.,
    7. Cheng, J.
    (2016). A small secreted virulence-related protein is essential for the necrotrophic interactions of Sclerotinia sclerotiorum with its host plants. PLoS Pathog. 12: e1005435.
    OpenUrl
  70. ↵
    1. Maere, S.,
    2. Heymans, K.,
    3. Kuiper, M.
    (2005). BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21: 3448–3449.
    OpenUrlCrossRefPubMed
  71. ↵
    1. Mahesh, S.K.,
    2. Liu, H.,
    3. Qiu, D.-w.
    (2012). The role of radical burst in plant defense responses to necrotrophic fungi. J. Integr. Agric. 11: 1305–1312.
    OpenUrl
  72. ↵
    1. Mantel, N.
    (1967). The detection of disease clustering and a generalized regression approach. Cancer Res. 27: 209–220.
    OpenUrlAbstract/FREE Full Text
  73. ↵
    1. Mao, G.,
    2. Meng, X.,
    3. Liu, Y.,
    4. Zheng, Z.,
    5. Chen, Z.,
    6. Zhang, S.
    (2011). Phosphorylation of a WRKY transcription factor by two pathogen-responsive MAPKs drives phytoalexin biosynthesis in Arabidopsis. Plant Cell 23: 1639–1653.
    OpenUrlAbstract/FREE Full Text
  74. ↵
    1. Mauricio, R.,
    2. Stahl, E.A.,
    3. Korves, T.,
    4. Tian, D.,
    5. Kreitman, M.,
    6. Bergelson, J.
    (2003). Natural selection for polymorphism in the disease resistance gene Rps2 of Arabidopsis thaliana. Genetics 163: 735–746.
    OpenUrlAbstract/FREE Full Text
  75. ↵
    1. Meena, M.,
    2. Prasad, V.,
    3. Zehra, A.,
    4. Gupta, V.K.,
    5. Upadhyay, R.S.
    (2015). Mannitol metabolism during pathogenic fungal-host interactions under stressed conditions. Front. Microbiol. 6: 1019.
    OpenUrl
  76. ↵
    1. Mendiburu, F.D.
    (2016). Agricolae: Statistical Procedures for Agricultural Research. https://cran.r-project.org/web/packages/agricolae/agricolae.pdf.
  77. ↵
    1. Mengiste, T.
    (2012). Plant immunity to necrotrophs. Annu. Rev. Phytopathol. 50: 267–294.
    OpenUrlCrossRefPubMed
  78. ↵
    1. Mengiste, T.,
    2. Chen, X.,
    3. Salmeron, J.,
    4. Dietrich, R.
    (2003). The BOTRYTIS SUSCEPTIBLE1 gene encodes an R2R3MYB transcription factor protein that is required for biotic and abiotic stress responses in Arabidopsis. Plant Cell 15: 2551–2565.
    OpenUrlAbstract/FREE Full Text
  79. ↵
    1. Michielse, C.B.,
    2. Becker, M.,
    3. Heller, J.,
    4. Moraga, J.,
    5. Collado, I.G.,
    6. Tudzynski, P.
    (2011). The Botrytis cinerea Reg1 protein, a putative transcriptional regulator, is required for pathogenicity, conidiogenesis, and the production of secondary metabolites. Mol. Plant Microbe Interact. 24: 1074–1085.
    OpenUrlCrossRefPubMed
  80. ↵
    1. Millet, Y.A.,
    2. Danna, C.H.,
    3. Clay, N.K.,
    4. Songnuan, W.,
    5. Simon, M.D.,
    6. Werck-Reichhart, D.,
    7. Ausubel, F.M.
    (2010). Innate immune responses activated in Arabidopsis roots by microbe-associated molecular patterns. Plant Cell 22: 973–990.
    OpenUrlAbstract/FREE Full Text
  81. ↵
    1. Moreau, M.,
    2. Tian, M.,
    3. Klessig, D.F.
    (2012). Salicylic acid binds NPR3 and NPR4 to regulate NPR1-dependent defense responses. Cell Res. 22: 1631–1633.
    OpenUrlCrossRefPubMed
  82. ↵
    1. Mou, Z.,
    2. Fan, W.,
    3. Dong, X.
    (2003). Inducers of plant systemic acquired resistance regulate NPR1 function through redox changes. Cell 113: 935–944.
    OpenUrlCrossRefPubMed
  83. ↵
    1. Mukhtar, M.S.,
    2. Nishimura, M.T.,
    3. Dangl, J.
    (2009). NPR1 in plant defense: It’s not over 'til it’s turned over. Cell 137: 804–806.
    OpenUrlCrossRefPubMed
  84. ↵
    1. Mulema, J.M.,
    2. Denby, K.J.
    (2012). Spatial and temporal transcriptomic analysis of the Arabidopsis thaliana-Botrytis cinerea interaction. Mol. Biol. Rep. 39: 4039–4049.
    OpenUrlCrossRefPubMed
  85. ↵
    1. Nawrath, C.,
    2. Heck, S.,
    3. Parinthawong, N.,
    4. Métraux, J.P.
    (2002). EDS5, an essential component of salicylic acid-dependent signaling for disease resistance in Arabidopsis, is a member of the MATE transporter family. Plant Cell 14: 275–286.
    OpenUrlAbstract/FREE Full Text
  86. ↵
    1. Nie, P.,
    2. Li, X.,
    3. Wang, S.,
    4. Guo, J.,
    5. Zhao, H.,
    6. Niu, D.
    (2017). Induced systemic resistance against Botrytis cinerea by Bacillus cereus AR156 through a JA/ET- and NPR1-dependent signaling pathway and activates PAMP-triggered immunity in Arabidopsis. Front. Plant Sci. 8: 238.
    OpenUrl
  87. ↵
    1. Nikolayeva, O.,
    2. Robinson, M.D.
    (2014). edgeR for differential RNA-seq and ChIP-seq analysis: an application to stem cell biology. Methods Mol. Biol. 1150: 45–79.
    OpenUrlCrossRefPubMed
  88. ↵
    1. Noda, J.,
    2. Brito, N.,
    3. González, C.
    (2010). The Botrytis cinerea xylanase Xyn11A contributes to virulence with its necrotizing activity, not with its catalytic activity. BMC Plant Biol. 10: 38.
    OpenUrlCrossRefPubMed
  89. ↵
    1. Okmen, B.,
    2. Doehlemann, G.
    (2014). Inside plant: biotrophic strategies to modulate host immunity and metabolism. Curr. Opin. Plant Biol. 20: 19–25.
    OpenUrlCrossRefPubMed
  90. ↵
    1. Oksanen, J.,
    2. Blanchet, F.G.,
    3. Friendly, M.,
    4. Kindt, R.,
    5. Legendre, P.,
    6. McGlinn, D.,
    7. Minchin, P.R.,
    8. O'Hara, R.B.,
    9. Simpson, G.L.,
    10. Solymos, P.,
    11. Stevens, M.H.H.,
    12. Szoecs, E.,
    13. Wagner, H.
    (2016). vegan: Community Ecology Package. https://cran.r-project.org/web/packages/vegan/vegan.pdf.
  91. ↵
    1. Palacio-López, K.,
    2. Beckage, B.,
    3. Scheiner, S.,
    4. Molofsky, J.
    (2015). The ubiquity of phenotypic plasticity in plants: a synthesis. Ecol. Evol. 5: 3389–3400.
    OpenUrlCrossRefPubMed
  92. ↵
    1. Pandey, S.P.,
    2. Roccaro, M.,
    3. Schön, M.,
    4. Logemann, E.,
    5. Somssich, I.E.
    (2010). Transcriptional reprogramming regulated by WRKY18 and WRKY40 facilitates powdery mildew infection of Arabidopsis. Plant J. 64: 912–923.
    OpenUrlCrossRefPubMed
  93. ↵
    1. Passarinho, P.A.,
    2. Van Hengel, A.J.,
    3. Fransz, P.F.,
    4. de Vries, S.C.
    (2001). Expression pattern of the Arabidopsis thaliana AtEP3/AtchitIV endochitinase gene. Planta 212: 556–567.
    OpenUrlCrossRefPubMed
  94. ↵
    1. Patel, T.K.,
    2. Williamson, J.D.
    (2016). Mannitol in plants, fungi, and plant-fungal interactions. Trends Plant Sci. 21: 486–497.
    OpenUrl
  95. ↵
    1. Pau, G.,
    2. Fuchs, F.,
    3. Sklyar, O.,
    4. Boutros, M.,
    5. Huber, W.
    (2010). EBImage--an R package for image processing with applications to cellular phenotypes. Bioinformatics 26: 979–981.
    OpenUrlCrossRefPubMed
  96. ↵
    1. Pearson, M.N.,
    2. Bailey, A.M.
    (2013). Viruses of Botrytis. Adv. Virus Res. 86: 249–272.
    OpenUrlCrossRefPubMed
  97. ↵
    1. Pedras, M.S.,
    2. Jha, M.,
    3. Okeola, O.G.
    (2005). Camalexin induces detoxification of the phytoalexin brassinin in the plant pathogen Leptosphaeria maculans. Phytochemistry 66: 2609–2616.
    OpenUrlPubMed
  98. ↵
    1. Pedras, M.S.,
    2. Minic, Z.,
    3. Jha, M.
    (2008). Brassinin oxidase, a fungal detoxifying enzyme to overcome a plant defense -- purification, characterization and inhibition. FEBS J. 275: 3691–3705.
    OpenUrlPubMed
  99. ↵
    1. Pedras, M.S.,
    2. Minic, Z.,
    3. Sarma-Mamillapalle, V.K.
    (2009). Synthetic inhibitors of the fungal detoxifying enzyme brassinin oxidase based on the phytoalexin camalexin scaffold. J. Agric. Food Chem. 57: 2429–2435.
    OpenUrlPubMed
  100. ↵
    1. Pedras, M.S.,
    2. Hossain, S.,
    3. Snitynsky, R.B.
    (2011). Detoxification of cruciferous phytoalexins in Botrytis cinerea: spontaneous dimerization of a camalexin metabolite. Phytochemistry 72: 199–206.
    OpenUrl
  101. ↵
    1. Pedras, M.S.,
    2. Gadagi, R.S.,
    3. Jha, M.,
    4. Sarma-Mamillapalle, V.K.
    (2007). Detoxification of the phytoalexin brassinin by isolates of Leptosphaeria maculans pathogenic on brown mustard involves an inducible hydrolase. Phytochemistry 68: 1572–1578.
    OpenUrlCrossRefPubMed
  102. ↵
    1. Penninckx, I.A.,
    2. Eggermont, K.,
    3. Terras, F.R.,
    4. Thomma, B.P.,
    5. De Samblanx, G.W.,
    6. Buchala, A.,
    7. Métraux, J.P.,
    8. Manners, J.M.,
    9. Broekaert, W.F.
    (1996). Pathogen-induced systemic activation of a plant defensin gene in Arabidopsis follows a salicylic acid-independent pathway. Plant Cell 8: 2309–2323.
    OpenUrlAbstract/FREE Full Text
  103. ↵
    1. Pieterse, C.M.,
    2. Van Loon, L.C.
    (2004). NPR1: the spider in the web of induced resistance signaling pathways. Curr. Opin. Plant Biol. 7: 456–464.
    OpenUrlCrossRefPubMed
  104. ↵
    1. Plaza, V.,
    2. Lagües, Y.,
    3. Carvajal, M.,
    4. Pérez-García, L.A.,
    5. Mora-Montes, H.M.,
    6. Canessa, P.,
    7. Larrondo, L.F.,
    8. Castillo, L.
    (2015). bcpmr1 encodes a P-type Ca(2+)/Mn(2+)-ATPase mediating cell-wall integrity and virulence in the phytopathogen Botrytis cinerea. Fungal Genet. Biol. 76: 36–46.
    OpenUrl
  105. ↵
    1. Pré, M.,
    2. Atallah, M.,
    3. Champion, A.,
    4. De Vos, M.,
    5. Pieterse, C.M.,
    6. Memelink, J.
    (2008). The AP2/ERF domain transcription factor ORA59 integrates jasmonic acid and ethylene signals in plant defense. Plant Physiol. 147: 1347–1357.
    OpenUrlAbstract/FREE Full Text
  106. ↵
    1. Pritchard, L.,
    2. Birch, P.R.
    (2014). The zigzag model of plant-microbe interactions: is it time to move on? Mol. Plant Pathol. 15: 865–870.
    OpenUrlCrossRefPubMed
  107. ↵
    1. Queitsch, C.,
    2. Sangster, T.A.,
    3. Lindquist, S.
    (2002). Hsp90 as a capacitor of phenotypic variation. Nature 417: 618–624.
    OpenUrlCrossRefPubMed
  108. ↵
    1. R Core Team
    (2011). R: A Language and Environment for Statistical Computing. (Vienna, Austria: R Foundation for Statistical Computing).
  109. ↵
    1. Rajniak, J.,
    2. Barco, B.,
    3. Clay, N.K.,
    4. Sattely, E.S.
    (2015). A new cyanogenic metabolite in Arabidopsis required for inducible pathogen defence. Nature 525: 376–379.
    OpenUrlCrossRefPubMed
  110. ↵
    1. Ren, D.,
    2. Liu, Y.,
    3. Yang, K.Y.,
    4. Han, L.,
    5. Mao, G.,
    6. Glazebrook, J.,
    7. Zhang, S.
    (2008). A fungal-responsive MAPK cascade regulates phytoalexin biosynthesis in Arabidopsis. Proc. Natl. Acad. Sci. USA 105: 5638–5643.
    OpenUrlAbstract/FREE Full Text
  111. ↵
    1. Robinson, M.D.,
    2. Smyth, G.K.
    (2008). Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics 9: 321–332.
    OpenUrlCrossRefPubMed
  112. ↵
    1. Robinson, M.D.,
    2. Oshlack, A.
    (2010). A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 11: R25.
    OpenUrlCrossRefPubMed
  113. ↵
    1. Robinson, M.D.,
    2. McCarthy, D.J.,
    3. Smyth, G.K.
    (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140.
    OpenUrlCrossRefPubMed
  114. ↵
    1. Rose, L.E.,
    2. Michelmore, R.W.,
    3. Langley, C.H.
    (2007). Natural variation in the Pto disease resistance gene within species of wild tomato (Lycopersicon). II. Population genetics of Pto. Genetics 175: 1307–1319.
    OpenUrlAbstract/FREE Full Text
  115. ↵
    1. Rowe, H.C.,
    2. Kliebenstein, D.J.
    (2007). Elevated genetic variation within virulence-associated Botrytis cinerea polygalacturonase loci. Mol. Plant Microbe Interact. 20: 1126–1137.
    OpenUrlCrossRefPubMed
  116. ↵
    1. Rowe, H.C.,
    2. Kliebenstein, D.J.
    (2008). Complex genetics control natural variation in Arabidopsis thaliana resistance to Botrytis cinerea. Genetics 180: 2237–2250.
    OpenUrlAbstract/FREE Full Text
  117. ↵
    1. Rowe, H.C.,
    2. Walley, J.W.,
    3. Corwin, J.,
    4. Chan, E.K.,
    5. Dehesh, K.,
    6. Kliebenstein, D.J.
    (2010). Deficiencies in jasmonate-mediated plant defense reveal quantitative variation in Botrytis cinerea pathogenesis. PLoS Pathog. 6: e1000861.
    OpenUrlCrossRefPubMed
  118. ↵
    1. Saga, H.,
    2. Ogawa, T.,
    3. Kai, K.,
    4. Suzuki, H.,
    5. Ogata, Y.,
    6. Sakurai, N.,
    7. Shibata, D.,
    8. Ohta, D.
    (2012). Identification and characterization of ANAC042, a transcription factor family gene involved in the regulation of camalexin biosynthesis in Arabidopsis. Mol. Plant Microbe Interact. 25: 684–696.
    OpenUrlCrossRefPubMed
  119. ↵
    1. Saleh, A.,
    2. Withers, J.,
    3. Mohan, R.,
    4. Marqués, J.,
    5. Gu, Y.,
    6. Yan, S.,
    7. Zavaliev, R.,
    8. Nomoto, M.,
    9. Tada, Y.,
    10. Dong, X.
    (2015). Posttranslational modifications of the master transcriptional regulator NPR1 enable dynamic but tight control of plant immune responses. Cell Host Microbe 18: 169–182.
    OpenUrl
  120. ↵
    1. Sato, M.,
    2. Tsuda, K.,
    3. Wang, L.,
    4. Coller, J.,
    5. Watanabe, Y.,
    6. Glazebrook, J.,
    7. Katagiri, F.
    (2010). Network modeling reveals prevalent negative regulatory relationships between signaling sectors in Arabidopsis immune signaling. PLoS Pathog. 6: e1001011.
    OpenUrlCrossRefPubMed
  121. ↵
    1. Schuhegger, R.,
    2. Nafisi, M.,
    3. Mansourova, M.,
    4. Petersen, B.L.,
    5. Olsen, C.E.,
    6. Svatos, A.,
    7. Halkier, B.A.,
    8. Glawischnig, E.
    (2006). CYP71B15 (PAD3) catalyzes the final step in camalexin biosynthesis. Plant Physiol. 141: 1248–1254.
    OpenUrlAbstract/FREE Full Text
  122. ↵
    1. Schumacher, J.,
    2. Pradier, J.M.,
    3. Simon, A.,
    4. Traeger, S.,
    5. Moraga, J.,
    6. Collado, I.G.,
    7. Viaud, M.,
    8. Tudzynski, B.
    (2012). Natural variation in the VELVET gene bcvel1 affects virulence and light-dependent differentiation in Botrytis cinerea. PLoS One 7: e47840.
    OpenUrlCrossRefPubMed
  123. ↵
    1. Schumacher, J.,
    2. Simon, A.,
    3. Cohrs, K.C.,
    4. Traeger, S.,
    5. Porquier, A.,
    6. Dalmais, B.,
    7. Viaud, M.,
    8. Tudzynski, B.
    (2015). The VELVET complex in the gray mold fungus Botrytis cinerea: Impact of BcLAE1 on differentiation, secondary metabolism, and virulence. Mol. Plant Microbe Interact. 28: 659–674.
    OpenUrl
  124. ↵
    1. Serrano, I.,
    2. Audran, C.,
    3. Rivas, S.
    (2016). Chloroplasts at work during plant innate immunity. J. Exp. Bot. 67: 3845–3854.
    OpenUrlCrossRefPubMed
  125. ↵
    1. Shannon, P.,
    2. Markiel, A.,
    3. Ozier, O.,
    4. Baliga, N.S.,
    5. Wang, J.T.,
    6. Ramage, D.,
    7. Amin, N.,
    8. Schwikowski, B.,
    9. Ideker, T.
    (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13: 2498–2504.
    OpenUrlAbstract/FREE Full Text
  126. ↵
    1. Sharma, H.C.,
    2. Pampapathy, G.,
    3. Dhillon, M.K.,
    4. Ridsdill-Smith, J.T.
    (2005). Detached leaf assay to screen for host plant resistance to Helicoverpa armigera. J. Econ. Entomol. 98: 568–576.
    OpenUrlCrossRefPubMed
  127. ↵
    1. Sheard, L.B., et al
    . (2010). Jasmonate perception by inositol-phosphate-potentiated COI1-JAZ co-receptor. Nature 468: 400–405.
    OpenUrlCrossRefPubMed
  128. ↵
    1. Shlezinger, N.,
    2. Minz, A.,
    3. Gur, Y.,
    4. Hatam, I.,
    5. Dagdas, Y.F.,
    6. Talbot, N.J.,
    7. Sharon, A.
    (2011). Anti-apoptotic machinery protects the necrotrophic fungus Botrytis cinerea from host-induced apoptotic-like cell death during plant infection. PLoS Pathog. 7: e1002185.
    OpenUrlCrossRefPubMed
  129. ↵
    1. Spanu, P.D., et al
    . (2010). Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism. Science 330: 1543–1546.
    OpenUrlAbstract/FREE Full Text
  130. ↵
    1. Staal, J.,
    2. Kaliff, M.,
    3. Dewaele, E.,
    4. Persson, M.,
    5. Dixelius, C.
    (2008). RLM3, a TIR domain encoding gene involved in broad-range immunity of Arabidopsis to necrotrophic fungal pathogens. Plant J. 55: 188–200.
    OpenUrlCrossRefPubMed
  131. ↵
    1. Staswick, P.E.,
    2. Tiryaki, I.,
    3. Rowe, M.L.
    (2002). Jasmonate response locus JAR1 and several related Arabidopsis genes encode enzymes of the firefly luciferase superfamily that show activity on jasmonic, salicylic, and indole-3-acetic acids in an assay for adenylation. Plant Cell 14: 1405–1415.
    OpenUrlAbstract/FREE Full Text
  132. ↵
    1. Stefanato, F.L.,
    2. Abou-Mansour, E.,
    3. Buchala, A.,
    4. Kretschmer, M.,
    5. Mosbach, A.,
    6. Hahn, M.,
    7. Bochet, C.G.,
    8. Metraux, J.P.,
    9. Schoonbeek, H.J.
    (2009). The ABC transporter BcatrB from Botrytis cinerea exports camalexin and is a virulence factor on Arabidopsis thaliana. Plant J. 58: 499–510.
    OpenUrlCrossRefPubMed
  133. ↵
    1. Stergiopoulos, I.,
    2. de Wit, P.J.
    (2009). Fungal effector proteins. Annu. Rev. Phytopathol. 47: 233–263.
    OpenUrlCrossRefPubMed
  134. ↵
    1. Strimmer, K.
    (2008). fdrtool: a versatile R package for estimating local and tail area-based false discovery rates. Bioinformatics 24: 1461–1462.
    OpenUrlCrossRefPubMed
  135. ↵
    1. ten Have, A.,
    2. Mulder, W.,
    3. Visser, J.,
    4. van Kan, J.A.
    (1998). The endopolygalacturonase gene Bcpg1 is required for full virulence of Botrytis cinerea. Mol. Plant Microbe Interact. 11: 1009–1016.
    OpenUrlCrossRefPubMed
  136. ↵
    1. Thaler, J.S.,
    2. Humphrey, P.T.,
    3. Whiteman, N.K.
    (2012). Evolution of jasmonate and salicylate signal crosstalk. Trends Plant Sci. 17: 260–270.
    OpenUrlCrossRefPubMed
  137. ↵
    1. Thines, B.,
    2. Katsir, L.,
    3. Melotto, M.,
    4. Niu, Y.,
    5. Mandaokar, A.,
    6. Liu, G.,
    7. Nomura, K.,
    8. He, S.Y.,
    9. Howe, G.A.,
    10. Browse, J.
    (2007). JAZ repressor proteins are targets of the SCF(COI1) complex during jasmonate signalling. Nature 448: 661–665.
    OpenUrlCrossRefPubMed
  138. ↵
    1. Thomma, B.P.,
    2. Eggermont, K.,
    3. Penninckx, I.A.,
    4. Mauch-Mani, B.,
    5. Vogelsang, R.,
    6. Cammue, B.P.A.,
    7. Broekaert, W.F.
    (1998). Separate jasmonate-dependent and salicylate-dependent defense-response pathways in Arabidopsis are essential for resistance to distinct microbial pathogens. Proc. Natl. Acad. Sci. USA 95: 15107–15111.
    OpenUrlAbstract/FREE Full Text
  139. ↵
    1. Thomma, B.P.,
    2. Nelissen, I.,
    3. Eggermont, K.,
    4. Broekaert, W.F.
    (1999a). Deficiency in phytoalexin production causes enhanced susceptibility of Arabidopsis thaliana to the fungus Alternaria brassicicola. Plant J. 19: 163–171.
    OpenUrlCrossRefPubMed
  140. ↵
    1. Thomma, B.P.,
    2. Eggermont, K.,
    3. Tierens, K.F.,
    4. Broekaert, W.F.
    (1999b). Requirement of functional ethylene-insensitive 2 gene for efficient resistance of Arabidopsis to infection by Botrytis cinerea. Plant Physiol. 121: 1093–1102.
    OpenUrlAbstract/FREE Full Text
  141. ↵
    1. Tsuda, K.,
    2. Mine, A.,
    3. Bethke, G.,
    4. Igarashi, D.,
    5. Botanga, C.J.,
    6. Tsuda, Y.,
    7. Glazebrook, J.,
    8. Sato, M.,
    9. Katagiri, F.
    (2013). Dual regulation of gene expression mediated by extended MAPK activation and salicylic acid contributes to robust innate immunity in Arabidopsis thaliana. PLoS Genet. 9: e1004015.
    OpenUrlCrossRefPubMed
  142. ↵
    1. Van der Hoorn, R.A.,
    2. De Wit, P.J.,
    3. Joosten, M.H.
    (2002). Balancing selection favors guarding resistance proteins. Trends Plant Sci. 7: 67–71.
    OpenUrlCrossRefPubMed
  143. ↵
    1. van Kan, J.A.
    (2006). Licensed to kill: the lifestyle of a necrotrophic plant pathogen. Trends Plant Sci. 11: 247–253.
    OpenUrlCrossRefPubMed
  144. ↵
    1. Van Valen, L.
    (1973). A new evolutionary law. Evol. Theory 1: 1–30.
    OpenUrl
  145. ↵
    1. Venables, W.N.,
    2. Ripley, B.D.
    (2002). Modern Applied Statistics with S. (New York: Springer).
  146. ↵
    1. Vlot, A.C.,
    2. Dempsey, D.A.,
    3. Klessig, D.F.
    (2009). Salicylic acid, a multifaceted hormone to combat disease. Annu. Rev. Phytopathol. 47: 177–206.
    OpenUrlCrossRefPubMed
  147. ↵
    1. Voegele, R.T.,
    2. Hahn, M.,
    3. Lohaus, G.,
    4. Link, T.,
    5. Heiser, I.,
    6. Mendgen, K.
    (2005). Possible roles for mannitol and mannitol dehydrogenase in the biotrophic plant pathogen Uromyces fabae. Plant Physiol. 137: 190–198.
    OpenUrlAbstract/FREE Full Text
  148. ↵
    1. von Saint Paul, V.,
    2. Zhang, W.,
    3. Kanawati, B.,
    4. Geist, B.,
    5. Faus-Kessler, T.,
    6. Schmitt-Kopplin, P.,
    7. Schäffner, A.R.
    (2011). The Arabidopsis glucosyltransferase UGT76B1 conjugates isoleucic acid and modulates plant defense and senescence. Plant Cell 23: 4124–4145.
    OpenUrlAbstract/FREE Full Text
  149. ↵
    1. Waddington, C.H.
    (1942). Canalization of development and the inheritance of acquired characters. Nature 150: 563–565.
    OpenUrlCrossRef
  150. ↵
    1. Waddington, C.H.
    (1959). Canalization of development and genetic assimilation of acquired characters. Nature 183: 1654–1655.
    OpenUrlCrossRefPubMed
  151. ↵
    1. Warnes, G.R.,
    2. Bolker, B.,
    3. Bonebakker, L.,
    4. Gentleman, R.,
    5. Liaw, W.H.A.,
    6. Lumley, T.,
    7. Maechler, M.,
    8. Magnusson, A.,
    9. Moeller, S.,
    10. Schwartz, M.,
    11. Venables, B.
    (2016). gplots: Various R Programming Tools for Plotting Data. https://CRAN.R-project.org/package=gplots.
  152. ↵
    1. Wasternack, C.
    (2007). Jasmonates: an update on biosynthesis, signal transduction and action in plant stress response, growth and development. Ann. Bot. 100: 681–697.
    OpenUrlCrossRefPubMed
  153. ↵
    1. Weinig, C.
    (2000). Plasticity versus canalization: population differences in the timing of shade-avoidance responses. Evolution 54: 441–451.
    OpenUrlCrossRefPubMed
  154. ↵
    1. Williamson, B.,
    2. Tudzynski, B.,
    3. Tudzynski, P.,
    4. van Kan, J.A.
    (2007). Botrytis cinerea: the cause of grey mould disease. Mol. Plant Pathol. 8: 561–580.
    OpenUrlCrossRefPubMed
  155. ↵
    1. Windram, O., et al
    . (2012). Arabidopsis defense against Botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis. Plant Cell 24: 3530–3557.
    OpenUrlAbstract/FREE Full Text
  156. ↵
    1. Wu, Y.,
    2. Zhang, D.,
    3. Chu, J.Y.,
    4. Boyle, P.,
    5. Wang, Y.,
    6. Brindle, I.D.,
    7. De Luca, V.,
    8. Després, C.
    (2012). The Arabidopsis NPR1 protein is a receptor for the plant defense hormone salicylic acid. Cell Rep. 1: 639–647.
    OpenUrlCrossRefPubMed
  157. ↵
    1. Xie, D.X.,
    2. Feys, B.F.,
    3. James, S.,
    4. Nieto-Rostro, M.,
    5. Turner, J.G.
    (1998). COI1: an Arabidopsis gene required for jasmonate-regulated defense and fertility. Science 280: 1091–1094.
    OpenUrlAbstract/FREE Full Text
  158. ↵
    1. Xu, L.,
    2. Liu, F.,
    3. Lechner, E.,
    4. Genschik, P.,
    5. Crosby, W.L.,
    6. Ma, H.,
    7. Peng, W.,
    8. Huang, D.,
    9. Xie, D.
    (2002). The SCF(COI1) ubiquitin-ligase complexes are required for jasmonate response in Arabidopsis. Plant Cell 14: 1919–1935.
    OpenUrlAbstract/FREE Full Text
  159. ↵
    1. Yan, S.,
    2. Dong, X.
    (2014). Perception of the plant immune signal salicylic acid. Curr. Opin. Plant Biol. 20: 64–68.
    OpenUrlCrossRefPubMed
  160. ↵
    1. Yang, L.,
    2. Li, B.,
    3. Zheng, X.Y.,
    4. Li, J.,
    5. Yang, M.,
    6. Dong, X.,
    7. He, G.,
    8. An, C.,
    9. Deng, X.W.
    (2015). Salicylic acid biosynthesis is enhanced and contributes to increased biotrophic pathogen resistance in Arabidopsis hybrids. Nat. Commun. 6: 7309.
    OpenUrlCrossRefPubMed
  161. ↵
    1. Yeh, Y.H.,
    2. Panzeri, D.,
    3. Kadota, Y.,
    4. Huang, Y.C.,
    5. Huang, P.Y.,
    6. Tao, C.N.,
    7. Roux, M.,
    8. Chien, H.C.,
    9. Chin, T.C.,
    10. Chu, P.W.,
    11. Zipfel, C.,
    12. Zimmerli, L.
    (2016). The Arabidopsis malectin-like/LRR-RLK IOS1 is critical for BAK1-dependent and BAK1-independent pattern-triggered immunity. Plant Cell 28: 1701–1721.
    OpenUrlAbstract/FREE Full Text
  162. ↵
    1. Yu, H.,
    2. Kim, P.M.,
    3. Sprecher, E.,
    4. Trifonov, V.,
    5. Gerstein, M.
    (2007). The importance of bottlenecks in protein networks: correlation with gene essentiality and expression dynamics. PLOS Comput. Biol. 3: e59.
    OpenUrlCrossRefPubMed
  163. ↵
    1. Zander, M.,
    2. Thurow, C.,
    3. Gatz, C.
    (2014). TGA transcription factors activate the salicylic acid-suppressible branch of the ethylene-induced defense program by regulating ORA59 expression. Plant Physiol. 165: 1671–1683.
    OpenUrlAbstract/FREE Full Text
  164. ↵
    1. Zarei, A.,
    2. Körbes, A.P.,
    3. Younessi, P.,
    4. Montiel, G.,
    5. Champion, A.,
    6. Memelink, J.
    (2011). Two GCC boxes and AP2/ERF-domain transcription factor ORA59 in jasmonate/ethylene-mediated activation of the PDF1.2 promoter in Arabidopsis. Plant Mol. Biol. 75: 321–331.
    OpenUrlCrossRefPubMed
  165. ↵
    1. Zhang, B.,
    2. Horvath, S.
    (2005). A general framework for weighted gene co-expression network analysis. Stat. Appl. Genet. Mol. Biol. 4: Article17.
  166. ↵
    1. Zhang, C.,
    2. He, Y.,
    3. Zhu, P.,
    4. Chen, L.,
    5. Wang, Y.,
    6. Ni, B.,
    7. Xu, L.
    (2015a). Loss of bcbrn1 and bcpks13 in Botrytis cinerea not only blocks melanization but also increases vegetative growth and virulence. Mol. Plant Microbe Interact. 28: 1091–1101.
    OpenUrl
  167. ↵
    1. Zhang, F., et al
    . (2015b). Structural basis of JAZ repression of MYC transcription factors in jasmonate signalling. Nature 525: 269–273.
    OpenUrlCrossRefPubMed
  168. ↵
    1. Zhang, W.,
    2. Kwon, S.T.,
    3. Chen, F.,
    4. Kliebenstein, D.J.
    (2016). Isolate dependency of Brassica rapa resistance QTLs to Botrytis cinerea. Front. Plant Sci. 7: 161.
    OpenUrl
  169. ↵
    1. Zhou, N.,
    2. Tootle, T.L.,
    3. Glazebrook, J.
    (1999). Arabidopsis PAD3, a gene required for camalexin biosynthesis, encodes a putative cytochrome P450 monooxygenase. Plant Cell 11: 2419–2428.
    OpenUrlAbstract/FREE Full Text
  170. ↵
    1. Züst, T.,
    2. Heichinger, C.,
    3. Grossniklaus, U.,
    4. Harrington, R.,
    5. Kliebenstein, D.J.,
    6. Turnbull, L.A.
    (2012). Natural enemies drive geographic variation in plant defenses. Science 338: 116–119.
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top

Table of Contents

Print
Download PDF
Email Article

Thank you for your interest in spreading the word on Plant Cell.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Plastic Transcriptomes Stabilize Immunity to Pathogen Diversity: The Jasmonic Acid and Salicylic Acid Networks within the Arabidopsis/Botrytis Pathosystem
(Your Name) has sent you a message from Plant Cell
(Your Name) thought you would like to see the Plant Cell web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Plastic Transcriptomes Stabilize Immunity to Pathogen Diversity: The Jasmonic Acid and Salicylic Acid Networks within the Arabidopsis/Botrytis Pathosystem
Wei Zhang, Jason A. Corwin, Daniel Copeland, Julie Feusier, Robert Eshbaugh, Fang Chen, Susana Atwell, Daniel J. Kliebenstein
The Plant Cell Nov 2017, 29 (11) 2727-2752; DOI: 10.1105/tpc.17.00348

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Plastic Transcriptomes Stabilize Immunity to Pathogen Diversity: The Jasmonic Acid and Salicylic Acid Networks within the Arabidopsis/Botrytis Pathosystem
Wei Zhang, Jason A. Corwin, Daniel Copeland, Julie Feusier, Robert Eshbaugh, Fang Chen, Susana Atwell, Daniel J. Kliebenstein
The Plant Cell Nov 2017, 29 (11) 2727-2752; DOI: 10.1105/tpc.17.00348
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • METHODS
    • Acknowledgments
    • AUTHOR CONTRIBUTIONS
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF

In this issue

The Plant Cell: 29 (11)
The Plant Cell
Vol. 29, Issue 11
Nov 2017
  • Table of Contents
  • Table of Contents (PDF)
  • Cover (PDF)
  • About the Cover
  • Index by author
View this article with LENS

More in this TOC Section

  • Alternative Crassulacean Acid Metabolism Modes Provide Environment-Specific Water-Saving Benefits in a Leaf Metabolic Model
  • ARADEEPOPSIS, an Automated Workflow for Top-View Plant Phenomics using Semantic Segmentation of Leaf States
  • Nonsense-Mediated RNA Decay Factor UPF1 Is Critical for Posttranscriptional and Translational Gene Regulation in Arabidopsis
Show more LARGE-SCALE BIOLOGY ARTICLES

Similar Articles

Our Content

  • Home
  • Current Issue
  • Plant Cell Preview
  • Archive
  • Teaching Tools in Plant Biology
  • Plant Physiology
  • Plant Direct
  • Plantae
  • ASPB

For Authors

  • Instructions
  • Submit a Manuscript
  • Editorial Board and Staff
  • Policies
  • Recognizing our Authors

For Reviewers

  • Instructions
  • Peer Review Reports
  • Journal Miles
  • Transfer of reviews to Plant Direct
  • Policies

Other Services

  • Permissions
  • Librarian resources
  • Advertise in our journals
  • Alerts
  • RSS Feeds
  • Contact Us

Copyright © 2021 by The American Society of Plant Biologists

Powered by HighWire