1. Introduction
Sweetpotato is widely cultivated in China and has been used as an important food source, feed material and industrial raw material because it is rich in carbohydrates, nutrients, and functional ingredients [
1]. The planting area and output of sweetpotato in China rank first in the world. However, the continuous planting of sweetpotato for many years has led to the occurrence and accumulation of a variety of diseases. Sweetpotato stem rot is one of the most serious sweetpotato diseases in China in recent years, which occurs in Zhejiang, Guangdong, Hebei, Henan and other provinces [
2,
3]. sweetpotato stem rot occurs seriously in land, which often causes serious seedling death and yield loss, and seriously affects the quality. sweetpotato stem and root rot can occur in all life stages of sweetpotato, such as seedling stage, growing period and even storage stage. The main feature of the disease is that dark brown and waterlogged disease spots are produced on the stem or petiole of sweetpotato. With the aggravation of the disease, the stem segment softens and dissociates, and the plant wilts and dies [
4].
Sweetpotato stem rot was first discovered and reported in the United States in 1974 [
5]. The pathogen causing sweetpotato stem rot was identified as Erwinia chrysanthemi, which was later classified as Pectobacterium according to 16s rDNA sequence analysis [
6]. With the emergence of new identification techniques, some species of the pathogen of stem rot are retained in the genus Pectinobacillus, while other species are classified into the new genus Dickeya [
7]. The identification of the pathogen of sweetpotato stem rot started late in China. Through molecular biology identification and sequencing identification, it was found that E.chrysanthemi and Dickeya dadantii were the main pathogens of sweetpotato in China [
3,
4,
8]. The pathogen of sweetpotato stem rot in Zhejiang Province was also identified as D.dadantii [
9]. D.dadantii has a wide host range and can infect more than 60 kinds of host plants , including ornamental plants such as Convolvulaceae, Compositae and Orchidaceae, and agricultural plants such as Solanaceae, Leguminosae, Gramineae, Dioscorea and Cruciferae, especially important food crops, such as sweetpotato, rice, corn, onion, eggplant, pepper, carrot, tomato and so on [
10]. Therefore, D.dadantii is considered to be one of the top ten plant pathogenic bacteria in molecular plant pathology [
11].
In recent years, sweetpotato stem rot has seriously occurred in growing areas in China, which poses a great threat to the development of sweetpotato industry and has become the bottleneck of the development of sweetpotato industry. At present, single chemical control, agricultural control and other means have been unable to control the harm of sweetpotato stem rot. Moreover, there are no varieties of sweetpotato with high resistance and immunity to stem rot [
9], so there is an urgent need to screen and cultivate sweetpotato resources resistant to stem rot. Therefore, the key to cultivate resistant varieties is to excavate plant disease resistance genes. As a powerful analytical tool, RNA-seq sequencing technology has been widely used in different kinds of plants. There have been many reports on transcriptome of plant pathogens under stress, from which a large number of stress response genes of pathogens have been identified, and their disease resistance mechanism has been discussed. Xiao et al used wheat resistant and susceptible materials for transcriptome sequencing and found that PR5 and PR14 were significantly up-regulated in resistant wheat [
12]. During the response of lettuce to Botrytis cinerea infection, RNA-seq data showed that the expression of genes related to phenylpropane and terpenoid synthesis was up-regulated, while the photosynthetic pathway was inhibited [
13]. A resistance gene SWC was isolated and identified from capsicum by RNA-seq technique. The protein encoded by this gene can recognize the effect factor AvrBs4 released by Xanthomonas [
14]. With the continuous development of RNA-seq sequencing technology, with the help of this technology, we can carry out more detailed and in-depth sequencing of plant transcriptome in response to pathogen infection. By studying the differential expression changes of plant genes, we can mine more and more reliable candidate genes for disease resistance and provide important information for further exploration of the molecular mechanism of plant disease resistance. Therefore, in this paper, using the medium-resistant sweetpotato variety Xushu 48 as material, the expression of stem rot response genes in different stages of stem rot infection was analyzed by transcriptome, and the sweetpotato stem rot resistance genes were excavated to provide theoretical basis for sweetpotato resistance breeding.
2. Materials and Methods
2.1. Plant materials and treatments
Xushu 48 is a fresh sweetpotato variety selected by Xuzhou Agricultural Science Research Institute in Xuhuai region of Jiangsu Province and is resistant to stem rot. sweetpotato stem rot pathogen (
Dickeya dadantii) was collected, isolated, identified and preserved by the Institute of Plant Protection and Microbiology of Zhejiang Academy of Agricultural Sciences, Dd1 and Dd2. The mixed strains of Dd1 and Dd2 were used to identify the resistance of sweetpotato to stem rot in laboratory[
15]. Pathogen inoculation: activate and preserve strain Dd1 and Dd2, pick up single colony on NA medium, culture at 28 ℃ for 36 h, rinse bacterial moss with sterile water on NA plate, and prepare bacterial suspension with concentration of 107/mL for resistance test and acupuncture inoculation[
9].
Resistance identification: the indoor inoculation resistance of vermiculite potato seedlings was identified by the method developed by the Institute of Plant Protection and Microbiology of Zhejiang Academy of Agricultural Sciences. Fill the 500ml Guangkou tissue culture bottle with vermiculite with fully wet sterile water, gently wipe out the wound from the stem base of 5cm with sandpaper with 3-4 leaves, and then cut it into vermiculite, so that part of the wound is not in vermiculite, the other part is exposed to the air, sprinkle the 1ml solution to the scratch of sweetpotato seedlings with a transfer gun, and then place it in a temperature-controlled and light-controlled culture room at 30 ℃. The air relative humidity is 80%. Aseptic water-drenched seedlings were used as control. Tissue culture bottle vermiculite cuttage each glass bottle 1 plant, each treatment 10 plants, set up 3 repeats. After being placed in the culture room, the incidence was observed and recorded day by day, and samples were taken at non-infection (T1), infection 3 days (T2), and infection 6 days (T3).
2.2. RNA isolation cDNA library preparation and RNA-Seq
Total RNA was extracted from three sweetpotato cultivars roots treated at 60℃ and CK using Trizol reagent kit (Invitrogen, Carlsbad, CA,USA) according to the manufacturer’s protocol. RNA quality was assessed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and checked using RNase free agarose gel electrophoresis. After total RNA was extracted, eukaryotic mRNA was enriched by Oligo (dT) beads, while prokaryotic mRNA was enriched by removing rRNA by Ribo-ZeroTM Magnetic Kit (Epicentre, Madison, WI, USA). Then the enriched mRNA was fragmented into short fragments using fragmentation buffer and reverse transcripted into cDNA with random primers. Second-strand cDNA were synthesized by DNA polymerase I, RNase H, dNTP and buffer. Then the cDNA fragments were purified with QiaQuick PCR extraction kit (Qiagen, Venlo, The Netherlands), end repaired, poly(A) added, and ligated to Illumina sequencing adapters. The ligation products were size selected by agarose gel electrophoresis, PCR amplified, and sequenced using Illumina HiSeq2500 by Gene Denovo Biotechnology Co. (Guangzhou, China).
2.3. Analysis of DEGs (differentially expressed genes)
RNAs differential expression analysis was performed by DESeq2 software [
16] between two different groups (and by edgeR [
17] between two samples). The genes/transcripts with the parameter of false discovery rate (FDR) ≤ 0.05 and absolute fold change ≥ 2 were considered as differentially expressed genes/transcripts[
18].
2.4. GO (Gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis
GO[
19] and KEGG[
20] pathways of these DEGs were analyzed to determine the significant abundance of GO terms and metabolic pathways. GOSeqR software package was used for GO enrichment analysis, and KEGG Orthology software was used for KEGG pathway analysis. Functional terms and pathways of DEGs were considered to be significantly different when the P < 0.05. The WEGO software was used to classify the functional classification of the GO annotations.
2.5. qRT-PCR (Real-time quantitative PCR) validation
A number of nine randomly selected unigenes were used for quantitative RT-PCR assays using QuantStudio TM 6 Flex Real-Time PCR System by Thermo Fisher Scientific. Total RNA was extracted from sweetpotato storage roots by Total RNA rapid extraction kit(GK3016,Shanghai Generay Biotech Co.,Ltd.,Shanghai,Chain)and was reverse-transcribed into cDNA using ReverTra Ace® qPCR RT Master Mix with gDNA Remover Kit (FSQ-301,TOYOBO Co.,Ltd.,Japan). The qRT-PCR reaction protocol was conducted using SYBR Green Real Time PCR Master Mix (10 μL), PCR Forward / Reverse Primer (10 μM, 0.5μL), cDNA template (1 μL), ddH2O (8 μL), and ARF (JX177359) was used as internal control (reference gene)[
21]. qRT-PCR primers of each unigene were designed using primer 3 Plus software are shown in
Supplementary Table S1. The data from three independent experiments were calculated with the 2−△△Ct method[
22] .
2.6. Statistical analysis
Means, standard errors and significance analysis of differences were calculated by SPSS 22.0 software. Tukey's multiple comparison test was used to determine the significance of differences between samples. The determination of physical and chemical indexes of samples and transcriptome sequencing carried out at least 3 biological repeats. The figures in this paper were drawn by origin 8.0 and photoshop C5 software.
4. Discussion
Like other crops, the production and cultivation of sweetpotato are infected by various pathogens, and stem rot is one of the common diseases of sweetpotato. The incidence of stem rot seriously affects the yield and nutritional quality of sweetpotato. Sweetpotato stem rot was first discovered and reported in the United States in 1974 [
5]. In recent years, the pathogen of sweetpotato stem rot has been isolated and identified as D.dadantii [
2]. However, there are few studies on the mechanism of stem rot resistance in sweetpotato, and there are few reports on the mining of stem rot resistance genes in sweetpotato. Therefore, in this study, transcriptome was used to determine the gene expression level of Xushu 48, moderately resistant to stem rot, in different stages of stem rot infection.
Plants rely on innate immunity to perceive and resist potential pathogens. Plant innate immune system consists of two related systems: Pathogen-associated molecular patterns-triggered immunity (PTI) and Effector-triggered immunity (ETI). PTI is the first layer of the plant immune system [
23]. Most invasive pathogens can be recognized by a transmembrane plant protein called Pattern recognition receptors (PRRs). After identifying the pathogen, PRRs transmits the immunity signal downstream through the Botrytis-induced kinase 1 (BIK1) [
24], Mitogen-activated protein kinase (MAPK) and Calcium-dependent protein kinase (CDPK) [
25], and then activates the corresponding immune response, such as the accumulation of reactive oxygen species (ROS) and callose. In the transcriptome data, we also found significant differentially expressed Mitogen-activated protein kinase genes and Calcium-dependent protein kinase genes, such as MAPK20 (Ibat.Brg.05A_G019550),MAPK18 (Ibat.Brg.06E_G014790) and CPK1 (Ibat.Brg.01E_G026440). In addition, GO annotation analysis showed that the number of DEGs enriched by catalytic activity (GO:0016194), oxidoreductase activity (GO:0003824) was more in the early stage of stem rot infection, while in the later stage of infection, the number of DEGs enriched in organonitrogen compound biosynthetic process (GO:1901566) and small molecule metabolic process (GO:0044281) was more. The expression pattern of these DEGs is consistent with the biological process of infection by plant resistant pathogens.
With the continuous evolution of pathogens, the effectors of some pathogens can avoid PTI into plants. Accordingly, plants have evolved the second layer of the immune system: ETI [
26,
27] . Plants have evolved a specific resistance gene: r gene. Most R genes encode a receptor that binds to leucine-rich repeats (NLR). NLR receptors can sense toxic proteins directly or indirectly in cells and trigger a series of immune responses to achieve plant disease resistance. Similarly, there are a large number of DEGs related to leucine-rich repeats and their related regulatory kinases in transcriptome data, such as LRR receptor-like serine/threonine-protein kinase (FEI1, Ibat.Brg.01D_G026840), plant intracellular Ras-group-related LRR protein (PIRL2, Ibat.Brg.07D_G004510). These candidate genes may be important regulatory genes for sweetpotato to resist the invasion of stem rot pathogen.
Transcription factor family also plays an important role in the process of plant resistance to pathogens. Among the transcription factors induced by MTI, WRKY transcription factor is the most prominent. 15 WRKY-TFs are strongly induced, including WRKY18, WRKY33 and WRKY40 [
28]. These three WRKY genes are also considered to be important node genes in the WRKY regulatory network [
29]. A total of 475 WRKY family members or transcripts were screened from transcriptome data, including WRKY33 (Ibat.Brg.15B_G007050) and WRKY40 (Ibat.Brg.S022350). In addition, bHLH [
30], ERF [
31], NAC [
32,
33], MYB [
34] and other transcription factor families have also been reported to be related to plant disease resistance or biological stress. A large number of transcription factors, such as bHLH137 (Ibat.Brg.04D_G021250), bHLH162 (Ibat.Brg.06B_G028170), ERF9 (Ibat.Brg.09A_G006830), ERF4 (Ibat.Brg.10F_G018750), MYB73 (Ibat.Brg.01F_G031560), and MYB2 (Ibat.Brg.02A_G015020), were also screened. These candidate transcription factors provide valuable genetic resources for studying the molecular mechanism of stem rot resistance of sweetpotato.
5. Conclusions
In this study, transcriptome sequencing analysis was carried out on the samples at different stages (T1, T2, T3) of stem rot pathogen infecting Xushu 48, and 44839 (T1-VS-T2), 81436 (T1-VS-T2) and 61932 (T2-VS-T3) differentially expressed genes (DEGs) were obtained. The DEGs were mainly concentrated in Alanine, aspartate and glutamate metabolism (ko00250), Carbon fixation in photosynthetic organisms (ko00710) and Amino sugar and nucleotide sugar metabolism (ko00520) and some candidate genes related to plant pathogen infection, such as receptor-like protein kinase (RLK5, RLK7), LRR receptor-like serine/threonine-protein kinase (SERK1), bHLH137, ERF9, MYB73 and NAC053, were screened. The results of this study provide genetic resources for the research on stem rot resistance of sweetpotato and provide theoretical basis for sweetpotato resistance breeding.
Author Contributions
Conceptualization, Y.-G.Z. and M.K.; methodology, W.-H.S., M.-Q.J., Q.L. and L.-R.C.; software, H.-H.J and M.K.; validation, H.-H.J. and ; formal analysis, L.Z.; investigation, M.-Q.J.; resources, Q.L. and M.K.; data curation, C.L.; writing—original draft preparation, C.L.; writing—review and editing, L.Z., W.-H.S. and M.K.; visualization, Z.-Y.Z and L.-R.C.; supervision, Z.-Y.Z.; project administration, Y.-G.Z.; funding acquisition, Q.L. and L.-R.C.. All authors have read and agreed to the published version of the manuscript.
Figure 1.
phenotype of different stages of stem rot infection of Xushu 48. The red arrow indicates that the sweetpotato stem is infected by stem rot. T1, T2 and T3 were uninfected, infected by stem rot for 3 days and infected by stem rot for 6 days, respectively..
Figure 1.
phenotype of different stages of stem rot infection of Xushu 48. The red arrow indicates that the sweetpotato stem is infected by stem rot. T1, T2 and T3 were uninfected, infected by stem rot for 3 days and infected by stem rot for 6 days, respectively..
Figure 2.
Transcriptome sample relationship analysis and RT-qPCR verification. Principal component (A) and heat map analysis(B) analysis of stem rot pathogen infecting stem segments of Xushu 48 at different stages. (C) Correlation analysis of FPKM in transcriptome data and fold change of RT-qPCR.
Figure 2.
Transcriptome sample relationship analysis and RT-qPCR verification. Principal component (A) and heat map analysis(B) analysis of stem rot pathogen infecting stem segments of Xushu 48 at different stages. (C) Correlation analysis of FPKM in transcriptome data and fold change of RT-qPCR.
Figure 3.
Total number of DEGs in stem rot pathogen infecting stem segments of Xushu 48 at different stage. (A) Venn diagram of DEG. The number in each circle represents the number of DEG in the corresponding group. The overlap of the circle represents the common number of DEG between groups. (B) Multi-point difference scatter diagrams of different comparison groups.Red and green dots represent differentially expressed genes up-regulated and down-regulated, respectively.
Figure 3.
Total number of DEGs in stem rot pathogen infecting stem segments of Xushu 48 at different stage. (A) Venn diagram of DEG. The number in each circle represents the number of DEG in the corresponding group. The overlap of the circle represents the common number of DEG between groups. (B) Multi-point difference scatter diagrams of different comparison groups.Red and green dots represent differentially expressed genes up-regulated and down-regulated, respectively.
Figure 4.
GO annotation analysis of DEGs in stem rot pathogen infecting stem segments of Xushu 48 at different stage. The orange column indicates the number of DEGs up-regulated. The blue column indicates the number of DEGs down-regulated.
Figure 4.
GO annotation analysis of DEGs in stem rot pathogen infecting stem segments of Xushu 48 at different stage. The orange column indicates the number of DEGs up-regulated. The blue column indicates the number of DEGs down-regulated.
Figure 5.
Figure 5. Top 20 significantly enriched GO terms. (A) TOP 20 GO terms with significant enrichment of DEGs compared with T1 and T2. (B) TOP 20 GO terms with significant enrichment of DEGs compared with T1 and T3. (C) TOP 20 GO terms with significant enrichment of DEGs compared with T2 and T3.
Figure 5.
Figure 5. Top 20 significantly enriched GO terms. (A) TOP 20 GO terms with significant enrichment of DEGs compared with T1 and T2. (B) TOP 20 GO terms with significant enrichment of DEGs compared with T1 and T3. (C) TOP 20 GO terms with significant enrichment of DEGs compared with T2 and T3.
Figure 6.
KEGG pathway enrichment analysis of DEGs in stem rot pathogen infecting stem segments of Xushu 48 at different stage. (A) Quantitative statistics of DEGs enriched in different metabolic pathways. Orange, blue, red, green and yellow columns represent the number of DEGs enriched in metabolism, genetic information processing, environmental information processing, organismal systems and cellular processes, respectively. The height of the column indicates the number of DEGs. Correlation analysis of metabolic pathway by DEGs in T1 and T2 (B), T1 and T3 (C) and T2 and T3 (D) comparison groups, respectively.
Figure 6.
KEGG pathway enrichment analysis of DEGs in stem rot pathogen infecting stem segments of Xushu 48 at different stage. (A) Quantitative statistics of DEGs enriched in different metabolic pathways. Orange, blue, red, green and yellow columns represent the number of DEGs enriched in metabolism, genetic information processing, environmental information processing, organismal systems and cellular processes, respectively. The height of the column indicates the number of DEGs. Correlation analysis of metabolic pathway by DEGs in T1 and T2 (B), T1 and T3 (C) and T2 and T3 (D) comparison groups, respectively.
Figure 7.
Analysis of expression trend of DEGs in stem rot pathogen infecting stem segments of Xushu 48 at different stage. (A) the significance analysis of DEGs in each module. (B) the number of DEGs in each module.
Figure 7.
Analysis of expression trend of DEGs in stem rot pathogen infecting stem segments of Xushu 48 at different stage. (A) the significance analysis of DEGs in each module. (B) the number of DEGs in each module.
Figure 8.
Quantitative analysis of differentially expressed transcription factors. The number on each column indicates the number of corresponding transcription factor family members or transcripts.
Figure 8.
Quantitative analysis of differentially expressed transcription factors. The number on each column indicates the number of corresponding transcription factor family members or transcripts.