Variation in shoot architecture traits and their relationship to canopy coverage and light interception in soybean (Glycine max).

BACKGROUND: In soybeans, faster canopy coverage (CC) is a highly desirable trait but a fully covered canopy is unfavorable to light interception at lower levels in the canopy with most of the incident radiation intercepted at the top of the canopy. Shoot architecture that influences CC is well studied in crops such as maize and wheat, and altering architectural traits has resulted in enhanced yield. However, in soybeans the study of shoot architecture has not been as extensive. RESULTS: This study revealed significant differences in CC among the selected soybean accessions. The rate of CC was found to decrease at the beginning of the reproductive stage (R1) followed by an increase during the R2-R3 stages. Most of the accessions in the study achieved maximum rate of CC between R2-R3 stages. We measured Light interception (LI), defined here as the ratio of Photosynthetically Active Radiation (PAR) transmitted through the canopy to the incoming PAR or the radiation above the canopy. LI was found to be significantly correlated with CC parameters, highlighting the relationship between canopy structure and light interception. The study also explored the impact of plant shape on LI and CO2 assimilation. Plant shape was characterized into distinct quantifiable parameters and by modeling the impact of plant shape on LI and CO2 assimilation, we found that plants with broad and flat shapes at the top maybe more photosynthetically efficient at low light levels, while conical shapes were likely more advantageous when light was abundant. Shoot architecture of plants in this study was described in terms of whole plant, branching and leaf-related traits. There was significant variation for the shoot architecture traits between different accessions, displaying high reliability. We found that that several shoot architecture traits such as plant height, and leaf and internode-related traits strongly influenced CC and LI. CONCLUSION: In conclusion, this study provides insight into the relationship between soybean shoot architecture, canopy coverage, and light interception. It demonstrates that novel shoot architecture traits we have defined here are genetically variable, impact CC and LI and contribute to our understanding of soybean morphology. Correlations between different architecture traits, CC and LI suggest that it is possible to optimize soybean growth without compromising on light transmission within the soybean canopy. In addition, the study underscores the utility of integrating low-cost 2D phenotyping as a practical and cost-effective alternative to more time-intensive 3D or high-tech low-throughput methods. This approach offers a feasible means of studying basic shoot architecture traits at the field level, facilitating a broader and efficient assessment of plant morphology.

Pathogen-driven coevolution across the CBP60 plant immune regulator subfamilies confers resilience on the regulator module.

Components of the plant immune signaling network need mechanisms that confer resilience against fast-evolving pathogen effectors that target them. Among eight Arabidopsis CaM-Binding Protein (CBP) 60 family members, AtCBP60g and AtSARD1 are partially functionally redundant, major positive immune regulators, and AtCBP60a is a negative immune regulator. We investigated possible resilience-conferring evolutionary mechanisms among the CBP60a, CBP60g and SARD1 immune regulatory subfamilies. Phylogenetic analysis was used to investigate the times of CBP60 subfamily neofunctionalization. Then, using the pairwise distance rank based on the newly developed analytical platform Protein Evolution Analysis in a Euclidean Space (PEAES), hypotheses of specific coevolutionary mechanisms that could confer resilience on the regulator module were tested. The immune regulator subfamilies diversified around the time of angiosperm divergence and have been evolving very quickly. We detected significant coevolutionary interactions across the immune regulator subfamilies in all of 12 diverse core eudicot species lineages tested. The coevolutionary interactions were consistent with the hypothesized coevolution mechanisms. Despite their unusually fast evolution, members across the CBP60 immune regulator subfamilies have influenced the evolution of each other long after their diversification in a way that could confer resilience on the immune regulator module against fast-evolving pathogen effectors.

A plant effector-triggered immunity signaling sector is inhibited by pattern-triggered immunity.

