Multiple immunity-related genes control susceptibility of Arabidopsis thaliana to the parasitic weed Phelipanche aegyptiaca

Plant growth and quantifying parasite attachment and development

The parasite attachment assay was conducted using a randomized incomplete block design due to large number of tested genotypes. For each experimental block, four mutant plant lines were randomly selected from the pool of all mutant lines in the collection. Approximately 25 A. thaliana seeds of each of the four mutant lines plus the wild type background ecotype were stratified in water at 4 °C for 2 days in then planted in Sunshine #1 potting mix (Sungro, Agawam, MA, USA). Plants were grown in a Conviron ATC40 growth chamber at 20 °C, 12-h light cycle and light intensity of 90 µmol m²s⁻¹ for 10 days. Eight Polyethylene (PE) bags with glass fiber grade A (GFA) paper (Whatman, Maidstone, UK) backings were made as previously described (Westwood, 2000) with dimensions of 26 cm × 9 cm (Fig. 1A). The 10-day-old A. thaliana seedlings were gently removed from the soil and the roots were washed with water until free of soil. The plants were transplanted such that the roots were positioned between the PE bag and the GFA paper and the hypocotyl extended from the top of the bag (Fig. 1A). One plant of each of the four mutant lines and one plant of the wild type background ecotype were randomly distributed into each of the PE bags and placed under a laboratory growth light with a 12-h light cycle for 12 days. P. aegyptiaca seeds were sterilized following previously described protocols (Westwood, 2000) and placed on wet GFA paper in a sealed Petri plate (Fisher Scientific, Waltham, MA, USA) for six days to condition the seeds (Westwood, 2000). The strigolactone analog rac-GR24 was then applied to the conditioned P. aegyptiaca seeds at a concentration of 2 mg/L to stimulate germination. The Petri dishes were resealed and stored in the dark for an additional 24 h. Conditioned and stimulated P. aegyptiaca seeds were inoculated onto the roots of the A. thaliana plants in the PE bags (12 days after transplanting from soil) using a fine-tipped paintbrush. By placing germinated P. aegyptiaca seeds immediately adjacent to developed roots, the impact of differential root development among the various tested Arabidopsis mutants was limited. The seeds were aligned to be within 0.3 mm of roots that had grown since the transplanting (lighter colored roots) and approximately 35 seeds were placed along the roots of each plant in each bag.

The attachment rate for each individual A. thaliana plant was quantified 15 days after inoculation. Every P. aegyptiaca seed was examined under a dissecting microscope (Zeiss Stemi SV11) at 12×g magnification and classified as either non-germinated, germinated but not attached, or attached. Not attached and attached germinated seeds were distinguished by probing with a 0.3 mm wide dissecting probe to test for adherence of the parasite radicle to the A. thaliana roots. The attached P. aegyptiaca plants were further classified as either early-stage attachment/haustoria connection, early-stage tubercle, or late-stage tubercle which correspond to the stages 3, 4.1 and 4.2 in the Parasitic Plant Genome Project datasets (Westwood et al., 2012; Yang et al., 2014) (Fig. 1D). The 0.3 mm dissecting probe was used to distinguish the three attached classifications as follows: an attached radicle thinner than 0.3 mm was categorized as an attached radicle/initial attachment (stage 3). A tubercle formed that is thicker than 0.3 mm but has no secondary roots was categorized as an early-stage tubercle (stage 4.1) and coincides with the completion of the vascular connection/feeding bridge between the parasite and the host plant. A tubercle with secondary roots longer than 0.3 mm was categorized as a late-stage tubercle (stage 4.2). The rate of attachment for each of the three stages was determined by calculating the ratio of all of the P. aegyptiaca seeds that had reached at least the designated stage (i.e., the rate of initial attachment included the count of late-stage tubercles plus early-stage tubercles) divided by the total number of germinated P. aegyptiaca seeds (i.e., the total number of all stages plus those seeds that germinated but did not attach to hosts).

Statistical analysis of attachment rates

