Next Article in Journal
Integrated Computational Approaches and Tools for Allosteric Drug Discovery
Next Article in Special Issue
An Overview of the Molecular Genetics of Plant Resistance to the Verticillium Wilt Pathogen Verticillium dahliae
Previous Article in Journal
Editorial of Special Issue “Adipokines 2.0”
Previous Article in Special Issue
Modulation of (Homo)Glutathione Metabolism and H2O2 Accumulation during Soybean Cyst Nematode Infections in Susceptible and Resistant Soybean Cultivars
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis of Rice Roots in Response to Root-Knot Nematode Infection

1
College of Plant Protection, Nematology Institute of Northern China, Shenyang Agricultural University, Shenyang 110866, China
2
Shaanxi Key Laboratory of Chinese Jujube, College of Life Science, Yan’an University, Yan’an, Shaanxi 716000, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2020, 21(3), 848; https://doi.org/10.3390/ijms21030848
Submission received: 6 December 2019 / Revised: 24 January 2020 / Accepted: 27 January 2020 / Published: 28 January 2020
(This article belongs to the Special Issue Plant Disease Resistance)

Abstract

:
Meloidogyne incognita and Meloidogyne graminicola are root-knot nematodes (RKNs) infecting rice (Oryza sativa L.) roots and severely decreasing yield, whose mechanisms of action remain unclear. We investigated RKN invasion and development in rice roots through RNA-seq transcriptome analysis. The results showed that 952 and 647 genes were differently expressed after 6 (invasion stage) and 18 (development stage) days post inoculation, respectively. Gene annotation showed that the differentially expressed genes were classified into diverse metabolic and stress response categories. Furthermore, phytohormone, transcription factor, redox signaling, and defense response pathways were enriched upon RKN infection. RNA-seq validation using qRT-PCR confirmed that CBL-interacting protein kinase (CIPK) genes (CIPK5, 8, 9, 11, 14, 23, 24, and 31) as well as brassinosteroid (BR)-related genes (OsBAK1, OsBRI1, D2, and D11) were altered by RKN infection. Analysis of the CIPK9 mutant and overexpressor indicated that the RKN populations were smaller in cipk9 and larger in CIPK9 OX, while more galls were produced in CIPK9 OX plant roots than the in wild-type roots. Significantly fewer numbers of second-stage infective juveniles (J2s) were observed in the plants expressing the BR biosynthesis gene D2 mutant and the BR receptor BRI1 activation-tagged mutant (bri1-D), and fewer galls were observed in bri1-D roots than in wild-type roots. The roots of plants expressing the regulator of ethylene signaling ERS1 (ethylene response sensor 1) mutant contained higher numbers of J2s and developed more galls compared with wild-type roots, suggesting that these signals function in RKN invasion or development. Our findings broaden our understanding of rice responses to RKN invasion and provide useful information for further research on RKN defense mechanisms.

1. Introduction

Rice (Oryza sativa) is one of the most important cereal crops. It is a major food source globally and thus has potential for improving global food security [1]. However, rice is vulnerable to injury from a range of pathogens and pests, such as those causing rice blast disease, viruses, insect pests, smut-causing pathogens, and root-knot nematodes (RKNs, Meloidogyne spp.). Disease occurrence severely reduces the production yield of rice.
RKNs include several species of Meloidogyne (e.g., Meloidogyne graminicola and Meloidogyne incognita) and are important soil pathogens that significantly influence the yield and quality of rice [2]. Among nematodes, RKNs and cyst nematodes (Heterodera spp. and Globodera spp.) are the most damaging to plants [3,4,5]. Economic losses of over $100 billion per year have been reported globally due to nematode infection [6]. Second-stage infective juveniles (J2) pierce root cell walls with their stylet at the elongation zone, migrate toward the tip, and then to the plant vascular cylinder to ingest cytoplasmic contents [7,8,9]. Plant cells respond to RKN infection by forming galls or root-knots, which are feeding sites for the parasites.
After feeding, nematodes progress through two developmental stages (J3 and J4) before becoming adults. Adult RKN females are pyriform and remain in the root to begin the next round of production, while males become vermiform and leave the root [10,11]. RKNs secrete effector proteins, which influence host functions, such as cell cycle, cell wall properties, metabolism, and defense, and hijack host cell physiology [12].
Plants have evolved chemical and physical defense mechanisms to protect against pathogen invasion and pathogen infection. Plant hormones such as salicylic acid, jasmonate, gibberellin, abscisic acid, brassinosteroids (BR), cytokinins, and auxin are important for basal immunity in rice [13,14,15,16,17,18]. Several previous studies have reported that plant–nematode interactions affect plant hormones and regulate response mechanisms [19,20]. Recently, sources of genetic resistance to RKNs have been found in the form of nematode resistance genes, which were isolated from different crop plants [21,22]. In rice, Has-1Og was reported to confer resistance to the cyst nematode Heterodera sacchari [23].
Transcriptome analysis is a useful tool to determine transcriptional changes in cells under different conditions [24]. Previous research has provided evidence to suggest that plant–nematode interactions affect the expression of genes associated with plant immune response in dicotyledons [25,26]. Nguyễn et al. (2014) demonstrated that some defense genes in rice were downregulated in the early stages of infection. This suggests that rice is a good model system for studying monocotyledon–nematode molecular interactions [27]. Moreover, increased expression of plant immunity genes in RKN-infected rice, which affects susceptibility to blast disease, has been previously tested using high-throughput sequencing [28]. The results demonstrated that RKN infection also affected oxidative stress response and hormone balance [29]. However, RKN resistance to molecular, physiological, and genetic mechanisms in rice remain poorly understood.
In the present study, the transcriptome profile of RKN-infected rice roots was analyzed, and the differentially expressed genes were annotated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Further genetic study using rice mutants suggested that our approach might be an efficient way to isolate genes associated with resistance to RKN. Overall, our findings provide a theoretical basis for studying the mechanisms underlying RKN defense in rice.

2. Results

2.1. Analysis of RKN Invasion and Development in Rice Roots

Most of the J2 penetrated into roots a few days post inoculation (3 dpi) and then migrated inside the root toward the feeding site. The J2s became sedentary (6 dpi) when they reached the feeding site after passing through the vascular cylinder of the roots. At 9 to 24 dpi, the cells of the feeding sites began to form root-knot galls around expanding giant cells, and the juveniles grew thicker, changing from a sausage-like shape to a globose shape by molting twice. At 27 dpi, the males emerged from the root, while the females continued to thicken and appeared pyriform (Figure 1).