Since signaling machineries for two modes of plant-induced immunity, pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), extensively overlap, PTI and ETI signaling likely interact. In an Arabidopsis quadruple mutant, in which four major sectors of the signaling network, jasmonate, ethylene, PAD4, and salicylate, are disabled, the hypersensitive response (HR) typical of ETI is abolished when the Pseudomonas syringae effector AvrRpt2 is bacterially delivered but is intact when AvrRpt2 is directly expressed in planta These observations led us to discovery of a network-buffered signaling mechanism that mediates HR signaling and is strongly inhibited by PTI signaling. We named this mechanism the ETI-Mediating and PTI-Inhibited Sector (EMPIS). The signaling kinetics of EMPIS explain apparently different plant genetic requirements for ETI triggered by different effectors without postulating different signaling machineries. The properties of EMPIS suggest that information about efficacy of the early immune response is fed back to the immune signaling network, modulating its activity and limiting the fitness cost of unnecessary immune responses.

The highly buffered Arabidopsis immune signaling network conceals the functions of its components.

Plant immunity protects plants from numerous potentially pathogenic microbes. The biological network that controls plant inducible immunity must function effectively even when network components are targeted and disabled by pathogen effectors. Network buffering could confer this resilience by allowing different parts of the network to compensate for loss of one another's functions. Networks rich in buffering rely on interactions within the network, but these mechanisms are difficult to study by simple genetic means. Through a network reconstitution strategy, in which we disassemble and stepwise reassemble the plant immune network that mediates Pattern-Triggered-Immunity, we have resolved systems-level regulatory mechanisms underlying the Arabidopsis transcriptome response to the immune stimulant flagellin-22 (flg22). These mechanisms show widespread evidence of interactions among major sub-networks-we call these sectors-in the flg22-responsive transcriptome. Many of these interactions result in network buffering. Resolved regulatory mechanisms show unexpected patterns for how the jasmonate (JA), ethylene (ET), phytoalexin-deficient 4 (PAD4), and salicylate (SA) signaling sectors control the transcriptional response to flg22. We demonstrate that many of the regulatory mechanisms we resolved are not detectable by the traditional genetic approach of single-gene null-mutant analysis. Similar to potential pathogenic perturbations, null-mutant effects on immune signaling can be buffered by the network.

Different Modes of Negative Regulation of Plant Immunity by Calmodulin-Related Genes.

Plant immune responses activated through the perception of microbe-associated molecular patterns, leading to pattern-triggered immunity, are tightly regulated. This results in low immune responses in the absence of pathogens and a rapid return to the resting state following an activation event. Here, we show that two CALMODULIN-LIKE genes, CML46 and CML47, negatively regulate salicylic acid accumulation and immunity in Arabidopsis (Arabidopsis thaliana). The double mutant cml46 cml47 is highly resistant to the pathogen Pseudomonas syringae pv maculicola (Pma). The effects of cml46 cml47 on Pma growth are genetically additive to that of cbp60a, a known negative regulator in the CALMODULIN-BINDING PROTEIN60 (CBP60) family. Transcriptome profiling revealed the effects of cbp60a and cml46 cml47 on both common and separate sets of genes, with the majorities of these differentially expressed genes being Pma responsive. CBP60g, a positive regulator of immunity in the CBP60 family, was found to be transcriptionally regulated by CBP60a, CML46, and CML47 Analysis of the flg22-induced mRNA levels of CBP60g in cbp60a and cml46 cml47 revealed that cml46 cml47 plants have higher induced expression while cbp60a plants retain elevated levels longer than wild-type plants. Assays for the effect of flg22 treatment on Pma growth showed that the effect is stronger in cml46 cml47 plants and lasts longer in cbp60a plants. Thus, the expression pattern of CBP60g is reflected in flg22-induced resistance to Pma.

Involvement of Adapter Protein Complex 4 in Hypersensitive Cell Death Induced by Avirulent Bacteria.