We used a generalized linear mixed model (SAS/STAT(R) 9.2 User’s Guide SE, 2016) to analyze different rates of parasitization at each of the three developmental stages. This model is based on the binomial distribution with the logit link. The SAS script and example files are uploaded to the Ag Data Commons (https://data.nal.usda.gov/dataset/data-multiple-immune-pathways-control-susceptibility-arabidopsis-thaliana-parasitic-weed-phelipanche-aegyptiaca). The one fixed effect of interest was treatment of the plants. The treatments are the different host plant genetic pathways which were compared to the respective wildtype in our model (Col-0, Ler-0, or Ws-2). Both bag and position were set as fixed effects to account for any systematic variability within both bag and position. The experiment was specified as a random block, which allows us to determine significance of the fixed effects above and beyond experiment to experiment variability. This also accounts for the fact that results within an experiment may be more related than results across experiments. Our response variable was the ratio of attachment stage total to overall germinated seed total (described above). We compared 95% and 99% confidence intervals on the odds ratio estimates at each of the three stages to determine which treatments are statistically significantly different from the wildtype. An odds ratio greater than one indicates that P. aegyptiaca is more successful on the mutant line relative to the wildtype. An odds ratio less than one indicates that P. aegyptiaca is less successful on the mutant line than on the wildtype.

Type III tests of fixed effects indicated that there was no significant variation dependent on the bag. However, there was a significant effect dependent on the position of the plant within the bag for early and late stage tubercle attachment rates (Table S1). Tukey–Kramer analysis of the least square means of the differences of position with adjustment for multiple comparisons showed that significant differences (p < 0.05) always included plants on the edge of the bag (positions 1 and 5). We hypothesize that these differences were due to the edge of the GFA paper drying out more unevenly than the center of the GFA paper. Therefore, it is essential to randomize the position of the plants within each bag, as was done here.

Time-lapse photography

A PE bag containing A. thaliana eco. Col-0 was inoculated with stimulated P. aegyptiaca seeds as described above. The bags were mounted on a metal hanger inside a 11.3-L Rubbermaid container that had a 70 mm diameter hole cut in the side for insertion of a macro camera lens and the top portion was covered with aluminum foil with a narrow opening so that the A. thaliana shoots could be exposed to light (12-h light-dark cycles) while the parasites were in darkness. The camera was a Sigma SD14 digital camera with a 70 mm Sigma DG macro lens. The camera flash provided the light source for the photography. The A. thaliana was watered through a tube that fed directly into the bag, so that the plant did not have to be disturbed during the time-lapse photography. For recording aboveground plant growth, a tomato (Solanum lycopersicum) was used as host because Arabidopsis does not support robust aboveground growth in our experimental setup. A Nikon D5100 camera with a Nikon DX AF-S Nikkor lens was used. Pictures were taken at 9-min intervals for 37 days over the course of the above and below ground segments. Sigma Photo Pro software (Sigma Corporation of America, Ronkonkoma, NY, USA) was used to convert the images from raw X3F format into JPEG. Adobe Premiere Pro CC 2017 was then used to assemble the images and for video editing to make the time-lapse video at 24 fps.

Quantitative RT-PCR

Three biological replicates of approximately 40 0.5 cm root sections were harvested from both inoculated and time-matched mock-inoculated A. thaliana eco. Col-0 plants at early attachment stage, early vascularization/early tubercle stage, and late tubercle stage (six total experimental conditions with three biological replicates and three technical replicates each). RNA was extracted using a Qiagen RNeasy kit. cDNA was synthesized using Superscript cDNA synthesis kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s recommended protocol. qRT-PCRs were performed by StepOnePlus Real-Time PCR System (Invitrogen, Carlsbad, CA, USA) with Power SYBR green master mix (Invitrogen, Carlsbad, CA, USA). The reaction condition was 2 min at 50 °C and 10 min at 95 °C followed by 40 cycles at 95 °C for 15 s, at 56 °C for 15 s and 72 °C for 15 s. Dissociation curves were evaluated to confirm the specificity. Relative quantification was calculated by 2^−∆∆Ct method using qBase+ software. Three housekeeping genes (SAND, UBQ10 and UFP) were used for the normalization. geNorm analysis revealed that geNorm M-value of all three genes (SAND, UBQ10 and UFP) were 0.762, 0.709 and 0.619, relatively (qBASE+, Biogazelle). The recommended cutoff value is one or less. Calibrated, normalized relative quantification values for each gene were compared between the mock-inoculated and Phelipanche-inoculated plants for each parasite developmental stage using a two-tailed, two-sample unequal variance t-test. The primers used are described in Table S2. All primers were previously designed in other studies: (Jacobs et al., 2011) (VSP2, MYB51, EXPPT1, BOI, SID2, WRKY53), (Paponov et al., 2008) (IAA13, IAA2, ACS6), (Nguyen et al., 2016) (ARR10, UBQ10), (Šašek et al., 2014) (SAND), (Tran, Chen & Wang, 2017) (UFP), (Pegadaraju et al., 2007) (PAD4), (Zhang et al., 2014) (PR1). All primers were validated to be specific to A. thaliana DNA and to not amplify any product from cDNA prepared from early tubercle stage P. aegyptiaca tissue.

A quantitative screen for host plant resistance and susceptibility to parasitic weeds

We developed a quantitative assay to measure the susceptibility of A. thaliana–P. aegyptiaca attachment and tubercle development based on previously described techniques (Westwood, 2000). A. thaliana plants grown in PE bags with GFA backing (Fig. 1A) allow for host roots to be directly inoculated with P. aegyptiaca. All P. aegyptiaca seeds were stimulated with the strigolactone analog rac-GR24 (Yoneyama et al., 2010), which is an effective germination stimulant for P. aegyptiaca (Fig. 1B). P. aegyptiaca attachment and development can be monitored in this growth system. Our analysis of time-lapse video of P. aegyptiaca (Video S1—available at https://figshare.com/articles/PhelipancheParasitization_mov/11894388) growth highlighted the acceleration in parasite growth following development of the floral meristem. Therefore, the challenge for our assay of parasite growth on mutants was to devise a method that could provide sensitive, reproducible measures of parasite success during a developmental stage characterized by relatively small changes in growth. The number of P. aegyptiaca seedlings at various growth stages can be counted under a dissecting microscope (Fig. 1C). Graphical representation of the typical attachment and tubercle development rates from inoculation of P. aegyptiaca onto wildtype A. thaliana is shown in Fig. 1D. All data from all parasite attachment experiments are in Table S3. Because of the substantial number of comparisons made in this analysis (three different stages for each of 46 mutant lines compared to wildtype), we primarily focus on differences that met the more stringent 99% confidence limit cutoff (see “Methods”) to minimize false positives. Based on this cutoff, only between 2 and 9 mutant lines showed altered susceptibility to P. aegyptiaca parasitization at any specific stage compared to wildtype A. thaliana (Table S4).

Host genetic pathways essential for susceptibility to P. aegyptiaca

Phelipanche aegyptiaca is a successful pathogen of A. thaliana (Goldwasser, Plakhine & Yoder, 2000); therefore, our assay of attachment and development rates is more suited to identify mutations that attenuate A. thaliana susceptibility to P. aegyptiaca parasitization than mutations that increase susceptibility. The 46 tested mutant lines include numerous disruptions in key signaling pathways (Table 1). It is important to note that many of these genes are involved in several different genetic pathways due to the substantial crosstalk among many of the signaling pathways. Examples include pad4 and eds1 which are involved in salicylic acid (SA) signaling, but also critical for several R gene-mediated responses (Cui et al., 2017; Rustérucci et al., 2001) and the substantial crosstalk between jasmonic acid (JA), SA and ethylene signaling (Koornneef & Pieterse, 2008).

For the stage of initial parasitization/haustoria development, no mutant lines were less susceptible than wildtype A. thaliana at the 99% confidence level (Fig. 2). This result demonstrates that none of the tested host endogenous signaling pathways are essential for P. aegyptiaca to form haustoria and attach to host roots. However, plants overexpressing Ethylene Response Factor 2 (ERF2), a transcription factor upregulated in the presence of ethylene that activates numerous immune responses (Catinot, Huang & Zimmerli, 2015), were significantly less susceptible than wildtype to P. aegyptiaca initial attachment at the 95% confidence level, suggesting that elevated host ethylene levels may attenuate initial attachment. The attachment rates were more variable within and among experiments for the early attachment stage compared to the two later stages measuring tubercle growth (Table S3).

Eight A. thaliana mutant lines were attenuated in the ability to support P. aegyptiaca early-stage tubercle formation at the 99% confidence limit (Fig. 3). Five of the less susceptible mutants involve knockouts in jasmonic acid biosynthesis or signaling (dde2, dde2/pad4, dde2/sid2, jar1/mlo, jar1/axr1) suggesting that P. aegyptiaca requires functional host JA signaling and biosynthesis to form successful tubercles. However, neither the jar1 nor jin1 single mutant is less susceptible to P. aegyptiaca tubercle formation. pad4 plants were also less susceptible to P. aegyptiaca tubercle development. PAD4 is involved in both SA signaling and TIR R-gene signaling. Because other mutants defective in SA accumulation (sid2 and NahG) and R-gene-mediated responses (sgt1, rar1 and rpm1) were not significantly less susceptible to parasite attachment, the reason pad4 plants are less susceptible at the early tubercle stage remains uncertain. Lastly, the mutants cpr5 and edr1 were substantially less susceptible to early-stage tubercle formation. Both of these mutants have been previously shown to constitutively express multiple plant defense pathways and be more resistant to several microbial biotrophic plant pathogens (Bowling et al., 1997; Christiansen et al., 2011; Tang, Christiansen & Innes, 2005; Yoshida et al., 2002).

Six of the mutants (jar1/axr1, jar1/mlo, cpr5, dde2/pad4, dde2 and dde2/sid2) identified as being less susceptible to early-stage tubercles also supported significantly fewer late-stage P. aegyptiaca tubercles at the 99% confidence limit (Fig. 4). Additionally, the jaz4/5/9 mutant, which is defective in JA responses because the JAZ4, JAZ5 and JAZ9 proteins are all transcriptional regulators triggered by jasmonate (Pauwels & Goossens, 2011; Yan et al., 2014), also supported fewer late-stage tubercles than WT A. thaliana. The JA biosynthesis mutant jar1 supported fewer late stage tubercles at the 95% confidence level giving further support to the hypothesis that several genes involved in host JA signaling and biosynthesis are essential for P. aegyptiaca parasitism.

In addition to the JA signaling and biosynthesis mutant lines, the DELLA quad mutant (gai/rga/rgl1/rgl2) and rdr6 mutant lines supported fewer late-stage tubercles at the 99% confidence level. The DELLA genes (GAI, RGR and RGLs) are negative growth regulators controlled by Gibberellic acid (GA), and this quad mutant is more resistant to the biotrophic pathogen Pseudomonas syringae and more susceptible to the necrotrophic pathogen Alterneria brassisicola (Navarro et al., 2008).

Host genetic pathways essential for resistance to P. aegyptiaca

Two tested mutants, mpk4 and pfd6-1, were more susceptible to initial P. aegyptiaca parasitization (Fig. 2). PFD6 is a subunit of the prefoldin complex, and pfd6 mutants are defective in microtubule function and formation (Gu et al., 2008) and is a putative hub of the A. thaliana immune system (Mukhtar et al., 2011). pfd6-1 mutants were also significantly more susceptible to supporting early-stage tubercle development by P. aegyptiaca (Fig. 3). We, therefore, hypothesize that PFD6 is an important component of the host plant immune response against plant parasitization.

None of the tested mutant lines supported significantly more late-stage P. aegyptiaca tubercles compared to wildtype A. thaliana at the 99% confidence level (Fig. 4). However, pfd6 did support significantly more late-stage tubers at the 95% confidence level. Additionally, the triple mutant arr1/arr10/arr11, which is deficient in cytokinin signaling (Mason et al., 2005), also supported more late-stage P. aegyptiaca tubercles at the 95% confidence level.

Parasitization by P. aegyptiaca alters the transcription of several immunity-related genes

We monitored transcript levels of 12 defense-and hormone signaling-associated genes of A. thaliana (Table S2) following P. aegyptiaca parasitization. Significant upregulation of the glucosinolate biosynthetic gene MYB51 and downregulation of the negative cell death regulator BOI were observed during initial haustorial attachment (Fig. 5). During early stage tubercle development, the pattern-trigged immunity (PTI) marker gene WRKY 53 is significantly upregulated. During late-stage tubercle development, only smaller magnitude changes in gene expression were observed, most notably significant down-regulation in the JA marker gene VSP2 and the SA marker gene SID2. Whether these changes are more directly associated with host defense responses or parasite virulence strategies remains unknown.

JA/SA signaling

Several host genes involved in JA biosynthesis and signaling are (e.g., dde2, jar1, jaz4/5/9) are required for full host susceptibility to development of parasite tubercles (Figs. 3 and 4). Several of these mutants also supported fewer early-stage haustorial attachments (Fig. 2), but not significantly as determined by our stringent cutoff (Table S4). Therefore, we conclude that aberrations in host endogenous JA signaling can compromise compatibility to P. aegyptiaca parasitization. Further supporting this hypothesis, JA signaling genes were previously identified as upregulated in response to compatible parasitization of Lotus japonicus by P. aegyptiaca (Hiraoka, Ueda & Sugimoto, 2008), A. thaliana by O. ramosa (Dos Santos et al., 2003), sorghum by Striga (Hiraoka & Sugimoto, 2008) and Medicago truncatula by Orobanche crenata (Die et al., 2007). The common interpretation of these studies is that hosts induce JA-related defense responses to deter the parasite, but all of these hosts are susceptible to the respective parasites. Exogenous jasmonate treatment was previously shown to limit P. aegyptiaca parasitism (Bar-Nun & Mayer, 2008), suggesting that general perturbations in JA signaling may attenuate the compatibility of host plants to root parasitization.

Jasmonic acid signaling plays a critical role in numerous components of plant response to pathogen attack, including antagonism to SA signaling (Gutjahr & Paszkowski, 2009; Kazan & Manners, 2008). A reasonable hypothesis is that the JA mutants are less susceptible to parasitism because of elevated SA levels. High levels of SA are often associated with reduced virulence of biotrophic plant pathogens (Thomma et al., 2001). Application of exogenous SA was previously shown to reduce parasitism of clover by the closely related parasite Orobanche minor (Kusumoto et al., 2007). However, several SA-deficient mutants (e.g., sid2 and NahG) were not more susceptible to parasite attachment and development. Additionally, a few mutants with known defects in SA accumulation and/or signaling (ndr1 and pad4) were actually less susceptible to early-and late-stage tubercle development at the 95% confidence interval. Taken together, these results suggest that perturbations to functional crosstalk between SA and JA signaling limit the ability of P. aegyptiaca to parasitize A. thaliana. Therefore, we hypothesize that wildtype-functional SA and JA signaling is required by the parasite, potentially because these pathways are precisely targeted and manipulated during parasitism to support the attachment and development of the parasite. This hypothesis is further supported by the observation that the DELLA quad mutant plants are significantly less susceptible to late-stage tubercle development. The DELLA quad mutant was previously shown to be much more susceptible to biotrophs, likely due to aberrations in SA–JA crosstalk (Navarro et al., 2008). The DELLA genes are repressed by GA and parasitized tissue likely contains low levels of GA based on the expression of the marker gene EXP-PT1 (Fig. 5).

Putative hubs and regulators of plant immunity

pfd6 plants were significantly more susceptible to P. aegyptiaca parasitization than wildtype host plants (Figs. 2–4). Even though our assay is biased toward identifying mutant host plants that are less susceptible to P. aegyptiaca parasitization, the increased susceptibility of pfd6 was one of the strongest and most consistent phenotypes identified. PFD6 is a putative immunity hub protein (Mukhtar et al., 2011) involved in microtubule dynamics (Gu et al., 2008). It is possible that microtubules are being deployed to prevent parasitization, and disruption of microtubule formation is sufficient to support increased parasitization. Several plant pathogens target microtubules as part of their attack (Hardham, 2013; Lee et al., 2012). Overexpression of PFD6 or genes involved in similar functions may be sufficient to increase the resistance of host plants against parasitic plants.

Obligate holoparasites, such as P. aegyptiaca, are analogous to biotrophic microbial plant pathogens because they rely on living host tissue for the completion of their lifecycle. Multiple mutants that have hyper-activated immunity and are resistant to several biotrophic plant pathogens (cpr5 and edr1) were substantially more resistant to P. aegyptiaca as well. This finding suggests that increased activation of known immune pathways that function in the control of biotrophic microbial plant pathogens is sufficient, in some instances, to also limit the severity of infestation by parasitic plants.

Auxin/cytokinin

The arr1/10/11 triple mutant supported more development of late-stage P. aegyptiaca tubercles than wildtype A. thaliana (Fig. 4) but only at the 95% confidence level. This finding suggests that host plant insensitivity to cytokinins leads to increased parasite attachment. The related parasite Phtheirospermum japonicum was recently shown to translocate cytokinins into host tissue to alter host development (Spallek et al., 2017). The reduced parasite attachment rates to axr1 plants (Figs. 3 and 4) potentially contradict the hypothesis that altered cytokinin signaling enhances host susceptibility to P. aegyptiaca. Though initially identified as an essential component of Arabidopsis response to auxin treatment, AXR1 is also involved in sensitivity to ethylene (Timpte et al., 1995), cytokinin (Li, Kurepa & Smalle, 2013) and methyl-jasmonate (Tiryaki & Staswick, 2002). axr1 mutant plants are expected to have lower sensitivity to cytokinins due to increased stability of the ARR5 protein, a negative regulator of cytokinin responses (Li, Kurepa & Smalle, 2013). We hypothesize that the reduced susceptibility to parasitization of the axr1 A. thaliana genotype is due to reduced sensitivity to ethylene (see below) as opposed to reduced sensitivity to cytokinins. Alternatively, reduced sensitivity to auxin signaling may explain the reduced susceptibility of axr1 mutants to P. aegyptiaca. Local auxin biosynthesis and response gene activation were recently implicated as essential for haustorial formation by the P. japonicum (Ishida et al., 2016). However, the auxin signaling mutants tir1 and aux1 were not differentially susceptible to formation of P. aegyptiaca tubercles.

Ethylene signaling

The ethylene signaling mutant lines, including ein2, did not support differential rates of P. aegyptiaca attachment and development compared to wildtype A. thaliana. Nevertheless, several observations suggest that ethylene signaling is an important part of host plant responses to P. aegyptiaca parasitism. First, host plants overexpressing Erf2, were substantially less susceptible to initial haustorial attachment by P. aegyptiaca (Fig. 2). Second, the dde2 mutant line (JA biosynthesis) exhibited substantially reduced susceptibility to P. aegyptiaca tubercle development, but the quad signaling mutant line (Tsuda et al., 2009) (dde2/pad4/sid2/ein2) was not less susceptible to parasitization. This result suggests that loss of function in SA signaling or ethylene signaling is sufficient to rescue the reduced susceptibility associated with defective host plant JA biosynthesis. Wildtype level P. aegyptiaca attachment rates were not rescued when infecting the double mutants dde2/pad4 and dde2/sid2 suggesting that the rescued attachment rate is due to the ein2 mutation in the dde2 background. We were not able to directly test this hypothesis because dde2/ein2 double mutant plants do not grow in the PE bag system for unknown reasons. Third, the resistant-to-parasitism mutant edr1 has perturbed ethylene signaling, specifically in the crosstalk of SA and ethylene signaling (Tang, Christiansen & Innes, 2005). However, the edr1 phenotype may be more predominant due to EDR1’s role as a negative regulator of plant immunity (see above). Given the high degree of crosstalk between JA and SA signaling with ethylene signaling (Koornneef & Pieterse, 2008), it will be challenging to determine the extent to which both of these host hormone-signaling pathways are independently required for successful P. aegyptiaca parasitization.

RNAi

RDR6 is involved in RNAi signaling and rdr6 mutant plants are less susceptible to late-stage P. aegyptiaca tubercle formation, suggesting that P. aegyptiaca may exploit host-derived RNAi pathways during parasitization. The distantly related parasitic plant Cuscuta campestris, for example, secretes sRNA molecules into host plant tissue to silence plant-immunity and auxin-signaling gene pathways (Shahid et al., 2018). However, other mutants deficient in RNAi signaling (dcl1and rdr1/2/6) are not less susceptible to P. aegyptiaca. Indeed, both dcl1 and rdr1/2/6 are among the most susceptible mutants tested in terms of late-stage tubercle development (Fig. 4). DCL1 is essential for most miRNA biogenesis (Bologna & Voinnet, 2014) and plays a role in the negative regulation in an RDR6-involved RNA-silencing pathway (Qu, Ye & Morris, 2008). It is possible that P. aegyptiaca exploits this RNA-silencing pathway during parasitization, hence dcl1 (negative regulator) mutants are highly susceptible and rdr6 mutants are significantly less susceptible to tubercle development. However, the differential susceptibility of the rdr6 and the rdr1/2/6 lines cannot be explained at this time. A more thorough examination is necessary to determine the role of RNAi in Arabidopsis–P. aegyptiaca interactions.

Pathways that appear not to be involved in susceptibility/resistance to P. aegyptiaca parasitization

It is important to note that the majority of tested mutant A. thaliana genotypes were not significantly different from wildtype in their susceptibility to P. aegyptiaca parasitization (Table S4). This result suggests that the ability of P. aegyptiaca to parasitize host plant tissue is robust, which is also supported by the relatively broad host range of P. aegyptiaca (Parker, 2013). Surprisingly, every tested A. thaliana mutant line was able to support some level of parasite attachment. In contrast, a recent similar screen of the obligate stem parasite Cuscuta reflexa on multiple tomato introgression lines revealed multiple genotypes where the host-parasite interaction was incompatible (Krause et al., 2018). We consider the lack of increased susceptibility of bak1-4 plants to be the most surprising result. We originally hypothesized that this mutant was the best candidate for a genotype more susceptible to parasitization because BAK1 is critical for many components of PTI (Chinchilla et al., 2009), which is a critical part of the plant immune response against many other classes of plant pests (Jones & Dangl, 2006) including the parasitic plant C. reflexa (Hegenauer et al., 2016). Future studies with the bak1-5 allele that specifically blocks innate immunity without impairing brassinosteroid signaling (Schwessinger et al., 2011) may be informative. mpk4 plants, which are deficient in map kinase signaling—another critical component of PTI (Rasmussen et al., 2012), were more susceptible to early attachment than wildtype (Fig. 2). But the increased susceptibility was not persistent through early and late stage P. aegyptiaca tubercle development. It is possible that P. aegyptiaca is already highly adept at subverting A. thaliana PTI. However, we show that the PTI marker gene WRKY53 is significantly upregulated in response to tubercle development (Fig. 5). Several of the borderline increased susceptibility genotypes, such as bak1-4, should be tested in the PE bag system using a related parasite such as O. minor that is largely incompatible on A. thaliana (Goldwasser, Plakhine & Yoder, 2000). Starting with an incompatible interaction will more readily enable discovery of host genes that underpin resistance to parasitism than the assay described here performed with the compatible P. aegyptiaca parasite.