2.2. Transcriptome Analysis to Identify Differentially Expressed Genes upon RKN Infection

Since the J2s were sedentary at 6 dpi and the galls formed between 9 and 24 dpi, RNA-sequencing data from root tissues after 6 and 18 dpi were selected to analyze transcriptome profiles and explore the mechanisms involved in invasion and gall formation. We found 952 and 647 differentially expressed genes (DEGs) after 6 and 18 dpi, respectively (Tables S1 and S2). Among them, 492 and 460 genes were up- and downregulated at 6 dpi, respectively, while 499 and 148 genes were up- and downregulated at 18 dpi, respectively, by RKN infection (Figure 2A,B). In addition, 15 genes were commonly upregulated, and 477 and 484 genes were uniquely upregulated at 6 and 18 dpi, respectively (Figure 2C). Similarly, 443 and 131 downregulated genes were unique to 6 and 18 dpi, respectively (Figure 2D). The 15 commonly upregulated genes included genes encoding actin-depolymerizing factor, myb-like DNA-binding domain-containing protein, lipase, and OsHKT2;1 sodium transporter, while the 17 commonly downregulated genes included genes encoding rhodanese-like domain-containing protein, dihydroflavonol-4-reductase, C3HC4-type domain-containing protein, and RING-H2 finger protein (Table S3).
Furthermore, the DEGs were classified according to GO and KEGG terms. 952 DEGs at 6 dpi were enriched in 20 GO terms, including response to wounding, MAPKKK activity, chitin activity, and chitin catabolic process (Figure 3A). Also, 647 DEGs at 18 dpi were classified by 20 GO terms, including symporter activity, response to wounding, response to chitin, monooxygenase activity, and lignin catabolic activity (Figure 3B). KEGG analysis showed that DEGs at 6 dpi displayed enrichment in diverse metabolic pathways including phenylpropanoid metabolism, phenylalanine metabolism, drug metabolism, and methane metabolism (Figure 3C). In addition, the DEGs at 18 dpi were classified as linked to phenylpropanoid metabolism, phenylalanine metabolism, starch and sucrose metabolism, and methane metabolism (Figure 3D). MapMan analysis was also performed to classify DEGs. This led to the identification of phytohormone- (auxin, BR, abscisic acid, ethylene, salicylic acid, and jasmonate) related, redox signaling, transcription factors (ERF, bZIP, WRKY, MYB, and DOF), and pathogen-related genes (Figure 4).

2.3. CIPK Genes in Rice Response to RKN

RNA-seq results showed that the expression levels of CIPK14(LOC_Os12g02200) and CIPK23(LOC_Os07g05620) were altered at 18 dpi (p < 0.05, Table S2). In addition, few CIPK genes (CIPK5(LOC_Os01g10890), CIPK8(LOC_Os01g35184), CIPK9(LOC_Os03g03510), CIPK11(LOC_Os01g60910), CIPK24(LOC_Os06g40370), and CIPK31(LOC_Os03g20380)) were listed in the RNA-seq data with altered expressions after RKN-infection, but these differences were not significant (p > 0.05, data not shown). To validate the RNA-seq data, qRT-PCR analysis was performed. CIPK5, 9, 14, 24, and 31 were significantly induced by RKN infection at 6 dpi compared with the control group, with CIPK9 exhibiting the highest induction (about 12-fold). Also, CIPK8, 14, and 31 were induced at 18 dpi by RKN infection, with CIPK8 exhibiting the highest induction (about 35-fold). In addition, CIPK11 and 23 were suppressed at 6 dpi, and CIPK9 and 11 were suppressed at 18 dpi compared with the control group (Figure 5A–H).
CIPK genes were differentially regulated by RKN infection. CIPK9 was the most highly induced gene at 6 dpi. Therefore, a cipk9 mutant, a CIPK9 overexpressor (CIPK9 OX), and a wild-type control were used to analyze their responses to RKN infection. After 6 dpi, a fewer J2 populations were observed in cipk9 mutants than in the wild-type control. However, no significant difference was observed between the wild-type and the CIPK9 OX roots (Figure 5I).
The number of galls was determined after 18 days of J2 inoculation, and the results showed that CIPK9 OX roots developed more galls than the wild-type and cipk9 mutant roots. The cipk9 mutant showed a reduction of galls of almost 32% compared with the wild-type control (Figure 5J). To further evaluate whether CIPK9 regulates RKN development in roots, the number of J2, sausage-shaped J2, and globose RKN in the roots of each plant were determined.
The results showed that the total number of nematodes was lower in the roots of cipk9 mutants and higher in those of CIPK9 OX plants. The cipk9 mutant roots contained lower numbers of J2 and sausage-shaped J2, while the CIPK9 OX roots contained higher numbers of J2 compared with the wild-type control. CIPK9 OX contained a similar number of sausage-shaped J2 to that of the wild-type plants, and the number of globose nematodes was similar in the roots of wild-type, cipk9, and CIPK9 OX plants (Figure 5K). Furthermore, the results showed that cipk9 mutants contained a lower percentage of J2, while CIPK9 OX contained a higher percentage of J2 than wild-type roots. Sausage-shaped J2 and globose RKNs were present in higher proportion in cipk9 mutants and in lower proportion in CIPK9 OX plants compared with the wild-type control (Figure 5L).

2.4. BR Signaling Regulation of RKN Invasion and Development