Plant immunity to avirulent bacterial pathogens is associated with subcellular membrane dynamics including fusion between the vacuolar and plasma membranes, resulting in hypersensitive cell death. Here, we report that ADAPTOR PROTEIN COMPLEX-4 (AP-4) subunits are involved in plant immunity associated with hypersensitive cell death. We isolated a mutant with a defect in resistance to an avirulent strain of Pseudomonas syringae pv. tomato (Pto) DC3000 avrRpm1 from a vacuolar protein sorting mutant library of Arabidopsis (Arabidopsis thaliana). The mutant was identical to gfs4-1, which has a mutation in the gene encoding the AP-4 subunit AP4B. Thus, we focused on AP4B and another subunit, AP4E. All of the mutants (ap4b-3, ap4b-4, ap4e-1, and ap4e-2) were defective in hypersensitive cell death and resistance to Pto DC3000 with the type III effector AvrRpm1 or AvrRpt2, both of which are recognized on the plasma membrane, while they showed slightly enhanced susceptibility to the type-III-secretion-deficient P. syringae strain hrcC On the other hand, both ap4b-3 and ap4b-4 showed no defect in resistance to Pto DC3000 with the type III effector AvrRps4, which is recognized in the cytosol and does not induce hypersensitive cell death. Upon infection with Pto DC3000 avrRpt2, the ap4b-3 and ap4b-4 leaf cells did not show fusion between vacuolar and plasma membranes, whereas the wild-type leaf cells did. These results suggest that AP-4 contributes to cell death-associated immunity, possibly via membrane fusion, after type III effector-recognition on the plasma membrane.

Pectin Biosynthesis Is Critical for Cell Wall Integrity and Immunity in Arabidopsis thaliana.

Plant cell walls are important barriers against microbial pathogens. Cell walls of Arabidopsis thaliana leaves contain three major types of polysaccharides: cellulose, various hemicelluloses, and pectins. UDP-D-galacturonic acid, the key building block of pectins, is produced from the precursor UDP-D-glucuronic acid by the action of glucuronate 4-epimerases (GAEs). Pseudomonas syringae pv maculicola ES4326 (Pma ES4326) repressed expression of GAE1 and GAE6 in Arabidopsis, and immunity to Pma ES4326 was compromised in gae6 and gae1 gae6 mutant plants. These plants had brittle leaves and cell walls of leaves had less galacturonic acid. Resistance to specific Botrytis cinerea isolates was also compromised in gae1 gae6 double mutant plants. Although oligogalacturonide (OG)-induced immune signaling was unaltered in gae1 gae6 mutant plants, immune signaling induced by a commercial pectinase, macerozyme, was reduced. Macerozyme treatment or infection with B. cinerea released less soluble uronic acid, likely reflecting fewer OGs, from gae1 gae6 cell walls than from wild-type Col-0. Although both OGs and macerozyme-induced immunity to B. cinerea in Col-0, only OGs also induced immunity in gae1 gae6. Pectin is thus an important contributor to plant immunity, and this is due at least in part to the induction of immune responses by soluble pectin, likely OGs, that are released during plant-pathogen interactions.

WRKY70 prevents axenic activation of plant immunity by direct repression of SARD1.

SARD1 is an activator of plant immunity that promotes production of the hormone salicylic acid (SA) and activation of defense gene expression. SARD1 itself is strongly inducible by infection. Here, we investigated the transcriptional control of SARD1. We used yeast one-hybrid assays to identify WRKY70. The WRKY70 binding site was defined using electrophoretic mobility shift assays, and its importance was investigated using an Arabidopsis thaliana protoplast system. The effect of wrky70 mutations was studied by measurements of pathogen growth, SA concentrations, and gene expression by RNA-seq. WRKY70 binds to a GACTTTT motif in the SARD1 promoter in yeast and Arabidopsis protoplasts. Plants with wrky70 mutations have elevated expression of SARD1 in the absence of pathogens, but not when infected. Expression profiling revealed that WRKY70 represses many pathogen-inducible genes in the absence of pathogens, yet is required for activation of many other pathogen-inducible genes in infected plants. The GACTTTT motif is enriched in the promoters of both these gene sets, and conserved in SARD1 orthologs within the Brassicaceae. WRKY70 represses SARD1 by binding the motif GACTTTT in the absence of pathogens. Conservation of the WRKY70 binding among the Brassicaceae suggests that WRKY70 repression of SARD1 is important for fitness.

Foundational and Translational Research Opportunities to Improve Plant Health.

