An Innate Immunity Pathway in the Moss Physcomitrella patens[OPEN]

A P. patens signaling pathway required for immunity triggered by PAMPs induces growth inhibition, a novel fluorescence burst, cell wall depositions, and accumulation of defense-related transcripts. MAP kinase (MPK) cascades in Arabidopsis thaliana and other vascular plants are activated by developmental cues, abiotic stress, and pathogen infection. Much less is known of MPK functions in nonvascular land plants such as the moss Physcomitrella patens. Here, we provide evidence for a signaling pathway in P. patens required for immunity triggered by pathogen associated molecular patterns (PAMPs). This pathway induces rapid growth inhibition, a novel fluorescence burst, cell wall depositions, and accumulation of defense-related transcripts. Two P. patens MPKs (MPK4a and MPK4b) are phosphorylated and activated in response to PAMPs. This activation in response to the fungal PAMP chitin requires a chitin receptor and one or more MAP kinase kinase kinases and MAP kinase kinases. Knockout lines of MPK4a appear wild type but have increased susceptibility to the pathogenic fungi Botrytis cinerea and Alternaria brassisicola. Both PAMPs and osmotic stress activate some of the same MPKs in Arabidopsis. In contrast, abscisic acid treatment or osmotic stress of P. patens does not activate MPK4a or any other MPK, but activates at least one SnRK2 kinase. Signaling via MPK4a may therefore be specific to immunity, and the moss relies on other pathways to respond to osmotic stress.