Despite CIPK genes, a BR co-receptor gene (OsBAK1(LOC_Os06g16330)) was induced by RKN infection at 18 dpi. To further examine BR signaling and biosynthesis gene response to RKN infection, the expression levels of two biosynthesis genes (D2(LOC_Os01g10040) and D11(LOC_Os04g39430)), the receptor gene OsBRI1 (LOC_Os01g52050), and the gene OsBAK1 were analyzed by qRT-PCR. The results indicated that D2, D11, and OsBRI1 were induced at 6 dpi, while OsBAK1 was induced at 18 dpi compared with wild-type roots (Figure 6A–D). Furthermore, the d2 mutant, the OsBRI1 activation tagged mutant bri1-D (same as the OsBRI1 overexpressor) [30], and wild-type control plants (Nipponbare) were used to analyze the effect of BR signaling on RKN invasion and development. After 6 dpi, fewer J2 populations were detected in d2 and bri1-D mutants compared with the wild-type plants (Figure 6E). At 18 dpi, bri1-D mutants developed significantly fewer galls compared with the d2 mutant and wild-type plants (Figure 6F).
To further analyze RKN development, the numbers of J2, sausage-shaped J2, and globose RKN were determined for wild-type, d2, and bri1-D mutant roots at 18 dpi. The results indicated that the total number of nematodes per gram of roots was lower in bri1-D and d2 mutants than in the wild-type roots. In detail, bri1-D and d2 roots contained more J2 populations but fewer sausage-shaped J2 populations than the wild-type roots. There were more globose RKNs in bri1-D roots and fewer globose RKNs in d2 mutant roots than in the wild-type ones (Figure 6G). Furthermore, the results showed that d2 and bri1-D mutant roots contained significantly more J2 than the wild-type roots. In addition, evidently fewer sausage-shaped J2 were found in d2 and bri1-D roots than in the wild-type ones. However, interestingly, there were more globose RKNs in the bir1-D mutant than in the wild-type and d2 mutant roots. There were fewer globose RKNs in the d2 mutant than in the wild-type roots (Figure 6H).

2.5. Activation of Ethylene Signaling during RKN Invasion and Development

MapMan analysis revealed that ethylene-related genes and ERFs were enriched. To analyze the role of ethylene signaling in RKN invasion and development, a mutant of ERS1, a negative regulator of ethylene signaling was used [31]. At 6 dpi, J2 populations were significantly more numerous in ers1 than in wild-type roots (Figure 7A); in addition, ers1 roots developed more galls after 18 days post J2 inoculation (Figure 7B).
Further examination revealed that the total number of nematodes was higher in the ers1 mutant than in the wild-type control. Among them, more numerous J2 populations and fewer globose RKN were found in ers1 compared with wild-type roots. In addition, ers1 and wild-type roots contained a similar number of sausage-shaped J2 (Figure 7C). Further analysis showed that the proportion of J2 was significantly higher and that of sausage-shaped J2 and globose RKN was significantly lower in ers1 roots compared with wild-type roots at 18 dpi (Figure 7D).

3. Discussion

The M. incognita–rice pathosystem can serve as a model to elucidate the interactions between RKNs and the host [27]. To identify different genes expressed in rice roots in both the invasion and the development phases of M. incognita, the root tissue was sampled at different time points and visualized by staining with acid fuchsin. Our histological examination demonstrated that most J2-stage RKNs invaded before 3 dpi, and feeding sites were successfully constructed at about 6 dpi. Furthermore, continual development of juveniles feeding on giant cells from 9 to 24 dpi, mature females beginning egg production, and males returning to the rhizosphere from 27 dpi were observed. In the RNA-seq-based transcriptome analysis, 952 and 647 DEGs were identified at 6 and 18 dpi in rice roots. At these two stages, similar numbers of induced genes were observed, but the number of downregulated genes was evidently lower at 18 dpi than at 6 dpi. However, only 15 and 17 commonly up- and downregulated genes, respectively, were found at these two stages, suggesting that the gene response differed in these stages.
DEGs from the two stages were commonly classified as linked to response to wound, lignin catabolic process, laccase activity, and response to chitin. Chitin is an essential component of the nematode egg shells and pharynx, and disturbing chitin synthesis or hydrolysis would lead to the failure of nematode embryos or molting [32]. This suggests that activating chitin catabolism by the induction of chitinase genes, plays a key role in host plant response to RKN invasion and gall formation. Nitrate transport and symporter activity were enriched at 6 and 18 dpi, respectively, implying that nutrient transport might be important in these stages.
At 6 dpi, the MAPK signaling term was identified. Previous research demonstrated that the MAPK gene family functions as a cohort during soybean (Glycine max) response to Heterodera glycines [33]. This suggests a key role for MAPK signaling in plant response to different species of plant parasitic nematodes. KEGG pathway analysis showed that C and N metabolism were enriched at 6 and 18 dpi, suggesting that the host plant actively fine-tunes nutrient metabolism in response to RKNs.
DEG annotation using MapMan suggested enrichment in genes related to phytohormones (auxin, BR, abscisic acid, ethylene, salicylic acid, and jasmonate), redox signaling, transcription factors (ERF, bZIP, WRKY, MYB, and DOF), and pathogens. Previous research demonstrated that in galls at 3 dpi infected with the RKN M. graminicola, genes for gibberellin and BR biosynthesis and signaling were activated, while those for salicylic acid and ethylene signaling were repressed [34]. Exogenous application of salicylic acid, jasmonate, or ethylene increased rice response to RKNs, with the strongest effect observed under jasmonate treatment [35]. In addition, ethylene was shown to promote plant defense against RKN by activating jasmonate biosynthesis [35]. Interestingly, the roots of BRI1 mutants infected with RKNs also showed a 30% reduction in gall numbers [17]. BR biosynthesis or BR receptor silencing increased the susceptibility of tomato plants to RKN, independent of salicylic acid and jasmonic acid signaling [36]. Furthermore, our recent work showed that BR signaling negatively regulated Arabidopsis resistance to RKN [37]. This result suggests that the response to RKNs is related to BR signaling.
To analyze phytohormone and RKN interactions, genes involved in BR signaling and BR biosynthesis were examined. The data indicated that BR signaling activation (bri1-D) or biosynthesis deficiency (d2) inhibited J2 M. incognita invasion of rice roots, suggesting the presence of diverse BR functions responding to different species of RKN. In addition, bri1-D, in which BR signaling was activated, developed fewer galls in the roots. There were fewer J2 and sausage-shaped J2 RKNs in d2 and bri1-D mutants than in the wild-type roots, while more and fewer globose RKNs were observed in bri1-D and d2 roots, respectively, than in the wild-type ones, suggesting that BR signaling also regulates RKN development in rice roots. In parallel, the expression of D2, D11, OsBRI1, and OsBAK1 was induced during RKN infection, implying a close connection between BR action and RKN invasion and development in rice roots. Ethylene signaling is known to positively regulate plant response to RKN invasion. In our analysis, ers1 mutants, in which ethylene signaling was activated, significantly promoted the invasion of J2 RKNs and the formation of galls, suggesting that activation of ethylene signaling inhibits rice root defense against M. incognita.
Calcium signaling plays a key role in plant response to environmental change. Previously, calcium/calmodulin-mediated defense signaling was reported to be a key mechanism underlying soybean cyst nematode resistance in wild Glycine soja [38]. Furthermore, calcium sensor Cbl10 together with its interacting partner CIPK6 regulate plant immunity in tomato plants [39]. In the RNA-seq data, we observed a significantly altered expression level of CIPK genes in rice roots. Further genetic study using CIPK9 mutants and the overexpressor revealed that CIPK9 negatively regulates rice defense to RKN invasion, as indicated by the fewer RKNs observed in the cipk9 mutant at 6 and 18 dpi. The CIPK9 mutation promoted RKN development, as indicated by the higher percentage of sausage-shaped J2 populations, while CIPK9 overexpression increased the total nematode population in roots and inhibited RKN development. CIPK9 overexpression in plant roots led to the development of more galls than in wild-type plants.
In the CIPK9 OX, bri1-D, and ers1 roots, the number of galls appears to be associated with the total RKN population, as shown by the observation that CIPK9 OX and ers1 contained a higher number of total RKNs and galls than bri1-D and the wild-type roots. However, the cipk9 mutant developed a similar number of galls as the wild-type control, even though cipk9 mutants contained significantly less numerous total RKN populations. This suggests that gall formation was regulated by diverse factors, and the details of the underlying mechanism require further analysis. Our findings will be important for further research on the defense mechanisms of rice against M. incognita invasion. Furthermore, the hormone and calcium signaling regulators identified herein might be useful for developing nematode-resistant breeding strategies.