Reader Comments | Submit a Comment The white paper reports the deliberations of a workshop focused on biotic challenges to plant health held in Washington, D.C. in September 2016. Ensuring health of food plants is critical to maintaining the quality and productivity of crops and for sustenance of the rapidly growing human population. There is a close linkage between food security and societal stability; however, global food security is threatened by the vulnerability of our agricultural systems to numerous pests, pathogens, weeds, and environmental stresses. These threats are aggravated by climate change, the globalization of agriculture, and an over-reliance on nonsustainable inputs. New analytical and computational technologies are providing unprecedented resolution at a variety of molecular, cellular, organismal, and population scales for crop plants as well as pathogens, pests, beneficial microbes, and weeds. It is now possible to both characterize useful or deleterious variation as well as precisely manipulate it. Data-driven, informed decisions based on knowledge of the variation of biotic challenges and of natural and synthetic variation in crop plants will enable deployment of durable interventions throughout the world. These should be integral, dynamic components of agricultural strategies for sustainable agriculture.

The µ Subunit of Arabidopsis Adaptor Protein-2 Is Involved in Effector-Triggered Immunity Mediated by Membrane-Localized Resistance Proteins.

Endocytosis has been suggested to be important in the cellular processes of plant immune responses. However, our understanding of its role during effector-triggered immunity (ETI) is still limited. We have previously shown that plant endocytosis, especially clathrin-coated vesicle formation at the plasma membrane, is mediated by the adaptor protein-2 (AP-2) complex and that loss of the µ subunit of AP-2 (AP2M) affects plant growth and floral organ development. Here, we report that AP2M is required for full-strength ETI mediated by the disease resistance (R) genes RPM1 and RPS2 in Arabidopsis. Reduced ETI was observed in an ap2m mutant plant, measured by growth of Pseudomonas syringae pv. tomato DC3000 strains carrying the corresponding effector genes avrRpm1 or avrRpt2 and by hypersensitive cell death response and defense gene expression triggered by these strains. In contrast, RPS4-mediated ETI and its associated immune responses were not affected by the ap2m mutation. While RPM1 and RPS2 are localized to the plasma membrane, RPS4 is localized to the cytoplasm and nucleus. Our results suggest that AP2M is involved in ETI mediated by plasma membrane-localized R proteins, possibly by mediating endocytosis of the immune receptor complex components from the plasma membrane.

BR-SIGNALING KINASE1 physically associates with FLAGELLIN SENSING2 and regulates plant innate immunity in Arabidopsis.

Pathogen-associated molecular pattern (PAMP)-trigged immunity (PTI) is the first defensive line of plant innate immunity and is mediated by pattern recognition receptors. Here, we show that a mutation in BR-SIGNALING KINASE1 (BSK1), a substrate of the brassinosteroid (BR) receptor BRASSINOSTEROID INSENSITIVE1, suppressed the powdery mildew resistance caused by a mutation in ENHANCED DISEASE RESISTANCE2, which negatively regulates powdery mildew resistance and programmed cell death, in Arabidopsis thaliana. A loss-of-function bsk1 mutant displayed enhanced susceptibility to virulent and avirulent pathogens, including Golovinomyces cichoracearum, Pseudomonas syringae, and Hyaloperonospora arabidopsidis. The bsk1 mutant also accumulated lower levels of salicylic acid upon infection with G. cichoracearum and P. syringae. BSK1 belongs to a receptor-like cytoplasmic kinase family and displays kinase activity in vitro; this kinase activity is required for its function. BSK1 physically associates with the PAMP receptor FLAGELLIN SENSING2 and is required for a subset of flg22-induced responses, including the reactive oxygen burst, but not for mitogen-activated protein kinase activation. Our data demonstrate that BSK1 is involved in positive regulation of PTI. Together with previous findings, our work indicates that BSK1 represents a key component directly involved in both BR signaling and plant immunity.

Dual regulation of gene expression mediated by extended MAPK activation and salicylic acid contributes to robust innate immunity in Arabidopsis thaliana.

