ORIGINAL RESEARCH article
Transcriptomic Analysis of Resistant and Susceptible Responses in a New Model Root-Knot Nematode Infection System Using Solanum torvum and Meloidogyne arenaria
- 1RIKEN Center for Sustainable Resource Science, Yokohama, Japan
- 2Central Region Agricultural Research Center, National Agriculture and Food Research Organization, Tsukuba, Japan
- 3INRES – Molecular Phytomedicine, Rheinische Friedrich-Wilhelms-University of Bonn, Bonn, Germany
- 4School of Life Science and Technology, Tokyo Institute of Technology, Yokohama, Japan
- 5NIBB Core Research Facilities, National Institute for Basic Biology, Okazaki, Japan
- 6RIKEN BioResource Research Center, Tsukuba, Japan
- 7Institute of Cellular and Molecular Botany, Rheinische Friedrich-Wilhelms-University of Bonn, Bonn, Germany
- 8Department of Entomology and Nematology, University of California, Davis, Davis, CA, United States
- 9Department of Biological Chemistry, College of Bioscience and Biotechnology, Chubu University, Kasugai, Japan
- 10Graduate School of Science, The University of Tokyo, Bunkyo, Japan
Root-knot nematodes (RKNs) are among the most devastating pests in agriculture. Solanum torvum Sw. (Turkey berry) has been used as a rootstock for eggplant (aubergine) cultivation because of its resistance to RKNs, including Meloidogyne incognita and M. arenaria. We previously found that a pathotype of M. arenaria, A2-J, is able to infect and propagate in S. torvum. In vitro infection assays showed that S. torvum induced the accumulation of brown pigments during avirulent pathotype A2-O infection, but not during virulent A2-J infection. This experimental system is advantageous because resistant and susceptible responses can be distinguished within a few days, and because a single plant genome can yield information about both resistant and susceptible responses. Comparative RNA-sequencing analysis of S. torvum inoculated with A2-J and A2-O at early stages of infection was used to parse the specific resistance and susceptible responses. Infection with A2-J did not induce statistically significant changes in gene expression within one day post-inoculation (DPI), but afterward, A2-J specifically induced the expression of chalcone synthase, spermidine synthase, and genes related to cell wall modification and transmembrane transport. Infection with A2-O rapidly induced the expression of genes encoding class III peroxidases, sesquiterpene synthases, and fatty acid desaturases at 1 DPI, followed by genes involved in defense, hormone signaling, and the biosynthesis of lignin at 3 DPI. Both isolates induced the expression of suberin biosynthetic genes, which may be triggered by wounding during nematode infection. Histochemical analysis revealed that A2-O, but not A2-J, induced lignin accumulation at the root tip, suggesting that physical reinforcement of cell walls with lignin is an important defense response against nematodes. The S. torvum-RKN system can provide a molecular basis for understanding plant-nematode interactions.
Plant-parasitic nematodes (PPNs) infect a broad range of commercially important crop families such as the Solanaceae (tomato, potato, pepper), Fabaceae (soybean, lucerne, lentils), Malvaceae (cotton), Amaranthaceae (sugar beet), and Poaceae (syn. Gramineae; rice, wheat, maize), causing an estimated annual loss of $80 billion USD (Jones et al., 2013; Sato et al., 2019). The most economically important group of PPNs are sedentary endoparasites, including root-knot nematodes (RKNs) and cyst nematodes (CNs) (Palomares-Rius et al., 2017). Sedentary endoparasites induce the formation of permanent feeding cells that provide specialized nutrient sources for nematodes (Bartlem et al., 2014; Siddique and Grundler, 2018; Smant et al., 2018). Infective second-stage juveniles (J2s) of RKNs (Meloidogyne spp.) predominantly invade near the root tip and then migrate intercellularly toward the apical meristematic region without crossing the endodermis, making a U-turn to enter the vascular cylinder where they induce several giant cells as a feeding site by stimulating the redifferentiation of root cells into multinucleate giant cells by repeated nuclear divisions without cytoplasmic division. After maturation, adult RKN females lay eggs in a gelatinous egg mass on or below the surface of the root (Sijmons et al., 1991; Abad et al., 2009; Escobar et al., 2015). In contrast, CNs move destructively through cells into the vascular cylinder, select a single cell, and form a syncytium as a feeding site by local dissolution of cell walls and protoplast fusion of neighboring cells. A CN female produces hundreds of eggs and its body forms a cyst that can protect the eggs for many years in the soil (Wyss and Zunke, 1986; Bohlmann and Sobczak, 2014; Bohlmann, 2015).
Nematicides have been commonly used to control PPNs in agriculture, but some nematicides such as methyl bromide and aldicarb are currently banned from use in many countries due to their negative effects on the environment and human health (Zasada et al., 2010; Brennan et al., 2020; Oka, 2020). It has therefore become important to understand the molecular mechanisms of plant immunity against PPNs to provide a foundation for the development of new environmentally friendly and effective control methods.
In general, the plant immune system is represented by two inter-related tiers (Jones and Dangl, 2006; Dodds and Rathjen, 2010). The first is governed by cell surface-localized pattern recognition receptors (PRRs) that perceive pathogen-associated molecular patterns (PAMPs), leading to pattern-triggered immunity (PTI) (Boutrot and Zipfel, 2017). Successful pathogens secrete effector molecules into the apoplast or directly into plant cells, which interfere with PTI, resulting in successful infection. Resistant plants recognize cell-invading effectors through recognition by intracellular nucleotide-binding domain leucine-rich repeat (NLR)-type immune receptors, which are encoded by resistance (R) genes. Similar mechanisms are also conserved in plant-PPN interactions. For example, the well-conserved nematode pheromone ascaroside has been identified as a PAMP (Manosalva et al., 2015), but the corresponding PRR has not yet been found. PPN genome sequence analyses identified a number of candidate virulence effectors (summarized in Mejias et al., 2019), and a handful of NLR protein-encoding R genes involved in PPN recognition have been well-studied and characterized, including tomato Mi-1.2, Mi-9, and Hero-A; potato Gpa2 and Gro1-4; pepper CaMi; and prune Ma (Milligan et al., 1998; van der Vossen et al., 2000; Ernst et al., 2002; Paal et al., 2004; Chen et al., 2007; Jablonska et al., 2007; Claverie et al., 2011; Kaloshian et al., 2011). Mi-1.2, Mi-9, CaMi, and Ma confer resistance against RKNs, whereas Hero-A, Gpa2, and Gro1-4 provide resistance against CNs.
Although the PPN perception mechanism is somewhat clearer at the molecular level, it is still largely unknown what kind of downstream responses are induced after the recognition of avirulent PPNs. It is also unclear what kind of host responses are induced after infection with virulent PPNs, leading to susceptibility and infestation. There are several difficulties in working on plant responses against PPNs. First and foremost, most model plants, such as Arabidopsis, are susceptible to PPNs and therefore cannot be used to study the cascade of responses leading to resistance. Second, PPNs migrate long-distances inside roots, inducing complicated responses as they go, triggered by mechanical stress and wounding, among others, making it difficult to pinpoint the key genes involved in resistance or susceptibility by transcriptome analyses. Some studies have used comparative transcriptomics using susceptible and resistant plants infected with a single genotype of nematode (Postnikova et al., 2015; Xing et al., 2017; Ye et al., 2017; Zhang H. et al., 2017). However, it is difficult to rule out the possibility that differences in gene expression were due to resistance or susceptibility rather than to differences in the genetic backgrounds of host plants. Lastly, susceptible responses such as the formation of feeding sites are induced in specific cells targeted by PPNs, and defense responses are likely to be induced in the cells directly impacted by PPN activity. Thus, cells responding to PPNs are rather limited, making the analysis technically challenging.
Here we have introduced Solanum torvum Sw “Torvum Vigor” to overcome these problems. S. torvum has been widely used as a rootstock of eggplant (aubergine) to prevent disease caused by PPNs, as well as the soil-borne pathogens Ralstonia solanacearum, Verticillium dahliae, and Fusarium oxysporum f. melongenae n. f. (Gousset et al., 2005; Yamaguchi et al., 2010; Bagnaresi et al., 2013; Yang et al., 2014; Uehara et al., 2017; García-Mendívil and Sorribas, 2019; García-Mendívil et al., 2019; Murata and Uesugi, 2021). S. torvum Sw “Torvum Vigor” is resistant to Meloidogyne arenaria pathotype A2-O (A2-O), but susceptible to M. arenaria pathotype A2-J (A2-J) (Uehara et al., 2017). By using S. torvum and avirulent or virulent isolates, we established an in vitro infection system and performed comparative transcriptome analyses to identify genes whose expressions were associated with either resistance or susceptibility by carefully collecting only root tips attacked by RKNs, which allowed us to detect gene expression only in cells directly affected by nematodes. In addition, observation of infected root tip morphology suggests that the success or failure of the immune system against PPNs is determined within a few days of invasion. Thus, we decided to focus on the transcriptional changes that occurred in the very early stages of infection, which has not been studied in previous transcriptomic analyses (Bagnaresi et al., 2013; Postnikova et al., 2015; Xing et al., 2017).
Comparative clustering analyses of gene expression identified a large number of novel genes, especially those involved in susceptibility through cell wall modification and transmembrane transport; resistance through lignin and isoprenoid biosynthesis and fatty acid metabolism; and suberin biosynthesis in mechanical wounding. Consistent with the transcriptional up-regulation of lignin biosynthetic genes from A2-O invasion, lignin is accumulated at the root tip of S. torvum infected with avirulent A2-O but not with virulent A2-J, suggesting that S. torvum reinforces the cell wall as a defense response against the avirulent RKN.
Materials and Methods
Plant Materials and Growth Conditions
Seeds of S. torvum cultivar “Torvum Vigor” were sown on half-strength Murashige-Skoog (MS) medium containing 1% sucrose and 0.5% Gelrite (Wako, Japan). Plants were grown in a controlled growth chamber under long-day photoperiods (16 hours light/8 hours dark) at 25°C.
Nematode Infection Assay
M. arenaria pathotypes A2-J and A2-O were propagated on Solanum lycopersicum cultivar “Micro-Tom” in a greenhouse. Nematode eggs were isolated from infected roots and then hatched at 25 °C. Freshly hatched J2s were collected and transferred to a Kimwipe filter (a folded Kimwipe tissue) placed on the top of a glass beaker filled with sterilized distilled water (SDW) containing 100 μg/ml streptomycin and 10 μg/ml nystatin. Only active J2s pass through the filter. Filtered J2s were surface sterilized with 0.002 % mercuric chloride, 0.002 % sodium azide, and 0.001 % Triton X-100 for 10 min, and then rinsed three times with SDW (Mitchum et al., 2004). Eleven-day-old S. torvum seedlings grown on the MS-Gelrite in 6-well plates were inoculated with 200–300 J2s resuspended in SDW. The plates were wrapped in aluminum foil for 2–3 days after inoculation to promote nematode infection. When mature giant cells were observed 18 days post-inoculation (DPI), we used the MS-Gelrite media without sucrose to prevent the formation of callus-like structures. The difference in the number of normal galls formed by A2-J or A2-O at 4 DPI was statistically tested using the Mann-Whitney U test with R software (v3.6.3). Nematodes resident in root tissues were stained with acid fuchsin 2–4 DPI (Bybd et al., 1983), photographed by light microscopy (Olympus BX51, Olympus, Japan), and the photomicrographs were processed using cellSens (Olympus, Japan). For the observation of giant cells and developing nematodes at 18 DPI, infection sites were fixed with glutaraldehyde and cleaned with benzyl-alcohol/benzyl-benzoate (BABB) (Cabrera et al., 2018). We observed BABB-cleaned samples by confocal laser scanning microscopy (TCS SP5, Leica Microsystems GmbH, Germany). Photomicrographs were processed using LAS X software (version 220.127.116.1179, Leica Microsystems GmbH, Germany).
RNA-Sequencing and de novo Assembly
S. torvum seedlings were grown on half-strength MS-Gelrite medium containing 1% sucrose. Eleven-day-old seedlings were treated with SDW as a mock infection or with 200–300 J2s of M. arenaria A2-J for susceptible infection or A2-O for resistant infection. Root tips attacked by the nematodes were checked under microscopy, and more than 50 root tips were cut (approximately 3–5 mm from the tip) and collected for each treatment (Figure 2A). Root tip samples were collected at 1, 2, and 3 DPI with four biological replicates. Whole shoot and root samples were collected at 1, 3, 6, and 9 DPI with four biological replicates. Root tip samples were used for de novo assembly and differential gene expression analyses, and whole shoot and root samples were used only for de novo assembly.
RNA-seq libraries were prepared from the collected samples using a high-throughput RNA-seq method (Kumar et al., 2012). The 85-bp paired-end reads for the root tip samples, and the 85-bp single-end reads for the whole shoot and root samples were sequenced on an Illumina NextSeq 500 platform (Illumina, CA, United States). The FASTX toolkit 0.0.13.2 (Hannonlab) was used for quality filtering. Low-quality nucleotides (Phred scores of <30) were removed from the 3′ ends, and short reads (<76bp) were excluded. Reads with at least 95% of nucleotides with Phred scores > 20 were kept and used for the downstream analyses (Supplementary Table 1A). Adaptor sequences were removed using custom scripts written in Perl (Kumar et al., 2012). Filtered reads were mapped to the genome assembly of M. arenaria A2-J (GenBank accession number JAEEAS010000000) or A2-O (GenBank accession number QEUI01000000) (Sato et al., 2018) using HISAT2 (version 2.1.0) (Kim et al., 2019) to exclude reads of nematode origin. Unmapped reads were used for de novo transcriptome assembly (Supplementary Figure 1A).
Three different transcriptome assemblers were used for de novo assembly: SOAPdenovo-Trans v1.03 (Xie et al., 2014), Velvet v1.2.10 (Zerbino and Birney, 2008)/Oases v0.2.09 (Schulz et al., 2012) and Trinity package v2.4.0 (Grabherr et al., 2011; Haas et al., 2013). Unmapped paired-end and single-end reads were normalized using Trinity and assembled independently (Mamrot et al., 2017). For assemblies of each dataset, SOAPdenovo-Trans was set at k-mer sizes: 21, 23, 25, 27, 29, 31, 35, 41, 51, 61, Velvet/Oases was set at k-mer sizes: 21, 23, 25, 27, 29, 31, 35, 41, 51, 61, and Trinity was set at k-mer size 25 (Supplementary Table 1B).
Assemblies were merged into a non-redundant dataset using the EvidentialGene pipeline1 as previously described (Nakasugi et al., 2014). Oases assembled scaffolds were split at gaps into contigs before merging with contigs from the other assemblies with the EvidentialGene tr2aacds pipeline. The tr2aacds pipeline produces ‘primary’ and ‘alternate’ sequences of non-redundant transcripts with ‘primary’ transcripts being the longest coding sequence for a predicted locus. Next, we used the evgmrna2tsa program from EvidentialGene to generate mRNA, coding, and protein sequences. BUSCO v3.0.2 (Simão et al., 2015) was applied for quantitative assessment of assembly completeness. This assembly and one previously reported for S. torvum by Yang et al. (2014) (GenBank accession number GBEG01000000) were compared to the Embryophyta odb9 dataset, which contains 1,440 BUSCO groups. The homology of the contigs from the final assembly was searched against the NCBI non-redundant database using BLASTX (BLAST+ v2.7.1) with an e-value threshold of 1E-05. We also compared the contigs with Arabidopsis genome annotation (TAIR10) (Berardini et al., 2015) using BLASTX at the e-value cutoff of 10. Results of the annotation are summarized in Supplementary Table 2.
Differential Expression Analysis
Quantification of S. torvum transcripts was performed with the mapping-based mode of Salmon v0.10.2 (Patro et al., 2017) by using reads that did not map to nematode genome assemblies. Quantified transcript-level abundance data were imported to R (v3.6.3) using tximport v.1.14.2 (Soneson et al., 2015) package, and differential gene expression analysis was carried out with the edgeR package (v3.28.1) (Robinson et al., 2010). Transcripts with very low expression values were filtered out with the “filterByExpr” function. Counts per million are shown in Supplementary Table 3. Differentially expressed genes (DEGs) (false discovery rate (FDR) ≤ 0.01, log-transformed fold change (logFC) ≥ 1 or logFC ≤ −1) were identified using the quasi-likelihood F-test by comparing expression levels during infection with A2-J or A2-O to SDW (mock) treatment at the same time point (Supplementary Table 4). We also compared gene expression levels between S. torvum infected with A2-J and A2-O, and all analyzed genes were listed in Supplementary Table 5.
Principal Component Analysis With Self-Organizing Map Clustering
To group genes by expression pattern, we applied the self-organizing map (SOM) clustering method on genes within the top 25 % of the coefficient of variation for expression across samples as previously described (Chitwood et al., 2013; Ranjan et al., 2014; Goto et al., 2020). Scaled expression values, representing the average principal component (PC) values among each gene in a cluster were used for multilevel three-by-three hexagonal SOM (Wehrens and Buydens, 2007). One hundred training iterations were used during the clustering, over which the α learning rate decreased from 0.0071 to 0.0061. The final assignment of genes to winning units formed the basis of the gene clusters. The results of SOM clustering were visualized in a principal component analysis (PCA) space where PC values were calculated based on gene expression across samples (Figure 4A).
Functional Annotation of Transcriptome and Gene Ontology Enrichment Analysis
We compared the contigs of our assembly with the NCBI non-redundant database using BLASTX (BLAST+ v2.7.1) with an e-value threshold of 1E-05. In addition, predicted amino acid sequences that begin with methionine were also annotated using InterProScan (v5.32-71.0) (Jones et al., 2014). BLASTX and InterProScan outputs were used for Blast2GO (v5.2.5) analysis to annotate the contigs with Gene Ontology (GO) terms (Gotz et al., 2008). GO enrichment analyses of the sets of genes induced by A2-O infection at 1 DPI or that were assigned to each cluster generated by SOM was performed by comparison with all genes using GO terms generated by Blast2GO at the FDR cutoff of 1E-04 (Gotz et al., 2008). We further used the “Reduce to most specific terms” option in Blast2GO to remove general GO terms and obtain only the most specific GO terms.
Histochemistry of Lignin Deposition
Lignin deposition during infection with A2-J or A2-O was visualized by phloroglucinol-HCl staining as previously described (Jensen, 1962). Eleven-day-old S. torvum seedlings were inoculated with nematodes or treated with SDW (mock) as described above. We collected root tips 3 DPI for lignin staining with phloroglucinol-HCl. Microphotographs were taken and processed as described above and combined manually.
Aliphatic Suberin Monomer Analysis
Quantification of aliphatic suberin was performed as described previously (Holbein et al., 2019). Eleven-day-old plants were treated with SDW (mock) or infected with A2-J or A2-O. At 4 DPI, root tips inoculated with nematodes were microscopically checked for infection, and more than 50 infected root tips were cut (approximately 3–5 mm from the tip) and collected for each treatment. To remove unbound lipids, samples were extracted in methanol for 24 h then in chloroform for 24 h, dried, and weighed. Samples were depolymerized and analyzed by gas chromatography-mass spectrometry (GC-MS) (Agilent 6890N-Agilent 5973N quadrupole mass-selective detector, Agilent Technologies, Germany) for monomer identification and for quantitative analysis based on an internal standard using an identical gas chromatography system coupled with a flame ionization detector as described previously (Franke et al., 2005).
A2-O Induces a Browning Response, and A2-J Induces Gall Formation in S. torvum
To understand the differential responses of S. torvum to M. arenaria A2-J and A2-O, we first established an in vitro infection system. Seedlings of S. torvum were grown in MS-Gelrite plates for 11 days and then inoculated with 200–300 J2s of A2-J or A2-O. At 4 DPI, more than 90 % of root tips infected with A2-J induced the formation of gall-like structures ranging in size. These galls are classified here as “normal” galls, while the rest produced brown pigments. Normal galls lacked obvious brown pigment accumulation and were further classified based on the width of the gall into small (shorter than 0.5 mm), medium (0.5–0.8 mm), and large (wider than 0.8 mm) (Figure 1A). In contrast, about 60 % of A2-O-infected root tips accumulated at least some brown pigment. Some of these brownish root tips also had an abnormal appearance due to the formation of balloon-like structures, and others had many localized and highly pigmented spots. There were a very few small gall-like structures formed after infection with A2-O, but far fewer and smaller than in root tips infected with A2-J (Figure 1B). RKN staining by acid fuchsin revealed that both A2-J and A2-O successfully invaded the roots (Figure 1C). Interestingly, host cells invaded by A2-O uniformly accumulated brownish pigments, suggesting that the surrounding tissue is strongly responding to, and highly correlated with A2-O infection, a response that was absent from A2-J infected roots. It is generally known that browning of plant tissue is related to enzymatic or non-enzymatic oxidation of phenolic substances (Mesquita and Queiroz, 2013), but the identity of the brown pigments synthesized upon infection with A2-O is unknown. By 18 DPI, A2-J had induced the formation of mature multinucleate giant cells and developed into fourth stage juveniles (Figure 1D). In contrast, A2-O did not induce the formation of giant cells nor develop past second stage juveniles. These results suggest that S. torvum rapidly induces defense responses against A2-O, which inhibits the maturation of A2-O and gall formation. In contrast, A2-J inhibits or evades the induction of defense responses, continues development, and induces gall formation.
Figure 1. M. arenaria A2-J induced gall formation in S. torvum upon infection, while M. arenaria A2-O infection resulted in the accumulation of brown pigment at root tips. (A) Microscopic examination of S. torvum root tips infected with A2-J or A2-O. Seedlings of S. torvum were infected with 200–300 M. arenaria A2-J or A2-O juveniles. The number of galls and brown root tips per plant were counted at 4 DPI. Normal galls without accumulation of brown pigment were binned by gall length into small (shorter than 0.5 mm), medium (0.5–0.8 mm), or large (wider than 0.8 mm). Scale bars indicate 1 mm. The average number of galls per plant (n = 25) is shown as a bar chart. Experiments were performed four times with similar results. (B) A2-J induced the formation of more galls than A2-O. The violin plot indicates variation in the number of galls per plant at 4 DPI. A p-value was calculated by the Mann-Whitney U test from the mean number of galls per plant. Similar results were obtained from four independent experiments. (C) Both A2-J and A2-O entered roots. Nematodes (N) in the roots were stained with acid fuchsin 2–4 DPI. The arrow indicates the accumulation of brown pigment. Scale bars indicate 200 μm. (D) S. torvum inhibited the growth of A2-O, but not A2-J. Nematodes and giant cells in the S. torvum roots 18 DPI were visualized by a method described previously (Cabrera et al., 2018). Nematodes and giant cells are indicated by N and asterisks, respectively. Scale bars indicate 200 μm.
RNA-Seq Analysis of S. torvum Root Tips Infected With A2-O or A2-J
RNA-seq analysis was performed to understand the differences in transcriptional regulation of the S. torvum response to infection by nematodes that induce an immune response or that are successful parasites. Eleven-day-old S. torvum seedlings grown on MS-Gelrite plates were inoculated with 200–300 surface-sterilized J2s of A2-J or A2-O, or treated with SDW (mock treatment) in vitro. Since there were clear morphological differences between the root tips infected with A2-J and A2-O after four days (Figure 1A), it should follow that the success or failure of infection is determined within a few days post-inoculation. We therefore decided to analyze the transcriptome at 1–3 DPI, corresponding to the early stages of infection. In addition, to detect gene expression in cells directly affected by the nematodes, we carefully collected infected root tips under a stereomicroscope (Figure 2A). Root tips were cut with precision forceps and flash-frozen with liquid nitrogen to preclude the induction of wound responses. More than 50 root tips were pooled for each treatment, and four biological replicates were used for the RNA-seq based transcriptome analyses. We also carried out RNA-seq of whole roots and shoots of S. torvum infected with A2-J or A2-O, or mock treatment (1, 3, 6, and 9 DPI) to improve the completeness of de novo transcriptome assembly. As a result, we obtained 218,024,788 paired-end reads from root tips and 341,297,551 single-end reads from whole shoots and roots after quality filtering (Supplementary Table 1A). After removing the reads derived from nematodes, we performed de novo assembly using multiple assemblers with a variety of k-mer sizes (Supplementary Figure 1A and Supplementary Table 1B). A set of non-redundant transcripts was generated by merging these multiple assemblies. The final assembly had 88,596 contigs with an N50 of 1,298 bp, an average size of 800.62 bp, and a total length of 70,931,593 bp (Table 1). We assessed the accuracy and completeness of the final assembly using BUSCO. The assembly included an estimated > 95 % of the assessed dataset, improving the current status of the transcriptome assembly of S. torvum (Yang et al., 2014; Supplementary Figure 1B) and provided a high-quality transcriptome assembly of S. torvum for further analyses.
Figure 2. Transcriptome analysis of DEGs in S. torvum infected with M. arenaria A2-J or A2-O. (A) Sampling of S. torvum root tips infected by A2-J or A2-O for RNA-seq analysis. Root tips (approximately 3–5 mm from the tip) were cut (scissors not to scale) and collected at 1, 2, and 3 DPI for RNA-seq. (B) Venn diagrams show the overlap of up-regulated genes (logFC ≥ 1, FDR ≤ 0.01) and down-regulated genes (logFC ≤ −1, FDR ≤ 0.01) for at least one time point after infection with A2-J or A2-O. (C,D) Infection with A2-O induced the expression of phytohormone marker genes and WRKY transcription factors. PR1 and PR2 are marker genes of salicylic acid (SA), LoxD and PI-II are marker genes for jasmonic acid (JA), and ACS and ACO are marker genes for ethylene (ET). The homologs of tomato WRKY genes WRKY31, WRKY33, WRKY39, and WRKY80 positively regulate immunity. The logFC values of the genes (compared to mock treatment) at 1, 2, and 3 DPI are shown. Heatmaps indicate relative gene expression levels. Gray color indicates values with no statistically significant differences at | logFC| ≥ 1 and FDR ≤ 0.01.
Infection With A2-O, but Not A2-J, Induces the Expression of Defense-Related Genes
Differential expression analysis showed that 1,220 genes were significantly up-regulated and 261 genes were down-regulated upon infection with A2-J, while 2,535 genes were up-regulated and 802 genes were down-regulated by infection with A2-O at at least one-time point during root tip infection, compared to the mock treatment (logFC ≥ 1 for up-regulation and logFC ≤ −1 for down-regulation, FDR ≤ 0.01, Figure 2B; Supplementary Figure 2; Supplementary Table 4). 1,029 genes were up-regulated, and 180 genes were down-regulated at at least one-time point in both A2-J and A2-O infected plants (Figure 2B). Previous studies showed that the expression of genes associated with the salicylic acid (SA), jasmonic acid (JA), and ethylene (ET) signaling pathways are induced in resistant plants infected with PPNs (Ye et al., 2017; Zhang H. et al., 2017; Bali et al., 2019; Du et al., 2020; Ghaemi et al., 2020), so we investigated the expression of marker genes for hormone biosynthesis, hormone signaling (Figure 2C), and defense responses (Figure 2D). Infection with A2-O significantly induced the expression of SA-regulated genes that encode basic pathogenesis-related protein 1 (PR1) and β-1,3-glucanase (PR2), as well as a JA biosynthesis gene (lipoxygenase D (LoxD)), the JA-regulated gene protease inhibitor II (PI-II), and two ET biosynthesis genes (1-aminocyclopropane-1-carboxylate synthase (ACS) and 1-aminocyclopropane-1-carboxylate oxidase (ACO)) (Farmer et al., 1992; Yan et al., 2013; Booker and DeLong, 2015; Zhang et al., 2018; Marhavý et al., 2019; Figure 2C). Infection with A2-J induced the expression of PR1 and PR2 to a lesser extent, but did not induce the expression of LoxD, PI-II, ACS, or ACO.
Tomato (S. lycopersicum) WRKY transcription factors are induced by pathogen infection and positively regulate defense responses against pathogens (Liu et al., 2014; Li et al., 2015a, b; Zhou et al., 2015). For example, SlWRKY31 and SlWRKY33 play important roles in resistance against fungal and oomycete pathogens, while SlWRKY39 is involved in resistance against bacterial pathogens (Huang et al., 2012; Sun et al., 2015). SlWRKY80 is required for Mi-1.2-mediated resistance against RKN Meloidogyne javanica (Atamian et al., 2012). We found that the S. torvum homologs of tomato WRKY31, 33, and 39 were specifically up-regulated upon infection with A2-O, but not with A2-J (Figure 2D). In the case of WRKY80, A2-J induced its expression at a later time point (3 DPI), while A2-O induced the WRKY80 homolog to much higher levels. These results suggest that A2-O, but not A2-J, strongly and specifically induces the expression of defense-related genes in S. torvum.
Infection With A2-O Rapidly Induces the Expression of Sesquiterpene Synthases and Class III Peroxidases
Infection with A2-O caused an increase in expression of more genes and to a greater extent than A2-J (Figures 2B, 3A and Supplementary Figure 2). Importantly, at 1 DPI, A2-J did not induce any statistically significant changes in the expression of any genes, whereas A2-O induced 204 genes, suggesting that infection with A2-O rapidly induces the expression of early responsive genes, which is prevented or avoided in A2-J infection. Since the speed of a defense response is one of the most important factors for successful immunity against pathogens, we hypothesized there must be important defense components among the 204 up-regulated genes. We therefore performed a GO enrichment analysis to identify significantly represented GO terms amongst the 204 up-regulated genes (FDR ≤ 1E-04, Supplementary Table 6). The list of enriched GO terms was further reduced using “Reduce to most specific terms” option in Blast2GO to remove general GO terms and obtain only the most specific terms (Table 2). Some GO terms that were significantly enriched among the 204 genes were related to the biosynthesis of isoprenoids (“farnesyl diphosphate catabolic process (GO:0045339)”, “sesquiterpene biosynthetic process (GO:0051762)”, and “terpenoid biosynthetic process (GO:0016114)”). To follow up on this result, we checked the expression of all the genes up-regulated by A2-O that are related to isoprenoid biosynthesis and found that A2-O infection induced the expression of genes encoding sesquiterpene synthases, such as viridiflorene synthase, vetispiradiene synthase, germacrene C synthase-like protein, and 5-epiaristolochene synthase. Several other enzymes involved in isoprenoid biosynthesis, such as xanthoxin dehydrogenase-like protein and UDP-glycosyltransferase 91C1 were also up-regulated (Figure 3B). Sesquiterpene synthases convert farnesyl diphosphate to sesquiterpenes such as germacrene C, 5-epiaristolochene, viridiflorene, and vetispiradiene. Because some isoprenoids have nematicidal activity (Ohri and Pannu, 2009), it is possible that the sesquiterpenes produced by S. torvum in response to infection with A2-O are nematicidal and contribute to suppressing A2-O infection. Other GO terms significantly enriched among the 204 up-regulated genes were related to oxidative stress (“hydrogen peroxide catabolic process (GO:0042744)” and “response to oxidative stress (GO:0006979)”). The most up-regulated genes by A2-O in the GO term group were class III peroxidases, which are involved in lignification, cell elongation, seed germination, and response to abiotic and biotic stresses (Cosio and Dunand, 2009; Shigeto and Tsutsumi, 2016; Supplementary Figure 3). The transcriptional up-regulation of the class III peroxidases is consistent with the fact that resistant tomato lines more strongly elevate peroxidase activity during RKN infection than susceptible lines (Zacheo et al., 1993).
Figure 3. Infection with A2-O, but not with A2-J induced rapid expression of genes involved in S. torvum isoprenoid biosynthesis. (A) Volcano plots showing the distribution of DEGs in S. torvum infected with M. arenaria A2-J or A2-O compared to mock treatment. The logarithms of the fold change of individual genes are plotted against the negative logarithms of their FDR. Red dots represent genes up-regulated (logFC ≥ 1) or down-regulated (logFC ≤ -1) with an FDR ≤ 0.01. (B) Infection with A2-O rapidly induced the expression of the genes involved in isoprenoid biosynthesis listed as “farnesyl diphosphate catabolic process (GO:0045339)”, “sesquiterpene biosynthetic process (GO:0051762)”, and “terpenoid biosynthetic process (GO:0016114)”. LogFC values of the genes (compared to mock treatment) at 1, 2, and 3 DPI are shown. In the heatmap, gray boxes indicate no statistically significant difference at | logFC| ≥ 1 and FDR ≤ 0.01.
Infection With A2-J Induces Genes Related to Cell Wall Modification and Transmembrane Transport
To identify the expression pattern of genes that are specific to A2-J or A2-O infection and common in both pathotypes, we clustered genes according to their transcript profiles by PCA with SOM clustering. SOM clustering grouped 6,502 genes into nine clusters based on their differential gene expression profiles after mock treatment or infection with either A2-J or A2-O (Figure 4 and Supplementary Table 7). 429 genes in Cluster 2 and 554 genes in Cluster 4 were specifically up-regulated after infection with A2-J. In contrast, 1,769 and 600 genes in Cluster 8 and 9, respectively, were specifically up-regulated after infection with A2-O. 1,000 genes in Cluster 7 were up-regulated after infection with either A2-J or A2-O (Figure 4B).
Figure 4. PCA and SOM clustering of gene expression in S. torvum. (A) Clustering of gene expression patterns in S. torvum treated with SDW (mock) or infected with A2-J or A2-O at 1, 2, and 3 DPI. The expression profiles of each gene were represented in a PCA space separated by PC1 and PC2, with SOM node memberships (cluster) represented by different colors and numbers. The percentage shown along the x or y-axis represents the percentage of variance explained by each component. A total of nine clusters with different expression patterns were defined. (B) Scaled expression profiles of genes by SOM clustering. The number of genes assigned to each SOM cluster is shown in parentheses. Each dot represents the value of the scaled expression of the gene. The upper and lower horizontal lines of the box correspond to the 25th and 75th percentiles. The horizontal bar in the box indicates the median. The upper and lower end of horizontal bars outside the box indicate 90th and 10th percentiles.
We once again used GO enrichment to identify functional terms enriched in the genes in each cluster (FDR ≤ 1E-04, Supplementary Table 8) by further filtering enriched GO terms using the “Reduce to most specific terms” option in Blast2GO. In Cluster 4 (specifically A2-J up-regulated genes), significantly enriched GO terms were related to cell wall remodeling, including “cell wall modification (GO:0042545)”, “cell wall organization or biogenesis (GO:0071554)”, and “pectin catabolic process (GO:0045490)” (Table 3 and Supplementary Table 8). This is consistent with the observation that the expansion of giant cells is associated with an increase in cell wall thickness (Bartlem et al., 2014; Bozbuga et al., 2018). GO terms associated with the significantly A2-J up-regulated genes include enzymes such as cellulose synthase-like protein, xyloglucan endotransglucosylase/hydrolase protein, and a non-catalytic subunit of a polygalacturonase isozyme (Supplementary Figure 4A). The transcriptional up-regulation of these enzymes is consistent with the presence of the common polysaccharides pectic homogalacturonan, xyloglucan, and pectic arabinan in the cell walls of giant cells (Bozbuga et al., 2018). Similarly, A2-J infection also activated the expression of COBRA-like protein, expansin, and LRR-RLK PXC1, which play important roles in cellulose deposition, loosening of cell walls, and secondary wall formation (Cosgrove, 2000; Brown et al., 2005; Wang et al., 2013; Kumar et al., 2016; Supplementary Figure 4A). Another GO term significantly enriched in Cluster 4 was “transmembrane transport (GO:0055085)” (Supplementary Table 8). Other significantly up-regulated genes encode sugar transporter ERD6-like protein (Hammes et al., 2005) and amino acid transporter family protein (Elashry et al., 2013). These transporters may promote the uptake of nutrients into giant cells or alter transportation through cells surrounding giant cells (Supplementary Figure 4B). We also performed GO enrichment analyses for Cluster 2, but no GO terms were enriched.
Table 3. GO enrichment analysis of Cluster 4 genes up-regulated specifically by infection with M. arenaria A2-Ja.
In addition to GO enrichment analyses, we looked for interesting genes whose expression were dramatically up-regulated in Clusters 2 and 4. These clusters included genes encoding chalcone synthase and a spermidine synthase that were specifically and highly expressed after infection with A2-J (Figure 5A). Chalcone synthase is the first enzyme of the flavonoid biosynthetic pathway (Dao et al., 2011), and spermidine synthase is a key enzyme involved in polyamine biosynthesis (Liu et al., 2007). In summary, A2-J infection significantly and specifically up-regulates genes related to cell wall modification and membrane transport, chalcone synthase, and spermidine synthase at an early phase of gall formation.
Figure 5. Up-regulated genes in Clusters 4, 7, 8, and 9. Highly up-regulated genes (logFC ≥ 4, FDR ≤ 0.01) in each cluster are listed. (A) Infection with A2-J specifically up-regulated the expression of chalcone synthase and spermidine synthase (Cluster 4). (B) Infection with A2-O dramatically up-regulated defense-related proteins, sesquiterpene synthase, fatty acid desaturase 2 (FAD2), ferulic acid 5-hydroxylase, and berberine bridge enzyme-like protein (Cluster 8 and 9). (C) Both A2-J and A2-O induced the expression of suberin biosynthetic enzymes (Cluster 7). LogFC values of the genes compared to mock treatment at 1, 2, and 3 DPI are shown. The heatmap gray to indicate values with no statistically significant differences at | logFC| ≥ 1 and FDR ≤ 0.01.
Infection With A2-O Induces Genes Related to Defense Responses
In Cluster 8 (genes specifically up-regulated by A2-O), GO terms that were significantly enriched were related to defense responses, including “defense response to fungus (GO: 0050832)”, “defense response to bacterium (GO: 0042742)”, “killing of cells of other organism (GO:0031640)”, and “regulation of salicylic acid biosynthetic process (GO: 0080142)”. In addition, GO terms involved in lignin biosynthesis, including “lignin biosynthetic process (GO:0009809)” was also overrepresented in Cluster 8 (Table 4 and Supplementary Table 8). In Cluster 9, the significantly enriched GO terms were related to biosynthesis of isoprenoids (“sesquiterpene biosynthetic process (GO: 0051762)”, “terpenoid biosynthetic process (GO: 0016114)”, and “farnesyl diphosphate catabolic process (GO:0045339)”) (Table 4 and Supplementary Table 8).
Table 4. GO enrichment analysis of Cluster 8 and 9 genes up-regulated specifically by infection with M. arenaria A2-Oa.
We also found that the genes that are highly expressed after infection with A2-O in Cluster 8 and 9 include (1) defense-related genes encoding chitinase, β-1,3-glucanase, and serine protease inhibitor, (2) sesquiterpene synthase, (3) fatty acid desaturase 2 (FAD2), (4) ferulic acid 5-hydroxylase (F5H) which is involved in lignin biosynthesis, (5) berberine bridge enzyme (BBE)-like protein, which is involved in oxidation of cinnamyl alcohol (Figure 5B). Fatty acids are major and essential components of all plant cells and are also precursors for a variety of plant metabolites, including signaling molecules and phytoalexins (Ohlrogge and Browse, 1995; Lim et al., 2017). FAD2 encodes Δ12-desaturase that catalyzes the conversion of oleic acid (C18:1) to linoleic acid (C18:2) (Ohlrogge and Browse, 1995). The Arabidopsis genome has only a single FAD2 gene (AT3G12120), but most other plant species carry multiple FAD2 homologs (Cao et al., 2013; Lee et al., 2020). The duplication of FAD2 genes in plants would have enabled the functional diversification of these enzymes, leading to divergent catalytic activities and the synthesis of novel metabolites. For example, recent studies have shown that tomato has non-canonical FAD2 family proteins that lack Δ12-desaturase activity (Jeon et al., 2020; Lee et al., 2020). In particular, ACET1a/b (Solyc12g100240 and Solyc12g100260) and FAD2-9 (Solyc12g100250) are non-canonical FAD2 involved in the biosynthesis pathway from linoleic acid to a phytoalexin, falcarindiol (Jeon et al., 2020). Falcarindiol has not only anti-bacterial and anti-fungal activities but also nematicidal activity to M. incognita and pinewood nematode Bursaphelenchus xylophilus (Liu et al., 2016). Infection with A2-O rapidly induced the expression of ACET1a/b and FAD2-9 (Figure 5B), suggesting that infection with A2-O rapidly activates a biosynthesis pathway similar to the falcarindiol pathway, but the production of falcarindiol by S. torvum needs to be experimentally confirmed in the future.
Infection With A2-O Induces Lignin Accumulation in S. torvum
The GO enrichment analysis of Cluster 8 revealed that “lignin biosynthetic process (GO:0009809)” was significantly enriched (Table 4 and Supplementary Table 8), and that the expression of F5H in Cluster 8 was very high and specifically induced after infection with A2-O (Figure 5B). These results suggest that the infection with A2-O transcriptionally activates lignin biosynthesis. Lignin is a phenylpropanoid polymer that is deposited predominantly in the secondary cell wall, making the cell wall rigid and impervious to water (Vanholme et al., 2010). Lignin polymer is synthesized via oxidative combinational coupling of lignin monomers (or monolignols), namely p-coumaryl alcohol, sinapyl alcohol, and coniferyl alcohol. The lignin subunits constituted by these monolignols are p-hydroxyphenyl (H), syringyl (S), and guaiacyl (G) groups, respectively. All of the monolignols are synthesized from phenylalanine through the general phenylpropanoid and monolignol-specific pathways (Vanholme et al., 2012; Figure 6A). Normally, lignin deposition occurs in the root endodermis of the differentiation zone and constitutes the Casparian strip, which functions as a physical barrier that prevents free diffusion of solutes and ions between the xylem and the soil (Robbins et al., 2014). However, biosynthesis and deposition of lignin can be induced in response to biotic stresses (Miedes et al., 2014; Mutuku et al., 2019), which prompted a closer examination of the expression patterns of genes involved in the lignin biosynthetic pathway whose expression was up-regulated by infection with either A2-J or A2-O (logFC ≥ 1, FDR ≤ 0.01) (Figure 6B). Infection with A2-O induced the expression of genes encoding phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumaroyl-CoA ligase (4CL), p-hydroxycinnamoyl-CoA:shikimate p-hydroxycinnamoyl transferase (HCT), caffeoyl-CoA O-methyltransferase (CCoAOMT), cinnamoyl-CoA reductase (CCR), F5H, caffeic acid O-methyltransferase (COMT), and cinnamyl alcohol dehydrogenase (CAD). Infection with A2-J also induced the expression of some of these genes, but to a much lesser extent than A2-O (Figure 6B). Phloroglucinol staining of infected roots allowed us to visualize the intensity and location of lignin accumulation (Figure 7). Infection with A2-O, but not with mock treatment, induced ectopic accumulation of lignin in root tips. With A2-J infection, the area of the root proximal to gall tissue was very slightly stained with phloroglucinol, but the gall itself had little or no detectable phloroglucinol staining. These differences in lignin staining intensity may reflect differences in the expression of lignin biosynthetic genes after infection with A2-J and A2-O (Figure 5B).
Figure 6. Infection with A2-O induced the expression of lignin biosynthetic genes. (A) Overview of the lignin biosynthesis pathway. PAL, phenylalanine ammonia-lyase; C4H, cinnamate 4-hydroxylase; 4CL, 4-coumaroyl-CoA ligase; HCT, p-hydroxycinnamoyl-CoA:shikimate p-hydroxycinnamoyl transferase; C3H, p-coumarate 3-hydroxylase; CCoAOMT, caffeoyl-CoA O-methyltransferase; CCR, cinnamoyl-CoA reductase; F5H, ferulate 5-hydroxylase; COMT, caffeic acid O-methyltransferase; CAD, cinnamyl alcohol dehydrogenase. (B) Expression patterns of genes involved in the lignin biosynthetic pathway whose expression was significantly up-regulated by infection with A2-J or A2-O (logFC ≥ 1, FDR ≤ 0.01). LogFC values of the genes compared to mock treatment at 1, 2, and 3 DPI are shown. The heatmap uses gray to indicate values with no statistically significant differences at | logFC| ≥ 1 and FDR ≤ 0.01.
Figure 7. Infection with A2-O but not with A2-J induced lignin deposition. Phloroglucinol staining of lignin (red) in S. torvum root tips treated with SDW (mock) or infected with A2-J or A2-O for three days. Scale bars indicate 0.5 mm. This experiment was performed four times with the same results.
Infection With A2-O or A2-J Induces the Expression of Genes Related to Suberin Biosynthesis
There was no enrichment of specific GO terms in Cluster 7 (up-regulated by the infection with A2-J or A2-O) (Supplementary Table 8). However, we found that both A2-J and A2-O strongly activate the expression of suberin biosynthetic genes, including aliphatic suberin feruloyl transferase (ASFT), cytochrome P450 86A1 (CYP86A1), cytochrome P450 86B1 (CYP86B1), glycerol-3-phosphate acyltransferase 5 (GPAT5), and β−ketoacyl−CoA synthase (KCS) (Figure 5C), but A2-O induced slightly higher expression of these genes than A2-J. Suberin is a cell wall component that restricts water loss, nutrient elution, and pathogen infection (Dean and Kolattukudy, 1976; Lulai et al., 2008; Leide et al., 2012; Froschel et al., 2020). It is normally deposited in the cell walls of endodermal cells, but not in the root tip (Barberon, 2017), and several reports showed that suberin synthesis is induced in wounded tissues (Dean and Kolattukudy, 1976; Bernards and Lewis, 1998). Considering that both A2-J and A2-O induced the expression of suberin biosynthetic genes, the expression of these genes was possibly up-regulated by a generalized wounding signal during nematode infection. We therefore investigated the expression levels of the wound-responsive genes Arabidopsis thaliana activation factor 2 (ATAF2), which encodes a wound-responsive NAC transcription factor (Collinge and Boller, 2001; Wang and Culver, 2012), and pathogen-related 4 (Stanford et al., 1989; Marhavý et al., 2019; Supplementary Figure 5). Not surprisingly, the expression patterns of wound-responsive genes were similar to those of the suberin biosynthetic genes. These results suggest that infection with A2-O may cause more wounding and damage than A2-J, but that the wounding signal generated during the earliest events in nematode infection, whether a resistant or susceptible reaction occurs later, may induce the expression of suberin biosynthetic genes.
Considering the suberin biosynthetic gene transcription data, we predicted that infection with A2-O or A2-J would induce the accumulation of suberin. To test our prediction, we examined and measured the chemical composition of aliphatic suberin at the nematode infection site, but total aliphatic suberin content did not increase after A2-J or A2-O infection (Welch’ s t-test with Bonferroni correction, p ≤ 0.05) (Figure 8A). However, the abundance of single monomers, ω-hydroxy acid (ω-OH-acid) C16 and α,ω-diacid C16, were significantly higher after infection with A2-O than mock treatment (p ≤ 0.05, Welch’s t-test followed by Bonferroni correction) (Figure 8B).
Figure 8. Infection with A2-O changed suberin monomer composition in the root tips. (A) Total aliphatic suberin and (B) suberin monomer composition in root tips treated with SDW (mock) or infected with A2-J or A2-O for four days. GC-MS was performed to measure the content of aliphatic suberin monomers in root tips. Bars indicate means ± SE of four biological replicates. Asterisks indicate significant differences compared with mock treatment (Welch’s t-test with Bonferroni correction, * p ≤ 0.05). There was no statistically significant difference at p ≤ 0.05 in the total amount of aliphatic suberin between samples (Welch’s t-test with Bonferroni correction).
Comparative Transcriptome Analyses for Host Defense and Susceptible Responses at the Early Stages of Nematode Infection
In this study we established an experimental system using a single cultivar of S. torvum and two pathotypes of RKNs M. arenaria A2-J, which is able to infect and establish a parasitic relationship with the plant, and A2-O, which induces a resistance response. To the best of our knowledge, this is the first comparative RNA-seq analyses in a single cultivar with virulent or avirulent nematodes. Using this experimental system, we were able to catalog changes in gene expression at the very early stage of infection (1–3 DPI) with a high degree of sensitivity compared to previous studies (Bagnaresi et al., 2013; Postnikova et al., 2015; Xing et al., 2017; Shukla et al., 2018). Because there are clear morphological differences in root tips infected with A2-J and A2-O within 4 DPI, the sum of host responses by this stage of infection likely has determined the outcome of infection. Interestingly, at 1 DPI, A2-O induced the expression of genes encoding class III peroxidases, fatty acid desaturases, and enzymes involved in isoprenoid biosynthesis (Figures 3B, 5B and Supplementary Figure 3), whereas A2-J had not induced any statistically significant changes in the expression of defense response genes (Figure 3A). These results suggest that A2-J initially evades recognition by S. torvum and/or actively suppresses the induction of transcriptional changes in the host. Our results also show that A2-J induces the expression of genes associated with susceptible responses that are related to gall formation at around 2–3 DPI (Supplementary Figure 4). The delay in expression of susceptibility-associated genes is consistent with the length of time it takes for nematodes to migrate through the plant’s vascular system to their target cells.
Sampling only infected root tips reduced background plant transcripts and enabled us to detect dynamic changes in gene expression that are directly associated with the establishment of a parasitic relationship or with active host defense. The range of calculated logFC values of root tip genes expressed at 3 DPI with A2-O was from −8 to 12.8 (that is, from 1/256 to 7131 times), which is quite broad for very early stages of infection. For example, in previous transcriptome analyses using S. torvum during infection with an avirulent isolate of M. incognita, chitinase expression in whole roots was only 2–4 times higher upon infection than in uninfected roots (Bagnaresi et al., 2013). In contrast, our analyses showed that the expression of chitinases in the root tips infected with A2-O at 3 DPI was several hundred times higher than in the mock treatment, clearly demonstrating the high sensitivity of this method for host-parasite transcriptomics (Figure 5B). This level of sensitivity allowed us to identify genes that had not previously been connected with resistance or susceptibility in this system. In addition, comparative transcriptome analyses of S. torvum combined with SOM clustering can be used to pinpoint genes that are associated with susceptibility or resistance.
Infection With A2-J Induces the Expression of Genes Related to Cell Wall Modification and Spermidine Synthase
A2-J specifically induced the expression of cell wall modification enzymes such as cellulose synthase-like protein, xyloglucan endotransglucosylase/hydrolase protein, and a non-catalytic subunit of the polygalacturonase isozyme (Table 3 and Supplementary Figure 3A). The cell walls of giant cells require thickening and loosening to allow cell expansion that increases the surface area of the plasma membrane, and to support nutrient uptake by the nematodes (Bartlem et al., 2014; Bozbuga et al., 2018; Meidani et al., 2019). The cell walls of giant cells and syncytial feeding sites induced by CNs contain high-ester pectic homogalacturonan, xyloglucan, and pectic arabinan (Davies et al., 2012; Zhang L. et al., 2017), suggesting that these polysaccharides are responsible for the flexible properties of feeding site cell walls. Our transcriptome results confirm the importance of cell wall modification enzymes during gall formation.
In addition to cell wall modification enzymes, infection with A2-J specifically induced the expression of spermidine synthase (Figure 5A). Interestingly, a virulence effector protein secreted from CNs, 10A06, functions through its interaction with Arabidopsis spermidine synthase 2 (Hewezi et al., 2010) by increasing spermidine concentrations, subsequently increasing polyamine oxidase activities. An increase in polyamine oxidase activity results in the induction of cellular antioxidant machinery in syncytia and disruption of SA-mediated defense signaling. Although there is no clear homolog of 10A06 in RKNs such as M. incognita, it is possible that there is a functional ortholog of 10A06 that up-regulates the expression of spermidine synthase. Potentially, transcriptional activation of the spermidine synthase gene by PPNs may turn out to be a common strategy for suppressing plant immunity.
Infection With A2-O Up-Regulates Genes Related to Sesquiterpene Biosynthesis and Fatty Acid Metabolism
Infection with A2-O, but not with A2-J, strongly induced genes encoding sesquiterpene synthases (Table 4 and Figure 5B). High levels of expression of these genes occurs at the very early stages of infection (1 DPI), which sets these genes apart from the other genes induced by A2-O infection (Figures 3B, 5B), suggesting the importance of sesquiterpenes as an early line of defense against PPNs. Capsidiol, a sesquiterpene, is the major phytoalexin produced in the Solanaceae plants Nicotiana spp. and Capsicum spp. in response to fungal and bacterial infection (Grosskinsky et al., 2011; Song et al., 2019). Capsidiol is toxic to many oospore and fungal pathogens, such as Phytophthora capsici and Botrytis cinerea (Stoessl et al., 1972; Ward et al., 1974), and suppresses the mobility of false root-knot nematode Nacobbus aberrans (Godinez-Vidal et al., 2010). We found that A2-O induces the expression of CYP71D7-like protein, the closest homologue to CYP71D20 from Nicotiana benthamiana, which converts 5-epiaristolochene to capsidiol (Ralston et al., 2001). A2-O also induces the expression of genes encoding 5-epiaristolochene synthase, which converts farnesyl diphosphate to 5-epiaristolochene (Vogeli and Chappell, 1988; Facchini and Chappell, 1992). Although capsidiol production is not common in Solanum spp., S. torvum may produce similar sesquiterpene derivatives that are toxic to PPNs. Another sesquiterpene that may be involved in resistance is the phytoalexin solavetivone, because A2-O induces the expression of genes encoding vetispiradiene synthase and premnaspirodiene oxygenase, which sequentially convert farnesyl diphosphate to solavetivone via vetispiradiene (Takahashi et al., 2007). Although the nematicidal activity of solavetivone has not yet been reported, genetic analyses indicate that the production of solavetivone is associated with resistance in potato against CN Globodera rostochiensis (Desjardins et al., 1997). It would be interesting to test the nematicidal activity of these sesquiterpenes and their production after A2-O infection in S. torvum.
We also found that infection with A2-O induces the expression of genes encoding the non-canonical FAD2 proteins ACET1a/b (Solyc12g100240 and Solyc12g100260) and FAD2-9 (Solyc12g100250), which are involved in the biosynthesis of falcarindiol, a modified fatty acid found in a variety of plants, including Solanaceae. It is toxic to RKNs and pinewood nematodes (Liu et al., 2016), but it is not known if falcarindiol is produced by S. torvum. In any case, the higher expression of FAD2 after A2-O infection suggests the importance of fatty acid metabolism in plant immunity against PPNs. We consistently found that A2-O infection results in the accumulation of linoleic acid (C18:2), while infection with A2-J reduced linoleic acid (C18:2) and palmitic acid (C16:0) (Supplementary Figure 6). Increases in linoleic acid elicits resistance to the fungal pathogen Colletotrichum gloeosporioides in avocado (Madi et al., 2003), and to B. cinerea in bean plants (Ongena et al., 2004). The importance of fatty acids in resistance is further supported by the fact that PPNs secrete fatty-acid- and retinol-binding family proteins as a mechanism to increase susceptibility (Prior et al., 2001; Iberkleid et al., 2013, 2015). Thus, the battle over fatty acid synthesis is likely to be important in S. torvum-PPN interactions.
Lignin Accumulation as a Defense Against PPNs
Because PPNs penetrate the cell wall and migrate within roots, reinforcement of cell walls by lignin accumulation has been implicated as an effective defense response to PPNs (Holbein et al., 2016; Sato et al., 2019). In fact, several studies showed that PPN infection induces more extensive lignin accumulation in resistant plants than in susceptible plants. For example, Veronico et al. (2018) performed a histochemical analysis of lignin in the roots of susceptible and resistant tomatoes infected with M. incognita, and found that accumulation of higher lignin levels in the root tissues (e.g., cortical cells) of resistant tomato than in susceptible varieties. Similarly, upon infection with the cereal CN Heterodera avenae, a resistant wheat cultivar gives a strong lignin accumulation response in the walls of cells affected by nematode infection, but a susceptible cultivar responded with only minor additional lignification, and only in cell walls mechanically damaged during nematode penetration (Andres et al., 2001). Since these studies used resistant and susceptible lines, it is certainly possible that differences in lignin accumulation are due to multiple differences in genetic background, and thus there would not necessarily be a direct relationship to resistance. In the present study, by using a single host cultivar with virulent and avirulent nematodes, we showed that susceptible and resistant responses are correlated with the levels of lignin. Infection with A2-O significantly up-regulated genes involved in lignin biosynthesis, and induced lignin accumulation near the root tip (Figures 6B, 7). In contrast, infection with A2-J induced lignin biosynthetic genes to some limited extent that did not result in additional lignin accumulation. We also found that infection with A2-O significantly up-regulated the expression of class III peroxidases, some of which are thought to be involved in the polymerization of lignin (Zacheo et al., 1993, 1997; Marjamaa et al., 2009). The importance of lignin accumulation in PPN resistance is also supported by the fact that plant defense inducers such as β-aminobutyric acid, thiamine, sclareol, and benzothiadiazole induce lignin biosynthesis and enhance resistance against PPNs (Fujimoto et al., 2015; Ji et al., 2015; Huang et al., 2016; Veronico et al., 2018). Thus, ectopic lignification may result in cell wall reinforcement that restricts RKN invasion and migration.
In contrast to the resistance interaction with A2-O, infection with A2-J did not induce lignin accumulation. One potential mechanism to explain this is that A2-J may alter the metabolic flux of lignin biosynthesis by specifically inducing the expression of the gene encoding chalcone synthase, which converts p-coumaroyl-CoA into naringenin chalcone (Figure 5A). The biosynthetic pathways leading to lignin and flavonoids diverge at the common intermediate p-coumaroyl CoA, so higher chalcone synthase expression should decrease lignin content. Indeed, silencing of the chalcone synthase gene increases lignin deposition in wheat (Eloy et al., 2017).
Activation of the phenylpropanoid biosynthesis pathway may also be involved in the production of anti-PPN secondary metabolites (Khanam et al., 2018). For example, trans-cinnamaldehyde is a phenylpropanoid compound that is highly toxic to RKNs and pinewood nematodes (Oka, 2001; Kong et al., 2007). Cinnamaldehyde can be synthesized from L-phenylalanine, and its biosynthesis requires enzymatic reactions mediated by PAL, 4CL, and CCR (Bang et al., 2016). Expression of genes encoding these enzymes was specifically induced in A2-O-infected S. torvum (Figure 6B). We also observed a drastic induction of a BBE-like protein homologous to AtBBE-like 13 and AtBBE-like 15 (Figure 5B), which can convert cinnamyl alcohol to cinnamaldehyde in Arabidopsis (Daniel et al., 2015). Thus, S. torvum may activate the production of cinnamaldehyde to suppress PPN infection.
Both A2-J and A2-O Induce the Expression of Suberin Synthetic Genes
Suberin is a lipid-phenolic biopolyester deposited in the cell walls of the root endodermis. Suberin lamellae and the lignin-based Casparian strip in the endodermis form the apoplastic diffusion barrier between vascular tissues and outer ground tissues to enforce selective absorption of water and nutrients (Alassimone et al., 2012; Vishwanath et al., 2015). The suberized endodermis may also play a role in defense against pathogens as the last line of defense before pathogens invade the vasculature (Thomas et al., 2007). Moreover, several reports showed that fungal or viral attack induces the deposition of suberin in the walls of cells in and around the site of penetration (Tippett and Hill, 1984; Kolattukudy and Espelie, 1989). In contrast, the virulent vascular pathogen, Verticillium longisporum reduces suberin deposition by suppressing the translation of genes involved in suberin biosynthesis (Froschel et al., 2020). These results may indicate the general importance of suberin in plant immunity. Infection with A2-J or A2-O induced the expression of genes involved in suberin biosynthesis at the root tip, with A2-O induction being slightly higher (Figure 5C). The total aliphatic suberin content tended to increase upon infection with A2-J and A2-O, but there was no statistically significant difference due to the large variation of the data (Figure 8A). The large variation in the data may be due to variation in the number of infected root tips, nematode aggressiveness, and/or the severity of wounding, which can vary by experiment and cause disparity in the proportions of suberized areas among the fragile root tips. Because infection with RKN M. incognita induces the expression of suberin biosynthetic genes and a patchy suberin accumulation pattern within the cell layer surrounding the infection site in later stages of infection (10 DPI) in susceptible Arabidopsis plants (Holbein et al., 2019), it is possible that the synthesis of large amounts of suberin takes longer in S. torvum that what was allowed for in this assay. Although there was no statistical difference in total aliphatic suberin content, infection with A2-O increased the accumulation of suberin monomers such as ω-OH-acid C16 and α, ω-diacid C16 (Figure 8). These results suggest that de novo suberin synthesis is induced at the site of infection. Considering that A2-J also induces the expression of suberin biosynthesis genes (Figure 5C), it is possible that de novo suberin synthesis is induced in both incompatible and compatible interactions between plants and RKNs. This would be consistent with the fact that infection with RKN M. incognita induces the expression of GPAT5, a suberin biosynthetic gene, as well as the patchy suberin accumulation pattern within the cell layer surrounding infection sites in susceptible Arabidopsis plants (Holbein et al., 2019). Because the expression pattern of the suberin biosynthetic genes was similar to that of wound-inducible genes (Supplementary Figure 5), the induction of suberin biosynthesis could be triggered by wounding during nematode infection, but this need to be confirmed experimentally. In addition, it is important to clarify the biological relevance of ectopic suberin accumulation at the site of infection in the resistance response against PPNs in the future. We may be able to assess the individual contribution of ectopic suberin accumulation to the resistance response by selectively expressing a suberin-degrading enzyme CDEF1 (cuticle destructing factor 1) only at the site of infection (Takahashi et al., 2010). However, transformation of S. torvum is technically challenging at the moment.
A Model for Resistant and Susceptible Responses to RKN Infection
In this study, we were able to examine detail the transcriptional reprogramming of S. torvum in response to infection with virulent and avirulent RKNs at early infection stages by comparative transcriptome analyses (Figure 9). Shortly after infection with A2-O, S. torvum rapidly induces the expression of genes encoding fatty acid desaturases and sesquiterpene synthases, which may act as the first, or at least an early line of defense. Infection with A2-O also induces the expression of SA, JA, and ET marker genes, and defense-related genes (WRKY transcription factors, chitinases, and β-1,3-glucanases), as well as the expression of lignin biosynthetic genes, leading to lignin accumulation in root tips. In contrast, infection with A2-J fails to alter gene expression at 1 DPI, but it does specifically induce the expression of genes encoding chalcone synthase and spermidine synthase, which might suppress lignin synthesis and the SA-mediated defense response (Hewezi et al., 2010). A2-J infection also induces the expression of genes related to cell wall modification and transmembrane transport, which may be important for the development and maturation of giant cells, and which in turn support nematode feeding, and are necessary to complete the parasitic life cycle. Interestingly, nematode invasion generally induces the expression of suberin synthetic genes, suggesting that the wounds caused by nematode entry and movement trigger suberin gene transcriptional activation. In the future, it will be important to independently assess the contribution of these genes to immunity and susceptibility against RKNs. It is also important to understand the molecular mechanisms and signal transduction pathways that decide the fate of infected plants and their nematode parasites, but will likely include whether or not plants have receptors that can recognize PPNs, and can transduce those signals to trigger defense responses, and/or whether or not PPNs produce virulence effectors that inhibit recognition, disrupt signal transduction, or suppress the defense response. As the next step, it is important to identify the specific effectors of A2-J and A2-O by comparative genomics, transcriptomics, and effectoromics, and to characterize their functions at the molecular level.
Figure 9. Schematic model of S. torvum responses to M. arenaria A2-O and A2-J. In the S. torvum/M. arenaria model, shortly after infection with avirulent RKN, plant genes encoding fatty acid desaturases and sesquiterpene synthases are induced to act as an early line of defense. Avirulent RKN infection also induces lignin accumulation in root tip, which may result in cell wall reinforcement that restricts RKN invasion and migration. Infection with virulent RKN results in a delayed gene expression response by specifically inducing chalcone synthase and spermidine synthase, which suppress multiple defense responses. Virulent RKN infection also results in induced expression of cell wall modification and transmembrane transport genes whose products contribute to giant cell formation for nematode feeding and completion of the life cycle.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) are as follows: NCBI (https://www.ncbi.nlm.nih.gov), PRJNA681478; Figshare (https://figshare.com), https://doi.org/10.6084/m9.figshare.14054354.
YK and KSh supervised the research. YK designed the experiments. YK and KSa carried out histopathological observation. YK and TU propagate the nematodes and YK collected samples for RNA-seq, lipid, and suberin analyses. NM generated RNA-seq libraries. TS carried out RNA-seq. YI performed quality filtering of raw RNA-seq reads. KSa carried out bioinformatic analyses, including assembly, annotation, gene expression analyses, and GO enrichment tests. PG helped in bioinformatic analyses. TB, KY, and SSh provided RKN genome assemblies. YK and KSa performed lignin staining. YS-S and HO performed lipid quantification. JH, RF, SSi, and FG carried out suberin monomer quantification. KSa, YK, and KSh wrote the manuscript. All authors read and approved the final manuscript.
This research was financially supported by Cabinet Office, Government of Japan, cross-ministerial Strategic Innovation Promotion Program (SIP), “Technologies for Creating Next-Generation Agriculture, Forestry and Fisheries” (funding agency, Bio-oriented Technology Research Advancement Institution, NARO) (to TU and YK), an NIBB Collaborative Research Program (20–421) (to YK), MEXT/JSPS KAKENHI Grant Number, JP16H06186 (to YK), JP16KT0037, JP19H02962, and JP20H02994 (to TU and YK), JP17H06172 and JP20H05909 (to KSh), and Grant-in-Aid for JSPS Fellows Grant Number JP19J00655 (to KSa).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank all members of the Shirasu lab for intensive and helpful discussions. We thank Akiko Ueno, Naomi Watanabe, Mamiko Kouzai, Yoko Nagai, and Kanako Hori for their support of this project. Hideaki Iwahori at Ryukoku University for sharing M arenaria A2-O. We thank Mitsuyasu Hasebe, Shoko Ohi, Miwako Matsumoto (NIBB), and Tomoko Shibata for the PacBio sequencing and the Data Integration and Analysis Facility at NIBB for computer resources. We also thank B. Favery for sharing the BABB cleaning protocols.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.680151/full#supplementary-material
CN, cyst nematode; DPI, day-post inoculation; RKN, root-knot nematode; PPN, plant-parasitic nematode.
- ^ http://arthropods.eugenes.org/genes2/about/EvidentialGene_trassembly_pipe.html (version 2017.12.21).
Abad, P., Castagnone-Sereno, P., Rosso, M. N., Engler, J. D. A., and Favery, B. (2009). “Invasion, Feeding and Development,” in Root-knot nematodes, eds R. N. Perry, M. Moens, and J. L. Starr, (Cambridge: CABI International Press), 163–181.
Alassimone, J., Roppolo, D., Geldner, N., and Vermeer, J. E. (2012). The endodermis — development and differentiation of the plant’s inner skin. Protoplasma 249, 433–443. doi: 10.1007/s00709-011-0302-5
Bagnaresi, P., Sala, T., Irdani, T., Scotto, C., Lamontanara, A., Beretta, M., et al. (2013). responses to the root-knot nematode Meloidogyne incognita. BMC Genomics 14:540. doi: 10.1186/1471-2164-14-540
Bali, S., Vining, K., Gleason, C., Majtahedi, H., Brown, C. R., and Sathuvalli, V. (2019). Transcriptome profiling of resistance response to Meloidogyne chitwoodi introgressed from wild species Solanum bulbocastanum into cultivated potato. BMC Genomics 20:907. doi: 10.1186/s12864-019-6257-1
Berardini, T. Z., Reiser, L., Li, D., Mezheritsky, Y., Muller, R., Strait, E., et al. (2015). The Arabidopsis information resource: Making and mining the “gold standard” annotated reference plant genome. Genesis 53, 474–485. doi: 10.1002/dvg.22877
Bohlmann, H. (2015). “Chapter Two - Introductory Chapter on the Basic Biology of Cyst Nematodes,” in Plant nematode interactions: a view on compatible interrelationships, vol 73, Advances in Botanical Research, eds C. Escobar and C. Fenoll, (Oxford, UK: Elsevier), 33–60.
Boutrot, F., and Zipfel, C. (2017). Function, discovery, and exploitation of plant pattern recognition receptors for broad-spectrum disease resistance. Annu. Rev. Phytopathol. 55, 257–286. doi: 10.1146/annurev-phyto-080614-120106
Bozbuga, R., Lilley, C. J., Knox, J. P., and Urwin, P. E. (2018). Host-specific signatures of the cell wall changes induced by the plant parasitic nematode. Meloid. Incognita. Sci. Rep. 8:17302. doi: 10.1038/s41598-018-35529-7
Brennan, R. J. B., Glaze-Corcoran, S., Wick, R., and Hashemi, M. (2020). Biofumigation: An alternative strategy for the control of plant parasitic nematodes. J. Integr. Agric. 19, 1680–1690. doi: 10.1016/S2095-3119(19)62817-0
Brown, D. M., Zeef, L. A., Ellis, J., Goodacre, R., and Turner, S. R. (2005). Identification of novel genes in Arabidopsis involved in secondary cell wall formation using expression profiling and reverse genetics. Plant Cell 17, 2281–2295. doi: 10.1105/tpc.105.031542
Cabrera, J., Olmo, R., Ruiz-Ferrer, V., Abreu, I., Hermans, C., Martinez-Argudo, I., et al. (2018). A Phenotyping Method of Giant Cells from Root-Knot Nematode Feeding Sites by Confocal Microscopy Highlights a Role for CHITINASE-LIKE 1 in Arabidopsis. Int. J. Mol. Sci. 19:429. doi: 10.3390/ijms19020429
Cao, S., Zhou, X. R., Wood, C. C., Green, A. G., Singh, S. P., Liu, L., et al. (2013). A large and functionally diverse family of Fad2 genes in safflower (Carthamus tinctorius L.). BMC Plant Biol. 13:5. doi: 10.1186/1471-2229-13-5
Chen, R., Li, H., Zhang, L., Zhang, J., Xiao, J., and Ye, Z. (2007). CaMi, a root-knot nematode resistance gene from hot pepper (Capsium annuum L.) confers nematode resistance in tomato. Plant Cell Rep. 26, 895–905. doi: 10.1007/s00299-007-0304-0
Chitwood, D. H., Maloof, J. N., and Sinha, N. R. (2013). Dynamic transcriptomic profiles between tomato and a wild relative reflect distinct developmental architectures. Plant Physiol. 162, 537–552. doi: 10.1104/pp.112.213546
Claverie, M., Dirlewanger, E., Bosselut, N., Van Ghelder, C., Voisin, R., Kleinhentz, M., et al. (2011). The Ma gene for complete-spectrum resistance to Meloidogyne species in Prunus is a TNL with a huge repeated C-terminal post-LRR region. Plant Physiol. 156, 779–792. doi: 10.1104/pp.111.176230
Collinge, M., and Boller, T. (2001). Differential induction of two potato genes, Stprx2 and StNAC, in response to infection by Phytophthora infestans and to wounding. Plant Mol. Biol. 46, 521–529. doi: 10.1023/a:1010639225091
Daniel, B., Pavkov-Keller, T., Steiner, B., Dordic, A., Gutmann, A., Nidetzky, B., et al. (2015). Oxidation of Monolignols by Members of the Berberine Bridge Enzyme Family Suggests a Role in Plant Cell Wall Metabolism. J. Biol. Chem. 290, 18770–18781. doi: 10.1074/jbc.M115.659631
Davies, L. J., Lilley, C. J., Paul Knox, J., and Urwin, P. E. (2012). Syncytia formed by adult female Heterodera schachtii in Arabidopsis thaliana roots have a distinct cell wall molecular architecture. New Phytol. 196, 238–246. doi: 10.1111/j.1469-8137.2012.04238.x
Desjardins, A. E., McCormick, S. P., Plaisted, R. L., and Brodie, B. B. (1997). Association between solavetivone production and resistance to Globodera rostochiensis in potato. J. Agric. Food. Chem. 45, 2322–2326.
Du, C., Jiang, J., Zhang, H., Zhao, T., Yang, H., Zhang, D., et al. (2020). Transcriptomic profiling of Solanum peruvianum LA3858 revealed a Mi-3-mediated hypersensitive response to Meloidogyne incognita. BMC Genomics 21:250. doi: 10.1186/s12864-020-6654-5
Elashry, A., Okumoto, S., Siddique, S., Koch, W., Kreil, D. P., and Bohlmann, H. (2013). The AAP gene family for amino acid permeases contributes to development of the cyst nematode Heterodera schachtii in roots of Arabidopsis. Plant Physiol. Biochem. 70, 379–386. doi: 10.1016/j.plaphy.2013.05.016
Eloy, N. B., Voorend, W., Lan, W., Saleme, M. L., Cesarino, I., Vanholme, R., et al. (2017). Silencing CHALCONE SYNTHASE in Maize Impedes the Incorporation of Tricin into Lignin and Increases Lignin Content. Plant Physiol. 173, 998–1016. doi: 10.1104/pp.16.01108
Ernst, K., Kumar, A., Kriseleit, D., Kloos, D. U., Phillips, M. S., and Ganal, M. W. (2002). The broad-spectrum potato cyst nematode resistance gene (Hero) from tomato is the only member of a large gene family of NBS-LRR genes with an unusual amino acid repeat in the LRR region. Plant J. 31, 127–136.
Escobar, C., Barcala, M., Cabrera, J., and Fenoll, C. (2015). “Chapter One - Overview of Root-Knot Nematodes and Giant Cells,” in Plant nematode interactions: a view on compatible interrelationships, vol 73, Advances in Botanical Research, eds C. Escobar and C. Fenoll, (Oxford, UK: Elsevier), 1–32.
Farmer, E. E., Johnson, R. R., and Ryan, C. A. (1992). Regulation of expression of proteinase inhibitor genes by methyl jasmonate and jasmonic Acid. Plant Physiol. 98, 995–1002. doi: 10.1104/pp.98.3.995
Franke, R., Briesen, I., Wojciechowski, T., Faust, A., Yephremov, A., Nawrath, C., et al. (2005). Apoplastic polyesters in Arabidopsis surface tissues – A typical suberin and a particular cutin. Phytochemistry 66, 2643–2658. doi: 10.1016/j.phytochem.2005.09.027
Froschel, C., Komorek, J., Attard, A., Marsell, A., Lopez-Arboleda, W. A., Le Berre, J., et al. (2020). Plant roots employ cell-layer-specific programs to respond to pathogenic and beneficial microbes. Cell Host. Microbe. 11:14. doi: 10.1016/j.chom.2020.11.014
Fujimoto, T., Mizukubo, T., Abe, H., and Seo, S. (2015). Sclareol induces plant resistance to root-knot nematode partially through ethylene-dependent enhancement of lignin accumulation. Mol. Plant Microbe Interact 28, 398–407. doi: 10.1094/MPMI-10-14-0320-R
García-Mendívil, H. A., and Sorribas, F. J. (2019). Fitness cost but no selection for virulence in Meloidogyne incognita after two consecutive crops of eggplant grafted onto Solanum torvum. Meloidogyne Incognita. 11:13092. doi: 10.1111/ppa.13092
Ghaemi, R., Pourjam, E., Safaie, N., Verstraeten, B., Mahmoudi, S. B., and Mehrabi, R. (2020). Molecular insights into the compatible and incompatible interactions between sugar beet and the beet cyst nematode. BMC Plant Biol. 20:483. doi: 10.1186/s12870-020-02706-8
Godinez-Vidal, D., Soto-Hernández, M., Rocha Sosa, M., Lozoya-Gloria, E., Rojas-Martínez, R. I., Guevara Olvera, L., et al. (2010). Capsidiol content in CM-334 pepper roots infected with Nacobbus aberrans and its effect on second stage juveniles. Nematropica 40, 227–237.
Goto, Y., Maki, N., Ichihashi, Y., Kitazawa, D., Igarashi, D., Kadota, Y., et al. (2020). Exogenous Treatment with Glutamate Induces Immune Responses in Arabidopsis. Mol. Plant Microbe Interact. 33, 474–487. doi: 10.1094/MPMI-09-19-0262-R
Gotz, S., Garcia-Gomez, J. M., Terol, J., Williams, T. D., Nagaraj, S. H., Nueda, M. J., et al. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36, 3420–3435. doi: 10.1093/nar/gkn176
Gousset, C., Collonnier, C., Mulya, K., Mariska, I., Rotino, G. L., Besse, P., et al. (2005). as a useful source of resistance against bacterial and fungal diseases for improvement of eggplant (S. melongena L.). Solanum torvum. Plant Sci. 7, 319–327. doi: 10.1016/j.plantsci.2004.07.034
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., and Amit, I. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Grosskinsky, D. K., Naseem, M., Abdelmohsen, U. R., Plickert, N., Engelke, T., Griebel, T., et al. (2011). Cytokinins mediate resistance against Pseudomonas. Plant Physiol. 11, 815–830. doi: 10.1104/pp.111.182931
Haas, B. J., Papanicolaou, A., Yassour, M., Grabherr, M., Blood, P. D., Bowden, J., et al. (2013). De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494–1512. doi: 10.1038/nprot.2013.084
Hammes, U. Z., Schachtman, D. P., Berg, R. H., Nielsen, E., Koch, W., McIntyre, L. M., et al. (2005). Nematode-induced changes of transporter gene expression in Arabidopsis roots. Mol. Plant Microbe Interact. 18, 1247–1257. doi: 10.1094/MPMI-18-1247
Hewezi, T., Howe, P. J., Maier, T. R., Hussey, R. S., Mitchum, M. G., Davis, E. L., et al. (2010). Arabidopsis spermidine synthase is targeted by an effector protein of the cyst nematode Heterodera schachtii. Plant Physiol. 152, 968–984. doi: 10.1104/pp.109.150557
Holbein, J., Franke, R. B., Marhavý, P., Fujita, S., Gorecka, M., Sobczak, M., et al. (2019). Root endodermal barrier system contributes to defence against plant-parasitic cyst and root-knot nematodes. Plant J. 100, 221–236. doi: 10.1111/tpj.14459
Huang, S., Gao, Y., Liu, J., Peng, X., Niu, X., Fei, Z., et al. (2012). Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum. Mol. Genet Genom. 287, 495–513. doi: 10.1007/s00438-012-0696-6
Huang, W. K., Ji, H. L., Gheysen, G., and Kyndt, T. (2016). Thiamine-induced priming against root-knot nematode infection in rice involves lignification and hydrogen peroxide generation. Mol. Plant Pathol. 17, 614–624. doi: 10.1111/mpp.12316
Iberkleid, I., Sela, N., and Brown Miyara, S. (2015). Meloidogyne javanica fatty acid- and retinol-binding protein (Mj-FAR-1) regulates expression of lipid-, cell wall-, stress- and phenylpropanoid-related genes during nematode infection of tomato. BMC Genomics 16:272. doi: 10.1186/s12864-015-1426-3
Iberkleid, I., Vieira, P., de Almeida Engler, J., Firester, K., Spiegel, Y., and Horowitz, S. B. (2013). Fatty acid-and retinol-binding protein, Mj-FAR-1 induces tomato host susceptibility to root-knot nematodes. PLoS One 8:e64586. doi: 10.1371/journal.pone.0064586
Jablonska, B., Ammiraju, J. S., Bhattarai, K. K., Mantelin, S., Martinez, de Ilarduya, O., et al. (2007). The Mi-9 gene from Solanum arcanum conferring heat-stable resistance to root-knot nematodes is a homolog of Mi-1. Plant Physiol. 143, 1044–1054. doi: 10.1104/pp.106.089615
Jeon, J. E., Kim, J. G., Fischer, C. R., Mehta, N., Dufour-Schroif, C., Wemmer, K., et al. (2020). A Pathogen-Responsive Gene Cluster for Highly Modified Fatty Acids in Tomato. Cell 180, 176–187e119. doi: 10.1016/j.cell.2019.11.037
Ji, H., Kyndt, T., He, W., Vanholme, B., and Gheysen, G. (2015). β-aminobutyric acid-induced resistance against root-knot nematodes in rice is based on increased basal defense. Mol. Plant Microbe Interact. 28, 519–533. doi: 10.1094/MPMI-09-14-0260-R
Jones, J. T., Haegeman, A., Danchin, E. G., Gaur, H. S., Helder, J., Jones, M. G., et al. (2013). Top 10 plant-parasitic nematodes in molecular plant pathology. Mol. Plant Pathol. 14, 946–961. doi: 10.1111/mpp.12057
Jones, P., Binns, D., Chang, H. Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kaloshian, I., Desmond, O. J., and Atamian, H. S. (2011). “Disease Resistance-Genes and Defense Responses During Incompatible Interactions,” in Genomics and Molecular Genetics of Plant-Nematode Interactions, eds J. Jones, G. Gheysen, and C. Fenoll, (Netherlands: Springer), 309–324.
Khanam, S., Bauters, L., Singh, R. R., Verbeek, R., Haeck, A., Sultan, S. M. D., et al. (2018). Mechanisms of resistance in the rice cultivar Manikpukha to the rice stem nematode Ditylenchus angustus. Mol. Plant Pathol. 19, 1391–1402. doi: 10.1111/mpp.12622
Kim, D., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4
Kong, J. O., Lee, S. M., Moon, Y. S., Lee, S. G., and Ahn, Y. J. (2007). Nematicidal Activity of Cassia and Cinnamon Oil Compounds and Related Compounds toward Bursaphelenchus xylophilus (Nematoda: Parasitaphelenchidae). J. Nematol. 39, 31–36.
Kumar, R., Ichihashi, Y., Kimura, S., Chitwood, D. H., Headland, L. R., Peng, J., et al. (2012). A High-Throughput Method for Illumina RNA-Seq Library Preparation. Front. Plant Sci. 3:202. doi: 10.3389/fpls.2012.00202
Lee, M. W., Padilla, C. S., Gupta, C., Galla, A., Pereira, A., Li, J., et al. (2020). The FATTY ACID DESATURASE2 Family in Tomato Contributes to Primary Metabolism and Stress Responses. Plant Physiol. 182, 1083–1099. doi: 10.1104/pp.19.00487
Leide, J., Hildebrandt, U., Hartung, W., Riederer, M., and Vogg, G. (2012). Abscisic acid mediates the formation of a suberized stem scar tissue in tomato fruits. New Phytol. 194, 402–415. doi: 10.1111/j.1469-8137.2011.04047.x
Li, J. B., Luan, Y. S., and Liu, Z. (2015a). Overexpression of SpWRKY1 promotes resistance to Phytophthora nicotianae and tolerance to salt and drought stress in transgenic tobacco. Physiol. Plant. 155, 248–266. doi: 10.1111/ppl.12315
Li, J. B., Luan, Y. S., and Liu, Z. (2015b). SpWRKY1 mediates resistance to Phytophthora infestans and tolerance to salt and drought stress by modulating reactive oxygen species homeostasis and expression of defense-related genes in tomato. Plant Cell Tiss. Organ. Cult. 123, 67–81. doi: 10.1007/s11240-015-0815-2
Liu, B., Hong, Y. B., Zhang, Y. F., Li, X. H., Huang, L., Zhang, H. J., et al. (2014). Tomato WRKY transcriptional factor SlDRW1 is required for disease resistance against Botrytis cinerea and tolerance to oxidative stress. Plant Sci. 227, 145–156. doi: 10.1016/j.plantsci.2014.08.001
Liu, G., Lai, D., Liu, Q. Z., Zhou, L., and Liu, Z. L. (2016). Identification of Nematicidal Constituents of Notopterygium incisum Rhizomes against Bursaphelenchus xylophilus and Meloidogyne incognita. Molecules 21:276. doi: 10.3390/molecules21101276
Liu, J.-H., Kitashiba, H., Wang, J., Ban, Y., and Moriguchi, T. (2007). Polyamines and their ability to provide environmental stress tolerance to plants. Plant Biotechnol. 24, 117–126. doi: 10.5511/plantbiotechnology.24.117
Madi, L., Wang, X., Kobiler, I., Lichter, A., and Prusky, D. (2003). Stress on avocado fruits regulates Δ9-stearoyl ACP desaturase expression, fatty acid composition, antifungal diene level and resistance to Colletotrichum gloeosporioides attack. Physiol. Mol. Plant Pathol. 62, 277–283. doi: 10.1016/S0885-5765(03)00076-6
Mamrot, J., Legaie, R., Ellery, S. J., Wilson, T., Seemann, T., Powell, D. R., et al. (2017). De novo transcriptome assembly for the spiny mouse (Acomys cahirinus). Sci. Rep. 7:8996. doi: 10.1038/s41598-017-09334-7
Manosalva, P., Manohar, M., von Reuss, S. H., Chen, S., Koch, A., Kaplan, F., et al. (2015). Conserved nematode signalling molecules elicit plant defenses and pathogen resistance. Nat. Commun. 6:7795. doi: 10.1038/ncomms8795
Marhavý, P., Kurenda, A., Siddique, S., Denervaud Tendon, V., Zhou, F., Holbein, J., et al. (2019). Single-cell damage elicits regional, nematode-restricting ethylene responses in roots. EMBO J. 38:972. doi: 10.15252/embj.2018100972
Meidani, C., Ntalli, N. G., Giannoutsou, E., and Adamakis, I. S. (2019). Cell Wall Modifications in Giant Cells Induced by the Plant Parasitic Nematode Meloidogyne incognita in Wild-Type (Col-0) and the fra2 Arabidopsis thaliana Katanin Mutant. Int. J. Mol. Sci. 20:21. doi: 10.3390/ijms20215465
Milligan, S. B., Bodeau, J., Yaghoobi, J., Kaloshian, I., Zabel, P., and Williamson, V. M. (1998). The root knot nematode resistance gene Mi from tomato is a member of the leucine zipper, nucleotide binding, leucine-rich repeat family of plant genes. Plant Cell 10, 1307–1319. doi: 10.1105/tpc.10.8.1307
Mitchum, M. G., Sukno, S., Wang, X., Shani, Z., Tsabary, G., Shoseyov, O., et al. (2004). The promoter of the Arabidopsis thaliana Cel1 endo-1,4-beta glucanase gene is differentially expressed in plant feeding cells induced by root-knot and cyst nematodes. Mol. Plant Pathol. 5, 175–181. doi: 10.1111/j.1364-3703.2004.00216.x
Mutuku, J. M., Cui, S., Hori, C., Takeda, Y., Tobimatsu, Y., Nakabayashi, R., et al. (2019). The Structural Integrity of Lignin Is Crucial for Resistance against Striga hermonthica Parasitism in Rice. Plant Physiol. 179, 1796–1809. doi: 10.1104/pp.18.01133
Nakasugi, K., Crowhurst, R., Bally, J., and Waterhouse, P. (2014). Combining transcriptome assemblies from multiple de novo assemblers in the allo-tetraploid plant Nicotiana benthamiana. PLoS One 9:e91776. doi: 10.1371/journal.pone.0091776
Ongena, M., Duby, F., Rossignol, F., Fauconnier, M. L., Dommes, J., and Thonart, P. (2004). Stimulation of the lipoxygenase pathway is associated with systemic resistance induced in bean by a nonpathogenic Pseudomonas. Mol. Plant Microbe Interact. 9, 1009–1018. doi: 10.1094/MPMI.2004.17.9.1009
Paal, J., Henselewski, H., Muth, J., Meksem, K., Menendez, C. M., Salamini, F., et al. (2004). Molecular cloning of the potato Gro1-4 gene conferring resistance to pathotype Ro1 of the root cyst nematode Globodera rostochiensis, based on a candidate gene approach. Plant J. 38, 285–297. doi: 10.1111/j.1365-313X.2004.02047.x
Palomares-Rius, J. E., Escobar, C., Cabrera, J., Vovlas, A., and Castillo, P. (2017). Anatomical alterations in plant tissues induced by plant-parasitic nematodes. Front. Plant Sci. 8:1987. doi: 10.3389/fpls.2017.01987
Patro, R., Duggal, G., Love, M. I., Irizarry, R. A., and Kingsford, C. (2017). Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419. doi: 10.1038/nmeth.4197
Postnikova, O. A., Hult, M., Shao, J., Skantar, A., and Nemchinov, L. G. (2015). Transcriptome analysis of resistant and susceptible alfalfa cultivars infected with root-knot nematode Meloidogyne incognita. PLoS One 10:e0118269. doi: 10.1371/journal.pone.0118269
Prior, A., Jones, J. T., Blok, V. C., Beauchamp, J., McDermott, L., Cooper, A., et al. (2001). A surface-associated retinol- and fatty acid-binding protein (Gp-FAR-1) from the potato cyst nematode Globodera pallida: lipid binding activities, structural analysis and expression pattern. Biochem. J. 356, 387–394. doi: 10.1042/0264-6021:3560387
Ralston, L., Kwon, S. T., Schoenbeck, M., Ralston, J., Schenk, D. J., Coates, R. M., et al. (2001). Cloning, heterologous expression, and functional characterization of 5-epi-aristolochene-1,3-dihydroxylase from tobacco (Nicotiana tabacum). Arch. Biochem. Biophys. 393, 222–235. doi: 10.1006/abbi.2001.2483
Ranjan, A., Ichihashi, Y., Farhi, M., Zumstein, K., Townsley, B., David-Schwartz, R., et al. (2014). De novo assembly and characterization of the transcriptome of the parasitic weed dodder identifies genes associated with plant parasitism. Plant Physiol. 166, 1186–1199. doi: 10.1104/pp.113.234864
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Sato, K., Kadota, Y., Gan, P., Bino, T., Uehara, T., Yamaguchi, K., et al. (2018). High-Quality Genome Sequence of the Root-Knot Nematode Meloidogyne arenaria Genotype A2-O. Genome Announc. 6:26. doi: 10.1128/genomeA.00519-18
Schulz, M. H., Zerbino, D. R., Vingron, M., and Birney, E. (2012). Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics 28, 1086–1092. doi: 10.1093/bioinformatics/bts094
Shukla, N., Yadav, R., Kaur, P., Rasmussen, S., Goel, S., Agarwal, M., et al. (2018). 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. Mol. Plant Pathol. 19, 615–633. doi: 10.1111/mpp.12547
Sijmons, P. C., Grundler, F. M. W., von Mende, N., Burrows, P. R., and Wyss, U. (1991). Arabidopsis thaliana as a new model host for plant-parasitic nematodes. Plant J. 1, 245–254. doi: 10.1111/j.1365-313X.1991.00245.x
Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., and Zdobnov, E. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212. doi: 10.1093/bioinformatics/btv351
Smant, G., Helder, J., and Goverse, A. (2018). Parallel adaptations and common host cell responses enabling feeding of obligate and facultative plant parasitic nematodes. Plant J. 93, 686–702. doi: 10.1111/tpj.13811
Song, N., Ma, L., Wang, W., Sun, H., Wang, L., Baldwin, I. T., et al. (2019). An ERF2-like transcription factor regulates production of the defense sesquiterpene capsidiol upon Alternaria alternata infection. J. Exp. Bot. 70, 5895–5908. doi: 10.1093/jxb/erz327
Sun, X. C., Gao, Y. F., Li, H. R., Yang, S. Z., and Liu, Y. S. (2015). Over-expression of SlWRKY39 leads to enhanced resistance to multiple stress factors in tomato. J. Plant Biol. 58, 52–60. doi: 10.1007/s12374-014-0407-4
Takahashi, K., Shimada, T., Kondo, M., Tamai, A., Mori, M., Nishimura, M., et al. (2010). Ectopic expression of an esterase, which is a candidate for the unidentified plant cutinase, causes cuticular defects in Arabidopsis thaliana. Plant Cell Physiol. 51, 123–131. doi: 10.1093/pcp/pcp173
Takahashi, S., Yeo, Y. S., Zhao, Y., O’Maille, P. E., Greenhagen, B. T., Noel, J. P., et al. (2007). Functional characterization of premnaspirodiene oxygenase, a cytochrome P450 catalyzing regio- and stereo-specific hydroxylations of diverse sesquiterpene substrates. J. Biol. Chem. 282, 31744–31754. doi: 10.1074/jbc.M703378200
Thomas, R., Fang, X., Ranathunge, K., Anderson, T. R., Peterson, C. A., and Bernards, M. A. (2007). Soybean root suberin: anatomical distribution, chemical composition, and relationship to partial resistance to Phytophthora sojae. Plant Physiol. 144, 299–311. doi: 10.1104/pp.106.091090
Uehara, T., Tateishi, Y., Kadota, Y., and Iwahori, H. (2017). Differences in parasitism of Meloidogyne incognita and two genotypes of M. arenaria on. J. Phytopathol. 94, 575–579. doi: 10.1111/jph.12594
van der Vossen, E. A., van der Voort, J. N., Kanyuka, K., Bendahmane, A., Sandbrink, H., Baulcombe, D. C., et al. (2000). Homologues of a single resistance-gene cluster in potato confer resistance to distinct pathogens: a virus and a nematode. Plant J. 23, 567–576.
Vanholme, R., Storme, V., Vanholme, B., Sundin, L., Christensen, J. H., Goeminne, G., et al. (2012). A systems biology view of responses to lignin biosynthesis perturbations in Arabidopsis. Plant Cell 24, 3506–3529. doi: 10.1105/tpc.112.102574
Veronico, P., Paciolla, C., Pomar, F., De Leonardis, S., Garcia-Ulloa, A., and Melillo, M. T. (2018). Changes in lignin biosynthesis and monomer composition in response to benzothiadiazole and root-knot nematode Meloidogyne incognita infection in tomato. J. Plant Physiol. 230, 40–50. doi: 10.1016/j.jplph.2018.07.013
Vishwanath, S. J., Delude, C., Domergue, F., and Rowland, O. (2015). Suberin: biosynthesis, regulation, and polymer assembly of a protective extracellular barrier. Plant Cell Rep. 34, 573–586. doi: 10.1007/s00299-014-1727-z
Vogeli, U., and Chappell, J. (1988). Induction of sesquiterpene cyclase and suppression of squalene synthetase activities in plant cell cultures treated with fungal elicitor. Plant Physiol. 88, 1291–1296. doi: 10.1104/pp.88.4.1291
Wang, J., Kucukoglu, M., Zhang, L., Chen, P., Decker, D., Nilsson, O., et al. (2013). The Arabidopsis LRR-RLK, PXC1, is a regulator of secondary wall formation correlated with the TDIF-PXY/TDR-WOX4 signaling pathway. BMC Plant Biol. 13:94. doi: 10.1186/1471-2229-13-94
Wang, X., and Culver, J. N. (2012). DNA binding specificity of ATAF2, a NAC domain transcription factor targeted for degradation by Tobacco mosaic virus. BMC Plant Biol. 12:157. doi: 10.1186/1471-2229-12-157
Ward, E. W. B., Unwin, C. H., and Stoessl, A. (1974). Postinfectional inhibitors from plants. XIII. Fungitoxicity of the phytoalexin, capsidiol, and related sesquiterpenes. Can. J. Bot. 52, 2481–2488. doi: 10.1139/b74-324
Xie, Y., Wu, G., Tang, J., Luo, R., Patterson, J., Liu, S., et al. (2014). SOAPdenovo-Trans: de novo transcriptome assembly with short RNA-Seq reads. Bioinformatics 30, 1660–1666. doi: 10.1093/bioinformatics/btu077
Xing, X., Li, X., Zhang, M., Wang, Y., Liu, B., Xi, Q., et al. (2017). Transcriptome analysis of resistant and susceptible tobacco (Nicotiana tabacum) in response to root-knot nematode Meloidogyne incognita infection. Biochem. Biophys. Res. Commun. 482, 1114–1121. doi: 10.1016/j.bbrc.2016.11.167
Yamaguchi, H., Fukuoka, H., Arao, T., Ohyama, A., Nunome, T., Miyatake, K., et al. (2010). Gene expression analysis in cadmium-stressed roots of a low cadmium-accumulating solanaceous plant. J. Exp. Bot. 131, 423–437. doi: 10.1093/jxb/erp313
Yan, L., Zhai, Q., Wei, J., Li, S., Wang, B., and Huang, T. (2013). Role of tomato lipoxygenase D in wound-induced jasmonate biosynthesis and plant immunity to insect herbivores. PLoS Genet 9:e1003964. doi: 10.1371/journal.pgen.1003964
Yang, X., Cheng, Y. F., Deng, C., Ma, Y., Wang, Z. W., Chen, X. H., et al. (2014). Comparative transcriptome analysis of eggplant (Solanum melongena L.) and turkey berry (Solanum torvum). BMC Genomics 15:412. doi: 10.1186/1471-2164-15-412
Ye, D. Y., Qi, Y. H., Cao, S. F., Wei, B. Q., and Zhang, H. S. (2017). Histopathology combined with transcriptome analyses reveals the mechanism of resistance to Meloidogyne incognita in Cucumis metuliferus. J. Plant Physiol. 212, 115–124. doi: 10.1016/j.jplph.2017.02.002
Zacheo, G., Bleve-Zacheo, T., and Melillo, M. T. (1997). “Biochemistry of Plant Defence Responses to Nematode Infection,” in Cellular and Molecular Aspects of Plant-Nematode Interactions. Developments in Plant Pathology, Vol. 10, eds C. Fenoll, F. M. W. Grundler, and S. A. Ohl, (Dordrecht: Springer).
Zasada, I. A., Halbrendt, J. M., Kokalis-Burelle, N., LaMondia, J., McKenry, M. V., and Noling, J. W. (2010). Managing nematodes without methyl bromide. Annu. Rev. Phytopathol. 48, 311–328. doi: 10.1146/annurev-phyto-073009-114425
Zhang, H., Kjemtrup-Lovelace, S., Li, C., Luo, Y., Chen, L. P., and Song, B. H. (2017). Comparative RNA-Seq Analysis Uncovers a Complex Regulatory Network for Soybean Cyst Nematode Resistance in Wild Soybean (Glycine soja). Sci. Rep. 7:9699. doi: 10.1038/s41598-017-09945-0
Zhang, L., Lilley, C. J., Imren, M., Knox, J. P., and Urwin, P. E. (2017). The Complex Cell Wall Composition of Syncytia Induced by Plant Parasitic Cyst Nematodes Reflects Both Function and Host Plant. Front. Plant Sci. 8:1087. doi: 10.3389/fpls.2017.01087
Zhang, P. J., He, Y. C., Zhao, C., Ye, Z. H., and Yu, X. P. (2018). Jasmonic Acid-Dependent Defenses Play a Key Role in Defending Tomato Against Bemisia tabaci Nymphs, but Not Adults. Front. Plant Sci. 9:1065. doi: 10.3389/fpls.2018.01065
Zhou, J., Wang, J., Zheng, Z., Fan, B., Yu, J. Q., and Chen, Z. (2015). Characterization of the promoter and extended C-terminal domain of Arabidopsis WRKY33 and functional analysis of tomato WRKY33 homologues in plant stress responses. J. Exp. Bot. 66, 4567–4583. doi: 10.1093/jxb/erv221
Keywords: Turkey berry, plant immunity, plant-parasitic nematode, lignin deposition, comparative transcriptomics
Citation: Sato K, Uehara T, Holbein J, Sasaki-Sekimoto Y, Gan P, Bino T, Yamaguchi K, Ichihashi Y, Maki N, Shigenobu S, Ohta H, Franke RB, Siddique S, Grundler FMW, Suzuki T, Kadota Y and Shirasu K (2021) Transcriptomic Analysis of Resistant and Susceptible Responses in a New Model Root-Knot Nematode Infection System Using Solanum torvum and Meloidogyne arenaria. Front. Plant Sci. 12:680151. doi: 10.3389/fpls.2021.680151
Received: 13 March 2021; Accepted: 30 April 2021;
Published: 25 May 2021.
Edited by:Daguang Cai, University of Kiel, Germany
Reviewed by:Mihir Kumar Mandal, Claflin University, United States
Yunliang Peng, Sichuan Academy of Agricultural Sciences, China
Copyright © 2021 Sato, Uehara, Holbein, Sasaki-Sekimoto, Gan, Bino, Yamaguchi, Ichihashi, Maki, Shigenobu, Ohta, Franke, Siddique, Grundler, Suzuki, Kadota and Shirasu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.