4. Materials and Methods

4.1. Nematode Culture

M. incognita were maintained on a nematode-susceptible tomato cultivar (L402, Liaoning Horticultural Seedling Co. Ltd.). Eggs were collected as described previously [40], with modifications. Briefly, roots of glasshouse-grown tomatoes were cut into pieces and then shook vigorously for 5 min with 10% commercial bleach. Pieces of root were poured through a 180 μm mesh, and eggs were collected on a 25 μm mesh screen. Eggs were quickly purified by centrifugation in 35% sucrose at 720 g for 10 min in a 100 mL centrifuge tube. The supernatant containing the eggs was subjected to another round of centrifugation for 5 min at 720 g in 10% bleach and then rinsed five times in sterile water. Eggs were hatched at 25 °C, and juveniles were collected after passing through 10 layers of Kimwipes tissues into sterile water.

4.2. Plant Growth and Nematode Inoculation Assay

O. sativa “Nipponbare”, “Dong jin”, cipk9: Ds and CIPK9 overexpressor revealed that CIPK9 regulates ammonium-dependent root growth [41]; bri1-D (BRI1 activation tagging line) in which BR receptor gene BRI1 was over expressed [42], d2 rice dwarf mutant with mutation of the gene for the BR biosynthesis enzyme [43], and ers1 mutant exhibiteing hyper-sensitive response to ethylene treatment [44] were used. The seeds were germinated for 6 days at 30 °C and transferred to plastic Ray Leach containers (SC7 Stubby, diameter, 3.81 cm; depth, 13.97 cm; Stuewe & Sons Inc. Corvallis, OR) that had been filled with sterile sand and further grown at 25 °C under a 16h light/8h dark photoperiod. The plants were inoculated with 1.5 mL of water containing 2000 freshly hatched J2s of M. incognita.
For RKN penetration and growth assays, 60 seeds were planted as described above. Ten-day-old seedlings were inoculated with 2000 J2s per seedling and maintained as described previously. Seedlings were gently removed from the container for staining at 3, 6, 9, 12, 15, 18, 21, 24, and 27 days post inoculation. Roots of seedlings were treated with 10% bleach for 1 min, flushed with water, and boiled for 1 min in acid fuchsin solution (3.5% acid fuchsin in 25% acetic acid). Then, the seedlings were transferred to a de-staining acidified glycerin solution (200 μL of 5 mol·L−1 HCl in 50 mL glycerol) that was heated to about 40 °C (not boiled) [45,46]. The number of nematodes inside the roots was counted using a stereoscope (Nikon SMZ800, Nikon, Tokyo, Japan).

4.3. RNA Extraction, Library Preparation, and Sequencing

For transcriptome analysis, the roots of Nipponbare rice were collected at days 6 and 18 post RKN inoculation (dpi). For each time point, three independent biological replicates and a control (plant in the same stage but without RKN inoculation) were used for RNA sequencing. Total RNA was extracted using Trizol reagent (Dingguo Biotechnology Co., Ltd., Beijing, China), following the manufacturer’s protocols. Total RNA quantity and purity were analyzed with the Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, Santa Clara, CA, USA) with RIN number > 7.0, and RNA sequencing was performed at LC Biotech Co., Hangzhou, China.
Approximately 10 µg of total RNA was used to isolate poly-(A) mRNA with poly-T oligo-attached magnetic beads (Invitrogen). Following purification, the mRNA was fragmented into small pieces using divalent cations under elevated temperature. Then. the cleaved RNA fragments were reverse-transcribed to create the final cDNA library in accordance with the protocol for the mRNA-Seq sample preparation kit (Illumina, San Diego, CA, USA). The average insert size for the paired-end libraries was 300 bp (±50 bp). Then, we performed paired-end sequencing on an Illumina HiSeq 2000/2500 (LC Sciences, Houston, TX, USA), following the manufacturer’s protocol. To validate the quality of the raw sequence data, forward and reverse reads were analyzed using fastQC (version 0.10.1). Sequencing saturation and 5′ and 3′ biases of RNA-seq data were assessed using RSeQC (version 2.4). Basic data processing consisted of splicing-aware alignment using Tophat (version 2.0.4). RNA-seq reads were aligned to the rice reference genome (https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Osativa) by Bowtie (version 2.0.0, Baltimore, MD, USA). Sequences were assembled using StringTie (version 1.3.0, Baltimore, MD, USA), with Ballgown used for differential expression. The original sequencing volume, effective sequencing volume, Q20, Q30, and GC content were determined (Table S4). The raw data were uploaded to the NCBI Sequence Read Arhcive (SAR) database. The BioProject ID is PRJNA589993, and the SRA ID is SRP230566.

4.4. Differential Gene Expression Analysis