Network robustness is a crucial property of the plant immune signaling network because pathogens are under a strong selection pressure to perturb plant network components to dampen plant immune responses. Nevertheless, modulation of network robustness is an area of network biology that has rarely been explored. While two modes of plant immunity, Effector-Triggered Immunity (ETI) and Pattern-Triggered Immunity (PTI), extensively share signaling machinery, the network output is much more robust against perturbations during ETI than PTI, suggesting modulation of network robustness. Here, we report a molecular mechanism underlying the modulation of the network robustness in Arabidopsis thaliana. The salicylic acid (SA) signaling sector regulates a major portion of the plant immune response and is important in immunity against biotrophic and hemibiotrophic pathogens. In Arabidopsis, SA signaling was required for the proper regulation of the vast majority of SA-responsive genes during PTI. However, during ETI, regulation of most SA-responsive genes, including the canonical SA marker gene PR1, could be controlled by SA-independent mechanisms as well as by SA. The activation of the two immune-related MAPKs, MPK3 and MPK6, persisted for several hours during ETI but less than one hour during PTI. Sustained MAPK activation was sufficient to confer SA-independent regulation of most SA-responsive genes. Furthermore, the MPK3 and SA signaling sectors were compensatory to each other for inhibition of bacterial growth as well as for PR1 expression during ETI. These results indicate that the duration of the MAPK activation is a critical determinant for modulation of robustness of the immune signaling network. Our findings with the plant immune signaling network imply that the robustness level of a biological network can be modulated by the activities of network components.

Putative Serine Protease Effectors of Clavibacter michiganensis Induce a Hypersensitive Response in the Apoplast of Nicotiana Species.

Clavibacter michiganensis subspp. michiganensis and sepedonicus cause diseases on solanaceous crops. The genomes of both subspecies encode members of the pat-1 family of putative serine proteases known to function in virulence on host plants and induction of hypersensitive responses (HR) on nonhosts. One gene of this family in C. michiganensis subsp. sepedonicus, chp-7, is required for triggering HR in Nicotiana tabacum. Here, further investigation revealed that mutation of the putative catalytic serine residue at position 232 to threonine abolished the HR induction activity of Chp-7, suggesting that enzymatic activity is required. Purified Chp-7 triggered an HR in N. tabacum leaves in the absence of the pathogen, indicating Chp-7 itself is the HR elicitor from C. michiganensis subsp. sepedonicus. Ectopic expression of chp-7 constructs in N. tabacum leaves revealed that Chp-7 targeted to the apoplast triggered an HR while cytoplasmic Chp-7 did not, indicating that Chp-7 induces the HR in the apoplast of N. tabacum leaves. Chp-7 also induced HR in N. sylvestris, a progenitor of N. tabacum, but not in other Nicotiana species tested. ChpG, a related protein from C. michiganensis subsp. michiganensis, also triggered HR in N. tabacum and N. sylvestris. Unlike Chp-7, ChpG triggered HR in N. clevelandii and N. glutinosa.

The receptor-like cytoplasmic kinase PCRK1 contributes to pattern-triggered immunity against Pseudomonas syringae in Arabidopsis thaliana.

In this paper we describe PATTERN-TRIGGERED IMMUNITY (PTI) COMPROMISED RECEPTOR-LIKE CYTOPLASMIC KINASE 1 (PCRK1) of Arabidopsis thaliana, an RLCK that is important for defense against the pathogen Pseudomonas syringae pv. maculicola ES4326 (Pma ES4326). We examined defense responses such as bacterial growth, production of reactive oxygen species (ROS) and callose deposition in pcrk1 mutant plants to determine the role of PCRK1 during pathogen infection. Expression of PCRK1 was induced following pathogen infection. Pathogen growth was significantly higher in pcrk1 mutant lines than in wild-type Col-0. Mutant pcrk1 plants showed reduced pattern-triggered immunity (PTI) against Pma ES4326 after pretreatment with peptides derived from flagellin (flg22), elongation factor-Tu (elf18), or an endogenous protein (pep1). Deposition of callose was reduced in pcrk1 plants, indicating a role of PCRK1 in activation of early immune responses. A PCRK1 transgene containing a mutation in a conserved lysine residue important for phosphorylation activity of kinases (K118E) failed to complement a pcrk1 mutant for the Pma ES4326 growth phenotype. Our study shows that PCRK1 plays an important role during PTI and that a conserved lysine residue in the putative kinase domain is important for PCRK1 function.

