Adaptation Mechanisms in the Evolution of Moss Defenses to Microbes

Bryophytes, including mosses, liverworts and hornworts are early land plants that have evolved key adaptation mechanisms to cope with abiotic stresses and microorganisms. Microbial symbioses facilitated plant colonization of land by enhancing nutrient uptake leading to improved plant growth and fitness. In addition, early land plants acquired novel defense mechanisms to protect plant tissues from pre-existing microbial pathogens. Due to its evolutionary stage linking unicellular green algae to vascular plants, the non-vascular moss Physcomitrella patens is an interesting organism to explore the adaptation mechanisms developed in the evolution of plant defenses to microbes. Cellular and biochemical approaches, gene expression profiles, and functional analysis of genes by targeted gene disruption have revealed that several defense mechanisms against microbial pathogens are conserved between mosses and flowering plants. P. patens perceives pathogen associated molecular patterns by plasma membrane receptor(s) and transduces the signal through a MAP kinase (MAPK) cascade leading to the activation of cell wall associated defenses and expression of genes that encode proteins with different roles in plant resistance. After pathogen assault, P. patens also activates the production of ROS, induces a HR-like reaction and increases levels of some hormones. Furthermore, alternative metabolic pathways are present in P. patens leading to the production of a distinct metabolic scenario than flowering plants that could contribute to defense. P. patens has acquired genes by horizontal transfer from prokaryotes and fungi, and some of them could represent adaptive benefits for resistance to biotic stress. In this review, the current knowledge related to the evolution of plant defense responses against pathogens will be discussed, focusing on the latest advances made in the model plant P. patens.

associated to present-day liverworts and mosses include bacteria adapted to aquatic, anaerobic and extreme drought environments, which is consistent with the transition of bryophytes from aquatic to terrestrial conditions (Tang et al., 2016). Due to their key evolutionary position, between unicellular green algae and flowering plants (angiosperms), bryophytes represent interesting organisms to explore the evolution of plant-microbe interactions across the green plant lineage. Recent advances in genome-sequencing projects, including basal angiosperms, gymnosperms, lycophytes, bryophytes, and charophyta algae, have facilitated the identification of gene acquisitions related to terrestrialization, including adaptation mechanisms to interacting microbes. Like other plants, bryophytes have associated microbial communities with beneficial activities. Microbial symbioses appear as one of the possible key innovations in land colonization by plants (Read et al., 2000). The symbiosis with arbuscular mycorrhizal fungi (AMF) and other fungal beneficial associations probably facilitated land colonization by plants improving the capture of nutrients (Delaux et al., 2015). Consistently, plant fossil evidence of fungal structures associated with early bryophytes in 460 and 400 million years old Ordovician sediments has been found (Redecker et al., 2000;Strullu-Derrien et al., 2014). AMF form mutualistic symbiosis with liverwort, improving plant growth and fitness through enhancing nutrient uptake (Humphreys et al., 2010). Liverworts associated with AMF exhibited significantly enhanced photosynthetic carbon gain, increased gametophytic growth and plant fitness, evidenced by greater production of asexual reproductive structures. The benefits of AMF in high CO 2 atmosphere of the Palaeozoic are amplified by increasing reproduction and growth-limiting rates of phosphorus uptake in a rootless basal land plant (Humphreys et al., 2010). Interestingly, symbiotic signaling genes were already present in algae, and in early land plants diversification of the transcription factors GRAS and half-ABC transporters took place via gene duplication events together with the acquisition of other downstream symbiotic genes (Delaux et al., 2015). The optimization and specificity of symbioses took place in liverworts and angiosperms by independent recruitments of symbiotic nutrient transport systems (Delaux et al., 2015). Rescue assays on flowering plant mutants with orthologs from bryophytes revealed that some plant components required for the establishment of symbiosis are conserved during land plant evolution (Delaux et al., 2013). Cyanobacterial associations with mosses, liverworts, and hornworts have also been reported (DeLuca et al., 2002;Meeks and Elhai, 2002). Interestingly, motile cyanobacterial (Nostoc sp.) filaments serving as infection units (hormogonia) require unknown moss signal(s) to move toward and colonize plant tissues (Bay et al., 2013). In the interactions with mosses, cyanobacterial colonization contributes to host nutrition by providing the plant with fixed nitrogen (Bay et al., 2013).
The moss Physcomitrella patens (P. patens) produces chemical compounds that increase swarming of Rhizobium leguminosarum, which is a process of rapid and coordinated movement of bacteria leading to rhizosphere colonization (Tambalo et al., 2014). Movement is probably related to increased expression of genes involved in flagellar filament formation and motility. Bacterial swarming behavior can be influenced by plant strigolactone, since extracts of a strigolactone-deficient moss mutant induce less swarming motility in R. leguminosarum. Accordingly, the existence of beneficial interactions between moss and rhizobia, suggest a potential ancient signal perception mechanism between these type of bacteria and plants (Tambalo et al., 2014). Studies of moss-associated bacteria in different ecosystems have shown the presence of diverse microbial populations with a high proportion of antagonistic isolates to fungal pathogens, suggesting that mosses harbor bacteria with antifungal activities that contribute to plant defense (Opelt et al., 2007a,b). Further research is needed to understand the adaptation mechanisms developed by bryophytes and beneficial microorganisms, including chemical communication and the biochemical function of the proteins involved in the establishment of the beneficial interaction.
Bryophytes are also colonized by microbial pathogens, including fungi, bacteria, oomycetes and viruses that cause disease evidenced by necrosis and tissue maceration (Ponce de León, 2011). Among bryophytes, P. patens has become an attractive model plant for studying plant-pathogen interactions and several moss defense strategies have been described (Ponce de León and Montesano, 2013). The sequenced genome and available genomic resources 1 (Rensing et al., 2008;Zimmer et al., 2013), allow the identification of genes and metabolic pathways involved in adaptation mechanisms adopted by mosses to defend themselves against microbial pathogens. In addition, targeted knockouts of genes with possible roles in defense can be generated in P. patens due to high rates of homologous recombination (Schaefer, 2002). The presence of one layer of cells in most tissues facilitates microscopic analysis of invasion and colonization processes as well as moss defense mechanisms (Ponce de León, 2011;Overdijk et al., 2016; Figure 1). This review summarizes the current knowledge on the adaptation mechanisms and the evolution of plant defenses against pathogens, focusing on the latest advances made in the model plant P. patens.