The DEGs between the two groups of samples were identified using Cuffdiff [47]. Only the genes with a log2 fold change ≥1 or ≤−1 and a p-value ≤ 0.05 were considered to be significant DEGs. DEGs associated with biological process GO terms [48] were determined as follows: Fisher’s exact test was used to classify the enrichment of the GO category, and the false discovery rate was calculated to correct the p values. The smaller the false discovery rate, the smaller the error in judging p [49]. Enrichment of GO terms among probe sets was found using one-tailed Fisher’s exact test based on 2 × 2 contingency tables, because it provides a measure of the significance of the enrichment as it increases, and Pathway analysis was performed against the Kyoto Encyclopedia of Genes and Genomes. For both 6 dpi and 18 dpi, the differentially modulated transcripts were compared and mapped by MapMan software (Version 3.6.0, Cologne, North Rhine-Westphalia, Germany), providing a metabolic overview of transcriptional changes.

4.5. qRT-PCR Analysis

For quantitative real-time PCR (qRT-PCR) analysis, total RNA was extracted from samples using Trizol reagent (Dingguo Biotechnology Co., Ltd., Beijing, China) with three biological replications. The cDNA was synthesized from 1 µg of total RNA using the PrimeScript RT reagent kit (Takara Bio, Tokyo, Japan), according to the manufacturer’s instructions. The primers used for these experiments are listed in Table S5. PCR reactions were set up in 96-well hard-shell PCR plates with 0.4 μM primers, using One Step SYBR PrimeScript RT-PCR kit (Takara Bio) in a 10 μL solution. Reaction conditions were as follows: denaturation at 95 °C for 3 min, 40 cycles at 95 °C for 10 s and at 58 °C for 30 s, heating from 60 °C to 95 °C at a rate of 1 °C per 4 s for melt curve analysis. Finally, relative gene expression levels were normalized to the levels of housekeeping gene ubiquitin and calculated using the 2−△△Ct method [50].

4.6. Statistical Analysis

All data were analyzed by ANOVA using SPSS version 22 (IBM, Armonk, NY, USA) for statistical significance. Index and density data complied with normality according to the Kolmogorov–Smirnov test and variance homogeneity as determined with the Kruskal–Wallis test. All data were analyzed using paired t-test and Duncan’s multiple comparison test.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/3/848/s1.

Author Contributions

X.Z., Y.X. and Y.Z. conceived and designed the experiments. X.Z. and D.Z. obtained the funding. Y.Z. carried out the experiments. X.Z., Y.Z., Y.X., L.S., and D.Z. analyzed the experiments. Y.X., X.Z., Y.D., D.X., and L.C. contributed to the materials and analysis tools. X.Z., Y.Z., and Y.X. wrote an initial version of the manuscript that was subsequently critically revised by Y.W., X.L., and H.F. All authors read and approved the final version of the paper.

Funding

This work was supported by a grant from the National Natural Science Foundation of China (31571985), National Sharing Service Platform for Parasite Resource (No. TDRC - 30), China Agriculture Research System (No. CARS-04-PS13), Educational Commission of Liaoning Province, China, and the Isolation and Identification of Nematicidal Compound and Analysis of Mechanism of the Fungi Snef5 (LSNQN201723).

Conflicts 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.

Abbreviations

RKNsroot-knot nematode
BRbrassinosteroid
J2second-stage infective juveniles
GOGene Ontology
KEGGKyoto Encyclopedia of Genes and Genomes
dpidays post inoculation
DEGsdifferentially expressed genes
qRT-PCRquantitative real-time PCR