Arabidopsis PECTIN METHYLESTERASEs contribute to immunity against Pseudomonas syringae.

Pectins, major components of dicot cell walls, are synthesized in a heavily methylesterified form in the Golgi and are partially deesterified by pectin methylesterases (PMEs) upon export to the cell wall. PME activity is important for the virulence of the necrotrophic fungal pathogen Botrytis cinerea. Here, the roles of Arabidopsis PMEs in pattern-triggered immunity and immune responses to the necrotrophic fungus Alternaria brassicicola and the bacterial hemibiotroph Pseudomonas syringae pv maculicola ES4326 (Pma ES4326) were studied. Plant PME activity increased during pattern-triggered immunity and after inoculation with either pathogen. The increase of PME activity in response to pathogen treatment was concomitant with a decrease in pectin methylesterification. The pathogen-induced PME activity did not require salicylic acid or ethylene signaling, but was dependent on jasmonic acid signaling. In the case of induction by A. brassicicola, the ethylene response factor, but not the MYC2 branch of jasmonic acid signaling, contributed to induction of PME activity, whereas in the case of induction by Pma ES4326, both branches contributed. There are 66 PME genes in Arabidopsis, suggesting extensive genetic redundancy. Nevertheless, selected pme single, double, triple and quadruple mutants allowed significantly more growth of Pma ES4326 than wild-type plants, indicating a role of PMEs in resistance to this pathogen. No decreases in total PME activity were detected in these pme mutants, suggesting that the determinant of immunity is not total PME activity; rather, it is some specific effect of PMEs such as changes in the pattern of pectin methylesterification.

The CALMODULIN-BINDING PROTEIN60 family includes both negative and positive regulators of plant immunity.

Two members of the eight-member CALMODULIN-BINDING PROTEIN60 (CBP60) gene family, CBP60g and SYSTEMIC ACQUIRED RESISTANCE DEFICIENT1 (SARD1), encode positive regulators of plant immunity that promote the production of salicylic acid (SA) and affect the expression of SA-dependent and SA-independent defense genes. Here, we investigated the other six family members in Arabidopsis (Arabidopsis thaliana). Only cbp60a mutations affected growth of the bacterial pathogen Pseudomonas syringae pv maculicola ES4326. In contrast to cbp60g and sard1 mutations, cbp60a mutations reduced pathogen growth, indicating that CBP60a is a negative regulator of immunity. Bacterial growth was increased by cbp60g only in the presence of CBP60a, while the increase in growth due to sard1 was independent of CBP60a, suggesting that the primary function of CBP60g may be to counter the repressive effect of CBP60a. In the absence of pathogen, levels of SA as well as of several SA-dependent and SA-independent pathogen-inducible genes were higher in cbp60a plants than in the wild type, suggesting that the enhanced resistance of cbp60a plants may result from the activation of immune responses prior to pathogen attack. CBP60a bound calmodulin, and the calmodulin-binding domain was defined at the C-terminal end of the protein. Transgenes encoding mutant versions of CBP60a lacking the ability to bind calmodulin failed to complement null cbp60a mutations, indicating that calmodulin-binding ability is required for the immunity-repressing function of CBP60a. Regulation at the CBP60 node involves negative regulation by CBP60a as well as positive regulation by CBP60g and SARD1, providing multiple levels of control over the activation of immune responses.

Arabidopsis lysin-motif proteins LYM1 LYM3 CERK1 mediate bacterial peptidoglycan sensing and immunity to bacterial infection.