PATHOGEN PERCEPTION AND KINASE ACTIVATION IN P. PATENS
Physcomitrella patens is infected by pathogens such as fungi, bacteria and oomycetes leading to tissue damage and disease (Lehtonen et al., 2009;Ponce de León and Montesano, 2013;Overdijk et al., 2016). Several of these pathogens infect mosses in nature, and some of them cause disease in important crops (Ponce de León, 2011). Similarly to angiosperms, P. patens perceives the presence of pathogens and activates a defense response evidenced by increased production of reactive oxygen species (ROS), induction of a hypersensitive-like response (HR), cell wall reinforcement, changes in hormones levels and expression of defense genes (Ponce de León et al., 2007León et al., , 2012 Figures 1, 2). Angiosperms have pattern recognition FIGURE 1 | Pathogen invasion and P. patens defense responses. P. patens has one layer of cells in most tissues facilitating microscopic analysis of pathogen colonization processes and moss defense mechanisms (Ponce de León, 2011;Overdijk et al., 2016). Leaves and rhizoids are colonized by different pathogens, which are easily visualized by different staining techniques (leaf; trypan blue and solophenyl flavine, rhizoid: solophenyl flavine) (Ponce de León and Montesano, 2013). Representative pictures of the oomycete P. irregulare and the fungus B. cinerea pathogen-infected tissues are shown. Plant defenses are activated after pathogen infection, including a HR-like response, cell wall fortification and chloroplast reorientation (red arrow), ROS accumulation, defense gene expression and increase in hormone levels (Ponce de León and Montesano, 2013). Representative pictures of defense responses are shown; HR-response, cell wall fortification (phenolic compounds accumulation, evidenced by stained cell walls) and ROS production in B. cinerea-infected leaves, callose deposition in P.c. carotovorum elictor treated protonemal cells (white arrow), and expression pattern of the defense gene Ppalpha-DOX and auxin levels in response to elicitors of P.c. carotovorum at different hours after treatment. The HR-like response was visualized by autofluorescent compounds, cell wall fortification by safranin-O staining, callose deposition by methyl blue staining, and intracellular ROS with 2 ,7 -dichlorodihydrofluorescein diacetate. "Graph showing auxin levels as originally published in Alvarez et al. (2016)." receptors (PRRs) at the plasma membrane to detect conserved pathogen-associated molecular patterns (PAMPs), and induce a defense response called PAMP-triggered immunity (PTI), providing protection against non-host pathogens Zipfel, 2009). Pathogens adapted to their host plants can deliver effector molecules into plant cells, which target key PTI components and inhibit plant defense (Boller and Felix, 2009;Boller and He, 2009). In turn, plants have a FIGURE 2 | Activation of defense responses in Arabidopsis and P. patens. (A) Angiosperms like Arabidopsis thaliana sense the presence of pathogen-associated molecular patterns (PAMPs) by plasma membrane (PM) pattern recognition receptors (PRRs), such as FLS2, EFR1, CERK1, and LYK5 (Couto and Zipfel, 2016). Pathogen recognition triggers calcium (Ca 2+ ) production and activates a MAPK cascade (Seybold et al., 2014). MAP kinase kinase kinases (MEKKs), MAP kinase kinase (MKKs), and MAP kinase (MPKs) are subsequently activated, leading to production of reactive oxygen species (ROS), and expression of defense genes (Meng and Zhang, 2013). The hormones salicylic acid (SA), jasmonic acid (JA) and its precursor cis-oxophytodienoic acid (OPDA), abscisic acid (ABA), auxins, and ethylene (ET) activate the expression of defense genes leading to the activation of defense mechanisms that involves pathogenesis-related (PR) proteins, incorporation of phenolics into the cell wall (CW), deposition of callose and lignin, activation of an hypersensitive response (HR) and systemic acquired resistance (SAR) (Glazebrook, 2005). In response to pathogen infection, oxylipins derived mainly from linoleic acid (18:2) and linolenic acid (18:3; C18 fatty acids) are synthesized . (B) P. patens lacks close homologs of the receptors FLS2 and EFR, while a functional CERK1 receptor perceives fungal chitin and bacterial peptidyl glycan (Bressendorff et al., 2016). At least one MEKK, one MKK, and two MPKs (MPK4a and MPK4b) participate in P. patens defense responses to fungal chitin (Bressendorff et al., 2016). ROS, SA, and auxin activate expression of defense genes, while JA is not produced (Ponce de León and Montesano, 2013;Reboledo et al., 2015). Further studies are needed to reveal the role of OPDA, ABA and ethylene (names in gray) in moss resistance against pathogens. After pathogen assault, P. patens activates several defense mechanisms, including the expression of genes encoding PR proteins, incorporation of phenolics into the CW, callose deposition and accumulation of pre-lignin compounds (Ponce de León and Montesano, 2013;Overdijk et al., 2016). An HR-like response and SAR are also activated in infected mosses, and oxylipins derived from C18 and polyunsaturated C20 fatty acids are synthesized producing a broader range of oxylipins with possible roles in plant defense (Ponce de León et al., 2012, 2015Winter et al., 2014). Orphan genes, some of them acquired by gene horizontal transfer from fungi and other microorganisms, could also play a role in moss defense against pathogens. Names in red or blue indicate their presence only in traqueophytes or P. patens, respectively. second layer of immune receptors encoded by resistance (R) genes to detect directly or indirectly the effector proteins leading to effector-triggered immunity (ETI), which is highly specific and often accompanied by a HR and systemic acquired resistance (SAR). Up to date, several PAMPs and the corresponding PRRs have been identified in angiosperms. In Arabidopsis thaliana (A. thaliana), the receptors FLS2 (for flagellin-sensing 2), EFR1 (for Elongation factor tu receptor 1), and LYK1/CERK1 (for LysM-containing receptor-like kinase1/chitin elicitor receptor kinase1), recognize bacterial flagellin, the elongation factor Tu, and the fungal chitin, respectively (Gómez-Gómez and Boller, 2000;Zipfel et al., 2006;Miya et al., 2007). LYK4 and LYK5 are also important for chitin signaling and immunity in A. thaliana (Wan et al., 2012;Cao et al., 2014). LYK5 was identified as the receptor that recognizes chitin in A. thaliana, forming a chitin inducible complex with CERK1 to induce plant immunity (Cao et al., 2014). Interestingly, P. patens lacks close homologs of the receptors FLS2 and EFR (Boller and Felix, 2009), which is in accordance to the insensibility of moss cells to flagellin (flg22) and the elongation factor Tu (Bressendorff et al., 2016). A functional CERK1 receptor was recently identified in P. patens and PAMPs such as fungal chitin and bacterial peptidyl glycan are perceived by this receptor (Bressendorff et al., 2016; Figure 2). Mutation in moss PpCERK1 leads to reduced defenses, including less expression of defense genes and less incorporation of phenolic compounds associated to cell wall defenses, showing that PTI is an ancient plant defense response (Bressendorff et al., 2016). Further studies are needed to understand chitin perception by P. patens, including analyses of molecular complexes related to PpCERK1 and putative LYK5-like receptors leading to the activation of moss immunity.
Other receptor-like kinases (RLKs) such as Catharanthus roseus RLK1-like (CrRLK1L) proteins are also linked to plantpathogen interactions, and orthologs have been found in P. patens, the liverwort Marchantia polymorpha (M. polymorpha) and the primitive vascular plant Selaginella moellendorffii (S. moellendorffii) (Galindo-Trigo et al., 2016). In A. thaliana two CrRLK1Ls have been directly associated with defense responses to biotic stress, and expression data in elicitortreated and pathogen-infected tissues suggest that additional CrRLK1Ls are involved (Lindner et al., 2012;Nissen et al., 2016). In A. thaliana, one of these CrRLK1Ls (Feronia) is phosphorylated upon treatment with flg22, and it is speculated that it may act as a co-receptor of FLS2 (Keinath et al., 2010). CrRLK1Ls could have a conserved role throughout the plant lineage recognizing cell wall changes caused by invading fungal hyphae, and relocating proteins involved in resistance toward the penetration site, such as MLO (Mildew Resistance Locus O) (Lindner et al., 2012). Interestingly, a closely related CpRLK1 protein was identified in a charophycean unicellular alga and functional characterization suggests that this protein could also be involved in cell wall stability sensing, specifically by detecting cell wall integrity during sexual reproduction (Galindo-Trigo et al., 2016). Other plant RLKs related to defenses to biotic stress are the RLK/Pelle genes, which expanded significantly with the establishment of land plants, suggesting an adaptation mechanism to fastevolving pathogens (Lehti-Shiu et al., 2009). Consistently, while P. patens and traqueohpytes have large RLK/Pelle gene families, only two predicted members with no extracellular domain were identified in the alga Chlamydomonas reinhardtii (Lehti-Shiu et al., 2009). P. patens has also several putative RLKs with high similarities to the rice receptor OsXa21 (Oryza sativa Xanthomonas resistance 21) as well as homologous proteins to the corresponding interacting proteins (PpXB3s) (Tanigaki et al., 2014). Some of these PpXB3s interacts with several PpRLKs in a yeast two hybrid system, suggesting that P. patens could have an OsXa21-like RLK-type sensing system to detect molecules from pathogens (Tanigaki et al., 2014).
Several putative intracellular receptor genes (R genes) have been identified in P. patens, including kinase-Nucleotide Binding Site (NBS)-Leucine Rich Repeat (LRR) receptors and Toll/interleukin-1 like Receptor (TIR)-NBS-LRR (Akita and Valkonen, 2002;Xue et al., 2012;Tanigaki et al., 2014). Since algae do not have NBS-LRR, TIR-NBS-LRR, or TIR-LRR homologs (Sarris et al., 2016), R-like genes appeared in basal land plants probably as an adaptation mechanism to detect pathogens and activate a defense response. The NBS-LRR containing protein family has expanded in flowering plant species and this increase is associated with polyploidy or ancient polyploidization events. Surprisingly, the repertoire of NBS-LRR and other putative R-genes is much smaller in the primitive vascular plant S. moellendorffii compared to P. patens (Sarris et al., 2016). Further studies are needed to understand if pathogen effectors could suppress moss defense responses, and whether they could be detected directly or indirectly by some of these R-like proteins activating an ETI in P. patens.
Calcium (Ca 2+ ) is a second messenger essential in both angiosperm PTI and ETI responses leading to ROS and nitric oxide production, and induction of plant defense genes (Zhang et al., 2014). An increase of cytosolic Ca 2+ is generated after pathogen perception, and Ca 2+ -dependent protein kinase (CDPK) and MAPK pathways modulate plant biotic stress signaling (Seybold et al., 2014). To our knowledge the role of Ca 2+ and CDPKs in moss-pathogen interactions has not been studied so far. MAPK cascades are involved in signaling multiple defense responses, including the biosynthesis and signaling of defense hormones, ROS generation, defense gene activation, phytoalexin biosynthesis, cell wall reinforcement, and HR cell death (Meng and Zhang, 2013). MAPK cascade typically consists of a MAP kinase kinase kinase (MEKK), a MAP kinase kinase (MKK), and a MPK (Rodriguez et al., 2010). In A. thaliana at least four MPKs (MPK3/4/6 and 11) are activated by PRRs (Fiil et al., 2009;Bethke et al., 2012). Activation of MPK4 requires MEKK1 and MKK1/MKK2 (Gao et al., 2008). It was recently discovered that in P. patens at least two MPKs are rapidly phosphorylated and activated in response to bacterial and fungal PAMPs (Bressendorff et al., 2016). P. patens MKK1b, MKK1c, MPK4a, and MPK4b transcript levels are rapidly induced after chitin treatment. By generating targeted knockout mutants in several kinases it was shown that at least one MEKK, one MKK, and two MPKs (MPK4a and MPK4b) are components of an immune pathway in P. patens, and are required for defense responses to the fungal chitin (Bressendorff et al., 2016). In addition, MPK4a mutants result in loss of basal resistance to pathogens, evidenced by reduced cell wall associated defenses and reduced expression of defense genes, showing that a basal defense signaling pathway is already present in P. patens. In contrast to A. thaliana, abiotic stress in P. patens does not activate MPK4a or any other MPK, suggesting that MPK4a is specific to immunity in basal land plants like mosses (Bressendorff et al., 2016).