References

  1. FAOSTAT Food and Agriculture Organization of the United Nations. FAOSTAT Database. 2002. Available online: http://www.fao.org/faostat/en/#data/QC (accessed on 30 May 2019).
  2. Castagnone-Sereno, P.; Danchin, E.G.; Perfus-Barbeoch, L.; Abad, P. Diversity and evolution of root-knot nematodes, genus Meloidogyne: New insights from the genomic era. Annu. Rev. Phytopathol. 2013, 51, 203–220. [Google Scholar] [CrossRef]
  3. Bridge, J.; Starr, J. Plant Nematodes of Agricultural Importance Acolor Handbook; Academic Press: San Diego, CA, USA, 2007; pp. 77–78. [Google Scholar] [CrossRef] [Green Version]
  4. Dutta, T.K.; Ganguly, A.K.; Gaur, H.S. Global status of rice root-knot nematode, Meloidogyne graminicola. Afr. J. Microbiol. Res. 2012, 6, 6016–6021. [Google Scholar]
  5. Williamson, V.M.; Gleason, C.A. Plant-nematode interactions. Curr. Opin. Plant Biol. 2003, 6, 327–333. [Google Scholar]
  6. Chitwood, D.J. Research on plant-parasitic nematode biology conducted by the united states department of agriculture-agricultural research service. Pest. Manag. Sci. 2003, 59, 748–753. [Google Scholar] [CrossRef]
  7. Bird, D.M.; Opperman, C.H.; Williamson, V.M. Plant Infection by Root-Knot Nematode. Cell Biology of Plant Nematode Parasitism. Plant Cell Monogr. 1970, 15, 1–13. [Google Scholar] [CrossRef]
  8. Berg, R.H.; Taylor, C. Plant Infection by Root-Knot Nematode of Cell Biology of Plant Nematode Parasitism; Springer: Berlin/Heidelberg, Germany, 2009; p. 113. [Google Scholar] [CrossRef]
  9. Jones, J.; Gheysen, G.; Fenoll, C. Genomics and Molecular Genetics of Plant–Nematode Interactions; Springer: Dordrecht, The Netherlands, 2011. [Google Scholar] [CrossRef]
  10. Bridge, J.P.R.; Peng, D. Nematode parasites of rice. In Plant Parasitic Nematodes in Subtropical and Tropical Agriculture; Luc, M., Sikora, R.A., Bridge, J., Eds.; CABI: Wallingford, UK, 2005. [Google Scholar] [CrossRef]
  11. Kyndt, T.; Fernandez, D.; Gheysen, G. Plant-parasitic nematode infections in rice: molecular and cellular insights. Annu. Rev. Phytopathol. 2014, 52, 135–153. [Google Scholar] [CrossRef]
  12. Hewezi, T.; Baum, T.J. Manipulation of plant cells by cyst and root-knot nematode effectors. MPMI 2013, 26, 9–16. [Google Scholar] [CrossRef] [Green Version]
  13. Bari, R.; Jones, J.D. Role of plant hormones in plant defence responses. Plant Mol. Biol. 2009, 69, 473–488. [Google Scholar] [CrossRef]
  14. De Vleesschauwer, D.; Gheysen, G.; Hofte, M. Hormone defense networking in rice: tales from a different world. Trends Plant Sci. 2013, 18, 555–565. [Google Scholar] [CrossRef]
  15. De Vleesschauwer, D.; Xu, J.; Hofte, M. Making sense of hormone-mediated defense networking: from rice to Arabidopsis. Front. Plant Sci. 2014, 5, 611. [Google Scholar] [CrossRef]
  16. Grant, M.R.; Jones, J.D. Hormone (dis)harmony moulds plant health and disease. Science 2009, 324, 750–752. [Google Scholar] [CrossRef]
  17. Nahar, K.; Kyndt, T.; Hause, B.; Hofte, M.; Gheysen, G. Brassinosteroids suppress rice defense against root-knot nematodes through antagonism with the jasmonate pathway. MPMI 2013, 26, 106–115. [Google Scholar] [CrossRef] [Green Version]
  18. Spoel, S.H.; Dong, X. Making sense of hormone crosstalk during plant immune responses. Cell Host Microbe 2008, 3, 348–351. [Google Scholar] [CrossRef] [Green Version]
  19. Gheysen, G.; Mitchum, M.G. How nematodes manipulate plant development pathways for infection. Curr. Opin. Plant Biol. 2011, 14, 415–421. [Google Scholar] [CrossRef]
  20. Ton, J.; Van Pelt, J.A.; Van Loon, L.C.; Pieterse, C.M. Differential effectiveness of salicylate-dependent and jasmonate/ethylene-dependent induced resistance in Arabidopsis. MPMI 2002, 15, 27–34. [Google Scholar] [CrossRef] [Green Version]
  21. Williamson, V.M. Plant nematode resistance genes. Curr. Opin. Plant Biol. 1999, 2, 327–331. [Google Scholar] [CrossRef]
  22. Williamson, V.M.; Kumar, A. Nematode resistance in plants: The battle underground. TIG 2006, 22, 396–403. [Google Scholar] [CrossRef]
  23. Lorieux, M.; Reversat, G.; Garcia Diaz, S.X.; Denance, C.; Jouvenet, N.; Ghesquiere, A. Linkage mapping of Hsa-1(Og), a resistance gene of African rice to the cyst nematode, Heterodera sacchari. TAG. Theoretical and applied genetics. Theor. Appl. Genet. 2003, 107, 691–696. [Google Scholar] [CrossRef]
  24. Wan, J.; Vuong, T.; Jiao, Y.; Joshi, T.; Zhang, H.; Nguyen, H.T. Whole-genome gene expression profiling revealed genes and pathways potentially involved in regulating interactions of soybean with cyst nematode (Heterodera glycines Ichinohe). BMC Genom. 2015, 16, 148. [Google Scholar] [CrossRef] [Green Version]
  25. Barcala, M.; Garcia, A.; Cabrera, J.; Casson, S.; Lindsey, K.; Favery, B.; García-Casado, G.; Solano, R.; Fenoll, C.; Escobar, C. Early transcriptomic events in microdissected Arabidopsis nematode-induced giant cells. Plant J. 2010, 61, 698–712. [Google Scholar] [CrossRef]
  26. Caillaud, M.C.; Dubreuil, G.; Quentin, M.; Perfus-Barbeoch, L.; Lecomte, P.; de Almeida Engler, J. Root-knot nematodes manipulate plant cell functions during a compatible interaction. Plant Physiol. 2008, 165, 104–113. [Google Scholar] [CrossRef]
  27. Nguyen, P.V.; Bellafiore, S.; Petitot, A.S.; Haidar, R.; Bak, A.; Fernandez, D. Meloidogyne incognita—Rrice (Oryza sativa) interaction: A new model system to study plant-root-knot nematode interactions in monocotyledons. Rice 2014, 7, 23. [Google Scholar] [CrossRef] [Green Version]
  28. Delteil, A.; Zhang, J.; Lessard, P.; Morel, J.-B. Potential candidate genes for improving rice disease resistance. Rice 2010, 3, 56. [Google Scholar] [CrossRef]
  29. Kyndt, T.; Zemene, H.Y.; Haeck, A.; Singh, R.; De Vleesschauwer, D.; Gheysen, G. Below-Ground Attack by the Root Knot Nematode Meloidogyne graminicola Predisposes Rice to Blast Disease. MPMI 2017, 30, 255–266. [Google Scholar] [CrossRef] [Green Version]
  30. Gao, Y.; Zhang, C.; Han, X.; Wang, Z.Y.; Ma, L.; Xuan, Y.H. Inhibition of OsSWEET11 function in mesophyll cells improves resistance of rice to sheath blight disease. Mol. Plant Pathol. 2018, 19, 2149–2161. [Google Scholar] [CrossRef] [Green Version]
  31. Yuan, P.; Zhang, C.; Wang, Z.Y.; Zhu, X.F.; Xuan, Y.H. RAVL1 Activates Brassinosteroids and Ethylene Signaling to Modulate Response to Sheath Blight Disease in Rice. Phytopathology 2018, 108, 1104–1113. [Google Scholar] [CrossRef] [Green Version]
  32. Chen, Q.; Peng, D. Nematode Chitin and Application. Adv. Exp. Med. Biol. 2019, 1142, 209–219. [Google Scholar] [CrossRef]
  33. McNeece, B.T.; Sharma, K.; Lawrence, G.W.; Lawrence, K.S.; Klink, V.P. The mitogen activated protein kinase (MAPK) gene family functions as a cohort during the Glycine max defense response to Heterodera glycines. Plant Physiol. Biochem. 2019, 137, 25–41. [Google Scholar] [CrossRef]
  34. Kyndt, T.; Denil, S.; Haegeman, A.; Trooskens, G.; Bauters, L.; Gheysen, G. Transcriptional reprogramming by root knot and migratory nematode infection in rice. N. Phytol. 2012, 196, 887–900. [Google Scholar] [CrossRef]
  35. Nahar, K.; Kyndt, T.; De Vleesschauwer, D.; Hofte, M.; Gheysen, G. The jasmonate pathway is a key player in systemically induced defense against root knot nematodes in rice. Plant Physiol. 2011, 157, 305–316. [Google Scholar] [CrossRef] [Green Version]
  36. Song, L.X.; Xu, X.C.; Wang, F.N.; Wang, Y.; Xia, X.J.; Yu, J.Q. Brassinosteroids act as a positive regulator for resistance against root-knot nematode involving RESPIRATORY BURST OXIDASE HOMOLOG-dependent activation of MAPKs in tomato. Plant Cell Environ. 2018, 41, 1113–1125. [Google Scholar] [CrossRef]
  37. Zhu, X.F.; Liu, Y.; Gai, X.T.; Zhou, Y.; Xia, Z.Y.; Xuan, Y.H. SNARE proteins SYP22 and VAMP727 negatively regulate plant defense. Plant Signal. Behav. 2019, 14, 1610300. [Google Scholar] [CrossRef]
  38. Zhang, H.; Kjemtrup-Lovelace, S.; Li, C.; Luo, Y.; Chen, L.P.; Song, B.H. Comparative RNA-Seq Analysis Uncovers a Complex Regulatory Network for Soybean Cyst Nematode Resistance in Wild Soybean (Glycine soja). Sci. Rep. 2017, 7, 9699. [Google Scholar] [CrossRef] [Green Version]
  39. De la Torre, F.; Gutierrez-Beltran, E.; Pareja-Jaime, Y.; Chakravarthy, S.; Martin, G.B.; del Pozo, O. The tomato calcium sensor Cbl10 and its interacting protein kinase Cipk6 define a signaling pathway in plant immunity. Plant Cell 2013, 25, 2748–2764. [Google Scholar] [CrossRef] [Green Version]
  40. Branch, C.; Hwang, C.-F.; Navarre, D.A.; Williamson, V.M. Salicylic acid is part of the Mi-1-mediated defense response to root-knot nematode in tomato. MPMI 2004, 17, 351–356. [Google Scholar] [CrossRef] [Green Version]
  41. Xuan, Y.H.; Kumar, V.; Han, X.; Kim, S.H.; Jeong, J.H.; Han, C.-D. CBL-INTERACTING PROTEIN KINASE 9 regulates ammonium-dependent root growth downstream of IDD10 in rice (Oryza sativa). Ann. Bot. 2019, 124, 947–960. [Google Scholar] [CrossRef]
  42. Jeong, D.H.; An, S.; Kang, H.G.; Moon, S.; Han, J.J.; Park, S.; Lee, H.S.; An, K.; An, G. T-DNA insertional mutagenesis foractivation tagging in rice. Plant Physiol. 2002, 130, 1636–1644. [Google Scholar]
  43. Hong, Z.; Ueguchi-Tanaka, M.; Umemura, K.; Uozu, S.; Fujioka, S.; Takatsuto, S.; Yoshida, S.; Ashikari, M.; Kitano, H.; Matsuoka, M. A rice brassinosteroid-deficient mutant, ebisu dwarf (d2), is caused by a loss of function of a new member of cytochrome P450. Plant Cell 2003, 15, 2900–2910. [Google Scholar]
  44. Ma, B.; Yin, C.C.; He, S.J.; Lu, X.; Zhang, W.K.; Lu, T.G.; Chen, S.Y. Ethylene-induced inhibition of root growth requires abscisic acid function in rice (Oryza sativa L.) seedlings. PLoS Genet. 2014, 10, e1004701. [Google Scholar]
  45. Bybd, D., Jr.; Kirkpatrick, T.; Barker, K. An improved technique for clearing and staining plant tissues for detection of nematodes. J. Nematol. 1983, 15, 142. [Google Scholar] [CrossRef]
  46. Thies, J.A.; Merrill, S.B.; Corley, E.L. Red food coloring stain: New, safer procedures for staining nematodes in roots and egg masses on root surfaces. J. Nematol. 2002, 34, 179. [Google Scholar]
  47. Trapnell, C.; Hendrickson, D.G.; Sauvageau, M.; Goff, L.; Rinn, J.L.; Pachter, L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat. Biotechnol. 2013, 31, 46–53. [Google Scholar]
  48. Ye, J.; Fang, L.; Zheng, H.; Zhang, Y.; Chen, J.; Wang, J. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34, W293–W297. [Google Scholar] [CrossRef]
  49. Dupuy, D.; Bertin, N.; Hidalgo, C.A.; Venkatesan, K.; Tu, D.; Lee, D. Genome-scale analysis of in vivo spatiotemporal promoter activity in Caenorhabditis elegans. Nat. Biotechnol. 2007, 25, 663–668. [Google Scholar] [CrossRef]
  50. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Figure 1. Meloidogyne incognita stained with acid fuchsin after infection of Nipponbare rice roots at 3 (A), 6 (B), 9 (C), 12 (D), 15 (E), 18 (F), 21 (G), 24 (H), and 27 (I) days post inoculation (dpi). J2: second-stage infective juveniles; one asterisk: sausage-shaped J2; two asterisks: globose root-knot nematode (RKN); ♀: female adult; ♂: male adult.