Recognition of microbial patterns by host pattern recognition receptors is a key step in immune activation in multicellular eukaryotes. Peptidoglycans (PGNs) are major components of bacterial cell walls that possess immunity-stimulating activities in metazoans and plants. Here we show that PGN sensing and immunity to bacterial infection in Arabidopsis thaliana requires three lysin-motif (LysM) domain proteins. LYM1 and LYM3 are plasma membrane proteins that physically interact with PGNs and mediate Arabidopsis sensitivity to structurally different PGNs from gram-negative and gram-positive bacteria. lym1 and lym3 mutants lack PGN-induced changes in transcriptome activity patterns, but respond to fungus-derived chitin, a pattern structurally related to PGNs, in a wild-type manner. Notably, lym1, lym3, and lym3 lym1 mutant genotypes exhibit supersusceptibility to infection with virulent Pseudomonas syringae pathovar tomato DC3000. Defects in basal immunity in lym3 lym1 double mutants resemble those observed in lym1 and lym3 single mutants, suggesting that both proteins are part of the same recognition system. We further show that deletion of CERK1, a LysM receptor kinase that had previously been implicated in chitin perception and immunity to fungal infection in Arabidopsis, phenocopies defects observed in lym1 and lym3 mutants, such as peptidoglycan insensitivity and enhanced susceptibility to bacterial infection. Altogether, our findings suggest that plants share with metazoans the ability to recognize bacterial PGNs. However, as Arabidopsis LysM domain proteins LYM1, LYM3, and CERK1 form a PGN recognition system that is unrelated to metazoan PGN receptors, we propose that lineage-specific PGN perception systems have arisen through convergent evolution.

An averaging model for analysis and interpretation of high-order genetic interactions.

While combinatorial genetic data collection from biological systems in which quantitative phenotypes are controlled by active and inactive alleles of multiple genes (multi-gene systems) is becoming common, a standard analysis method for such data has not been established. The currently common approaches have three major drawbacks. First, although it is a long tradition in genetics, modeling the effect of an inactive allele (a null mutant allele) contrasted against that of the active allele (the wild-type allele) is not suitable for mechanistic understanding of multi-gene systems. Second, a commonly-used additive model (ANOVA with interaction) mathematically fails in estimation of interactions among more than two genes when the phenotypic response is not linear. Third, interpretation of higher-order interactions defined by an additive model is not intuitive. I derived an averaging model based on algebraic principles to solve all these problems within the framework of a general linear model. In the averaging model: the effect of the active allele is contrasted against the effect of the inactive allele for easier mechanistic interpretations; there is mathematical stability in estimation of higher-order interactions even when the phenotypic response is not linear; and interpretations of higher-order interactions are intuitive and consistent-interactions are defined as the mean effects of the last active genes added to the system. Thus, the key outcomes of this study are development of the averaging model, which is suitable for analysis of multi-gene systems, and a new, intuitive, and mathematically and interpretationally consistent definition of a genetic interaction, which is central to the averaging model.

Decomposition of dynamic transcriptomic responses during effector-triggered immunity reveals conserved responses in two distinct plant cell populations.

Rapid plant immune responses in the appropriate cells are needed for effective defense against pathogens. Although transcriptome analysis is often used to describe overall immune responses, collection of transcriptome data with sufficient resolution in both space and time is challenging. We reanalyzed public Arabidopsis time-course transcriptome data obtained after low-dose inoculation with a Pseudomonas syringae strain expressing the effector AvrRpt2, which induces effector-triggered immunity in Arabidopsis. Double-peak time-course patterns are prevalent among thousands of upregulated genes. We implemented a multi-compartment modeling approach to decompose the double-peak pattern into two single-peak patterns for each gene. The decomposed peaks reveal an "echoing" pattern: the peak times of the first and second peaks correlate well across most upregulated genes. We demonstrated that the two peaks likely represent responses of two distinct cell populations that respond either cell autonomously or indirectly to AvrRpt2. Thus, the peak decomposition has extracted spatial information from the time-course data. The echoing pattern also indicates a conserved transcriptome response with different initiation times between the two cell populations despite different elicitor types. A gene set highly overlapping with the conserved gene set is also upregulated with similar kinetics during pattern-triggered immunity. Activation of a WRKY network via different entry-point WRKYs can explain the similar but not identical transcriptome responses elicited by different elicitor types. We discuss potential benefits of the properties of the WRKY activation network as an immune signaling network in light of pressure from rapidly evolving pathogens.