REACTIVE OXYGEN SPECIES AND HYPERSENSITIVE RESPONSE CELL DEATH IN P. PATENS DEFENSE
After pathogen recognition, one of the first plant responses involves the production of ROS, which are directly toxic to pathogens and are involved in cell wall reinforcement, HR induction and signaling leading to gene expression (Torres et al., 2006). Several enzymes are involved in ROS production during plant-pathogen interactions, including peroxidases and NADPH oxidases (Torres et al., 2002;Bindschedler et al., 2006). Like in angiosperms, several lines of evidence suggest that ROS participate in P. patens defense against pathogens. Exposure of P. patens tissues to the fungal elicitor chitosan leads to rapid production of ROS, increased expression of alternative oxidase (PpAOX), NADPH-oxidase (PpNOX), and increased peroxidase (PpPrx34) activity (Lehtonen et al., 2009(Lehtonen et al., , 2012a. Rapid ROS production was also evidenced in single moss leaf cells in contact with B. cinerea hyphae (Ponce de León et al., 2012). Other proteins related to oxidative stress such as cysteine peroxiredoxin, several thioredoxins, and copper/zinc superoxide dismutases are secreted to the apoplast after chitosan treatment (Lehtonen et al., 2014). The expression of a gene encoding a PpTSPO1 protein involved in mitochondrial tetrapyrrole transport is also induced by chitosan (Lehtonen et al., 2012a). Peroxidases, thioredoxins, and glutathione S-transferases transcript levels increase after treatment with bacterial culture filtrate of Pectobacterium carotovorum subsp. carotovorum (P.c. carotovorum) containing elicitors such as harpin and plant cell wall degrading enzymes secreted by this bacteria (Alvarez et al., 2016). The liverwort Dumortiera hirsute produces superoxide radicals by the action of apoplastic peroxidases and tyrosinases leading to an extracellular oxidative burst during rehydration (Li et al., 2010). Liverwort peroxidases may also protect plant tissues against pathogens by producing ROS and melanin-like compounds (Li et al., 2010). Functional analyses in P. patens knockouts have shown that PpPrx34 and PpTSPO1 are important for resistance against pathogens (Lehtonen et al., 2009(Lehtonen et al., , 2012a, indicating that ROS production is a common defense mechanisms in mosses and angiosperms. Chloroplasts are a rich source of ROS and these organelles are repositioned in P. patens cells toward the penetration site of several pathogens (Oliver et al., 2009;Lehtonen et al., 2012b;Reboledo et al., 2015; Figure 1). Although the precise role of chloroplast repositioning is at present unknown, plastidic ROS sources as well as the production of specific chloroplastic molecules could be part of the defense response to biotic stress. Interestingly, Nicotiana benthamiana chloroplasts send out dynamic tubular extensions called stromules during plant defense responses, which form connections with the nucleus (Caplan et al., 2015). These stromules transport prodefense signals such as H 2 O 2 into the nucleus during plant defense. Like in angiosperms, P. patens induces accumulation of cytoplasm and repositioning of the nucleus to pathogen contact sites, and rearrangement of the actin cytoskeleton after Phytophthora capsici (P. capsici) infection (Overdijk et al., 2016). Actin cytoskeleton is important for stromules formation and probably facilitates vesicle delivery to the infection site (Takemoto et al., 2006;Caplan et al., 2015). Further studies are needed to understand the role of chloroplasts repositioning during pathogen infection, including the relation with the pathogen penetration site, nucleus movement and cytoskeleton rearrangement.
The HR response contributes to plant resistance to biotrophic pathogens by confining the pathogen and limiting its growth, while necrotrophic pathogens actively stimulate the HR to enhance tissues colonization (Govrin and Levine, 2000;Jones and Dangl, 2006). P. patens cells infected or treated with elicitors of HR-inducing pathogens such as Botrytis cinerea (B. cinerea), Colletotrichum gloeosporioides (C. gloeosporioides), and P.c. carotovorum, induce a HR-like reaction (Ponce de León et al., 2007;Reboledo et al., 2015). Typical programmed cell death (PCD) hallmarks are present in elicitor-treated or pathogeninfected moss tissues, including chloroplasts breakdown, accumulation of autofluorescent compounds, cytoplasmic shrinkage, fragmentations of the nucleus and nuclease activities (Ponce de León et al., 2007, 2012Lawton and Saidasan, 2009).
In plants and algae, ROS generated during pathogen infection induce autophagy, which is a regulated process involved in recycling intracellular components (Pérez-Pérez et al., 2012). Autophagy is important in plant innate immunity and is involved in HR cell death (Lenz et al., 2011). A. thaliana autophagy deficient mutants (atg) develop spreading necrosis upon Alternaria brassicicola infection, and fungal toxin fumonisin B1 treatment (Lenz et al., 2011). However, atg mutant plants are more resistant to Pseudomonas syringae pv. tomato infection, suggesting that different mechanisms related to autophagy contribute to resistance after infection with different pathogens (Lenz et al., 2011). Like in A. thaliana, P. patens atg5 mutant shows early senescence under nutrient deficient conditions, which was partly explained by amino acid imbalance due to lack of cytoplasmic degradation by autophagy (Mukae et al., 2015). Although further studies are needed to understand the involvement of autophagy in moss defense against pathogens, preliminary results suggest that atg5 P. patens mutant plants are more susceptible than wild type plants to Sordariomycetes sp. fungus infection (Bressendorff, 2012).

FUNCTIONAL EVOLUTION OF HORMONAL PATHWAYS AND THEIR ROLE IN PLANT DEFENSE
The plant hormones salicylic acid (SA), jasmonic acid (JA), and ethylene have a central role in defense against pathogens in angiosperm. Additionally, other hormones such as abscisic acid (ABA), auxins, cytokinins, gibberellins, and brassinosteroids have emerged as regulators in plant immunity (Denancé et al., 2013). These hormones function in complex networks to regulate many aspects of plant resistance. In bryophytes, only few studies have related some of these hormones with defense against pathogens. However, genome-wide analyses of multiple species have been useful to gain insights into the origin and evolution of hormonal signaling components across the plant lineage. While algal genomes have the orthologs of only some components of the SA, JA, ABA signaling pathways, land plants have all the components, suggesting these signaling pathways might have developed during land colonization by plants (Wang et al., 2015). Orthologs of the Non-expresser of Pathogenesis-Related (NPR) genes have been found in P. patens (Wang et al., 2015). Studies in A. thaliana have identified NPR3 and NPR4 as SA receptors, while NPR1 is a key regulator in SAR, which is an important long distance signaling mechanism present in flowering plants that provides broad spectrum and long-lasting resistance to secondary infections throughout the plant (Fu et al., 2012). NPR3 and NPR4 are negative regulators of immunity, acting as substrate adaptors for the recruitment of NPR1 to an E3-ubiquitin ligase, which leads to NPR1 degradation by the proteasome. When NPR1 is not degraded and accumulates in the cytoplasm, it activates transcription factors such as WRKYs and TGAs, turning on defense genes encoding pathogenesis related (PR) proteins such as PR-1, PR-2, and PR-5 (Fu et al., 2012). SA levels increase in P. patens after B. cinerea infection and treatments with SA induce the expression of the defense gene phenylalanine ammonia lyase (PAL.1), suggesting that SA is involved in moss defense (Ponce de León et al., 2012). Consistently, SA application to P. patens colonies increases resistance to P. wasabiae (Andersson et al., 2005). Furthermore, the moss Amblystegium serpens activates a SAR-like response upon inoculation with Pythium irregulare (P. irregulare) or treatment with β-1,3 glucan increasing resistance against subsequent P. irregulare infections (Winter et al., 2014). This finding suggests that a SAR-like response exists before the divergence of mosses and vascular plants. Further research is needed to identify if SA and/or other signaling molecules participate in moss systemic responses. The use of SA-deficient transgenic moss plants, overexpressing a bacterial salicylate hydroxylase (encoded by NahG) that degrades SA, will help to elucidate if SA is involved in local and systemic moss defense.
Upon pathogen attack, JA-isoleucine accumulates rapidly in angiosperms promoting the interaction between the receptor COI1 and the repressor proteins JAZ, leading to the ubiquitination and subsequent degradation of JAZs through the proteasome (Chini et al., 2007). The degradation of JAZ relieves the repression on MYC transcription factors allowing the expression of JA-responsive genes. P. patens encode the orthologs of all JA signaling components, although moss and liverworts do not produce the oxygenated fatty acid (oxylipin) JA (Ponce de Yamamoto et al., 2015). The precursor 12-oxo-phytodienoic acid (OPDA) accumulates in P. patens tissues infected with P. irregulare, P. debaryanum (P. debaryanum) and B. cinerea (Oliver et al., 2009;Ponce de León et al., 2012). Similarly as in A. thaliana, OPDA could act as a signaling molecule in bryophytes leading to the induction of defense related genes (Taki et al., 2005). Surprisingly, P. patens and A. thaliana respond similarly to OPDA or JA, showing reduced rhizoid length and moss colony size, and reduced growth of seedlings and roots, respectively (Ponce de . Furthermore, jasmonates (JA and methyl jasmonate), and OPDA induce the expression of PAL.2, suggesting that all components needed for sensing JA and transducing the signal are present in basal land plants (Oliver et al., 2009). Since genes encoding proteins homologous to the JA-isoleucine receptor COI1 are present in the P. patens genome, other oxylipins than JA-isoleucine could be detected resulting in the activation of a defense response (Ponce de León et al., 2012). Interestingly, the binding sites are well conserved between COI orthologs of P. patens, the liverwort M. polymorpha and A. thaliana COI1 (Wang et al., 2015). The similarities between COI1 and the auxin receptor TIR1, at the amino acid sequence level and their mode of actions, as well as the fact that they probably originated from gene duplication events in bryophytes, suggest that JA and auxin pathways may have evolved in parallel from a common ancestor that later diverged by successive mutations to bind different hormones (Chico et al., 2008;Wang et al., 2010). Jasmonates could have evolved later in plant evolution as a result of selection pressures by insects and novel pathogens (Groen and Whiteman, 2014).
Abscisic acid plays a role in priming of callose biosynthesis after pathogen recognition and it induces the expression of defense genes in traqueophytes . While several studies have demonstrated the involvement of ABA in moss defense against abiotic stress (Cuming and Stevenson, 2015), little is known so far about the role played by ABA in defense against biotic stress. P. patens plants infected with B. cinerea show a small increase of ABA, however, it cannot be discarded that it is produced by the pathogen (Ponce de León et al., 2012). In P. patens, ABA induces the synthesis of defense proteins such as a RPM1-like R protein, a lipoxygenase (LOX) involved in the production of oxylipins, an intracellular pathogenesisrelated protein, a N-hydroxycinnamoyl/benzoyl transferase involved in the production of phytoalexin, a hydroxyprolinerich glycoprotein and a proline-rich protein involved in cell wall reinforcement, an ascorbate peroxidase and a peroxiredoxin, suggesting that ABA could be involved in moss defense against pathogens (Wang et al., 2010).
Bryophytes produce ethylene and the orthologs of the ethylene sensing and signaling components are present in several land plant genomes (Wang et al., 2015). P. patens has two putative genes encoding enzymes involved in the synthesis of the ethylene precursor, 1-aminocyclopropane 1 carboxylic acid (ACC) (Rensing et al., 2008), although their activities need to be confirmed. At least one of the seven ETR-like ethylene receptors binds ethylene in P. patens and this receptor is needed for a full ethylene water submergence moss response (Yasumura et al., 2012), suggesting that ethylene signaling was already present in mosses (Wang et al., 2015). Until present, there are only few lines of evidence suggesting that ethylene participates in moss defenses against pathogens. In P. patens, chitosan treatment and P. capsici infection induce the expression of the ethylene responsive transcription factor ERF2 and ERF5, respectively (Bressendorff et al., 2016;Overdijk et al., 2016), and treatment with ACC induces the expression of moss defense genes (Ponce de León and Montesano, 2013). In addition, several putative P. patens AP2/ERF encoding transcripts increase after bacterial elicitor treatment (Alvarez et al., 2016). Similarly, increased production of ethylene in A. thaliana leads to a cascade of transcription factors, including downstream ERF-like transcription factors that control the expression of several defense genes (Broekaert et al., 2006). Further studies using moss mutants in key components of the ethylene pathway, including ETRlike receptors will help to elucidate their role in defense against pathogens.
Auxin is produced in algae, and charophyte genomes have several components of auxin signaling and transport machinery, including ARFs, AUX/IAAs, and PINs, suggesting that AUX signaling originated in algae (Wang et al., 2015). In P. patens and M. polymorpha transcriptomes orthologs of all the signaling and transport core components were found, except for the liverwort Auxin-Binding Protein1, probably due to low coverage of the transcriptome (Wang et al., 2015). In P. patens, auxin levels increase in P.c. carotovorum elicitor-treated and Pythiuminfected tissues (Mittag et al., 2015;Alvarez et al., 2016). Furthermore, the use of a moss reporter line harboring the auxin inducible promoter GH3 (indole-3-acetic acid amino acid conjugate synthetases) fused to β-glucuronidase, shows that auxin signaling is induced in response to P.c. carotovorum elicitors and the pathogens P. irregulare, P. debaryanum, and C. gloeosporoides (Mittag et al., 2015;Reboledo et al., 2015;Alvarez et al., 2016). Like in A. thaliana where auxin induces the expression of genes and activates defense, auxin has been associated to P. patens defense against pathogens (Mittag et al., 2015). High auxin levels present in P. patens GH3 double knockout plants increase resistance against P. debaryanum (Mittag et al., 2015), and auxin induces the expression of the defense gene PAL.1 (Reboledo et al., 2015).
Orthologs of all the key cytokinin signaling components were identified in the charophyte and land plants genomes, suggesting that this hormonal signaling pathway originated in algae (Wang et al., 2015). For brassinosteroids, several of the signaling transduction components are present in algal genomes, although the key regulator BKI1 (Brassinosteroid Insensitive1Kinase Inhibitor) is found in angiosperms and not in algae, bryophytes, and lycophytes, suggesting that the brassinosteroids signaling pathway originated before the emergence of angiosperms (Wang et al., 2015). The roles played by cytokinin and brassinosteroids in moss defense have not been studied so far. Apparently, P. patens does not produce gibberellin, although it has part of the gibberellin biosynthetic pathway from geranylgeranyl diphosphate to entkaurenoic acid via ent-kaurene (Miyazaki et al., 2015). In angiosperms, the gibberellin signal is perceived and transduced by the GID1 gibberellin receptor/DELLA repressor pathway. A comparative analysis of gibberellin signaling components has shown that land plant genomes, except M. polymorpha probably due to low coverage of the transcriptome, encode the orthologs of all the core components of gibberellin signaling (Wang et al., 2015). While GID1 and DELLA do not interact with each other in P. patens, the lycophyte S. moellendorffii produces gibberellin and has functional GID1-DELLA signaling components (Yasumura et al., 2007). These results suggest that gibberellin signaling originated after the evolutionary split of bryophytes and traqueophytes. In contrast to flowering plants, P. patens has a bifunctional ent-kaurene synthase (KS) which has both ent-diphosphate synthase and ent-KS activities (Hayashi et al., 2006). P. patens mutant in this bifunctional CPS/KS has a phenotype with defective development that is restored with ent-kaurene or ent-kaurenoic acid (Miyazaki et al., 2015). The homologous diterpene cyclase (HpDTC1) from the moss Hypnum plumaeforme (H. plumaeforme) catalyzes a two-step cyclization reaction of geranylgeranyl diphosphate to syn-pimara-7,15-diene, which is a precursor for phytoalexin momilactones (Okada et al., 2016). HpDTC1 gene expression and momilactone accumulation is induced by chitosan and B. cinerea infection, suggesting a defensive role for momilactone biosynthesis in this moss (Okada et al., 2016). Interestingly, cis-OPDA but not JA enhanced HpDTC1 expression and momilactone accumulation, indicating that this oxylipin could regulate momilactone biosynthesis in H. plumaeforme (Okada et al., 2016). The tetracyclic diterpene 16α-hydroxykaurane (Von Schwartzenberg et al., 2004), is produced in P. patens tissues treated with chitosan, and transcripts levels encoding enzymes associated to kaurenes production increase, suggesting that this type of metabolites are also involved in moss defense (Lawton and Saidasan, 2009). Taken together, some of these hormonal signaling pathways represent innovations that probably played a role during the adaptation of terrestrial plants to microbial pathogens.

PHENYLPROPANOID AND OXYLIPINS PATHWAYS IN MOSS ADAPTATION TO BIOTIC STRESS
Genome analysis has revealed that P. patens has large gene families involved in protection against stresses (Rensing et al., 2008). Examples are the expansion of heat shock proteins 70, photoprotective early light-induced proteins, and enzymes involved in the phenylpropanoid pathway such as PALs and chalcone synthases (CHSs) (Lang et al., 2005;Rensing et al., 2008;Wolf et al., 2010). Interestingly, metabolic genes have been maintained in excess in P. patens following genome duplication, which took place about 45 million years ago, and this could be related to a more generalist metabolic gene complement that allow mosses to grow in more diverse habitats (Rensing et al., 2007). In addition, plant colonization required improved defense mechanisms to biotic challengers, including the production of novel or more specialized secondary metabolites. The phenylpropanoid pathway is probably one of the most relevant adaptation mechanisms that have evolved in early land plants to cope with abiotic stresses and to avoid invasion by already diversified microbial soil communities (Emiliani et al., 2009). In P. patens, the phenylpropanoid pathway is activated in response to B. cinerea, C. gloeosporioides, two Pythium species, P. capsici, chitin, and secreted elicitors of the pathogenic bacteria P.c. carotovorum (Ponce de León et al., 2007León et al., , 2012Oliver et al., 2009;Reboledo et al., 2015;Alvarez et al., 2016;Bressendorff et al., 2016;Overdijk et al., 2016). Several genes from the shikimate and phenylpropanoid pathways are induced after P.c. carotovorum elicitor treatment, including genes encoding a chorismate synthase, a chorismate mutase, PALs, a cinnamic acid 4-hydroxylase, a 4-coumarate CoA Ligase, CHSs, CHIs, a flavanone 3-hydroxylase and a cinnamyl alcohol dehydrogenase (Alvarez et al., 2016). The phenylpropanoid metabolism generates in angiosperms an enormous array of secondary metabolites with different roles in defense against pathogens, such as SA, flavonoids, isoflavonoids, phytoalexins, and lignins (Dixon and Paiva, 1995). Metabolites from the phenylpropanoid pathway have also been identified in P. patens, including caffeic acid, the cinnamate conjugate, caffeoylquinic acid, and the intermediate of lignin biosynthesis, chlorogenic acid (Erxleben et al., 2012). Cinnamic acid increases in P.c. carotovorum elicitor-treated moss tissues and genes encoding Dirigent proteins involved in monolignol coupling are also induced by these elicitors as well as by several pathogens (Ponce de León et al., 2012;Reboledo et al., 2015;Alvarez et al., 2016). This suggests that cinnamic acid-derived and monolignol-containing compounds are produced in P. patens after pathogen assault. The presence of lignin in bryophytes is still controversial and P. patens could use cell wall-bound phenolics that resemble lignin for defensive purposes. Lignin-like polymers have also been identified in cell walls of algae, other mosses, liverworts and hornworts where they are proposed to principally serve as a defense mechanism against microorganisms and UV radiation (Ligrone et al., 2008). Callose is a highmolecular weight β-(1,3)-glucan polymer involved in plant cell wall defenses and it is usually associated to phenolic compounds, polysaccharides and antimicrobial proteins forming cell wall appositions (papilla) (Luna et al., 2011). Cell wall appositions were identified in P. patens tissues inoculated with P. debaryanum, P. capsici, P. infestans (P. infestans) and Apiospora montagnei (Oliver et al., 2009;Lehtonen et al., 2012b;Overdijk et al., 2016). Similarly, papilla formation was observed in the moss Funaria hygrometrica inoculated with the fungus Atradidymella muscivora and Coniochaeta velutina (Davey et al., 2009(Davey et al., , 2010. Papilla-like structures are capable of blocking P. infestans and P. capsici penetration and further colonization, which are defensive barriers also observed in angiosperms (Overdijk et al., 2016).
Physcomitrella patens has two CHIs, which are enhancer of flavonoid production genes encoding type IV CHI proteins with no CHI activity (Ngaki et al., 2012), and the conversion of flavanone to flavone by flavone synthase seems to be absent in mosses (Koduri et al., 2010). Both CHI genes are induced in response to P.c. carotovorum elicitor treatment (Alvarez et al., 2016), suggesting that they participate in moss defense. Since flavonoids are important defense compounds in plant responses to pathogens and insects, further analyses are needed to characterize enzymatic activities of the pathway and to have a broader metabolic scenario in P. patens tissues under biotic stress conditions.
Physcomitrella patens has a high proportion of proteins involved in metabolism, and alternative metabolic pathways have been discovered in this moss (Lang et al., 2005). Interestingly, P. patens has a lipoxygenase with a hydroperoxidase and fatty acid chain-cleaving lyase activity, which uses C18-fatty acids (linolenic and linoleic acid) and C20-fatty acids (20:4 arachidonic and 20:5 eicosapentaenoic acids) as substrates yielding a broader range of oxylipins than angiosperms (Figure 2; Senger et al., 2005;Wichard et al., 2005;Anterola et al., 2009). The high abundance of long and very long chain fatty acids as well as the presence of oxylipins derived from C20:4 and C20:5 represent a metabolic difference between mosses and angiosperms that may provide a metabolic advantage to adapt to abiotic stress (Mikami and Hartmann, 2004), and probably to protect themselves against pathogen invasion. Mosses could have evolved the ability to produce C20-derived oxylipins for defensive purposes against pathogens, although further research is needed to determine if some of them have antimicrobial activities or induce defense genes. Like in angiosperms, pathways leading to oxylipin biosynthesis are activated in P. patens after pathogen assault or elicitor treatment, evidenced by increased endogenous levels of free unsaturated fatty acid and induced expression of several lipases, LOXs, allene oxide synthases, OPDA reductases and alpha-Dioxygenases (alpha-DOX) (Lehtonen et al., 2014;Ponce de León et al., 2015;Alvarez et al., 2016;Bressendorff et al., 2016). These evidences show that different LOXs and alpha-DOX pathways are activated after pathogen assault, leading to the production of different oxylipins that could play different roles in moss defense. Interestingly, LOX-derived oxylipins, synthesized from C20 and C18 polyunsaturated fatty acids, are also produced in algae, where they function in defense against an algal pathogen (Bouarab et al., 2004). Oxilipin profiling after pathogen assault and the generation of multiple knockouts plants in the different enzymes of these pathways will help to elucidate their role in P. patens defense. Like in angiosperms, moss tissues treated with pathogens increase production of oxylipins derived from the alpha-DOX, which reduce cellular damage caused by elicitors of P.c. carotovorum (Machado et al., 2015).

MOSS DEFENSE GENES AND GENE ACQUISITION BY HORIZONTAL TRANSFER
Different genes encode proteins that are involved in plant resistance against pathogens, and generally their expression is up-regulated when plants are attacked by bacteria, fungi, oomycetes, or viruses. In angiosperms, the expression of genes encoding PR proteins increases at the infection site and during SAR (van Loon et al., 2006). Several PR encoding genes are induced in P. patens after P.c. carotovorum elicitor treatment or pathogen infection, including PR-1, PR-10 and PR-3 (chitinase) (Ponce de León et al., 2007;Oliver et al., 2009;Alvarez et al., 2016). Moreover, thaumatin-like protein (PR-5) and chitinases are secreted to the extracellular space when moss plants grown in liquid medium are treated with a fungal elicitor (Lehtonen et al., 2014). Transcript levels of genes encoding transcription factors, betaine-aldehyde dehydrogenase and proteins involved in modification or reinforcement of the cell wall such as Dirigent proteins and xyloglucan:xyloglucosyl transferases also increase in moss tissues treated with bacterial secreted elicitors (Alvarez et al., 2016).
It was recently shown that chitin oligosaccharides induce the expression of P. patens chitinase, which inhibits the growth of the fungus Trichoderma viride, suggesting that this enzyme could play a role in moss defense (Kobaru et al., 2016). Moss lines overexpressing PpPR-10 are more resistant to P. irregulare infection and accumulate unknown cell wall polymers, probably reinforcing the cell walls and increasing resistance to pathogens (Castro et al., 2016). PpPR-10 binds cytokinin (Gonneau et al., 2001), however, the biological relevance of this binding capacity in plant defense against pathogens is currently unknown. Coexpression of PpPR-10 with moss genes encoding enzymes of the phenylpropanoid pathway, together with the fact that some PR-10 proteins from other plants are related to phenylpropanoids biosynthesis, suggest that PpPR-10 could also affect phenolic compounds repertoire in moss tissues (Castro et al., 2016). Other proteins with possible roles in defense are also secreted to the apoplast in chitosan-treated moss tissues, including a D-mannose binding lectin protein (Lehtonen et al., 2014), which could act by recognizing pathogen surfaces and activating a moss defense response similarly as in angiosperms (Hwang and Hwang, 2011).
Several genes present in P. patens and associated to defense were acquired by horizontal gene transfer, providing an adaptive advantage that probably played a role in defense against microbial colonization in basal land plants. Examples include the bifunctional CPS/KS enzyme, which probably came from fungi by horizontal gene transfer, and unifunctional CPS and KS evolved gradually in angiosperms (Mukherjee, 2015). PAL was probably transferred horizontally from soil bacteria (Emiliani et al., 2009), and the major intrinsic protein acquired from bacteria probably facilitates the transport of glycerol in P. patens (Gustavsson et al., 2005). The acquisition of glycerol transporters in plants could have been important for symbiotic fungi and bacteria, representing an effective mechanism for transferring carbon sources from the plant to microbes (Wei et al., 2004;Gustavsson et al., 2005). Several genes acquired in P. patens by horizontal transfer from bacteria and fungi are directly or indirectly related to plant defense responses and stress tolerance. A glutamate-cysteine ligase acquired from bacteria is involved in glutathione formation, which is important for plant disease resistance (Yue et al., 2012). Other acquired genes encoding proteins with possible roles in pathogen resistance include a killer toxin Protein, a bifunctional nuclease, and an acid phosphatase (Yue et al., 2012). These acquired genes could have expanded the metabolic capabilities of P. patens to better adapt to pathogen infection as well as to beneficial microbes. P. patens has 13% orphan genes (Zimmer et al., 2013), and some of them could represent innovative adaptive strategies to biotic stress. Consistently, some of these orphan genes were identified in a library of P.c. carotovorum elicitortreated P. patens tissues (Alvarez et al., 2016), suggesting that they could play a role in defense against pathogens. Due to the high number of orphan genes expressed in P. patens during abiotic stresses, species-or lineage-specific genes could have played an important role in the acquisition of abiotic stress tolerance in mosses (Khraiwesh et al., 2015). Further studies during biotic stress will help to elucidate if orphan genes could also play a role in moss resistance against pathogens.

CONCLUSION
As early land-growing plants, bryophytes like P. patens adapted cellular and molecular responses to interact with beneficial microbes and to cope with microbial pathogens during their life cycle. Comparative studies between beneficial and pathogenic interactions open the possibility to gain insights into the molecular mechanisms involved during coevolution, and how basal land plants have evolved the ability to distinguish between these different types of microbes. Further research is needed to identify and characterize more moss PRRs and possible targets of pathogen effector proteins that inhibit PTI or ETI mediated signaling in bryophytes. Studies conducted so far have shown that several defense mechanisms against pathogen infection are conserved among moss and angiosperms. Some of these mechanisms emerged when plants colonized land, probably influencing fitness and conferring adaptive advantages to biotic stressors. Certain precursor stages of defense mechanisms were already present in the algal ancestor of plants, while several moss genes associated to defense were acquired with terrestrialization, either via horizontal gene transfer from microbes or by duplication events and probably further diversification and functional specialization. The expansion of metabolic gene families, the production of novel metabolites, and the presence of orphan genes with possible roles in P. patens defense, suggest that 470 million years of divergent evolution lead to the development of alternative defense strategies for survival in extant moss and angiosperms. The efficient transient transformation and generation of P. patens knockouts in several genes using CRISPR-Cas9 open new perspectives for functional analysis of gene families that participate in moss defense (Lopez-Obando et al., 2016). Furthermore, genome wide expression analysis, generation of knockout moss plants, and comparative studies with species representing different plant lineages, will provide new insights into the evolutionary dynamics, diversity and origin of genes, including moss-specific genes, involved in defense against biotic stressors and adaptation mechanisms that played a crucial role in colonization by land plant ancestors.

AUTHOR CONTRIBUTIONS
IPDL conceived the review and wrote the manuscript. MM participated in the discussions and helped to draft and write the manuscript.

ACKNOWLEDGMENTS
This work was supported by Programa de Desarrollo de las Ciencias Básicas (PEDECIBA), Uruguay. Both authors read and approved the final version of this manuscript.