Activated PRR signals are transmitted and amplified within minutes by MAP kinase cascades. These involve activation by phosphorylation of a MAP kinase kinase kinase (MEKK), a MAP kinase kinase (MKK), and a MAP kinase (MPK), which phosphorylates substrate proteins including transcription factors . At least four Arabidopsis MPKs (MPK3/4/6 and 11) are activated by PRRs (Fiil et al., 2009;Bethke et al., 2012), and Arabidopsis MPK1 and 2 are also implicated in immunity to the fungal pathogen Botrytis cinerea (Weiberg et al., 2013). In addition to their roles in Arabidopsis immunity, MPK3/4/6 have been implicated in responses to abiotic stress including osmotic, salt, and cold stresses (Droillard et al., 2004;Teige et al., 2004;Rodriguez et al., 2010). However, the contributions of single MPKs in Arabidopsis immunity or stress responses remain unclear due to the pleiotropic phenotypes of some of their lossof-function mutants and to their potential redundancies (Petersen et al., 2000;Ren et al., 2008;Zhang et al., 2012).
The Arabidopsis genome encodes 20 MPKs, including MPK4, which functions in a cascade(s) with MEKK1 and MKK1 and 2 (Ichimura et al., 1998;Qiu et al., 2008). In response to pathogens or PAMPs, this cascade permits the WRKY33 transcription factor to enhance the expression of defensive genes including PAD3 required for synthesis of antimicrobial camalexin (Suarez- Rodriguez et al., 2007;Qiu et al., 2008). The Arabidopsis mpk4 loss-of-function mutants exhibit dwarfism and ectopically express defense responses but are suppressed by mutation of SUMM2, which encodes a resistance (R) protein (Zhang et al., 2012). This indicates that the MPK4 pathway is a target of pathogen effectors and is monitored by SUMM2, in line with the guard hypothesis of effector-triggered immunity in plant-pathogen evolution (Jones and Dangl, 2006;Zhang et al., 2012).
The moss Physcomitrella patens is a bryophyte, an early branch of land plants (Zimmer et al., 2007). P. patens responds to pathogens, including fungi, oomycetes, and bacteria, by activating defenses including production of reactive oxygen species (ROS), cell wall fortifications, hormonal perturbations, and expression of defenserelated genes (Lehtonen et al., 2012;Mittag et al., 2015;Oliver et al., 2009;Ponce de León et al., 2007. However, PTI components have not been identified in P. patens or other bryophytes. The P. patens genome encodes only eight MPKs (Dóczi et al., 2012). They represent a basal set of pathways whose functions may be assessed by their targeted knockouts (KOs) or knockins (KIs) via homologous recombination, a genetic approach unavailable in vascular plants. Comparative analyses of such moss lines may elucidate aspects of the complexity of MPK activation in flowering plants and how plant immunity and stress responses have evolved.
To initiate such analyses, we examined if P. patens MPKs respond to PAMPs or abiotic stresses by assaying MPK phosphorylation and in-gel kinase activity. This showed that at least two MPKs are rapidly phosphorylated and activated in response to bacterial and fungal PAMPs. Targeted gene KOs and KIs identified components of this pathway, including a PRR, at least one MEKK and one MKK, and two MPKs required for responses to the fungal PAMP chitin and resistance to two fungal pathogens. In contrast, osmotic stress did not result in detectable MPK phosphorylation or activation in the moss but did result in rapid activation of at least one SnRK2 kinase of a canonical abscisic acid response pathway.

Phosphorylation and Activation of P. patens MPKs in Response to Chitin
We first assayed MPK phosphorylation in P. patens protein extracts using an anti-p44/42-ERK antibody (a-pTEpY) following treatment with bacterial and fungal PAMPs. MPK phosphorylation was undetectable following treatment with bacterial elongation factor Tu peptide (elf18; Roux et al., 2011), conserved flagellin peptides (Pseudomonas flg22; Gómez-Gómez and Boller, 2000; Salmonella flgII-28; Clarke et al., 2013) or full-length flagellins from Pseudomonas or Borrelia (Supplemental Figures 1A and 1B). This is in contrast to the elicitation patterns of at least three MPKs in Arabidopsis by these PAMPs, excluding Borrelia flagellin and Salmonella flgII-28 used as negative controls for Arabidopsis. The flagellin and elf18 insensitivity of P. patens is in line with the absence of close homologs of their respective Arabidopsis receptors FLS2 and EFR in P. patens (Boller and Felix, 2009). In contrast, treatment with the fungal PAMP chitin or its derivative chitosan induced phosphorylation and activation, assayed by phosphorylation of the general kinase substrate MBP, of two P. patens MPKs within 1 min (Supplemental Figures 1A and 1C). These findings are in line with those of Lehtonen et al. (2012Lehtonen et al. ( , 2014 showing that chitosan triggers responses in P. patens including a ROS burst, defense transcript accumulation, and secretion of pathogen response proteins. In addition, these MPKs were rapidly phosphorylated in response to bacterial peptidyl glycans (PGNs; Supplemental Figure 1D), which have molecular patterns similar to chitin. This response is in line with studies showing that the Arabidopsis and rice (Oryza sativa) chitin receptors are also involved in PGN detection (Willmann et al., 2011;Ao et al., 2014). We conclude that at least two P. patens MPKs are activated in response to fungal and bacterial PAMPs.

The Moss Chitin Receptor CERK1
P. patens encodes four homologs of the Arabidopsis chitin receptor At CERK1 (Zhang et al., 2009), and the mRNA levels of two of these moss genes (Pp1s364 10V6.1 and Pp1s41 280V6.1) were enhanced by chitosan treatment (Supplemental Figures 2A and 2B and Supplemental Data Set 1). This may be relevant as stimuli that affect the activities of regulatory proteins may alter the levels of their corresponding mRNAs (Hruz et al., 2008). Targeted KO lines of these two P. patens CERK1 homologs were generated (Supplemental Figures 3A to 3C) and tested for PTI responses. Plants were transferred on cellophane sheets from normal media to media containing chitin to assess their responses. In contrast to the wild-type and KO lines of Pp1s41 280V6.1, chitin did not induce apparent cell wall modification readily stainable with Toluidine blue in Pp1s364 10V6.1 KO lines (hereafter designated Dcerk1; Figures  1A and 1B). This modification includes ethanol-insoluble, brownish red, cell wall-associated depositions in wild-type tissues, which become visible after 8 h of chitin treatment (Supplemental Figures  4A and 4B). Interestingly, chitin treatment of wild-type, but not Dcerk1 lines, led to complete growth inhibition within 1 min and to a strong fluorescence burst visible after 2 h, which peaked by 4 h and decayed by 8 h (Supplemental Figures 4A, 4C, and 4D, and Supplemental Movies 1 and 2). In addition, chitin, chitosan, chitin hexamer, or bacterial peptidyl glycan also did not induce phosphorylation of the two MPKs in Dcerk1 KO lines ( Figures 1C and 1D). This effect was not due to a general impairment of PTI, as the MPKs were phosphorylated in these lines treated with bacterial cell-free culture filtrate from the harpin negative strain Pectobacterium wasabiae (ex Pectobacterium carotovorum strain SCC3193) ( Figure  1E; Nykyri et al., 2012).
Dcerk1 KO lines were also reduced in two other responses to chitin compared with the wild type. First, Dcerk1 tissues failed to exhibit enhanced conductivity ( Figure 1F) using a standard assay for immunity-associated electrolyte leakage in Arabidopsis (Mackey et al., 2003). Second, Dcerk1 lines were compromised in the accumulation of mRNAs of the defense genes Phenylalanine Ammonia Lyase4 (PAL4) and Chalcone Synthase (CHS) ( Figure  1G; Oliver et al., 2009). Finally, expression of Pp CERK1 complemented the Arabidopsis cerk1 mutant (Miya et al., 2007), as measured by chitin-induced MPK phosphorylation and accumulation of WRKY22 transcripts (Supplemental Figures 5A and  5B). We conclude that Pp CERK1 is required for moss responses to chitin and bacterial peptidyl glycans and is therefore an ortholog of the Arabidopsis chitin receptor At CERK1.

Two P. patens MEKKs Function Redundantly in Responses to Chitin
Arabidopsis PRR signaling via MAP kinase cascades includes MEKK1-MKK1/MKK2-MPK4 as well as MEKK1-MKK4/MKK5-MPK3/MPK6 (Asai et al., 2002;Qiu et al., 2008). The P. patens genome encodes 22 MEKK homologs (Sun et al., 2014). From the phylogenetic relationships with 21 MEKKs from Arabidopsis (Supplemental Figure 6A and Supplemental Data Set 2), two P. patens homologs of At MEKK1 were chosen as candidates and are here referred to as Pp MEKK1a (Pp1s136 5V6.3) and Pp MEKK1b (Pp1s31 252V6.1). The mRNA levels of these two genes increase in response to chitosan, although Pp MEKK1a displayed a stronger response (Supplemental Figure 6B). Targeted KO lines of these two MEKK genes (D mekk1a and D mekk1b;Supplemental Figures 3D and 3E) were generated and tested for responses to chitin.
Dmekk1a exhibited strongly reduced chitin induced cell wall associated depositions compared with the wild type and Dmekk1b (Figures 2A and 2B). To assess potential redundancy of MEKK1a and MEKK1b, we generated the double KO Dmekk1a/b (Supplemental Figure 3F). This line, like the Dcerk1 KO line, was unable to induce cell wall depositions in response to chitin (Figures 2A  and 2B). Chitin-induced MPK phosphorylation was also reduced in both Dmekk1b and Dmekk1a and absent in Dmekk1a/b ( Figure 2C).
The same trend was seen in these lines on the levels of mRNAs of the defense-related PAL4 and CHS genes ( Figure 2D). In addition, Dmekk1a/b failed to induce the fluorescent burst when exposed to chitin, similarly to Dcerk1 (Supplemental Figure 7A and Supplemental Movie 3). We conclude that Pp MEKK1a and b have redundant functions in these chitin responses. Interestingly, the growth of Dmekk1a/b was inhibited within 1 min after chitin exposure as seen in wild-type moss (Supplemental Figure 7B and Supplemental Movie 4). This chitin response is therefore independent of MEKK1a and b.

Redundancy among MKKs in Response to Chitin
The P. patens genome encodes seven MKK homologs, three of which are most closely related to Arabidopsis MKK1, 2, and 6. Since At MKK1 and 2 are strongly implicated in immunity, we designated the three P. patens homologs Pp MKK1a (Pp1s32 219v6.1), Pp MKK1b (Pp1s106 83v6.1), and Pp MKK1c (Pp1s50 83v6.1) (Supplemental Figure 8A and Supplemental Data Set 3). The mRNAs of Pp MKK1a and b were slightly induced in response to chitosan, while Pp MKK1c mRNA was not (Supplemental Figure  8B). Targeted KO lines of these MKK genes were generated (Supplemental Figures 3G to 3I) and tested for responses to chitin. Toluidine blue staining revealed that chitin induced cell wall modification was slightly reduced in the two Dmkk1a and c KO lines compared with the wild type, but not in the Dmkk1b line ( Figures 3A and 3B). Chitin-induced MPK phosphorylation also appeared slightly reduced in Dmkk1a and b lines and further reduced in the Dmkk1c lines ( Figure 3C). In addition, the levels of mRNAs of the defense-related PAL4 and CHS showed a trend of being less induced following chitin treatment in Dmkk1a and c (D) Immunoblot assay for MAP kinase phosphorylation of wild type and Dcerk1-1 KO line in response to no treatment (N.T.), water, 0.01% acetic acid (AA), 100 mg/mL chitin (Chi), 100 mg/mL chitosan (Cho), or 100 mg/mL hexa-N-acetylchitohexaose (CH6) for 5 min. (E) Immunoblot assay for MAP kinase phosphorylation of wild type and Dcerk1-1 KO line in response to no treatment (N.T.), water, Luria broth (LB), 100 mg/mL chitin (Chi), 200 mg/mL peptidyl glycan from Staphylococcus aureus (PGN), or cell-free culture filtrate (CFC) from the harpin negative strain P. wasabiae SCC3193. (F) Tissue ion leakage assay of the wild type and Dcerk1-1 treated with water or chitin for times specified. Data are presented as percent of total conductivity as an average of four replicates each with five plants. Error bars represent SD. (G) RT-qPCR of PAL4 and CHS transcript levels in the wild type and Dcerk1-1 as fold change relative to the wild type untreated at times specified after 100 mg/mL chitin treatment. Data points represent an average of three independent experiments (three biological replicates) each with three technical replicates for each PCR reaction. Error bars represent SE of the means. compared with the wild type and Dmkk1b ( Figure 3D). Although these results do not pinpoint the involvement of a single moss MKK in chitin responses, they indicate that KOs of individual MKKs, in particular mkk1a and c, reduce responses to chitin. Combinatorial, multiple KOs of MKK1a, b, and c would be required to fully delineate their involvement in PAMP signal transduction.

MPK4a and b Activation in Response to Chitin and B. cinerea Inoculation
A phylogeny of MPKs from 10 plant species (Dóczi et al., 2012) identified eight putative MPKs in P. patens. Two of these homologs group with the Arabidopsis group B MPKs MPK4/5/ 11/12/13 (Supplemental Figure 9A and Supplemental Data Set 4; (C) MAP kinase phosphorylation in the same moss lines treated with 100 mg/mL chitin and assayed by immunoblotting with anti-phospho-p44/42 MPK antibody (a-pTEpY). Loading control is Ponceau staining of total protein.
(D) Quantitative RT-PCR of PAL4 and CHS transcript levels in the wild type and Dmekk1a-1, Dmekk1b-1, and Dmekk1a/b as fold change relative to the wild type untreated at times specified after 100 mg/mL chitin treatment. Data points represent an average of three independent experiments, each with three technical replicates. Error bars represent SE of the means. Ichimura et al., 2002). Arabidopsis MPK3/4/6/11 have all been shown to be phosphorylated and activated upon PAMP treatment (Gómez-Gómez and Boller, 2000;Bethke et al., 2012). Since the two P. patens MPKs are most similar to At MPK4 (Supplemental Figure 9A), we refer to them as Pp MPK4a (Pp1s149 39V6.1) and Pp MPK4b (Pp1s59 325V6.1). Pp MPK4a mRNA accumulated in response to chitosan within 15 min, peaked 8-fold at 2 h, and returned to basal levels by 8 h (Supplemental Figure 9B). In contrast, Pp MPK4b transcripts were only marginally affected by chitosan. Since the sizes of the two P. patens MPKs (C) MAP kinase phosphorylation in the same moss lines treated with 100 mg/mL chitin and assayed by immunoblotting with anti-phospho-p44/42 MPK antibody (a-pTEpY). Loading control is Ponceau staining of total protein.
(D) Quantitative RT-PCR of PAL4 and CHS transcript levels in the wild type and Dmkk1a-2, Dmkk1b-1, and Dmkk1c-2 as fold change relative to the wild type untreated at times specified after 100 mg/mL chitin treatment. Data points represent an average of three independent experiments each with three technical replicates. Error bars represent SE of the means. phosphorylated and activated in response to chitin (Figures 1C to 1E;Supplemental Figures 1C and 1D) match those predicted for MPK4a (42.8 kD) and MPK4b (43.5 kD), we created stable KI lines of these two MPK genes containing C-terminal GFP reporter fusions. We also created KO lines of MPK4a (Dmpk4a) but were unable to obtain stable knockout lines of MPK4b despite several attempts using different constructs. The recombined loci were verified in the three lines (Dmpk4a, MPK4a-GFP, and MPK4b-GFP; Supplemental Figures 3K to 3N).
Wild-type and selected Dmpk4a, MPK4a-GFP, and MPK4b-GFP lines were treated with chitin and protein extracts subjected to immunoblot analysis of MPK activation loop phosphorylation with a-pTEpY and to in-gel kinase assays with MBP as substrate. Both the 42.8-and 43.5-kD MPKs became phosphorylated and activated in the wild type, but the 42.8-kD band could not be detected in the Dmpk4a or MPK4a-GFP lines ( Figure 4A). In contrast, the 43.5-kD band was absent in MPK4b-GFP lines ( Figure 4A). In addition, phosphorylated MPK4a-GFP and MPK4b-GFP fusion proteins with expected sizes of ;70 kD (;43-kD MPK with ;27-kD GFP fusion) appeared after chitin treatment, although these bands were masked in the kinase assay by an unspecific band in all samples ( Figure 4A). However, kinase assays of the ;70-kD MPK-GFP fusions immunoprecipitated with an anti-GFP antibody confirmed that MPK4a-GFP and MPK4b-GFP could be activated by chitin treatment ( Figure 4B). Thus, MPK4a and MPK4b represent the 42.8-and 43.5-kD bands, respectively, and are the two MPKs phosphorylated and activated in response to chitin and other PAMPs tested ( Figures 1C to 1E).
We examined the localization of the MPK4a-GFP and MPK4b-GFP fusions as their chitin-inducible activities indicated that they were functional. Moss plants have three cell types that extend by tip growth: (1) chloronema with many large chloroplasts and transverse cell walls perpendicular to the growth axis; (2) caulonema with fewer, smaller chloroplasts and oblique transverse cell walls; and (3) rhizoids lacking chloroplasts that develop from the mature gametophyte. The MPK4a-GFP signal was strong in the apical cells of caulonemal air filaments and rhizoids in fully developed plants and less strong but readily detectable in filamentous protonemal tissue at the edge of the plant consisting of both chloronema and caulonema (Supplemental Figure 10A). When grown on media overlayed with cellophane, which promotes filamentous growth of protonema, the MPK4a-GFP signal was strongest in newly formed apical tip cells of protonemal tissue, whereas a 35S promoter-driven GFP control showed uniform expression throughout the plant (Supplemental Figure 10B). MPK4a-GFP localized to the cytoplasm and nucleus, and this pattern did not appear to change significantly following chitin treatment (Supplemental Figure 10C). This pattern of expression is similar to that of At MPK4-GFP in Arabidopsis (Andreasson et al., 2005). In contrast to MPK4a-GFP, MPK4b-GFP was undetectable. This may reflect the ;20-fold lower transcript level of MPK4b compared with MPK4a in untreated plants (Supplemental Figure 9C).
Since chitin is a component of fungal cell walls, we examined if MPK4a and MPK4b are activated in response to the fungal pathogen B. cinerea. To this end, we assayed the kinase activities and phosphorylation patterns of immunoprecipitated MPK4a-GFP and MPK4b-GFP following inoculation with B. cinerea spores (Supplemental Figure 11). MPK4a-GFP was clearly activated and phosphorylated in response to B. cinerea. In contrast, MPK4b-GFP activation and phosphorylation were less marked. The signals detected after inoculation with fungal spores were weaker than those detected after chitin treatment ( Figure 4B). This is presumably because even high spore concentrations only trigger a fraction of the PRR receptors compared with those triggered by soluble chitin or other PAMPs.

Kinase Activation in Response to Osmotic Stress
In addition to their PAMP responsiveness, At MPK3/4/6 are phosphorylated in response to abiotic stresses including hyperosmolarity and salt (Ichimura et al., 2000;Droillard et al., 2004;Teige et al., 2004). P. patens also accumulates specific transcripts and proteins in response to osmotic stress (Cuming et al., 2007;Lunde et al., 2007;Wang et al., 2008). To examine if Pp MPK4a or Pp MPK4b are activated in response to osmotic stress, we treated lines with 500 mM NaCl or 800 mM mannitol and assayed kinase activity. We did not detect activation or phosphorylation of Pp MPK4a or Pp MPK4b, or of any other MPKs, in response to these treatments ( Figure 5A).
However, both treatments activated a kinase of ;40 kD, which phosphorylated MBP. Since MBP is a substrate of many kinases, including stress-activated SnRKs , we performed in-gel kinase assays with the SnRK-specific synthetic peptide P3 (Halfter et al., 2000). This demonstrated that the P3 phosphorylation patterns by the ;40-kD kinase were the same as those of MBP in response to the NaCl and mannitol ( Figure 5A). Importantly, P3 was not phosphorylated by the MPKs activated in response to chitin, and a-pTEpY immunoblotting only detected MPK phosphorylation in response to chitin and not to NaCl or mannitol ( Figure 5A). It thus appears that one or more SnRK2s, but not MPKs, are involved in early responses to osmotic stress in P. patens. P. patens encodes four SnRK2 group III kinases with sizes of 38.5 to 40 kD (Fujita et al., 2013). The osmoregulatory phytohormone abscisic acid (ABA) activates both MPKs and SnRK2s in Arabidopsis (Lu et al., 2002;Boudsocq et al., 2004;Ortiz-Masia et al., 2007;Xing et al., 2008). We therefore assayed kinase activity and found that ABA activated the same P. patens SnRK2 kinase(s) of ;40 kD as did NaCl treatment, but not activities with sizes of MPKs ( Figure 5B). This SnRK2 activation by NaCl is dependent on ABA since it was strongly reduced in the ABA-insensitive moss line D2-1 ( Figure 5B; Tougane et al., 2010). Since the kinase assay of NaCl-treated moss indicated activation of SnRK2s, we measured the transcript levels of the four moss SnRK2 group III kinases in response to NaCl treatment. Transcript levels of all four SnRK2s responded similarly in the wild type and Dmpk4a (Supplemental Figure 12), indicating that MPK4a is not involved in signaling osmotic changes. Interestingly, SnRK2a transcripts accumulated to higher levels than those of the other three. We therefore generated a targeted KO of this kinase (Dsnrk2a; Supplemental Figure  3J). The levels of SnRK2 activated by NaCl or ABA were markedly less in Dsnrk2a lines compared with the wild type ( Figures 5C and 5D).

Characterization of Moss MPK4a in Immunity
Analysis of Dmpk4 KO lines showed that they are markedly reduced in chitin-induced cell wall-associated depositions assayed by Toluidine blue staining over a time course up to 48 h (Figures 4C and 4D; Supplemental Figure 13). To assess whether the levels of pathogen defense-related transcripts were altered in Dmpk4a, we measured the response to chitosan of five transcripts involved in defense against pathogens and/or are induced in response to chitosan in both Arabidopsis and P. patens: PAL4, CHS, ERF2, a-DOX, and LOX7 (Oliver et al., 2009;Lehtonen et al., 2012;Machado et al., 2015;Overdijk et al., 2016). Although the response to chitosan (Supplemental Figures 14A and 14B) is slower than to chitin ( Figure 4E), all these transcripts accumulated less in Dmpk4a compared with the wild type (Supplemental Figures 14A to 14E). These data implicate P. patens MPK4a in the transduction of signals from chitosan perception to the activation of defense genes.
In Arabidopsis, the contribution of PTI to resistance is often studied using pathogens that are unable to overcome PTI (Heese et al., 2007). While such host-pathovar systems have yet to be developed for P. patens, the broad host range necrotrophic fungus B. cinerea proliferates on P. patens (Ponce de León et al., 2007León et al., , 2012 and delivers effectors and secretes metabolites aimed at circumventing PTI in vascular plant hosts (Frías et al., 2011;Weiberg et al., 2013). We therefore examined the susceptibility of the Dmpk4a mutant to B. cinerea. After 2 d, symptoms of susceptibility characterized by browning protonema and stems became visible in both the wild type and Dmpk4a (Supplemental Figures 15A to 15H). We also stained tissues for the production of ROS and phenolic compounds that accumulate during B. cinerea infection (Ponce de León et al., 2007León et al., , 2012. While ROS production was detected in infected tissues, clear differences could not be seen between the wild type and Dmpk4a (Supplemental Figures  15K to 15N). Although phenolic compounds were incorporated in (E) Quantitative RT-PCR of PAL4 and CHS transcript levels in the wild type and Dmpk4a-5 as fold change relative to the wild type untreated at times specified after 100 mg/mL chitin treatment. Data points represent an average of three independent experiments, each with three technical replicates. Error bars represent SE of the means. moss cell walls in contact with B. cinerea hyphae, especially around the point of hyphal penetration, no clear differences were noted between the wild type and Dmpk4a (Supplemental Figures  15I and 15J).
In contrast to such qualitative analyses, host cell death is a hallmark of B. cinerea infection and can be quantified by Evans blue staining (Ponce De León et al., 2012). Using this assay, significantly more staining was quantified 2 d after inoculation in two independent Dmpk4a lines than in the wild type ( Figure 6A). To substantiate these observations, we assessed susceptibility to another necrotrophic fungus, Alternaria brassicicola. In this assay, 14-d-old wild-type and Dmpk4a lines were inoculated with a suspension of A. brassicicola spores and the moss plants harvested after 4 d for spore counting (van Wees et al., 2003). A. brassicicola produced significantly more spores on the independent Dmpk4a lines compared with the wild type ( Figure 6B). Taken together, these results indicate that although the fungal pathogens B. cinerea and A. brassicicola establish infections on wild-type moss, the absence of MPK4a results in further loss of basal resistance.

DISCUSSION
Many studies have shown that MAP kinase cascades contribute to PTI in vascular plants (Rasmussen et al., 2012). Much less is known about PTI in nonvascular plants including the bryophyte moss model P. patens. We provide a primary example of how the activation of a specific MPK (Pp MPK4a) in response to the fungal PAMP chitin is reduced in targeted gene knockout lines of a chitin receptor or an upstream MEKK or MPKK. In addition, KO of moss MPK4a results in reduced accumulation of defense-related transcripts and increased susceptibility to pathogens. MPK4a is not activated by ABA or osmotic stress and therefore appears to function primarily in PTI. Apart from their disease susceptibility phenotype, Dmpk4a KO lines appear indistinguishable from wildtype moss. Loss of Pp MPK4a activity is therefore not associated with the kinds of pleiotropic effects detected in MPK KO lines seen in Arabidopsis.
More specifically, our experiments showed that two P. patens MPKs are rapidly activated within 1 min in response to PAMPs, including chitin, chitin derivatives, bacterial peptidyl glycan, and bacterial culture filtrates. These moss MPKs were not activated in response to the bacterial PAMP elf18 peptide, flg22 peptides from P. syringae or Salmonella, or full-length flagellin proteins from P. syringae or Borrelia. This is in line with the lack of close homologs of the Arabidopsis EF2 and FLS2 receptors in P. patens (Boller and Felix, 2009).
The rapid and robust MPK activation by chitin and derivatives prompted us to try to identify components of this signaling pathway using targeted gene KOs and KIs by homologous recombination. Of the two moss homologs of the Arabidopsis chitin receptor At CERK1, one KO mutant (Dcerk1) failed to activate the MPKs in response to chitin. In addition, expression of Pp CERK1 restored chitin sensitivity in the Arabidopsis cerk1 mutant, indicating that Pp CERK1 is an ortholog of At CERK1. Furthermore, Dcerk1 loss of function in the moss reduced ion leakage and the induced accumulation of defense related mRNAs compared with the wild type. Loss of chitin sensitivity in Dcerk1 lines versus the wild type was also easily quantified by staining of brownish-red cell wall-associated depositions with Toluidine blue. Time-lapse photography revealed that chitin treatment of wild-type, but not Dcerk1 lines, led to growth inhibition within minutes. Molecular (A) In-gel kinase activity assay with MBP or synthetic peptide P3 as substrates, and immunoblot analysis of MPK phosphorylation with a-pTEpY of extracts from wild type treated with water, 100 mg/mL chitin (Chi), 500 mM NaCl, or 800 mM mannitol for indicated times. Loading control is Ponceau staining. (B) In-gel kinase assay with MBP and immunoblotting of MPK phosphorylation with a-pTEpY of extracts of the wild type or the ABA-insensitive line (D2-1) treated with 100 mg/mL chitin (Chi), 500 mM NaCl, or 10 mM ABA. Loading control is Ponceau staining of total protein.
(C) In-gel kinase assay with MBP of extracts of the wild type and Dsnrk2a-1 untreated (N.T.) or treated with water, 500 mM NaCl, or 10 mM ABA for 10 min.
(D) Quantitation of kinase activities in extracts relative to the wild type of two Dsnrk2a lines treated with 500 mM NaCl for 10 min. Data points represent an average of three biological replicates with SD as error bars. components involved in PAMP-triggered growth inhibition have been identified in Arabidopsis (Li et al., 2009), and studies of their P. patens homologs may provide tools to study tip growth in the moss (Vidali and Bezanilla, 2012). In addition, chitin treatment of wild-type moss, but not Dcerk1, led to a strong fluorescence burst by 4 h, which decayed by 8 h. The identities of these fluorescent molecules and the brownish-red depositions remain to be determined. We note that the chitin-induced depositions resemble those seen around hyphal penetrations during B. cinerea infection.
We used the same targeted KO approach to identify components of kinase cascades functioning downstream of CERK1 in moss. This revealed that KO of the MAP kinase kinase kinase MEKK1a had markedly reduced chitin sensitivity, while the chitin sensitivity of the KO of MEKK1b was slightly reduced. However, as chitin responsiveness was reduced in a double KO Dmekk1a/b line to the level of the chitin receptor KO Dcerk1, MEKK1a and b may have partially redundant functions in responses to chitin. However, it is also possible that the absence of MEKK1a in Dmekk1a lines leads to ectopic recruitment of MEKK1b to the pathway. Such effects have been documented in other systems (Madhani and Fink, 1998). Interestingly, chitin induced the same rapid growth inhibition in Dmekk1a/b as in the wild type. Since chitin failed to inhibit the growth of Dcerk1, this inhibitory response to chitin is independent of the kinase cascade in which MEKK1a and b function.
Targeted KOs of three MKK genes revealed that KO of MKK1a and c caused a modest reduction in chitin sensitivity. This indicates that MKK1a and c may be partially redundant with other MKK genes in responses to chitin, with the above caveat of potential ectopic recruitment of other MKK s in the Dmkk1a KO line. However, in Arabidopsis, At MPKK1 and 2 act redundantly in signal transduction downstream of PRR activation; therefore, redundancy at this level of the MPK signaling cascade is expected .
A recent phylogenic analysis could not establish clear orthologous relationships between MPKs from Arabidopsis and basal land plants including P. patens (Dóczi et al., 2012). Since we and others have shown that the Arabidopsis group B MAP kinase MPK4 is involved in immunity, we created KO and KI lines of the two P. patens group B MPKs MPK4a and MPKb. Analysis of KI lines expressing Pp MPK4a or b GFP fusions demonstrated that these two kinases are the two enzymes activated in response to chitin and other PAMPs. While we were unsuccessful in generating Dmpk4b KO lines, Dmpk4a lines were strongly compromised in (1) accumulation of chitin-induced cell wall-associated depositions, (2) expression of defense-related genes, and, most significantly, (3) basal resistance against necrotrophic fungi.
Interestingly, we did not detect phosphorylation or activation of P. patens MPKs in response to osmotic stress treatments or ABA. This is in contrast to recent reports of MPK functions in ABA and osmotic stress responses in Arabidopsis (Montillet et al., 2013;Danquah et al., 2015). We also found no difference between the wild type and Dmpk4a in the accumulation of abiotic stress responsive transcripts upon NaCl treatment or in the survival of the wild type and Dmpk4a on high salt media (S. Bressendorff, unpublished data). These results indicate that Pp MPK4a primarily functions in PTI signaling and thus differs from its close homologs in vascular plants, which are also activated by abiotic stress (Ichimura et al., 2000;Droillard et al., 2004;Teige et al., 2004;González Besteiro et al., 2011). We conclude that MPK signaling in response to the abiotic stresses examined here could either have been lost in the P. patens lineage or have evolved subsequently in vascular lineages.
While ABA and osmotic stress treatments did not induce the activity of kinases of similar size of moss MPKs using MBP as substrate ($42 kD), they did induce the activity of one or more kinases of ;40 kD within 5 min. This activity was presumably not that of an MPK, as it was not detected with a-pTEpY antibodies. We showed that this activity corresponds to that of Pp SnRK2a, and likely of other redundant SnRK2s, as seen in Arabidopsis (Fujii et al., 2007). This activity was strongly reduced in an ABA insensitive moss line expressing a constitutively active PP2C phosphatase (Tougane et al., 2010). Such ectopic PP2C activity deactivates SnRK2s (Yoshida et al., 2006;Boudsocq et al., 2007;Vlad et al., 2010). This is in accordance with the findings of Komatsu et al. (2013) on the activation of one or more SnRK2s of <40 kD by ABA in P. patens and of D'Souza and Johri (2002) on an ;38-kD kinase activated by NaCl and ABA in the closely related moss Funaria hygrometrica. We showed that this kinase activity phosphorylates the SnRK-specific substrate peptide P3 (Halfter et al., 2000) using a novel in-gel kinase assay with acrylamide cross-linked acryloylated P3. In addition, this osmotic stress and ABA-inducible SnRK2 activity was strongly reduced in the targeted knockout lines of Dsnrk2a. These results indicate that SnRK2a, but not MPKs, functions in a canonical ABA signaling pathway in P. patens. Saruhashi et al. (2015) recently reported similar findings of SnRK2 and not MPK activity induced in (A) Evans blue staining 2 d after spraying with water (white bars) or 2 3 10 5 B. cinerea spores per milliliter (gray). Staining measured as OD 600 per mg dry weight (DW). Treated samples are presented as an average of four independent biological replicates with SD as error bars. Analysis of variance followed by Tukey's test determined statistical differences indicated by two asterisks (P < 0.01) between the wild type and the two independent knockout lines (Dmpk4a-1 and 2). The experiment was repeated twice with similar results. (B) A. brassicicola spore counts 4 d after inoculation with ;2500 spores per P. patens plant. Each data point is an average of four pools of 16 plants per genotype. Analysis of variance followed by Tukey's test determined statistical differences indicated by three asterisks (P < 0.001) between the wild type and the two knockout lines. The experiment was repeated twice with similar results. P. patens in response to ABA and hyperosmotic stress. They further found that this SnRK2 activity was dependent on a Raf-like kinase.
The apparent phosphorylation of moss MPK4b in response to chitin and PGN indicates that it is also involved in PAMP signaling, perhaps redundantly with MPK4a. While our inability to produce Dmpk4b lines currently obscures such a function, it suggests that MPK4b is also required for metabolism and/or development in the moss. Future studies using kinase chemical genetics in P. patens could elucidate MPK4b functions (Brodersen et al., 2006).
Pp MPK4a is a close homolog of At MPK4 but, in apparent contrast to the dwarfed Arabidopsis mpk4 mutant, the moss Dmpk4a mutant appears phenotypically normal under all conditions tested. While this suggests that Pp MPK4a and At MPK4 have different functions, both are responsive to pathogens and PAMPs and At MPK4 is implicated in the regulation of PTI (Qiu et al., 2008). One explanation for their different mutant phenotypes is that At MPK4 and orthologs in vascular plants have evolved additional functions. A second explanation relates to the evolution of host-pathogen interactions. Successful pathogens deliver effectors to modify host proteins and suppress immune responses (Bent and Mackey, 2007). Plants have developed mechanisms to detect pathogen effectors via cytoplasmic immune receptors (termed resistance or R proteins) that typically trigger a rapid, localized host cell death known as the hypersensitive response (DeYoung and Innes, 2006). R proteins may thus guard effector host targets or guardees (Jones and Dangl, 2006). Zhang et al. (2012) recently reported that the absence of At MPK4 in the mpk4 mutant triggers immunity via an R protein. Thus, Arabidopsis MPK4 or its pathway may be targeted by pathogen effectors to block PTI and thereby prevent defense gene activation. If so, it is possible that moss MPK4a remains unguarded in the moss. This difference may have enabled us to assess the phenotype and function of Pp MPK4 more readily than has been the case for At MPK4. In a similar vein, our failure to isolate Dmpk4b mutants may be because Pp MPK4b performs essential functions or because it is in fact guarded. Although some 18 R genes are found in P. patens (Xue et al., 2012), nothing is currently known of their functions.

Generation of P. patens Mutants
Cloning primers are listed in Supplemental Table 1. The Dmpk4a transformation vector was constructed by cloning a PCR product of the 59 upstream region of MPK4a into the NotI and BamHI sites of P. patens transformation vector pMBL10a. The 39 downstream region of MPK4a was then cloned in the EcoRV site of this vector. The KO transformation vectors for P. patens CERK1, 1s41 280V6.1, MEKK1a, MEKK1b, MKK1a, MKK1b, MKK1c, and SnRK2a were initially constructed by adding USER cloning sites on both sides of the NptII selection cassette in vector pMBL6 (Nour-Eldin et al., 2006). The respective genes were cloned into the resulting USER cloning vector pMBLU by PCR amplification of their 59 upstream and 39 downstream regions with USER compatible primers. These two PCR products were combined in a single, four-fragment USER cloning step into pMBLU cut with PacI and AsiSI and nicked with Nt.BbvCI (Jacobsen et al., 2011). GFP and a nos terminator were cloned in the pMBLU KpnI site to make the USER cloning vector pMBLU-GFP for tagging proteins with a C-terminal GFP. P. patens MPK4a-GFP and MPK4b-GFP were cloned by PCR amplification of their genomic regions and of 39 regions downstream of the genes with USER compatible primers. To make the double KO line Dmekk1a/b, pMBLU was modified to contain an HtpII selection cassette (hygromycin resistance) by cutting pMBLU with SalI and KpnI and inserting HtpII amplified from pUNI33. P. patens MEKK1b was cloned into the resulting vector pMBLU-hyg and transformed into the single KO line Dmekk1a-1 to create Dmekk1a/b. The pro35S:GFP line was made for random integration by cloning the 35S promoter from pEarly103 into the pMBLU-GFP by cutting with PacI and nicked with Nt.BbvCI. P. patens protoplasts were transformed with 30 mg of linearized vectors using polyethylene glycol as previously described (Liu and Vidali, 2011). Stable transformants were selected by transferring the protoplasts on cellophane to media containing 50 mg/mL G418 or 30 mg/mL hygromycin B for 2 weeks, followed by transfer to media without selection for 2 weeks, and finally by transfer onto media with selection for a further 2 weeks (Schaefer et al., 1991).
Complementation of Arabidopsis cerk1 with P. patens CERK1 P. patens CERK1 was cloned into pENTR/D-TOPO, modified with InFusion Overhangs, to be fused with a Gateway destination vector. First, the cDNA corresponding to CERK1 was isolated by RT-PCR using RevertAid reverse transcriptase (Thermo Scientific) from 0.5 mg of total P. patens RNA. The open reading frame of CERK1 was obtained by PCR using Phusion High-Fidelity DNA Polymerase (Thermo Scientific) and gene-specific primers. pENTR was digested with AscI and NotI and the CERK1 cDNA cloned into the vector using In-Fusion HD Cloning Plus (Clontech). The plasmid obtained was then transferred to Gateway destination vector pGWB514 with a C-terminal HA tag using Gateway LR Clonase II Enzyme mix (Life Technologies) according with the manufacturer's specifications. Escherichia coli XL blue cells were transformed with these plasmid constructs that were then transformed into Agrobacterium tumefaciens Agl1 electro-competent cells. Arabidopsis cerk1 sil1-2 (Gabi-kat line) was transformed using Agrobacterium floral dipping as previously described (Clough and Bent, 1998).

Treatments for Kinase Assays
Unless otherwise stated, all treatments were applied for 5 min by spraying 3 mL of the solution used onto 16 moss plants grown for 14 d on a Petri dish (9 cm) with BCDAT media overlaid with cellophane. Chitin oligosaccharide (Yaizu Suisankagaku Industry) was diluted to 100 mg/mL in water. Hexa-Nacetylchitohexaose (Sigma-Aldrich) was diluted to 100 mg/mL in water. Chitosan (Sigma-Aldrich), 75 to 85% deacetylated and low molecular weight (50 to 190 kD), was solubilized at 100 mg/mL in 0.010% acetic acid followed by adjusting to pH 5.5 with NaOH. Acetic acid (0.010%) adjusted to pH 5.5 was used as a control. Full-length flagellins from Pseudomonas aeruginosa (Invitrogen) and from Borrelia (Alpha Diagnostic), and flg22 and flgII-28 purified peptides (Schafer-N) were sprayed on P. patens and applied to Arabidopsis in liquid culture at the indicated concentrations. Then, 1 3 10 7 B. cinerea spores per milliliter were sprayed onto plants. Abscisic acid (Sigma-Aldrich) was diluted to a concentration of 10 mM in 0.01% (w/w) ethanol. Peptidyl glycan, PGN-Staphylococcus aureus (Sigma-Aldrich), was applied by gentle shaking in a suspension of 200 mg/mL and Xanthomonas campestris PGN, obtained as previously described (Erbs et al., 2008), was similarly applied at 100 mg/mL. NaCl was applied at 500 mM and mannitol at 800 mM. Cell-free culture filtrates of P. wasabiae SCC3193 were made by growing bacteria in LB broth overnight, removing bacterial cells by centrifugation (10 min at 4000g), and filter sterilizing the supernatant. Detached leaves from three independent 14-d-old T1 complementing lines expressing Pp CERK1 in the Arabidopsis cerk1 mutant background were preincubated for 2 h in 10 mM MgCl 2 . Chitin oligosaccharide (Yaizu Suisankagaku Industry) was added to a final concentration of 200 mg/mL for 10 min to the treated samples.

Toluidine Blue Staining
Wild-type and mutant moss were grown for 14 d on BCDAT overlaid with cellophane and then transferred on the cellophane onto BCDAT plates supplemented with 300 mg/mL chitin. After 48 h the five plants per sample were weighed in microfuge tubes and plants stained in 0.05% Toluidine blue (Sigma-Aldrich) in citrate-citric acid buffer (50 mM, pH 3.5) with 50% ethanol for 10 min. The plants were quickly washed twice in water and then washed twice in 96% ethanol and left in 96% ethanol overnight. Representative plants were then photographed. The Toluidine blue staining was quantified by destaining in 10% acetic acid, 86% ethanol at 60°C for 4 h, and the OD of the supernatant was measured at 600 nm as previously described (Beasley et al., 1974) and data presented as OD 600 /mg fresh weight. Treated samples represent an average of three samples of five plants each with SD as error bars. Nontreated controls were treated equally except plants were not exposed to chitin.

Electrolyte Leakage Assay
Wild-type and Dcerk1-1 moss were grown for 14 d on BCDAT overlaid with cellophane. Five plants were then briefly washed in milli-Q water and transferred to culture tubes containing either 5 mL water or 5 mL 1 mg/mL chitin. The conductivity was measured at indicated times, and at the end of the time series (34 h) total conductivity of each sample was established by freezing the samples in liquid nitrogen, thawing them, and measuring conductivity. Electrolyte leakage is presented as the percentage of total conductivity for each sample as an average of four samples with each five plants. Error bars represent SD.

Evans Blue Staining for Detection of Cell Death
Cell death quantitation by Evans blue staining was similar to the assay of Oliver et al. (2009). Four plates of each line, each with 16 plants, were sprayed with 3 mL of 2 3 10 5 B. cinerea spores per milliliter and one plate sprayed with water. Two days later, plants were incubated in 0.1% Evans blue (Bie and Berntsen) in 0.53 PBS for 2 h, washed four times in water, then destained in 50% (v/v) methanol, 1% SDS at 60°C for 30 min, and the OD 600 measured. Plants were dried over night at 70°C and weighed. Each data point is presented as OD 600 /mg dry weight as an average of four samples of 16 Plants each. Error bars represent SD. An analysis of variance followed by a Tukey's test was applied to assess significant differences. The experiment was repeated twice with similar results.

A. brassicicola Spore Count Assay
The assay was adapted from van Wees et al. (2003). A. brassicicola spores were harvested in water from 14-d-old plates and counted in a hemocytometer. Moss plants were inoculated by placing 5 mL of 2 3 10 5 A. brassicicola spores/mL on top of each plant. Four plates of each line, each with 16 plants, were inoculated, and 4 d later the plants were transferred to a tube containing 8 mL 0.01% (v/v) Tween 20 (Sigma-Aldrich) and shaken vigorously. Spores in the suspension were counted in a hemocytometer. Each data point is an average of four pools of 16 plants per genotype. Analysis of variance followed by a Tukey's test assessed differences between the wild-type and two Dmpk4a lines. The experiment was repeated twice with similar results.

RNA Extraction and RT-qPCR
Three 14-d-old moss plants grown on media overlaid with cellophane were sampled at each time point after spraying with 500 mM NaCl, 100 mg/mL chitosan (in 0.010% acetic acid followed by adjusting to pH 5.5 with NaOH), or 100 mg/mL chitin and flash-frozen in liquid nitrogen. Arabidopsis plants were grown for 14 d on soil, and detached leaves were preincubated for 2 h in 10 mM MgCl 2 before chitin was added to a final concentration of 200 mg/mL for 30 and 60 min and the samples flash-frozen in liquid nitrogen. RNA was extracted with NucleoSpin RNA plant (Macherey-Nagel) according to the manufacturer's protocol. RNA concentrations were measured in a Nano Drop 1000 (Thermo Scientific) and adjusted to 5 ng/mL. Quantitative PCR used Brilliant II SYBR green one-step kit (Agilent Technologies), with 10 pmol of each primer and 12.5 ng total RNA in 10 mL. Reactions were run on a CFX 96 thermocycler (Bio-Rad) with three technical replicates for each PCR reaction, and means of normalized expression relative to untreated wild type calculated with b-TUBULIN as internal normalization control for the P. patens samples (Zhang et al., 2015) and UBQ10 for the Arabidopsis samples. Primers used for qPCR are in Supplemental Table 1.

In-Gel Kinase Assay
Three 14-d-old plants grown on media overlaid with cellophane were flashfrozen at given times after treatment. Protein was extracted in Lacus buffer. Forty micrograms of protein was resolved on 13% SDS-PAGE with 14 mM MBP (Millipore) or 14 mM P3 (N-terminal acryolated by 6-amino hexanoic acid peptide sequence, ALARAASAAALARRRALARAASAAALARRR; Schafer-N) polymerized in separating gels as substrate (Zhang and Klessig, 1997). Gels were washed six times for 3 h in 25 mM Tris-HCl (pH 7.5), 0.5 mM DTT, 5 mM NaF, 0.1 mM Na 3 VO 4 , 0.5 mg/mL BSA, and 0.1% Triton X-100. Proteins were denatured by rinsing gels twice for 30 min at room temperature in 25 mM Tris-HCl (pH 7.5), 6 M guanidine hydrochloride, and 5 mM DTT, followed by incubation overnight with five changes of 25 mM Tris-HCl (pH 7.5), 1 mM DTT, 5 mM NaF, and 0.1 mM Na 3 VO 4 . After gel equilibration for 1 h in kinase activity buffer (25 mM Tris-HCl, pH 7.5, 2 mM EGTA, 12 mM MgCl 2 , 1 mM DTT, and 0.1 mM Na 3 VO 4 ), the assay was performed for 1 h at room temperature in that buffer with 100 nM ATP plus 100 mCi [g-32 P]ATP. The reaction was stopped with 5% TCA and 1% sodium pyrophosphate, and unincorporated [g-32 P]ATP was removed by washing five times for 15 min in the same buffer. Gels were dried and phosphorylation determined by phosphor imaging.

Time-Lapse Movies of Growth Inhibition and Fluorescent Burst
Wild-type and Dcerk1-1 or wild-type and Dmekk1a/b moss plants were grown in close proximity on BCD media overlaid with cellophane for 10 d and then transferred on cellophane onto BCD plates supplemented with 300 mg/mL chitin. Pictures were taken every 2nd or 4th minute using a Sony A7S camera. Growth inhibition (Supplemental Movies 1 and 4) was recorded in white light, while the fluorescent burst (Supplemental Movies 2 and 3) was recorded using GFP settings (excitation <488 nm and detection >510 nm).

Confocal Microscopy
Confocal images were acquired using a Zeiss LSM 700 and analyzed with ZEN imaging software (version 2011) with GFP settings (excitation 488 nm and detection 300 to 610 nm). The 14-d-old MPK4a-GFP plants were mounted on a glass slide embedded in BCDAT with 1% agar and supplemented with 300 mg/mL chitin.

Staining of B. cinerea Infection Symptoms
The 14-d-old moss plants of wild type and Dmpk4a-1 were sprayed with 2 3 10 5 B. cinerea spores per milliliter. Toluidine blue staining of cell wall modifications was conducted 1 d after infection by staining tissue with 0.05% Toluidine blue in citrate-citric acid buffer (50 mM, pH 3.5). Cell death was visualized by staining for 2 h with 0.1% Evans blue (Sigma-Aldrich) in 0.53 PBS and washed four times in water to remove excess unbound dye. Intracellular ROS production was analyzed by incubating moss tissues with 10 mM 29,79-dichlorodihydrofluorescein diacetate for 15 min in 0.1 M phosphate buffer (pH 7.5) in the dark.

Accession Numbers
A list of accession numbers is provided in Supplemental Table 1.

Supplemental Data
Supplemental Figure 1. MAP kinase activities in response to PAMP treatments.
Supplemental Figure 2. Phylogeny of CERK1 homologs and mRNA levels of P. patens CERK1 homologs.