Figure 1. Meloidogyne incognita stained with acid fuchsin after infection of Nipponbare rice roots at 3 (A), 6 (B), 9 (C), 12 (D), 15 (E), 18 (F), 21 (G), 24 (H), and 27 (I) days post inoculation (dpi). J2: second-stage infective juveniles; one asterisk: sausage-shaped J2; two asterisks: globose root-knot nematode (RKN); ♀: female adult; ♂: male adult.
Ijms 21 00848 g001
Figure 2. Differentially expressed genes (DEGs) induced by root-knot nematodes at 6 and 18 dpi. (A) Cluster heat map showing the level of expression of differential genes; (B) Up- and downregulated genes; (C) Venn diagrams showing intersection of genes that were upregulated at either 6 or 18 dpi; (D) Venn diagrams showing intersection of genes that were downregulated at either 6 or 18 dpi.
Figure 2. Differentially expressed genes (DEGs) induced by root-knot nematodes at 6 and 18 dpi. (A) Cluster heat map showing the level of expression of differential genes; (B) Up- and downregulated genes; (C) Venn diagrams showing intersection of genes that were upregulated at either 6 or 18 dpi; (D) Venn diagrams showing intersection of genes that were downregulated at either 6 or 18 dpi.
Ijms 21 00848 g002
Figure 3. The 20 most enriched Gene Ontology (GO) terms (A and C) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (B and D) for Nipponbare rice roots infected by Meloidogyne incognita at 6 dpi (A and B) and 18 dpi (C and D). “Rich factor” indicates the ratio between the number of DEGs and the number of genes in this pathway. The higher the rich factor, the greater the degree of enrichment.
Figure 3. The 20 most enriched Gene Ontology (GO) terms (A and C) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (B and D) for Nipponbare rice roots infected by Meloidogyne incognita at 6 dpi (A and B) and 18 dpi (C and D). “Rich factor” indicates the ratio between the number of DEGs and the number of genes in this pathway. The higher the rich factor, the greater the degree of enrichment.
Ijms 21 00848 g003
Figure 4. DEG analysis of Nipponbare rice roots infected by M. incognita using MapMan (A) 6dpi; (B) 18dpi.
Figure 4. DEG analysis of Nipponbare rice roots infected by M. incognita using MapMan (A) 6dpi; (B) 18dpi.
Ijms 21 00848 g004
Figure 5. Altered disease resistance of cipk9 and CIPK9 OX plants against RKNs. (AH) CIPK gene family expression. Error bars indicate the SD between biological repeats (n = 3). (I) Numbers of second-stage infective juveniles (J2) inside roots on the cipk9 and CIPK9 OX plants at 6 d after RKN inoculation (SD, n = 16). (J) Numbers of root galls in the cipk9 and CIPK9 OX plants at 18 d after RKN inoculation (SD, n = 16). (K,L) Numbers and percentages of different larval instars on the cipk9 and CIPK9 OX plants at 18 d after RKN inoculation (SD, n = 16); * p < 0.05, ** p < 0.01.
Figure 5. Altered disease resistance of cipk9 and CIPK9 OX plants against RKNs. (AH) CIPK gene family expression. Error bars indicate the SD between biological repeats (n = 3). (I) Numbers of second-stage infective juveniles (J2) inside roots on the cipk9 and CIPK9 OX plants at 6 d after RKN inoculation (SD, n = 16). (J) Numbers of root galls in the cipk9 and CIPK9 OX plants at 18 d after RKN inoculation (SD, n = 16). (K,L) Numbers and percentages of different larval instars on the cipk9 and CIPK9 OX plants at 18 d after RKN inoculation (SD, n = 16); * p < 0.05, ** p < 0.01.
Ijms 21 00848 g005
Figure 6. Altered disease resistance of bri-D and d2 plants against RKNs. (AD) D2, D11, OsBRI1, and OsBAK1 expression. Error bars indicate the SD between biological replicates (n = 3). (E) Numbers of second-stage infective juveniles (J2) inside the roots of Bir-D and d2 mutants at 6 d after RKN inoculation (SD, n = 16). (F) Numbers of root galls on the Bir-D and d2 mutants at 18 d after RKN inoculation (SD, n = 16). (G,H) Numbers and percentages of different larval instars on the Bir-D and d2 mutant sat 18 d after RKN inoculation (SD, n = 16); * p < 0.05, ** p < 0.01.
Figure 6. Altered disease resistance of bri-D and d2 plants against RKNs. (AD) D2, D11, OsBRI1, and OsBAK1 expression. Error bars indicate the SD between biological replicates (n = 3). (E) Numbers of second-stage infective juveniles (J2) inside the roots of Bir-D and d2 mutants at 6 d after RKN inoculation (SD, n = 16). (F) Numbers of root galls on the Bir-D and d2 mutants at 18 d after RKN inoculation (SD, n = 16). (G,H) Numbers and percentages of different larval instars on the Bir-D and d2 mutant sat 18 d after RKN inoculation (SD, n = 16); * p < 0.05, ** p < 0.01.
Ijms 21 00848 g006
Figure 7. Altered disease resistance of ers1 plants against RKNs. (A) Numbers of J2 in the roots of the ers1 mutant at 6 d after RKN inoculation (SD, n = 16). (B) Numbers of root galls in the ers1 mutant at 18 d after RKN inoculation (SD, n = 16). (C,D) Numbers and percentages of different larval instars in the ers1 mutant at 18 d after RKN inoculation (SD, n = 16); * p < 0.05, ** p < 0.01.
Figure 7. Altered disease resistance of ers1 plants against RKNs. (A) Numbers of J2 in the roots of the ers1 mutant at 6 d after RKN inoculation (SD, n = 16). (B) Numbers of root galls in the ers1 mutant at 18 d after RKN inoculation (SD, n = 16). (C,D) Numbers and percentages of different larval instars in the ers1 mutant at 18 d after RKN inoculation (SD, n = 16); * p < 0.05, ** p < 0.01.
Ijms 21 00848 g007

