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.

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.

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.

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.

Environmental Association Identifies Candidates for Tolerance to Low Temperature and Drought.

Barley (Hordeum vulgare ssp. vulgare) is cultivated from the equator to the Arctic Circle. The wild progenitor species, Hordeum vulgare ssp. spontaneum, occupies a relatively narrow latitudinal range (∼30 - 40° N) primarily at low elevation (< 1,500 m). Adaptation to the range of cultivation has occurred over ∼8,000 years. The genetic basis of adaptation is amenable to study through environmental association. An advantage of environmental association in a well-characterized crop is that many loci that contribute to climatic adaptation and abiotic stress tolerance have already been identified. This provides the opportunity to determine if environmental association approaches effectively identify these loci of large effect. Using published genotyping from 7,864 SNPs in 803 barley landraces, we examined allele frequency differentiation across multiple partitions of the data and mixed model associations relative to bioclimatic variables. Using newly generated resequencing data from a subset of these landraces, we tested for linkage disequilibrium (LD) between SNPs queried in genotyping and SNPs in neighboring loci. Six loci previously reported to contribute to adaptive differences in flowering time and abiotic stress in barley and six loci previously identified in other plant species were identified in our analyses. In many cases, patterns of LD are consistent with the causative variant occurring in the immediate vicinity of the queried SNP. The identification of barley orthologs to well-characterized genes may provide a new understanding of the nature of adaptive variation and could permit a more targeted use of potentially adaptive variants in barley breeding and germplasm improvement.

Review: Plant immune signaling from a network perspective.

For predictive biology of a large and complex network, important mechanistic information consists of network topology, signal-convergence rules, and signal dynamics. In practice, the key to enabling predictive modeling of a complex network is reducing it sufficiently to allow modeling without omitting important factors affecting network behavior. Here I argue that the plant immune signaling network must have high levels of resilience and tunability based on the fundamental facts that plants are evolutionarily disadvantaged relative to microbial pathogens and that unnecessary immune response is bad for plants. By reducing the middle part of the immune signaling network to a four-sector network, we previously showed that the middle part indeed has high levels of resilience and tunability and from what signaling strategies the network properties emerge. The ability to comprehensively reconstitute the reduced network enabled a reductionist approach to a resilient network, which was crucial for revealing signaling mechanisms concealed by network resilience. However, this four-sector network may be reduced too much to make a predictive model with relatively simple signal-convergence rules and probably requires deconvolution to generate a highly predictive network model with a set of simple signal-convergence rules.

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.

Quantification of Plant Cell Death by Electrolyte Leakage Assay.

We describe a protocol to measure the electrolyte leakage from plant tissues, resulting from loss of cell membrane integrity, which is a common definition of cell death. This simple protocol is designed to measure the electrolyte leakage from a tissue sample over a time course, so that the extent of cell death in the tissue can be monitored dynamically. In addition, it is easy to handle many tissue samples in parallel, which allows a high level of biological replication. Although the protocol is exemplified by cell death in Arabidopsis in response to pathogen challenge, it is easily applicable to other types of plant cell death.

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.

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.

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.

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.

Network Reconstitution for Quantitative Subnetwork Interaction Analysis.

A fundamental task in systems biology is to quantify the contributions of the systems' parts and their interactions. Here I describe a powerful concept and tool for this purpose: network reconstitution. Genotypes of an organism that represent all possible combinations of the subnetworks in question will be quantitatively phenotyped. The quantitative phenotype data is analyzed using an R script to obtain estimates for single subnetwork contributions and their interactions.

Nup82 functions redundantly with Nup136 in a salicylic acid-dependent defense response of Arabidopsis thaliana.

The nuclear pore complex (NPC) comprises more than 30 nucleoporins (Nups). NPC mediates macromolecular trafficking between the nucleoplasm and the cytoplasm, but specific roles of individual Nups are poorly understood in higher plants. Here, we show that the novel nucleoporin unique to angiosperm plants (designated as Nup82) functions in a salicylic acid-dependent defense in a redundant manner with Nup136, which is a component of the nuclear basket in the NPC. Arabidopsis thaliana Nup82 had a similar amino acid sequence to the N-terminal half of Nup136 and a Nup82-GFP fusion was localized on the nuclear envelope. Immunoprecipitation and bimolecular fluorescence complementation analyses revealed that Nup82 interacts with the NPC components Nup136 and RAE1. The double knockout mutant nup82 nup136 showed severe growth defects, while the single knockout mutant nup82 did not, suggesting that Nup82 functions redundantly with Nup136. nup82 nup136 impaired benzothiadiazole (an analog of salicylic acid)-induced resistance to the virulent bacteria Pseudomonas syringae pv. tomato DC3000. Furthermore, transcriptome analysis of nup82 nup136 indicates that deficiency of Nup82 and Nup136 causes noticeable downregulation of immune-related genes. These results suggest that Nup82 and Nup136 are redundantly involved in transcriptional regulation of salicylic acid-responsive genes through nuclear transport of signaling molecules.

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.

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.

Identification of differentially expressed genes between developing seeds of different soybean cultivars.

Soybean is a major source of protein and oil and a primary feedstock for biodiesel production. Research on soybean seed composition and yield has revealed that protein, oil and yield are controlled quantitatively and quantitative trait loci (QTL) have been identified for each of these traits. However, very limited information is available regarding the genetic mechanisms controlling seed composition and yield. To help address this deficiency, we used Affymetrix Soybean GeneChips® to identify genes that are differentially expressed between developing seeds of the Minsoy and Archer soybean cultivars, which differ in seed weight, yield, protein content and oil content. A total of 700 probe sets were found to be expressed at significantly different (defined as having an adjusted p-value below or equal to 0.05 and an at least 2-fold difference) levels between the two cultivars at one or more of the three developmental stages and in at least one of the two years assayed. Comparison of data from soybeans collected in two different years revealed that 97 probe sets were expressed at significantly different levels in both years. Functional annotations were assigned to 78% of these 97 probe sets based on the SoyBase Affymetrix™ GeneChip® Soybean Genome Array Annotation. Genes involved in receptor binding/activity and protein binding are overrepresented among the group of 97 probe sets that were differentially expressed in both years assayed. Probe sets involved in growth/development, signal transduction, transcription, defense/stress response and protein and lipid metabolism were also identified among the 97 probe sets and their possible implications in the regulation of agronomic traits are discussed. As the Minsoy and Archer soybean cultivars differ with respect to seed size, yield, protein content and lipid content, some of the differentially expressed probe sets identified in this study may thus play important roles in controlling these traits. Others of these probe sets may be involved in regulation of general seed development or metabolism. All microarray data and expression values after GCRMA are available at the Gene Expression Omnibus (GEO) at NCBI (http://www.ncbi.nlm.nih.gov/geo), under accession number GSE21598.

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.