Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Root Transcriptome Analysis of Wild Peanut Reveals Candidate Genes for Nematode Resistance

Abstract

Wild peanut relatives (Arachis spp.) are genetically diverse and were adapted to a range of environments during the evolution course, constituting an important source of allele diversity for resistance to biotic and abiotic stresses. The wild diploid A. stenosperma harbors high levels of resistance to a variety of pathogens, including the root-knot nematode (RKN) Meloidogyne arenaria, through the onset of the Hypersensitive Response (HR). In order to identify genes and regulators triggering this defense response, a comprehensive root transcriptome analysis during the first stages of this incompatible interaction was conducted using Illumina Hi-Seq. Overall, eight cDNA libraries were produced generating 28.2 GB, which were de novo assembled into 44,132 contigs and 37,882 loci. Differentially expressed genes (DEGs) were identified and clustered according to their expression profile, with the majority being downregulated at 6 DAI, which coincides with the onset of the HR. Amongst these DEGs, 27 were selected for further qRT-PCR validation allowing the identification of nematode-responsive candidate genes that are putatively related to the resistance response. Those candidates are engaged in the salycilic (NBS-LRR, lipocalins, resveratrol synthase) and jasmonic (patatin, allene oxidase cyclase) acids pathways, and also related to hormonal balance (auxin responsive protein, GH3) and cellular plasticity and signaling (tetraspanin, integrin, expansin), with some of them showing contrasting expression behavior between Arachis RKN-resistant and susceptible genotypes. As these candidate genes activate different defensive signaling systems, the genetic (HR) and the induced resistance (IR), their pyramidding in one genotype via molecular breeding or transgenic strategy might contribute to a more durable resistance, thus improving the long-term control of RKN in peanut.

Introduction

Peanut (Arachis hypogaea) is the fourth most important oil crop in the world and is a major source of edible oil and digestible protein, especially in Asia, Africa and Americas. This crop also has a significant role in sustainable agriculture in terms of global food security and nutrition, fuel and energy, sustainable fertilization, and enhanced agricultural productivity as a rotation crop [1]. The narrow genetic base of the cultivated genepool together with the tetraploid and complex nature of the genome are a challenge to breeding efforts to develop cultivars with multiple resistances, high quality and yield [2]. Diversely, peanut diploid wild relatives (Arachis spp.) are important sources of resistance genes due to their higher genetic diversity, and because they have been selected during evolution in a range of environments and biotic stresses, constituting a rich source of new alleles to be introgressed into the cultivated species [3].

Peanut yield is worldwide affected by fungi, bacteria, virus and the root-knot nematode (RKN) Meloidogyne arenaria, which causes substantial yield losses at relatively low population densities [4]. The current approach to manage M. arenaria in the crop include chemical control, adequate cultural practices and most recently, the use of resistant cultivars, which once integrated with good pest management practices, can reduce the risk of pest outbreaks and the expense and environmental toxicity of nematicides [5]. Up to now, three peanut nematode resistant cultivars (COAN, Nematan and Tifguard) have been developed, all by the introgression of a chromosome region from the wild diploid species A. cardenasii [6, 7], for which recently, tightly linked genetic markers have been developed [4, 8]. Whilst quantitative sources of resistance to RKN have been found in peanut, strong resistance probably derived from a small number of genes, such as displayed by Hypersensitive Reaction (HR), has only been identified in diploid species [9, 10]. Therefore, the identification of novel genes related to plant defense against RKN in wild Arachis species is desirable, to broaden the sources of resistances deployed in peanut, and also secure greater durability of the resistance.

Resistance to RKN in various species is often mediated by HR, which occurs soon after invasion and is characterized by the generation of reactive oxygen species (ROS) in plant root cells associated with the presence of second stage juveniles (J2), with the Mi-1-mediated resistance in tomato being the most well characterized [11]. In previous studies [10], we showed that the wild peanut relative A. stenosperma is resistant to M. arenaria, with its mechanism of resistance operating in at least two stages: firstly, only few nematodes intrude the resistant roots, and then, those that successfully penetrated, were killed by HR in the surrounding cells. Transcriptional profiling studies in this species have identified a number of genes potentially related to this defense response [1214]. Nevertheless, the molecular events and the functional role of these genes are yet to be clarified. In this regard, high-throughput mRNA sequencing (RNA-Seq) can be a complementary strategy to predict the roles and interactions of genes, and help the elucidation of complex signaling pathways activated in response to nematode parasitism. This technique enables a comprehensive transcriptome survey of induced changes in gene expression, and has already been successfully applied for the identification of genes responsive to diverse biotic and abiotic stresses, or during specific developmental stages in peanut [2, 1517].

In this work, a genome-wide overview of gene expression during the first stages of the incompatible A. stenosperma/M. arenaria interaction was obtained using RNA-Seq, aiming to better understand the ways in which the resistant host transcriptome respond to the attack of this endoparasite and identify novel genes potentially involved with this nematode resistance. This approach provided new insights into the defense mechanisms used by this species to avoid infection and the ultimate development of the HR response. The understanding of these mechanisms and genes underlying this response will facilitate their transference to more adapted cultivars through the use of tightly linked genetic markers or via plant transformation, leading to improved peanut varieties with higher and more durable resistance to RKN.

Material and Methods

Plant material and libraries construction

A. stenosperma (accession V10309) plants challenged with M. arenaria race 1 were obtained as described before [14]. Shortly, plants were maintained under greenhouse conditions at the University of California, Riverside, USA, and approximately 20,000 juveniles of M. arenaria (J2) were resuspended in 1 ml of fresh deionised water and pipetted into soil depressions around each root of the 4-week-old plants. Control roots were collected at day zero and inoculated roots were collected at 3, 6 and 9 days after inoculation (DAI) in liquid nitrogen for total RNA extraction, using the lithium chloride modified protocol [18], and then purified using Invisorb Plant RNA Mini Kit (Invitek, Germany). Equal amounts of total RNA per collecting point (control and 3, 6 and 9 DAI) were pooled from five different plants, forming four combined samples with two independent biological replicates each. For the paired-end cDNA libraries construction and sequencing, services of FASTERIS (www.fasteris.com) were used employing the mRNA-seq and TruSeq (TM) SBS v5 protocols (Illumina, USA) on a Hi-Seq 2000 sequencing system.

Sequence processing and analysis

