Expression-based network biology identifies immune-related functional modules involved in plant defense
© Tully et al.; licensee BioMed Central Ltd. 2014
Received: 13 August 2013
Accepted: 27 May 2014
Published: 3 June 2014
Plants respond to diverse environmental cues including microbial perturbations by coordinated regulation of thousands of genes. These intricate transcriptional regulatory interactions depend on the recognition of specific promoter sequences by regulatory transcription factors. The combinatorial and cooperative action of multiple transcription factors defines a regulatory network that enables plant cells to respond to distinct biological signals. The identification of immune-related modules in large-scale transcriptional regulatory networks can reveal the mechanisms by which exposure to a pathogen elicits a precise phenotypic immune response.
We have generated a large-scale immune co-expression network using a comprehensive set of Arabidopsis thaliana (hereafter Arabidopsis) transcriptomic data, which consists of a wide spectrum of immune responses to pathogens or pathogen-mimicking stimuli treatments. We employed both linear and non-linear models to generate Arabidopsis immune co-expression regulatory (AICR) network. We computed network topological properties and ascertained that this newly constructed immune network is densely connected, possesses hubs, exhibits high modularity, and displays hallmarks of a “real” biological network. We partitioned the network and identified 156 novel modules related to immune functions. Gene Ontology (GO) enrichment analyses provided insight into the key biological processes involved in determining finely tuned immune responses. We also developed novel software called OCCEAN (One Click Cis-regulatory Elements ANalysis) to discover statistically enriched promoter elements in the upstream regulatory regions of Arabidopsis at a whole genome level. We demonstrated that OCCEAN exhibits higher precision than the existing promoter element discovery tools. In light of known and newly discovered cis-regulatory elements, we evaluated biological significance of two key immune-related functional modules and proposed mechanism(s) to explain how large sets of diverse GO genes coherently function to mount effective immune responses.
We used a network-based, top-down approach to discover immune-related modules from transcriptomic data in Arabidopsis. Detailed analyses of these functional modules reveal new insight into the topological properties of immune co-expression networks and a comprehensive understanding of multifaceted plant defense responses. We present evidence that our newly developed software, OCCEAN, could become a popular tool for the Arabidopsis research community as well as potentially expand to analyze other eukaryotic genomes.
KeywordsPathogen and pathogen-mimic stimuli Linear and non-linear models Transcriptional regulatory network Network topology GO terms OCCEAN software Immune-related functional modules
Agriculture, in struggling to meet global food demands of a rapidly growing population, has been plagued with plant diseases that thwart crop production and account for a global annual average yield loss of 16 percent [1–3]. In light of the importance of agriculture to humans, it is essential to elucidate and understand the mechanisms of plant immunity at the molecular level. As with any host-pathogen conflict, plants and their disease agents are in an evolutionary arms race. When the host mounts a defense reaction, the pathogen develops new strategies to evade the defensive mechanisms, which causes the continuation of this cycle, ad infinitum[4, 5]. Our present understanding of plant immune systems reveals two primary means by which plants recognize their invaders. The recognition of Microbial-Associated Molecular Patterns (MAMPs) is the first line of defense and its associated MAMPs-Triggered Immunity (MTI) is highly efficient at repelling most pathogens. If pathogens are able to breach this first line of defense through the production of effector proteins, plants deploy Effector-Triggered Immunity (ETI) [4, 6]. Both modes of defense cause massive transcriptional reprogramming which involves complex signal transduction networks and cross-talk that is mediated by plant hormones, including: Salicylic Acid (SA), Jasmonic Acid (JA) and Ethylene (ET) [7, 8]. The Arabidopsis genome encodes approximately 1,922 transcription factors (TFs) that are implicated in diverse biological processes including the regulation of immune signaling pathways [9–11]. However, the global organization of TF-DNA interactions remains elusive. Moreover, what remains mostly unknown are the precise mechanisms by which the plant cell integrates innumerable synergistic and antagonistic immune transcriptomic signals to orchestrate fine-tuned and pathogen-specific defense responses.
In recent years, systems biology approaches, specifically network analyses that integrate experimentally derived information with computational modeling, have emerged as powerful tools for studying complex traits in diverse species [5, 12–15]. Network biology provides comprehensive analyses of the system’s components (nodes) and the relationships among them (edges). In addition, the analyses of the topological properties of the network can provide further understanding of the hierarchical organization of a complex biological system and contribute to the overall interpretation of biological complexity. In any eukaryotic cell, thousands of genes and their products orchestrate their transcriptional and translational activities to ensure the proper execution of cellular functions [13, 15, 16]. It has become evident that biological functions can be accomplished by functional modules that are embedded within the interaction networks including transcriptional gene regulatory networks [17, 18]. Comprehensive high-throughput analyses of gene expression can allow for the identification of gene clusters that are highly correlated in expression levels across multiple samples in any given cellular state [13, 16]. Generally, it is thought that genes in the same co-expression sub-network are often enriched with similar functional annotations. Additionally, the metric to measure the co-expression falls into one of two major categories: correlation coefficients or mutual information measures . Finally, finding common cis-regulatory elements (transcription factor binding sites) can aid in the identification of co-regulated gene clusters and characterization of transcriptional regulatory networks [13, 14, 17].
In the current study we employed a systems biology approach and report the construction of a large-scale Arabidopsis immune co-expression regulatory (AICR) network. We tested five diverse algorithms, i.e. PCC (Pearson Correlation Coefficient), ARACNE multiplicative (Algorithm for the Reconstruction of Accurate Cellular Networks), ARACNE additive, CLR (Context Likelihood of Relatedness), and MRNET (Minimum Redundancy NETwork) with different thresholds, which yielded 15 pairs of experimental networks along with their respective random networks [19–25]. We employed network biology analyses and determined that ARACNE multiplicative network (5,147 nodes and 38,610 edges), with threshold of 0.8 exhibits properties of a “true” network as it possesses the scale-freeness attribute (degree distribution follows a power law). Next, we partitioned the AICR network and predicted 156 functional modules containing at least six members with the largest module encompassing 178 nodes. Subsequently, we analyzed functional annotations of genes within each module and calculated enrichment of specific Gene Ontology terms to evaluate the biological significance of functional modules. To establish a causal relationship between co-expression and co-regulation, we sought to identify common cis-regulatory elements. First we developed a new, comprehensive software interface for cis-regulatory elements discovery in Arabidopsis (named OCCEAN - One Click Cis-regulatory Elements ANalysis). We demonstrated that OCCEAN boasts higher capacity (can process the entire genome scale; over 30 million characters in a single run) and features higher precision than MEME (Multiple EM for Motif Elicitation). We identified several statistically enriched novel cis-regulatory elements in our dataset. Finally, we evaluated and discussed key immune-related functional modules.
Results and discussion
Construction of the Arabidopsis immune co-expression regulatory (AICR) network
Comparison of multiple algorithms used to generate Arabidopsis Immune Co-expression (AICR) Network
PCC K = 3
PCC K = 10
PCC K = 20
PCC K = 50
PCC K = 100
PCC K = 250
Global topological properties of the AICR network
Biological implications of topological properties of the AICR network
Topological properties of the network can also be utilized to prioritize genes for further functional analysis. Highly connected and highly centered nodes are likely to play crucial roles in maintaining network integrity and controlling information flow throughout the plant immune system. Specifically, we identified 190 highly connected nodes (Hub50 encompassing at least 50 connections in the AICR). To characterize the significance of these highly connected proteins, we combined Hub50 of the AICR with biophysical interactions from plant-pathogen immune network, version 1′ (PPIN-1; Figure 1) . The PPIN-1 was constructed using effectors (virulence factors) from two pathogens, Arabidopsis immune proteins, and ~8,000 other Arabidopsis proteins covering almost one-third of the genome . Pathogens utilize effector molecules to rewire the host network in a manner conducive to pathogen proliferation and dispersal . PPIN-1 discovered 165 effector interacting proteins (effector targets) and it was also demonstrated that these effector targets are highly interconnected host proteins . Despite the unavailability of effector targets from the entire Arabidopsis genome in PPIN-1, we found that the AICR network contains 51 effector targets and three proteins (At3g63210, At4g11890 and At3g07780) are Hub50 in the AICR network (Figure 1). Given that transcriptional regulation does not necessarily coincide with biophysical interactions, this overlap is potentially significant and biologically meaningful. This also suggests that pathogen effectors target key regulatory proteins in both transcriptional as well as protein-protein interactions networks.
In addition, we selected ninety nodes (“top 90”) with the highest betweeness centrality (Figure 4). Among them, we found a number of previously characterized immune regulators or defense-related proteins, encompassing various levels of signal transduction flow: resistance proteins (RPP4, uncharacterized TIR-NBS-LRR) [42, 43], kinases mediating early signaling events (MAPKKK13, WAK) , transcriptional executors (ANAC072, WRKY54, MYB66, WER1) [45–47], secretory proteins assisting with folding and modification of newly synthesized proteins (TRAP, CNX3, Sec14p, cyclophilin, UDP-Glycosyltransferase, DnaJ family, KMS1, GlcNAc1pUT2) [48, 49], enzymes involved in production of phytohormones and antimicrobial compounds (AOC2, ILL3, PAL2) and pathogenesis-related proteins (PR1) [50, 51], hormonal response modulators (PYL6, IAA28) [52, 53] and finally, components of ubiquitin-mediated protein degradation (UBQ4, an F-box protein that interacts with SKP1) [54, 55]. Identification of these key immune regulators serves as a proof of concept for our analytical approach and justifies further analyses on suite of “top 90” proteins selected among thousands of other proteins in the AICR network.
Network clustering and module annotations
Another essential feature of a majority of biological networks is their ability to naturally organize into modules [14, 17, 18]. Discovery of such functional modules within biological networks is imperative for understanding principles of cellular organization and functions. To identify functional modules (subnetworks) that are comprised of highly interconnected sets of genes within the AICR network, we employed a fast agglomerative algorithm (FAG-EC) . This new, low time complexity algorithm operates based on a local variable, edge clustering coefficients, making it ideal to partition a relatively large network [37, 56]. Overall, we identified 156 immune-related functional modules containing at least six members with the largest module encompassing 178 nodes (Additional file 11: Table S3). Furthermore, we subjected the AICR network to a module size distribution analysis. We computed the frequency and module size and demonstrated that the distribution of the module size follows a power law distribution with an exponential truncation (Additional file 12: Figure S8), which is common for several previously described “real-world” networks [14, 18]. The module size distribution property further indicates that the AICR network possesses a modular structure.
Significantly enriched GO terms in the ten largest immune-related modules with sizes ranging from 38 to 178 genes
GO:0000166 ~ nucleotide binding
GO:0006796 ~ phosphate metabolic process
GO:0004672 ~ protein kinase activity
GO:0008219 ~ cell death
GO:0009536 ~ plastid
GO:0004672 ~ protein kinase activity
GO:0019748 ~ secondary metabolic process
GO:0009536 ~ plastid
GO:0055114 ~ oxidation reduction
GO:0005886 ~ plasma membrane
GO:0006793 ~ phosphorus metabolic process
GO:0009628 ~ response to abiotic stimulus
GO:0005618 ~ cell wall
GO:0009507 ~ chloroplast
GO:0051186 ~ cofactor metabolic process
GO:0000166 ~ nucleotide binding
GO:0009725 ~ response to hormone stimulus
GO:0006793 ~ phosphorus metabolic process
GO:0031224 ~ intrinsic to membrane
GO:0006793 ~ phosphorus metabolic process
GO:0032555 ~ purine ribonucleotide binding
GO:0005886 ~ plasma membrane
Development of OCCEAN software and identification of cis-regulatory elements
Inferring biological relevance of immune-related modules
Diverse kinases account for the major functional category of Module 1, encompassing 46 out of 178 proteins representing receptor-like kinases (RLKs), mitogen-activated protein kinases (MAPKs), calcium-dependent protein kinases (CDPKs), wall-associated kinases (WAKs), etc. Typically, RLKs perceive MAMPs such as flg22, elf18, chitin and OGs and trigger MTI . Subsequently, a MAPK cascade is activated and downstream signaling ensues including phytohormonal crosstalk [8, 27, 44, 67]. These events result in the activation of a wide-spectrum of immune responses, such as induced biosynthesis of phytoalexins and other antimicrobial secondary metabolites such as glucosinolates, calcium influx, as well as production and accumulation of reactive oxygen species (ROS) [27, 68].
At the second, more powerful layer of the immune response, plants usually deploy NLR receptors (nucleotide-binding domain and leucine-rich repeat receptors; a major class of R proteins). The NLR receptors directly recognize specific effectors or indirectly detect effector activities and trigger ETI [1, 2, 6, 7]. We discovered 8 NLR genes in the kinase sub-cluster of Module 1. These include proteins conferring resistance to diverse pathogens such as Pseudomonas syringae, Albugo candida, Hyaloperonospora arabidopsidis and various Fusarium races.
The second most abundant GO category in Module 1 was immune-related genes (count: 32). This functional sub-cluster contains various defense-related transcription factors such as WRKY15, WRKY33, MYB113 and ANAC061 [9, 47], four members of the RING-H2 finger protein family, a heat shock transcription factor HSF A-4a , an ethylene-responsive transcription factor ERF105 and, intriguingly, three scarecrow-like GRAS family transcription factor genes, known to be primarily implicated in plant root development [69–71]. Collectively, transcription factor genes account for over 40% of this sub-cluster’s gene count. In addition, we identified a small group of genes strongly associated with response to fungal cell wall elicitor chitin, such as two auxin-regulated SAUR genes, one lectin-like, Enhanced Disease Susceptibility 1 (EDS1) and senescence-associated protein SAG102 [72–76].
Another enriched GO category in Module 1 was membrane-associated proteins including both plasma membrane- and organellar membrane-resident peptides; not surprisingly, the largest class of these membrane proteins is transporters. In plants, uptake and translocation of nutrients play essential roles in physiological processes including plant growth, nutrition, signal transduction, and development [77–79]. Transport processes are also critical for the reallocation of resources and it was previously shown that defense signals cause realignment of transport activities to redirect resources toward immune responses and protect tissues of high value . Of 31 genes assigned to this sub-cluster, seven were transport proteins (22.6% of this sub-cluster’s gene count). The prevailing class was sugar (galactose and sucrose) transporters, in agreement with the fact that this class of proteins constitutes a key component for carbon partitioning at the whole plant level and is involved in both symbiotic and parasitic plant-fungi interactions . Additionally, we identified CNI1 (Carbon-Nitrogen Insensitive 1) that is a key regulator of the Carbon/Nitrogen response for growth phase transition in Arabidopsis seedlings , as well as ATP Binding Cassette (ABC) and inorganic phosphate transporters that were previously postulated to be required for organ growth, nutrition, development, and stress responses . Last but not least, we also identified two other well-known immune regulators: Syntaxin 121/PEN1 and MLO2 genes, required for resistance against barley powdery mildew, Blumeria graminis sp. hordei and a fungal pathogen Colletotrichum higginsianum[84, 85].
Another essential aspect of large-scale transcriptomic analyses is to link co-expression with co-regulation, i.e., to determine presence of common cis-regulatory elements among genes of the same module. Several known cis-regulatory elements were discovered in the promoters of Module 1 genes (Additional file 14: Table S5). Two of the enriched elements, TCATGG and CATGGA, overlap with the octadecameric CArG box sequence 5′-CTTACCTTTCATGGATTA-3′, identified in the APETALA3 promoter . APETALA3, a class B homeotic organ identity gene, was originally discovered as the central regulator of petal development in Arabidopsis flowers, but later also shown to be involved in control of the floral meristem proliferation, regulation of flowering time, and other plant reproductive processes [87–89]. Inhibition of reproductive processes is an expected outcome of an immune stimulation, given the recent report describing a transition from growth to defense following immune stimuli treatments in Arabidopsis . Involvement of a developmental regulator in defense responses is not unprecedented and was previously shown for a MYB-related gene ASYMMETRIC LEAVES 1 (AS1) .
Jasmonic Acid (JA) and Ethylene (ET) are two plant hormones antagonistic to SA  and down-regulation of major JA/ET signaling proteins will promote SA signaling by suppressing the JA/ET pathway . Consistent with this observation, we identified an enriched element, ATCTTG that resembles the binding site for Ethylene-Insensitive 3 (EIN3) and three EIN3-LIKE (EIL) proteins  as well as two hexamers GTCGTC and CGTCGT, which overlap with the core binding site of JASE2 motif in two JA biosynthetic genes OPR1/OPR2. Recently, EIN3 and EIL1 were shown to negatively regulate MAMP-triggered immunity via a direct binding and down-regulation of the SA biosynthetic gene SID2. Given that OPRs are implicated in JA biosynthesis, we expect a fine-tuned interplay between SA and JA/ET signaling pathways . In addition, several novel cis-regulatory elements were discovered in our study (Additional file 14: Table S5), and the identity of the cognate transcription factors and their functional relevance will be subject to future studies. Very recently, it was demonstrated that numerous transcription factors can recognize secondary elements, which in some cases are completely sequence-unrelated to the primary element . In conjunction with this study, our analysis of Module 1 shows that ARR14 (a MYB family transcription factor) not only recognizes the known primary cis-regulatory element AGATA/TCG but also can specifically bind to a previously unknown cis-regulatory element AGATCT. Thus, it’s likely that transcription factor(s) have extended list of binding specificities beyond the currently known cis-regulatory elements. We propose that the above described transcription factors and additional unknown regulatory proteins coordinate the gene regulation in this module.
Our data provide further insights into the transcriptional regulatory mechanisms repressing additional putative negative regulators of plant defense upon treatments with pathogen-mimicking stimuli.
In this study, we used a systems-level network biology approach to construct genome-wide Arabidopsis immune co-expression network and demonstrated that this network shares properties of a ‘real-world network’. Topological properties-based partitioning allowed us to unravel 156 distinct immune-related functional modules. We demonstrated that nodes in the same module are not only highly correlated at the expression level but also densely connected to each other. Detailed analyses of two key immune-related modules provided a systems-level understanding of how plant cells coordinate distinct immune signals to orchestrate fine-tuned and pathogen-specific defense responses. Our novel approach to discover cis-regulatory elements using OCCEAN is an effective method of solving the issue of finding novel motifs within a sequence set. OCCEAN has advantages over other programs of the same purpose, such as APPLES (Analysis of Plant Promoter-Linked Elements) and MEME [63, 97]. APPLES requires finding organisms of a specified relational distance for comparison, which can burden the user, and MEME has the statistical risk of discarding actual existing motifs. These problems are avoided in our solution whilst continuing to maintain a fair amount of focus on client-side simplicity. In addition, OCCEAN has the capacity to be expanded for analyses of other eukaryotes genomes, such as fly and human. Our systems-level approach to examine cis-regulatory elements (the putative transcription factor binding sites) in the promoters of the co-expressed genes made it possible to link co-expression to co-regulation of genes in the same module.
Data download, selection criteria for differentially expressed (DE) genes dataset and promoter sequence acquisition
We utilized available transcriptomic data of transcriptional responses extracted from 271 microarray experiments representing nine major immune-related studies (Additional file 1: Table S1 and Additional file 2: Supporting methods) [26–34]. Priority was given to well-referenced studies, employing the Affymetrix ATH1 GeneChip array, encompassing the broadest possible spectrum of plant defense responses upon pathogen infections or pathogen-mimicking stimuli treatment (Additional file 2: Supporting methods). Lists of probes showing significant up- or down-regulation in each experimental condition were compiled, using criteria for significance employed in the respective original study (Additional file 3: Table S2). This led us to the identification of a list 8,377 genes differentially expressed between all treatments . For each of these genes, we acquired 1000 bp upstream of the transcriptional start site from TAIR Version 10 at http://www.arabidopsis.org. These upstream regions were searched for putative transcription factor binding sites.
Network construction, topological properties, network clustering and Gene Ontology (GO)
The microarray data presented in Additional file 3: Table S2 was used to construct a gene co-expression network using both linear and non-linear models. In the linear model, Pearson Correlation Coefficients (PCC) was measured based on the mutual k-nearest neighbor method of Ruan et al. [18, 97] with some modifications. In contrast to Ruan et al. [18, 97], the k-list for each gene included k-1 other genes, and at least one (and possibly more) gene sharing the smallest of the k PCCs. For example, for a k of 10, the k-list might include 9 genes in order of decreasing PCC and one or more genes having the next lower PCC value. This was done in consideration of the fact that all genes of equal PCC are equally valid as connected nodes. Random networks were constructed by randomly assigning gene expression values (0 = no change, +1 = up-regulated, -1 = down-regulated) to the genes from the immune transcriptional profiling experiments. The network based on PCC-threshold was inferred by first calculating the PCC for every gene-gene pair in the dataset yielding 8,377 × 8,377 PCC values. These values were tested against 0.9 threshold and only the PCC values that were greater than or equal to 0.9 were selected to indicate that there is an edge between this gene pairs. Other values were discarded. R was used to calculate the non-linear mutual information (MI) . First the MI was calculated using parmigene package were the MI for every gene pair was calculated resulting in 8377 × 8377 MI values . Then these values where used by the ARACNE multiplicative (Algorithm for the Reconstruction of Accurate Cellular Networks), ARACNE additive, CLR (Context Likelihood of Relatedness), and MRNET (Minimum Redundancy NETwork) algorithms, using the same package, to infer a weighted adjacency matrix [21–25]. 80% and 90% of the maximum MI values in each weight adjacency matrix were chosen as the threshold to determine the existence of an edge between each pair of genes. If the value is greater than or equal to the threshold, this indicates an edge between a given pair of genes.
A modified and parallelized Cytoscape version 2.8.4 [37, 38] and the Clusterviz version 1.2 (FAG-EC algorithm) [37, 38, 56, 100] and ClusterMaker version 1.10 [100, 101] plugins were used to visualize this network, calculate its parameters and network topological properties. Network partitioning was performed using FAG-EC algorithm [56, 100] and the Markov Clustering Algorithm . For each gene in the network sub-clusters, a Gene Ontology [57, 58] molecular function and biological process was assigned from the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 .
BLAST-based generation of common cis-regulatory elements Sequences obtained from TAIR  (1000 bp 5′ of the transcription start site) were incorporated into a local BLAST [103, 104] database using the formatdb utility that is included in NCBI’s package download of BLAST. Each sequence in the local BLAST database was used in turn as a query sequence to find subsequences shared in the promoter regions of a number of genes. The resulting common subsequences were imported into a relational database. Similar sequences with at least 75% identity and an E-value ≤ 1× 10-4 were scanned using a sliding window search in order to extract all contained 6-mers. All 6-mers having no more than four consecutive single nucleotide repeats were deemed putative transcription factor binding sites and retained for further computational analysis. This process has been incorporated into our web-accessible tool, OCCEAN.
Statistical enrichment of cis-regulatory elements using modified version of POBO
Putative transcription factor binding sites results from BLAST were analyzed for statistically significant over representation/under representation using the POBO (a promoter bootstrapping program) tool . POBO is an exceptional tool for statistical sequence significance analysis, but is inhibited by the constraint that it can perform analysis only on one cis-regulatory element at a time. As our sequence sets were high in number, manual entry would be infeasible. Therefore, a wrapper program was written to take each individual sequence, create a cluster file consisting of the set of promoter identifiers/sequences from the original experimentally found gene promoter dataset that contain the sequence, and run an instance of POBO with the sequence and cluster file as input. The program then extracted the desired information from the resulting HTML files from all of the individual runs and put them in a single spreadsheet for analysis. Running POBO locally required downloading the source code from http://ekhidna.biocenter.helsinki.fi/poxo/pobo/ and setting up a local MySQL database with all known Arabidopsis gene promoters corresponding to ~33,000 genes as the background data. 1 kb promoter sequence for all ~33,000 genes were obtained from the TAIR10 genome release, and used as background for the POBO analysis. The POBO results were converted to spreadsheet format for further analysis. True positives (nTPs), false positives (nFPs) and precision for MEME as well as OCCEAN were calculated as described in .
Development of One-Click Cis-Element ANalysis (OCCEAN) software for genome-wide identification of cis-regulatory elements
mutual k-nearest neighbor
One click cis-regulatory elements ANalysis
Microbial-associated molecular patterns
Pearson correlation coefficient
Fast agglomerate algorithm- edge clustering
Multiple EM for motif elicitation
Arabidopsis thaliana map
A database of plant cis-acting regulatory DNA elements
The arabidopsis gene regulatory information server
A promoter bootstrapping program
Analysis of plant promoter-linked elements, nTPs, True positives
Algorithm for the reconstruction of accurate cellular networks)
Context likelihood of relatedness
Minimum redundancy NETwork
Arabidopsis immune co-expression regulatory.
The authors are greatly thankful to Dr. Purushotham Bangalore and Mr. Larry Owen for technical support. The authors thank Drs. Karolina Mukhtar and Michelle Amaral, Ms. Cassandra Garbutt and Ms. Jessica Lopez for critical reading of the manuscript and valuable suggestions.
- Mukhtar MS: Engineering NLR immune receptors for broad-spectrum disease resistance. Trends Plant Sci. 2013, 18 (9): 469-472. 10.1016/j.tplants.2013.08.005.PubMed
- Pajerowska-Mukhtar KM, Emerine DK, Mukhtar MS: Tell me more: roles of NPRs in plant immunity. Trends Plant Sci. 2013, 18 (7): 402-411. 10.1016/j.tplants.2013.04.004.PubMed
- Wulff BB, Horvath DM, Ward ER: Improving immunity in crops: new tactics in an old game. Curr Opin Plant Biol. 2011, 14 (4): 468-476. 10.1016/j.pbi.2011.04.002.PubMed
- Jones JD, Dangl JL: The plant immune system. Nature. 2006, 444 (7117): 323-329. 10.1038/nature05286.PubMed
- Mukhtar MS, Carvunis AR, Dreze M, Epple P, Steinbrenner J, Moore J, Tasan M, Galli M, Hao T, Nishimura MT, Pevzner SJ, Donovan SE, Ghamsari L, Santhanam B, Romero V, Poulin MM, Gebreab F, Gutierrez BJ, Tam S, Monachello D, Boxem M, Harbort CJ, McDonald N, Gai L, Chen H, He Y, Consortium EUF, Vandenhaute J, Roth FP, Hill DE, et al: Independently evolved virulence effectors converge onto hubs in a plant immune system network. Science. 2011, 333 (6042): 596-601. 10.1126/science.1203659.PubMed CentralPubMed
- Nishimura MT, Dangl JL: Arabidopsis and the plant immune system. Plant J Cell Mole Biol. 2010, 61 (6): 1053-1066. 10.1111/j.1365-313X.2010.04131.x.
- Mukhtar MS, Nishimura MT, Dangl J: NPR1 in plant defense: it’s not over ‘til it’s turned over. Cell. 2009, 137 (5): 804-806. 10.1016/j.cell.2009.05.010.PubMed
- Spoel SH, Dong X: How do plants achieve immunity? Defence without specialized immune cells. Nat Rev Immunol. 2012, 12 (2): 89-100. 10.1038/nri3141.PubMed
- Eulgem T, Somssich IE: Networks of WRKY transcription factors in defense signaling. Curr Opin Plant Biol. 2007, 10 (4): 366-371. 10.1016/j.pbi.2007.04.020.PubMed
- Guo A, He K, Liu D, Bai S, Gu X, Wei L, Luo J: DATF: a database of Arabidopsis transcription factors. Bioinformatics. 2005, 21 (10): 2568-2569. 10.1093/bioinformatics/bti334.PubMed
- Yilmaz A, Mejia-Guerra MK, Kurz K, Liang X, Welch L, Grotewold E: AGRIS: the arabidopsis gene regulatory information server, an update. Nucleic Acids Res. 2011, 39 (Database issue): D1118-D1122.PubMed CentralPubMed
- Arabidopsis Interactome Mapping C: Evidence for network evolution in an Arabidopsis interactome map. Science. 2011, 333 (6042): 601-607.
- Carrera J, Rodrigo G, Jaramillo A, Elena SF: Reverse-engineering the Arabidopsis thaliana transcriptional network under changing environmental conditions. Genome Biol. 2009, 10 (9): R96-10.1186/gb-2009-10-9-r96.PubMed CentralPubMed
- Mao L, Van Hemert JL, Dash S, Dickerson JA: Arabidopsis gene co-expression network and its functional modules. BMC Bioinformatics. 2009, 10: 346-10.1186/1471-2105-10-346.PubMed CentralPubMed
- Vidal M, Cusick ME, Barabasi AL: Interactome networks and human disease. Cell. 2011, 144 (6): 986-998. 10.1016/j.cell.2011.02.016.PubMed CentralPubMed
- Ma S, Bohnert HJ: Integration of Arabidopsis thaliana stress-related transcript profiles, promoter structures, and cell-specific expression. Genome Biol. 2007, 8 (4): R49-10.1186/gb-2007-8-4-r49.PubMed CentralPubMed
- Heyndrickx KS, Vandepoele K: Systematic identification of functional plant modules through the integration of complementary data sources. Plant Physiol. 2012, 159 (3): 884-901. 10.1104/pp.112.196725.PubMed CentralPubMed
- Ruan J, Perez J, Hernandez B, Lei C, Sunter G, Sponsel VM: Systematic identification of functional modules and cis-regulatory elements in Arabidopsis thaliana. BMC Bioinformatics. 2011, 12 (Suppl 12): S2-10.1186/1471-2105-12-S12-S2.PubMed CentralPubMed
- Song L, Langfelder P, Horvath S: Comparison of co-expression measures: mutual information, correlation, and model based indices. BMC Bioinformatics. 2012, 13: 328-10.1186/1471-2105-13-328.PubMed CentralPubMed
- Ruan J, Zhang W: Identifying network communities with a high resolution. Physical Rev E Stat Nonlinear Soft Matter Physics. 2008, 77 (1 Pt 2): 016104-
- Margolin AA, Nemenman I, Basso K, Wiggins C, Stolovitzky G, Dalla Favera R, Califano A: ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics. 2006, 7 (Suppl 1): S7-10.1186/1471-2105-7-S1-S7.PubMed CentralPubMed
- Zoppoli P, Morganella S, Ceccarelli M: TimeDelay-ARACNE: Reverse engineering of gene networks from time-course data by an information theoretic approach. BMC Bioinformatics. 2010, 11: 154-10.1186/1471-2105-11-154.PubMed CentralPubMed
- Olsen C, Meyer PE, Bontempi G: On the impact of entropy estimation on transcriptional regulatory network inference based on mutual information. EURASIP J Bioinformatics Systems Biol. 2009, 2009 (1): 308959-
- Meyer PE, Kontos K, Lafitte F, Bontempi G: Information-theoretic inference of large transcriptional regulatory networks. EURASIP J Bioinformatics Systems Biol. 2007, 2007: 79879-
- Sales G, Romualdi C: Parmigene–a parallel R package for mutual information estimation and gene network reconstruction. Bioinformatics. 2011, 27 (13): 1876-1877. 10.1093/bioinformatics/btr274.PubMed
- Chandran D, Inada N, Hather G, Kleindt CK, Wildermuth MC: Laser microdissection of Arabidopsis cells at the powdery mildew infection site reveals site-specific processes and regulators. Proc Nat Acad Sci USA. 2010, 107 (1): 460-465. 10.1073/pnas.0912492107.PubMed CentralPubMed
- Denoux C, Galletti R, Mammarella N, Gopalan S, Werck D, De Lorenzo G, Ferrari S, Ausubel FM, Dewdney J: Activation of defense response pathways by OGs and Flg22 elicitors in Arabidopsis seedlings. Mole Plant. 2008, 1 (3): 423-445. 10.1093/mp/ssn019.
- Eulgem T, Weigman VJ, Chang HS, McDowell JM, Holub EB, Glazebrook J, Zhu T, Dangl JL: Gene expression signatures from three genetically separable resistance gene signaling pathways for downy mildew resistance. Plant Physiol. 2004, 135 (2): 1129-1144. 10.1104/pp.104.040444.PubMed CentralPubMed
- Goda H, Sasaki E, Akiyama K, Maruyama-Nakashita A, Nakabayashi K, Li W, Ogawa M, Yamauchi Y, Preston J, Aoki K, Kiba T, Takatsuto S, Fujioka S, Asami T, Nakano T, Kato H, Mizuno T, Sakakibara H, Yamaguchi S, Nambara E, Kamiya Y, Takahashi H, Hirai MY, Sakurai T, Shinozaki K, Saito K, Yoshida S, Shimada Y: The AtGenExpress hormone and chemical treatment data set: experimental design, data evaluation, model data analysis and data access. Plant J Cell Mole Biol. 2008, 55 (3): 526-542. 10.1111/j.1365-313X.2008.03510.x.
- Ramonell K, Berrocal-Lobo M, Koh S, Wan J, Edwards H, Stacey G, Somerville S: Loss-of-function mutations in chitin responsive genes show increased susceptibility to the powdery mildew pathogen Erysiphe cichoracearum. Plant Physiol. 2005, 138 (2): 1027-1036. 10.1104/pp.105.060947.PubMed CentralPubMed
- Thilmony R, Underwood W, He SY: Genome-wide transcriptional analysis of the Arabidopsis thaliana interaction with the plant pathogen Pseudomonas syringae pv. tomato DC3000 and the human pathogen Escherichia coli O157:H7. Plant J Cell Mole Biol. 2006, 46 (1): 34-53. 10.1111/j.1365-313X.2006.02725.x.
- Truman W, de Zabala MT, Grant M: Type III effectors orchestrate a complex interplay between transcriptional networks to modify basal defence responses during pathogenesis and resistance. Plant J Cell Mole Biol. 2006, 46 (1): 14-33. 10.1111/j.1365-313X.2006.02672.x.
- Wang D, Amornsiripanitch N, Dong X: A genomic approach to identify regulatory nodes in the transcriptional network of systemic acquired resistance in plants. PLoS Pathogens. 2006, 2 (11): e123-10.1371/journal.ppat.0020123.PubMed CentralPubMed
- Zipfel C, Kunze G, Chinchilla D, Caniard A, Jones JD, Boller T, Felix G: Perception of the bacterial PAMP EF-Tu by the receptor EFR restricts Agrobacterium-mediated transformation. Cell. 2006, 125 (4): 749-760. 10.1016/j.cell.2006.03.037.PubMed
- Altay G, Emmert-Streib F: Inferring the conservative causal core of gene regulatory networks. BMC Systems Biol. 2010, 4: 132-10.1186/1752-0509-4-132.
- Luo F, Yang Y, Zhong J, Gao H, Khan L, Thompson DK, Zhou J: Constructing gene co-expression networks and predicting functions of unknown genes by random matrix theory. BMC Bioinformatics. 2007, 8: 299-10.1186/1471-2105-8-299.PubMed CentralPubMed
- Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, Christmas R, Avila-Campilo I, Creech M, Gross B, Hanspers K, Isserlin R, Kelley R, Killcoyne S, Lotia S, Maere S, Morris J, Ono K, Pavlovic V, Pico AR, Vailaya A, Wang PL, Adler A, Conklin BR, Hood L, Kuiper M, Sander C, Schmulevich I, Schwikowski B, Warner GJ, et al: Integration of biological networks and gene expression data using Cytoscape. Nat. Protoc. 2007, 2 (10): 2366-2382. 10.1038/nprot.2007.324.PubMed CentralPubMed
- Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T: Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011, 27 (3): 431-432. 10.1093/bioinformatics/btq675.PubMed CentralPubMed
- Seebacher J, Gavin AC: SnapShot: Protein-protein interaction networks. Cell. 2011, 144 (6): 1000-10.1016/j.cell.2011.02.025. 1000 e1001PubMed
- Batada NN, Reguly T, Breitkreutz A, Boucher L, Breitkreutz BJ, Hurst LD, Tyers M: Stratus not altocumulus: a new view of the yeast protein interaction network. PLoS Biol. 2006, 4 (10): e317-10.1371/journal.pbio.0040317.PubMed CentralPubMed
- Dangl JL, Jones JD: Plant pathogens and integrated defence responses to infection. Nature. 2001, 411 (6839): 826-833. 10.1038/35081161.PubMed
- van der Biezen EA, Freddie CT, Kahn K, Parker JE, Jones JD: Arabidopsis RPP4 is a member of the RPP5 multigene family of TIR-NB-LRR genes and confers downy mildew resistance through multiple signalling components. Plant J Cell Mole Biol. 2002, 29 (4): 439-451. 10.1046/j.0960-7412.2001.01229.x.
- Knoth C, Ringler J, Dangl JL, Eulgem T: Arabidopsis WRKY70 is required for full RPP4-mediated disease resistance and basal defense against Hyaloperonospora parasitica. MPMI. 2007, 20 (2): 120-128. 10.1094/MPMI-20-2-0120.PubMed
- Rasmussen MW, Roux M, Petersen M, Mundy J: MAP kinase cascades in arabidopsis innate immunity. Front Plant Sci. 2012, 3: 169-PubMed CentralPubMed
- Li J, Besseau S, Toronen P, Sipari N, Kollist H, Holm L, Palva ET: Defense-related transcription factors WRKY70 and WRKY54 modulate osmotic stress tolerance by regulating stomatal aperture in Arabidopsis. New Phytol. 2013, 200 (2): 457-472. 10.1111/nph.12378.PubMed CentralPubMed
- Besseau S, Li J, Palva ET: WRKY54 and WRKY70 co-operate as negative regulators of leaf senescence in Arabidopsis thaliana. J Exp Bot. 2012, 63 (7): 2667-2679. 10.1093/jxb/err450.PubMed CentralPubMed
- Jensen MK, Hagedorn PH, de Torres-Zabala M, Grant MR, Rung JH, Collinge DB, Lyngkjaer MF: Transcriptional regulation by an NAC (NAM-ATAF1,2-CUC2) transcription factor attenuates ABA signalling for efficient basal defence towards Blumeria graminis f. sp. hordei in Arabidopsis. Plant J Cell Mole Biol. 2008, 56 (6): 867-880. 10.1111/j.1365-313X.2008.03646.x.
- Bassham DC, Brandizzi F, Otegui MS, Sanderfoot AA: The secretory system of Arabidopsis. Arabidopsis Book/Am Soc Plant Biol. 2008, 6: e0116-
- Mishiba K, Nagashima Y, Suzuki E, Hayashi N, Ogata Y, Shimada Y, Koizumi N: Defects in IRE1 enhance cell death and fail to degrade mRNAs encoding secretory pathway proteins in the Arabidopsis unfolded protein response. Proc Nat Acad Sci USA. 2013, 110 (14): 5713-5718. 10.1073/pnas.1219047110.PubMed CentralPubMed
- Tierens KF, Thomma BP, Brouwer M, Schmidt J, Kistner K, Porzel A, Mauch-Mani B, Cammue BP, Broekaert WF: Study of the role of antimicrobial glucosinolate-derived isothiocyanates in resistance of Arabidopsis to microbial pathogens. Plant Physiol. 2001, 125 (4): 1688-1699. 10.1104/pp.125.4.1688.PubMed CentralPubMed
- Damon C, Dmitrieva J, Muhovski Y, Francis F, Lins L, Ledoux Q, Luwaert W, Marko IE, Mauro S, Ongena M, Thonart P, Veys P, Portetelle D, Twizere JC, Vandenbol M: Interaction network of antimicrobial peptides of Arabidopsis thaliana, based on high-throughput yeast two-hybrid screening. Plant Physiol Biochem PPB/Societe francaise de physiologie vegetale. 2012, 58: 245-252.
- Boatwright JL, Pajerowska-Mukhtar K: Salicylic acid: an old hormone up to new tricks. Mole Plant Pathol. 2013, 14 (6): 623-634. 10.1111/mpp.12035.
- Pajerowska-Mukhtar KM, Wang W, Tada Y, Oka N, Tucker CL, Fonseca JP, Dong X: The HSF-like transcription factor TBF1 is a major molecular switch for plant growth-to-defense transition. CB. 2012, 22 (2): 103-112.PubMed CentralPubMed
- Goritschnig S, Zhang Y, Li X: The ubiquitin pathway is required for innate immunity in Arabidopsis. Plant J Cell Mole Biol. 2007, 49 (3): 540-551. 10.1111/j.1365-313X.2006.02978.x.
- Trujillo M, Ichimura K, Casais C, Shirasu K: Negative regulation of PAMP-triggered immunity by an E3 ubiquitin ligase triplet in Arabidopsis. CB. 2008, 18 (18): 1396-1401.PubMed
- Wang J, Li M, Deng Y, Pan Y: Recent advances in clustering methods for protein interaction networks. BMC Genomics. 2010, 11 (Suppl 3): S10-10.1186/1471-2164-11-S3-S10.
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralPubMed
- Gaudet P, Livstone MS, Lewis SE, Thomas PD: Phylogenetic-based propagation of functional annotations within the Gene Ontology consortium. Brief Bioinform. 2011, 12 (5): 449-462. 10.1093/bib/bbr042.PubMed CentralPubMed
- Jiao X, Sherman BT, da Huang W, Stephens R, Baseler MW, Lane HC, Lempicki RA: DAVID-WS: a stateful web service to facilitate gene/protein list analysis. Bioinformatics. 2012, 28 (13): 1805-1806. 10.1093/bioinformatics/bts251.PubMed CentralPubMed
- Higo K, Ugawa Y, Iwamoto M, Korenaga T: Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999, 27 (1): 297-300. 10.1093/nar/27.1.297.PubMed CentralPubMed
- O’Connor TR, Dyreson C, Wyrick JJ: Athena: a resource for rapid visualization and systematic analysis of Arabidopsis promoter sequences. Bioinformatics. 2005, 21 (24): 4411-4413. 10.1093/bioinformatics/bti714.PubMed
- Palaniswamy SK, James S, Sun H, Lamb RS, Davuluri RV, Grotewold E: AGRIS and AtRegNet. a platform to link cis-regulatory elements and transcription factors into regulatory networks. Plant Physiol. 2006, 140 (3): 818-829. 10.1104/pp.105.072280.PubMed CentralPubMed
- Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS: MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009, 37 (Web Server issue): W202-W208.PubMed CentralPubMed
- Kankainen M, Holm L: POBO, transcription factor binding site verification with bootstrapping. Nucleic Acids Res. 2004, 32 (Web Server issue): W222-W229.PubMed CentralPubMed
- Bulow L, Brill Y, Hehl R: AthaMap-assisted transcription factor target gene identification in Arabidopsis thaliana. Database J Biol Databases Curation. 2010, 2010: baq034-
- Newman MA, Sundelin T, Nielsen JT, Erbs G: MAMP (microbe-associated molecular pattern) triggered immunity in plants. Front Plant Sci. 2013, 4: 139-PubMed CentralPubMed
- Spoel SH, Dong X: Making sense of hormone crosstalk during plant immune responses. Cell Host Microbe. 2008, 3 (6): 348-351. 10.1016/j.chom.2008.05.009.PubMed
- Mersmann S, Bourdais G, Rietz S, Robatzek S: Ethylene signaling regulates accumulation of the FLS2 receptor and is required for the oxidative burst contributing to plant immunity. Plant Physiol. 2010, 154 (1): 391-400. 10.1104/pp.110.154567.PubMed CentralPubMed
- Meng X, Xu J, He Y, Yang KY, Mordorski B, Liu Y, Zhang S: Phosphorylation of an ERF transcription factor by Arabidopsis MPK3/MPK6 regulates plant defense gene induction and fungal resistance. Plant Cell. 2013, 25 (3): 1126-1142. 10.1105/tpc.112.109074.PubMed CentralPubMed
- Nakano T, Suzuki K, Fujimura T, Shinshi H: Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006, 140 (2): 411-432. 10.1104/pp.105.073783.PubMed CentralPubMed
- Gobbato E, Marsh JF, Vernie T, Wang E, Maillet F, Kim J, Miller JB, Sun J, Bano SA, Ratet P, Mysore KS, Dénarié J, Schultze M, Oldroyd GE: A GRAS-type transcription factor with a specific function in mycorrhizal signaling. CB. 2012, 22 (23): 2236-2241.PubMed
- Schon M, Toller A, Diezel C, Roth C, Westphal L, Wiermer M, Somssich IE: Analyses of wrky18 wrky40 plants reveal critical roles of SA/EDS1 signaling and indole-glucosinolate biosynthesis for Golovinomyces orontii resistance and a loss-of resistance towards Pseudomonas syringae pv. tomato AvrRPS4. MPMI. 2013, 26 (7): 758-767. 10.1094/MPMI-11-12-0265-R.PubMed
- Wiermer M, Feys BJ, Parker JE: Plant immunity: the EDS1 regulatory node. Curr Opin Plant Biol. 2005, 8 (4): 383-389. 10.1016/j.pbi.2005.05.010.PubMed
- Anai T, Kono N, Kosemura S, Yamamura S, Hasegawa K: Isolation and characterization of an auxin-inducible SAUR gene from radish seedlings. DNA Seq J DNA Sequencing Mapping. 1998, 9 (5–6): 329-333.
- Leyser HM, Pickett FB, Dharmasiri S, Estelle M: Mutations in the AXR3 gene of Arabidopsis result in altered auxin response including ectopic expression from the SAUR-AC1 promoter. Plant J Cell Mole Biol. 1996, 10 (3): 403-413. 10.1046/j.1365-313x.1996.10030403.x.
- Li Y, Liu ZB, Shi X, Hagen G, Guilfoyle TJ: An auxin-inducible element in soybean SAUR promoters. Plant Physiol. 1994, 106 (1): 37-43. 10.1104/pp.106.1.37.PubMed CentralPubMed
- Tommasini R, Vogt E, Fromenteau M, Hortensteiner S, Matile P, Amrhein N, Martinoia E: An ABC-transporter of Arabidopsis thaliana has both glutathione-conjugate and chlorophyll catabolite transport activity. Plant J Cell Mole Biol. 1998, 13 (6): 773-780. 10.1046/j.1365-313X.1998.00076.x.
- Guo WJ, Nagy R, Chen HY, Pfrunder S, Yu YC, Santelia D, Frommer WB, Martinoia E: SWEET17, a facilitative transporter, mediates fructose transport across the tonoplast of Arabidopsis roots and leaves. Plant Physiol. 2014, 164 (2): 777-789. 10.1104/pp.113.232751.PubMed CentralPubMed
- Shkolnik-Inbar D, Adler G, Bar-Zvi D: ABI4 downregulates expression of the sodium transporter HKT1;1 in Arabidopsis roots and affects salt tolerance. Plant J Cell Mole Biol. 2013, 73 (6): 993-1005. 10.1111/tpj.12091.
- Nour-Eldin HH, Andersen TG, Burow M, Madsen SR, Jorgensen ME, Olsen CE, Dreyer I, Hedrich R, Geiger D, Halkier BA: NRT/PTR transporters are essential for translocation of glucosinolate defence compounds to seeds. Nature. 2012, 488 (7412): 531-534. 10.1038/nature11285.PubMed
- Doidy J, Grace E, Kuhn C, Simon-Plas F, Casieri L, Wipf D: Sugar transporters in plants and in their interactions with fungi. Trends Plant Sci. 2012, 17 (7): 413-422. 10.1016/j.tplants.2012.03.009.PubMed
- Sato T, Maekawa S, Yasuda S, Sonoda Y, Katoh E, Ichikawa T, Nakazawa M, Seki M, Shinozaki K, Matsui M, Goto DB, Ikeda A, Yamaguchi J: CNI1/ATL31, a RING-type ubiquitin ligase that functions in the carbon/nitrogen response for growth phase transition in Arabidopsis seedlings. Plant J Cell Mole Biol. 2009, 60 (5): 852-864. 10.1111/j.1365-313X.2009.04006.x.
- Kang J, Park J, Choi H, Burla B, Kretzschmar T, Lee Y, Martinoia E: Plant ABC Transporters. Arabidopsis book / Am Soc Plant Biol. 2011, 9: e0153-
- Consonni C, Bednarek P, Humphry M, Francocci F, Ferrari S, Harzen A, Ver Loren van Themaat E, Panstruga R: Tryptophan-derived metabolites are required for antifungal defense in the Arabidopsis mlo2 mutant. Plant Physiol. 2010, 152 (3): 1544-1561. 10.1104/pp.109.147660.PubMed CentralPubMed
- Lorek J, Griebel T, Jones AM, Kuhn H, Panstruga R: The role of Arabidopsis heterotrimeric G-protein subunits in MLO2 function and MAMP-triggered immunity. MPMI. 2013, 26 (9): 991-1003. 10.1094/MPMI-03-13-0077-R.PubMed
- Tilly JJ, Allen DW, Jack T: The CArG boxes in the promoter of the Arabidopsis floral organ identity gene APETALA3 mediate diverse regulatory effects. Development. 1998, 125 (9): 1647-1657.PubMed
- Mara CD, Huang T, Irish VF: The Arabidopsis floral homeotic proteins APETALA3 and PISTILLATA negatively regulate the BANQUO genes implicated in light signaling. Plant Cell. 2010, 22 (3): 690-702. 10.1105/tpc.109.065946.PubMed CentralPubMed
- Krizek BA, Meyerowitz EM: The Arabidopsis homeotic genes APETALA3 and PISTILLATA are sufficient to provide the B class organ identity function. Development. 1996, 122 (1): 11-22.PubMed
- Purugganan MD, Suddith JI: Molecular population genetics of floral homeotic loci. Departures from the equilibrium-neutral model at the APETALA3 and PISTILLATA genes of Arabidopsis thaliana. Genetics. 1999, 151 (2): 839-848.PubMed CentralPubMed
- Nurmberg PL, Knox KA, Yun BW, Morris PC, Shafiei R, Hudson A, Loake GJ: The developmental selector AS1 is an evolutionarily conserved regulator of the plant immune response. Proc Nat Acad Sci USA. 2007, 104 (47): 18795-18800. 10.1073/pnas.0705586104.PubMed CentralPubMed
- Denance N, Sanchez-Vallet A, Goffner D, Molina A: Disease resistance or growth: the role of plant hormones in balancing immune responses and fitness costs. Front Plant Sci. 2013, 4: 155-PubMed CentralPubMed
- Solano R, Stepanova A, Chao Q, Ecker JR: Nuclear events in ethylene signaling: a transcriptional cascade mediated by ETHYLENE-INSENSITIVE3 and ETHYLENE-RESPONSE-FACTOR1. Genes Dev. 1998, 12 (23): 3703-3714. 10.1101/gad.12.23.3703.PubMed CentralPubMed
- He Y, Gan S: Identical promoter elements are involved in regulation of the OPR1 gene by senescence and jasmonic acid in Arabidopsis. Plant Mole Biol. 2001, 47 (5): 595-605. 10.1023/A:1012211011538.
- Chen H, Xue L, Chintamanani S, Germain H, Lin H, Cui H, Cai R, Zuo J, Tang X, Li X, Guo H, Zhou JM: Ethylene insensitive3 and ethylene insensitive3-like1 repress salicylic acid induction deficient2 expression to negatively regulate plant innate immunity in Arabidopsis. Plant Cell. 2009, 21 (8): 2527-2540. 10.1105/tpc.108.065193.PubMed CentralPubMed
- Chehab EW, Kim S, Savchenko T, Kliebenstein D, Dehesh K, Braam J: Intronic T-DNA insertion renders Arabidopsis opr3 a conditional jasmonic acid-producing mutant. Plant Physiol. 2011, 156 (2): 770-778. 10.1104/pp.111.174169.PubMed CentralPubMed
- Franco-Zorrilla JM, Lopez-Vidriero I, Carrasco JL, Godoy M, Vera P, Solano R: DNA-binding specificities of plant transcription factors and their potential to define target genes. Proc Nat Acad Sci USA. 2014, 111 (6): 2367-2372. 10.1073/pnas.1316278111.PubMed CentralPubMed
- Baxter L, Jironkin A, Hickman R, Moore J, Barrington C, Krusche P, Dyer NP, Buchanan-Wollaston V, Tiskin A, Beynon J, Denby K, Ott S: Conserved noncoding sequences highlight shared components of regulatory networks in dicotyledonous plants. Plant Cell. 2012, 24 (10): 3949-3965. 10.1105/tpc.112.103010.PubMed CentralPubMed
- Lamesch P, Berardini TZ, Li D, Swarbreck D, Wilks C, Sasidharan R, Muller R, Dreher K, Alexander DL, Garcia-Hernandez M, Karthikeyan AS, Lee CH, Nelson WD, Ploetz L, Singh S, Wensel A, Huala E: The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2012, 40 (Database issue): D1202-D1210.PubMed CentralPubMed
- Team RC: R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. 2013, Vienna, Austria, http://www.R-project.org,
- Saito R, Smoot ME, Ono K, Ruscheinski J, Wang PL, Lotia S, Pico AR, Bader GD, Ideker T: A travel guide to cytoscape plugins. Nat Met. 2012, 9 (11): 1069-1076. 10.1038/nmeth.2212.
- Morris JH, Apeltsin L, Newman AM, Baumbach J, Wittkop T, Su G, Bader GD, Ferrin TE: clusterMaker: a multi-algorithm clustering plugin for Cytoscape. BMC Bioinformatics. 2011, 12: 436-10.1186/1471-2105-12-436.PubMed CentralPubMed
- Enright AJ, Van Dongen S, Ouzounis CA: An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res. 2002, 30 (7): 1575-1584. 10.1093/nar/30.7.1575.PubMed CentralPubMed
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mole Biol. 1990, 215 (3): 403-410. 10.1016/S0022-2836(05)80360-2.
- Mount DW: Using the Basic Local Alignment Search Tool (BLAST). CSH Protoc. 2007, 2007: pdb top17-PubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.