Share and Cite

MDPI and ACS Style

Zhou, Y.; Zhao, D.; Shuang, L.; Xiao, D.; Xuan, Y.; Duan, Y.; Chen, L.; Wang, Y.; Liu, X.; Fan, H.; et al. Transcriptome Analysis of Rice Roots in Response to Root-Knot Nematode Infection. Int. J. Mol. Sci. 2020, 21, 848. https://doi.org/10.3390/ijms21030848

AMA Style

Zhou Y, Zhao D, Shuang L, Xiao D, Xuan Y, Duan Y, Chen L, Wang Y, Liu X, Fan H, et al. Transcriptome Analysis of Rice Roots in Response to Root-Knot Nematode Infection. International Journal of Molecular Sciences. 2020; 21(3):848. https://doi.org/10.3390/ijms21030848

Chicago/Turabian Style

Zhou, Yuan, Di Zhao, Li Shuang, Dongxue Xiao, Yuanhu Xuan, Yuxi Duan, Lijie Chen, Yuanyuan Wang, Xiaoyu Liu, Haiyan Fan, and et al. 2020. "Transcriptome Analysis of Rice Roots in Response to Root-Knot Nematode Infection" International Journal of Molecular Sciences 21, no. 3: 848. https://doi.org/10.3390/ijms21030848

APA Style

Zhou, Y., Zhao, D., Shuang, L., Xiao, D., Xuan, Y., Duan, Y., Chen, L., Wang, Y., Liu, X., Fan, H., & Zhu, X. (2020). Transcriptome Analysis of Rice Roots in Response to Root-Knot Nematode Infection. International Journal of Molecular Sciences, 21(3), 848. https://doi.org/10.3390/ijms21030848

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop