Received 6 February 2016; accepted 19 June 2016; published 22 June 2016
While many adaptive or learned strategies have been acquired through biological evolution, nowhere are the learning processes more acute or dynamic than at mutualistic interfaces  . The constant and rapidly changing evolutionary tensions between intimately associated organisms may be most developmentally and metabolically dynamic in the sessile parasitic plants where chemical and tactile cues direct lifecycle phases  -  . The host-derived cues that mediate both acute and chronic developmental decision-making events are known as xenognosins, from the Greek words for “other” (xenos) and “knowledge” (gnosis).
Studies on Striga asiatica, an obligate hemi-parasitic angiosperm commonly known as witchweed, have defined the chemistry underlying host/parasite transition strategies in this noxious weed. Distinct classes of xenognosins provide the necessary and sufficient stimulus to initiate at least three separate developmental transitions, each of which is crucial to the successful development of the parasite  -  (Figure 1). Seeds of S. asiatica can remain viable under soil conditions for up to twenty years, germinating only in the presence of an SXSg (Sorghum Xenognosin for Striga germination) which includes strigolactones and/or dihydrosorgoleone derived from the roots of prospective host plants, including agriculturally important crops like maize and sorghum (Figure 1(a))   . Subsequent development of the host attachment organ, the haustorium, recommits the root apical meristem in a process known as semagenesis  -  . Reactive oxygen species (ROS) are actively exuded at the parasite root tip such that contact with host cell walls oxidatively liberates p-benzoquinone (pBQ) products. These pBQs are necessary and sufficient to initiate haustorial development; arrest of vegetative growth, induction of radial swelling of the root tip, and initiation of the ectopic root (haustorial) hairs (Figure 1(b)). Given the limited resources of the parasite for vegetative growth (3 - 5 d), this process, which takes about 24 h, prevents premature commitment. Following attachment, additional xenognosins appear necessary to initiate the cells of the shoot apical meristem (SAM). While specific endogenous factors that initiate this transition have yet to be defined formally, the exogenous addition of cytokinins (Figure 1(c)) can initiate SAM development  . Endogenous hormones derived by tapping the host vascular system  may allow the plant to complete its lifecycle, transitioning back to vegetative growth to reach seed set and, ultimately, senescence (Figure 1(d) & Figure 1(e)).
Figure 1. The life-cycle of S. asiatica is divided into stages that are either dependent ((a)-(c), red letters) or independent ((d) and (e), green letters) of xenognosins: (a) germination is initiated by strigolactones or SXSg; (b) haustorial organogenesis by p-benzoquinones; (c) shoot apical meristem development by cytokinins; (d) & (e) Flowering, seed set, and seed dispersal are host independent.
Many of these xenognosins and their analogues are now readily available, so these transitions can be easily and effectively managed under laboratory conditions. Most importantly, these developmental transitions are so highly synchronized across the population that it becomes possible to explore the molecular elements associated with surmounting the checkpoints on development. Not only will an analysis of the molecular events inform downstream regulatory and functional elements associated with specific plant developmental transitions    but the mechanisms of xenognosin regulation and how these cues impact existing machinery will provide insight into the functional evolution and successful strategies of host and parasite adaptation.
Here we combine this xenognosin-directed control over development with transcriptional profiling to identify specific elements associated with the commitment to haustorial organogenesis in S. asiatica. The identified genes suggest that semagenesis and haustorial organogenesis are connected to a quorum sensing-like process which facilities the synthesis and integration of biotic signals derived from neighboring plants to manage general growth and development. We also present evidence in support of existing models that the functional origins of these processes lie in the existing defense strategies of plants, most notably the oxidative burst response  .
2. Materials and Methods
CM-H2DCFDA was purchased from Molecular Probes (Carlsbad, CA). Murashige and Skoog (MS) media was purchased from Caisson Labs (North Logan, UT). All other chemicals were obtained from Sigma Aldrich (St. Louis, MO). Columbia-0 (Col-0) wild-type Arabidopsis thaliana seeds were purchased from Lehle Seeds (Round Rock, TX). Striga asiatica seeds were obtained from the U.S. Department of Agriculture (Beltsville, MD) and all experiments were conducted under quarantine conditions for noxious plants. These experiments were conducted under the auspices of the USDA quarantine license awarded to Emory University.
2.1. S. asiatica Germination and Haustorium Development
Seeds of Striga asiatica were pre-treated and germinated as previously described using 10−9 M Strigol  . Germination was evaluated after 36 - 48 hours and at this point experiments were initiated. Assays for haustorium development were performed as previously described  .
2.2. S. asiatica Plant Cultures
Plant cultures were obtained by transferring one day old seedlings of Striga asiatica to Magenta boxes containing 1xMS medium supplemented with either: 1 mg/L 6-BA and 0.1 mg/L IAA or 1 mg/L 6-BA and 0.5 mg/L IAA. Seedlings were grown at 23˚C with a 16 h photoperiod/day. Regenerated shoots were visible after one week, and the resulting Striga plants typically flowered within 4 - 6 weeks under these conditions.
2.3. SSH Libraries of S. asiatica Expression in Response to DMBQ
Forward and reverse subtracted cDNA libraries were generated from mRNA extracted from day-old S. asiatica seedlings either treated with 10 μM DMBQ for 6 h (pre-haustorium formation) or untreated. Suppression Subtractive Hybridization (SSH) was performed using a Clontech SSH Kit. These two cDNA libraries were used as both TESTER and DRIVER to generate a complete library of targets both up and down regulated by DMBQ exposure. Sequenced hits were identified by BLAST analysis and added to the NCBI database. These hits were utilized as the source of primers for subsequent relative transcription-PCR analyses at subsequent time points in haustorial organogenesis using actin as an internal standard for amplification.
2.4. NADPH Oxidase Cloning
Total RNA was extracted from germinated seedlings of S. asiatica with the RNeasy Plant Mini Kit (Qiagen). Integrity of the isolation was analyzed by electrophoresis in a 8% formaldehyde/1.5% agarose gel. cDNA was generated using Superscript TM III Transcriptase (Invitrogen) at 50˚C. Based on the conserved NADPH and FAD binding domains of the Respiratory burst oxidases from Arabidopsis thaliana (RbohA) and Nicotiana tabacum (NtRbohA), several degenerate primers were prepared and investigated. For example, the primer pair forward 5’-GGCAYCCITTYTCWATYACITC-3’ and reverse: 5’-GGHGTIGCWCCDATICCNARWC-3’ as well as several other degenerate primer pairs successfully cloned genes from 1.5-day-old Striga seedling cDNA. PCR products were isolated and three distinct gene sequences (SaNOX1, SaNOX2, and SaNOX3) were cloned with the help of TOPO TA Cloning Kit (Invitrogen). Whole cDNA sequences of the three genes were obtained with the help of SMART™ RACE cDNA Amplification Kit (Clontech). Based on cDNA sequences, three genomic sequences were obtained from S. asiatica genomic DNA, extracted using the DNeasy Plant Mini Kit (Qiagen). The promoter region for each putative SaNOX homologue was identified by the thermal asymmetric interlaced (TAIL)-PCR method  .
2.5. Evaluating Localization and Regulation of SaNOX Expression
To analyze the expression pattern of these genes, RT-PCR was performed on 3 distinct sets of tissues: 1) roots, 2) shoots, and 3) leaves and flowers. RNA and cDNA were obtained as described above. In order to exclude genomic DNA contamination, RNA was pretreated with DNase. The RT‑PCR studies utilized two primers for each gene (SaNOX1-3) and two primers for actin as an internal standard. SaNOX1, 2, and 3 expression levels were evaluated by 28, 30, and 30 cycles respectively with an annealing temperature of 58°C. The RT-PCR experiments for each gene were repeated in triplicate with different cDNA samples. Gene specific primers used in RT-PCR were: SaNOX1 reverse: 5’-ctgcaccggacgatgactatcttagc-3’; forward: 5’-ctgctatatcatccataacgcctttg-3’ or reverse: 5’-cagatcttccgagggacgaatccgtaaaat-3’; forward: 5’-gccatgttgaatttgacggtctcggcag-3’. SaNOX2 reverse: 5’-ttgcctaagccatttgaccgcctca-3’; forward: 5’-cctaactgccttatgtgaatgctgagg-3’. SaNOX3 reverse: 5’-cgagctattggcatttcgtgttgagc-3’; forward: 5’-ccttgtggctttgacatgtgcagagcc-3’. Actin reverse: 5’-caggctgttctctccctttat-3’; forward: 5’tccgatccagacactgtactt-3’. Changes in the expression of SaNOX1, 2, or 3 in responses to quinone treatments were evaluated in a similar manner. For Northern analysis, Total RNA from specific Striga tissue was extracted as above. Approximately, 10 μg of each RNA sample was transferred to a nylon membrane, and probed with the P32 labeled cDNA sequence for SaNOX1. Samples are washed 3 times and exposed to X-ray film (Kodak).
2.6. Cloning SaNOX Promoters and Transient Arabidopsis Transformation
Isolated promoters for SaNOX1-3 were cut by SalI and XbaI, and ligated into the PBI101 vector to generate the promoter-directed β-glucuronidase (GUS) reporters. Transient expression of the reporter constructs in Arabidopsis thaliana Columbia-0 (wildtype) seedlings was facilitated via vacuum infiltration of Agrobacterium tumefaciens GV3101 strains bearing the appropriate promoter:reporter construct. Transformed seedlings were selected on 1xMS plates with Kanamycin (50 μg/ml) and subjected to GUS staining  .
2.7. GUS Assay
Arabidopsis transformants are vacuum infiltrated for 5 minutes in the staining solution (50 mM sodium phosphate buffer, pH 7.0, 0.2% triton-X-100, 1 mM X-Gluc), incubated at 37˚C, and fixed in 75% ethyl alcohol. The pictures were taken with a Canon digital camera.
3.1. Developmental Phases Can Be Synchronized by Xenognosins in S. asiatica Seedlings
Our approach depends on the synchronized population-wide response of S. asiatica to xenognosin exposure. We therefore began by considering exposure times and conditions for each of the xenognosin-mediated developmental transitions: root apical meristem initiation (germination), haustorial organogenesis, and shoot apical meristem (SAM) initiation. Previous studies established that a 6 h exposure to sub-micromolar concentrations of SXSg results in near quantitative germination of S. asiatica   . Similarly a 6 h exposure to 10 μM 2,6-di- methoxy-p-benzoquinone (DMBQ), a previously established haustorial inducer, induces population-wide terminal commitment to haustorial organogenesis  .
A specific xenognosin for shoot apical meristem initiation has yet to be identified. However, exogenous cytokinins have been shown to induce germination, haustorial organogenesis, and shoot apical meristem development in S. asiatica  . In these prior studies, the arrested dormancy of the shoot apical meristem was observed by the shedding of the seed coat, which we used as a simple reporter. In S. asiatica seedlings in which haustoria were pre-formed, a 4 h exposure to 10 μM trans-zeatin resulted in 50% seed coat shedding (t1/2) within 24 hours (Figure 2). This effect was not isolated to trans-zeatin as 10 μM treatments with 6-benzyl adenine (6-BA) or kinetin also induced population-wide seed coat shedding with t1/2 exposures of 6 h and 10 h, respectively.
Figure 2. Cytokinins induce shoot apical meristem development in S. asiatica. Pre-haustoria are formed on seedlings of S. asiatica (See Methods). These seedlings are rinsed in triplicate then placed in buffer either (■) with or (●) without 10 μM trans-zeatin at room temperature. Seedlings are scored for shedding of the seed coat, as an indicator of SAM initiation, at the indicated time point. Results are the average of three trials with error equal to +/− the standard deviation of the mean. n = 50 seedlings.
3.2. Difference Library Construction and Analysis of Regulated Genes
Knowing the exposure time required for developmental commitment allowed us to optimize library construction. Suppression subtractive hybridization (SSH) permits the identification of differential gene expression between two populations and has been exploited for studies ranging from metagenomic analyses to the responses of a single individual to exogenously added molecular signals   . In this process, mRNA is isolated from the two populations to be compared and reverse transcribed into cDNA carrying distinct adaptor sequences. dsDNA hybrids formed between members of the two libraries are removed, identifying the differential expression of specific genes.
The phenotypes that score for the initiation of germination and apical meristem initiation (seed coat shedding) provide simple yes/no readouts. The timing of the various events during haustorial organogenesis can be divided into 3 distinct stages: (I) pre-commitment (<2 h), (II) transitioning (2 - 6 h), and (III) post-commitment (8 - 24 h). Stage II is further characterized by an arrest in root elongation and Stage III by (III.a) early commitment/radial swelling, (III.b) haustorial hair initiation (16 h) and (III.c) complete organogenesis. These transitions provide well-defined intermediate time points for library construction, offering improved resolution of the changes in gene expression associated with this developmental transition. We therefore selected haustorial organogenesis as the first xenognosin-regulated transition for evaluation by SSH.
Focusing on Stage II, when mRNA is extracted from day-old S. asiatica seedlings treated with 10 μM DMBQ for 6 h and compared with untreated seedlings, 34 unique hits are identified. Estimated fold-change analysis of these products via relative transcription-Polymerase Chain Reaction (rt-PCR), with actin as an internal standard, indicated that 5 of these products show <2-fold change in expression and are therefore excluded from further analysis. Initial assignments of the remaining 29 products is accomplished via BLAST analysis utilizing the NCBI database, and organized into separate groups based on their proposed functions: 1) Cell wall expansion and vascular tissue development; 2) plant defense; 3) nutrient metabolism and transport; and 4) hormone regulation (Table 1). In addition, two previously identified expansins of S. asiatica  as well as a homolog of the quinone oxidoreductase QR2 from Triphysaria versicolor  , a facultative root parasite in the Orobanche (broomrape) family, are identified in this library. The observed changes in expression appear specific to xenognostic quinones as the non-inducing di-tert-butyl-benzoquinone has no discernable effect on gene expression and tetrafluorobenzoquinone (TFBQ), a reversible inhibitor of haustorium development, blocks the differential expression mediated by DMBQ  .
The differential expression of these genes over the 24 h period of haustorium development, relative to un-
Table 1. Sequences obtained from the DMBQ-induced haustorium library.
a: NCBI ID of deposited full-length sequence determined from cDNA; b: Matching sequence identified through BLAST analysis; c: Time seeds were harvested for library construction; d: Expression-levels are based on relative transcription PCR to untreated seedlings with actin used as internal standard to normalize results; e: Time point for initial SSH library construction. ND = Not Determined.
treated (DMBQ-free) seedlings of similar age by the same rt-PCR analysis method, is summarized in Table 1. Specifically, these results are categorized and discussed in terms of the individual developmental stages:
I. Pre-commitment (0.5 h). Of the 29 genes identified in our library, 12 (40%) are differentially expressed within 0.5h after DMBQ exposure. This group includes two known expansins (3 and 4)  , as well as a newly identified expansin (2). Up-regulation of these proteins may be correlated with remodeling of the root tip. However, some of the most dramatic changes in expression at this early time point are associated with ROS production. This group includes the down regulation of an NADPH oxidase (10) as well as peroxidases 13 and 15. The significant up-regulation (>7-fold) of a QR2 homologue (17) is particularly noteworthy given previous studies in T. versicolor that implicate QR2 as a potential receptor for xenognostic p-benzoquinones  .
Also apparent are changes in the expression of genes associated with metabolism and transport of nutrients as well as hormone regulation. A putative auxin regulatory protein, 29, supports phytohormone regulation early during haustorial organogenesis. A role for auxin at this early time point is also consistent with the arrest of root elongation observed upon DMBQ exposure prior to commitment. Phytohormone studies in T. versicolor have also implicated a role for auxin signaling in haustorial organogenesis  . As cytokinins induce haustorium development, the developmental transition may also be regulated by both of these often antagonistic hormones  .
II. Transitioning (2 - 6 h). Differential regulation of the remaining 17 (≈ 60%) genes in the library largely occurs throughout the period of commitment to haustorium development (2 - 6 h). The expansins (2-4) remain up regulated, while Elongation factor 1 (5) and α-tubulin (6) are both significantly down regulated, suggesting some relationship between the rate of protein synthesis and microtubule assembly in regulating haustorium development   . We note that α-tubulin knockdowns in A. thaliana are characterized by radial expansion of the root tip, arrested root elongation, and ectopic root hair formation adopting a structure morphologically similar to the haustorium  . Regarding plant defense responses, the NADPH oxidase (10), peroxidases (13 and 15), and QR2 (17) remain differentially regulated throughout this period as well. Calmodulin (11) expression increases during this time period implicating Ca2+-mediated signaling events in commitment and/or development. Some genes associated with metabolism and transport of nutrients as well as hormone regulation (19, 21, 24, and 29) remains differentially regulated during this time period. Two cytokinin associated response regulators (26 and 27) are also upregulated during this period.
III. Post-commitment (16 - 24 h). Finally, the rt-PCR assays of library members at 16 and 24h post-DMBQ exposure show no significant variation in the expression relative to the 6h treatment. These results are consistent with the molecular elements required for haustorium development being largely in place after 6h, coupling expression to the clock for terminal commitment  .
3.3. Library Guided Identification of Genes for Semagenic ROS Production
The connection to defense and wounding is the most dominant and striking feature of all three stages. ROS production is central to the semagenesis mechanism and likely an early and important element in the functional emergence of parasitic angiosperms. Production is tightly regulated by the presence of xenognostic quinines  , which may be due to a reduction in the expression of a ROS-producing catalyst. Likely candidates for this catalyst in plants include xanthine oxidases  , peroxidases  , lipoxygenases  , polyamine oxidases  , mitochondrial NADH dehydrogenases  , oxalate oxidases  , and NADPH oxidases  -  . The SSH library identified two peroxidases (13 and 15) and one NADPH oxidase (10) down regulated by DMBQ exposure during pre-commitment (0.5 h, stages I) and transition (2 - 6 h, stage II).
Several possible inhibitors for ROS synthesis were evaluated including Salicylhydroxamic acid (SHAM), AEBCF, Apocynin, Phenylarsine oxide (PAO), Diphenylene iodonium (DPI), Oxypurinol, Allopurinol, Rotenone, MDL7257, and Nordihydroguaiartic acid  -  . Of these, only diphenylene iodonium (DPI) and phenylarsine oxide (PAO) inhibits ROS production in S. asiatica seedlings (Supplementary Figure S1), most consistent with a superoxide (O2•−) generating NADPH oxidase (NOX), like 10, as the catalyst  -  . The putative Striga asiatica NADPH oxidase 1 (10) is referred to as SaNOX1 for the remainder of the discussion.
Degenerate oligonucletoide primers for the conserved regions of NOX proteins successfully identified two additional putative NOX homologues: SaNOX2 and SaNOX3 (accession numbers DQ431678 and DQ431677)   . Distinct promoter sites for each of these genes were obtained using thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR)  . Exon-intron analyses of the full-length genomic sequences, as well as structural predictions based on the mRNA, are available in the Supplementary Information and summarized below (Supplementary Figures S2-S4).
NADPH oxidases in plants are often referred to as respiratory burst oxidase homologue (Rboh) proteins due to their initial discovery as the source of the oxidant in this immune response   . Consistent with known Rboh proteins, analysis of the SaNOX1 sequence supports the presence of well-conserved structural features including (from N à C): 1) 2 Ca2+-binding EF-hand domains, 2) 6 transmembrane helices, 3) corresponding histidines for heme-binding on helices III and V, 4) an FAD-binding domain, and 5) an NADPH-binding domain   . In contrast, SaNOX2 and SaNOX3 appear to be N-terminal truncations that begin in the loop between transmembrane helix I and II. All three proteins have distinct, homologous sequences with unique exon-intron patterns.
Northern blot and PCR-based expression analysis of SaNOX1-3 in S. asiatica plants grown in culture confirm SaNOX2 & SaNOX3 are transcribed throughout the plant, while SaNOX1 expression is limited to the roots (Figure 3(a) and Figure 3(b)). In order to further refine SaNOX1-3 localization, Arabidopsis thaliana was employed as a heterologous host expression system. Transient expression via Agrobacterium-mediated vacuum infiltration was used to identify specific sites of SaNOX expression. SaNOX2::GUS and SaNOX3::GUS transformants accumulate stain throughout the plant (Figure 3(c) and Figure 3(d)) while in the SaNOX1::GUS staining is limited to the root tip (Figure 3(e) and Figure 3(f)), the corresponding site of ROS production in S. asiatica.
Prior to the commitment to haustorial organogenesis, the removal of pBQs restores semagenic ROS production  , and as shown in Figure 4, the expression of SaNOX1 as well. Additional haustorial inducing pBQs, like the unsubstituted p-benzoquinone or methyl-pBQ also attenuate ROS production as well as SaNOX1 expression, further supporting a connection to semagenesis (data not shown).
Figure 3. Tissue localization and expression of SaNOX1-3. Tissues harvested from germinated seedlings (root) or mature cultured S. asiatica plants were (a) screened for SaNOX1-3 localization of expression by relative transcription PCR using actin as a control; (b) Expression of SaNOX1 and SaNOX3 were confirmed by Northern analyses using 18s rRNA as a control. Lane labels are R: Root, L + S: Leaf + Stem; and F: Flower. (c)-(f) A. thaliana mutant lines expressing β-glucoronidase downstream of the promoters for SaNOX2 (c), SaNOX3 (d), and SaNOX1 (e) & (f) were generated to establish expression localization of each product.
Figure 4. SaNOX1 expression is reversible prior to commitment to haustorial organogenesis. SaNOX1 expression was evaluated via RT-PCR at the indicated time points following the addition of 10 μM DMBQ. After 2 h of DMBQ exposure seedlings were washed in triplicate and placed in either buffer (grey) or returned to a solution of 10 μM DMBQ (white). Auxin was used as an internal standard for normalization.
The successful integration of a eukaryotic parasite with its host presents many of the same synchronization challenges faced during multicellular tissue development. The temporal and spatial control over cell fate defines tissue architecture, and the mechanisms underlying these precise responses appear to be similarly constrained between two distinct organisms. Defining the temporal and spatial organization of the chemistry that makes this control possible has significant practical and theoretical consequences for plant development as well as the adaptations necessary for both pathogenic and mutualistic symbioses.
Here we have focused our attention on the haustorium, the host attachment organ whose development is common across the various manifestations of the parasitic strategy in plants   . In S. asiatica, as well as other parasitic plants, the detection of simple p-benzoquinones is both necessary and sufficient to mobilize the genetic elements necessary for the rapid functional transition to the haustorium. Our initial differential expression library identified 29 genes that were regulated during commitment to haustorial organogenesis. These genes fall into the general categories of cell wall expansion and vascular tissue development, nutrient metabolism and transport, hormone regulation, and cellular defense (Table 1).
The haustorium provides the vascular conduit between host and parasite in the parasitic plants, and its development begins with the radial swelling and ectopic root hair formation necessary for host attachment. The genes associated with cell wall expansion (e.g., expansins) as well as vascularization (e.g., chorismate mutase) are reasonable candidates for differential expression during haustorial organogenesis. The down-regulation of α- tubulin expression is of particular interest given that A. thaliana knockdowns of this protein give root phenotypes similar to that observed in haustorial organogenesis: radial expansion of the root tip, arrested root elongation, and ectopic root hair formation   . The general up-regulations of genes associated with nutrient metabolism and transport are consistent with the principal function of the haustorium: tapping host resources. The regulated flux of auxins and cytokinins is probably the most expected of a developmental transition and consistent with previous studies of haustorial organogenesis  .
Conversely, the changes in defense gene expression seem contradictory. On one hand, the grafting of two tissues through the haustorium must involve the regulation of xenognostic processes, while at the same time constitutive ROS production in the tissue undergoing haustorial organogenesis is central to the chemical mechanism of semagenesis. The SSH library and subsequent PCR analyses confirmed the down regulation of two putative sources of semagenic ROS, a peroxidase (15) and an NADPH oxidase (SaNOX1) (10). Small molecule inhibitors further implicated the NADPH oxidase (SaNOX1) as the source of ROS at the root tip that was regulated in response to DMBQ. Furthermore, heterologous expression in A. thaliana confirmed SaNOX1 expression only at the root tip, the site of haustorial development.
Like other respiratory burst oxidase homologues (Rboh) in plants, structural predictions of SaNOX1 indicate all the features required for calcium-regulated ROS production are conserved in this protein. Based on these predictions, as well as its observed regulation by pBQ concentrations, we suspect that SaNOX1 provides the source of semagenic ROS. Direct support for the redox activity of SaNOX1, however, that may require expression in a heterologous host as knockout/knockdown experiments are currently limited in S. asiatica.
As for SaNOX2 and SaNOX3, while they appear to be N-terminal truncations of these Rboh proteins, they retain all of the features associated with ROS production and may be well expressed in their corresponding tissues. The apparent absence of EF-hand domains for calcium regulation is unusual in plants but is not uncommon in other eukaryotes  . The variations in their sequence as well as exon-intron utilization argue against their arising simply from recent gene duplication events. Indeed comparative exon-intron analysis between the eight Rboh proteins from Arabidopsis thaliana confirm distinct exon-intron patterns for several of those proteins as well, although none display truncations like those in SaNOX2 or SaNOX3   . Such unusual truncations may reflect the accelerated rate at which mutations appear to accumulate among parasitic angiosperms, relative to their non-parasitic relatives  . These increased rates of mutation may well provide an adaptive advantage to parasites in dealing with the evolution of resistance in hosts.
In addition to their well-established roles in defense, ROS production has also been implicated in plant growth and development including root hair formation and pollen tube growth  -  . In this context, the redox events associated with semagenesis may not be unique to haustorial organogenesis, but more generally associated with cell-cell signaling and development. Increasing evidence implicates the semagenic quinones as not only regulating haustorial development in the parasitic plants  , but more broadly in root meristem growth and development. The addition of semagenic quinones impacts development of both root architecture and notably root hair density in A. thaliana and N. tabaccum as well as the production of ROS   .
Root growth and architecture is sensitive to the density and identity of neighboring roots (self and non-self)   and the chemistry exploited by S. asiatica in semagenesis could serve a more general quorum sensing- type role in plant roots, controlling morphological development (phenotype) based on their proximal neighbors. In this context, understanding the integration of accumulating pBQs via semagenesis in S. asiatica may provide greater insight into how plants ‘decide’ to organize the architecture of their root systems. Indeed, the strigolactones, which were originally discovered as germination stimulants for Striga spp. (Figure 1(a)), have now been shown to function as regulators of branching in both plants and arbuscular mychorizzal fungi   . These findings bolster the model that xenognosin-regulated transitions in parasites likely arose from pre-existing regulatory pathways present among non-parasites. Given the tight regulation of these transitions, these systems now serve as valuable model systems for defining the molecular machinery responsible for signal integration, response to environmental inputs, and decision-making in critical developmental commitments. Such information can be used to regulate the “decision making” processes and drive productive/lethal transitions in plants. This information then is critical to reading and controlling the molecular code of plant environmental sensing and adaptation to an ever-changing environmental landscape.
We thank Eric Weeks (Emory Department of Physics) for support with confocal microscopy.
Figure S1. Effect of ROS inhibitors on CM-H2DCFDA assay. Two day old seedlings of Striga asiatica are incubated with the indicated concentration of inhibitor for 2 hours then loaded with 10 μM CM-H2DCFDA for 5 minutes and imaged for ROS production (See Methods). The ROS source for each inhibitor is indicated in parentheses behind each inhibitor. (a) □ Salicylhydroxamic acid (peroxidases)  , о AEBCF (Animal not plant NOX proteins), rApocynin (Animal not plant NOX)  ; (b) Phenylarsineoxide (PAO) (NOX)  ; (c) Diphenyleneiodonium (DPI) (NOX)  ; (d) ■ Oxypurinol (Xanthine oxidases)  , ● Allopurinol (Xanthine oxidases)  , ▲ Rotenone (Mitochondrial complex I)  , ? MDL 72,527 (polyamine oxidase)  , ♦Nordihydroguaiartic acid (lipoxygenase)  . Fluorescence intensity expressed relative to Striga seedlings treated with CM-H2DCFDA only (set at 100). Results are the average of three trials with error equal to +/− the standard deviation of the mean.
Figure S2. Proposed structure of SaNOX1-3 based on BLAST analysis, CLUSTALW alignments, and published literature on the structure of gp91phox and the Arabidopsis thaliana respiratory burst oxidase A (Rboh A). Both of these proteins produce superoxide anions but only Rboh A is regulated be calcium due to the presence of the Ca2+-binding EF-hand domains which are also conserved in SaNOX1 (underlined) but not SaNOX 2 or 3. In addition the first transmembrane helix found in SaNOX1, Rboh A, and gp91phox (H1 and bold text) is absent in SaNOX2 and 3. Near the C-terminus of the proteins are the highly conserved FAD and NADPH binding domains (BD) indicated by the single line or double lines respectively.
Figure S3. Results of SPIDEY analysis using mRNA and genomic sequences for: SaNOX1, 2, and 3 respectively (http://www.ncbi.nlm.nih.gov/IEB/Research/Ostell/Spidey/).
Figure S4. Exon-intron analysis of SaNOX1, 2 and 3. Genomic sequences of SaNOX1-3 are organized by exons and introns, based on the SPIDEY analysis above (See Supplementary Table 1). The first exons of SaNOX2 and 3 are positioned to reflect the exon of SaNOX1 to which they show the greatest similarity (based on CLUSTALW analysis). Position is indicated by vertical dashed lines Sequences are drawn to scale.
*Authors contributed equally to this work.
 Hillesland, K.L. and Stahl, D.A. (2010) Rapid Evolution of Stability and Productivity at the Origin of a Microbial Mutualism. Proceedings of the National Academy of Sciences of the United States of America, 107, 2124-2129.
 Kim, D., Kocz, R., Boone, L., Keyes, W.J. and Lynn, D.G. (1998) On Becoming a Parasite: Evaluating the Role of Wall Oxidases in Parasitic Plant Development. Chemistry & Biology, 5, 103-117.
 Keyes, W.J., Malley, R.C.O., Kim, D. and Lynn, D.G. (2000) Signaling Organogenesis in Parasitic Angiosperms: Xenognosin Generation, Perception and Response. Journal of Plant Growth Regulation, 19, 217-231.
 Palmer, A.G., Gao, R., Maresh, J., Erbil, W.K. and Lynn, D.G. (2004) Chemical Biology of Multi-Host/Pathogen Interactions: Chemical Perception and Metabolic Complementation. Annual Review of Phytopathology, 42, 439-464.
 Lynn, D.G. and Chang, M. (1990) Phenolic Signals in Cohabitation: Implications for Plant Development. Annual Reviews of Plant Physiology and Plant Molecular Biology, 41, 497-526.
 Kim, D., Kocz, R., Boone, L., Keyes, W.J. and Lynn, D.G. (1998) On Becoming a Parasite: Evaluating the Role of Wall Oxidases in Parasitic Plant Development. Chemistry and Biology, 5, 103-117.
 Keyes, W.J., Palmer, A.G., Erbil, W.K., Taylor, J.V., Apkarian, R.P., Weeks, E.R., et al. (2007) Semagenesis and the Parasitic Angiosperm Striga asiatica. Plant Journal, 51,707-716.
 Palmer, A.G., Liu, Y., Adkins, S.M., Zhang, X., Wu, I., Chen, X., et al. (2008) The Molecular Language of Semagenesis. Plant Signaling and Behavior, 3, 560-561.
 Worsham, A.D., Moreland, D.E. and Klingman, G.C. (1959) Stimulation of Striga asiatica (Witchweed) Seed Germination by 6-Substituted Purines. Science, 130, 1654-1656.
 Smith, C.E., Ruttledge, T., Zeng, Z., O’Malley, R.C. and Lynn, D.G. (1996) A Mechanism for Inducing Plant Development: The Genesis of a Specific Inhibitor. Proceedings of the National Academy of Sciences of the United States of America, 93, 6986-6991.
 Tomilov, A., Tomilova, N., Shin, D.H., Jamison-McClung, D., Torres, M., Reagan, R., McGray, H., Horning, T., Truong, R., Nava, A.J., Nava, A. and Yoder, J.I. (2006) Chemical Signaling between Plants: Mechanistic Similarities between Allelopathy and Host Plant Recognition by Parasitic ANGIOSPERMS. In: Dicke, M. and Takken, W., Eds., Chemical Ecology: From Genes to Ecosystem, Springer, Berlin, 55-69.
 Esfandiari, E., Jin, Z., Abdeen, A., Griffiths, J.S., Western, T.L. and Haughn, G.W. (2013) Identification and Analysis of an Outer-Seed-Coat-Specific Promoter from Arabidopsis thaliana. Plant Molecular Biology, 81, 93-104.
 Tague, B.W. and Mantis, J. (2006) In Planta Agrobacterium-Mediated Transformation by Vacuum Infiltration. In: Salinas, J. and Sanchez-Serrano, J.J., Eds., Arabidopsis Protocols, Humana Press, New York, 215-223.
 Galbraith, E.A., Antonopoulos, D.A. and White, B.A. (2008) Application of Suppressive Subtractive Hybridization to Uncover the Metagenomic Diversity of Environmental Samples. In: Martin, C.C., Ed., Environmental Genomics, Humana Press, New York, 295-334.
 Mahalingam, R., Gomez-Buitrago, A., Eckardt, N., Shah, N., Guevara-Garcia, A., Day, P., et al. (2003) Characterizing the Stress/Defense Transcriptome of Arabidopsis. Genome Biology, 4, R20.
 Matvienko, M., Wojtowicz, A., Wrobel, R., Jamison, D., Goldwasser, Y. and Yoder, J.I. (2001) Quinone Oxidoreductase Message Levels Are Differentially Regulated in Parasitic and Non-Parasitic Plants Exposed to Allelopathic Quinones. The Plant Journal, 25, 375-387.
 Bandaranayake, P.C.G., Filappova, T., Tomilov, A., Tomilova, N.B., Jamison-McClung, D., Ngo, Q., et al. (2010) A Single-Electron Reducing Quinone Oxidoreductase Is Necessary to Induce Haustorium Development in the Root Parasitic Plant Triphysaria. The Plant Cell, 22, 1404-1419.
 Tomilov, A.A., Tomilova, N.B., Abdallah, I. and Yoder, J.I. (2005) Localized Hormone Fluxes and Early Haustorium Development in the Hemiparasitic Plant Triphysaria versicolor. Plant Physiology, 138, 1469-1480.
 Bao, Y., Kost, B. and Chua, N.H. (2001) Reduced Expression of Alpha-Tubulin Genes in Arabidopsis thaliana Specifically Affects Root Growth and Morphology, Root Hair Development and Root Gravitropism. The Plant Journal, 28, 145-157.
 Ursin, V.M., Irvine, J.M., Hiatt, W.R. and Shewmaker, C.K. (1991) Developmental Analysis of Elongation Factor-1 Alpha Expression in Transgenic Tobacco. The Plant Cell, 3, 583-591.
 Montillet, J.L., Cacas, J.L., Garnier, L., Montané, M.H., Douki, T., Bessoule, J.J., et al. (2004) The Upstream Oxylipin Profile of Arabidopsis thaliana: A Tool to Scan for Oxidative Stresses. Plant Journal, 40, 439-451.
 Jarvinen, A., Grigorenko, N., Khomutov, A.R., Hyvonen, M.T., Uimari, A., Vepsalainen, J., et al. (2005) Metabolic Stability of α-Methylated Polyamine Derivatives and Their Use as Substitutes for the Natural Polyamines. The Journal of Biological Chemistry, 280, 6595-6601.
 Moller, I.M. (2001) PLANT MITOCHONDRIA AND OXIDATIVE STRESS: Electron Transport, NADPH Turnover, and Metabolism of Reactive Oxygen Species. Annual Review of Plant Physiology and Plant Molecular Biology, 52, 561-591.
 Gucciardo, S., Wisniewski, J.P., Brewin, N.J. and Bornemann, S. (2007) A Germin-Like Protein with Superoxide Dismutase Activity in Pea Nodules with High Protein Sequence Identity to a Putative Rhicadhesin Receptor. Journal of Experimental Botany, 58, 1161-1171.
 Foreman, J., Demidchik, V., Bothwell, J.H.F., Mylona, P., Miedema, H., Torres, M.A., et al. (2003) Reactive Oxygen Species Produced by NADPH Oxidase Regulate Plant Cell Growth. Nature, 422, 442-446.
 Jabs, T., Tschope, M., Colling, C., Hahlbrock, K. and Scheel, D. (1997) Elicitor-Stimulated Ion Fluxes and O2- from the Oxidative Burst Are Essential Components in Triggering Defense Gene Activation and Phytoalexin Synthesis in Parsley. Proceedings of the National Academy of Sciences of the United States of America, 94, 4800-4805.
 Rea, G., de Pinto, M.C., Tavazza, R., Biondi, S., Gobbi, V., Ferrante, P., et al. (2004) Ectopic Expression of Maize Polyamine Oxidase and Pea Copper Amine Oxidase in the Cell Wall of Tobacco Plants. Plant Physiology, 134, 1414-1426.
 Jaquet, V., Scapozza, L., Clark, R.A., Krause, K.-H. and Lambeth, J.D. (2009) Small-Molecule NOX Inhibitors: ROS-Generating NADPH Oxidases as Therapeutic Targets. Antioxidants & Redox Signaling, 11, 2535-2552.
 Abramov, A.Y., Jacobson, J., Wientjes, F., Hothersall, J., Canevari, L. and Duchen, M.R. (2005) Expression and Modulation of an NADPH Oxidase in Mammalian Astrocytes. The Journal of Neuroscience, 25, 9176-9184.
 Das, D.K., Engelman, R.M., Clement, R., Otani, H., Prasad, M.R. and Rao, P.S. (1987) Role of Xanthine Oxidase Inhibitor as Free Radical Scavenger: A Novel Mechanism of Action of Allopurinol and Oxypurinol in Myocardial Salvage. Biochemical and Biophysical Research Communications, 148, 314-319.
 Le Cabec, V. and Maridonneau-Parini, I. (1995) Complete and Reversible Inhibition of NADPH Oxidase in Human Neutrophils by Phenylarsine Oxide at a Step Distal to Membrane Translocation of the Enzyme Subunits. Journal of Biological Chemistry, 270, 2067-2073.
 Matsunaga, T., Kudo, J., Takahashi, K., Dohmen, K., Hayashida, K., Okamura, S., et al. (1996) Rotenone, a Mitochondrial NADH Dehydrogenase Inhibitor, Induces Cell Surface Expression of CD13 and CD38 and Apoptosis in HL-60 Cells. Leukemia & Lymphoma, 20, 487-494.
 Pivovarov, A.S. and Egido-Villareal, W. (1995) The Influence of an Inhibitor of Lipoxygenases on the Modulation of the Plasticity of Cholinoreceptors by 15-Hete. Neuroscience and Behavioral Physiology, 25, 483-487.
 Potocky, M., Jones, M.A., Bezvoda, R., Smirnoff, N. and Zársky, V. (2007) Reactive Oxygen Species Produced by NADPH Oxidase Are Involved in Pollen Tube Growth. New Phytologist, 174, 742-751.
 Torres, M.A., Onouchi, H., Hamada, S., Machida, C., Hammond-kosack, K.E. and Jones, J.D.G. (1998) Six Arabidopsis thaliana Homologues of the Human Respiratory Burst Oxidase (gp91phox). The Plant Journal, 14, 365-370.
 Yoshioka, H., Numata, N., Nakajima, K., Katou, S., Kawakita, K., Rowland, O., et al. (2003) Nicotiana benthamiana gp91phox Homologs NbrbohA and NbrbohB Participate in H2O2 Accumulation and Resistance to Phytophthora infestans. The Plant Cell, 15, 706-718.
 Barkman, T.J., McNeal, J.R., Lim, S.-H., Coat, G., Croom, H.B., Young, N.D., et al. (2007) Mitochondrial DNA Suggests at Least 11 Origins of Parasitism in Angiosperms and Reveals Genomic Chimerism in Parasitic Plants. BMC Evolutionary Biology, 7, 248.
 Bromham, L., Cowman, P.F. and Lanfear, R. (2013) Parasitic Plants Have Increased Rates of Molecular Evolution across All Three Genomes. BMC Evolutionary Biology, 13, 126.
 Palmer, A.G., Chen, M.C., Kinger, N.P. and Lynn, D.G. (2009) Parasitic Angiosperms, Semagenesis and General Strategies for Plant-Plant Signaling in the Rhizosphere. Pest Management Science, 65, 512-519.
 Bhatt, M.V., Khandelwal, A. and Dudley, S.A. (2011) Kin Recognition, Not Competitive Interactions, Predicts Root Allocation in Young Cakile edentula Seedling Pairs. New Phytologist, 189, 1135-1142.
 Besserer, A., Puech-Pagès, V., Kiefer, P., Gomez-Roldan, V., Jauneau, A., Roy, S., et al. (2006) Strigolactones Stimulate Arbuscular Mycorrhizal Fungi by Activating Mitochondria. PLoS Biology, 4, 1239-1247.
 Dun, E.A., Brewer, P.B. and Beveridge, C.A. (2009) Strigolactones: Discovery of the Elusive Shoot Branching Hormone. Trends in Plant Science, 14, 364-372.