- Research article
- Open Access
Transcriptomic profiling of Solanum peruvianum LA3858 revealed a Mi-3-mediated hypersensitive response to Meloidogyne incognita
BMC Genomics volume 21, Article number: 250 (2020)
The Mi-1 gene was the first identified and cloned gene that provides resistance to root-knot nematodes (RKNs) in cultivated tomato. However, owing to its temperature sensitivity, this gene does not meet the need for breeding disease-resistant plants that grow under high temperature. In this study, Mi-3 was isolated from the wild species PI 126443 (LA3858) and was shown to display heat-stable resistance to RKNs. However, the mechanism that regulates this resistance remains unknown.
In this study, 4760, 1024 and 137 differentially expressed genes (DEGs) were enriched on the basis of pairwise comparisons (34 °C vs. 25 °C) at 0 (before inoculation), 3 and 6 days post-inoculation (dpi), respectively. A total of 7035 DEGs were identified from line LA3858 in the respective groups under the different soil temperature treatments. At 3 dpi, most DEGs were enriched in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to plant biotic responses, such as “plant-pathogen interaction” and “plant hormone signal transduction”. Significantly enriched DEGs were found to encode key proteins such as R proteins and heat-shock proteins (HSPs). Moreover, other DEGs were found to participate in Ca2+ signal transduction; the production of ROS; DEGs encoding transcription factors (TFs) from the bHLH, TGA, ERF, heat-shock transcription factor (HSF) and WRKY families were highly expressed, which contribute to be involved into the formation of phytohormones, such as salicylic acid (SA), jasmonic acid (JA) and ethylene (ET), the expression of most was upregulated at 3 dpi at the 25 °C soil temperature compared with the 34 °C soil temperature.
Taken together, the results of our study revealed reliable candidate genes from wild materials LA3858, that are related to Mi-3-mediate resistance to Meloidogyne incognita. A large number of vital pathways and DEGs were expressed specifically in accession LA3858 grown at 34 °C and 25 °C soil temperatures at 3 dpi. Upon infection by RKNs, pattern-recognition receptors (PRRs) specifically recognized conserved pathogen-associated molecular patterns (PAMPs) as a result of pathogen-triggered immunity (PTI), and the downstream defensive signal transduction pathway was likely activated through Ca2+ signal channels. The expression of various TFs was induced to synthesize phytohormones and activate R proteins related to resistance, resulting in the development of effector-triggered immunity (ETI). Last, a hypersensitive response in the roots occurred, which was probably induced by the accumulation of ROS.
Members of the genus Meloidogyne, which represent major pests worldwide, have a substantial negative influence on the development of various plant species . The J2 stage (the second stage of juveniles) is the main infection stage of RKNs; during this stage, giant cells (GCs) form in the roots of plants and absorb nutrients from the roots for their own growth and reproduction . Once established, J2-stage RKNs undergo three successive molts to become adult females . The development of genetic resistance is an effective method for reducing yield losses caused by RKN infection.
ROS have a substantial influence on reactions to biotic and abiotic stresses. Phytohormones such as JA and SA not only function to regulate plant growth but also are involved in plant defense signaling pathways . Genes that encode nucleotide-binding site–leucine-rich repeats (NBS-LRRs) are the predominant members of the R gene family, accounting for approximately 80% of the more than 140 cloned R genes . In potato, Gpa2 is an R gene that encodes the GPA2 protein, which depends on the recognition specificity afforded by both amino acid 187 and the Gpa2 LRR domain and provides resistance against two Globodera pallida nematode populations (D383 and D372) . Additionally, the Rhg1 gene from soybean, the Me3 and Me4 genes from pepper, and the Mi-1 gene from tomato provide resistance against specific strains of nematodes.
The Mi-1 gene, which has been mapped to the short arm of chromosome 6 and whose product contains a putative coiled-coil domain preceding the nucleotide-binding site (NBS) , is the major effective R gene against RKN species in tomato (Solanum lycopersicum) . The Mi genomic region contains three homologous genes, which are referred to as Mi-1.1, Mi-1.2 and Mi-1.3. Of these genes, only Mi-1.2 provides resistance against RKNs, including Meloidogyne incognita, Meloidogyne javanica and Meloidogyne arenaria . Mi-1 is an effective genetic resource for use against RKNs; however, Mi-1-mediated resistance is inactive at soil temperatures greater than 28 °C . Thus, additional heat-stable genes that provide resistance against RKNs must be identified to overcome this barrier. Recently, Mi-3 has attracted increased attention because of its heat-stable characteristics. In Solanum peruvianum LA3858, Mi-3 is located on the short arm of chromosome 12 and contains a 600-kb contig between the Mi-3-flanking markers TG180 and NR18, corresponding to a genetic distance of approximately 7.2 cM . Although it originates from a wild species, Mi-3, whose product functions effectively when temperatures reach 32 °C, could be a valuable source of resistance for cultivated tomato. However, self-incompatibility and distant hybridization incompatibility are the primary barriers preventing this gene from being finely mapped and cloned .
This study was designed to investigate via RNA sequencing (RNA-seq) Mi-3-mediated resistance in plants grown at two different soil temperatures (34 °C and 25 °C) following inoculation with M. incognita. Our goal was to identify key DEGs at the transcriptional level from the perspectives of PTI and ETI. We analyzed the plant defense response pathways with which these DEGs were significantly involved. Last, the Mi-3 gene-mediated disease response was characterized via a functional analysis of the proteins encoded by the DEGs in combination with analyses of the differences in DEG expression trends between the resistant (25 °C) and susceptible (34 °C) lines.
RKN disease evaluation under different temperature treatments
According to the infection results, the Moneymaker line displayed relatively consistent susceptibility to M. incognita at both soil temperatures (Fig. 1a). According to the resistance index results, all the seedlings were rated as S and HS (Additional file 1). In accession LA3858, all 8 plants at the soil temperature of 25 °C were immune to infection (rated as HR). At the soil temperature of 32 °C, very small galls had developed on the roots; however, all plants were rated as HR. When the plants were subjected to a soil temperature of 34 °C, gall formation on the roots was enhanced, and the resistance was rated as S and HS. According to a one-way ANOVA, although the number of egg masses from the Moneymaker line was not significant (P > 0.05), the root gall numbers on accession LA3858 plants in the 34 °C soil temperature treatment were obviously greater than those on plants in the other two soil treatments (P ≤ 0.05), which indicated that any resistance to RKNs of LA3858 will be completely absent (HS) at 34 °C (Fig. 1b).
Illumina sequencing and alignment to the reference genome
RNA-seq data were generated from 18 samples of the wild species LA3858 at different stages following inoculation (0, 3 and 6 dpi). After sequencing a total of 89 billion fragments of clean reads, we obtained approximately 49 M reads for each sample after aligning them to the reference genome (SL 3.0, ftp://ftp.sgn.cornell.edu/genomes/Solanum_lycopersicum/assembly/build_3.00/) via TopHat2 (version 22.214.171.124). Approximately 46 M clean reads per sample were obtained for subsequent analysis after the rRNA sequences, adapter sequences and low-quality reads were filtered and removed. The Q20 values (base calling error probability = 99%) of the 18 samples were greater than 98%. The expression profiles of 35,768 genes were ultimately used for further analysis (Fig. 2).
DEGs observed at 25 °C and 34 °C soil temperatures
The DEG analysis revealed that 5921 DEGs were enriched (P ≤ 0.05) in the three groups; the expression of 2802, 904 and 100 of these genes was upregulated, and that of 1958, 120 and 37 DEGs was downregulated in 34 °C vs. 25 °C at 0 (before inoculation), 3 and 6 dpi (HS0 vs. HR0, HS3 vs. HR3 and HS6 vs. HR6), respectively (Fig. 3). The overlapping genes among these groups are shown in Fig. 4.
GO enrichment analysis revealed the enrichment of the DEGs with respect to the following three categories: cellular component, molecular function and biological process. Most of the DEGs were enriched (P ≤ 0.05) in the cellular component category and were involved in the “cell”, “cell part”, “membrane”, and “organelle” terms. The significantly enriched terms (P ≤ 0.05) in the molecular function category were “binding”, “catalytic activity” and “transporter activity”, and the enriched terms (P ≤ 0.05) in the biological process category included “biological regulation”, “cellular process”, “metabolic process”, and “response to stimulus and signaling”, which were related to disease resistance  (Figure S1). Notably, in the HS3 vs. HR3 comparison, many key KEGG pathways related to biotic stress were significantly enriched (P ≤ 0.05), including “plant-pathogen interaction” (16 DEGs), “plant hormone signal transduction” (16 DEGs), and “brassinosteroid biosynthesis” (3) (Fig. 5). However, in the HS0 vs. HR0 and HS6 vs. HR6 comparisons, enrichment of these key pathways was not obvious.
Trends of DEGs in the same lines (25 °C and 34 °C)
To analyze the DEGs in the same lines at different time points (0, 3 and 6 dpi), trend analysis was used to discover DEG expression patterns. When LA3858 was subjected to the 25 °C soil temperature treatment, most DEGs were enriched in profiles 1 and 6 (P ≤ 0.05) (Fig. 6a). The KEGG analysis revealed a unique increasing expression trend (Fig. 6b), as profile 6 exhibited a high enrichment of DEGs. Several key pathways related to disease resistance were also significantly enriched, such as “plant hormone signal transduction” (31 DEGs) and “plant-pathogen interaction” (25 DEGs) (Fig. 6c). In 34 °C line, most DEGs were significantly enriched in profiles 2 and 4 (P ≤ 0.05) (Fig. 7a). Similarly, because DEGs enriched in profile 2 exhibited a specific downregulation expression trend (Fig. 7b), profile 2 received increased attention. In addition to the “plant-pathogen interaction” (9 DEGs) pathway, the “phenylpropanoid biosynthesis” (15 DEGs) and “flavonoid biosynthesis” (4 DEGs) pathways, whose metabolites often play an active role in regulating the response to biotic stimulus, were also enriched in profile 2 (Fig. 7c).
Gene expression under different temperature treatments at 3 dpi
At the important time period of 3 dpi, in terms of the key pathways, “plant-pathogen interaction” and “plant hormone signal transduction” were the primary ones identified.
In the “plant-pathogen interaction” pathway, the expression levels of 2 genes that encode calcium-dependent protein kinases (CDPKs), 3 genes that encode respiratory burst oxidases (RBOHs), and 3 genes that encode LRR receptor-like serine/threonine-protein kinases (FLS2s) were upregulated. Additionally, the expression levels of 2 genes that encode disease resistance proteins (RPs; RPM1 and PRS2) and 2 genes that encode HSPs (HSP90s) were also upregulated in the plants grown at 25 °C soil temperature. The other genes encode enhanced disease susceptibility 1 protein (EDS1), which is involved in programmed cell death (PCD), and WRKY TFs. In the “plant hormone signal transduction” pathway, 16 genes were enriched significantly (P ≤ 0.05), with 11 involved in the auxin response, including those encoding auxin-responsive proteins (IAAs), auxin response factors (ARFs), auxin-responsive GH3 family members and SAUR proteins; the expression of 10 of these genes was upregulated. Other genes encode a serine/threonine-protein kinase (SRK2), which is involved in the abscisic acid pathway, a brassinosteroid (BR) signaling kinase (BSK), which is involved in BR biosynthesis, and 2 TGA TFs, which are involved with SA (Table 1).
The trend analysis revealed that, in 25 °C line, DEGs were enriched significantly (P ≤ 0.05) in profile 6 class. In the “plant-pathogen interaction” pathway; in addition to the above mentioned protein-coding genes, the genes encoding RPM1-interacting protein (RIN4) and pathogen-induced protein kinase (PIK1) were highly enriched. Compared with the 34 °C vs. 25 °C, a large number of genes were involved in other key processes, such as those associated with cytokinin (CK), gibberellin (GA), ET and JA in the “plant hormone signal transduction” pathway (Table 2). In 34 °C line, DEGs in profile 2 were highly enriched (P ≤ 0.05). The expression levels of 9 DEGs that encoded FLS2, RBOH, RPS2, HSP90 and WRKY, which involved into pathway “plant-pathogen interaction”, decreased at 3 dpi (Table 3).
Analysis of the hub genes from the coexpression network during incompatible interactions
The genes related to the regulation of the resistance mechanism of LA3858 at different soil temperatures were further investigated. After clustering the module genes according to the standards mentioned above, we selected a total of 17,184 genes for the construction of a scale-free coexpression network. Thirteen coexpression modules were constructed (Fig. 8). Of these modules, a total of 6 (i.e., bisque4, brown, darkmagenta, darkorange, darkorange 2 and pink) were selected because of the specificity of the expression pattern at each time period. Among these 6 modules, which had been subjected to KEGG analysis, the darkorange and pink modules attracted our attention because the expression of the DEGs in both modules tended increase at 3 dpi (Fig. 9), and the “plant-pathogen interaction” pathway was also significantly enriched (P ≤ 0.05) in both modules (Figure S2). With respect to genes, Pearson correlation coefficients ≥0.8 were filtered to establish DEGs coexpression network to reveal hub genes whose expression is induced during Mi-3-mediated resistance (Fig. 10). In the darkorange module, 6 hub genes encode histidine decarboxylase (HDC), 3 genes encode calcium-binding protein (CML) and one encodes an ET-responsive transcription factor (ERF). In the pink module, CDPK and RBOH were also found to be encoded by hub genes. These results suggested that Ca2+ channels may play a key role in signal transduction during nematode infection , at the same time, the level of ROS may also contribute to the regulation of resistance, which is consistent with our previous results obtained from incompatible interactions  (Additional file 2).
Key TFs that function during infection in Mi-3-mediated resistance
RNA-seq analysis revealed a total of 832 TFs. Five main TF categories related to disease resistance in plants were identified: WRKYs, heat-shock transcription factors (HSFs), bHLHs, TGAs, and ERFs . On the basis of the significant expression levels (P ≤ 0.05), we compared the detailed expression trends of the DEGs that encode TFs from the above 5 families between 34 °C and 25 °C (Additional file 3). There were more enriched DEGs that encode TFs from these families at 0 days (before inoculation) than at the other two stages (3 dpi and 6 dpi). This was most likely due to the increase in soil temperature, and the DEGs that encode TFs involved in the response to abiotic stress were significantly enriched. At 3 dpi, the expression of nearly all the significantly enriched DEGs of the TFs was upregulated, and among those that encode TFs, the number of DEGs encoding WRKY TFs was the largest. Almost no DEGs encoding TFs were enriched at 6 dpi, which indicated again that Mi-3-mediated resistance occurs mainly during the early stages of nematode infection (Fig. 11).
Validation of RNA-seq data via RT-qPCR
To validate the RNA-seq results, 15 DEGs with different expression patterns were selected for qRT-PCR analysis via gene-specific primers (Additional file 4). All of these candidate genes were screened from the 34 °C vs. 25 °C comparison at 3 dpi, and there was a strong positive correlation coefficient (R2 = 0.9404) between the qRT-PCR results and the RNA-seq data, suggesting that the RNA-seq data were reproducible and reliable (Fig. 12).
In recent years, the Mi-3 gene has attracted increased attention because of its heat-stable characteristics, as this gene remains active when soil temperatures reach 32 °C. In the present study, the results of disease resistance tests demonstrated that accession LA3858 was susceptible to M. incognita when the plants were growing at a soil temperature of 34 °C compared with a normal soil temperature (25 °C). We further elucidated the regulatory networks and models of the Mi-3 gene via RNA-seq and a WGCNA (Figure S3). And RNA-seq was used to verify the transcriptomic profiles of LA3858 in response to M. incognita, and the reliability of the RNA-seq data was verified by the significant positive correlation detected between the qRT-PCR data and the RNA-seq data (Additional file 5). GO analysis revealed that most of the DEGs were enriched in the “cell wall”, “external encapsulating structure”, “cell division” and other key terms at 3 dpi between 34 °C and 25 °C, which may be related to GC formation . DEGs were significantly enriched in each of the identified KEGG pathways, which were shown to be associated with plant immunity in a previous study. The initial stage of infection was a key period for RKN penetration and migration in tomato ; thus, on the basis of this timing, 3 dpi was selected as the primary stage for analysis.
The PTI defense system is formed by Mi-3-mediated disease resistance
The resistance formed in response to external stimuli in plants is referred to as the plant natural immune system and can be divided into two levels, the first of which is called PTI. Plant-surface PRRs specifically recognize the conserved molecular structure of PAMPs of the pathogen, eventually inducing PTI . flg22 is by far the most researched type of PAMP, which is specifically recognized by FLS2 . In this work, LA3858 plants growing at 25 °C soil temperature were resistant to M. incognita, and compared with those in plants growing at a 34 °C soil temperature, the expression of Solyc02g072480.3-, Solyc04g014650.3- and Solyc06g006020.2-encoded FLS2 was upregulated at 3 dpi. And Solyc06g006020.2 were having an increasing trend in 25 °C line (profile 6) at 3 dpi and 6 dpi; however, Solyc02g072480.3 and another DEG (Solyc02g072393.1) also encoding FLS2 tended to decrease in 34 °C line (profile 2) at both 3 and 6 dpi (Tables 1, 2 and 3). This most likely indicated that in plants with Mi-3 at 34 °C soil temperature, losing resistance was related to the obstruction of developing PTI. After PRRs recognize PAMP molecules, plants can induce rapid defense responses in a short period of time, which can involve activating mitogen-activated protein kinase (MAPK) signaling pathways, increasing ROS levels, and activating SA and JA signaling pathways .
Increases in Ca2+ signal channels and ROS in response to J2 inoculation
After PTI develops, Ca2+ signal channels are instantly activated by MAPKs, and both CMLs and CDPKs are important components of these channels [22, 23]. The expression of three DEGs encoding CPK2, CPK16 and CML41 was upregulated at 3 dpi in the 34 °C vs. 25 °C comparison. Previous research has suggested that CDPKs are closely related to ROS levels during pathogen damages. RBOHs, also known as plant NADPH oxidases, are involved in the production of ROS. RBOH N-terminal domains contain two EF-hands, which are involved in calcium binding . Via protein phosphorylation, RBOHs are activated to produce ROS . Consistent with this function, in the present study, RBOHs were encoded by 3 DEGs whose expression was upregulated at 25 °C compared with 34 °C at 3 dpi. Therein, Solyc06g068680.3, which encodes RBOHD, was also significantly enriched (P ≤ 0.05) in 25 °C line according to the results of the trend analysis (profile 6), and its expression was downregulated in 34 °C line (profile 2) at both 3 and 6 dpi (Tables 1, 2 and 3). It was also found that this DEG is a hub gene in the pink module and that its regulates the production of ROS together with that of Solyc12g005030.2, which encodes another CDPK (Additional file 2). Apparently, RBOHD-dependent ROS production may play a positive role during M. incognita infection. According to a previous study, early production of ROS could be a key biological barrier against disease progression in plant development . In the process of Mi-3-mediated resistance, we speculate that, owing to the rapid accumulation of ROS, M. incognita at the J2 stage cannot survive the relatively high oxygen concentration in the environment. Furthermore, cell damage and PCD resulting from ROS may be two reasons why nematode feeding sites (NFSs) could not be established, which prevented nutrient extraction and prevented RKNs from forming GCs .
Differential expression of WRKY TFs between 34 °C and 25 °C soil temperatures
When Ca2+ signal channels are instantaneously activated, in addition to ROS levels, the expression of downstream TFs is also regulated . The results indicated that the expression of many more TFs whose encoding genes were significantly enriched was upregulated at 3 dpi. WRKY TFs are extremely important in plant growth and disease resistance and are widely studied in tomato . The expression of six key genes encoding WRKY TFs was significantly upregulated at 25 °C compared with 34 °C at 3 dpi. Among these TF genes, the expression of Solyc04g051540.3, Solyc12g006170.2, Solyc10g009550.3, Solyc08g067340.3, Solyc01g095630.3 and Solyc03g095770.3, which encode WRKY13, 20, 30, 40, 41 and 70, respectively, differed by approximately 5.10-, 2.58-, 9.23-, 4.70-, 3.46- and 4.91-fold (Additional file 3). In addition to WKRY20, which was most likely related to drought tolerance , other five WRKYs all might contribute to plants disease defense. In a previous study, WRKY13 was shown to mediate disease resistance to RKNs by activating SA-dependent pathways and by suppressing JA-dependent pathways in rice . Overexpression of WRKY30 in rice increased resistance to rice sheath blight fungus by upregulating the expression of the JA synthesis-related genes LOX and AOS2 and a series of pathogenesis-related (PR) genes . The expression of GmWRKY40 was strongly induced in soybean following infection with Phytophthora sojae, and silencing GmWRKY40 was shown to increase susceptibility . WRKY41 is likely to be a key regulator in the cross talk of SA and JA pathways. The overexpression of WRKY41 leads to an increasing of PR5 expression, for enhancing the resistance to Pto wild-type . Other studies have shown that, in tomato, silencing of SlWRKY70 attenuates Mi-1-mediated resistance against RKNs . In short, WRKYs are often involved in various defense mechanisms associated with plant growth and stress responses. Despite their involvement in the regulation of expression of some key PR genes and R genes, many WRKYs may also contribute to JA- and SA-defense signaling processes [35, 36].
Phytohormone regulation during M. incognita infection
Phytohormones have simple molecular structures and low cellular levels . SA, JA and ET play crucial roles in resistance to biotrophic and necrotrophic pathogens, and TGAs, MYC2 and ERFs are key TFs involved in the biosynthesis of the above three phytohormones . In this research, Solyc10g009290.1 was found to encode MYC2, which, as a member of the basic helix-loop-helix (bHLH) TF family, controls JA-dependent responses and was highly expressed in 25 °C line (profile 6). 2 and 3 DEGs encoding TGA and ERF TFs, respectively, were highly expressed at 3 dpi between 34 °C and 25 °C. Moreover, the expression of several DEGs that encoded ethylene receptors (ETRs), EIN3-binding F-box proteins (EBF1/2), and ethylene-insensitive protein 3 (EIN3) was upregulated in 25 °C line (profile 6), which contributed to the biosynthesis of the ET (Table 2). In addition, Solyc03g093610.1, which encodes the ET-responsive TF (ERF1), was a hub gene in the darkorange module (Additional file 2). Previous reports have shown that obstacles to the synthesis or accumulation of SA result in consistent increases in the number of phenotypes associated with susceptibility to soybean cyst nematode (SCN) . JA can reduce the damage caused by RKNs, and a recent report revealed that, in JA-overexpressing transgenic tomato plants, proteinase inhibitor II (PI-II) translation levels gradually increased with time after inoculation with RKN; this increase was related to the production of JA, which acts at the initial site of infection, inhibiting nematode invasion and propagation . ERF1 integrates defense signals from the ET and JA pathways and induces the expression of downstream defense-related genes . In short, SA and JA defense signaling processes are important components of the first-layer defense mechanism of plants after the development of PTI . According to previous research, both signaling pathways are likely to be involved in Mi-3-mediated resistance to M. incognita.
Roles of induced disease resistance genes related to Mi-3
When a host plant recognizes a pathogen and produces an immune response, the pathogen will produce a substance—an effector—to inhibit the host’s recognition. To prevent additional damage, plants directly or indirectly recognize effectors via R proteins, and initiate a rapid and violent hypersensitive response at the infection site, which constitutes the second level of ETI . Most R genes encode proteins that include a unique NBS-LRR domain. Members of the NBS-LRR protein family generally participate in downstream signal transduction during plant-pathogen interactions .
Two genes (Solyc08g007250.2 and Solyc04g007090.1) were found to encode the disease RPs RPS2 and RPM1, which attracted our attention because of their significant increase (P ≤ 0.05) in expression at 25 °C compared with 34 °C at 3 dpi. It was suggested that both RPM1 and RPS2 are members of the NBS-LRR protein family. A previous report showed that these two R proteins recognize the AvrRpm1 type III effector avirulence protein involved in protecting plants against pathogens . The Solyc04g007090.1 gene encodes RPM1 and exhibited 2.58- and 2.71-FCs in expression levels in 25 °C line (profile 6) at 3 and 6 dpi, respectively; this gene was not obviously enriched in 34 °C line (profile 2). In addition, the gene Solyc02g037540.2, which encodes another RPS2, showed 2.57- and 2.71-fold increases in 25 °C line (profile 6) at 3 and 6 dpi, respectively. And showed downward trend in 34 °C line (profile 2). Therefore, RPM1 and RPS2 may be key proteins during infection. Additionally, Solyc09g059430.3 was found to encode RIN4, a negative regulator of basal defense responses, whose expression was upregulated in 25 °C line (profile 6) along with that of PRS2 and RPM1 (Tables 1, 2 and 3). The biological functions of these two R proteins are associated with RIN4 . RIN4 physically interacts with RPM1 and indirectly interacts with RPS2. Moreover, the RPM1 and RPS2 proteins work together to maintain RIN4 expression levels during pathogen infection . However, the biological mechanisms of these three proteins, which are regulated by Mi-3, in response to RKNs require further study.
HSPs may play a positive role against nematode infection
In most cases, disease resistance mediated by both RPS2 and RPM1 is associated with RAR1, SGT1 and HSPs. RAR1 interacts with the N-terminal half of HSP90, which contains an ATPase domain, and HSP90 specifically interacts with SGT1. In Arabidopsis, the interaction of these three proteins regulates the stability and function of RPM1, which is an HSP-client protein . In this study, the expression of Solyc05g010670.3 and Solyc06g036290.3, each of which encodes an HSP90 (Hsp83 and HSP83A), was upregulated at 25 °C compared with 34 °C at 3 dpi; therein, the expression of Solyc05g010670.3 obviously tended to decrease in 34 °C line (profile 2), and Solyc06g036290.3 had a significant increasing trend in 25 °C line (profile 6) at 3 and 6 dpi (Tables 1, 2 and 3). Additionally, another 2 DEGs were found to encode HSP90 (HSP83A, HSC80), and the expression both tended to increase in 25 °C line (profile 6). According to a previous study, genes that encode certain HSPs activated by HSFs perceive biotic stress signals via Ca2+ channels, which are highly expressed at the soybean Rhg1 locus and are involved in resistance to SCN . In addition, silencing of tomato HSP90 and SGT1 led to a reduction in Mi-1 protein levels, which reduced resistance to M. incognita , demonstrating that HSP90 contributes to the resistance process in plants and that DEGs encoding HSPs in our study might have a positive effect during infection of LA3858 by nematodes. Moreover, the production of HSP can also used as a standard for detecting levels of ROS in plants [48, 49].
In summary, this is the first report on the resistance and susceptibility of Mi-3 at the transcriptional level under different soil temperatures. A soil temperature of 32 °C is likely to be limiting for Mi-3 activity, which is consistent with the previous research. When the soil temperature was less than 32 °C, after M. incognita infected the roots, a rapid response (PTI) occurred, and downstream signal transduction was most likely triggered through Ca2+ channels via MAPKs and other signaling pathway components. Key defense-related TFs were subsequently triggered, such as HSFs, TGAs, ERFs, bHLHs and WRKYs, which activated disease RPs and downstream defense pathways, such as the SA, JA and ET pathways. Notably, HSPs likely associated with R proteins such as RPM1 and RPS2 in the development of ETI during infection. Last, the formation of a hypersensitive response in the roots was likely induced by ROS, resulting in cell damage and PCD. These actions prevented the establishment of root NFSs and might constitute a key mechanism for specific resistance to M. incognita.
Plant growth conditions and nematode assays
Seeds of S. peruvianum accession LA3858 were obtained from the Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, and seeds of the susceptible breeding line Moneymaker were obtained from the Northeast Agricultural University Tomato Research Institute. At 20–25 days after the seeds were sown in plots, when the seedlings reached a height of 12–15 cm (the second-leaf stage), they were transplanted into 10-cm-diameter plastic pots that each contained 50% turf soil and 50% roseite. The plants were maintained at a 24 ± 2 °C temperature, at 60% relative humidity and under a 14-h light/10-h dark photoperiod . The M. incognita strain was obtained from the Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences. We obtained infected roots and stripped the egg masses, which were cleaned with 1% NaOCl, after which the eggs were hatched in distilled water to obtain J2-stage individuals for inoculation .
Heat treatment for phenotypic identification and transcriptomic analysis
Approximately 60 plants of accession LA3858 and Moneymaker were grown at different soil temperatures (25 °C, 32 °C and 34 °C) for five days before inoculation. Afterward, 8 plants of both lines in each treatment were randomly picked for gall and egg mass counting 45–50 days after inoculation with 2000 J2-stage M. incognita nematodes per plant. For transcriptomic analysis of accession LA3858, high-temperature assays were performed in a greenhouse in which the soil temperature was greater than 32 °C (34 °C) five days before inoculation. In the other group, the plants were grown at normal temperature (25 °C). Afterward, the roots of plants from the two groups were collected at 0 days (before inoculation) and at 3 and 6 days post-inoculation (dpi) with 2000 J2-stage M. incognita nematodes per plant . At each time point, the roots of three replicates were collected. For each biological replicate, inoculated tissue was collected from 3 random plants, and the tissue samples were pooled together to obtain sufficient root tissue materials for RNA-seq and qRT-PCR verification. All root tissue samples were stored in liquid nitrogen and then transferred to − 80 °C conditions.
Disease score on root systems
The roots were removed from the soil, washed for several minutes, cleaned with a NaOCl solution, soaked in distilled water for 15 min, and then dyed with acid fuchsin. Finally, the galls and egg masses were counted under a microscope to assess the resistance. The root gall index (percentage of roots with galls) and the resistance index (rated according to the root gall index) were the standards used for evaluating disease resistance, and each was divided into 6 levels as follows: 0 – no galls, 1–1-10% of roots with galls, 2–10-20%, 3–20-50%, 4–50-80%, and 5–80-100% for the former and immune (I), highly resistant (HR), resistant (R), moderately resistant (MR), susceptible (S), and highly susceptible (HS) for the latter .
mRNA library construction and sequencing
In total, eighteen independent mRNA libraries from the roots of LA3858 plants in 2 treatments (25 °C soil temperature and 34 °C soil temperature) and at three time points (0 days, 3 days, 6 days) with three biological replicates per treatment were sequenced . After the total RNA was isolated via TRIzol™ reagent (No: 15596026, Thermo Fisher Scientific), oligo (dT) beads were used to enrich the eukaryotic mRNA, while prokaryotic mRNA was enriched by the removal of rRNA with a Ribo-Zero™ Magnetic Kit (Epicentre). The enriched mRNA was then fragmented into short fragments by fragmentation buffer and reverse transcribed into cDNA with random primers. Second-strand cDNA was synthesized with DNA polymerase I, RNase H, dNTPs and buffer. The cDNA fragments were then purified with a QiaQuick PCR Extraction Kit, end repaired, polyadenylated and then ligated to Illumina sequencing adapters. The appropriate sizes of the ligation products were selected via agarose gel electrophoresis, amplified by PCR, and sequenced by an Illumina HiSeq™ 2500 instrument by Gene Denovo Biotechnology Co. (Guangzhou, China). All the raw read data have been deposited into the NCBI Sequence Read Archive database (PRJNA494774).
Analysis of DEGs
To identify DEGs across the 18 samples, the edgeR package (http://www.rproject.org/) was used. We considered genes significantly differentially expressed when their log2|fold change (FC)| was > 1 and their false discovery rate (FDR) was < 0.05 in a particular comparison [55, 56]. Gene expression pattern analysis was then used to cluster genes with similar expression patterns for multiple samples (at least 3 at a specific time point, space, or treatment dose). To examine the expression pattern of the DEGs, the expression data for each sample (in the order of treatment) were normalized to 0, log2(v1/v0) transformed, and log2(v2/v0) transformed, after which they were clustered via Short Time-series Expression Miner (STEM) software . The clustered profiles with P-values ≤0.05 were considered significant. The DEGs in the profiles were then subjected to Gene Ontology (GO) functional analysis and KEGG pathway enrichment analysis. In this article, trend analysis was performed by clustering gene expression patterns of the characteristics of multiple continuous samples (at least 3). The gene sets that met certain biological characteristics were then selected from the clustering results. And trend analysis was together used to discuss DEGs that are in the same line and are associated with pairwise comparisons.
Weighted gene coexpression network analysis (WGCNA)
WGCNA is a systems biology method for describing correlation patterns among genes across multiple samples. Genes that express similar patterns can be clustered and analyzed for associations between modules and specific traits or phenotypes. In this study, coexpression networks were constructed via the WGCNA (v1.47) package in R . After the low-quality samples and samples that had an unstable effect on the results (genes that were not expressed in more than half of the samples and samples in which more than half of the genes were not expressed) were removed, the gene expression values were subjected to a WGCNA to construct coexpression modules via the automatic network construction function blockwiseModules with the default settings, with the exceptions that the power was 6, the TOMType was unsigned, the mergeCutHeight was 0.8, and the minModuleSize was 50 [59, 60].
To identify biologically significant modules, eigengenes were used to calculate the correlation coefficients with samples or sample traits. The intramodular connectivity (function softConnectivity) of each gene was calculated, and genes with a high connectivity tended to be hub genes (q. weighted < 0.01 as a cutoff) , which may act as key factors that regulate a large subset of genes to perform biological functions together. The networks were subsequently visualized by Cytoscape 3.3.0.
qRT-PCR and validation of RNA-seq results
qRT-PCR was performed to validate the accuracy of the RNA-seq results. First-strand cDNA was synthesized via a RevertAid First Strand cDNA Synthesis Kit (K1621) from Thermo Scientific. Each sample included three biological replicates. We used the NCBI BLAST program to design the primers used for the unigenes. Relative quantitative data were calculated according to the ΔΔCT method: normalization (ΔCT = CT (sample) – CT (GAPDH)); ΔΔCT = ΔCT (sample A) -ΔCT (sample B); relative quantification = 2-ΔΔCT .
Availability of data and materials
We have deposited our data in Sequence Read Archive (SRA) (http://www.ncbi.nlm.nih.gov/sra/), the accession number for our submissions are: PRJNA494774.
Calcium-dependent protein kinase
Differentially expressed genes
LRR receptor-like serine/threonine-protein kinases
Heat shock transcription factors
Heat shock proteins
Nicotinamide adenine dinucleotide phosphate
Nucleotide-binding site-leucine-rich repeat
Nematode feeding sites
Pathogen-associated molecular patterns
Programmed cell death
Respiratory burst oxidases
Reactive oxygen species
Soybean cyst nematode
Trudgill DL. Origins of root-knot nematodes (Meloidogyne spp., Nematoda) in relation to their cultural control [J]. Phytoparasitica. 1995;23(3):191–4.
Schaff JE, Nielsen DM, Smith CP, et al. Comprehensive transcriptome profiling in tomato reveals a role for glycosyltransferase in Mi-mediated nematode resistance [J]. Plant Physiol. 2007;144(2):1079–92..
Das S, Ehlers JD, Close TJ, et al. Transcriptional profiling of root-knot nematode induced feeding sites in cowpea (Vigna unguiculata L. Walp.) using a soybean genome array [J]. BMC Genomics. 2010;11(1):480.
Meldau S, Ullman-Zeunert L, Govind G, et al. MAPK-dependent JA and SA signalling in Nicotiana attenuata affects plant growth and fitness during competition with conspecifics [J]. BMC Plant Biol. 2012;12(1):213.
Kochetov AV, Glagoleva AY, Strygina KV, et al. Differential expression of NBS-LRR-encoding genes in the root transcriptomes of two Solanum phureja genotypes with contrasting resistance to Globodera rostochiensis [J]. BMC Plant Biol. 2017;17(2):251.
Sacco MA, Koropacka K, Grenier E, et al. The cyst nematode SPRYSEC protein RBP-1 elicits Gpa2-and RanGAP2-dependent plant cell death [J]. PLoS Pathog. 2009;5(8): e1000564.
Seah S, Yaghoobi J, Rossi M, et al. The nematode-resistance gene, Mi-1, is associated with an inverted chromosomal segment in susceptible compared to resistant tomato [J]. Theor Appl Genet. 2004;108(8):1635–42.
Devran Z, Göknur A, Mesci L. Development of molecular markers for the Mi-1 gene in tomato using the KASP genotyping assay [J]. Hortic Environ Biotechnol. 2016;57(2):156–60.
Reddy YS, Sellaperumal C, Prasanna HC, et al. Screening of tomato genotypes against root-knot nematode and validation of Mi 1 gene linked markers [J]. Proc Natl Acad Sci India Sect B. 2018;88(1):65–72.
Jablonska B, Ammiraju JSS, Bhattarai KK, et al. The Mi-9 gene from Solanum arcanum conferring heat-stable resistance to root-knot nematodes is a homolog of Mi-1[J]. Plant Physiol. 2007;143(2):1044–54.
Yaghoobi J, Yates JL, Williamson VM. Fine mapping of the nematode resistance gene Mi-3 in Solanum peruvianum and construction of a S. lycopersicum DNA contig spanning the locus [J]. Mol Gen Genomics. 2005;274(1):60–9.
Yaghoobi J, Kaloshian I, Wen Y, et al. Mapping a new nematode resistance locus in Lycopersicon peruvianum [J]. Theor Appl Genet. 1995;91(3):457–64.
Maboreke HR, Feldhahn L, Bönn M, et al. Transcriptome analysis in oak uncovers a strong impact of endogenous rhythmic growth on the interaction with plant-parasitic nematodes [J]. BMC Genomics. 2016;17(1):627.
Bali S, Vining K, Gleason C, et al. Transcriptome profiling of resistance response to Meloidogyne chitwoodi introgressed from wild species Solanum bulbocastanum into cultivated potato [J]. BMC Genomics. 2019;20(1):907.
Melillo MT, Leonetti P, Leone A, et al. ROS and NO production in compatible and incompatible tomato-Meloidogyne incognita interactions [J]. Eur J Plant Pathol. 2011;130(4):489–502.
Depuydt S, Hardtke CS. Hormone signalling crosstalk in plant growth regulation [J]. Curr Biol. 2011;21(9):R365–73.
Wang X, Cheng C, Zhang K, et al. Comparative transcriptomics reveals suppressed expression of genes related to auxin and the cell cycle contributes to the resistance of cucumber against Meloidogyne incognita [J]. BMC Genomics. 2018;19(1):583.
Shukla N, Yadav R, Kaur P, et al. Transcriptome analysis of root-knot nematode (Meloidogyne incognita)-infected tomato (Solanum lycopersicum) roots reveals complex gene expression profiles and metabolic networks of both host and nematode during susceptible and resistance responses [J]. Mol Plant Pathol. 2018;19(3):615–33.
Hong CY, Zheng JL, Chen TY, et al. PFLP-intensified disease resistance against bacterial soft rot through the MAPK pathway in PAMP-triggered immunity [J]. Phytopathology. 2018;108(12):1467–74.
Smith JM, Heese A. Rapid bioassay to measure early reactive oxygen species production in Arabidopsis leave tissue in response to living Pseudomonas syringae [J]. Plant Methods. 2014;10(1):6.
Dong S, Yin W, Kong G, et al. Phytophthora sojae avirulence effector Avr3b is a secreted NADH and ADP-ribose pyrophosphorylase that modulates plant immunity [J]. PLoS Pathog. 2011;7(11): e1002353.
Zhang K, Han YT, Zhao FL, et al. Genome-wide identification and expression analysis of the CDPK gene family in grape, Vitis spp [J]. BMC Plant Biol. 2015;15(1):164.
Davies LJ, Brown CR, Elling AA. Calcium is involved in the R Mc1 (blb)-mediated hypersensitive response against Meloidogyne chitwoodi in potato [J]. Plant Cell Rep. 2015;34(1):167–77.
Liu Y, He C. Regulation of plant reactive oxygen species (ROS) in stress responses: learning from AtRBOHD [J]. Plant Cell Rep. 2016;35(5):995–1007.
Zhang JZ, Zhao K, Ai XY, et al. Involvements of PCD and changes in gene expression profile during self-pruning of spring shoots in sweet orange (Citrus sinensis)[J]. BMC Genomics. 2014;15(1):892.
Li R, Rashotte AM, Singh NK, et al. Integrated signaling networks in plant responses to sedentary endoparasitic nematodes: a perspective [J]. Plant Cell Rep. 2015;34(1):5–22.
Wang JP, Xu YP, Munyampundu JP, et al. Calcium-dependent protein kinase (CDPK) and CDPK-related kinase (CRK) gene families in tomato: genome-wide identification and functional analyses in disease resistance [J]. Mol Gen Genomics. 2016;291(2):661–76.
Huang S, Gao Y, Liu J, et al. Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum [J]. Mol Gen Genomics. 2012;287(6):495–513.
Ning W, Zhai H, Yu J, et al. Overexpression of Glycine soja WRKY20 enhances drought tolerance and improves plant yields under drought stress in transgenic soybean [J]. Mol Breed. 2017;37(2):19.
Nguyễn PV, Bellafiore S, Petitot AS, et al. Meloidogyne incognita-rice (Oryza sativa) interaction: a new model system to study plant-root-knot nematode interactions in monocotyledons [J]. Rice. 2014;7(1):1–13.
Peng X, Hu Y, Tang X, et al. Constitutive expression of rice WRKY30 gene increases the endogenous jasmonic acid accumulation, PR gene expression and resistance to fungal pathogens in rice [J]. Planta. 2012;236(5):1485–98.
Cui X, Yan Q, Gan S, et al. GmWRKY40, a member of the WRKY transcription factor genes identified from Glycine max L., enhanced the resistance to Phytophthora sojae [J]. BMC Plant Biol. 2019;19(1):1–15.
Higashi K, Ishiga Y, Inagaki Y, et al. Modulation of defense signal transduction by flagellin-induced WRKY41 transcription factor in Arabidopsis thaliana [J]. Mol Gen Genomics. 2008;279(3):303–12.
Atamian HS, Eulgem T, Kaloshian I. SlWRKY70 is required for Mi-1-mediated resistance to aphids and nematodes in tomato [J]. Planta. 2012;235(2):299–309.
Rushton PJ, Somssich IE, Ringler P, et al. WRKY transcription factors [J]. Trends Plant Sci. 2010;15(5):247–58.
Tao Z, Liu H, Qiu D, et al. A pair of allelic WRKY genes play opposite roles in rice-bacteria interactions [J]. Plant Physiol. 2009;151(2):936–48.
Zhang Y, Li D, Zhang H, et al. Tomato histone H2B monoubiquitination enzymes SlHUB1 and SlHUB2 contribute to disease resistance against Botrytis cinerea through modulating the balance between SA-and JA/ET-mediated signaling pathways [J]. BMC Plant Biol. 2015;15(1):252.
Matthews BF, Beard H, Brewer E, et al. Arabidopsis genes, AtNPR1, AtTGA2 and AtPR-5, confer partial resistance to soybean cyst nematode (Heterodera glycines) when overexpressed in transgenic soybean roots [J]. BMC Plant Biol. 2014;14(1):96.
Ma KW, Ma W. Phytohormone pathways as targets of pathogens to facilitate infection [J]. Plant Mol Biol. 2016;91(6):713–25.
Cao FY, DeFalco TA, Moeder W, et al. Arabidopsis Ethylene Response Factor 8 (ERF8) has dual functions in ABA signaling and immunity [J]. BMC Plant Biol. 2018;18(1):211.
Kang YJ, Kim KH, Shim S, et al. Genome-wide mapping of NBS-LRR genes and their association with disease resistance in soybean [J]. BMC Plant Biol. 2012;12(1):139.
Geng X, Shen M, Kim JH, et al. The Pseudomonas syringae type III effectors AvrRpm1 and AvrRpt2 promote virulence dependent on the F-box protein COI1[J]. Plant Cell Rep. 2016;35(4):921–32.
Axtell MJ, Staskawicz BJ. Initiation of RPS2-specified disease resistance in Arabidopsis is coupled to the AvrRpt2-directed elimination of RIN4[J]. Cell. 2003;112(3):369–77.
Mackey D, Belkhadir Y, Alonso JM, et al. Arabidopsis RIN4 is a target of the type III virulence effector AvrRpt2 and modulates RPS2-mediated resistance [J]. Cell. 2003;112(3):379–89.
Takahashi A, Casais C, Ichimura K, et al. HSP90 interacts with RAR1 and SGT1 and is essential for RPS2-mediated disease resistance in Arabidopsis [J]. Proc Natl Acad Sci. 2003;100(20):11777–82.
Kandoth PK, Ithal N, Recknor J, et al. The soybean Rhg1 locus for resistance to the soybean cyst nematode Heterodera glycines regulates the expression of a large number of stress-and defense-related genes in degenerating feeding cells [J]. Plant Physiol. 2011;155(4):1960–75.
Van Ooijen G, Lukasik E, Van Den Burg HA, et al. The small heat shock protein 20 RSI2 interacts with and is required for stability and function of tomato resistance protein I-2[J]. Plant J. 2010;63(4):563–72.
Hubert DA, Tornero P, Belkhadir Y, et al. Cytosolic HSP90 associates with and modulates the Arabidopsis RPM1 disease resistance protein [J]. EMBO J. 2003;22(21):5679–89.
Barcala M, García A, Cubas P, et al. Distinct heat-shock element arrangements that mediate the heat shock, but not the late-embryogenesis induction of small heat-shock proteins, correlate with promoter activation in root-knot nematode feeding cells [J]. Plant Mol Biol. 2008;66(1–2):151–64.
Bagnaresi P, Sala T, Irdani T, et al. Solanum torvum responses to the root-knot nematode Meloidogyne incognita [J]. BMC Genomics. 2013;14(1):540.
Hussey RS. A comparison of methods of collecting inocula of Meloidogyne spp., including a new technique [J]. Plant Dis Rep. 1973;57:1025–8.
Kong LA, Wu DQ, Huang WK, et al. Large-scale identification of wheat genes resistant to cereal cyst nematode Heterodera avenae using comparative transcriptomic analysis [J]. BMC Genomics. 2015;16(1):801.
Taylor AL, Sasser JN. Biology, identification and control of root-knot nematodes (Meloidogyne spp.)[J]. Raleigh: North Carolina State University Graphics; 1978.
Schroeder A, Mueller O, Stocker S, et al. The RIN: an RNA integrity number for assigning integrity values to RNA measurements [J]. BMC Mol Biol. 2006;7(1):3.
Zhou R, Yu X, Zhao T, et al. Physiological analysis and transcriptome sequencing reveal the effects of combined cold and drought on tomato leaf [J]. BMC Plant Biol. 2019;19(1):377.
Zhao T, Liu W, Zhao Z, et al. Transcriptome profiling reveals the response process of tomato carrying Cf-19 and Cladosporium fulvum interaction [J]. BMC Plant Biol. 2019;19(1):1–12.
Ernst J, Bar-Joseph Z. STEM: a tool for the analysis of short time series gene expression data [J]. BMC Bioinformatics. 2006;7(1):191.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis [J]. BMC Bioinformatics. 2008;9(1):559.
Botía JA, Vandrovcova J, Forabosco P, et al. An additional k-means clustering step improves the biological features of WGCNA gene co-expression networks [J]. BMC Syst Biol. 2017;11(1):47.
Moschen S, Higgins J, Di Rienzo JA, et al. Network and biosignature analysis for the integration of transcriptomic and metabolomic data to characterize leaf senescence process in sunflower [J]. BMC Bioinformatics. 2016;17(5):174.
Xu P, Yang J, Liu J, et al. Identification of glioblastoma gene prognosis modules based on weighted gene co-expression network analysis [J]. BMC Med Genet. 2018;11(1):1–12.
Köppel R, Bucher TB. Duplex real-time PCR for the determination of wasabi (Eutrema wasabi) contents in horseradish (Armoracia rusticana) products applying the ΔΔct-method [J]. Eur Food Res Technol. 2016;242(7):1111–5.
We thank the Gene Denovo at Guangzhou for its assistance in related bioinformatics analysis.
This work was supported by the Major Research Plan (2016YFD0101703) and National Outstanding Youth Foundation of China (Grant No. 31501777).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Disease resistance statistics of LA3858 and Moneymaker infected with M. incognita after treatment with different soil temperatures.
List of hub genes filtered from darkorange and pink modules by WGCNA.
List of transcription factors significantly enriched in the 34 °C vs. 25 °C comparison at 0 (before inoculation), 3 and 6 dpi.
List of qRT-PCR primers used in this study.
Verification of the correlation results between the RNA-seq and qRT-PCR data
Level 2 GO terms identified for DEGs in 34 °C vs. 25 °C comparison at 0 (before inoculation), 3 and 6 dpi.
Top 20 significantly enriched pathways in the darkorange and pink modules according to WGCNA.
Description of all the analyses performed in this work.
About this article
Cite this article
Du, C., Jiang, J., Zhang, H. et al. Transcriptomic profiling of Solanum peruvianum LA3858 revealed a Mi-3-mediated hypersensitive response to Meloidogyne incognita. BMC Genomics 21, 250 (2020). https://doi.org/10.1186/s12864-020-6654-5
- RNA sequencing
- Root-knot nematode
- Soil temperature