The de novo contig assembly of A. stenosperma reads was performed by FASTERIS using Velvet assembler (http://www.ebi.ac.uk/~zerbino/velvet/) that takes in short read sequences removes errors and produces high quality unique contigs, followed by its additional module specialized in transcriptome assembly, OASES, which allows ambiguity assembly, such as alternative splicing or polymorphism (http://www.ebi.ac.uk/~zerbino/oases/). To validate the de novo assembly and estimate the number of assembled reads, a BWA mapping (http://bio-bwa.sourceforge.net/) was carried on with the de novo contigs for the best assembly selected, using 1 million pairs of each library. A complete mapping was only carried for the best selected assembly. The coverage of each transcripts of the best oases assembly was computed using the BEDtools software (http://code.google.com/p/bedtools).

Differential expression analysis of contigs was performed with DESeq package (http://www.bioconductor.org/biocLite.R) using the contig count matrix built from the Bedtools coverage output. Scaling factor for each given sample in each performed pairwise treatment-control was computed using the “estimateSizeFactors” function. Dispersion estimate call was performed with the “estimateDispersions” function from DESeq package for null model evaluation with few replicates with parameters method = "pooled", fitType = "parametric" and sharingMode = "fit-only". The read counts for each contig in each treatment-control condition were fitted to a negative binomial (NB) distribution via the “nbinomTest”- algorithm that involves borrowing information from count data for the complete set of genes. Transcript abundance for each gene was then inferred from the NB mean, in combination with a normalisation obtained by matching the distribution of read counts over a subset of genes which are likely not to be Differentially Expressed Genes (DEGs) according to the model fit. Thus, the null hypothesis corresponding to no differential expression is that the transcript abundance is the same in both conditions. P-values for each estimate of False Discovery Rates (FDRs) were added using the Benjamini-Hochberg algorithm.

For each time point (3, 6 and 9 DAI), the number of DEGs in the infected libraries was compared to the control samples to estimate the Fold Change (FC). Genes were considered as significant DEGs when their relative gene expression levels showed at least 4-FC (log2FC ≥2.0 or ≤ -2.0) difference between inoculated and control samples, with P-value<0.05. Clustering analysis of the significant DEGs based on common expression patterns was conducted using the Multiexperiment Viewer (MeV) (http://www.tm4.org/mev.html). K-means clustering was performed using TM4: MeV 4.7 software (http://www.tm4.org/mev.html) and Pearson’s correlation coefficient.

Similarity search and functional annotation

Functional annotation of contigs was performed by sequence similarity searches using BLASTX against NCBI’s non-redundant sequence database, and similarities were considered significant with E-values ≤1e-7 (S1 Table). InterProScan (http://www.ebi.ac.uk/Tools/InterProScan/) was employed to perform protein domain and motif searches. Gene ontology (GO) terms were assigned by Blast2GO (https://www.blast2go.com/). A. stenosperma putative Transcription Factors (TFs) were identified from the functional BLASTX annotation, and classified according to their TF family using KEGG (http://www.genome.jp/kegg/).

For the MapMan ontology, all A. stenosperma DEGs were compared to Glycine max Affymetrix chip (http://soybase.org/AffyChip/), using BLASTN (E-values ≤1e-7) default parameters. Only A. stenosperma DEGs with identity hits higher than 70% were used for the analysis, and to compile a MapMan ontology (http://mapman.gabipd.org/) classification. Pictorial representations for the biotic stress pathways were uploaded from the MapMan website.

Expression analysis by qRT-PCR

For qRT-PCR analysis, the Platinum® SYBR® Green qPCR Super Mix-UDG w/ROX kit (Invitrogen, USA) was used according to manufacturer's recommendations on ABI 7300 Real-Time PCR System (Applied Biosystem Foster City, USA). A. stenosperma RNA samples were reassembled in pools of five random individuals with two independent biological replicates formed per collecting point (control and 3, 6 and 9 DAI). A total of 2 μg of each RNA pool was treated with 2 U of DNase (Fermentas, Germany) and reverse transcribed using the Super Script II enzyme and oligo(dT) 20 primer (Invitrogen, USA), according to manufacturer’s instructions. Specific primer pairs were designed for the selected candidate genes (S2 Table) and PCR reactions carried out as previously described [14]. Average cycle threshold (Cq) values were estimated using the online real-time PCR Miner tool [19] and normalized to two reference genes (As60S and AsGAPDH), as previously established [18]. Expression ratios of mRNA transcripts at 3, 6 and 9 DAI, relative to control were determined and statistically tested using REST 2009 ver. 2.0.13 software [20]. For qRT-PCR analysis using contrasting genotypes, a new bioassay was conducted in greenhouse conditions essentially as described above using M. arenaria resistant A. stenosperma (V10309) and susceptible A. hypogaea cv Runner. Total RNA from control and 6DAI inoculated roots were assembled in pools of three individuals to form two biological replicates per treatment.

Results and Discussion

RNA-sequencing and de novo assembly

Non-inoculated and M. arenaria race 1 inoculated (3, 6, 9 DAI) A. stenosperma roots were collected for total RNA extraction. Two biological replicates per collecting point were subsequently formed by pooling total RNA samples. The eight cDNA libraries constructed were pair-ended sequenced using Illumina Hi-Seq 2000, with two technical replicates for each library, in two different sequencing runs. A total of 141.3 million raw reads was produced, constituting 28.2 GB of cDNA sequences. An overview on the data of the RNA-seq experiment and contig mapping in the best de novo assembly is shown in Table 1 and S1 Fig. This Transcriptome Shotgun Assembly project (TSA) has been deposited at DDBJ/EMBL/GenBank under the accession GDBK00000000 (PRJNA284674). The version described in this paper is the first version, GDBK01000000.

thumbnail
Table 1. Summary of the A. stenosperma Illumina Hi-Seq sequencing data.

https://doi.org/10.1371/journal.pone.0140937.t001

The de novo transcript assembly consisted of 37,882 loci which include a combined total of 44,132 contigs (S1 Table). In total 3,548 of the loci had multiple contigs consisting of different combinations of splice variants or different regions of the same loci due to partial assembly. Most contigs were 200–300 bp in length, with more than 37,000 contigs common among all libraries (S1 Fig). Overall, 34,489 sequences (78.14%) showed similarity to proteins in the NCBI databank using BLASTX algorithm (E value < e-7) (S1 Table). The number of loci assembled for A. stenosperma with gene structure found herein, is comparable to those described in other tropical legumes such as soybean (43,367) [21], chickpea (42,141) [22], pigeonpea (48,680) [23] and common bean (39,572) [24].

To date, publicly available transcriptome data for A. stenosperma comprise 6,264 ESTs and 7,723 contigs originated from a Transcriptome Shotgun Assembly (TSA) both produced by this group [25, 26]. In this study, the sequencing of root samples during the early stages of the A. stenosperma/RKN interaction, associated to the robustness of the Illumina Hi-Seq technique, enabled a 3-fold increase on the number of loci identified for this species, which is known to harbor resistance to a number of pests but lacks a reference genome sequenced. This study will, therefore, expand considerably the transcriptomic resources available to be used for gene discovery, molecular genetics and functional genomics in the genus Arachis.

Functional annotation

The annotation of the A. stenosperma contigs and their assignment into functional categories based on conserved PFAM domain predictions and gene ontology (GO) allowed the classification of almost 30% (13,157) of the 44,132 contigs in the three main GO categories: biological process (41.59%), molecular function (33.74%) and cellular component (24.66%) (Fig 1). This distribution was in accordance with previous analysis on peanut and wild Arachis transcriptome [16, 2628], with proteins related to metabolic processes (15%) and cellular processes (12%) enriched in the biological process, whilst cell (8%) and cell parts (7%) were the most highly represented in the cellular components category. Interestingly, within the molecular function category, binding proteins (14%) and proteins related to catalytic activity (11%) were the most enriched, with only 207 contigs annotated in the response to stimulus category (0.5%) (Fig 1). This included genes with putative defense function regulated during RKN infection, such as peroxidases, chitinases, extensins and proteinase inhibitors, perhaps playing a global response to pathogen invasion or specific response to the nematode, including serine–threonine protein kinases and glycosyltransferases [29]. Genes related to regulation of hormones, such as auxin repressed proteins (ARP) and cytokinin oxidase, previously related to A. stenosperma response to M. arenaria were also identified [14, 30].

thumbnail
Fig 1. Gene Ontology (GO) terms of A. stenosperma transcriptome during the interaction with M. arenaria.

https://doi.org/10.1371/journal.pone.0140937.g001

Transcription factors distribution

Transcriptional factors (TFs) are of foremost importance, as gene induction by stress occurs primarily at the level of transcription. TFs regulation of the temporal and spatial expression of these genes is an important part of downstream signaling pathways to initiate protective stress responses, as seen in peanut and other legumes [3133]. A. stenosperma TFs here identified were classified into 23 families according to KEGG, and those related to biotic stress response were amongst the most abundant, including MYB (22%), EREBP (13%), bZIP (12%), HD-ZIP (10%), Heat shock (9%), Zinc finger (8%), CCAAT binding (7%) and WRKY (2%) (Fig 2).

thumbnail
Fig 2. Distribution of A. stenosperma Transcription Factors families during the interaction with M. arenaria.

https://doi.org/10.1371/journal.pone.0140937.g002

The three most expressed TF families in this study MYB, EREBP and bZIP, which are implicated in many cell physiological and biochemical processes, have also been described as being involved in different plant-pathogen interactions, through the regulation of pathogen resistance and cell death [33]. Thus, it is possible that the TFs highly expressed here have closely overlapping functions, and might activate or repress genes through cis-acting sequences that respond to the nematode infection.

Clustering analyses of the differentially expressed genes

Differential expression of A. stenosperma contigs was achieved with DESeq package that performs variance estimation and differential expression of the raw read counts, using the contig count matrix built from the Bedtools coverage output. Therefore, at each sampling point, the number of differentially expressed genes (DEGs) in the infected libraries was compared to the control (Fold Change-FC). Genes were considered to be differentially expressed when relative gene expression levels showed at least 4-FC difference (log2FC ≥2.0 or ≤ -2.0) between inoculated and control samples, with and P-value<0.05. In total, 958 (2.17%) significant DEGs were identified among the 44,132 contigs in at least one of the sampling points (3, 6 or 9 DAI), indicating that, the great majority of contigs were evenly expressed between inoculated and control roots during the assay time course.

A Heat Map representation of the clustering analysis of the global expression of these DEGs indicated a general downregulation during the initial time points of the M. arenaria infection (Fig 3A). Numbers of DEGs greatly differed at each stage of this plant-pathogen interaction and point up to a considerable imbalance of gene expression during nematode infection. At 3 DAI, the majority of DEGs (268) were dowregulated in A. stenosperma infected roots, whilst 165 were upregulated. The profile changed more drastically at 6 DAI, when the number of downregulated DEGs (473) is 52-fold higher than that of upregulated genes (9). At subsequent stage (9 DAI), most of the genes (91) was still repressed, and only few (13) were induced. This global gene downregulation at early stages of RKN interaction is in agreement with the expression profiles previously reported for peanut and A. stenosperma [30, 34], as well as for other legumes.

thumbnail
Fig 3. HeatMap of A. stenosperma Differentially Expressed Genes (DEGs).

A) Distribution of A. stenosperma DEGs significantly (FC ≥2.0 or ≤ -2.0 and p˂0.05) up and downregulated during M. arenaria infection with in at least one of the sampling points (3, 6 or 9 DAI). Normalised values, relative to control, are shown in a red–green scale; the brightest red being the most upregulated and the brightest green, the most downregulated genes; B) DEGs K-means clustering analysis, with the pink lines indicating the representative expression profile of the Cluster; x- axes represent days after nematode inoculation (DAI) and y-axes, log2FC of mRNA levels of DEGs in the sampling points (3, 6 or 9 DAI) normalized to control.

https://doi.org/10.1371/journal.pone.0140937.g003

The expression trends at the early stages of M. arenaria/A. stenosperma interaction were corroborated by subsequent MeV clustering analysis (P< 0.05 and log2FC ≥2.0 or ≤ -2.0), which gathered the 958 significant DEGs into three distinct clusters based on common expression patterns (Fig 3A and 3B). Cluster I grouped 284 DEGs downregulated at 3 DAI, with a slight upregulation at 6 and 9 DAI (Fig 3B). Contrastingly, Cluster II showed 44 DEGs which were induced only at 3 DAI, followed by a steady downregulation at the two other points. The larger group, Cluster III, was composed of 630 DEGs displaying an upregulation only at 3 DAI, followed by a strong downregulation at 6 and 9 DAI. Representatives of each cluster were selected for further validation of their expression profile using qRT-PCR.

Regardless the trend of the expression pattern (up or down), the great majority of genes showed differential regulation at 6 DAI (482), followed by 3 DAI (433) and 9 DAI (104), with 61 showing significant regulation in more than one sampling point. Most commonly regulated genes were identified between 6 and 9 DAI, with few in common to 3 DAI, as most DEGs were downregulated at 6 and 9 DAI, but upregulated at 3DAI (Fig 3A and 3B).

Several studies, including those in Arachis, have identified genes responsive to nematode infection, in both, the compatible and incompatible interactions, which regulation vary dramatically due to, among other factors, the time elapsed since pathogen perception [30, 34 35]. In the resistant A. stenosperma, the great majority of genes are repressed in response to the infection, which is in accordance to studies on other plant-RKN interactions [3537], and on the soybean cyst nematode [38]. This suggests that the downregulation of a significant number of genes might be necessary for the onset of the specific HR in A. stenosperma, as in other resistant genotypes, which is the main mechanism of resistance to M. arenaria [10, 39].

MapMan ontology

In order to obtain an overview of the metabolic pathways or other processes underlined by the 958 DEGs, we used the MapMan tool [40]. For that, due to the lack of an Affymetrix chip containing Arachis spp. sequences, only A. stenosperma DEGs displaying sequence identity higher than 70% with the G. max Affymetrix chip sequences were considered for the analysis. A total of 405 genes (42.27%) with sufficient homology to the G. max genome probes were annotated into 20 MapMan Bins. A diagram displaying the expression of these DEGs into biotic stress pathways and processes is shown in Fig 4. Using MapMan annotation software to organize and display our data sets in the context of biological pathways we found 345 genes related to biotic stress responses among distinct DEGs, with a large number of genes involved in signaling (114), including 39 TFs. Interestingly, those TF genes revealed in MapMan analysis belong to TF families, as MYB, EREBP and b-ZIP (Fig 4), that are also the most abundant families previously identified in the TF distribuition analysis (Fig 2). Also, genes related to pathogen recognition and defense response (9), cell wall modification (58), hormonal metabolism (36) and proteolysis (67), enzymes of the phenylpropanoid pathway and others related to secondary metabolites (22), which are well known to be involved in the defense response to nematode infection, could be mapped (Fig 4).

thumbnail
Fig 4. MapMan visualization of A. stenosperma genes involved in the response to M. arenaria infection.

Overview of DEGs expression patterns (log2FC of mRNA levels) in infected roots relative to control. Dots show the different paralogous genes encoding proteins that are related to a certain defense response step. Red dots indicate upregulation and green dots downregulation.

https://doi.org/10.1371/journal.pone.0140937.g004

Overall, MapMan ontology enabled a genome-wide expression outline of A. stenosperma genes responding to M. arenaria infection, by identifying pathways that are inlvolved in the main steps leading to the HR response. Innitially, pathogen recognition by R genes (Fig 4) triggers Effector-Triggered Immunity (ETI), eliciting reactive oxygen intermediates and ion fluxes (redox state species and peroxidases) and also enzymes that protect the cell from the attack by reactive oxygen (e.g. Glutathione-S-transferase). Following, the upregulation of genes engaged in signal transduction pathways (e.g. ERF, b-ZIP, MYB) promote the induction of a large set of genes commonly known as defense related genes (e.g. Heat Shock Proteins, Proteolitic proteins), Pathogenesis related proteins (PR) and enzymes involved in the synthesis of anti-microbial compounds called phytoalexins (secondary metabolites). Genes involved in phytohormone signaling pathway directly related to plant defense response (Jasmonic acid (JA), Salycilic Acid (SA)) and others that play a role in plant imunity (ABA and Auxin) were also upregulated (Fig 4). Some of the proteins and transcription factors regulated here are also known to be activated by environmental factors (abiotic stress), whilst PAMP-triggered immunity (PTI) activated by pathogen associated molecular patterns (PAMPs), activate cell wall enzymes, induce proteolysis and transcriptional reprogramming as a common plant immune responses (Fig 4).

Genes expression profiling

From the 958 DEGs identified, representatives of each cluster (5%) identified in MeV cluster (Fig 2) were selected as candidate genes for expression profile validation by qRT-PCR. Within each cluster, candidate genes were chosen based on their fold-change level (FC), (positive- or negatively regulated), assuming that those with the highest fold-change ratio would be more likely to have a biological role in this plant-nematode interaction, and also based on their homology in BLASTX to proteins with a functional role in plant defense mechanisms. From those, 27 candidate genes showed significant (P<0.05) differencial expression in at least on time point by RT-qPCR analysis (Table 2 and S2 Table). Among these candidates, representatives of different steps of the plant defense response, previously identified in MapMan analysis (Fig 4), were found, including: (i) genes associated to the perception of PAMPs such as glucan binding proteins (AsBger) that perceives beta-glucans from microbes; (ii) enzymes that degrade pathogen cell wall such as chitinases (AsCHI2); (iii) resistance genes such as NBS-LRR proteins (AsTIR-NBS-LRR) that recognize pathogen effectors and (iv) TFs previously identified in silico as the most abundant and related to plant responses to stresses (AsMYB25; AsWRKY49; AsERF6) (Fig 2). In addition, genes codifying for (v) key enzymes on the jasmonic and salicylic acids metabolic pathways, which are triggered as plant responses to pathogen attack, such as ureases (AsUreD) and Response Regulator (AsPRR37), (vi) enzymes related to Reactive Oxygen Response (ROS), such as allene oxidase (AsAOC3) and (vii) enzymes related to hormonal balance (AsGH3.1) and (ix) cell wall composition and trafficking, such as two cell wall hydrolases (AsCWAH and AsCWAH2), were also found to be significantly differentially regulated.

thumbnail
Table 2. Differentially expressed candidate genes in A. stenosperma roots infected with M.arenaria.

https://doi.org/10.1371/journal.pone.0140937.t002

All 27 candidate genes analyzed by qRT-PCR showed specificity of transcript amplification, with high amplification efficiencies (S2 Table). Their mRNA levels in A. stenosperma roots inoculated with M. arenaria relative to control are shown in Fig 5. The majority of candidate genes (17) displayed a gradual increase in their expression as the nematode infection progressed, achieving a positive regulation at 9 DAI (Fig 5). AsAOC, showed an unique profile, as it also increased its expression along the assay but was still negatively regulated at 9 DAI. Their expression profile corresponded to genes previously classified in Cluster I (Fig 3). The opposite expression profile, consisting of a gradual decrease on gene expression along the assay, was observed only for two genes (AsALKBH2 and AsTAT; Fig 5) representatives of Clusters II (Fig 3). For the remaining eight candidate genes (AsAUX/IAA, AsBerg, AsBTB, AsCHI2, AsCOX1, AsERF6, AsKEL and AsTMV), 6 DAI was the inflection point of the curve representing their expression profiles (U- and V-shaped) (Fig 5), with (AsAUX/IAA and AsKEL) being representatives of Cluster III.

thumbnail
Fig 5. A. stenosperma candidate genes expression profile by qRT-PCR.

Relative mRNA levels of 27 candidate genes in A. stenosperma roots inoculated with M. arenaria collected at 3, 6 and 9 days (DAI), relative to control. Bars represent the standard deviation of two biological replicates. Significantly (P < 0.05) up- or downregulated genes are indicated *.

https://doi.org/10.1371/journal.pone.0140937.g005

As a follow-up analysis, the expression profiles of these candidate genes by qRT-PCR and in silico were compared. For that, graphs displaying the transcript abundance for each contig inferred in silico (log2 transformed relative mRNA levels) and the relative quantification obtained for the corresponding transcript by qRT-PCR (log2 transformed) were built (Fig 6). An additional set of nine genes previously identified in our studies as interesting candidates for nematode resistance [14, 26] was also included to broaden the analysis (AsMG13, AsU-BOX, AsLIP, AsDC1, AsPN, AsINT, AsTET, AsRS2 and AsEXLB) (S3 Table). When primers designed for these additional candidate genes [14] (S3 Table) amplified by electronic PCR more than one contig in our RNA-Seq database, only the contig with the expression profile most similar to the qRT-PCR was considered (Fig 6).

thumbnail
Fig 6. Comparisons between RNA-seq and qRT–PCR expression profile.

Log2 transformed relative mRNA levels of DEGs in roots of A. stenosperma inoculated with M. arenaria collected at 3, 6 and 9 DAI relative to control, analysed by qRT-PCR (left Y-axis) and by RNA-Seq right Y-axis.

https://doi.org/10.1371/journal.pone.0140937.g006

Overall, 17 DEGs showed a comparable expression profile between both analyses, with qRT-PCR results showing to be positively correlated to RNA-Seq data analysis (Pearson correlation coefficients r = 0.68). This set of genes exhibited distinctive responses to nematode challenge, and include those related to phytohormones signaling, TFs, cell wall structure and trafficking, which participate in plant resistance to nematodes, including the onset of the HR response, the resistance mechanism described for A. stenosperma against M. arenaria [10].

The association of their expression profiles in response to nematode infection with their primary functions in defense might contribute to unravel their putative roles during this incompatible plant-nematode interaction, as described below.

Genes involved in the salycilic acid pathway.

After perception of pathogen or microbe–associated molecular patterns (PAMPS), which are sensed through leucine-rich repeat (LRR) receptors, plants induce stomatal closure and produce innate immunity responses. If the plant has a specific resistance gene (R) that recognizes the corresponding avirulence (avr) gene from the pathogen, a rapid defense mechanism known as the HR occurs to prevent infection. In general, salycilic acid (SA) activates the defense response to biotrophic and hemi-biotrophic pathogens, induces Systemic Acquired Resistance (SAR), and triggers the expression of SAR-associated pathogenesis related genes (PRs). SA has already been described as playing an important role in decreasing susceptibility to RKN, as it can cause accumulation of Reactive Oxygen Species (ROS), which is observed during the early stages of incompatible interactions and it is one of the components of the HR reaction [4143].

From the relatively few endoparasite nematode resistance (Nem-R) genes isolated so far, the vast majority fall into the NBS-LRR class of R-genes, including the MI gene from tomato, which confers resistance to M. incognita and M. javanica and five Me genes in pepper, that control different species of Meloidogyne [39, 44, 45]. In the A. stenosperma/ M. arenaria interaction, several typical cell features, including densely stained cytoplasm and altered organelle structures were observed in the central cylinder, indicating a R-gene mediated HR of infested host cells [10]. In the present study, two members of the NBS-LRR class of R-genes, AsTIR-NBS-LRR and AsMG13 showed an increased level of expression in A. stenosperma roots, in concert with the onset of the HR response, what suggests their involvement in the early stages of defense response in this incompatible interaction (Fig 6). R-genes might detect directly the pathogen avirulence effector or recognize changes in the host protein elicited by the pathogen effector (guard hypothesis) [46], and can be used in combinations to be deployed in cultivars to maximize their durability against pathogen changes. The identification of new distinct R genes is, therefore, an important topic in most crop improvement programmes, and AsTIR-NBS-LRR and AsMG13 constitute, thus, interesting candidate genes for further functional characterization.

High levels of SA that often lead to HR are regulated by other proteins, such as lipocalins, which are involved in modulating tolerance to oxidative stress [47]. These proteins have already shown to act as scavengers of potentially harmful molecules known to be induced by extreme temperature stress and excessive light [48]. In our study, the lipocalin AsLIP gene expression is induced after nematode perception (3 DAI) following by a sharp reduction on its expression along the assay (Fig 6), which might contribute to the outbreak of the HR in this species. In our previous studies, another known ROS plant scavenger, a metallothionein homolog (AsMET) also had its expression diminished during the A. stenosperma/M. arenaria interaction [30], reinforcing the potential role of these molecules in this defense mechanism.

Salicylic acid is a known secondary metabolism inducing agent, including the phenylpropanoid pathway that, on its turn, produces a variety of structural and defense-related phenolics including resveratrol [49]. Resveratrol Synthase (RS) mediates the synthesis of resveratrol and is induced by different stresses, [48] constituting an important regulator for initiation of HR-related cell death [50]. Peanut is one of the few crop plants that is able to synthesize resveratrol and the accumulation of RS and its mRNA has been previously observed in wild Arachis in response to treatment with ultraviolet light [51]. In our previous studies [14], the RS homolog AsRS2, showed upregulation along the nematode infection (3, 6 and 9 DAI), whilst here, similar trend was observed in silico at at 6 and 9 DAI (Fig 6), both coinciding with the occurrence of HR in A. stenosperma. Together with other genes involved in the SA pathway described here, AsRS2 is suggested to be involved in the HR response of A. stenosperma to M. arenaria.

In addition to this immediate, gene-a-gene contest, plants also employ general resistance mechanisms induced after the HR to combat secondary infections from a broad spectrum of pathogens, known as systemic acquired resistance (SAR). This mechanism is commonly associated with activation of genes that code for defense related proteins, including chitinases, that have been shown to degrade chitinous components of eggs of M. incognita causing a significant reduction on their development and reproduction rate [52]. Here, we found a chitinase homolog (AsCHI2) upregulated at 6 and 9 DAI (Fig 6), which we suggest as having a part on the SAR response of A. stenosperma to the pathogen.

Genes involved in the jasmonic acid pathway.

Besides the genetic resistance mediated by R-genes, induced resistance (IR), activated under the influence of metabolites of plant pathogens and various biotic and abiotic factors called elicitors also occurs. Jasmonic acid (JA) and ethylene are important for the induction of nonspecific disease resistance, and play an important role in IR activating genes that are involved in plant defense in response to wounding and pathogen attack [53]. Jasmonates are closely associated with the production of secondary metabolites, such as flavonoids, and can also trigger synthesis of a wide range of defensive proteins in roots, which are toxic to nematodes in vitro, and have been implicated in plant defense against these parasites [54]. One of these compounds are Patatin-like proteins, that play an essential role in HR and cell death during plant-pathogen interactions, as these proteins mediate the production of JA and defense responses in some species [55]. Previous works have shown the upregulation of patatin in response to M. arenaria infection in Arachis spp. [14, 34]. Here, the induction of AsPN was confirmed both in silico and in RT-qPCR (Fig 6), at 6 and 9 DAI, the time points corresponding to the induction and aftermath of the HR in A. stenosperma, which reinforces its role in this defense response.

Another protein involved in mounting resistance responses via the accumulation of JA is allene oxidase cyclase (AsAOC3), which has been demonstrated to be a functional enzyme involved in the biosynthesis of JA and related compounds [56]. In this study, AsAOC3 showed a good correlation between the in silico and qRT-PCR expression profiles (Fig 6). However, its role in this plant-nematode interaction is not yet clear, because regardless showing a trend of increased expression along the whole assay, this candidate gene was downregulated at all time points analysed (Fig 6).

Genes related to hormonal balance, cellular plasticity and signaling.

Auxin signaling, which is generally recognized to be involved in plant growth and development, is also involved in the complicated network of plant–nematode interactions, with evidences that RKN interfere in auxin biosynthesis in feeding cells, thus, altering local hormone balance in order to form the giant cell complex [57]. Here, an Auxin Responsive Protein (AsAUX/IAA) was downregulated in response to M. incognita infection, in both qRT-PCR and in silico analysis, especially at 6 DAI (Fig 6). Moreover, previous analyses by macroarray and Northern blotting, showed that AsARP, a gene repressed by the presence of auxin, was strongly upregulated in A. stenosperma roots infected with M. arenaria [30]. These findings, the downregulation of an auxin-induced gene, AsAUX/IAA, and the upregulation of an auxin-repressed gene, AsARP, corroborate our previous hypothesis that the lack of auxin accumulation could lead to a failure in giant cell formation in incompatible host-RKN interactions, being part of the suite of defense genes regulated by hormones harbored by A. stenosperma.

Another constituent of the complex network of auxin signaling that regulates stress adaptation responses in plants is GH3), which encodes an IAA-amido synthetase that maintains auxin homeostasis by conjugating excess IAA (indole-3-acetic acid) to amino acids [58]. GH3 also inhibits expansin production, which on its turn is induced by IAA [59]. Expansins are cell wall-associated proteins that contribute for the cell wall plasticity and may also contribute to the plant vulnerability to pathogens. Therefore, the suppression of expansin might contribute to plant resistance to some pathogens. Considering that GH3 is induced by multiple pathogens, its mediated role in nematode resistance might be due, at least partly, to inhibiting the expression of expansins by suppressing auxin signaling [59].

In our study, AsGH3.1, which is an enzyme that prevents the free accumulation of auxin, showed a steady increase on its expression immediately after nematode perception (3 to 6 DAI) (Fig 6). At the same interaction period, an expansin homolog (AsEXLB) displayed a notable decrease on its expression, specially at 6 DAI, when HR is initially observed in the A.stenosperma/M.arenaria interaction [10] (Fig 6). This antagonist gene expression between GH3.1 and expansin homologs has previously been described in rice infected with Xanthomonas oryzae, and proved to be a mechanism which contributed to an increased resistance to the pathogen [58]. In this study, we suggest that AsGH3.1 and AsEXLB may also be interligated and contribute to the host defense response. Therefore, we hypothesize that M. arenaria initially (3 DAI) induces auxin production by the plant, in order to weaken the cell wall as an attempt to establish the giant cell [35]. However, the resistant A. stenosperma in turn, defends itself from auxin-facilitated nematode infection by preventing the loosening of its cell wall via suppressing auxin accumulation, as indicated by AsAUX/IAA and AsARP expression behavior discussed above. This is through the activation of AsGH3.1 by the pathogen that, in its turn, inactivates IAA, by conjugating it to amino acids, preventing cell wall-loosening and nematode penetration, and therefore the establishment of the giant cell by the nematode.

To further support the role of membrane structural proteins in A. stenosperma/M.arenaria interaction, two other membrane-associated genes, coding for Tetraspanin (AsTET) and Integrin (AsINT) also displayed a compatible profile of increased expression in both, in silico and qRT-PCR (Fig 6). Tetraspanin is a membrane spanning protein containing a large extracellular loop, which provides interactive binding sites and reportedly has a role in development, pathogenesis and immune response in animal and plants [60]. The main suggested function for Tetraspanin is cell-to-cell signaling, which is also one of the main roles attributed to Integrins, that play a critical part in cell structure, migration and anchoring cells to the extracellular matrices [61]. A well known Integrin-like protein is the non-race specific disease resistance NDR1 gene in plants, a plasma membrane-localized protein that plays an essential role in resistance mediated by the NBS-LRR proteins, as it is a key signaling component of these responses during pathogen infection [62]. In our study, both AsTET and AsINT showed a steady increase in their expression along plant-nematode interaction, with upregulation after 6 DAI (Fig 6). This contributes to the assumption of their role in A. stenosperma in cell signaling and triggering the plant defense response against nematode parasitism.

Besides structural proteins, cell wall enzymes are thought to be important components of the mechanisms by which plants sense pathogens, including endo-ß-mannanases [63], which overexpression has shown to increase resistance to root fungus in transgenic tobacco [64]. Here, AsMAN was upregulated during the whole period of the nematode-plant interaction analysed (Fig 6). Altough it is well known that carbohydrate moieties on nematode surface are specific elicitors, functioning in the determinative stages of nematode-plant interactions [65], no indicative of mannan eliciting resistance by RKN has yet being described.

Transcription factors and DNA binding elements.

The responses to pathogen invasion require large-scale transcriptional reprogramming, including those in some TF families. In this study, two defense related TFs, AsMYB25 and AsWRKY49, were shown to be upregulated in response to infection in both in silico and RT-qPCR analyses (Fig 6). These candidates are members of MYB and WRKY TF families identified as some of the most abundant in this study (Fig 2), suggesting their putative role in this incompatible plant-nematode interaction. Here, AsWRKY49 increased its expression along the assay (Fig 6) and might be involved in the host response either by directly increasing the expression of PRs, or through the mediation of the NPR1 protein [66]. Likewise, MYB genes have been identified in other plant species displaying resistance to endosedentary nematodes [67], which contributes to the assumption that some representatives of this TF family in A. stenosperma may be involved in the response to M. arenaria parasitism.

Additionally, DC1 domain-containing zinc finger proteins have been demonstrated to play an important role in regulating the defense responses to biotic stresses, as they function as positive regulators of plant cell death programme (PCD) and SA-dependent defense responses [68]. In this study, an A. stenosperma DC1 homolog, AsDC1 is induced at 6 and 9 DAI in response to the nematode infection both in silico and in qRT-PCR. Comparably, similar expression trend was observed for another protein containing a AsU-box, that through ubiquitin ligase activity, may act as regulators in defense against pathogen and PCD in plants [69].

Expression profiling of candidate genes in resistant and susceptible Arachis genotypes.

In order to verify whether the 17 candidate genes selected above (DEGs) (Fig 6) are involved in global responses to RKN infestation in any Arachis species or in specific responses by RKN-resistant A. stenosperma, their expression profiles were analysed in two contrasting Arachis genotypes for M. arenaria resistance: A. stenosperma (resistant) and A. hypogaea cv Runner (susceptible). For that, qRT-PCR analysis was conducted with samples from a new bioassay including control and inoculated roots of the two contrasting Arachis genotypes, using primers designed for A. stenosperma (S2 Table).

For 10 DEGs analyzed, single peaked melting curves were generated, indicating a specific amplicon in both A. hypogaea and A. stenosperma. Contrasting expression behaviors (up or downregulation) between the two genotypes were observed for eight genes (AsOAC, AsAUX/IAA, AsCHI2, AsMAN, AsMYB25, AsTIR-NBS-LRR, AsWRKY49 and AsLIP) (Fig 7). The two remaining DEGs (AsDC1 and AsEXLB), whereas showing a similar response to RKN infection (downregulation for AsDC1 and upregulatuion for AsEXLB), differed considerably in their fold change magnitude in each genotype (Fig 7).

thumbnail
Fig 7. A. stenosperma and A. hypogaea qRT–PCR expression profile.

Relative mRNA levels of 10 candidate genes in A. stenosperma and A. hypogaea roots inoculated with M. arenaria relative to control. Bars represent the standard deviation of two biological replicates. Significantly (P < 0.05) up- or downregulated genes are indicated by*.

https://doi.org/10.1371/journal.pone.0140937.g007

The discongruity in expression behavior of DEGs representing PRs, ETI receptors, and jasmonic acid (JA)- and salicylic acid (SA)-related proteins in RKN-resistant and susceptible genotypes following M. arenaria inoculation, coupled with their putative role in HR response discussed throghout this study, reinforce their involvement in the specific responses by RKN-resistant A. stenosperma, as in contrast to transcripts involved in global responses to RKN infestation in other genotypes.

Current knowledge of signaling mechanisms that operate during plant infestation by parasitic nematodes are limited, especially in wild species. In the present work, a deep transcriptome analysis allowed the identification of differentially expressed genes during the early stages (3 to 9 DAI) of the incompatible A. stenosperma/M. areanaria interaction. These genes, either inlvolved in the PAMP-triggered immunity (PTI) or Effector-Triggered Immunity (ETI), were further validated by qRT-PCR, corroborating their potential as genes of interest. We found that both JA-dependent (patatin, allene oxidase cyclase) and SA-dependent (R genes, lipocalin, resveratrol synthase, chitinase) defense genes and their regulators (MYB, WRKY) were triggered in A. stenosperma roots when infected with M. arenaria. Additionaly, proteins that are involved in the regulation of defense responses (DC1, U-BOX), cell signaling (mannanases, tetraspanin, integrin), cell plasticity (expansin) and hormonal balance (auxin responsive, GH3) were also differentially expressed in infected roots, with a number of them showing divergent expression profile in contrasting genotypes for M. arenaria resistance, reinforcing their role in specific responses by RKN-resistant A. stenosperma. Most of these genes constitute interesting candidates for further functional analysis, and their knockdown or gain-of-function effect in A. rhizogenes-induced hairy roots in composite peanut plants [70], will contribute to decipher their potential role in A. stenosperma defense mechanisms, including HR and SAR against M. arenaria.

The fact that genes participating in both defense hormone pathways (SA and JA) were identified in this plant-nematode interaction probably indicates the absence of a conflict between these signaling systems, as observed in the Mi mediated resistance in tomato [71] and other plant species [72, 73], where cross-talk between different hormone pathways has been shown to optimize the defense response against antagonists. These findings encourage the use of different signaling pathways mediating defense in peanut plants to optimize their defensive responses with both genetic (HR) and induced resistance (IR), which can be used in single or cumulative properties, and contribute to the reduction of its susceptibility to nematodes and possibly, also to other pests. Pyramiding or alternating resistance genes against Meloidogyne spp. has already proved to be effective in other crops [45], and once used to produce improved peanut cultivars, may also contribute to a more durable resistance and improve the long-term control of root-knot nematodes in the species.

Conclusions

Parasitic RKN evolved sophisticated strategies to infect plant causing worldwide yield losses. Sources of host resistance against the most damaging nematode species in peanut (M. arenaria) have been identified in its wild relative A. stenosperma and is due to the onset of the Hypersensitive Response. The present study suggests putative genes that trigger changes in specific signalling pathways either inlvolved in the PAMP-triggered immunity (PTI) or Effector-Triggered Immunity (ETI), and also in structural arrangements involved in A. stenosperma resistance to RKN. The deployment of individual or combined resistance genes by marker assisted selection or via transgenic expression can offer several benefits for nematode control in an integrated management strategy, reducing risks to the environment and to human health, accessibility for food, and the possibility of achieving durable, broad-spectrum nematode resistance.

Supporting Information

S1 Table. A. stenosperma contigs and their best protein homology at NCBI (BLASTX).

https://doi.org/10.1371/journal.pone.0140937.s001

(XLSX)

S2 Table. Primer sequences of the differentially expressed genes (DEGs).

https://doi.org/10.1371/journal.pone.0140937.s002

(DOCX)

S3 Table. Additional differentially expressed candidate genes in A. stenosperma roots infected with M. arenaria.

https://doi.org/10.1371/journal.pone.0140937.s003

(DOCX)

S1 Fig. A.stenosperma ILLUMINA HI-Seq 2000 statistics.

https://doi.org/10.1371/journal.pone.0140937.s004

(TIF)

Acknowledgments

José Francisco Montenegro Valls for providing the germplasm; Ana Zotta Mota for the bioinformatic support; Phillip A. Roberts for the nematode bioassays at Riverside University; Leandro Mesquita for green house assistance.

Author Contributions

Conceived and designed the experiments: PMG ACMB. Performed the experiments: PMG ACMB ACGA LAG ACQM MS CVM SCMLB TNO. Analyzed the data: OBS RCT DJB PMG ACMB LAG. Contributed reagents/materials/analysis tools: PMG ACMB DJB SCMLB. Wrote the paper: PMG ACMB ACGA.

References

  1. 1. Guo B, Chen C, Chu Y, Holbrook C, Ozias-Akins P, Stalker HT. Advances in genetics and genomics for sustainable peanut production. Sustainable Agriculture and New Biotechnologies. 2011:341–67.
  2. 2. Pandey MK, Monyo E, Ozias-Akins P, Liang X, Guimarães P, Nigam SN, et al. Advances in Arachis genomics for peanut improvement. Biotechnology Advances. 2012;30(3):639–51. pmid:22094114
  3. 3. Leal-Bertioli S, Jose A, Alves-Freitas D, Moretzsohn M, Guimaraes P, Nielen S, et al. Identification of candidate genome regions controlling disease resistance in Arachis. BMC Plant Biology. 2009;9:112. pmid:19698131
  4. 4. Burow M, Starr J, Park C-H, Simpson C, Paterson A. Introgression of homeologous quantitative trait loci (QTLs) for resistance to the root-knot nematode (Meloidogyne arenaria (Neal) Chitwood) in an advanced backcross-QTL population of peanut (Arachis hypogaea L.). Molecular Breeding. 2014;34(2):393–406.
  5. 5. Danchin EGJ, Arguel M-J, Campan-Fournier A, Perfus-Barbeoch L, Magliano M, Rosso M-N, et al. Identification of Novel Target Genes for Safer and More Specific Control of Root-Knot Nematodes from a Pan-Genome Mining. PLoS Pathog. 2013;9(10):e1003745. pmid:24204279
  6. 6. Tillman B, Stalker HT. Peanut. In: Vollmann J, Rajcan I, editors. Oil Crops: Springer New York; 2010. p. 287–315.
  7. 7. Holbrook C, Stalker H. Peanut breeding and genetic resources. Plant Breeding Reviews. 2003;22:297–355.
  8. 8. Nagy E, Chu Y, Guo Y, Khanal S, Tang S, Li Y, et al. Recombination is suppressed in an alien introgression in peanut harboring Rma, a dominant root-knot nematode resistance gene. Molecular Breeding. 2010;26(2):357–70.
  9. 9. Bendezu IF, Starr JL. Mechanism of Resistance to Meloidogyne arenaria in the Peanut Cultivar COAN. Journal of Nematology. 2003;35(1):115–8. pmid:19265984
  10. 10. Proite K, Carneiro R, Falcão R, Gomes A, Leal-Bertioli S, Guimarães P, et al. Post-infection development and histopathology of Meloidogyne arenaria race 1 on Arachis spp. Plant Pathology. 2008;57(5):974–80.
  11. 11. Abad P, Williamson VM. Plant Nematode Interaction: A Sophisticated Dialogue. In: Jean-Claude K, Michel D, editors. Advances in Botanical Research: Academic Press; 2010. p. 147–92.
  12. 12. Proite K, Leal-Bertioli S, Bertioli D, Moretzsohn M, Silva F, Martins N, et al. ESTs from a wild Arachis species for gene discovery and marker development. BMC Plant Biol. 2007;7:7. pmid:17302987
  13. 13. Guimarães P, Brasileiro A, Proite K, de Araújo A, Leal-Bertioli S, Pic-Taylor A, et al. A Study of Gene Expression in the Nematode Resistant Wild Peanut Relative, Arachis stenosperma, in Response to Challenge with Meloidogyne arenaria. Tropical Plant Biol. 2010;3(4):183–92.
  14. 14. Morgante CV, Brasileiro ACM, Roberts PA, Guimaraes LA, Araujo ACG, Fonseca LN, et al. A survey of genes involved in Arachis stenosperma resistance to Meloidogyne arenaria race 1. Functional Plant Biology. 2013;40(12):1298–309.
  15. 15. Zhang J, Liang S, Duan J, Wang J, Chen S, Cheng Z, et al. De novo assembly and Characterisation of the Transcriptome during seed development, and generation of genic-SSR markers in Peanut (Arachis hypogaea L.). BMC Genomics. 2012;13(1):90.
  16. 16. Chen X, Zhu W, Azam S, Li H, Zhu F, Li H, et al. Deep sequencing analysis of the transcriptomes of peanut aerial and subterranean young pods identifies candidate genes related to early embryo abortion. Plant Biotechnology Journal. 2013;11(1):115–27. pmid:23130888
  17. 17. Yin D, Wang Y, Zhang X, Li H, Lu X, Zhang J, et al. De Novo Assembly of the Peanut (Arachis hypogaea L.) Seed Transcriptome Revealed Candidate Unigenes for Oil Accumulation Pathways. PLoS ONE. 2013;8(9):e73767. pmid:24040062
  18. 18. Morgante C, Guimaraes P, Martins A, Araujo A, Leal-Bertioli S, Bertioli D, et al. Reference genes for quantitative reverse transcription-polymerase chain reaction expression studies in wild and cultivated peanut. BMC Research Notes. 2011;4(1):339.
  19. 19. Zhao S, Fernald RD. Comprehensive Algorithm for Quantitative Real-Time Polymerase Chain Reaction. Journal of computational biology: a journal of computational molecular cell biology. 2005;12(8):1047–64.
  20. 20. Pfaffl M, Tichopad A, Prgomet C, Neuvians T. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper—Excel-based tool using pair-wise correlations. Biotechnology Letters. 2004;26(6):509–15. pmid:15127793
  21. 21. Kim KH, Kang YJ, Kim DH, Yoon MY, Moon J-K, Kim MY, et al. RNA-Seq Analysis of a Soybean Near-Isogenic Line Carrying Bacterial Leaf Pustule-Resistant and -Susceptible Alleles. DNA Research. 2011;18(6):483–97. pmid:21987089
  22. 22. Hiremath PJ, Farmer A, Cannon SB, Woodward J, Kudapa H, Tuteja R, et al. Large-scale transcriptome analysis in chickpea (Cicer arietinum L.), an orphan legume crop of the semi-arid tropics of Asia and Africa. Plant Biotechnology Journal. 2011;9(8):922–31. pmid:21615673
  23. 23. Varshney RK, Chen W, Li Y, Bharti AK, Saxena RK, Schlueter JA, et al. Draft genome sequence of pigeonpea (Cajanus cajan), an orphan legume crop of resource-poor farmers. Nat Biotech. 2012;30(1):83–9.
  24. 24. Kalavacharla V, Liu Z, Meyers B, Thimmapuram J, Melmaiee K. Identification and analysis of common bean (Phaseolus vulgaris L.) transcriptomes by massively parallel pyrosequencing. BMC Plant Biology. 2011;11(1):135.
  25. 25. Proite K, Leal-Bertioli S, Bertioli D, Moretzsohn M, da Silva F, Martins N, et al. ESTs from a wild Arachis species for gene discovery and marker development. BMC Plant Biology. 2007;7(1):7.
  26. 26. Guimaraes P, Brasileiro A, Morgante C, Martins A, Pappas G, Silva O, et al. Global transcriptome analysis of two wild relatives of peanut under drought and fungi infection. BMC Genomics. 2012;13(1):387.
  27. 27. Wu N, Matand K, Wu H, Li B, Li Y, Zhang X, et al. De novo next-generation sequencing, assembling and annotation of Arachis hypogaea L. Spanish botanical type whole plant transcriptome. Theoretical and Applied Genetics. 2013;126(5):1145–9. pmid:23338522
  28. 28. Brasileiro AM, Morgante C, Araujo AG, Leal-Bertioli SM, Silva A, Martins AQ, et al. Transcriptome Profiling of Wild Arachis from Water-Limited Environments Uncovers Drought Tolerance Candidate Genes. Plant Mol Biol Rep. 2015:1–17.
  29. 29. Schaff JE, Nielsen DM, Smith CP, Scholl EH, Bird DM. Comprehensive Transcriptome Profiling in Tomato Reveals a Role for Glycosyltransferase in Mi-Mediated Nematode Resistance. Plant Physiology. 2007;144(2):1079–92. pmid:17434994
  30. 30. Guimarães P, Brasileiro A, Proite K, de Araújo A, Leal-Bertioli S, Pic-Taylor A, et al. A Study of Gene Expression in the Nematode Resistant Wild Peanut Relative, Arachis stenosperma, in Response to Challenge with Meloidogyne arenaria. Tropical Plant Biol. 2010:1–10.
  31. 31. Chen N, Yang Q, Pan L, Chi X, Chen M, Hu D, et al. Identification of 30 MYB transcription factor genes and analysis of their expression during abiotic stress in peanut (Arachis hypogaea L.). Gene. 2014;533(1):332–45. pmid:24013078
  32. 32. Hong L, Hu B, Liu X, He CY, Yao Y, Li XL, et al. Molecular cloning and expression analysis of a new stress-related AREB gene from Arachis hypogaea. Biol Plant. 2013;57(1):56–62.
  33. 33. Udvardi MK, Kakar K, Wandrey M, Montanari O, Murray J, Andriankaja A, et al. Legume Transcription Factors: Global Regulators of Plant Development and Response to the Environment. Plant Physiology. 2007;144(2):538–49. pmid:17556517
  34. 34. Tirumalaraju SV, Jain M, Gallo M. Differential gene expression in roots of nematode-resistant and -susceptible peanut (Arachis hypogaea) cultivars in response to early stages of peanut root-knot nematode (Meloidogyne arenaria) parasitization. Journal of Plant Physiology. 2011;168(5):481–92. pmid:20863592
  35. 35. Portillo M, Cabrera J, Lindsey K, Topping J, Andrés MF, Emiliozzi M, et al. Distinct and conserved transcriptomic changes during nematode-induced giant cell development in tomato compared with Arabidopsis: a functional role for gene repression. New Phytologist. 2013;197(4):1276–90. pmid:23373862
  36. 36. Das S, Ehlers J, Close T, Roberts P. Transcriptional profiling of root-knot nematode induced feeding sites in cowpea (Vigna unguiculata L. Walp.) using a soybean genome array. BMC Genomics. 2010;11(1):480.
  37. 37. Jammes F, Lecomte P, de Almeida-Engler J, Bitton F, Martin-Magniette M-L, Renou JP, et al. Genome-wide expression profiling of the host response to root-knot nematode infection in Arabidopsis. The Plant Journal. 2005;44(3):447–58. pmid:16236154
  38. 38. Klink V, Overall C, Alkharouf N, MacDonald M, Matthews B. A time-course comparative microarray analysis of an incompatible and compatible response by Glycine max (soybean) to Heterodera glycines (soybean cyst nematode) infection. Planta. 2007;226(6):1423–47. pmid:17653570
  39. 39. Fuller VL, Lilley CJ, Urwin PE. Nematode resistance. New Phytologist. 2008;180(1):27–44. pmid:18564304
  40. 40. Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, et al. mapman: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. The Plant Journal. 2004;37(6):914–39. pmid:14996223
  41. 41. Nandi B, Kundu K, Banerjee N, Babu SPS. Salicylic acid-induced suppression of Meloidogyne incognita infestation of okra and cowpea. Nematology. 2003;5(5):747–52.
  42. 42. Branch C, Hwang C-F, Navarre DA, Williamson VM. Salicylic Acid Is Part of the Mi-1-Mediated Defense Response to Root-Knot Nematode in Tomato. Molecular Plant-Microbe Interactions. 2004;17(4):351–6. pmid:15077667
  43. 43. Melillo MT, Leonetti P, Bongiovanni M, Castagnone-Sereno P, Bleve-Zacheo T. Modulation of reactive oxygen species activities and H2O2 accumulation during compatible and incompatible tomato–root-knot nematode interactions. New Phytologist. 2006;170(3):501–12. pmid:16626472
  44. 44. Williamson VM, Kumar A. Nematode resistance in plants: the battle underground. Trends in Genetics. 2006;22(7):396–403. pmid:16723170
  45. 45. Djian-Caporalino C, Palloix A, Fazari A, Marteu N, Barbary A, Abad P, et al. Pyramiding, alternating or mixing: comparative performances of deployment strategies of nematode resistance genes to promote plant resistance efficiency and durability. BMC Plant Biology. 2014;14(1):53.
  46. 46. van der Hoorn RAL, Kamoun S. From Guard to Decoy: A New Model for Perception of Plant Pathogen Effectors. The Plant Cell. 2008;20(8):2009–17. pmid:18723576
  47. 47. Zinovieva SV, Vasyukova NI, Udalova ZV, Gerasimova NG, Ozeretskovskaya OL. Involvement of salicylic acid in induction of nematode resistance in plants. Biol Bull Russ Acad Sci. 2011;38(5):453–8.
  48. 48. Charron J-B, Ouellet F, Houde M, Sarhan F. The plant Apolipoprotein D ortholog protects Arabidopsis against oxidative stress. BMC Plant Biology. 2008;8(1):86.
  49. 49. Kiselev KV, Dubrovina AS, Isaeva GA, Zhuravlev YN. The effect of salicylic acid on phenylalanine ammonia-lyase and stilbene synthase gene expression in Vitis amurensis Cell Culture. Russ J Plant Physiol. 2010;57(3):415–21.
  50. 50. Chang X, Heene E, Qiao F, Nick P. The Phytoalexin Resveratrol Regulates the Initiation of Hypersensitive Cell Death in Vitis Cell. PLoS ONE. 2011;6(10):e26405. pmid:22053190
  51. 51. Lopes R, Silveira D, Gimenes M, Vasconcelos P, Alves Rd, Silva J, et al. Characterization of resveratrol content in ten wild species of section Arachis, genus Arachis. Genetic Resources and Crop Evolution. 2013;60(8):2219–26.
  52. 52. Gan Z, Yang J, Tao N, Liang L, Mi Q, Li J, et al. Cloning of the gene Lecanicillium psalliotae chitinase Lpchi1 and identification of its potential role in the biocontrol of root-knot nematode Meloidogyne incognita. Appl Microbiol Biotechnol. 2007;76(6):1309–17. pmid:17665191
  53. 53. Isayenkov S, Mrosk C, Stenzel I, Strack D, Hause B. Suppression of Allene Oxide Cyclase in Hairy Roots of Medicago truncatula Reduces Jasmonate Levels and the Degree of Mycorrhization with Glomus intraradices. Plant Physiology. 2005;139(3):1401–10. pmid:16244141
  54. 54. Soriano I, Riley I, Potter M, Bowers W. Phytoecdysteroids: A Novel Defense Against Plant-Parasitic Nematodes. J Chem Ecol. 2004;30(10):1885–99. pmid:15609826
  55. 55. Camera SL, Balagué C, Göbel C, Geoffroy P, Legrand M, Feussner I, et al. The Arabidopsis Patatin-Like Protein 2 (PLP2) Plays an Essential Role in Cell Death Execution and Differentially Affects Biosynthesis of Oxylipins and Resistance to Pathogens. Molecular Plant-Microbe Interactions. 2009;22(4):469–81. pmid:19271961
  56. 56. Riemann M, Haga K, Shimizu T, Okada K, Ando S, Mochizuki S, et al. Identification of rice Allene Oxide Cyclase mutants and the function of jasmonate for defence against Magnaporthe oryzae. The Plant Journal. 2013;74(2):226–38. pmid:23347338
  57. 57. Doyle EA, Lambert KN. Meloidogyne javanica Chorismate Mutase 1 Alters Plant Cell Development. Molecular Plant-Microbe Interactions. 2003;16(2):123–31. pmid:12575746
  58. 58. Ding X, Cao Y, Huang L, Zhao J, Xu C, Li X, et al. Activation of the Indole-3-Acetic Acid–Amido Synthetase GH3-8 Suppresses Expansin Expression and Promotes Salicylate- and Jasmonate-Independent Basal Immunity in Rice. The Plant Cell Online. 2008;20(1):228–40.
  59. 59. Fu J, Liu H, Li Y, Yu H, Li X, Xiao J, et al. Manipulating Broad-Spectrum Disease Resistance by Suppressing Pathogen-Induced Auxin Accumulation in Rice. Plant Physiology. 2011;155(1):589–602. pmid:21071600
  60. 60. Wang F, Vandepoele K, Van Lijsebettens M. Tetraspanin genes in plants. Plant Science. 2012;190(0):9–15.
  61. 61. Zhu W, Wei W, Fu Y, Cheng J, Xie J, Li G, et al. A Secretory Protein of Necrotrophic Fungus Sclerotinia sclerotiorum That Suppresses Host Resistance. PLoS ONE. 2013;8(1):e53901. pmid:23342034
  62. 62. Knepper C, Savory EA, Day B. Arabidopsis NDR1 Is an Integrin-Like Protein with a Role in Fluid Loss and Plasma Membrane-Cell Wall Adhesion. Plant Physiology. 2011;156(1):286–300. pmid:21398259
  63. 63. Schröder R, Atkinson RG, Redgwell RJ. Re-interpreting the role of endo-β-mannanases as mannan endotransglycosylase/hydrolases in the plant cell wall. Annals of Botany. 2009;104(2):197–204. pmid:19454593
  64. 64. Xu X, Zhang Y, Meng Q, Meng K, Zhang W, Zhou X, et al. Overexpression of a Fungal β-Mannanase from Bispora sp. MEY-1 in Maize Seeds and Enzyme Characterization. PLoS ONE. 2013;8(2):e56146. pmid:23409143
  65. 65. Davies KG, Curtis RHC. Cuticle Surface Coat of Plant-Parasitic Nematodes. Annual Review of Phytopathology. 2011;49(1):135–56.
  66. 66. Singh KB, Foley RC, Oñate-Sánchez L. Transcription factors in plant defense and stress responses. Current Opinion in Plant Biology. 2002;5(5):430–6. pmid:12183182
  67. 67. Uehara T, Sugiyama S, Matsuura H, Arie T, Masuta C. Resistant and Susceptible Responses in Tomato to Cyst Nematode are Differentially Regulated by Salicylic Acid. Plant and Cell Physiology. 2010;51(9):1524–36. pmid:20660227
  68. 68. Hwang IS, Choi DS, Kim NH, Kim DS, Hwang BK. The pepper cysteine/histidine-rich DC1 domain protein CaDC1 binds both RNA and DNA and is required for plant cell death and defense response. New Phytologist. 2014;201(2):518–30. pmid:24117868
  69. 69. Zeng L-R, Qu S, Bordeos A, Yang C, Baraoidan M, Yan H, et al. Spotted leaf11, a Negative Regulator of Plant Cell Death and Defense, Encodes a U-Box/Armadillo Repeat Protein Endowed with E3 Ubiquitin Ligase Activity. The Plant Cell. 2004;16(10):2795–808. pmid:15377756
  70. 70. Chu Y, Guimarães LA, Wu CL, Timper P, Holbrook CC, Ozias-Akins P. A Technique to Study Meloidogyne arenaria Resistance in Agrobacterium rhizogenes-Transformed Peanut. Plant Disease. 2014;98(10):1292–9.
  71. 71. Zinovieva SV, Vasyukova NI, Udalova ZV, Gerasimova NG. The Participation of salicylic and jasmonic acids in genetic and induced resistance of tomato to Meloidogyne incognita (Kofoid and White, 1919). Biol Bull Russ Acad Sci. 2013;40(3):297–303.
  72. 72. Zambounis A, Kalamaki M, Tani E, Paplomatas E, Tsaftaris A. Expression Analysis of Defense-Related Genes in Cotton (Gossypium hirsutum) after Fusarium oxysporum f. sp. vasinfectum Infection and Following Chemical Elicitation using a Salicylic Acid Analog and Methyl Jasmonate. Plant Mol Biol Rep. 2012;30(1):225–34.
  73. 73. Thaler JS, Humphrey PT, Whiteman NK. Evolution of jasmonate and salicylate signal crosstalk. Trends in Plant Science.17(5):260–70. pmid:22498450